G4XTankSteppingAction.cc

Go to the documentation of this file.

#include "G4XTankSteppingAction.h"
#include "G4XTankTrackingAction.h"
#include "G4XTankSimulator.h"

#include <G4Step.hh>

#include <G4VProcess.hh>
#include <sstream>

using namespace G4XTankSimulatorAG;
using namespace std;

int G4XTankSteppingAction::fgParticleId = 0;

G4XTankSteppingAction::G4XTankSteppingAction() 
{
}

G4XTankSteppingAction::~G4XTankSteppingAction() 
{
}

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

  fgParticleId = theStep->GetTrack()->GetDefinition()->GetPDGEncoding();

  if (fgParticleId == 0) {

    if (theStep->GetTrack()->GetNextVolume() != 0) {
            
      if (theStep->GetTrack()->GetNextVolume()->GetName() == "top" || 
          theStep->GetTrack()->GetNextVolume()->GetName() == "side" ||
          theStep->GetTrack()->GetNextVolume()->GetName() == "bottom") {

        G4XTankTrackingAction::fNumBounces++;
        
      } 
    }    
  }



  // Particle data must be taken when they reach the muon counter's depth

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

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

  // Check if the particle entered the ground level
  if( prePVName == "ground"){

    // Convert the depth to milimeters
    double depth = fCounterDepth*CLHEP::mm;

    // Check if the particle reached the counter's depth
    // The last condition (z > depth*2) is to avoid problems with neutrinos 
    if(z <= depth &&  z0 > depth  && z > depth*2) {

      G4int    trackid   = theTrack->GetTrackID();
      
/*
      G4cerr << "\n===========================================\n";
      G4cerr << "        NEW PARTICLE ARRIVED TO COUNTER DEPTH of " << depth  << " mm !!!!            \n";
      G4cerr <<  particleType->GetParticleName() << G4endl;
      G4cerr <<  particleType->GetParticleName() << G4endl;
      G4cerr <<  "PDG code: " << particleType->GetPDGEncoding()  << G4endl;
      G4cerr << "prePV: " << prePVName << G4endl;
      G4cerr << "trackID " << trackid << G4endl;
      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 );

    }

  }



}