IonizationAndScintillation.h

Go to the documentation of this file.

#ifndef LARG4_IONIZATIONANDSCINTILLATION_H
#define LARG4_IONIZATIONANDSCINTILLATION_H

#include <memory>
#include <string>

#include "larsim/LegacyLArG4/ISCalculation.h"

class G4Step;
class TH1F;
class TH2F;

namespace CLHEP {
  class HepRandomEngine;
}

namespace detinfo {
  class DetectorPropertiesData;
}

namespace larg4 {

  // The Ionization and Scintillation singleton
  class IonizationAndScintillation {
  public:
    static IonizationAndScintillation* CreateInstance(
      detinfo::DetectorPropertiesData const& detProp,
      CLHEP::HepRandomEngine& engine);
    static IonizationAndScintillation* Instance();

    // Method to reset the internal variables held in the ISCalculation
    // This method should be called at the start of any G4Step
    void Reset(const G4Step* step);

    double EnergyDeposit() const { return fISCalc->EnergyDeposit(); }
    double VisibleEnergyDeposit() const { return fISCalc->VisibleEnergyDeposit(); }
    double NumberIonizationElectrons() const { return fISCalc->NumberIonizationElectrons(); }
    double NumberScintillationPhotons() const { return fISCalc->NumberScintillationPhotons(); }
    double StepSizeLimit() const { return fISCalc->StepSizeLimit(); }

  private:
    IonizationAndScintillation(detinfo::DetectorPropertiesData const& detProp,
                               CLHEP::HepRandomEngine& engine);

    std::unique_ptr<larg4::ISCalculation>
      fISCalc;                    
    std::string fISCalculator;    
    G4Step const* fStep{nullptr}; 
    int fStepNumber{-1};          
    int fTrkID{-1};               

    TH1F* fElectronsPerStep{nullptr};   
    TH1F* fStepSize{nullptr};           
    TH1F* fPhotonsPerStep{nullptr};     
    TH1F* fEnergyPerStep{nullptr};      
    TH1F* fElectronsPerLength;          
    TH1F* fPhotonsPerLength;            
    TH1F* fElectronsPerEDep;            
    TH1F* fPhotonsPerEDep;              
    TH2F* fElectronsVsPhotons{nullptr}; 
    CLHEP::HepRandomEngine& fEngine;    
  };

} // namespace larg4

#endif // LARG4_IONIZATIONANDSCINTILLATION