#include <ISCalculationCorrelated.h>

Inheritance diagram for larg4::ISCalculationCorrelated:

Public Member Functions
	ISCalculationCorrelated (detinfo::DetectorPropertiesData const &detProp)

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 Member Functions
double	EscapingEFraction (double const dEdx)

double	FieldCorrection (double const EF, double const dEdx)

Private Attributes
double	fStepSize
	maximum step to take More...

double	fEfield
	value of electric field from LArProperties service More...

double	fWion
	W_ion (23.6 eV) == 1/fGeVToElectrons. More...

double	fWph
	W_ph (19.5 eV) More...

double	fScintPreScale
	scintillation pre-scale from LArProperties service More...

double	fRecombA
	from LArG4Parameters service More...

double	fRecombk
	from LArG4Parameters service More...

double	fModBoxA
	from LArG4Parameters service More...

double	fModBoxB
	from LArG4Parameters service More...

double	fLarqlChi0A
	from LArG4Parameters service More...

double	fLarqlChi0B
	from LArG4Parameters service More...

double	fLarqlChi0C
	from LArG4Parameters service More...

double	fLarqlChi0D
	from LArG4Parameters service More...

double	fLarqlAlpha
	from LArG4Parameters service More...

double	fLarqlBeta
	from LArG4Parameters service More...

bool	fUseModBoxRecomb
	from LArG4Parameters service More...

bool	fUseModLarqlRecomb
	from LArG4Parameters service 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 30 of file ISCalculationCorrelated.h.

Constructor & Destructor Documentation

larg4::ISCalculationCorrelated::ISCalculationCorrelated ( detinfo::DetectorPropertiesData const & detProp )

explicit

Definition at line 39 of file ISCalculationCorrelated.cxx.

   {
     std::cout << "LegacyLArG4/ISCalculationCorrelated Initialize." << std::endl;
     art::ServiceHandle<sim::LArG4Parameters const> lgpHandle;
     const detinfo::LArProperties* larp = lar::providerFrom<detinfo::LArPropertiesService>();
 
     double density = detProp.Density(detProp.Temperature());
     fEfield = detProp.Efield();
     fScintPreScale = larp->ScintPreScale();
 
     // ionization work function
     fWion = 1. / lgpHandle->GeVToElectrons() * 1e3; // MeV
 
     // ion+excitation work function (\todo: get from LArG4Parameters or LArProperties?)
     fWph = 19.5 * 1e-6; // MeV
 
     // the recombination coefficient is in g/(MeVcm^2), but
     // we report energy depositions in MeV/cm, need to divide
     // Recombk from the LArG4Parameters service by the density
     // of the argon we got above.
     fRecombA = lgpHandle->RecombA();
     fRecombk = lgpHandle->Recombk() / density;
     fModBoxA = lgpHandle->ModBoxA();
     fModBoxB = lgpHandle->ModBoxB() / density;
     fLarqlChi0A = lgpHandle->LarqlChi0A();
     fLarqlChi0B = lgpHandle->LarqlChi0B();
     fLarqlChi0C = lgpHandle->LarqlChi0C();
     fLarqlChi0D = lgpHandle->LarqlChi0D();
     fLarqlAlpha = lgpHandle->LarqlAlpha();
     fLarqlBeta  = lgpHandle->LarqlBeta();
     fUseModBoxRecomb = lgpHandle->UseModBoxRecomb();
     fUseModLarqlRecomb = lgpHandle->UseModLarqlRecomb();
 
     // determine the step size using the voxel sizes
     art::ServiceHandle<sim::LArVoxelCalculator const> lvc;
     double maxsize =
       std::max(lvc->VoxelSizeX(), std::max(lvc->VoxelSizeY(), lvc->VoxelSizeZ())) * CLHEP::cm;
 
     fStepSize = 0.1 * maxsize;
   }

Member Function Documentation

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

virtual

Implements larg4::ISCalculation.

Definition at line 95 of file ISCalculationCorrelated.cxx.

   {
     fEnergyDeposit = step->GetTotalEnergyDeposit() / CLHEP::MeV;
 
     // calculate total quanta (ions + excitons)
     double Nq = fEnergyDeposit / fWph;
 
     // Get the recombination factor for this voxel - Nucl.Instrum.Meth.A523:275-286,2004
     // R = A/(1 + (dE/dx)*k)
     // dE/dx is given by the voxel energy deposition, but have to convert it to MeV/cm
     // from GeV/voxel width
     // A = 0.800 +/- 0.003
     // k = (0.097+/-0.001) g/(MeVcm^2)
     // the dx depends on the trajectory of the step
     // k should be divided by the density as our dE/dx is in MeV/cm,
     // the division is handled in the constructor when we set fRecombk
     // B.Baller: Add Modified Box recombination - ArgoNeuT result submitted to JINST
 
     G4ThreeVector totstep = step->GetPostStepPoint()->GetPosition();
     totstep -= step->GetPreStepPoint()->GetPosition();
     double dx = totstep.mag() / CLHEP::cm;
     double dEdx = (dx == 0.0) ? 0.0 : fEnergyDeposit / dx;
     double EFieldStep = EFieldAtStep(fEfield, step);
 
     // Guard against spurious values of dE/dx. Note: assumes density of LAr
     if (dEdx < 1.) dEdx = 1.;
 
     // calculate the recombination survival fraction
     double recomb = 0.;
     if (fUseModBoxRecomb) {
       if (dx) {
         double Xi = fModBoxB * dEdx / EFieldStep;
         recomb = log(fModBoxA + Xi) / Xi;
       }
       else
         recomb = 0;
     }
     else {
       recomb = fRecombA / (1. + dEdx * fRecombk / EFieldStep);
     }
 
     if(fUseModLarqlRecomb){ //Use corrections from LArQL model
       recomb += EscapingEFraction(dEdx)*FieldCorrection(EFieldStep, dEdx); //Correction for low EF
     }
 
     // using this recombination, calculate number of ionization electrons
     fNumIonElectrons = (fEnergyDeposit / fWion) * recomb;
 
     // calculate scintillation photons
     fNumScintPhotons = Nq - fNumIonElectrons;
 
     // apply the scintillation pre-scaling (normally this is already folded into
     // the particle-specific scintillation yields)
     fNumScintPhotons *= fScintPreScale;
 
     MF_LOG_DEBUG("ISCalculationCorrelated")
       << " Electrons produced for " << fEnergyDeposit << " MeV deposited with " << recomb
       << " recombination: " << fNumIonElectrons;
     MF_LOG_DEBUG("ISCalculationCorrelated") << "number photons: " << fNumScintPhotons;
 
     return;
   }

double larg4::ISCalculationCorrelated::EscapingEFraction ( double const dEdx )

private

Definition at line 158 of file ISCalculationCorrelated.cxx.

                                                                     { //LArQL chi0 function = fraction of escaping electrons
     return fLarqlChi0A/(fLarqlChi0B+exp(fLarqlChi0C+fLarqlChi0D*dEdx));
   }

double larg4::ISCalculationCorrelated::FieldCorrection	(	double const	EF,
		double const	dEdx
	)

private

Definition at line 162 of file ISCalculationCorrelated.cxx.

                                                                                    { //LArQL f_corr function = correction factor for electric field dependence
     return exp(-EF/(fLarqlAlpha*log(dEdx)+fLarqlBeta));
   }

void larg4::ISCalculationCorrelated::Reset ( )

virtual

Implements larg4::ISCalculation.

Definition at line 83 of file ISCalculationCorrelated.cxx.

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

double larg4::ISCalculationCorrelated::StepSizeLimit ( ) const

inlinevirtual

Implements larg4::ISCalculation.

Definition at line 37 of file ISCalculationCorrelated.h.

     {
       return fStepSize;
     }

Member Data Documentation

double larg4::ISCalculationCorrelated::fEfield

private

value of electric field from LArProperties service

Definition at line 44 of file ISCalculationCorrelated.h.

double larg4::ISCalculationCorrelated::fLarqlAlpha

private

from LArG4Parameters service

Definition at line 56 of file ISCalculationCorrelated.h.

double larg4::ISCalculationCorrelated::fLarqlBeta

private

from LArG4Parameters service

Definition at line 57 of file ISCalculationCorrelated.h.

double larg4::ISCalculationCorrelated::fLarqlChi0A

private

from LArG4Parameters service

Definition at line 52 of file ISCalculationCorrelated.h.

double larg4::ISCalculationCorrelated::fLarqlChi0B

private

from LArG4Parameters service

Definition at line 53 of file ISCalculationCorrelated.h.

double larg4::ISCalculationCorrelated::fLarqlChi0C

private

from LArG4Parameters service

Definition at line 54 of file ISCalculationCorrelated.h.

double larg4::ISCalculationCorrelated::fLarqlChi0D

private

from LArG4Parameters service

Definition at line 55 of file ISCalculationCorrelated.h.

double larg4::ISCalculationCorrelated::fModBoxA

private

from LArG4Parameters service

Definition at line 50 of file ISCalculationCorrelated.h.

double larg4::ISCalculationCorrelated::fModBoxB

private

from LArG4Parameters service

Definition at line 51 of file ISCalculationCorrelated.h.

double larg4::ISCalculationCorrelated::fRecombA

private

from LArG4Parameters service

Definition at line 48 of file ISCalculationCorrelated.h.

double larg4::ISCalculationCorrelated::fRecombk

private

from LArG4Parameters service

Definition at line 49 of file ISCalculationCorrelated.h.

double larg4::ISCalculationCorrelated::fScintPreScale

private

scintillation pre-scale from LArProperties service

Definition at line 47 of file ISCalculationCorrelated.h.

double larg4::ISCalculationCorrelated::fStepSize

private

maximum step to take

Definition at line 43 of file ISCalculationCorrelated.h.

bool larg4::ISCalculationCorrelated::fUseModBoxRecomb

private

from LArG4Parameters service

Definition at line 58 of file ISCalculationCorrelated.h.

bool larg4::ISCalculationCorrelated::fUseModLarqlRecomb

private

from LArG4Parameters service

Definition at line 59 of file ISCalculationCorrelated.h.

double larg4::ISCalculationCorrelated::fWion

private

W_ion (23.6 eV) == 1/fGeVToElectrons.

Definition at line 45 of file ISCalculationCorrelated.h.

double larg4::ISCalculationCorrelated::fWph

private

W_ph (19.5 eV)

Definition at line 46 of file ISCalculationCorrelated.h.

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

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation