PurgMagSteppingAction Class Reference

#include <PurgMagSteppingAction.hh>

Inheritance diagram for PurgMagSteppingAction:

Public Member Functions
	PurgMagSteppingAction (PurgMagDetectorConstruction *)
	~PurgMagSteppingAction ()
void	UserSteppingAction (const G4Step *)
Public Member Functions inherited from G4UserSteppingAction
	G4UserSteppingAction ()
virtual	~G4UserSteppingAction ()
void	SetSteppingManagerPointer (G4SteppingManager *pValue)

Additional Inherited Members
Protected Attributes inherited from G4UserSteppingAction
G4SteppingManager *	fpSteppingManager

Detailed Description

Definition at line 52 of file PurgMagSteppingAction.hh.

Constructor & Destructor Documentation

PurgMagSteppingAction::PurgMagSteppingAction

(

PurgMagDetectorConstruction *

det

)

Definition at line 53 of file PurgMagSteppingAction.cc.

:Detector(det)
{ }

PurgMagSteppingAction::~PurgMagSteppingAction

(

)

Definition at line 59 of file PurgMagSteppingAction.cc.

{ }

Member Function Documentation

void PurgMagSteppingAction::UserSteppingAction ( const G4Step *  aStep )

virtual

Reimplemented from G4UserSteppingAction.

Definition at line 64 of file PurgMagSteppingAction.cc.

References python.hepunit::cm, PurgMagAnalysisManager::fill_Tuple_Electrons(), PurgMagAnalysisManager::fill_Tuple_Gamma(), PurgMagAnalysisManager::fill_Tuple_Positrons(), G4DynamicParticle::GetDefinition(), G4Track::GetDynamicParticle(), PurgMagAnalysisManager::getInstance(), G4Track::GetKineticEnergy(), PurgMagDetectorConstruction::GetMeasureVolume(), G4Track::GetMomentum(), G4Track::GetNextVolume(), G4ParticleDefinition::GetParticleName(), G4StepPoint::GetPhysicalVolume(), G4Track::GetPosition(), G4Step::GetPostStepPoint(), G4Step::GetPreStepPoint(), G4StepPoint::GetProcessDefinedStep(), G4VProcess::GetProcessName(), G4Step::GetTrack(), PurgMagDetectorConstruction::GetWorld(), python.hepunit::MeV, CLHEP::Hep3Vector::x(), CLHEP::Hep3Vector::y(), and CLHEP::Hep3Vector::z().

{ 
  //Collection at SSD in N-tuples. Electrons and photons separated
  //Prestep point in World, next volume MeasureVolume, process transportation
 PurgMagAnalysisManager* analysis = PurgMagAnalysisManager::getInstance();  
  if ((aStep->GetPreStepPoint()->GetPhysicalVolume() == Detector->GetWorld())&&
      (aStep->GetTrack()->GetNextVolume() == Detector->GetMeasureVolume())&&
      //(aStep->GetTrack()->GetMomentumDirection().z()>0.)&& // only particles with positive momentum
      (aStep->GetPostStepPoint()->GetProcessDefinedStep()->GetProcessName() == "Transportation"))
    {
      G4double gx, gy, gz, ge, gpx, gpy, gpz, ex, ey, ez, ee; 
      G4double epx, epy, epz, px, py, pz, pe, ppx, ppy, ppz;
   
         
      // Electrons 
      if(aStep->GetTrack()->GetDynamicParticle()->GetDefinition()->GetParticleName()
      == "e-")
    {//Position 
         ex = (aStep->GetTrack()->GetPosition().x())/cm;
         ey = (aStep->GetTrack()->GetPosition().y())/cm; 
         ez = (aStep->GetTrack()->GetPosition().z())/cm;
     // Energy
     ee = (aStep->GetTrack()->GetKineticEnergy())/MeV;
     // Momentum
     epx = aStep->GetTrack()->GetMomentum().x();
     epy = aStep->GetTrack()->GetMomentum().y();
     epz = aStep->GetTrack()->GetMomentum().z();
      
     // Fill N-tuple electrons
     analysis->fill_Tuple_Electrons(ex, ey, ez, ee, epx, epy, epz);
      }

      // Photons
      if (aStep->GetTrack()->GetDynamicParticle()->GetDefinition()->
          GetParticleName() == "gamma")
    {
      
      // Position
          gx = (aStep->GetTrack()->GetPosition().x())/cm;
      gy = (aStep->GetTrack()->GetPosition().y())/cm;
      gz = (aStep->GetTrack()->GetPosition().z())/cm;
      
      // Energy
      ge = (aStep->GetTrack()->GetKineticEnergy())/MeV;
      
      // Momentum
      gpx = aStep->GetTrack()->GetMomentum().x();
      gpy = aStep->GetTrack()->GetMomentum().y();
      gpz = aStep->GetTrack()->GetMomentum().z();

      // Fill N-tuple photons
      analysis->fill_Tuple_Gamma(gx, gy, gz, ge, gpx, gpy, gpz);
    }


      // Positrons
      if (aStep->GetTrack()->GetDynamicParticle()->GetDefinition()->GetParticleName() == "e+")
    {

      // Position
          px = (aStep->GetTrack()->GetPosition().x())/cm;
      py = (aStep->GetTrack()->GetPosition().y())/cm;
      pz = (aStep->GetTrack()->GetPosition().z())/cm;
      
      // Energy
      pe = (aStep->GetTrack()->GetKineticEnergy())/MeV;
      
      // Momentum
      ppx = aStep->GetTrack()->GetMomentum().x();
      ppy = aStep->GetTrack()->GetMomentum().y();
      ppz = aStep->GetTrack()->GetMomentum().z();
      
      // Fill Ntuple positrons
      analysis->fill_Tuple_Positrons(px, py, pz, pe, ppx, ppy, ppz);
    }
    }
}

---

The documentation for this class was generated from the following files:

• PurgMagSteppingAction.hh
• PurgMagSteppingAction.cc