MdGroundPropagatorAG/SteppingAction.cc

Go to the documentation of this file.

#include "SteppingAction.h"

#include <fwk/CentralConfig.h>
#include <utl/ErrorLogger.h>
#include <utl/Reader.h>
#include <utl/MathConstants.h>

//G4 classes
#include <G4SteppingManager.hh>
#include <G4Track.hh>
#include <G4Step.hh>
#include <G4StepPoint.hh>
#include <G4TrackStatus.hh>
#include <G4VPhysicalVolume.hh>
#include <G4ParticleDefinition.hh>
#include <G4ParticleTypes.hh>

#include <string>
using std::string;

using namespace fwk;
using namespace utl;
using namespace GroundPropagatorAG;

SteppingAction::SteppingAction  ( double depth ){
fDepth = depth;
}

SteppingAction::~SteppingAction (  ){;}

void SteppingAction::UserSteppingAction ( const G4Step * theStep )
{

  G4Track  * theTrack = theStep->GetTrack (  );

  G4StepPoint       * thePostPoint = theStep->GetPostStepPoint (  );
  G4VPhysicalVolume * thePostPV = thePostPoint->GetPhysicalVolume (  );
  G4ParticleDefinition * particleType = theTrack->GetDefinition (  );

  G4String process = "Undefined";
  if ( thePostPoint->GetProcessDefinedStep() != NULL)
   process = thePostPoint->GetProcessDefinedStep()->GetProcessName();

  if ( !thePostPV ) {//the post step volume is out of the world

    G4double x =  theStep->GetPostStepPoint()->GetPosition().x();
    G4double y =  theStep->GetPostStepPoint()->GetPosition().y();
    G4double z =  theStep->GetPostStepPoint()->GetPosition().z();
    G4VPhysicalVolume *thePrePV = theStep->GetPreStepPoint (  )->GetPhysicalVolume (  );
    string  prePVName  = string (thePrePV->GetName());
    string  postPVName = "OutOfWorld";

    //Check if it a particle entering the counter level
    if( prePVName == "SoilPhys"&& fabs(fDepth+z)< 1e-5 ) {
      G4int    trackid   = theTrack->GetTrackID();
      /*
      G4cerr << "\n===========================================\n";
      G4cerr << "        NEW PARTICLE ARRIVED TO COUNTER DEPTH!!!!            \n";
      G4cerr <<  particleType->GetParticleName() << G4endl;
      G4cerr <<  "PDG code: " << particleType->GetPDGEncoding()  << G4endl;
      G4cerr << "prePV: " << prePVName << G4endl;
      G4cerr << "trackID " << trackid << G4endl;
      G4cerr << "fDepth " << fDepth/CLHEP::m << " (m)" << std::endl;
      G4cerr << "X:" << x/CLHEP::m << " m" << G4endl;
      G4cerr << "Y:" << y/CLHEP::m << " m" << G4endl;
      G4cerr << "Z:" << z/CLHEP::m << " m" << G4endl;
      G4cerr << "===========================================\n";
      */                             
      fPData.Pname   = particleType->GetParticleName();
      fPData.PDGCode = particleType->GetPDGEncoding();
      fPData.TrackID = trackid;
      fPData.x       = x;
      fPData.y       = y;
      fPData.z       = z;
      fPData.px      = theStep->GetPreStepPoint()->GetMomentumDirection().x();
      fPData.py      = theStep->GetPreStepPoint()->GetMomentumDirection().y();
      fPData.pz      = theStep->GetPreStepPoint()->GetMomentumDirection().z();
      fPData.KinE    = theStep->GetPreStepPoint()->GetKineticEnergy();
      fPData.Time    = theStep->GetPreStepPoint()->GetGlobalTime();

      fPartColl.push_back( fPData );

    }
    return;//out of world
  }

 return;

}