71 if ((aStep->GetPreStepPoint()->GetPhysicalVolume() ==
Detector->
GetWorld())&&
74 (aStep->GetPostStepPoint()->GetProcessDefinedStep()->GetProcessName() ==
"Transportation"))
76 G4double gx, gy, gz, ge, gpx, gpy, gpz, ex, ey, ez, ee;
77 G4double epx, epy, epz, px, py, pz, pe, ppx, ppy, ppz;
81 if(aStep->GetTrack()->GetDynamicParticle()->GetDefinition()->GetParticleName()
84 ex = (aStep->GetTrack()->GetPosition().x())/
cm;
85 ey = (aStep->GetTrack()->GetPosition().y())/
cm;
86 ez = (aStep->GetTrack()->GetPosition().z())/
cm;
88 ee = (aStep->GetTrack()->GetKineticEnergy())/
MeV;
90 epx = aStep->GetTrack()->GetMomentum().x();
91 epy = aStep->GetTrack()->GetMomentum().y();
92 epz = aStep->GetTrack()->GetMomentum().z();
99 if (aStep->GetTrack()->GetDynamicParticle()->GetDefinition()->
100 GetParticleName() ==
"gamma")
104 gx = (aStep->GetTrack()->GetPosition().x())/
cm;
105 gy = (aStep->GetTrack()->GetPosition().y())/
cm;
106 gz = (aStep->GetTrack()->GetPosition().z())/
cm;
109 ge = (aStep->GetTrack()->GetKineticEnergy())/
MeV;
112 gpx = aStep->GetTrack()->GetMomentum().x();
113 gpy = aStep->GetTrack()->GetMomentum().y();
114 gpz = aStep->GetTrack()->GetMomentum().z();
122 if (aStep->GetTrack()->GetDynamicParticle()->GetDefinition()->GetParticleName() ==
"e+")
126 px = (aStep->GetTrack()->GetPosition().x())/
cm;
127 py = (aStep->GetTrack()->GetPosition().y())/
cm;
128 pz = (aStep->GetTrack()->GetPosition().z())/
cm;
131 pe = (aStep->GetTrack()->GetKineticEnergy())/
MeV;
134 ppx = aStep->GetTrack()->GetMomentum().x();
135 ppy = aStep->GetTrack()->GetMomentum().y();
136 ppz = aStep->GetTrack()->GetMomentum().z();
const G4VPhysicalVolume * GetMeasureVolume() const
void fill_Tuple_Gamma(G4double, G4double, G4double, G4double, G4double, G4double, G4double)
const PurgMagDetectorConstruction * Detector
const G4VPhysicalVolume * GetWorld() const
void fill_Tuple_Positrons(G4double, G4double, G4double, G4double, G4double, G4double, G4double)
Uncomment #include to switch to ROOT or XML output file.
void fill_Tuple_Electrons(G4double, G4double, G4double, G4double, G4double, G4double, G4double)
static PurgMagAnalysisManager * getInstance()