gar::rosim::ISCalculationNEST Class Reference

#include <ISCalculationNEST.h>

Inheritance diagram for gar::rosim::ISCalculationNEST:

Public Member Functions
	ISCalculationNEST (CLHEP::HepRandomEngine &engine)
virtual	~ISCalculationNEST ()
void	Initialize ()
void	Reset ()
void	CalculateIonizationAndScintillation (const gar::sdp::EnergyDeposit *dep)
double	EnergyDeposit () const
int	NumberIonizationElectrons () const
int	NumberScintillationPhotons () const
double	StepSizeLimit () const
Public Member Functions inherited from gar::rosim::ISCalculation
	ISCalculation ()
virtual	~ISCalculation ()

Private Attributes
NestAlg *	fNest
	the fast optical simulation process More...
double	fStepSize
	maximum step to take More...
CLHEP::HepRandomEngine &	fEngine
	random engine More...

Additional Inherited Members
Protected Attributes inherited from gar::rosim::ISCalculation
double	fEnergyDeposit
	total energy deposited in the step More...
int	fNumIonElectrons
	number of ionization electrons for this step More...
int	fNumScintPhotons
	number of scintillation photons for this step More...

Detailed Description

Definition at line 20 of file ISCalculationNEST.h.

Constructor & Destructor Documentation

gar::rosim::ISCalculationNEST::ISCalculationNEST

(

CLHEP::HepRandomEngine &

engine

)

Definition at line 18 of file ISCalculationNEST.cxx.

    : ISCalculation()
    , fNest(0)
    , fEngine(engine)
    {
      return;
    }

gar::rosim::ISCalculationNEST::~ISCalculationNEST ( )

virtual

Definition at line 27 of file ISCalculationNEST.cxx.

    {
      if(fNest) delete fNest;
      
      return;
    }

Member Function Documentation

void gar::rosim::ISCalculationNEST::CalculateIonizationAndScintillation ( const gar::sdp::EnergyDeposit *  dep )

virtual

Implements gar::rosim::ISCalculation.

Definition at line 60 of file ISCalculationNEST.cxx.

    {
      // get a const representation of the track for this step
      //const G4Track track(*(step->GetTrack()));
      
        //fNest->CalculateIonizationAndScintillation(track, *step);
      
      // compare the energy deposition of this step to what is in the fNest object
      if(fNest->EnergyDeposition() != dep->Energy() / CLHEP::GeV)
        MF_LOG_WARNING("ISCalculationNest")
        << "NEST and G4 step depositions do not agree!\n"
        << fNest->EnergyDeposition()
        << " vs "
        << dep->Energy() / CLHEP::GeV;
    
      // Nest uses Geant units, GArSoft assumes energy is in units of GeV here
      fEnergyDeposit   = fNest->EnergyDeposition() / CLHEP::GeV;
      fNumIonElectrons = fNest->NumberIonizationElectrons();
      fNumScintPhotons = fNest->NumberScintillationPhotons();
      
      return;
    }

double gar::rosim::ISCalculationNEST::EnergyDeposit ( ) const

inlinevirtual

Implements gar::rosim::ISCalculation.

Definition at line 30 of file ISCalculationNEST.h.

{ return fEnergyDeposit;   }

void gar::rosim::ISCalculationNEST::Initialize ( )

virtual

Implements gar::rosim::ISCalculation.

Definition at line 35 of file ISCalculationNEST.cxx.

    {
      // \todo should ideally make the yield factor passed to the NestAlg ctor a parameter
      if(!fNest) fNest = new NestAlg(1., fEngine);
      
      // Set the step size to small value if NEST is chosen, per Matthew Szydagis,
      // "because without delta rays, the yields are wrong.  The ICARUS model that is
      // in GArSoft uses a fudge factor to compensate, but NEST is "purer" -- no
      // fudge factor. "
      fStepSize = 0.05 * CLHEP::micrometer;
      
      return;
    }

int gar::rosim::ISCalculationNEST::NumberIonizationElectrons ( ) const

inlinevirtual

Implements gar::rosim::ISCalculation.

Definition at line 31 of file ISCalculationNEST.h.

{ return fNumIonElectrons; }

int gar::rosim::ISCalculationNEST::NumberScintillationPhotons ( ) const

inlinevirtual

Implements gar::rosim::ISCalculation.

Definition at line 32 of file ISCalculationNEST.h.

{ return fNumScintPhotons; }

void gar::rosim::ISCalculationNEST::Reset ( )

virtual

Implements gar::rosim::ISCalculation.

Definition at line 50 of file ISCalculationNEST.cxx.

    {
      fEnergyDeposit   = 0.;
      fNumIonElectrons = 0.;
      fNumScintPhotons = 0.;
      
      return;
    }

double gar::rosim::ISCalculationNEST::StepSizeLimit ( ) const

inlinevirtual

Implements gar::rosim::ISCalculation.

Definition at line 33 of file ISCalculationNEST.h.

{ return fStepSize;        }

Member Data Documentation

CLHEP::HepRandomEngine& gar::rosim::ISCalculationNEST::fEngine

private

random engine

Definition at line 39 of file ISCalculationNEST.h.

NestAlg* gar::rosim::ISCalculationNEST::fNest

private

the fast optical simulation process

Definition at line 37 of file ISCalculationNEST.h.

double gar::rosim::ISCalculationNEST::fStepSize

private

maximum step to take

Definition at line 38 of file ISCalculationNEST.h.

---

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

• garsoft/ReadoutSimulation/ISCalculationNEST.h
• garsoft/ReadoutSimulation/ISCalculationNEST.cxx