gar::garg4::EnergyDepositAction Class Reference

list of energy deposits from Geant4 More...

#include <EnergyDepositAction.h>

Inheritance diagram for gar::garg4::EnergyDepositAction:

Public Member Functions
	EnergyDepositAction (CLHEP::HepRandomEngine *engine, fhicl::ParameterSet const &pset)
virtual	~EnergyDepositAction ()
void	reconfigure (fhicl::ParameterSet const &pset)
void	BeginOfEventAction (const G4Event *)
void	EndOfEventAction (const G4Event *)
void	PreTrackingAction (const G4Track *)
void	PostTrackingAction (const G4Track *)
void	SteppingAction (const G4Step *)
std::string	GetVolumeName (const G4Track *track)
std::vector< gar::sdp::EnergyDeposit > const &	EnergyDeposits () const

Private Member Functions
void	AddEnergyDeposition (const G4Step *step)

Private Attributes
double	fEnergyCut
	be included in the list. More...
std::vector< std::string >	fVolumeName
	volume we will record energy depositions in More...
std::string	fMaterialMatchString
	Energy deposition will be recorded for materials that match this. More...
std::vector< gar::sdp::EnergyDeposit >	fDeposits
	energy fDeposits More...

Detailed Description

list of energy deposits from Geant4

Definition at line 35 of file EnergyDepositAction.h.

Constructor & Destructor Documentation

gar::garg4::EnergyDepositAction::EnergyDepositAction	(	CLHEP::HepRandomEngine *	engine,
		fhicl::ParameterSet const &	pset
	)

Definition at line 38 of file EnergyDepositAction.cxx.

      {
        this->reconfigure(pset);
      }

gar::garg4::EnergyDepositAction::~EnergyDepositAction ( )

virtual

Definition at line 47 of file EnergyDepositAction.cxx.

      {
        // Delete anything that we created with "new'.
      }

Member Function Documentation

void gar::garg4::EnergyDepositAction::AddEnergyDeposition ( const G4Step *  step )

private

Definition at line 147 of file EnergyDepositAction.cxx.

      {
        // get the track id for this step
        auto trackID  = ParticleListAction::GetCurrentTrackID();

        // first check that we are not dealing with ionization, ie delta rays
        // if we have one of those, set the trackID to be the ID of the parent
        //if(trackID < 0)
        //  MF_LOG_VERBATIM("EnergyDepositAction")
        //  << "TrackID "
        //  << trackID
        //  << " was created by process "
        //  << step->GetTrack()->GetCreatorProcess()->GetProcessName();
        //else
        //  MF_LOG_VERBATIM("EnergyDepositAction")
        //  << "TrackID "
        //  << trackID
        //  << " is a primary particle ";

        // the step mid point is used for the position of the deposit
        auto midPoint = 0.5 * (step->GetPreStepPoint()->GetPosition() +
        step->GetPostStepPoint()->GetPosition() );
        float time = step->GetPreStepPoint()->GetGlobalTime();

        fDeposits.emplace_back(trackID,
          time,
          step->GetTotalEnergyDeposit() * CLHEP::MeV / CLHEP::GeV,
          midPoint.x() / CLHEP::cm,
          midPoint.y() / CLHEP::cm,
          midPoint.z() / CLHEP::cm,
          step->GetStepLength() / CLHEP::cm,
          (trackID > 0));

          return;
        }

void gar::garg4::EnergyDepositAction::BeginOfEventAction

(

const G4Event *

)

Definition at line 67 of file EnergyDepositAction.cxx.

      {
        // Clear any previous information.
        fDeposits.clear();
      }

void gar::garg4::EnergyDepositAction::EndOfEventAction

(

const G4Event *

)

Definition at line 184 of file EnergyDepositAction.cxx.

        {
          // sort the EnergyDeposit lists in each EnergyDeposits object
          std::sort(fDeposits.begin(), fDeposits.end());
        }

std::vector<gar::sdp::EnergyDeposit> const& gar::garg4::EnergyDepositAction::EnergyDeposits ( ) const

inline

Definition at line 57 of file EnergyDepositAction.h.

{ return fDeposits; }

std::string gar::garg4::EnergyDepositAction::GetVolumeName

(

const G4Track *

track

)

Definition at line 191 of file EnergyDepositAction.cxx.

        {
          std::string VolName = track->GetVolume()->GetName();
          VolName.erase(VolName.length()-3, 3);
          return VolName;
        }

void gar::garg4::EnergyDepositAction::PostTrackingAction

(

const G4Track *

)

Definition at line 79 of file EnergyDepositAction.cxx.

      {
      }

void gar::garg4::EnergyDepositAction::PreTrackingAction

(

const G4Track *

)

Definition at line 74 of file EnergyDepositAction.cxx.

      {
      }

void gar::garg4::EnergyDepositAction::reconfigure

(

fhicl::ParameterSet const &

pset

)

Definition at line 53 of file EnergyDepositAction.cxx.

      {
        fEnergyCut  = pset.get< double >("EnergyCut");
        fVolumeName = pset.get< std::vector<std::string> >("GArVolumeName");
        fMaterialMatchString = pset.get< std::string >("GArMaterialMatchString");

        // std::cout << "EnergyDepositAction: Name of the Volumes to track for the GArTPC" << std::endl;
        // for(unsigned int i = 0; i < fVolumeName.size(); i++) std::cout << fVolumeName.at(i) << " ";
        // std::cout << std::endl;

        return;
      }

void gar::garg4::EnergyDepositAction::SteppingAction

(

const G4Step *

step

)

Definition at line 85 of file EnergyDepositAction.cxx.

      {
        MF_LOG_DEBUG("EnergyDepositAction")
        << "EnergyDepositAction::SteppingAction";

        art::ServiceHandle<geo::GeometryGAr> geo;
        TGeoManager *geomanager = geo->ROOTGeoManager();

        // Get the pointer to the track
        G4Track *track = step->GetTrack();

        const CLHEP::Hep3Vector &start = step->GetPreStepPoint()->GetPosition();
        const CLHEP::Hep3Vector &stop  = track->GetPosition();

        // If it's a null step, don't use it.
        if(start == stop) return;
        if(step->GetTotalEnergyDeposit() == 0) return;

        // check that we are in the correct material to record a hit
        std::string VolumeName   = this->GetVolumeName(track);

        if( std::find( fVolumeName.begin(), fVolumeName.end(), VolumeName ) == fVolumeName.end() )
        return;

        // check the material
        auto pos = 0.5 * (start + stop);
        TGeoNode *node = geomanager->FindNode(pos.x()/CLHEP::cm, pos.y()/CLHEP::cm, pos.z()/CLHEP::cm);//Node in cm

        if(!node){
          MF_LOG_DEBUG("EnergyDepositAction")
          << "Node not found in "
          << pos.x() << " mm "
          << pos.y() << " mm "
          << pos.z() << " mm";
          return;
        }

        std::string volmaterial = node->GetMedium()->GetMaterial()->GetName();
        if ( ! std::regex_match(volmaterial, std::regex(fMaterialMatchString)) ) return;

        // only worry about energy depositions larger than the minimum required
        if(step->GetTotalEnergyDeposit() * CLHEP::MeV / CLHEP::GeV > fEnergyCut ){

          MF_LOG_DEBUG("EnergyDepositAction")
          << "In volume "
          << VolumeName
          << " step size is "
          << step->GetStepLength() / CLHEP::cm
          << " cm and deposited "
          << step->GetTotalEnergyDeposit()
          << " MeV of energy with a minimum of "
          << fEnergyCut
          << " required.";

          // save this deposition
          this->AddEnergyDeposition(step);

        } // end if enough energy to worry about this step

      }// end of EnergyDepositAction::SteppingAction

Member Data Documentation

std::vector<gar::sdp::EnergyDeposit> gar::garg4::EnergyDepositAction::fDeposits

private

energy fDeposits

Definition at line 68 of file EnergyDepositAction.h.

double gar::garg4::EnergyDepositAction::fEnergyCut

private

be included in the list.

The minimum energy in GeV for a deposit to

Definition at line 63 of file EnergyDepositAction.h.

std::string gar::garg4::EnergyDepositAction::fMaterialMatchString

private

Energy deposition will be recorded for materials that match this.

Definition at line 66 of file EnergyDepositAction.h.

std::vector<std::string> gar::garg4::EnergyDepositAction::fVolumeName

private

volume we will record energy depositions in

Definition at line 65 of file EnergyDepositAction.h.

---

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

• garsoft/GArG4/EnergyDepositAction.h
• garsoft/GArG4/EnergyDepositAction.cxx