larg4::ISCalculationNEST Class Reference

#include <ISCalculationNEST.h>

Inheritance diagram for larg4::ISCalculationNEST:

Public Member Functions
	ISCalculationNEST (CLHEP::HepRandomEngine &engine)
void	Reset ()
void	CalculateIonizationAndScintillation (const G4Step *step)
double	StepSizeLimit () const
Public Member Functions inherited from larg4::ISCalculation
virtual	~ISCalculation ()=default
double	EnergyDeposit () const
double	NumberIonizationElectrons () const
double	NumberScintillationPhotons () const
double	VisibleEnergyDeposit () const
double	EFieldAtStep (double fEfield, const G4Step *step) const

Private Attributes
NestAlg	fNest
	the fast optical simulation process More...
double	fStepSize
	maximum step to take More...

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

Detailed Description

Definition at line 22 of file ISCalculationNEST.h.

Constructor & Destructor Documentation

larg4::ISCalculationNEST::ISCalculationNEST

(

CLHEP::HepRandomEngine &

engine

)

Definition at line 19 of file ISCalculationNEST.cxx.

    : // \todo should ideally make the yield factor passed to the NestAlg ctor a parameter
    fNest(1., engine)
  {
    // 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 LArSoft uses a fudge factor to compensate, but NEST is "purer" -- no
    // fudge factor. "
    fStepSize = 0.05 * CLHEP::micrometer;
  }

Member Function Documentation

void larg4::ISCalculationNEST::CalculateIonizationAndScintillation ( const G4Step *  step )

virtual

Implements larg4::ISCalculation.

Definition at line 42 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() != step->GetTotalEnergyDeposit() / CLHEP::MeV)
      mf::LogWarning("ISCalculationNest")
        << "NEST and G4 step depositions do not agree!\n"
        << fNest.EnergyDeposition() << " vs " << step->GetTotalEnergyDeposit() / CLHEP::MeV;

    // Nest uses Geant units, LArSoft assumes energy is in units of MeV here
    fEnergyDeposit = fNest.EnergyDeposition() / CLHEP::MeV;
    fNumIonElectrons = fNest.NumberIonizationElectrons();
    fNumScintPhotons = fNest.NumberScintillationPhotons();
    fVisibleEnergyDeposition = 0.; //Not implimented for NEST.
  }

void larg4::ISCalculationNEST::Reset ( )

virtual

Implements larg4::ISCalculation.

Definition at line 32 of file ISCalculationNEST.cxx.

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

double larg4::ISCalculationNEST::StepSizeLimit ( ) const

inlinevirtual

Implements larg4::ISCalculation.

Definition at line 29 of file ISCalculationNEST.h.

    {
      return fStepSize;
    }

Member Data Documentation

NestAlg larg4::ISCalculationNEST::fNest

private

the fast optical simulation process

Definition at line 35 of file ISCalculationNEST.h.

double larg4::ISCalculationNEST::fStepSize

private

maximum step to take

Definition at line 36 of file ISCalculationNEST.h.

---

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

• ISCalculationNEST.h
• ISCalculationNEST.cxx