Message: Re: 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

More Re: dE [ GetTotalEnergyDeposit() ] value remains constant at 0.0332 MeV 

Forum: Physics List
Re: Question dE [ GetTotalEnergyDeposit() ] value remains constant at 0.0332 MeV (Ashwath Rao)
Date: 06 Apr, 2012
From: Ashwath Rao <Ashwath Rao>

hello everyone, Also kindly tell me where to add:

G4double eDeposit = step-> GetTotalEnergyDeposit()

I want to obtain energy deposit after every gamma event.

Thanking all.

Regards, Ashwath.

On Thu, 05 Apr 2012 10:38:37 GMT, Ashwath Rao wrote:

> 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 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: "Re: 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 ]