How to get track starting point

#1

Hello.

I made detector and I need to get coordinates where secondary track start. I want to know how can I bring coordinates for separate each components.

Below is code I tried.

if(aTrack->GetTrackID() != 1)

{

if(aTrack->GetCreatorProcess()->GetProcessName() != “proptonInelastic”)

{

G4double positionZ = step->GetPreStepPoint()->GetPosition().z();

G4double r = positionZ;
}
{

#2

Hi,

If you want proton inelastic process, you should make sure the correct process’s name.
proptonInelastic => protonInelastic

Cheers,
Ye

#3

Thank you.

I correct process name proptonInelastic -> protonInelastic and I get result.

However, still some problem is remain. I measure how many time the process occur two method. One is using my code(already fixed) and another is count manually. but in detector, result from two method is different. result from code is have more event than second method.

#4

I found cause of problem and solve it.

It bring next step starting point not interaction point.
I change my code
G4double positionZ = step->GetPreStepPoint()->GetPosition().z();
to
G4double pz = aTrack->GetVertexPosition().z();

then it bring interaction point.

1 Like