DetectorClocksServiceStandard_service.cc

Go to the documentation of this file.

#include "lardata/DetectorInfoServices/DetectorClocksServiceStandard.h"
#include "lardataalg/DetectorInfo/DetectorClocksStandardTriggerLoader.h"
#include "TFile.h"
#include "art/Framework/IO/Root/RootDB/SQLite3Wrapper.h"
#include "fhiclcpp/make_ParameterSet.h"

//-----------------------------------------------------------------------------------------
detinfo::DetectorClocksServiceStandard::DetectorClocksServiceStandard
  (fhicl::ParameterSet const& pset, art::ActivityRegistry &reg)
  : fClocks(std::make_unique<detinfo::DetectorClocksStandard>(pset))
{
  
  reg.sPreProcessEvent.watch (this, &DetectorClocksServiceStandard::preProcessEvent);
  reg.sPostOpenFile.watch    (this, &DetectorClocksServiceStandard::postOpenFile);
  reg.sPreBeginRun.watch     (this, &DetectorClocksServiceStandard::preBeginRun);
  
}

//------------------------------------------------------------------
void detinfo::DetectorClocksServiceStandard::reconfigure(fhicl::ParameterSet const& pset)
//------------------------------------------------------------------
{
  fClocks->Configure(pset);

}

//------------------------------------------------------------
void detinfo::DetectorClocksServiceStandard::preProcessEvent(const art::Event& evt)
//------------------------------------------------------------
{
  detinfo::setDetectorClocksStandardTrigger(*fClocks, evt);
}

//------------------------------------------------------
void detinfo::DetectorClocksServiceStandard::preBeginRun(art::Run const& run)
//------------------------------------------------------
{

  fClocks->ApplyParams();

}


//---------------------------------------------------------------
void detinfo::DetectorClocksServiceStandard::postOpenFile(const std::string& filename)
//---------------------------------------------------------------
{

  // Method inheriting from DetectorProperties
  
  if(!fClocks->InheritClockConfig()) return;

  // The only way to access art service metadata from the input file
  // is to open it as a separate TFile object.  Do that now.
  
  if(!filename.empty()) {
    
    TFile* file = TFile::Open(filename.c_str(), "READ");
    if(file != 0 && !file->IsZombie() && file->IsOpen()) {

      std::vector<std::string> cfgName = fClocks->ConfigNames();
      std::vector<double> cfgValue = fClocks->ConfigValues();
      
      // Open the sqlite datatabase.
      
      art::SQLite3Wrapper sqliteDB(file, "RootFileDB");
      
      // Loop over all stored ParameterSets.

      std::vector<size_t> config_count(detinfo::kInheritConfigTypeMax,0);
      std::vector<double> config_value(detinfo::kInheritConfigTypeMax,0);

      sqlite3_stmt * stmt = 0;
      sqlite3_prepare_v2(sqliteDB, "SELECT PSetBlob from ParameterSets;", -1, &stmt, NULL);

      while (sqlite3_step(stmt) == SQLITE_ROW) {

        fhicl::ParameterSet ps;
        fhicl::make_ParameterSet(reinterpret_cast<char const *>(sqlite3_column_text(stmt, 0)), ps);

        if(!fClocks->IsRightConfig(ps)) continue;

        for(size_t i=0; i<detinfo::kInheritConfigTypeMax; ++i) {

          double value_from_file = ps.get<double>(cfgName.at(i).c_str());
          
          if(!(config_count.at(i)))

            config_value.at(i) = value_from_file;

          else if(config_value.at(i) != value_from_file)

            throw cet::exception(__FUNCTION__) << Form("\033[95mFound historical value disagreement for %s ... %g != %g",
                                                       cfgName.at(i).c_str(),
                                                       config_value.at(i),
                                                       value_from_file)
                                               << "\033[00m" << std::endl;
          config_count.at(i) +=1;

        }
        
      }

      // Override parameters 

      for(size_t i=0; i<detinfo::kInheritConfigTypeMax; ++i)

        if(config_count.at(i) && cfgValue.at(i) != config_value.at(i)) {

          std::cout  << Form("\033[93mOverriding configuration parameter %s ... %g (fcl) => %g (data file)\033[00m",
                             cfgName.at(i).c_str(),
                             cfgValue.at(i),
                             config_value.at(i))
                     << std::endl;
          
          fClocks->SetConfigValue(i,config_value.at(i));
          
        }
    }
    
    // Close file.
    if(file != 0) {
      if(file->IsOpen())
        file->Close();
      delete file;
    }
  }

  // Reset parameters
  fClocks->ApplyParams();

}

DEFINE_ART_SERVICE_INTERFACE_IMPL(detinfo::DetectorClocksServiceStandard, detinfo::DetectorClocksService)