DriftCache.cxx

Go to the documentation of this file.

////////////////////////////////////////////////////////////////////////
/// \brief   Hold drift constants from current run
/// \author  bckhouse@caltech.edu
/// \date
/////////////////////////////////////////////////////////////////////////
#include "Calibrator/func/DriftCache.h"

#include "messagefacility/MessageLogger/MessageLogger.h"

#ifdef NOVACMAKE
#include "Database/cxx/include/Table.h"
#else
#include "Database/Table.h"
#endif

#include <cassert>
#include "cetlib_except/exception.h"

namespace calib
{
  //......................................................................
  DriftCache::DriftCache()
    : fTableLoaded(false),
      fUseEpochs(false),
      fDet(0),
      fFactorIdx(-1),
      fRun(0),
      fCSVFile(""),
      fTag(""),
      fEpochTag(""),
      fEpoch(""),      
      fTbl(0)
  {
  }

  //......................................................................
  DriftCache::~DriftCache()
  {
    if (fTbl) delete fTbl;
  }

  //......................................................................
  double DriftCache::GetFactor(int plane, int cell)
  {
    if (!fTableLoaded) Load();
    
    geo::OfflineChan chan(plane,cell);
    int tchan = 12*plane + int(cell/32) + 1;
    
    //**********
    nova::dbi::Row* row = fTbl->GetVldRow(tchan, fRun);
    if (!row) {
      std::cout << "DriftCache no factor found for channel " << tchan << std::endl;
      std::abort();
    }
    double factor;
    if(!row->Col(fFactorIdx).Get(factor)) std::abort();
    return factor;
  }

  //......................................................................
  void DriftCache::SetRun(int det, int run)
  {
    fDet = det;
    
    fRun =run;

    if (fTbl) delete fTbl;
    
    fTbl = new nova::dbi::Table("Calibration/tables/DriftCorrVld.xml");
    fTbl->SetDetector(fDet);
    fTbl->SetVerbosity(100);
    fTbl->SetMinTSVld(fRun);
    fTbl->SetMaxTSVld(fRun);
    fTbl->SetDataTypeMask(nova::dbi::kDataOnly);
    fTbl->SetTag(fTag);
    fFactorIdx = fTbl->GetColIndex("corr");
    fTableLoaded = false;
    
  }

  //......................................................................
  void DriftCache::GetEpochFromDB()
  {
    // first get epoch
    if (fUseEpochs) {
      nova::dbi::Table tbl("Calibration/tables/CalibDriftEpochs.xml");
      tbl.SetDataTypeMask(nova::dbi::kDataOnly);
      tbl.SetTag(fEpochTag);
      tbl.SetDetector(fDet);
      tbl.SetMinTSVld(fRun);
      tbl.SetMaxTSVld(fRun);
      tbl.SetVerbosity(100);
      assert(tbl.Load());
      assert(tbl.NRow() == 1);
      
      int epochIdx = tbl.GetColIndex("epoch");
      tbl.GetRow(0)->Col(epochIdx).Get(fEpoch);
    }
  }

  //......................................................................
  std::string DriftCache::GetCSVPath()
  {

    // if configured to provide own file
    if(!fCSVFile.empty()){   
      return fCSVFile.c_str();
    }
    
    // necessary components of path
    const char* ups_path = getenv("CALIBCSVS_CSV_PATH");
    assert(ups_path);
    const char* detStr = "fd";

    if (fUseEpochs) {
      GetEpochFromDB();
      assert(!fEpoch.empty());
    }
    
    std::string fullPath = std::string(ups_path) + "/" + fTag + "/calib_drift." + detStr + "." + fTag + "." + fEpoch + ".csv";
    //  Return final path
    
    return fullPath;    
  }
  
  //......................................................................
  void DriftCache::Load()
  {
    if (fTableLoaded) return;
    
    mf::LogInfo("DriftCache")
      << "Initializing drift cache ...\n";

    const std::string constsCSVSource = GetCSVPath();
    
    GetEpochFromDB(); // will access dB if configured to do so

    if (fTbl->LoadFromCSV(constsCSVSource) <= 0){
      std::cout<< "DriftCache Load failed to read from CSV file " <<std::endl;
      std::abort();
    }
    if (fTbl->NRow() == 0){
      std::cout<< "DriftCache Load found 0 rows in CSV file " <<std::endl;
      std::abort();
    }
    
    fTableLoaded = true;
    fTbl->FillChanRowMap();
  }
} // end namespace calib
////////////////////////////////////////////////////////////////////////