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()