Message: Running IAEA phase space files as particle generator Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Running IAEA phase space files as particle generator 

Forum: Run Management
Date: 01 Jun, 2016
From: Brett Cunningham <Brett Cunningham>

Hello Geant4 community,

I am currently simulating a mammography x-ray tube, and to aid in computational efficiency downstream I am collecting phase space files at different locations. To test my writing and reading of those files, I'm running a 10^7 original event simulation whereby I'm collecting about 10^3 photons at my phase space plane. However, when I read the file back in as the particle generator I get, what I would consider, some unexpected results.

To read the phase space file back in I'm using the command /run/beamOn, but the number of events I have to use to get 10^3 photons out from the file is 10^7 - i.e. /run/beamOn 10000000 when reading in the phase space file results in all 10^3 from the phase space file being loaded in. This is not what I would expect, what I would expect is when using a phase space file with 10^3 photons a command of /run/beamOn 1000, and NOT /run/beamOn 10000000, would load the entire phase space file, but it doesn't.

So my question is, is this typical behavior when reading in phase space files as particle generators? or am I missing something totally obvious? or is there some other command I need to be using with a phase space file?

Another somewhat related issue, when I recycle photons using theIAEAReader->SetTimesRecycled(9) all it does is multiply the energy of the particle from the phase space file by the number of times it was recycled. i.e. if I had SetTimesRecycled(9) to recycle 10 times and I were to detect the particles using a detection volume, I would see one entry of 170keV from a 17keV phase space photon.

Thankyou in advance to anyone that may be able to help, this is one of the last hurdles I need to get over to start collecting data.

Brett Cunningham

P.S. I haven't posted any code as it is essentially the same as that found in the IAEA "HowTo" guide.

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

1 Note: Re: Running IAEA phase space files as particle generator   (Miguel Cortes-Giraldo - 08 Jun, 2016)
2 Note: Re: Running IAEA phase space files as particle generator   (Miguel Cortes-Giraldo - 08 Jun, 2016)
1 None: Re: Running IAEA phase space files as particle generator   (Brett Cunningham - 14 Jun, 2016)
2 None: Re: Running IAEA phase space files as particle generator   (Juliana Cristina Martins - 17 Jan, 2018)
1 Warning: Re: Running IAEA phase space files as particle generator   (Miguel Cortes-Giraldo - 24 Jan, 2018)
... 1 Message(s)
 Add Message Add Message
to: "Running IAEA phase space files as particle generator"

 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 ]