#include "ISCalculationSeparate.h"

Inheritance diagram for larg4::ISCalculationSeparate:

Public Member Functions
	ISCalculationSeparate ()

virtual	~ISCalculationSeparate ()

void	Initialize (const detinfo::LArProperties larp, const detinfo::DetectorProperties detp, const sim::LArG4Parameters lgp, const spacecharge::SpaceCharge sce)

void	Reset ()

void	CalculateIonizationAndScintillation (sim::SimEnergyDeposit const &edep)

double	EnergyDeposit () const

double	NumberIonizationElectrons () const

double	NumberScintillationPhotons () const

double	EFieldAtStep (double efield, sim::SimEnergyDeposit const &edep)

	ISCalculationSeparate (CLHEP::HepRandomEngine &)

virtual	~ISCalculationSeparate ()

void	Initialize ()

void	Reset ()

void	CalculateIonizationAndScintillation (const G4Step *step)

double	StepSizeLimit () const

double	VisibleEnergyDeposit () const

double	EFieldAtStep (double fEfield, const G4Step *step) const

Protected Attributes
double	fVisibleEnergyDeposition

Private Member Functions
double	EFieldAtStep (double efield, float x, float y, float z)

void	CalculateIonization (float e, float ds, float x, float y, float z)

void	CalculateIonization (sim::SimEnergyDeposit const &edep)

void	CalculateScintillation (float e, int pdg)

void	CalculateScintillation (sim::SimEnergyDeposit const &edep)

Private Attributes
double	fRecombA
	from LArG4Parameters service More...

double	fRecombk
	from LArG4Parameters service More...

double	fModBoxA
	from LArG4Parameters service More...

double	fModBoxB
	from LArG4Parameters service More...

bool	fUseModBoxRecomb
	from LArG4Parameters service More...

double	fScintYieldFactor
	scintillation yield factor More...

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...

geo::Vector_t	fEfieldOffsets

const detinfo::LArProperties *	fLArProp

const spacecharge::SpaceCharge *	fSCE

const detinfo::DetectorProperties *	fDetProp

const sim::LArG4Parameters *	fLArG4Prop

double	fStepSize
	maximum step to take More...

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

double	fGeVToElectrons
	conversion factor from LArProperties service More...

bool	fScintByParticleType
	from LArProperties service More...

G4EmSaturation *	fEMSaturation
	pointer to EM saturation More...

Detailed Description

Definition at line 26 of file ISCalculationSeparate.h.

Constructor & Destructor Documentation

larg4::ISCalculationSeparate::ISCalculationSeparate ( )

Definition at line 28 of file ISCalculationSeparate.cxx.

29 {

30 }

larg4::ISCalculationSeparate::~ISCalculationSeparate ( )

virtual

Definition at line 33 of file ISCalculationSeparate.cxx.

Referenced by ISCalculationSeparate().

34 {

35 }

larg4::ISCalculationSeparate::ISCalculationSeparate ( CLHEP::HepRandomEngine & )

Definition at line 27 of file ISCalculationSeparate.cxx.

References ~ISCalculationSeparate().

28 {

29 }

virtual larg4::ISCalculationSeparate::~ISCalculationSeparate ( )

virtual

Member Function Documentation

void larg4::ISCalculationSeparate::CalculateIonization	(	float	e,
		float	ds,
		float	x,
		float	y,
		float	z
	)

private

Definition at line 78 of file ISCalculationSeparate.cxx.

References detinfo::DetectorProperties::Efield(), EFieldAtStep(), fDetProp, fEnergyDeposit, fLArG4Prop, fModBoxA, fModBoxB, fNumIonElectrons, fRecombk, fUseModBoxRecomb, sim::LArG4Parameters::GeVToElectrons(), LOG_DEBUG, sim::LArG4Parameters::RecombA(), x, y, and z.

Referenced by CalculateIonization(), and CalculateIonizationAndScintillation().

                                                                             {
     double recomb = 0.;
     double dEdx   = e/ds;
     double EFieldStep = EFieldAtStep(fDetProp->Efield(),x,y,z);
     
     // Guard against spurious values of dE/dx. Note: assumes density of LAr
     if(dEdx < 1.) dEdx = 1.;
     
     if(fUseModBoxRecomb) {
       if(ds){
         double Xi = fModBoxB * dEdx / EFieldStep;
         recomb = log(fModBoxA + Xi) / Xi;
       }
       else 
         recomb = 0;
     }
     else{
       recomb = fLArG4Prop->RecombA() / (1. + dEdx * fRecombk / EFieldStep);
     }
     
     
     // 1.e-3 converts fEnergyDeposit to GeV
     fNumIonElectrons = fLArG4Prop->GeVToElectrons() * 1.e-3 * e * recomb;
 
     LOG_DEBUG("ISCalculationSeparate") 
       << " Electrons produced for " << fEnergyDeposit 
       << " MeV deposited with "     << recomb 
       << " recombination: "         << fNumIonElectrons << std::endl; 
   }

void larg4::ISCalculationSeparate::CalculateIonization ( sim::SimEnergyDeposit const & edep )

private

Definition at line 111 of file ISCalculationSeparate.cxx.

References CalculateIonization(), sim::SimEnergyDeposit::Energy(), sim::SimEnergyDeposit::MidPointX(), sim::SimEnergyDeposit::MidPointY(), sim::SimEnergyDeposit::MidPointZ(), and sim::SimEnergyDeposit::StepLength().

                                                                                 {
     CalculateIonization(edep.Energy(),edep.StepLength(),
                         edep.MidPointX(),edep.MidPointY(),edep.MidPointZ());
   }

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

virtual

Implements larg4::ISCalculation.

Definition at line 86 of file ISCalculationSeparate.cxx.

References EFieldAtStep(), fEfield, fEMSaturation, fEnergyDeposit, fGeVToElectrons, fModBoxA, fModBoxB, fNumIonElectrons, fNumScintPhotons, fRecombA, fRecombk, fScintByParticleType, fScintYieldFactor, fUseModBoxRecomb, larg4::ISCalculation::fVisibleEnergyDeposition, and LOG_DEBUG.

   {
 
     fEnergyDeposit = step->GetTotalEnergyDeposit()/CLHEP::MeV;
 
     // 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 recomb = 0.;
     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.;
 
     if(fUseModBoxRecomb) {
       if (dx){
         double Xi = fModBoxB * dEdx / EFieldStep;
         recomb = log(fModBoxA + Xi) / Xi;
       }
       else 
         recomb = 0;
     } 
     else{
       recomb = fRecombA/(1. + dEdx * fRecombk / EFieldStep);
     }
 
 
     // 1.e-3 converts fEnergyDeposit to GeV
     fNumIonElectrons = fGeVToElectrons * 1.e-3 * fEnergyDeposit * recomb;
 
     LOG_DEBUG("ISCalculationSeparate") << " Electrons produced for " << fEnergyDeposit 
                                        << " MeV deposited with "     << recomb 
                                        << " recombination: "         << fNumIonElectrons;
 
     // Now do the scintillation
     G4MaterialPropertiesTable* mpt = step->GetTrack()->GetMaterial()->GetMaterialPropertiesTable();
     if( !mpt) 
       throw cet::exception("ISCalculationSeparate") << "Cannot find materials property table"
                                                     << " for this step! "
                                                     << step->GetTrack()->GetMaterial() << "\n";
 
     // if not doing the scintillation by particle type, use the saturation
     double scintYield = mpt->GetConstProperty("SCINTILLATIONYIELD");
 
     if(fScintByParticleType){
 
       LOG_DEBUG("ISCalculationSeparate") << "scintillating by particle type";
 
       // Get the definition of the current particle
       G4ParticleDefinition *pDef = step->GetTrack()->GetDynamicParticle()->GetDefinition();
       //G4MaterialPropertyVector *Scint_Yield_Vector = NULL;
 
       // Obtain the G4MaterialPropertyVectory containing the
       // scintillation light yield as a function of the deposited
       // energy for the current particle type
       
       // Protons
       if(pDef == G4Proton::ProtonDefinition()){
         scintYield = mpt->GetConstProperty("PROTONSCINTILLATIONYIELD");
       }
       // Muons
       else if(pDef == G4MuonPlus::MuonPlusDefinition() ||
               pDef == G4MuonMinus::MuonMinusDefinition()){
         scintYield = mpt->GetConstProperty("MUONSCINTILLATIONYIELD");
       }
       // Pions
       else if(pDef == G4PionPlus::PionPlusDefinition() ||
               pDef == G4PionMinus::PionMinusDefinition()){
         scintYield = mpt->GetConstProperty("PIONSCINTILLATIONYIELD");
       }
       // Kaons
       else if(pDef == G4KaonPlus::KaonPlusDefinition() ||
               pDef == G4KaonMinus::KaonMinusDefinition()){
         scintYield = mpt->GetConstProperty("KAONSCINTILLATIONYIELD");
       }
       // Alphas
       else if(pDef == G4Alpha::AlphaDefinition()){
         scintYield = mpt->GetConstProperty("ALPHASCINTILLATIONYIELD");
       }
       // Electrons (must also account for shell-binding energy
       // attributed to gamma from standard PhotoElectricEffect)
       else if(pDef == G4Electron::ElectronDefinition() ||
               pDef == G4Gamma::GammaDefinition()){
         scintYield = mpt->GetConstProperty("ELECTRONSCINTILLATIONYIELD");
       }       
       // Default for particles not enumerated/listed above
       else{
         scintYield = mpt->GetConstProperty("ELECTRONSCINTILLATIONYIELD");
       }
 
       // If the user has not specified yields for (p,d,t,a,carbon)
       // then these unspecified particles will default to the 
       // electron's scintillation yield
            
       // Throw an exception if no scintillation yield is found
       if (!scintYield) 
         throw cet::exception("ISCalculationSeparate") << "Request for scintillation yield for energy "
                                                       << "deposit and particle type without correct "
                                                       << "entry in MaterialPropertiesTable\n"
                                                       << "ScintillationByParticleType requires at "
                                                       << "minimum that ELECTRONSCINTILLATIONYIELD is "
                                                       << "set by the user\n";
 
       fNumScintPhotons =  scintYield * fEnergyDeposit;
     }
     else if(fEMSaturation){
       // The default linear scintillation process
       //fEMSaturation->SetVerbose(1);
       fVisibleEnergyDeposition = fEMSaturation->VisibleEnergyDepositionAtAStep(step);
       fNumScintPhotons = fScintYieldFactor * scintYield * fVisibleEnergyDeposition;
       //fNumScintPhotons = fScintYieldFactor * scintYield * fEMSaturation->VisibleEnergyDepositionAtAStep(step);
       //I need a dump here
       //mf::LogInfo("EMSaturation") <<"\n\nfEMSaturation VisibleEnergyDepositionAtAStep(step): "<<fEMSaturation->VisibleEnergyDepositionAtAStep(step)<<"\n" <<"BirksCoefs: \n";
      // fEMSaturation->DumpBirksCoefficients();
       //mf::LogInfo("EMSaturation")<<"\n" <<"G4Birks: \n";
       //fEMSaturation->DumpG4BirksCoefficients();
       //mf::LogInfo("EMSaturation")<<"fScintYieldFactor: "<<fScintYieldFactor <<"\nscintYield: "<<scintYield<<"\nfEMVisAtStep: "<<fEMSaturation->VisibleEnergyDepositionAtAStep(step)<<"\nfNumScintPhotons: "<<fNumScintPhotons<<"\n";
       //mf::LogInfo("EMSaturation")<<"fTotalEnergyDeposit: "<< step->GetTotalEnergyDeposit()/CLHEP::MeV<<"\nfNumIonElectrons: "<<fNumIonElectrons<<"\n";
     }
     else{
       fNumScintPhotons = fScintYieldFactor * scintYield * fEnergyDeposit;
       fVisibleEnergyDeposition = 0.0; //This is set to zero because I have not made a correct implimentation of this value for anything but EMSaturation.
     }
 
     LOG_DEBUG("ISCalculationSeparate") << "number photons: " << fNumScintPhotons 
                                        << " energy: "        << fEnergyDeposit/CLHEP::MeV
                                        << " saturation: " 
                                        << fEMSaturation->VisibleEnergyDepositionAtAStep(step)
                                        << " step length: "   << step->GetStepLength()/CLHEP::cm;
 
 
     return;
   }

void larg4::ISCalculationSeparate::CalculateIonizationAndScintillation ( sim::SimEnergyDeposit const & edep )

Definition at line 168 of file ISCalculationSeparate.cxx.

References CalculateIonization(), CalculateScintillation(), sim::SimEnergyDeposit::Energy(), and fEnergyDeposit.

Referenced by larg4::ISCalculationAna::analyze(), spacecharge::ShiftEdepSCE::produce(), phot::PhotonLibraryPropagation::produce(), and detsim::SimDriftElectrons::produce().

   {
     fEnergyDeposit = edep.Energy();
     CalculateIonization(edep);
     CalculateScintillation(edep);
   }

void larg4::ISCalculationSeparate::CalculateScintillation	(	float	e,
		int	pdg
	)

private

Definition at line 117 of file ISCalculationSeparate.cxx.

References detinfo::LArProperties::AlphaScintYield(), e, detinfo::LArProperties::ElectronScintYield(), fLArProp, fNumScintPhotons, fScintYieldFactor, detinfo::LArProperties::KaonScintYield(), LOG_DEBUG, detinfo::LArProperties::MuonScintYield(), detinfo::LArProperties::PionScintYield(), detinfo::LArProperties::ProtonScintYield(), detinfo::LArProperties::ScintByParticleType(), and detinfo::LArProperties::ScintYield().

Referenced by CalculateIonizationAndScintillation(), and CalculateScintillation().

   {
     double scintYield = fLArProp->ScintYield(true);
     if(fLArProp->ScintByParticleType()){
 
       LOG_DEBUG("ISCalculationSeparate") << "scintillating by particle type";
 
       switch(pdg) {
 
       case 2212:
         scintYield = fLArProp->ProtonScintYield(true);
         break;
       case 13:
       case -13:
         scintYield = fLArProp->MuonScintYield(true);
         break;
       case 211:
       case -211:
         scintYield = fLArProp->PionScintYield(true);
         break;
       case 321:
       case -321:
         scintYield = fLArProp->KaonScintYield(true);
         break;
       case 1000020040:
         scintYield = fLArProp->AlphaScintYield(true);
         break;
       case 11:
       case -11:
       case 22:
         scintYield = fLArProp->ElectronScintYield(true);
         break;
       default:
         scintYield = fLArProp->ElectronScintYield(true);
 
       }
 
       fNumScintPhotons =  scintYield * e;
     }
     else
       fNumScintPhotons = fScintYieldFactor * scintYield * e;
     
   }

void larg4::ISCalculationSeparate::CalculateScintillation ( sim::SimEnergyDeposit const & edep )

private

Definition at line 162 of file ISCalculationSeparate.cxx.

References CalculateScintillation(), sim::SimEnergyDeposit::Energy(), and sim::SimEnergyDeposit::PdgCode().

   {
     CalculateScintillation(edep.Energy(),edep.PdgCode());
   }

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

inherited

Definition at line 30 of file ISCalculation.cxx.

Referenced by larg4::ISCalculation::VisibleEnergyDeposit().

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

double larg4::ISCalculationSeparate::EFieldAtStep	(	double	efield,
		sim::SimEnergyDeposit const &	edep
	)

Definition at line 175 of file ISCalculationSeparate.cxx.

References sim::SimEnergyDeposit::MidPointX(), sim::SimEnergyDeposit::MidPointY(), and sim::SimEnergyDeposit::MidPointZ().

Referenced by CalculateIonization(), and CalculateIonizationAndScintillation().

   {
     return EFieldAtStep(efield,
                         edep.MidPointX(),edep.MidPointY(),edep.MidPointZ());
   }

double larg4::ISCalculationSeparate::EFieldAtStep	(	double	efield,
		float	x,
		float	y,
		float	z
	)

private

Definition at line 181 of file ISCalculationSeparate.cxx.

References spacecharge::SpaceCharge::EnableSimEfieldSCE(), fEfieldOffsets, fSCE, spacecharge::SpaceCharge::GetEfieldOffsets(), x, and y.

   {
     double EField = efield;
     if (fSCE->EnableSimEfieldSCE())
       {
         fEfieldOffsets = fSCE->GetEfieldOffsets(geo::Point_t{x,y,z});
         EField = std::sqrt( (efield + efield*fEfieldOffsets.X())*(efield + efield*fEfieldOffsets.X()) +
                             (efield*fEfieldOffsets.Y()+efield*fEfieldOffsets.Y()) +
                             (efield*fEfieldOffsets.Z()+efield*fEfieldOffsets.Z()) );
       }
     return EField;
   }

double larg4::ISCalculationSeparate::EnergyDeposit ( ) const

inline

Definition at line 39 of file ISCalculationSeparate.h.

Referenced by larg4::ISCalculationAna::analyze().

39 { return fEnergyDeposit; }

larg4::ISCalculationSeparate::fEnergyDeposit

double fEnergyDeposit

total energy deposited in the step

Definition: ISCalculationSeparate.h:57

void larg4::ISCalculationSeparate::Initialize ( )

virtual

Implements larg4::ISCalculation.

Definition at line 37 of file ISCalculationSeparate.cxx.

References detinfo::DetectorProperties::Density(), fEfield, fEMSaturation, fEnergyDeposit, fGeVToElectrons, fModBoxA, fModBoxB, fNumIonElectrons, fNumScintPhotons, fRecombA, fRecombk, fScintByParticleType, fScintYieldFactor, fStepSize, fUseModBoxRecomb, sim::LArG4Parameters::GeVToElectrons(), max, sim::LArG4Parameters::ModBoxA(), sim::LArG4Parameters::ModBoxB(), sim::LArG4Parameters::RecombA(), sim::LArG4Parameters::Recombk(), Reset(), detinfo::LArProperties::ScintByParticleType(), sim::LArG4Parameters::UseModBoxRecomb(), sim::LArVoxelCalculator::VoxelSizeX(), sim::LArVoxelCalculator::VoxelSizeY(), and sim::LArVoxelCalculator::VoxelSizeZ().

   {
     art::ServiceHandle<sim::LArG4Parameters> lgpHandle;
     const detinfo::LArProperties* larp = lar::providerFrom<detinfo::LArPropertiesService>();
     const detinfo::DetectorProperties* detprop = lar::providerFrom<detinfo::DetectorPropertiesService>();
 
     double density       = detprop->Density(detprop->Temperature());
     fEfield       = detprop->Efield();
     fScintByParticleType = larp->ScintByParticleType();
     fGeVToElectrons      = lgpHandle->GeVToElectrons();
     
     // \todo get scintillation yield from LArG4Parameters or LArProperties
     fScintYieldFactor  = 1.;
 
     // 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;
     fUseModBoxRecomb     = lgpHandle->UseModBoxRecomb();  
 
     // Use Birks Correction in the Scintillation process    
     fEMSaturation = G4LossTableManager::Instance()->EmSaturation();
 
     // determine the step size using the voxel sizes
     art::ServiceHandle<sim::LArVoxelCalculator> lvc;
     double maxsize = std::max(lvc->VoxelSizeX(), std::max(lvc->VoxelSizeY(), lvc->VoxelSizeZ())) * CLHEP::cm;
 
     fStepSize = 0.1 * maxsize;
 
     return;
   }

void larg4::ISCalculationSeparate::Initialize	(	const detinfo::LArProperties *	larp,
		const detinfo::DetectorProperties *	detp,
		const sim::LArG4Parameters *	lgp,
		const spacecharge::SpaceCharge *	sce
	)

Definition at line 38 of file ISCalculationSeparate.cxx.

References detinfo::DetectorProperties::Density(), fDetProp, fLArG4Prop, fLArProp, fModBoxA, fModBoxB, fRecombA, fRecombk, fSCE, fScintYieldFactor, fUseModBoxRecomb, sim::LArG4Parameters::ModBoxA(), sim::LArG4Parameters::ModBoxB(), sim::LArG4Parameters::RecombA(), sim::LArG4Parameters::Recombk(), Reset(), detinfo::DetectorProperties::Temperature(), and sim::LArG4Parameters::UseModBoxRecomb().

Referenced by spacecharge::ShiftEdepSCE::beginJob(), larg4::ISCalculationAna::beginJob(), detsim::SimDriftElectrons::beginJob(), and phot::PhotonLibraryPropagation::produce().

   {
     fLArProp = larp;
     fSCE = sce;
     fDetProp = detp;
     fLArG4Prop = lgp;
     
     // \todo get scintillation yield from LArG4Parameters or LArProperties
     fScintYieldFactor  = 1.;
 
     // 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             = (double)lgp->RecombA();
     fRecombk             = (double)lgp->Recombk()/detp->Density(detp->Temperature());
     fModBoxA             = (double)lgp->ModBoxA();
     fModBoxB             = (double)lgp->ModBoxB()/detp->Density(detp->Temperature());
     fUseModBoxRecomb     = (bool)lgp->UseModBoxRecomb();  
 
     this->Reset();
     
     return;
   }

double larg4::ISCalculationSeparate::NumberIonizationElectrons ( ) const

inline

Definition at line 40 of file ISCalculationSeparate.h.

Referenced by larg4::ISCalculationAna::analyze(), spacecharge::ShiftEdepSCE::produce(), and detsim::SimDriftElectrons::produce().

40 { return fNumIonElectrons; }

larg4::ISCalculationSeparate::fNumIonElectrons

double fNumIonElectrons

number of ionization electrons for this step

Definition: ISCalculationSeparate.h:58

double larg4::ISCalculationSeparate::NumberScintillationPhotons ( ) const

inline

Definition at line 41 of file ISCalculationSeparate.h.

Referenced by larg4::ISCalculationAna::analyze(), spacecharge::ShiftEdepSCE::produce(), and phot::PhotonLibraryPropagation::produce().

41 { return fNumScintPhotons; }

larg4::ISCalculationSeparate::fNumScintPhotons

double fNumScintPhotons

number of scintillation photons for this step

Definition: ISCalculationSeparate.h:59

void larg4::ISCalculationSeparate::Reset ( )

virtual

Implements larg4::ISCalculation.

void larg4::ISCalculationSeparate::Reset ( )

virtual

Implements larg4::ISCalculation.

Definition at line 67 of file ISCalculationSeparate.cxx.

References fEnergyDeposit, fNumIonElectrons, and fNumScintPhotons.

Referenced by larg4::ISCalculationAna::analyze(), Initialize(), spacecharge::ShiftEdepSCE::produce(), phot::PhotonLibraryPropagation::produce(), and detsim::SimDriftElectrons::produce().

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

double larg4::ISCalculationSeparate::StepSizeLimit ( ) const

inlinevirtual

Implements larg4::ISCalculation.

Definition at line 33 of file ISCalculationSeparate.h.

33 { return fStepSize; }

larg4::ISCalculationSeparate::fStepSize

double fStepSize

maximum step to take

Definition: ISCalculationSeparate.h:37

double larg4::ISCalculation::VisibleEnergyDeposit ( ) const

inlineinherited

Definition at line 32 of file ISCalculation.h.

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

32 { return fVisibleEnergyDeposition; }

larg4::ISCalculation::fVisibleEnergyDeposition

double fVisibleEnergyDeposition

Definition: ISCalculation.h:42

Member Data Documentation

const detinfo::DetectorProperties* larg4::ISCalculationSeparate::fDetProp

private

Definition at line 66 of file ISCalculationSeparate.h.

Referenced by CalculateIonization(), and Initialize().

double larg4::ISCalculationSeparate::fEfield

private

value of electric field from LArProperties service

Definition at line 38 of file ISCalculationSeparate.h.

Referenced by CalculateIonizationAndScintillation(), and Initialize().

geo::Vector_t larg4::ISCalculationSeparate::fEfieldOffsets

private

Definition at line 61 of file ISCalculationSeparate.h.

Referenced by EFieldAtStep().

G4EmSaturation* larg4::ISCalculationSeparate::fEMSaturation

private

pointer to EM saturation

Definition at line 47 of file ISCalculationSeparate.h.

Referenced by CalculateIonizationAndScintillation(), and Initialize().

double larg4::ISCalculationSeparate::fEnergyDeposit

private

total energy deposited in the step

Definition at line 57 of file ISCalculationSeparate.h.

Referenced by CalculateIonization(), CalculateIonizationAndScintillation(), Initialize(), and Reset().

double larg4::ISCalculationSeparate::fGeVToElectrons

private

conversion factor from LArProperties service

Definition at line 39 of file ISCalculationSeparate.h.

Referenced by CalculateIonizationAndScintillation(), and Initialize().

const sim::LArG4Parameters* larg4::ISCalculationSeparate::fLArG4Prop

private

Definition at line 67 of file ISCalculationSeparate.h.

Referenced by CalculateIonization(), and Initialize().

const detinfo::LArProperties* larg4::ISCalculationSeparate::fLArProp

private

Definition at line 64 of file ISCalculationSeparate.h.

Referenced by CalculateScintillation(), and Initialize().

double larg4::ISCalculationSeparate::fModBoxA

private

from LArG4Parameters service

Definition at line 51 of file ISCalculationSeparate.h.

Referenced by CalculateIonization(), CalculateIonizationAndScintillation(), and Initialize().

double larg4::ISCalculationSeparate::fModBoxB

private

from LArG4Parameters service

Definition at line 52 of file ISCalculationSeparate.h.

Referenced by CalculateIonization(), CalculateIonizationAndScintillation(), and Initialize().

double larg4::ISCalculationSeparate::fNumIonElectrons

private

number of ionization electrons for this step

Definition at line 58 of file ISCalculationSeparate.h.

Referenced by CalculateIonization(), CalculateIonizationAndScintillation(), Initialize(), and Reset().

double larg4::ISCalculationSeparate::fNumScintPhotons

private

number of scintillation photons for this step

Definition at line 59 of file ISCalculationSeparate.h.

Referenced by CalculateIonizationAndScintillation(), CalculateScintillation(), Initialize(), and Reset().

double larg4::ISCalculationSeparate::fRecombA

private

from LArG4Parameters service

Definition at line 49 of file ISCalculationSeparate.h.

Referenced by CalculateIonizationAndScintillation(), and Initialize().

double larg4::ISCalculationSeparate::fRecombk

private

from LArG4Parameters service

Definition at line 50 of file ISCalculationSeparate.h.

Referenced by CalculateIonization(), CalculateIonizationAndScintillation(), and Initialize().

const spacecharge::SpaceCharge* larg4::ISCalculationSeparate::fSCE

private

Definition at line 65 of file ISCalculationSeparate.h.

Referenced by EFieldAtStep(), and Initialize().

bool larg4::ISCalculationSeparate::fScintByParticleType

private

from LArProperties service

Definition at line 45 of file ISCalculationSeparate.h.

Referenced by CalculateIonizationAndScintillation(), and Initialize().

double larg4::ISCalculationSeparate::fScintYieldFactor

private

scintillation yield factor

Definition at line 55 of file ISCalculationSeparate.h.

Referenced by CalculateIonizationAndScintillation(), CalculateScintillation(), and Initialize().

double larg4::ISCalculationSeparate::fStepSize

private

maximum step to take

Definition at line 37 of file ISCalculationSeparate.h.

Referenced by Initialize().

bool larg4::ISCalculationSeparate::fUseModBoxRecomb

private

from LArG4Parameters service

Definition at line 53 of file ISCalculationSeparate.h.

Referenced by CalculateIonization(), CalculateIonizationAndScintillation(), and Initialize().

double larg4::ISCalculation::fVisibleEnergyDeposition

protectedinherited

Definition at line 42 of file ISCalculation.h.

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

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

larsim/v06_53_00/source/larsim/IonizationScintillation/ISCalculationSeparate.h
larsim/v06_53_00/source/larsim/IonizationScintillation/ISCalculationSeparate.cxx

Public Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation