|Message: Re: primary particle number||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 Sun, 31 Jan 2010 11:56:05 GMT, riazi wrote: > Dear all > I need to count the number of primary protons which can reach to a fantom or detector after traveling a beamline. would you please help me to do this?. > > thanks a lot for your time and consideration riazi > As explained in Application Developers Guide (5.1.2), each tracked particle has an ID, the primary particle being ID=1
For instance in TestEm1, SteppingAction::UserSteppingAction() you will find :
// step length of primary particle G4int ID = aStep->GetTrack()->GetTrackID(); G4double steplen = aStep->GetStepLength(); if (ID == 1) ...
(TestEm1 in examples/extended/electromagnetic)
To know where the particle is :
// get volume of the current step G4VPhysicalVolume* volume = aStep->GetPreStepPoint()->GetTouchableHandle()->GetVolume(); if (volume == detector->GetAbsorber()) ...
More generaly, the attached text (in FAQ), may help.