Message: dE [ GetTotalEnergyDeposit() ] value remains constant at 0.0332 MeV Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question dE [ GetTotalEnergyDeposit() ] value remains constant at 0.0332 MeV 

Forum: Physics List
Date: 05 Apr, 2012
From: Ashwath Rao <Ashwath Rao>

Respected Experts,
  I have been trying to simulate interaction process of NaI(Tl) crystal with gamma rays (662 keV). While doing so, I observe that the total energy deposited by the gamma (as displayed on the terminal) is 0.0332 MeV, and it stays constant throughout the simulation.
I am using geant4.9.3.p01
My physics list contains the following :
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 AddTransportation();

 ConstructEM();
 G4PhotoElectricEffect*thePhotoElectricEffect= new G4PhotoElectricEffect();
 G4ComptonScattering * theComptonScattering = new G4ComptonScattering();

 ConstructGeneral();
 pmanager->AddProcess(theDecayProcess);
 pmanager->SetProcessOrdering(theDecayProcess,idxPostStep);
 pmanager->SetProcessOrdering(theDecayProcess, idxAtRest);

 ConstructOp(); [for scintillation]
  SetVerboseLevel(1);
  theScintillationProcess->SetTrackSecondariesFirst(true);
  theScintillationProcess->SetScintillationYieldFactor(0.3);
////////////
  pmanager->AddProcess(theScintillationProcess);
  pmanager->SetProcessOrderingToLast(theScintillationProcess, idxAtRest);
  pmanager->SetProcessOrderingToLast(theScintillationProcess, idxPostStep);
///////////
 if   (particle==G4OpticalPhoton::OpticalPhoton())
{  pmanager->AddDiscreteProcess (new G4OpBoundaryProcess()); }

 SetCuts();
double cut_em=  0.010*mm;
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

Fluctuations/variations in the total energy deposited will only give me the defined NaI spectrum with photopeak at 662 keV.(as desired).

where can i find the error in the codes? and why is it coming a constant value always? I request for guidance to all of you. Thanking all in advance. Regards, Ashwath.

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

1 More: Re: dE [ GetTotalEnergyDeposit() ] value remains constant at 0.0332 MeV   (Ashwath Rao - 06 Apr, 2012)
(_ None: Re: dE [ GetTotalEnergyDeposit() ] value remains constant at 0.0332 MeV   (michel maire - 06 Apr, 2012)
(_ Agree: Re: dE [ GetTotalEnergyDeposit() ] value remains constant at 0.0332 MeV   (Ashwath Rao - 06 Apr, 2012)
(_ Feedback: Re: dE [ GetTotalEnergyDeposit() ] value remains constant at 0.0332 MeV   (Ashwath Rao - 06 Apr, 2012)
 Add Message Add Message
to: "dE [ GetTotalEnergyDeposit() ] value remains constant at 0.0332 MeV"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]