Difference between revisions of "HRRL Simulations"

From New IAC Wiki
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.


HRRL ExperimentCell World 1.jpg

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


Go back