OpDigiProperties.h

Go to the documentation of this file.

// \file OpDigiProperties.h
//
// \brief service containing information about optical digitization
//
// \author bjpjones@mit.edu
// \author kazuhiro@nevis.columbia.edu
//

#ifndef OPDET_OPDIGIPROPERTIES_H
#define OPDET_OPDIGIPROPERTIES_H

// LArSoft includes
#include "lardataobj/OpticalDetectorData/OpticalTypes.h"
#include "larcore/Geometry/Geometry.h"

// ART includes
#include "fhiclcpp/ParameterSet.h"
#include "art/Framework/Services/Registry/ActivityRegistry.h"
#include "art/Framework/Services/Registry/ServiceHandle.h"
#include "art/Framework/Services/Registry/ServiceMacros.h"

// ROOT includes
#include "TF1.h"

namespace opdet
{
    class OpDigiProperties {
    public:

      OpDigiProperties(fhicl::ParameterSet const& pset, art::ActivityRegistry& reg);
      ~OpDigiProperties();

      void reconfigure(fhicl::ParameterSet const& p);

      double SampleFreq()   { return fSampleFreq; }
      double TimeBegin()    { return fTimeBegin;  }
      double TimeEnd()      { return fTimeEnd;    }

      std::vector<double> WaveformInit(std::string WaveformFile);

      optdata::TimeSlice_t GetTimeSlice(double time_ns);

      double QE()              const { return fQE;              }
      double DarkRate()        const { return fDarkRate;        }
      double PedFlucRate()     const { return fPedFlucRate;     }
      optdata::ADC_Count_t PedFlucAmp()      const { return fPedFlucAmp;      }
      optdata::ADC_Count_t SaturationScale() const { return fSaturationScale; }
      optdata::ADC_Count_t ADCBaseline() const { return fADCBaseline; }
      double WFTimeConstant()  const { return fWFTimeConstant;  }
      double WFPowerFactor()   const { return fWFPowerFactor;   }
      double LowGainMean  () const { return fLowGainMean;  }
      double HighGainMean () const { return fHighGainMean; }
      double LowGainMean  (optdata::Channel_t ch) const;
      double HighGainMean (optdata::Channel_t ch) const;
      
      double GainSpread        () const { return fGainSpread;         }
      double GainSpread_PMT2PMT() const { return fGainSpread_PMT2PMT; }
      double LowGain(optdata::Channel_t ch) const;
      double HighGain(optdata::Channel_t ch) const;

      std::vector<double> SinglePEWaveform() const { return fWaveform;        } 
      std::vector<double> HighGainArray()    const { return fHighGainArray;   }
      std::vector<double> LowGainArray()     const { return fLowGainArray;    }
      std::vector<double> GainSpreadArray()  const { return fGainSpreadArray; }
      std::vector<optdata::ADC_Count_t> PedMeanArray() const {return fPedMeanArray;}

      double GetSPEArea();
      double GetSPECumulativeArea();
      double GetSPEAmplitude();
      double GetSPECumulativeAmplitude();
      
    private:

      double fSampleFreq;
      double fTimeBegin;
      double fTimeEnd;
      double fPERescale;

      std::vector<double> GenEmpiricalWF(std::string WaveformFile);
      std::vector<double> GenAnalyticalWF();
      void GenerateWaveform();
      void FillGainArray();
      void FillPedMeanArray();

      bool   fUseEmpiricalGain;
      bool   fUseEmpiricalShape;
      double fQE;
      double fPedFlucRate;
      optdata::ADC_Count_t fPedFlucAmp;
      double fDarkRate;
      optdata::ADC_Count_t fSaturationScale;
      optdata::ADC_Count_t fADCBaseline;
      double fADCBaseSpread;
      double fLowGainMean;
      double fHighGainMean;
      double fVoltageAmpForSPE;
      double fWFTimeConstant;
      double fWFPowerFactor;
      double fWFLength;
      double fGainSpread;
      double fGainSpread_PMT2PMT;
      TF1* fAnalyticalSPE;

      std::string fHighGainFile;
      std::string fLowGainFile;
      std::string fWaveformFile;
      std::string fGainSpreadFile;
      std::vector<double> fWaveform;
      std::vector<double> fLowGainArray;
      std::vector<double> fHighGainArray;
      std::vector<double> fGainSpreadArray;
      std::vector<optdata::ADC_Count_t> fPedMeanArray;
      art::ServiceHandle<geo::Geometry> fGeometry;
 
    }; // class OpDigiProperties

} //namespace opdet

DECLARE_ART_SERVICE(opdet::OpDigiProperties, LEGACY)

#endif