|Message: Re: Energy deposited by secondary particles||Not Logged In (login)|
Click on the Forum title, e.g. on the "Forums by Category" page, to read a sequence of postings to the Forum and its threads all in one page. If you are only interested in one thread or the thread following a specific posting, click the thread or the posting, which takes you to a smaller page, which contains only the part you are interested in and may be easier to navigate.
Messages are "chained" if there are only replies at the first level, i.e. 1/1.html, 1/1/1.html etc. In case of "chained" messages the message number is replaced by the icon and there is no indentation.
Inline: Display the subject line only or also the text of the posting(s); for the choice "All" the "Outline" choices are switched off.
|1||0||1||no text / full text of posting|
|2||1||All||text for level 1 only / text for All postings|
Outline: Choose the depth of the posting thread, successive toggle controls provide increasing detail.
|1||2||1||2 levels / 1 level (original posting)|
|2||3||2||3 levels / 2 levels|
|3||3||All||3 levels / all levels (all postings)|
On Thu, 10 Apr 2014 09:05:53 GMT, Ivan Brilkov wrote:
> I'm currently working on similar problem. May be you should see example > named "radioprotection": there are code that collects energies of > secondary particles in root histograms. Look at RemSimSteppingAction.cc. > > // Retrieve the secondary particles > G4TrackVector* fSecondary = steppingManager -> GetfSecondary(); > // Retrieve the info about the generation of secondary particles in the phantom and > // in the vehicle > G4String volumeName = (*fSecondary)[lp1] -> GetVolume() -> GetName(); > G4String secondaryParticleName = (*fSecondary)[lp1]->GetDefinition() -> GetParticleName(); > G4double secondaryParticleKineticEnergy = (*fSecondary)[lp1] -> GetKineticEnergy(); > G4String process = (*fSecondary)[lp1]-> GetCreatorProcess()-> GetProcessName(); > > In another example this code is used for getting Energy deposit, but > it's applicable only for primary particles: > > // collect energy and track length step by step > G4double edep = aStep->GetTotalEnergyDeposit(); > > If you find a solution of this problem, please post here. Ivan >
I also came across the method you suggested where we could simply find the information of the kinetic energy of the secondary particles. Hopefully there is a way to find the total energy deposited for the secondary particles if this is possible for the primary particles. If I find a way I will let you know!