SPIM PhotElectricEffect Lab

The objective of this lab is to evaluate the implementatin of the Photoelectric effect in GEANT4.

The photoelectric effect is a discrete process by which a photon ejects electrons from the surface of a metal.

Here the process is generalized to mean the ejection of an electron due to a collision with an incident photon (I call this process direct ionization). The incident photon is absorbed in the process and the ejected electron carries away the excess energy. This differs from Compton scattering where the photon is not absorbed but rather scatters from the electron.

To accomplish the above objective you will choose an element from the table below which has at least 4 energy levels.

http://brems.physics.isu.edu/~tforest/NucSim/Day9/BindingEnergyTable.pdf

Then you will use GEANT to create histograms of the photon energy lost as simulated by 2 different implementations of the Photoelectric effect. This is done in the PhysicsList. The two processes are added in the physics list for gammas such that you have either

    pmanager->AddDiscreteProcess(new G4PhotoElectricEffect);

or

     pmanager->AddDiscreteProcess(new G4LowEnergyPhotoElectric);

you can't have both at the same time.

Our goal will be to see which one of the above is a more realistic implementation of the Photoelectric effect.

1.) Step 1 is to create a target made from one of the elements above which is long enough to almost guarantee a photoelectric event for each incident photon.

I chose Argon. You will choose something else.

I edited the Detector Construction code to have the following

  //Argon gas
 G4Material* ArgonGas = 
   new G4Material("ArgonGas", z=18., a=39.948*g/mole, density= 1.784*mg/cm3);

and I made my target length long

fTargetLength = 15. * cm;

2.) Now I make sure the two processes are in my Physicslists

be sure the following header files are included in the PhysicsList process definition.

1. include "G4PhotoElectricEffect.hh"
2. include "G4LowEnergyPhotoElectric.hh"

and have at least one physics process "un-commented" out.

   if (particleName == "gamma") {
     // gamma         

pmanager->AddDiscreteProcess(new G4PhotoElectricEffect);

     //pmanager->AddDiscreteProcess(new G4LowEnergyPhotoElectric);
     
   } else if (particleName == "e-") {

3.) Now you need to alter the SteppingVerbose code so it spits out the photon KE lost. Mine looked like this

    if( fTrack->GetDefinition()->GetPDGEncoding()==22 && fStep->GetPostStepPoint()->GetProcessDefinedStep()->GetProcessName() == "phot" && fTrack->GetVolume()->GetName() =="Target")

// if( fTrack->GetDefinition()->GetPDGEncoding()==22 && fStep->GetPostStepPoint()->GetProcessDefinedStep()->GetProcessName() == "LowEnPhotoElec" && fTrack->GetVolume()->GetName() =="Target")

   {
     
     //      G4cout  << "  Photon "  ;
     //G4cout  << 
     outfile <<

// fTrack->GetKineticEnergy() << " " fStep->GetTotalEnergyDeposit()<< " " << fTrack->GetPosition().x()<< " " << fTrack->GetPosition().y()<< " " << fTrack->GetPosition().z()<< " " << fTrack->GetMomentum().x() << " " << fTrack->GetMomentum().y() << " " << fTrack->GetMomentum().z() << " " << G4endl;

   }

Note there are 2 if statements above , one is commented out, which need to change according to your entry in the physics list.

I used G4out and set /tracking/versose 2 to check that I was printing out the right thing for each simulation event (above it is commented out). If you comment our the line with "oufile" and un-comment out the two lines with G4cout then what was printing to a file will print to the terminal window.

4.) You now have the infrastructure to start your investigation.

ie: use /gun/energy to change the energy of the incident photon to cover the range of electron binding energies in the atom you selected.

you can create root trees containing the photon kinetic energy lost in the target and compare those with the electron binding energies of your chosen atom.

Is there a distribution of photon energies lost or is just one specific photon energy lost?

5.) As in the previous lab you will write your results up in TeX, include all neccessary figures, and reach a conclusion describing which entry to use in your physicslist and what are the limitations for its use.

 Note: My grading scheme is as follows

0 points /10 if you use a target atom with less than 3 energy levels (ie H, He, Li ...). This means all atoms should have "A" of Sodium(Na) or above.

8 points/10 if you supply a pdf file with ionization spectrum for the two different physics list entries

9 points /10 if you insert a reference and tabulated values for the binding energy of your atom (along with uncertainty in that binding energy)

10 points /10 if you show via a referenced paper that the sampled binding energies have been cross-section weighted.