SPIM PhotElectricEffect Lab - Revision history

Foretony at 18:32, 31 March 2025

2025-03-31T18:32:44Z

Foretony at 17:23, 31 March 2025

2025-03-31T17:23:37Z

Foretony at 17:18, 31 March 2025

2025-03-31T17:18:07Z

Foretony at 17:17, 31 March 2025

2025-03-31T17:17:47Z

Foretony at 02:37, 17 March 2025

2025-03-17T02:37:46Z

Foretony at 02:35, 17 March 2025

2025-03-17T02:35:34Z

Foretony at 02:35, 17 March 2025

2025-03-17T02:35:01Z

Foretony at 21:51, 22 February 2019

2019-02-22T21:51:17Z

Foretony at 20:20, 7 October 2015

2015-10-07T20:20:56Z

Foretony at 22:39, 24 October 2013

2013-10-24T22:39:33Z

← Older revision		Revision as of 18:32, 31 March 2025
Line 117:		Line 117:
	binding energy) and compare that value with the GEANT4 prediction.		binding energy) and compare that value with the GEANT4 prediction.

−	10 points /10 if you are able to implement the ~~G4PenelopePhotoElectricModel~~() and ~~describe what happens to~~ the ~~ejected electron in this model~~.	+	10 points /10 if you are able to implement the G4EmLivermorePhysics() and G4EmPenelopePhysics() models. Compare the energy levels among the three models and the Measured/Mosley values.

@@ Line 53: / Line 53: @@
 out the photon KE lost.  Mine looked like this
-     if( fTrack->GetDefinition()->GetPDGEncoding()==22 &&
+  if(step->GetTrack()->GetParticleDefinition()->GetParticleName()=="gamma" && step->GetPostStepPoint()->GetProcessDefinedStep()->GetProcessName()=="phot")
-                fStep->GetPostStepPoint()->GetProcessDefinedStep()->GetProcessName() ==  "phot"
-              && fTrack->GetVolume()->GetName() =="Target")
+        /*
+          G4cout
-      //      G4cout  << "  Photon "  ;
+          << "stepping:cc "
-      //G4cout  <<
+          << " \n\t Particle=" << step->GetTrack()->GetParticleDefinition()->GetParticleName()
-      outfile <<
+          << " \n \t Process Name=" << step->GetPostStepPoint()->GetProcessDefinedStep()->GetProcessName()
-           //   fTrack->GetKineticEnergy() << "    "
+          << "\n \t ParentID=  " << step->GetTrack()->GetParentID()
-            fStep->GetTotalEnergyDeposit()<< "    "
+          << "\n\t volume=" << step->GetTrack()->GetVolume()->GetName()
-              << fTrack->GetPosition().x()<< "    "
+          << "\n \t x,y,z=" << step->GetTrack()->GetPosition().x()<< ","
-              << fTrack->GetPosition().y()<< "    "
+          <<step->GetTrack()->GetPosition().y() << ","
-             << fTrack->GetPosition().z()<< "    "
+          << step->GetTrack()->GetPosition().z() << ","
-             << fTrack->GetMomentum().x() << "    "
+          << "\n \t Px,Py,Pz=" << step->GetTrack()->GetMomentum().x() << ", "
-             << fTrack->GetMomentum().y() << "    "
+          << step->GetTrack()->GetMomentum().y() << ", "
-             << fTrack->GetMomentum().z() << "    "
+          << step->GetTrack()->GetMomentum().z() << ", "
-            << G4endl;<br>
+          << "\n \t dz=" << step->GetStepLength()
+          << " \n \t Delta E=" <<  step->GetTotalEnergyDeposit()
+          << " \n\t edep=" << edep
+          << G4endl;
 I used G4out and set /tracking/versose 2 to check that I was

@@ Line 14: / Line 14: @@
 To accomplish the above objective you will choose an element from
 the table below which has at least 4 energy levels.
 https://wiki.iac.isu.edu/index.php?title=TF_SPIM_e-gamma#Moseley.27s_Law

← Older revision		Revision as of 17:17, 31 March 2025
Line 17:		Line 17:
	https://wiki.iac.isu.edu/index.php/TF_SPIM_e-gamma#Moseley.27s_Law		https://wiki.iac.isu.edu/index.php/TF_SPIM_e-gamma#Moseley.27s_Law

		+	https://wiki.iac.isu.edu/index.php?title=TF_SPIM_e-gamma#Moseley.27s_Law

	Then you will use GEANT to create histograms of the photon energy		Then you will use GEANT to create histograms of the photon energy

@@ Line 1: / Line 1: @@
 The objective of this lab is to evaluate the implementation of the
-Photoelectric effect in GEANT4.
+Photoelectric (photo) effect in GEANT4.
 The photoelectric effect is a discrete process by which a photon

← Older revision		Revision as of 02:35, 17 March 2025
Line 114:		Line 114:

	10 points /10 if you are able to implement the G4PenelopePhotoElectricModel() and describe what happens to the ejected electron in this model.		10 points /10 if you are able to implement the G4PenelopePhotoElectricModel() and describe what happens to the ejected electron in this model.
−
−
−
−	~~Adding the Penelope model~~
−	~~<pre>~~
−
−	~~G4PhotoElectricEffect* thePhotoElectricEffect = new G4PhotoElectricEffect();~~
−
−
−	~~G4PenelopePhotoElectricModel* thePEPenelopeModel = new G4PenelopePhotoElectricModel();~~
−	~~thePEPenelopeModel->SetHighEnergyLimit(PenelopeHighEnergyLimit);~~
−	~~thePhotoElectricEffect->AddEmModel(0,thePEPenelopeModel);~~
−	~~ph->RegisterProcess(thePhotoElectricEffect, particle);~~
−	~~</pre>~~

@@ Line 35: / Line 35: @@
     //Argon gas
-   G4Material* ArgonGas =
+   //G4Material* ArgonGas =
-    new G4Material("ArgonGas", z=18., a=39.948*g/mole, density= 1.784*mg/cm3);
+  //  new G4Material("ArgonGas", z=18., a=39.948*g/mole, density= 1.784*mg/cm3);
 and I made my target length long
@@ Line 42: / Line 43: @@
   fTargetLength  = 15. * cm;
-.) Then make sure that only one process is in the PhysicsList for the gamma particle
+.) Then make sure that only one process is in the physics.cc for the gamma particle
 <pre>
+     RegisterPhysics (new G4EmStandardPhysics());
-    if (particleName == "gamma") {
+       //    RegisterPhysics(new G4EmLivermorePhysics());
-      // gamma
+       //    RegisterPhysics(new G4EmPenelopePhysics());
 </pre>
-.)  Now you need to alter the SteppingVerbose code so it writes
+.)  Now you need to alter the stepping.cc code so it writes
 out the photon KE lost.  Mine looked like this

← Older revision		Revision as of 21:51, 22 February 2019
Line 117:		Line 117:

	10 points /10 if you are able to implement the G4PenelopePhotoElectricModel() and describe what happens to the ejected electron in this model.		10 points /10 if you are able to implement the G4PenelopePhotoElectricModel() and describe what happens to the ejected electron in this model.
		+
		+
		+
		+	Adding the Penelope model
		+	<pre>
		+
		+	G4PhotoElectricEffect* thePhotoElectricEffect = new G4PhotoElectricEffect();
		+
		+
		+	G4PenelopePhotoElectricModel* thePEPenelopeModel = new G4PenelopePhotoElectricModel();
		+	thePEPenelopeModel->SetHighEnergyLimit(PenelopeHighEnergyLimit);
		+	thePhotoElectricEffect->AddEmModel(0,thePEPenelopeModel);
		+	ph->RegisterProcess(thePhotoElectricEffect, particle);
		+	</pre>

@@ Line 20: / Line 20: @@
 Then you will use GEANT to create histograms of the photon energy
 lost as predicted by GEANT4's
-Photoelectric effect.  Your Physics list should only have the process
+Photoelectric effect.  Your Physics list should only have the G4PhotoElectricEffecr physics process.
 you should comment out all other physics processes for a gamma particle.