Inheritance diagram for lightana::DriftEstimatorPMTRatio:

Classes
struct	Config

Public Member Functions
	DriftEstimatorPMTRatio (art::ToolConfigTable< Config > const &config)

double	GetDriftPosition (std::vector< double > PE_v) override

double	GetPropagationTime (double drift) override

double	PEToPropagationTime (std::vector< double > PE_v)

Private Member Functions
double	Interpolate (double val)

Private Member Functions inherited from lightana::DriftEstimatorBase
virtual	~DriftEstimatorBase () noexcept=default

Private Attributes
std::string	fCalibrationFile

double	fVGroupVUV

double	fVGroupVIS

double	fVGroupVUV_I

double	fDriftDistance

double	fVISLightPropTime

double	fKinkDistance

std::vector< double >	fPMTRatioCal

std::vector< double >	fDriftCal

int	fNCalBins

double	fPMTRatio_MinVal

double	fPMTRatio_MaxVal

opdet::sbndPDMapAlg	fPDSMap

Detailed Description

Definition at line 36 of file DriftEstimatorPMTRatio_tool.cc.

Constructor & Destructor Documentation

lightana::DriftEstimatorPMTRatio::DriftEstimatorPMTRatio ( art::ToolConfigTable< Config > const & config )

explicit

Definition at line 100 of file DriftEstimatorPMTRatio_tool.cc.

     : fCalibrationFile { config().CalibrationFile() },
     fVGroupVUV { config().VGroupVUV() },
     fVGroupVIS { config().VGroupVIS() }
   {
 
     // Open input file
     std::string file_name;
     cet::search_path sp("FW_SEARCH_PATH");
     if ( !sp.find_file(fCalibrationFile, file_name) )
       throw cet::exception("DriftEstimatorPMTRatio") << "Calibration file " <<
           fCalibrationFile << " not found in FW_SEARCH_PATH\n";
 
     TFile* input_file = TFile::Open(file_name.c_str(), "READ");
     TProfile * hProf_Calibration = (TProfile*)input_file->Get("PMTRatioCalibrationProfile");
 
     //Fill calibration variables
     fNCalBins = hProf_Calibration->GetNbinsX();
     for (int ix=1; ix<=fNCalBins; ix++){
       fPMTRatioCal.push_back( hProf_Calibration->GetBinCenter(ix) );
       fDriftCal.push_back( hProf_Calibration->GetBinContent(ix) );
     }
     fPMTRatio_MinVal = hProf_Calibration->GetBinCenter(1);
     fPMTRatio_MaxVal = hProf_Calibration->GetBinCenter(fNCalBins);
 
     input_file->Close();
 
     geo::GeometryCore const& geom = *(lar::providerFrom<geo::Geometry>());
 
     fDriftDistance = geom.TPC(0.0).DriftDistance();
     fVISLightPropTime = fDriftDistance/fVGroupVIS;
     fKinkDistance = 0.5*fDriftDistance*(1-fVGroupVUV/fVGroupVIS);
 
     fVGroupVUV_I = 1./fVGroupVUV;
   }

Member Function Documentation

double lightana::DriftEstimatorPMTRatio::GetDriftPosition ( std::vector< double > PE_v )

overridevirtual

Implements lightana::DriftEstimatorBase.

Definition at line 136 of file DriftEstimatorPMTRatio_tool.cc.

                                                                        {
 
     double coatedPE=0, uncoatedPE=0;
 
     for(size_t oc=0; oc<PE_v.size(); oc++){
       if( fPDSMap.pdType(oc)=="pmt_coated" ) coatedPE+=PE_v[oc];
       else if( fPDSMap.pdType(oc)=="pmt_uncoated" ) uncoatedPE+=PE_v[oc];
     }
 
     double pmtratio=1.;
     if(coatedPE!=0)
       pmtratio = 4*uncoatedPE/coatedPE;
 
     double drift_distance;
     if(pmtratio<=fPMTRatioCal[0])
       drift_distance=fDriftCal[0];
     else if(pmtratio>=fPMTRatioCal[fNCalBins-1])
       drift_distance=fDriftCal[fNCalBins-1];
     else
       drift_distance=Interpolate(pmtratio);
 
     return drift_distance;
    }

double lightana::DriftEstimatorPMTRatio::GetPropagationTime ( double drift )

overridevirtual

Implements lightana::DriftEstimatorBase.

Definition at line 160 of file DriftEstimatorPMTRatio_tool.cc.

                                                                {
 
     // drift is here the X coordinate
     // cathode: x=0 cm, PDS: x=200 cm
     if(std::abs(drift) > fKinkDistance)
       return (fDriftDistance-std::abs(drift)) * fVGroupVUV_I ;
     else
       return std::abs(drift) * fVGroupVUV_I + fVISLightPropTime;
   }

double lightana::DriftEstimatorPMTRatio::Interpolate ( double val )

private

Definition at line 177 of file DriftEstimatorPMTRatio_tool.cc.

                                                       {
 
     size_t upix = std::upper_bound(fPMTRatioCal.begin(), fPMTRatioCal.end(), val)-fPMTRatioCal.begin();
 
     double slope = ( fDriftCal[upix]-fDriftCal[upix] ) / ( fPMTRatioCal[upix]-fPMTRatioCal[upix-1] );
 
     return fDriftCal[upix-1] + slope * ( val - fPMTRatioCal[upix-1] );
   }

double lightana::DriftEstimatorPMTRatio::PEToPropagationTime ( std::vector< double > PE_v )

Definition at line 170 of file DriftEstimatorPMTRatio_tool.cc.

                                                                           {
 
     double _drift = GetDriftPosition(PE_v);
 
     return GetPropagationTime(_drift);
   }