ISCalcSeparate.h

Go to the documentation of this file.

#ifndef IS_ISCALCSEPARATE_H
#define IS_ISCALCSEPARATE_H

#include <vector>
#include "larcoreobj/SimpleTypesAndConstants/geo_vectors.h"

// forward declaration
namespace detinfo { class LArProperties; class DetectorProperties; }
namespace sim { class SimEnergyDeposit; class LArG4Parameters; }
namespace spacecharge { class SpaceCharge; }

namespace larg4 {

 class ISCalcSeparate{

 public:

   ISCalcSeparate();
   virtual ~ISCalcSeparate();

   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 { return fEnergyDeposit;       }
   double NumberIonizationElectrons()  const { return fNumIonElectrons;     }
   double NumberScintillationPhotons() const { return fNumScintPhotons;     }

   double EFieldAtStep(double efield, sim::SimEnergyDeposit const& edep); //value of field with any corrections for this step  

 private:

   double                fRecombA;             
   double                fRecombk;             
   double                fModBoxA;             
   double                fModBoxB;             
   bool                  fUseModBoxRecomb;     
   double                fGeVToElectrons;      

   double                fScintYieldFactor;    

   double fEnergyDeposit;   
   double fNumIonElectrons; 
   double fNumScintPhotons; 
   
   geo::Vector_t fEfieldOffsets;
   double EFieldAtStep(double efield, float x, float y, float z);

   const detinfo::LArProperties*      fLArProp;   
   const spacecharge::SpaceCharge*    fSCE;
   const detinfo::DetectorProperties* fDetProp;
   const sim::LArG4Parameters*        fLArG4Prop;
   
   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);
   
 };
}
#endif // LARG4_ISCALCULATIONSEPARATE_H