Difference between revisions of "HRRL Simulations"
| Line 16: | Line 16: | ||
[[Image:HRRL_ExperimentCell_World_1.jpg | 200 px]] | [[Image:HRRL_ExperimentCell_World_1.jpg | 200 px]] | ||
| + | == Run 1 == | ||
I then ran 1 million events in which an incident 7 MeV electron hit the 2 mm thick Tungsten with the physics processes | I then ran 1 million events in which an incident 7 MeV electron hit the 2 mm thick Tungsten with the physics processes | ||
| Line 32: | Line 33: | ||
I only kept event in which had a momentum component towards the ceiling | I only kept event in which had a momentum component towards the ceiling | ||
| + | The energy distribution (in MeV) of the gammas headed towards the roof is shown below. Of the 1 million electrons incident on the 2 mm thick Tungsten target, only 70,000 gammas were headed towards the cieling. This does not mean that they hit it, they could have hit the wall. | ||
| + | |||
| + | [[Image:Run_1_Gamma_E_dist.gif]] | ||
[http://www.iac.isu.edu/mediawiki/index.php/HRRL Go back] | [http://www.iac.isu.edu/mediawiki/index.php/HRRL Go back] | ||
Revision as of 11:31, 11 July 2008
This document describes the simulations performed to estimate the radiological background in the event that a 1 nA electron current from the HRRL enters the experimental cell.
Radiation monitors: Ludlum Model 45-8.
Low energy gamma cutoff = 60 keV
Tungsten SImulation
I created a world volume filled with Air to represent the experimental cell.
I created a 2mm thick tungsten target that is 30 cm x 30 cm in area.
The image below shows several electrons hitting the tungsten foil, then scattering in air. The red lines are electrons and the green ones are photons.
Run 1
I then ran 1 million events in which an incident 7 MeV electron hit the 2 mm thick Tungsten with the physics processes
if (particleName == "gamma") {
// gamma
pmanager->AddDiscreteProcess(new G4PhotoElectricEffect);
pmanager->AddDiscreteProcess(new G4ComptonScattering);
pmanager->AddDiscreteProcess(new G4GammaConversion);
} else if (particleName == "e-") {
//electron
pmanager->AddProcess(new G4MultipleScattering,-1, 1,1);
pmanager->AddProcess(new G4eIonisation, -1, 2,2);
pmanager->AddProcess(new G4eBremsstrahlung, -1, 3,3);
I only kept event in which had a momentum component towards the ceiling
The energy distribution (in MeV) of the gammas headed towards the roof is shown below. Of the 1 million electrons incident on the 2 mm thick Tungsten target, only 70,000 gammas were headed towards the cieling. This does not mean that they hit it, they could have hit the wall.
