larg4::ISCalculationSeparate Class Reference

#include "ISCalculationSeparate.h"

Inheritance diagram for larg4::ISCalculationSeparate:

Public Member Functions
	ISCalculationSeparate (CLHEP::HepRandomEngine &)
virtual	~ISCalculationSeparate ()
void	Initialize ()
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
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...
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...
bool	fScintByParticleType
	from LArProperties service More...
double	fScintYieldFactor
	scintillation yield factor More...
G4EmSaturation *	fEMSaturation
	pointer to EM saturation More...

Detailed Description

Definition at line 23 of file ISCalculationSeparate.h.

Constructor & Destructor Documentation

larg4::ISCalculationSeparate::ISCalculationSeparate

(

CLHEP::HepRandomEngine &

)

Definition at line 27 of file ISCalculationSeparate.cxx.

  {
  }

larg4::ISCalculationSeparate::~ISCalculationSeparate ( )

virtual

Definition at line 32 of file ISCalculationSeparate.cxx.

  {
  }

Member Function Documentation

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

virtual

Implements larg4::ISCalculation.

Definition at line 86 of file ISCalculationSeparate.cxx.

References larg4::ISCalculation::EFieldAtStep(), fEfield, fEMSaturation, larg4::ISCalculation::fEnergyDeposit, fGeVToElectrons, fModBoxA, fModBoxB, larg4::ISCalculation::fNumIonElectrons, larg4::ISCalculation::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;
  }

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

inherited

Definition at line 30 of file ISCalculation.cxx.

Referenced by CalculateIonizationAndScintillation(), and 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::ISCalculation::EnergyDeposit ( ) const

inlineinherited

Definition at line 29 of file ISCalculation.h.

References larg4::ISCalculation::fEnergyDeposit.

Referenced by larg4::IonizationAndScintillation::Reset().

{ return fEnergyDeposit;   }

void larg4::ISCalculationSeparate::Initialize ( )

virtual

Implements larg4::ISCalculation.

Definition at line 37 of file ISCalculationSeparate.cxx.

References detinfo::DetectorProperties::Density(), fEfield, fEMSaturation, fGeVToElectrons, fModBoxA, fModBoxB, fRecombA, fRecombk, fScintByParticleType, fScintYieldFactor, fStepSize, fUseModBoxRecomb, sim::LArG4Parameters::GeVToElectrons(), max, sim::LArG4Parameters::ModBoxA(), sim::LArG4Parameters::ModBoxB(), sim::LArG4Parameters::RecombA(), sim::LArG4Parameters::Recombk(), 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;
  }

double larg4::ISCalculation::NumberIonizationElectrons ( ) const

inlineinherited

Definition at line 30 of file ISCalculation.h.

References larg4::ISCalculation::fNumIonElectrons.

Referenced by larg4::IonizationAndScintillation::Reset().

{ return fNumIonElectrons; }

double larg4::ISCalculation::NumberScintillationPhotons ( ) const

inlineinherited

Definition at line 31 of file ISCalculation.h.

References larg4::ISCalculation::fNumScintPhotons.

Referenced by larg4::IonizationAndScintillation::Reset().

{ return fNumScintPhotons; }

void larg4::ISCalculationSeparate::Reset ( )

virtual

Implements larg4::ISCalculation.

Definition at line 75 of file ISCalculationSeparate.cxx.

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

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

    return;
  }

double larg4::ISCalculationSeparate::StepSizeLimit ( ) const

inlinevirtual

Implements larg4::ISCalculation.

Definition at line 33 of file ISCalculationSeparate.h.

{ return fStepSize;            }

double larg4::ISCalculation::VisibleEnergyDeposit ( ) const

inlineinherited

Definition at line 32 of file ISCalculation.h.

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

{ return fVisibleEnergyDeposition; }

Member Data Documentation

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

G4EmSaturation* larg4::ISCalculationSeparate::fEMSaturation

private

pointer to EM saturation

Definition at line 47 of file ISCalculationSeparate.h.

Referenced by CalculateIonizationAndScintillation(), and Initialize().

double larg4::ISCalculation::fEnergyDeposit

protectedinherited

total energy deposited in the step

Definition at line 39 of file ISCalculation.h.

Referenced by larg4::ISCalculationNEST::CalculateIonizationAndScintillation(), CalculateIonizationAndScintillation(), larg4::ISCalculation::EnergyDeposit(), larg4::ISCalculationNEST::Reset(), 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().

double larg4::ISCalculationSeparate::fModBoxA

private

from LArG4Parameters service

Definition at line 42 of file ISCalculationSeparate.h.

Referenced by CalculateIonizationAndScintillation(), and Initialize().

double larg4::ISCalculationSeparate::fModBoxB

private

from LArG4Parameters service

Definition at line 43 of file ISCalculationSeparate.h.

Referenced by CalculateIonizationAndScintillation(), and Initialize().

double larg4::ISCalculation::fNumIonElectrons

protectedinherited

number of ionization electrons for this step

Definition at line 40 of file ISCalculation.h.

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

double larg4::ISCalculation::fNumScintPhotons

protectedinherited

number of scintillation photons for this step

Definition at line 41 of file ISCalculation.h.

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

double larg4::ISCalculationSeparate::fRecombA

private

from LArG4Parameters service

Definition at line 40 of file ISCalculationSeparate.h.

Referenced by CalculateIonizationAndScintillation(), and Initialize().

double larg4::ISCalculationSeparate::fRecombk

private

from LArG4Parameters service

Definition at line 41 of file ISCalculationSeparate.h.

Referenced by CalculateIonizationAndScintillation(), 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 46 of file ISCalculationSeparate.h.

Referenced by CalculateIonizationAndScintillation(), 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 44 of file ISCalculationSeparate.h.

Referenced by 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:

• /cvmfs/larsoft.opensciencegrid.org/products/larsim/v07_09_00/source/larsim/LArG4/ISCalculationSeparate.h
• /cvmfs/larsoft.opensciencegrid.org/products/larsim/v07_09_00/source/larsim/LArG4/ISCalculationSeparate.cxx