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
double	EnergyDeposit () const
double	NumberIonizationElectrons () const
double	NumberScintillationPhotons () const
double	VisibleEnergyDeposit () const
double	EFieldAtStep (double fEfield, const G4Step *step) const

Protected Attributes
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

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

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.

References fStepSize.

    : // \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 40 of file ISCalculationNEST.cxx.

References NestAlg::CalculateIonizationAndScintillation(), NestAlg::EnergyDeposition(), larg4::ISCalculation::fEnergyDeposit, fNest, larg4::ISCalculation::fNumIonElectrons, larg4::ISCalculation::fNumScintPhotons, larg4::ISCalculation::fVisibleEnergyDeposition, NestAlg::NumberIonizationElectrons(), NestAlg::NumberScintillationPhotons(), and track.

  {
    // 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.
  }

double larg4::ISCalculation::EFieldAtStep ( double  fEfield, const G4Step *  step  ) const

inherited

Definition at line 18 of file ISCalculation.cxx.

Referenced by larg4::ISCalculationSeparate::CalculateIonizationAndScintillation(), larg4::ISCalculationCorrelated::CalculateIonizationAndScintillation(), and larg4::ISCalculation::VisibleEnergyDeposit().

  {
    auto const* SCE = lar::providerFrom<spacecharge::SpaceChargeService>();
    if (!SCE->EnableSimEfieldSCE()) return efield;
    geo::Point_t midPoint{
      (step->GetPreStepPoint()->GetPosition() + step->GetPostStepPoint()->GetPosition()) * 0.5 /
      CLHEP::cm};
    auto const eFieldOffsets = SCE->GetEfieldOffsets(midPoint);
    return efield * std::hypot(1 + eFieldOffsets.X(), eFieldOffsets.Y(), eFieldOffsets.Z());
  }

double larg4::ISCalculation::EnergyDeposit ( ) const

inlineinherited

Definition at line 22 of file ISCalculation.h.

References larg4::ISCalculation::fEnergyDeposit.

{ return fEnergyDeposit; }

double larg4::ISCalculation::NumberIonizationElectrons ( ) const

inlineinherited

Definition at line 23 of file ISCalculation.h.

References larg4::ISCalculation::fNumIonElectrons.

{ return fNumIonElectrons; }

double larg4::ISCalculation::NumberScintillationPhotons ( ) const

inlineinherited

Definition at line 24 of file ISCalculation.h.

References larg4::ISCalculation::fNumScintPhotons.

{ return fNumScintPhotons; }

void larg4::ISCalculationNEST::Reset ( )

virtual

Implements larg4::ISCalculation.

Definition at line 31 of file ISCalculationNEST.cxx.

References larg4::ISCalculation::fEnergyDeposit, larg4::ISCalculation::fNumIonElectrons, larg4::ISCalculation::fNumScintPhotons, and larg4::ISCalculation::fVisibleEnergyDeposition.

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

double larg4::ISCalculationNEST::StepSizeLimit ( ) const

inlinevirtual

Implements larg4::ISCalculation.

Definition at line 28 of file ISCalculationNEST.h.

{ return fStepSize; }

double larg4::ISCalculation::VisibleEnergyDeposit ( ) const

inlineinherited

Definition at line 25 of file ISCalculation.h.

References larg4::ISCalculation::EFieldAtStep(), and larg4::ISCalculation::fVisibleEnergyDeposition.

{ return fVisibleEnergyDeposition; }

Member Data Documentation

double larg4::ISCalculation::fEnergyDeposit

protectedinherited

total energy deposited in the step

Definition at line 32 of file ISCalculation.h.

Referenced by larg4::ISCalculationSeparate::CalculateIonizationAndScintillation(), CalculateIonizationAndScintillation(), larg4::ISCalculationCorrelated::CalculateIonizationAndScintillation(), larg4::ISCalculation::EnergyDeposit(), larg4::ISCalculationSeparate::Reset(), Reset(), and larg4::ISCalculationCorrelated::Reset().

NestAlg larg4::ISCalculationNEST::fNest

private

the fast optical simulation process

Definition at line 31 of file ISCalculationNEST.h.

Referenced by CalculateIonizationAndScintillation().

double larg4::ISCalculation::fNumIonElectrons

protectedinherited

number of ionization electrons for this step

Definition at line 33 of file ISCalculation.h.

Referenced by larg4::ISCalculationSeparate::CalculateIonizationAndScintillation(), CalculateIonizationAndScintillation(), larg4::ISCalculationCorrelated::CalculateIonizationAndScintillation(), larg4::ISCalculation::NumberIonizationElectrons(), larg4::ISCalculationSeparate::Reset(), Reset(), and larg4::ISCalculationCorrelated::Reset().

double larg4::ISCalculation::fNumScintPhotons

protectedinherited

number of scintillation photons for this step

Definition at line 34 of file ISCalculation.h.

Referenced by larg4::ISCalculationSeparate::CalculateIonizationAndScintillation(), CalculateIonizationAndScintillation(), larg4::ISCalculationCorrelated::CalculateIonizationAndScintillation(), larg4::ISCalculation::NumberScintillationPhotons(), larg4::ISCalculationSeparate::Reset(), Reset(), and larg4::ISCalculationCorrelated::Reset().

double larg4::ISCalculationNEST::fStepSize

private

maximum step to take

Definition at line 32 of file ISCalculationNEST.h.

Referenced by ISCalculationNEST().

double larg4::ISCalculation::fVisibleEnergyDeposition

protectedinherited

Definition at line 35 of file ISCalculation.h.

Referenced by larg4::ISCalculationSeparate::CalculateIonizationAndScintillation(), CalculateIonizationAndScintillation(), Reset(), and larg4::ISCalculation::VisibleEnergyDeposit().

---

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

• larsim/v09_43_00/source/larsim/LegacyLArG4/ISCalculationNEST.h
• larsim/v09_43_00/source/larsim/LegacyLArG4/ISCalculationNEST.cxx