StepPrinter.cc

Go to the documentation of this file.

//~~~ SLIC ~~~//
#include "StepPrinter.hh" 
 
namespace slic {
 
void StepPrinter::stepping(const G4Step *step) {
  // Get the track associated with this step
  auto track{step->GetTrack()};
 
  auto trackID{track->GetTrackID()}; 
  if ((trackID_ > 0) && (trackID != trackID_)) return;
 
  // Get the particle name.
  auto particleName{track->GetParticleDefinition()->GetParticleName()};
 
  // Get the energy of the particle
  auto energy{step->GetPostStepPoint()->GetTotalEnergy()};
 
  // Get the volume the particle is in.
  auto volume{track->GetVolume()->GetName()};
 
  // Get the next volume
  auto nextVolume{track->GetNextVolume()->GetName()};
 
  // Get the region
  auto region{track->GetVolume()->GetLogicalVolume()->GetRegion()->GetName()};
 
  std::cout << "*******************************\n"
            << "*   Step " << track->GetCurrentStepNumber() << "\n"
            << "********************************\n"
            << "\tEnergy of " << particleName << " : " << energy << "\n"
            << "\tTrack ID: " << track->GetTrackID() << "\n"
            << "\tStep #: " << track->GetCurrentStepNumber() << "\n"
            << "\tParticle currently in " << volume << "\n"
            << "\tRegion: " << region << "\n"
            << "\tNext volume: " << nextVolume << "\n"
            << "********************************\n"
            << std::endl;
}
} // namespace slic
 
DECLARE_PLUGIN(slic, StepPrinter)