Message: Re: How score Kinetic Energy? Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Re: How score Kinetic Energy? 

Forum: Biasing and Scoring
Re: Question How score Kinetic Energy? (lu luo)
Date: 05 May, 2014
From: Mahdi <Mahdi>

Dear all; I have such a same problem. I want to simulate energy spectrum of X and UV rays that generated via atomic deexcitation caused by an alpha particle source interaction with air or any other materials. I want this spectrum in world region. the spatial distribution is not important to me.It is some thing like f8:p tally in MCNP/X.

for this problem can I use the example Testem4 or RE03 ?

I used TestEm4 example and modified its physics list to consider alpha particle. Now I have a histogram of energy spectrum but I cant filter this spectrum to be just for gamma particle not electrons and etc...

And in RE03; can I use /score/filter/particleWithKineticEnergy and define bins to have the gamma rays spectrum.because of in this example the filtering of particles sems to be easy.

Could you please advise me abot this?

Inline Depth:
 0 0
 1 1
Add message: (add)

None Re: How score Kinetic Energy?  by Tsukasa Aso <Tsukasa Aso>,   11 May, 2014
Re: Question Re: How score Kinetic Energy? (Mahdi)
Hi,

   In the command-line scorer, there is no filter to find a particle generation.
i.e. If you count the tracks such as current, the count is incremented every steps.
If you want score the track at only generated point, you need to ensure that the particle step is 1.  

  Unfortunately, the command-line scorer does not have such filter.
Therefore you have to create it by yourself. 
The G4PSMinKinEAtGeneration class in the digits_hits directory in Geant4 source code will be helpful.

  Regards,
  Tsukasa Aso

None Re: How score Kinetic Energy?  by michel maire <michel maire>,   07 May, 2014
Re: Question Re: How score Kinetic Energy? (Mahdi)
On Mon, 05 May 2014 20:08:19 GMT, Mahdi wrote:
> Dear all; I have such a same problem. I want to simulate energy spectrum
> of X and UV rays that generated via atomic deexcitation caused by an
> alpha particle source interaction with air or any other materials. I
> want this spectrum in world region. the spatial distribution is not
> important to me.It is some thing like f8:p tally in MCNP/X.
> 
> for this problem can I use the example Testem4 or RE03 ?
> 
> I used TestEm4 example and modified its physics list to consider alpha
> particle. Now I have a histogram of energy spectrum but I cant filter
> this spectrum to be just for gamma particle not electrons and etc...

 In steppingAction, you can do
 G4ParticleDefinition* particle = aStep->GetTrack()->GetGefinition();
 if (particle == G4Gamma::Gamma()) fEventAction->addEdep(EdepStep);

 But you will be surprised by the result.

 Add Message Add Message
to: "Re: How score Kinetic Energy?"

 Subscribe Subscribe

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