Difference between revisions of "HRRL Simulations"
Jump to navigation
Jump to search
Line 15: | Line 15: | ||
[[Image:HRRL_ExperimentCell_World_1.jpg | 200 px]] | [[Image:HRRL_ExperimentCell_World_1.jpg | 200 px]] | ||
+ | |||
+ | I then ran 1 million events 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 | ||
+ | |||
+ | |||
[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:27, 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.
I then ran 1 million events 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