calib::LifetimeCalibProtoDUNE Class Reference

#include <LifetimeCalibProtoDUNE.h>

Inheritance diagram for calib::LifetimeCalibProtoDUNE:

Public Member Functions
	LifetimeCalibProtoDUNE ()
	LifetimeCalibProtoDUNE (fhicl::ParameterSet const &pset)
	LifetimeCalibProtoDUNE (LifetimeCalibProtoDUNE const &)=delete
virtual	~LifetimeCalibProtoDUNE ()=default
bool	Configure (fhicl::ParameterSet const &pset)
bool	Update (uint64_t ts=0)
virtual double	GetLifetime () override
virtual double	GetLifetimeLow () override
virtual double	GetLifetimeHigh () override
void	SetIsMC (bool v)
void	SetLifetime (double val)
void	SetLifetimeLow (double val)
void	SetLifetimeHigh (double val)
void	SetUseCondb (bool v)
void	SetLifetimeFileName (std::string f)
void	SetInterpolate (bool v)
Public Member Functions inherited from calib::LifetimeCalib
	LifetimeCalib (const LifetimeCalib &)=delete
	LifetimeCalib (LifetimeCalib &&)=delete
LifetimeCalib &	operator= (const LifetimeCalib &)=delete
LifetimeCalib &	operator= (LifetimeCalib &&)=delete
virtual	~LifetimeCalib ()=default

Protected Member Functions
bool	LoadLifetime ()
Protected Member Functions inherited from calib::LifetimeCalib
	LifetimeCalib ()=default

Protected Attributes
bool	fUseCondbLifetime
bool	fLifetimeLoaded
bool	fIsMC
bool	fInterpolate
uint64_t	fCurrentTS
std::string	fLifetimeFileName
std::string	fLifetimeDBTag
std::map< int, LifetimePurMon_t >	fLifetimePurMon

Detailed Description

Definition at line 37 of file LifetimeCalibProtoDUNE.h.

Constructor & Destructor Documentation

calib::LifetimeCalibProtoDUNE::LifetimeCalibProtoDUNE

(

)

Definition at line 31 of file LifetimeCalibProtoDUNE.cxx.

{
  fIsMC = true;
  fLifetimeLoaded = false;
  fInterpolate = false;
  fCurrentTS = 0;
  fLifetimeFileName="";
  fLifetimeDBTag="";
}

calib::LifetimeCalibProtoDUNE::LifetimeCalibProtoDUNE

(

fhicl::ParameterSet const &

pset

)

Definition at line 43 of file LifetimeCalibProtoDUNE.cxx.

{
  fIsMC = true;
  fLifetimeLoaded = false;
  fInterpolate = false;
  fCurrentTS = 0;
  fLifetimeFileName="";
  fLifetimeDBTag="";
  Configure(pset);
}

calib::LifetimeCalibProtoDUNE::LifetimeCalibProtoDUNE ( LifetimeCalibProtoDUNE const &  )

delete

virtual calib::LifetimeCalibProtoDUNE::~LifetimeCalibProtoDUNE ( )

virtualdefault

Member Function Documentation

bool calib::LifetimeCalibProtoDUNE::Configure

(

fhicl::ParameterSet const &

pset

)

Definition at line 57 of file LifetimeCalibProtoDUNE.cxx.

{  
  fIsMC             = pset.get<bool>("IsMC");
  fUseCondbLifetime = pset.get<bool>("UseCondbLifetime");
  fInterpolate      = pset.get<bool>("Interpolate");
  fLifetimeFileName = pset.get<std::string>("LifetimeFileName");
  fLifetimeDBTag    = pset.get<std::string>("LifetimeDBTag");
  return true;
}

double calib::LifetimeCalibProtoDUNE::GetLifetime ( )

overridevirtual

Implements calib::LifetimeCalib.

Definition at line 83 of file LifetimeCalibProtoDUNE.cxx.

{
  if (!fLifetimeLoaded) this->LoadLifetime();

  if (fLifetimePurMon.find(1) == fLifetimePurMon.end()) {
    mf::LogError("LifetimeCalibProtoDUNE") << "channel not found!";
    return 0.;
  }

  return fLifetimePurMon[1].center;
}

double calib::LifetimeCalibProtoDUNE::GetLifetimeHigh ( )

overridevirtual

Implements calib::LifetimeCalib.

Definition at line 96 of file LifetimeCalibProtoDUNE.cxx.

{
  if (!fLifetimeLoaded) this->LoadLifetime();

  if (fLifetimePurMon.find(1) == fLifetimePurMon.end()) {
    mf::LogError("LifetimeCalibProtoDUNE") << "channel not found!";
    return 0.;
  }

  return fLifetimePurMon[1].high;
}

double calib::LifetimeCalibProtoDUNE::GetLifetimeLow ( )

overridevirtual

Implements calib::LifetimeCalib.

Definition at line 109 of file LifetimeCalibProtoDUNE.cxx.

{
  if (!fLifetimeLoaded) this->LoadLifetime();

  if (fLifetimePurMon.find(1) == fLifetimePurMon.end()) {
    mf::LogError("LifetimeCalibProtoDUNE") << "channel not found!";
    return 0.;
  }

  return fLifetimePurMon[1].low;
}

bool calib::LifetimeCalibProtoDUNE::LoadLifetime ( )

protected

Definition at line 122 of file LifetimeCalibProtoDUNE.cxx.

{
  if (!fUseCondbLifetime) return true;

  if (fLifetimeLoaded) return true;
  
  nutools::dbi::Table LifetimePurMonTable;
  LifetimePurMonTable.SetDetector("pdunesp");
  LifetimePurMonTable.SetTableName("lifetime_purmon");
  LifetimePurMonTable.SetTableType(nutools::dbi::kConditionsTable);
  LifetimePurMonTable.SetDataTypeMask(nutools::dbi::kDataOnly);
  if (fIsMC) LifetimePurMonTable.SetDataTypeMask(nutools::dbi::kMCOnly);

  int centerIdx = LifetimePurMonTable.AddCol("center", "double");
  int lowIdx = LifetimePurMonTable.AddCol("low", "double");
  int highIdx = LifetimePurMonTable.AddCol("high", "double");
  
  if (!fInterpolate) {
    LifetimePurMonTable.SetMinTSVld(fCurrentTS); // only load one previous lifetime
    LifetimePurMonTable.SetMaxTSVld(fCurrentTS);
  }
  else {
    // load lifetime IOV: runtime +- 1 days for interpolation
    LifetimePurMonTable.SetMinTSVld(fCurrentTS-1*86400.);
    LifetimePurMonTable.SetMaxTSVld(fCurrentTS+1*86400.);
  }
  LifetimePurMonTable.SetTag(fLifetimeDBTag);
  
  LifetimePurMonTable.SetVerbosity(100);

  bool readOk = false;
  if (!fLifetimeFileName.empty())
    readOk = LifetimePurMonTable.LoadFromCSV(fLifetimeFileName);
  else
    readOk = LifetimePurMonTable.Load();
  
  if (! readOk) {
   mf::LogError("LifetimeCalibProtoDUNE") << "Load from lifetime calib database table failed.";
   throw cet::exception("LifetimeCalibProtoDUNE") << "Failed to query lifetime database. Please check URL: https://dbdata0vm.fnal.gov:9443/dune_con_prod/app" << "\n";

   return false; //std::abort();
  }

  if (LifetimePurMonTable.NRow() == 0) {
    mf::LogError("LifetimeCalibProtoDUNE") << "Number of rows in lifetime calib table is 0.  This should never be the case!";
    return false;
  }

  std::vector<double> loaded_tv;
  std::vector<double> loaded_center;
  std::vector<double> loaded_low;
  std::vector<double> loaded_high;
  loaded_tv.clear();
  loaded_center.clear();
  loaded_low.clear();
  loaded_high.clear();

  if (LifetimePurMonTable.NRow() == 1) fInterpolate = false;

  nutools::dbi::Row* row;
  uint64_t chan;
  
  for (int i=0; i<LifetimePurMonTable.NRow(); ++i) {

    LifetimePurMon_t lifetime;
    row = LifetimePurMonTable.GetRow(i);
    chan = row->Channel(); // One channel only, start with 1

    if (chan != 1) {
      mf::LogError("LifetimeCalibProtoDUNE") << "Channel numuber in lifetime calib table is not  1.  This should never be the case!";
      return false;
    }

    row->Col(centerIdx).Get(lifetime.center);
    row->Col(lowIdx).Get(lifetime.low);
    row->Col(highIdx).Get(lifetime.high);

    loaded_tv.push_back(row->VldTime());
    loaded_center.push_back(lifetime.center);
    loaded_low.push_back(lifetime.low);
    loaded_high.push_back(lifetime.high);
    
    if (!fInterpolate) fLifetimePurMon[chan] = lifetime;
  }

  if (fInterpolate && loaded_tv.size()>1) {

    LifetimePurMon_t interpolate_lifetime;

    // find the previous and following lifetime for fCurrentTS
    int t0_idx = -1; 
    int t1_idx = -1;
    double temp_0 = 1e10;
    double temp_1 = 1e10; 
    for (size_t j=0; j<loaded_tv.size(); j++) {
      if ( fCurrentTS - loaded_tv[j] >= 0 && fCurrentTS - loaded_tv[j] < temp_0 ) {
        temp_0 = fCurrentTS- loaded_tv[j];
        t0_idx = j;
      }

      if ( loaded_tv[j]-fCurrentTS >= 0 && loaded_tv[j]-fCurrentTS < temp_1 ) {
        temp_1 = loaded_tv[j]-fCurrentTS;
        t1_idx = j;
      }
    }
    
    if (t0_idx == -1) { // only lifetime measurement after fCurrentTS during selected IOV
      // use the following one only
      interpolate_lifetime.center = loaded_center[t1_idx];
      interpolate_lifetime.low = loaded_low[t1_idx];
      interpolate_lifetime.high = loaded_high[t1_idx];
    }
    if (t1_idx == -1) { // only lifetime measurement before fCurrentTS during selected IOV
      // use the previous one only
      interpolate_lifetime.center = loaded_center[t0_idx];
      interpolate_lifetime.low = loaded_low[t0_idx];
      interpolate_lifetime.high = loaded_high[t0_idx];
    }
    if (t0_idx != -1 && t1_idx != -1) {
      if (t0_idx == t1_idx) { // ideally, this is not possible
        interpolate_lifetime.center = loaded_center[t0_idx];
        interpolate_lifetime.low = loaded_low[t0_idx];
        interpolate_lifetime.high = loaded_high[t0_idx];
      }
      else {
        // do linear interpolation using previous one and following one
        // refer: https://en.wikipedia.org/wiki/Linear_interpolation
        interpolate_lifetime.center = ( loaded_center[t0_idx]*(loaded_tv[t1_idx]-fCurrentTS) + loaded_center[t1_idx]*(fCurrentTS-loaded_tv[t0_idx]) ) / (loaded_tv[t1_idx]-loaded_tv[t0_idx]);
        interpolate_lifetime.low = ( loaded_low[t0_idx]*(loaded_tv[t1_idx]-fCurrentTS) + loaded_low[t1_idx]*(fCurrentTS-loaded_tv[t0_idx]) ) / (loaded_tv[t1_idx]-loaded_tv[t0_idx]);
        interpolate_lifetime.high = ( loaded_high[t0_idx]*(loaded_tv[t1_idx]-fCurrentTS) + loaded_high[t1_idx]*(fCurrentTS-loaded_tv[t0_idx]) ) / (loaded_tv[t1_idx]-loaded_tv[t0_idx]);
      }
    }

    fLifetimePurMon[chan] = interpolate_lifetime;
  } // end if fInterpolate

  fLifetimeLoaded = true;
  return true;
}

void calib::LifetimeCalibProtoDUNE::SetInterpolate ( bool  v )

inline

Definition at line 59 of file LifetimeCalibProtoDUNE.h.

{ fInterpolate = v; }

void calib::LifetimeCalibProtoDUNE::SetIsMC ( bool  v )

inline

Definition at line 53 of file LifetimeCalibProtoDUNE.h.

{ fIsMC = v; }

void calib::LifetimeCalibProtoDUNE::SetLifetime

(

double

val

)

void calib::LifetimeCalibProtoDUNE::SetLifetimeFileName ( std::string  f )

inline

Definition at line 58 of file LifetimeCalibProtoDUNE.h.

{ fLifetimeFileName=f; }

void calib::LifetimeCalibProtoDUNE::SetLifetimeHigh

(

double

val

)

void calib::LifetimeCalibProtoDUNE::SetLifetimeLow

(

double

val

)

void calib::LifetimeCalibProtoDUNE::SetUseCondb ( bool  v )

inline

Definition at line 57 of file LifetimeCalibProtoDUNE.h.

{ fUseCondbLifetime = v; }

bool calib::LifetimeCalibProtoDUNE::Update

(

uint64_t

ts = 0

)

Definition at line 68 of file LifetimeCalibProtoDUNE.cxx.

{
  
  if (fLifetimeLoaded && ts != fCurrentTS) {
    fLifetimePurMon.clear();
    fLifetimeLoaded = false;
  }

  fCurrentTS = ts;
  // all done! 

  return true;
}

Member Data Documentation

uint64_t calib::LifetimeCalibProtoDUNE::fCurrentTS

protected

Definition at line 68 of file LifetimeCalibProtoDUNE.h.

bool calib::LifetimeCalibProtoDUNE::fInterpolate

protected

Definition at line 67 of file LifetimeCalibProtoDUNE.h.

bool calib::LifetimeCalibProtoDUNE::fIsMC

protected

Definition at line 66 of file LifetimeCalibProtoDUNE.h.

std::string calib::LifetimeCalibProtoDUNE::fLifetimeDBTag

protected

Definition at line 71 of file LifetimeCalibProtoDUNE.h.

std::string calib::LifetimeCalibProtoDUNE::fLifetimeFileName

protected

Definition at line 70 of file LifetimeCalibProtoDUNE.h.

bool calib::LifetimeCalibProtoDUNE::fLifetimeLoaded

protected

Definition at line 65 of file LifetimeCalibProtoDUNE.h.

std::map<int,LifetimePurMon_t> calib::LifetimeCalibProtoDUNE::fLifetimePurMon

protected

Definition at line 73 of file LifetimeCalibProtoDUNE.h.

bool calib::LifetimeCalibProtoDUNE::fUseCondbLifetime

protected

Definition at line 64 of file LifetimeCalibProtoDUNE.h.

---

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

• dunecalib/dunecalib/Calib/LifetimeCalibProtoDUNE.h
• dunecalib/dunecalib/Calib/LifetimeCalibProtoDUNE.cxx