ISCalculationCorrelated.h

Go to the documentation of this file.

// \file  ISCalculationCorrelated.h
// \brief Interface to algorithm class for a specific calculation of
//        ionization electrons and scintillation photons, based on
//        simple microphysics arguments to establish an anticorrelation
//        between these two quantities.
//
//        To enable this in simulation, change LArG4Parameters variable
//        in your fhicl file:
//
//        services.LArG4Parameters.IonAndScintCalculator: "Correlated"
//
// \author wforeman @ iit.edu
// Modified: Adding corrections for low electric field (LArQL model)
// Mar 2021 by L. Paulucci and F. Marinho
#ifndef LARG4_ISCALCULATIONCORRELATED_H
#define LARG4_ISCALCULATIONCORRELATED_H

#include "larsim/LegacyLArG4/ISCalculation.h"
namespace detinfo {
  class DetectorPropertiesData;
}

// forward declarations
class G4Step;

namespace larg4 {

  class ISCalculationCorrelated : public ISCalculation {
  public:
    explicit ISCalculationCorrelated(detinfo::DetectorPropertiesData const& detProp);

    void Reset();
    void CalculateIonizationAndScintillation(const G4Step* step);
    double StepSizeLimit() const { return fStepSize; }

  private:
    double fStepSize;            
    double fEfield;              
    double fWion;                
    double fWph;                 
    double fScintPreScale;       
    double fRecombA;             
    double fRecombk;             
    double fModBoxA;             
    double fModBoxB;             
    double fLarqlChi0A;          
    double fLarqlChi0B;          
    double fLarqlChi0C;          
    double fLarqlChi0D;          
    double fLarqlAlpha;          
    double fLarqlBeta;           
    double fEllipsModBoxA;       
    double fEllipsModBoxB;       
    double fEllipsModBoxR;       
    bool fUseEllipsModBoxRecomb; 
    bool fUseModBoxRecomb;       
    bool fUseModLarqlRecomb;     

    double EscapingEFraction(
      double const dEdx); //LArQL chi0 function = fraction of escaping electrons
    double FieldCorrection(
      double const EF,
      double const dEdx); //LArQL f_corr function = correction factor for electric field dependence
  };
}
#endif // LARG4_ISCALCULATIONCORRELATED_H