#include "OpDigiProperties.h"

Public Member Functions
	OpDigiProperties (fhicl::ParameterSet const &pset, art::ActivityRegistry &reg)

	~OpDigiProperties ()

void	reconfigure (fhicl::ParameterSet const &p)

double	SampleFreq ()
	Returns sample frequency in MHz. More...

double	TimeBegin ()
	Returns window start time in us ... with respect to MC photon T0. More...

double	TimeEnd ()
	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
	Returns quantum efficiency. More...

double	DarkRate () const
	Returns rate of dark noise. More...

double	PedFlucRate () const
	Returns rate of pedestal fluctuation. More...

optdata::ADC_Count_t	PedFlucAmp () const
	Returns amplitude of pedestal fluctuation. More...

optdata::ADC_Count_t	SaturationScale () const
	Returns the saturation scale of the electronics. More...

optdata::ADC_Count_t	ADCBaseline () const
	Returns the ADCBaseline set mean value. More...

double	WFTimeConstant () const
	Returns WF time constant used in analytical model. More...

double	WFPowerFactor () const
	Returns WF power factor used in analytical model. More...

double	LowGainMean () const
	Returns set mean gain value for LOW gain. More...

double	HighGainMean () const
	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
	Returns set value for intrinsic gain spread (common to HIGH and LOW) More...

double	GainSpread_PMT2PMT () const
	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 >	SinglePEWaveform () const
	Returns a vector of double which represents a binned SPE waveform. More...

std::vector< double >	HighGainArray () const
	Returns an array of HIGH gain. More...

std::vector< double >	LowGainArray () const
	Returns an array of LOW gain. More...

std::vector< double >	GainSpreadArray () const
	Returns an array of gain spread. More...

std::vector< optdata::ADC_Count_t >	PedMeanArray () const
	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 >	fGeometry

Detailed Description

Definition at line 29 of file OpDigiProperties.h.

Constructor & Destructor Documentation

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

Definition at line 24 of file OpDigiProperties_service.cc.

References reconfigure().

     : fAnalyticalSPE(0)
   {
     this->reconfigure(pset);
   }

opdet::OpDigiProperties::~OpDigiProperties ( )

Definition at line 32 of file OpDigiProperties_service.cc.

   {
     
   }

Member Function Documentation

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

inline

Returns the ADCBaseline set mean value.

Definition at line 67 of file OpDigiProperties.h.

References fADCBaseline.

67 { return fADCBaseline; }

opdet::OpDigiProperties::fADCBaseline

optdata::ADC_Count_t fADCBaseline

Definition: OpDigiProperties.h:130

double opdet::OpDigiProperties::DarkRate ( ) const

inline

Returns rate of dark noise.

Definition at line 59 of file OpDigiProperties.h.

References fDarkRate.

59 { return fDarkRate; }

opdet::OpDigiProperties::fDarkRate

double fDarkRate

Definition: OpDigiProperties.h:128

void opdet::OpDigiProperties::FillGainArray ( )

private

Definition at line 209 of file OpDigiProperties_service.cc.

References fGainSpread, fGainSpread_PMT2PMT, fGainSpreadArray, fGainSpreadFile, fGeometry, fHighGainArray, fHighGainFile, fHighGainMean, fLowGainArray, fLowGainFile, fLowGainMean, fUseEmpiricalGain, and geo::GeometryCore::NOpChannels().

Referenced by reconfigure().

                                       {
     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 203 of file OpDigiProperties_service.cc.

References fADCBaseline, fADCBaseSpread, fGeometry, fPedMeanArray, and geo::GeometryCore::NOpChannels().

Referenced by reconfigure().

                                          {
     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

inline

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

Definition at line 82 of file OpDigiProperties.h.

References fGainSpread.

82 { return fGainSpread; }

opdet::OpDigiProperties::fGainSpread

double fGainSpread

Definition: OpDigiProperties.h:138

double opdet::OpDigiProperties::GainSpread_PMT2PMT ( ) const

inline

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

Definition at line 84 of file OpDigiProperties.h.

References fGainSpread_PMT2PMT, HighGain(), and LowGain().

84 { return fGainSpread_PMT2PMT; }

opdet::OpDigiProperties::fGainSpread_PMT2PMT

double fGainSpread_PMT2PMT

Definition: OpDigiProperties.h:139

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

inline

Returns an array of gain spread.

Definition at line 97 of file OpDigiProperties.h.

References fGainSpreadArray.

97 { return fGainSpreadArray; }

opdet::OpDigiProperties::fGainSpreadArray

std::vector< double > fGainSpreadArray

Definition: OpDigiProperties.h:150

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

private

Definition at line 342 of file OpDigiProperties_service.cc.

References DEFINE_ART_SERVICE, e, fAnalyticalSPE, fChargeNormalized, fSampleFreq, fVoltageAmpForSPE, fWFLength, fWFPowerFactor, and fWFTimeConstant.

Referenced by GenerateWaveform().

                                                      {
     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 307 of file OpDigiProperties_service.cc.

References fChargeNormalized, fSampleFreq, and fWFLength.

Referenced by GenerateWaveform().

   {
     // 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 292 of file OpDigiProperties_service.cc.

References fUseEmpiricalShape, fWaveform, fWaveformFile, GenAnalyticalWF(), and GenEmpiricalWF().

Referenced by reconfigure().

   {
     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 48 of file OpDigiProperties_service.cc.

References fWaveform.

Referenced by PedMeanArray().

   {
     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 39 of file OpDigiProperties_service.cc.

References fWaveform.

Referenced by PedMeanArray().

   {
     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 69 of file OpDigiProperties_service.cc.

References fWaveform.

Referenced by PedMeanArray().

   {
     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 57 of file OpDigiProperties_service.cc.

References fWaveform.

Referenced by PedMeanArray().

   {
     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 104 of file OpDigiProperties_service.cc.

References fSampleFreq, fTimeBegin, fTimeEnd, and max.

Referenced by TimeEnd().

   {
     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 99 of file OpDigiProperties_service.cc.

References fGainSpreadArray, and fHighGainArray.

Referenced by GainSpread_PMT2PMT().

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

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

inline

Returns an array of HIGH gain.

Definition at line 93 of file OpDigiProperties.h.

References fHighGainArray.

93 { return fHighGainArray; }

opdet::OpDigiProperties::fHighGainArray

std::vector< double > fHighGainArray

Definition: OpDigiProperties.h:149

double opdet::OpDigiProperties::HighGainMean ( ) const

inline

Returns set mean gain value for HIGH gain.

Definition at line 75 of file OpDigiProperties.h.

References fHighGainMean, and LowGainMean().

75 { return fHighGainMean; }

opdet::OpDigiProperties::fHighGainMean

double fHighGainMean

Definition: OpDigiProperties.h:133

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 89 of file OpDigiProperties_service.cc.

References fHighGainArray.

   {
     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 94 of file OpDigiProperties_service.cc.

References fGainSpreadArray, and fLowGainArray.

Referenced by GainSpread_PMT2PMT().

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

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

inline

Returns an array of LOW gain.

Definition at line 95 of file OpDigiProperties.h.

References fLowGainArray.

95 { return fLowGainArray; }

opdet::OpDigiProperties::fLowGainArray

std::vector< double > fLowGainArray

Definition: OpDigiProperties.h:148

double opdet::OpDigiProperties::LowGainMean ( ) const

inline

Returns set mean gain value for LOW gain.

Definition at line 73 of file OpDigiProperties.h.

References fLowGainMean.

Referenced by HighGainMean().

73 { return fLowGainMean; }

opdet::OpDigiProperties::fLowGainMean

double fLowGainMean

Definition: OpDigiProperties.h:132

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 83 of file OpDigiProperties_service.cc.

References fLowGainArray.

   {
     return fLowGainArray[ch];
   }

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

inline

Returns amplitude of pedestal fluctuation.

Definition at line 63 of file OpDigiProperties.h.

References fPedFlucAmp.

63 { return fPedFlucAmp; }

opdet::OpDigiProperties::fPedFlucAmp

optdata::ADC_Count_t fPedFlucAmp

Definition: OpDigiProperties.h:127

double opdet::OpDigiProperties::PedFlucRate ( ) const

inline

Returns rate of pedestal fluctuation.

Definition at line 61 of file OpDigiProperties.h.

References fPedFlucRate.

61 { return fPedFlucRate; }

opdet::OpDigiProperties::fPedFlucRate

double fPedFlucRate

Definition: OpDigiProperties.h:126

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

inline

Returns an array of generated pedestal mean value per channel.

Definition at line 99 of file OpDigiProperties.h.

References fPedMeanArray, GetSPEAmplitude(), GetSPEArea(), GetSPECumulativeAmplitude(), and GetSPECumulativeArea().

99 {return fPedMeanArray;}

opdet::OpDigiProperties::fPedMeanArray

std::vector< optdata::ADC_Count_t > fPedMeanArray

Definition: OpDigiProperties.h:151

double opdet::OpDigiProperties::QE ( ) const

inline

Returns quantum efficiency.

Definition at line 57 of file OpDigiProperties.h.

References fQE.

Referenced by opdet::FlashHypothesisCreator::CreateFlashHypothesesFromSegment(), cosmic::BeamFlashTrackMatchTaggerAlg::GetMIPHypotheses(), and opdet::SimPhotonCounterAlg::InitializeCounters().

57 { return fQE; }

opdet::OpDigiProperties::fQE

double fQE

Definition: OpDigiProperties.h:125

void opdet::OpDigiProperties::reconfigure ( fhicl::ParameterSet const & p )

Definition at line 113 of file OpDigiProperties_service.cc.

References fADCBaseline, fADCBaseSpread, fChargeNormalized, fDarkRate, fGainSpread, fGainSpread_PMT2PMT, fGainSpreadFile, fGeometry, fHighGainArray, fHighGainFile, fHighGainMean, FillGainArray(), FillPedMeanArray(), fLowGainArray, fLowGainFile, fLowGainMean, fPedFlucAmp, fPedFlucRate, fPedMeanArray, fPERescale, fQE, fSampleFreq, fSaturationScale, fTimeBegin, fTimeEnd, fUseEmpiricalGain, fUseEmpiricalShape, fVoltageAmpForSPE, fWaveformFile, fWFLength, fWFPowerFactor, fWFTimeConstant, GenerateWaveform(), fhicl::ParameterSet::get(), and geo::GeometryCore::NOpChannels().

Referenced by OpDigiProperties().

   {
     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();
       
     return;
   }

double opdet::OpDigiProperties::SampleFreq ( )

inline

Returns sample frequency in MHz.

Definition at line 38 of file OpDigiProperties.h.

References fSampleFreq.

38 { return fSampleFreq; }

opdet::OpDigiProperties::fSampleFreq

double fSampleFreq

Definition: OpDigiProperties.h:112

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

inline

Returns the saturation scale of the electronics.

Definition at line 65 of file OpDigiProperties.h.

References fSaturationScale.

65 { return fSaturationScale; }

opdet::OpDigiProperties::fSaturationScale

optdata::ADC_Count_t fSaturationScale

Definition: OpDigiProperties.h:129

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

inline

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

Definition at line 91 of file OpDigiProperties.h.

References fWaveform.

91 { return fWaveform; }

opdet::OpDigiProperties::fWaveform

std::vector< double > fWaveform

Definition: OpDigiProperties.h:146

double opdet::OpDigiProperties::TimeBegin ( )

inline

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

Definition at line 40 of file OpDigiProperties.h.

References fTimeBegin.

40 { return fTimeBegin; }

opdet::OpDigiProperties::fTimeBegin

double fTimeBegin

Definition: OpDigiProperties.h:113

double opdet::OpDigiProperties::TimeEnd ( )

inline

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

Definition at line 42 of file OpDigiProperties.h.

References fTimeEnd, GetTimeSlice(), and WaveformInit().

42 { return fTimeEnd; }

opdet::OpDigiProperties::fTimeEnd

double fTimeEnd

Definition: OpDigiProperties.h:114

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 178 of file OpDigiProperties_service.cc.

References fPERescale.

Referenced by TimeEnd().

   {
     // 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

inline

Returns WF power factor used in analytical model.

Definition at line 71 of file OpDigiProperties.h.

References fWFPowerFactor.

71 { return fWFPowerFactor; }

opdet::OpDigiProperties::fWFPowerFactor

double fWFPowerFactor

Definition: OpDigiProperties.h:136

double opdet::OpDigiProperties::WFTimeConstant ( ) const

inline

Returns WF time constant used in analytical model.

Definition at line 69 of file OpDigiProperties.h.

References fWFTimeConstant.

69 { return fWFTimeConstant; }

opdet::OpDigiProperties::fWFTimeConstant

double fWFTimeConstant

Definition: OpDigiProperties.h:135

Member Data Documentation

optdata::ADC_Count_t opdet::OpDigiProperties::fADCBaseline

private

Definition at line 130 of file OpDigiProperties.h.

Referenced by ADCBaseline(), FillPedMeanArray(), and reconfigure().

double opdet::OpDigiProperties::fADCBaseSpread

private

Definition at line 131 of file OpDigiProperties.h.

Referenced by FillPedMeanArray(), and reconfigure().

TF1* opdet::OpDigiProperties::fAnalyticalSPE

private

Definition at line 140 of file OpDigiProperties.h.

Referenced by GenAnalyticalWF().

bool opdet::OpDigiProperties::fChargeNormalized

private

Definition at line 147 of file OpDigiProperties.h.

Referenced by GenAnalyticalWF(), GenEmpiricalWF(), and reconfigure().

double opdet::OpDigiProperties::fDarkRate

private

Definition at line 128 of file OpDigiProperties.h.

Referenced by DarkRate(), and reconfigure().

double opdet::OpDigiProperties::fGainSpread

private

Definition at line 138 of file OpDigiProperties.h.

Referenced by FillGainArray(), GainSpread(), and reconfigure().

double opdet::OpDigiProperties::fGainSpread_PMT2PMT

private

Definition at line 139 of file OpDigiProperties.h.

Referenced by FillGainArray(), GainSpread_PMT2PMT(), and reconfigure().

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

private

Definition at line 150 of file OpDigiProperties.h.

Referenced by FillGainArray(), GainSpreadArray(), HighGain(), and LowGain().

std::string opdet::OpDigiProperties::fGainSpreadFile

private

Definition at line 145 of file OpDigiProperties.h.

Referenced by FillGainArray(), and reconfigure().

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

private

Definition at line 152 of file OpDigiProperties.h.

Referenced by FillGainArray(), FillPedMeanArray(), and reconfigure().

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

private

Definition at line 149 of file OpDigiProperties.h.

Referenced by FillGainArray(), HighGain(), HighGainArray(), HighGainMean(), and reconfigure().

std::string opdet::OpDigiProperties::fHighGainFile

private

Definition at line 142 of file OpDigiProperties.h.

Referenced by FillGainArray(), and reconfigure().

double opdet::OpDigiProperties::fHighGainMean

private

Definition at line 133 of file OpDigiProperties.h.

Referenced by FillGainArray(), HighGainMean(), and reconfigure().

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

private

Definition at line 148 of file OpDigiProperties.h.

Referenced by FillGainArray(), LowGain(), LowGainArray(), LowGainMean(), and reconfigure().

std::string opdet::OpDigiProperties::fLowGainFile

private

Definition at line 143 of file OpDigiProperties.h.

Referenced by FillGainArray(), and reconfigure().

double opdet::OpDigiProperties::fLowGainMean

private

Definition at line 132 of file OpDigiProperties.h.

Referenced by FillGainArray(), LowGainMean(), and reconfigure().

optdata::ADC_Count_t opdet::OpDigiProperties::fPedFlucAmp

private

Definition at line 127 of file OpDigiProperties.h.

Referenced by PedFlucAmp(), and reconfigure().

double opdet::OpDigiProperties::fPedFlucRate

private

Definition at line 126 of file OpDigiProperties.h.

Referenced by PedFlucRate(), and reconfigure().

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

private

Definition at line 151 of file OpDigiProperties.h.

Referenced by FillPedMeanArray(), PedMeanArray(), and reconfigure().

double opdet::OpDigiProperties::fPERescale

private

Definition at line 115 of file OpDigiProperties.h.

Referenced by reconfigure(), and WaveformInit().

double opdet::OpDigiProperties::fQE

private

Definition at line 125 of file OpDigiProperties.h.

Referenced by QE(), and reconfigure().

double opdet::OpDigiProperties::fSampleFreq

private

Definition at line 112 of file OpDigiProperties.h.

Referenced by GenAnalyticalWF(), GenEmpiricalWF(), GetTimeSlice(), reconfigure(), and SampleFreq().

optdata::ADC_Count_t opdet::OpDigiProperties::fSaturationScale

private

Definition at line 129 of file OpDigiProperties.h.

Referenced by reconfigure(), and SaturationScale().

double opdet::OpDigiProperties::fTimeBegin

private

Definition at line 113 of file OpDigiProperties.h.

Referenced by GetTimeSlice(), reconfigure(), and TimeBegin().

double opdet::OpDigiProperties::fTimeEnd

private

Definition at line 114 of file OpDigiProperties.h.

Referenced by GetTimeSlice(), reconfigure(), and TimeEnd().

bool opdet::OpDigiProperties::fUseEmpiricalGain

private

Definition at line 123 of file OpDigiProperties.h.

Referenced by FillGainArray(), and reconfigure().

bool opdet::OpDigiProperties::fUseEmpiricalShape

private

Definition at line 124 of file OpDigiProperties.h.

Referenced by GenerateWaveform(), and reconfigure().

double opdet::OpDigiProperties::fVoltageAmpForSPE

private

Definition at line 134 of file OpDigiProperties.h.

Referenced by GenAnalyticalWF(), and reconfigure().

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

private

Definition at line 146 of file OpDigiProperties.h.

Referenced by GenerateWaveform(), GetSPEAmplitude(), GetSPEArea(), GetSPECumulativeAmplitude(), GetSPECumulativeArea(), and SinglePEWaveform().

std::string opdet::OpDigiProperties::fWaveformFile

private

Definition at line 144 of file OpDigiProperties.h.

Referenced by GenerateWaveform(), and reconfigure().

double opdet::OpDigiProperties::fWFLength

private

Definition at line 137 of file OpDigiProperties.h.

Referenced by GenAnalyticalWF(), GenEmpiricalWF(), and reconfigure().

double opdet::OpDigiProperties::fWFPowerFactor

private

Definition at line 136 of file OpDigiProperties.h.

Referenced by GenAnalyticalWF(), reconfigure(), and WFPowerFactor().

double opdet::OpDigiProperties::fWFTimeConstant

private

Definition at line 135 of file OpDigiProperties.h.

Referenced by GenAnalyticalWF(), reconfigure(), and WFTimeConstant().

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

/cvmfs/larsoft.opensciencegrid.org/products/larana/v07_04_02/source/larana/OpticalDetector/OpDigiProperties.h
/cvmfs/larsoft.opensciencegrid.org/products/larana/v07_04_02/source/larana/OpticalDetector/OpDigiProperties_service.cc

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation