opdet::OpDigiProperties Class Reference

#include <OpDigiProperties.h>

Public Member Functions
	OpDigiProperties (fhicl::ParameterSet const &pset)
double	SampleFreq () const noexcept
	Returns sample frequency in MHz. More...
double	TimeBegin () const noexcept
	Returns window start time in us ... with respect to MC photon T0. More...
double	TimeEnd () const noexcept
	Returns window end time in us ... with respect to MC photon T0. More...
std::vector< double >	WaveformInit (std::string WaveformFile)
optdata::TimeSlice_t	GetTimeSlice (double time_ns)
double	QE () const noexcept
	Returns quantum efficiency. More...
double	DarkRate () const noexcept
	Returns rate of dark noise. More...
double	PedFlucRate () const noexcept
	Returns rate of pedestal fluctuation. More...
optdata::ADC_Count_t	PedFlucAmp () const noexcept
	Returns amplitude of pedestal fluctuation. More...
optdata::ADC_Count_t	SaturationScale () const noexcept
	Returns the saturation scale of the electronics. More...
optdata::ADC_Count_t	ADCBaseline () const noexcept
	Returns the ADCBaseline set mean value. More...
double	WFTimeConstant () const noexcept
	Returns WF time constant used in analytical model. More...
double	WFPowerFactor () const noexcept
	Returns WF power factor used in analytical model. More...
double	LowGainMean () const noexcept
	Returns set mean gain value for LOW gain. More...
double	HighGainMean () const noexcept
	Returns set mean gain value for HIGH gain. More...
double	LowGainMean (optdata::Channel_t ch) const
	Returns generated LOW gain value for input channel (PMT-to-PMT spread applied) More...
double	HighGainMean (optdata::Channel_t ch) const
	Returns generated HIGH gain value for input channel (PMT-to-PMT spread applied) More...
double	GainSpread () const noexcept
	Returns set value for intrinsic gain spread (common to HIGH and LOW) More...
double	GainSpread_PMT2PMT () const noexcept
	Returns set value for PMT-to-PMT gain spread. More...
double	LowGain (optdata::Channel_t ch) const
	Generate & return LOW gain value for an input channel using mean & spread for this channel. More...
double	HighGain (optdata::Channel_t ch) const
	Generate & return HIGH gain value for an input channel using mean & spread for this channel. More...
std::vector< double > const &	SinglePEWaveform () const noexcept
	Returns a vector of double which represents a binned SPE waveform. More...
std::vector< double > const &	HighGainArray () const noexcept
	Returns an array of HIGH gain. More...
std::vector< double > const &	LowGainArray () const noexcept
	Returns an array of LOW gain. More...
std::vector< double > const &	GainSpreadArray () const noexcept
	Returns an array of gain spread. More...
std::vector< optdata::ADC_Count_t > const &	PedMeanArray () const noexcept
	Returns an array of generated pedestal mean value per channel. More...
double	GetSPEArea ()
	Utility function ... To be verified (Kazu 08/05/13) More...
double	GetSPECumulativeArea ()
	Utility function ... To be verified (Kazu 08/05/13) More...
double	GetSPEAmplitude ()
	Utility function ... To be verified (Kazu 08/05/13) More...
double	GetSPECumulativeAmplitude ()
	Utility function ... To be verified (Kazu 08/05/13) More...

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

Private Attributes
double	fSampleFreq
double	fTimeBegin
double	fTimeEnd
double	fPERescale
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
bool	fChargeNormalized
std::vector< double >	fLowGainArray
std::vector< double >	fHighGainArray
std::vector< double >	fGainSpreadArray
std::vector< optdata::ADC_Count_t >	fPedMeanArray
art::ServiceHandle< geo::Geometry const >	fGeometry

Detailed Description

Definition at line 27 of file OpDigiProperties.h.

Constructor & Destructor Documentation

opdet::OpDigiProperties::OpDigiProperties

(

fhicl::ParameterSet const &

pset

)

Definition at line 25 of file OpDigiProperties_service.cc.

    : fAnalyticalSPE(0)
  {
    fSampleFreq        = p.get< double       >("SampleFreq"       );
    fTimeBegin         = p.get< double       >("TimeBegin"        );
    fTimeEnd           = p.get< double       >("TimeEnd"          );
    fWaveformFile      = p.get< std::string  >("WaveformFile"     );
    fPERescale         = p.get< double       >("PERescale"        );

    // PMT properties
    fQE              = p.get<double>("QE");
    fDarkRate        = p.get<double>("DarkRate");
    fSaturationScale = p.get<optdata::ADC_Count_t>("SaturationScale");

    // Shaper properties
    fUseEmpiricalGain= p.get<bool        >("UseEmpiricalGain");
    fHighGainFile    = p.get<std::string >("HighGainFile");
    fLowGainFile     = p.get<std::string >("LowGainFile");
    fGainSpreadFile  = p.get<std::string >("GainSpreadFile");

    fHighGainMean       = p.get<double   >("HighGainMean");
    fLowGainMean        = p.get<double   >("LowGainMean");
    fGainSpread          = p.get<double   >("GainSpread");
    fGainSpread_PMT2PMT  = p.get<double   >("GainSpread_PMT2PMT");

    // Digitization ped fluc
    fPedFlucRate     = p.get<double>("PedFlucRate");
    fPedFlucAmp      = p.get<optdata::ADC_Count_t>("PedFlucAmp");
    fADCBaseline     = p.get<optdata::ADC_Count_t>("ADCBaseline");
    fADCBaseSpread   = p.get<double>("ADCBaseSpread");

    // WF related stuff
    fUseEmpiricalShape = p.get< bool         >("UseEmpiricalShape");
    fWFLength          = p.get< double       >("WFLength"         );

    fWaveformFile      = p.get< std::string  >("WaveformFile"     );
    fChargeNormalized      = p.get< bool  >("WaveformChargeNormalized", false);

    // Option 2: WF from analytical function
    fWFPowerFactor     = p.get< double       >("WFPowerFactor"    );
    fWFTimeConstant    = p.get< double       >("WFTimeConstant"   );
    fVoltageAmpForSPE  = p.get< double       >("VoltageAmpForSPE" );

    // Generate the SPE waveform (i.e. fWaveform)
    GenerateWaveform();

    // Fill gain array
    FillGainArray();

    // Fill baseline mean
    FillPedMeanArray();

    // Report
    std::string msg(Form("%-10s ... %-10s ... %-10s ... %-10s\n","Ch. Number","Pedestal","High Gain","Low Gain"));
    for(unsigned int i=0;i<fGeometry->NOpChannels();++i) {
      msg+=Form("%-10d ... %-10d ... %-10g ... %-10g\n",i,fPedMeanArray[i],fHighGainArray[i],fLowGainArray[i]);
    }
    mf::LogInfo(__FUNCTION__)<<msg.c_str();
  }

Member Function Documentation

optdata::ADC_Count_t opdet::OpDigiProperties::ADCBaseline ( ) const

inlinenoexcept

Returns the ADCBaseline set mean value.

Definition at line 62 of file OpDigiProperties.h.

{ return fADCBaseline; }

double opdet::OpDigiProperties::DarkRate ( ) const

inlinenoexcept

Returns rate of dark noise.

Definition at line 54 of file OpDigiProperties.h.

{ return fDarkRate;        }

void opdet::OpDigiProperties::FillGainArray ( )

private

Definition at line 193 of file OpDigiProperties_service.cc.

                                      {
    if(fUseEmpiricalGain)
      {
        // Fill fron user's text files.
        mf::LogWarning("OpDigiProperties")<<"Using empirical table of gain for each PMT...";
        std::string FullPath;
        cet::search_path sp("FW_SEARCH_PATH");

        if( !sp.find_file(fHighGainFile, FullPath) )
          throw cet::exception("OpDigiProperties") << "Unable to find high gain spread file in " << sp.to_string() << "\n";

        mf::LogWarning("OpDigiProperties")<<"OpDigiProperties opening high gain spread file at " << FullPath.c_str();
        std::ifstream HighGainFile(FullPath.c_str());
        if(HighGainFile.is_open()) {
          std::string line;
          while ( HighGainFile.good() ){
            getline(HighGainFile, line);
            fHighGainArray.push_back(strtod(line.c_str(),NULL));
          }
        }else throw cet::exception("OpDigiProperties")<<"Unable to open file!\n";

        FullPath="";
        if( !sp.find_file(fLowGainFile, FullPath) )
          throw cet::exception("OpDigiProperties") << "Unable to find low gain spread file in " << sp.to_string() << "\n";

        mf::LogWarning("OpDigiProperties")<<"OpDigiProperties opening low gain spread file at " << FullPath.c_str();
        std::ifstream LowGainFile(FullPath.c_str());
        if(LowGainFile.is_open()) {
          std::string line;
          while ( LowGainFile.good() ){
            getline(LowGainFile, line);
            fLowGainArray.push_back(strtod(line.c_str(),NULL));
          }
        }else throw cet::exception("OpDigiProperties")<<"Unable to open file!\n";

        FullPath="";
        if( !sp.find_file(fGainSpreadFile, FullPath) )
          throw cet::exception("OpDigiProperties") << "Unable to find low gain spread file in " << sp.to_string() << "\n";

        mf::LogWarning("OpDigiProperties")<<"OpDigiProperties opening low gain spread file at " << FullPath.c_str();
        std::ifstream GainSpreadFile(FullPath.c_str());
        if(GainSpreadFile.is_open()) {
          std::string line;
          while ( GainSpreadFile.good() ){
            getline(GainSpreadFile, line);
            fGainSpreadArray.push_back(strtod(line.c_str(),NULL));
          }
        }else throw cet::exception("OpDigiProperties")<<"Unable to open file!\n";

      }
    else{
      // Generate a set of gake gain mean & spread.
      std::string txt;
      txt+=     "    Generating gain for each pmt.\n";
      txt+=Form("        High gain mean: %g ADC/p.e.\n", fHighGainMean);
      txt+=Form("        Low  gain mean: %g ADC/p.e.\n", fLowGainMean);
      txt+=Form("        PMT-to-PMT gain spread : %g \n", fGainSpread_PMT2PMT);
      txt+=Form("        Intrinsic gain spread  : %g \n", fGainSpread);
      mf::LogWarning("OpDigiProperties")<<txt.c_str();
      for(unsigned int i=0;i<fGeometry->NOpChannels();++i) {
        fLowGainArray.push_back(CLHEP::RandGauss::shoot(fLowGainMean,fLowGainMean*fGainSpread_PMT2PMT));
        fHighGainArray.push_back(CLHEP::RandGauss::shoot(fHighGainMean,fHighGainMean*fGainSpread_PMT2PMT));
        fGainSpreadArray.push_back(fGainSpread);
      }
    }

    //
    // Note:
    // Check for # entries. These exceptions ensures we have enough # of elements.
    // When a user access the elements by a channel number using LowGainMean(Channel_t ch) etc.,
    // those functions do not check if the given channel number is valid or not to keep a high
    // performance of the code. If there's an invalid memory access error in run-time, then
    // it must mean the user provided an invalid channel number and not due to insufficient
    // vector elements filled in this function.
    //
    if(fLowGainArray.size()<fGeometry->NOpChannels())
      throw cet::exception("OpDigiProperties")<<"Low gain missing for some channels!\n";
    if(fHighGainArray.size()<fGeometry->NOpChannels())
      throw cet::exception("OpDigiProperties")<<"High gain missing for some channels!\n";
    if(fGainSpreadArray.size()<fGeometry->NOpChannels())
      throw cet::exception("OpDigiProperties")<<"Gain spread missing for some channels!\n";
  }

void opdet::OpDigiProperties::FillPedMeanArray ( )

private

Definition at line 187 of file OpDigiProperties_service.cc.

                                         {
    for(unsigned int i=0;i<fGeometry->NOpChannels();++i)
      fPedMeanArray.push_back((optdata::ADC_Count_t)(CLHEP::RandGauss::shoot((double)fADCBaseline,fADCBaseSpread)+0.5));
  }

double opdet::OpDigiProperties::GainSpread ( ) const

inlinenoexcept

Returns set value for intrinsic gain spread (common to HIGH and LOW)

Definition at line 77 of file OpDigiProperties.h.

{ return fGainSpread;         }

double opdet::OpDigiProperties::GainSpread_PMT2PMT ( ) const

inlinenoexcept

Returns set value for PMT-to-PMT gain spread.

Definition at line 79 of file OpDigiProperties.h.

{ return fGainSpread_PMT2PMT; }

std::vector<double> const& opdet::OpDigiProperties::GainSpreadArray ( ) const

inlinenoexcept

Returns an array of gain spread.

Definition at line 92 of file OpDigiProperties.h.

{ return fGainSpreadArray; }

std::vector< double > opdet::OpDigiProperties::GenAnalyticalWF ( )

private

Definition at line 326 of file OpDigiProperties_service.cc.

                                                     {
    mf::LogWarning("OpDigiProperties")<<"    OpDigiProperties using analytical function for WF generation.";
    //
    // Generate waveform from analytical form
    //
    if(!fAnalyticalSPE) {// Create analytical function if not yet created
      fAnalyticalSPE = new TF1("_analyticalSPE",
                               "10^(22)*x^[1]*[0]*exp(-x/[2])/TMath::Factorial([1])",
                               0,fWFLength*2);             // x is time in micro-seconds
      fAnalyticalSPE->SetParameter(0,fVoltageAmpForSPE);             // voltage amplitude for s.p.e.
      fAnalyticalSPE->SetParameter(1,fWFPowerFactor);                // power factor (no unit)
      fAnalyticalSPE->SetParameter(2,fWFTimeConstant);               // time constant in us
    }
    //
    // Define a waveform vector
    //
    // Size of WF = time width [us] * frequency [MHz]
    std::vector<double> PEWaveform(int( fWFLength * fSampleFreq), 0.0);
    double SamplingDuration = 1./fSampleFreq; // in micro seconds
    double MaxAmp=0;
    double Charge=0;
    for(unsigned short i = 0; i<PEWaveform.size(); ++i){
      double Value=fAnalyticalSPE->Integral(   i   * SamplingDuration,
                                              (i+1) * SamplingDuration) / SamplingDuration;
      PEWaveform[i]=Value;
      if(PEWaveform[i]>MaxAmp) MaxAmp=PEWaveform[i];
      Charge+=Value;
    }

    // rescale
    if(MaxAmp<=0) throw cet::exception("OpDigiProperties_module")<<"Waveform amplitude <=0!\n";

    if(!fChargeNormalized) {
      for(unsigned short i=0; i<PEWaveform.size(); i++) {
        PEWaveform[i]=PEWaveform[i]/MaxAmp;
        if(PEWaveform[i]<1.e-4) PEWaveform[i]=0;
      }
    }
    else for(unsigned short i=0; i<PEWaveform.size(); i++){ PEWaveform[i]=PEWaveform[i]/Charge; }

    return PEWaveform;
  }

std::vector< double > opdet::OpDigiProperties::GenEmpiricalWF ( std::string  WaveformFile )

private

Definition at line 291 of file OpDigiProperties_service.cc.

  {
    // Read in waveform vector from text file
    std::ifstream WaveformFile (fWaveformFile.c_str());
    std::string line;

    mf::LogWarning("OpDigiProperties")<<"OpDigiProperties opening OpDet waveform at " << fWaveformFile.c_str();

    // Read in each line and place into waveform vector
    std::vector<double> PEWaveform;
    if (WaveformFile.is_open())
      {
        double MaxAmp=0;
        double Charge=0;
        int    NSample=0;
        while ( WaveformFile.good() && NSample<int(fWFLength*fSampleFreq))
          {
            getline (WaveformFile, line);
            double Amp=strtod(line.c_str(),NULL);
            PEWaveform.push_back(Amp);
            if(Amp>MaxAmp) MaxAmp=Amp;
            NSample++;
            Charge+=Amp;
          }
        // rescale
        if(MaxAmp<=0) throw cet::exception("OpDigiProperties_module")<<"Waveform amplitude <=0!\n";
        if(!fChargeNormalized)for(unsigned short i=0; i<PEWaveform.size(); i++){ PEWaveform[i]=PEWaveform[i]/MaxAmp; }
        else for(unsigned short i=0; i<PEWaveform.size(); i++){ PEWaveform[i]=PEWaveform[i]/Charge; }
      }
    else throw cet::exception("No Waveform File") << "Unable to open file\n";

    WaveformFile.close();
    return(PEWaveform);
  }

void opdet::OpDigiProperties::GenerateWaveform ( )

private

Definition at line 276 of file OpDigiProperties_service.cc.

  {
    std::string FullPath;

    if(fUseEmpiricalShape){
      cet::search_path sp("FW_SEARCH_PATH");
      if( !sp.find_file(fWaveformFile, FullPath) )

        throw cet::exception("OpDigiProperties") << "Unable to find PMT waveform file in " << sp.to_string() << "\n";

      fWaveform   = GenEmpiricalWF(FullPath);

    }else fWaveform = GenAnalyticalWF();
  }

double opdet::OpDigiProperties::GetSPEAmplitude

(

)

Utility function ... To be verified (Kazu 08/05/13)

Definition at line 95 of file OpDigiProperties_service.cc.

  {
    double AmpSoFar=0;
    for(size_t i=0; i!=fWaveform.size(); ++i)
      if(fWaveform.at(i)>AmpSoFar) AmpSoFar=fWaveform.at(i);
    return AmpSoFar;
  }

double opdet::OpDigiProperties::GetSPEArea

(

)

Utility function ... To be verified (Kazu 08/05/13)

Definition at line 86 of file OpDigiProperties_service.cc.

  {
    double SPEArea=0;
    for(size_t i=0; i!=fWaveform.size(); ++i)
      SPEArea+=fWaveform.at(i);
    return SPEArea;
  }

double opdet::OpDigiProperties::GetSPECumulativeAmplitude

(

)

Utility function ... To be verified (Kazu 08/05/13)

Definition at line 116 of file OpDigiProperties_service.cc.

  {
    double AmpSoFar=0, Cumulative=0;
    for(size_t i=0; i!=fWaveform.size(); ++i)
      {
        Cumulative += fWaveform.at(i);
        if(Cumulative>AmpSoFar) AmpSoFar=Cumulative;
      }
    return AmpSoFar;
  }

double opdet::OpDigiProperties::GetSPECumulativeArea

(

)

Utility function ... To be verified (Kazu 08/05/13)

Definition at line 104 of file OpDigiProperties_service.cc.

  {
    double Cumulative=0, SPEArea=0;
    for(size_t i=0; i!=fWaveform.size(); ++i)
      {
        Cumulative += fWaveform.at(i);
        SPEArea    += Cumulative;
      }
    return SPEArea;
  }

optdata::TimeSlice_t opdet::OpDigiProperties::GetTimeSlice

(

double

time_ns

)

Convert the given time into time-slice number. Input time should be in ns unit and measurd w.r.t. MC photon T0

Definition at line 151 of file OpDigiProperties_service.cc.

  {
    if( time_ns/1.e3 > (fTimeEnd-fTimeBegin)) return std::numeric_limits<optdata::TimeSlice_t>::max();

    else return optdata::TimeSlice_t((time_ns/1.e3-fTimeBegin)*fSampleFreq);
  }

double opdet::OpDigiProperties::HighGain

(

optdata::Channel_t

ch

)

const

Generate & return HIGH gain value for an input channel using mean & spread for this channel.

Definition at line 146 of file OpDigiProperties_service.cc.

  {
    return CLHEP::RandGauss::shoot(fHighGainArray[ch],fGainSpreadArray[ch]*fHighGainArray[ch]);
  }

std::vector<double> const& opdet::OpDigiProperties::HighGainArray ( ) const

inlinenoexcept

Returns an array of HIGH gain.

Definition at line 88 of file OpDigiProperties.h.

{ return fHighGainArray;   }

double opdet::OpDigiProperties::HighGainMean ( ) const

inlinenoexcept

Returns set mean gain value for HIGH gain.

Definition at line 70 of file OpDigiProperties.h.

{ return fHighGainMean; }

double opdet::OpDigiProperties::HighGainMean

(

optdata::Channel_t

ch

)

const

Returns generated HIGH gain value for input channel (PMT-to-PMT spread applied)

Definition at line 136 of file OpDigiProperties_service.cc.

  {
    return fHighGainArray[ch];
  }

double opdet::OpDigiProperties::LowGain

(

optdata::Channel_t

ch

)

const

Generate & return LOW gain value for an input channel using mean & spread for this channel.

Definition at line 141 of file OpDigiProperties_service.cc.

  {
    return CLHEP::RandGauss::shoot(fLowGainArray[ch],fGainSpreadArray[ch]*fLowGainArray[ch]);
  }

std::vector<double> const& opdet::OpDigiProperties::LowGainArray ( ) const

inlinenoexcept

Returns an array of LOW gain.

Definition at line 90 of file OpDigiProperties.h.

{ return fLowGainArray;    }

double opdet::OpDigiProperties::LowGainMean ( ) const

inlinenoexcept

Returns set mean gain value for LOW gain.

Definition at line 68 of file OpDigiProperties.h.

{ return fLowGainMean;  }

double opdet::OpDigiProperties::LowGainMean

(

optdata::Channel_t

ch

)

const

Returns generated LOW gain value for input channel (PMT-to-PMT spread applied)

Definition at line 130 of file OpDigiProperties_service.cc.

  {
    return fLowGainArray[ch];
  }

optdata::ADC_Count_t opdet::OpDigiProperties::PedFlucAmp ( ) const

inlinenoexcept

Returns amplitude of pedestal fluctuation.

Definition at line 58 of file OpDigiProperties.h.

{ return fPedFlucAmp;      }

double opdet::OpDigiProperties::PedFlucRate ( ) const

inlinenoexcept

Returns rate of pedestal fluctuation.

Definition at line 56 of file OpDigiProperties.h.

{ return fPedFlucRate;     }

std::vector<optdata::ADC_Count_t> const& opdet::OpDigiProperties::PedMeanArray ( ) const

inlinenoexcept

Returns an array of generated pedestal mean value per channel.

Definition at line 94 of file OpDigiProperties.h.

{return fPedMeanArray;}

double opdet::OpDigiProperties::QE ( ) const

inlinenoexcept

Returns quantum efficiency.

Definition at line 52 of file OpDigiProperties.h.

{ return fQE;              }

double opdet::OpDigiProperties::SampleFreq ( ) const

inlinenoexcept

Returns sample frequency in MHz.

Definition at line 33 of file OpDigiProperties.h.

{ return fSampleFreq; }

optdata::ADC_Count_t opdet::OpDigiProperties::SaturationScale ( ) const

inlinenoexcept

Returns the saturation scale of the electronics.

Definition at line 60 of file OpDigiProperties.h.

{ return fSaturationScale; }

std::vector<double> const& opdet::OpDigiProperties::SinglePEWaveform ( ) const

inlinenoexcept

Returns a vector of double which represents a binned SPE waveform.

Definition at line 86 of file OpDigiProperties.h.

{ return fWaveform;        }

double opdet::OpDigiProperties::TimeBegin ( ) const

inlinenoexcept

Returns window start time in us ... with respect to MC photon T0.

Definition at line 35 of file OpDigiProperties.h.

{ return fTimeBegin;  }

double opdet::OpDigiProperties::TimeEnd ( ) const

inlinenoexcept

Returns window end time in us ... with respect to MC photon T0.

Definition at line 37 of file OpDigiProperties.h.

{ return fTimeEnd;    }

std::vector< double > opdet::OpDigiProperties::WaveformInit

(

std::string

WaveformFile

)

DEPRECATED as far as Kazu is concerned for UBooNE PMT simulation. Any comment? –Kazu 08/05/2013

Definition at line 162 of file OpDigiProperties_service.cc.

  {
    // Read in waveform vector from text file
    std::ifstream WaveformFile (fWaveformFile.c_str());
    std::string line;

    mf::LogInfo("OpDigiProperties")<<"OpDigiProperties opening OpDet waveform at " << fWaveformFile.c_str();

    // Read in each line and place into waveform vector
    std::vector<double> PEWaveform;
    if (WaveformFile.is_open())
      {
        while ( WaveformFile.good() )
          {
            getline (WaveformFile, line);
            PEWaveform.push_back( fPERescale * strtod( line.c_str(), NULL ) );
          }
      }
    else throw cet::exception("OpDigiProperties") << "No Waveform File: Unable to open file\n";

    WaveformFile.close();
    return(PEWaveform);
  }

double opdet::OpDigiProperties::WFPowerFactor ( ) const

inlinenoexcept

Returns WF power factor used in analytical model.

Definition at line 66 of file OpDigiProperties.h.

{ return fWFPowerFactor;   }

double opdet::OpDigiProperties::WFTimeConstant ( ) const

inlinenoexcept

Returns WF time constant used in analytical model.

Definition at line 64 of file OpDigiProperties.h.

{ return fWFTimeConstant;  }

Member Data Documentation

optdata::ADC_Count_t opdet::OpDigiProperties::fADCBaseline

private

Definition at line 125 of file OpDigiProperties.h.

double opdet::OpDigiProperties::fADCBaseSpread

private

Definition at line 126 of file OpDigiProperties.h.

TF1* opdet::OpDigiProperties::fAnalyticalSPE

private

Definition at line 135 of file OpDigiProperties.h.

bool opdet::OpDigiProperties::fChargeNormalized

private

Definition at line 142 of file OpDigiProperties.h.

double opdet::OpDigiProperties::fDarkRate

private

Definition at line 123 of file OpDigiProperties.h.

double opdet::OpDigiProperties::fGainSpread

private

Definition at line 133 of file OpDigiProperties.h.

double opdet::OpDigiProperties::fGainSpread_PMT2PMT

private

Definition at line 134 of file OpDigiProperties.h.

std::vector<double> opdet::OpDigiProperties::fGainSpreadArray

private

Definition at line 145 of file OpDigiProperties.h.

std::string opdet::OpDigiProperties::fGainSpreadFile

private

Definition at line 140 of file OpDigiProperties.h.

art::ServiceHandle<geo::Geometry const> opdet::OpDigiProperties::fGeometry

private

Definition at line 147 of file OpDigiProperties.h.

std::vector<double> opdet::OpDigiProperties::fHighGainArray

private

Definition at line 144 of file OpDigiProperties.h.

std::string opdet::OpDigiProperties::fHighGainFile

private

Definition at line 137 of file OpDigiProperties.h.

double opdet::OpDigiProperties::fHighGainMean

private

Definition at line 128 of file OpDigiProperties.h.

std::vector<double> opdet::OpDigiProperties::fLowGainArray

private

Definition at line 143 of file OpDigiProperties.h.

std::string opdet::OpDigiProperties::fLowGainFile

private

Definition at line 138 of file OpDigiProperties.h.

double opdet::OpDigiProperties::fLowGainMean

private

Definition at line 127 of file OpDigiProperties.h.

optdata::ADC_Count_t opdet::OpDigiProperties::fPedFlucAmp

private

Definition at line 122 of file OpDigiProperties.h.

double opdet::OpDigiProperties::fPedFlucRate

private

Definition at line 121 of file OpDigiProperties.h.

std::vector<optdata::ADC_Count_t> opdet::OpDigiProperties::fPedMeanArray

private

Definition at line 146 of file OpDigiProperties.h.

double opdet::OpDigiProperties::fPERescale

private

Definition at line 110 of file OpDigiProperties.h.

double opdet::OpDigiProperties::fQE

private

Definition at line 120 of file OpDigiProperties.h.

double opdet::OpDigiProperties::fSampleFreq

private

Definition at line 107 of file OpDigiProperties.h.

optdata::ADC_Count_t opdet::OpDigiProperties::fSaturationScale

private

Definition at line 124 of file OpDigiProperties.h.

double opdet::OpDigiProperties::fTimeBegin

private

Definition at line 108 of file OpDigiProperties.h.

double opdet::OpDigiProperties::fTimeEnd

private

Definition at line 109 of file OpDigiProperties.h.

bool opdet::OpDigiProperties::fUseEmpiricalGain

private

Definition at line 118 of file OpDigiProperties.h.

bool opdet::OpDigiProperties::fUseEmpiricalShape

private

Definition at line 119 of file OpDigiProperties.h.

double opdet::OpDigiProperties::fVoltageAmpForSPE

private

Definition at line 129 of file OpDigiProperties.h.

std::vector<double> opdet::OpDigiProperties::fWaveform

private

Definition at line 141 of file OpDigiProperties.h.

std::string opdet::OpDigiProperties::fWaveformFile

private

Definition at line 139 of file OpDigiProperties.h.

double opdet::OpDigiProperties::fWFLength

private

Definition at line 132 of file OpDigiProperties.h.

double opdet::OpDigiProperties::fWFPowerFactor

private

Definition at line 131 of file OpDigiProperties.h.

double opdet::OpDigiProperties::fWFTimeConstant

private

Definition at line 130 of file OpDigiProperties.h.

---

The documentation for this class was generated from the following files:

• larana/larana/OpticalDetector/OpDigiProperties.h
• larana/larana/OpticalDetector/OpDigiProperties_service.cc