SiStripSummaryBuilder.cc

Go to the documentation of this file.

#include "CondTools/SiStrip/plugins/SiStripSummaryBuilder.h"
#include "CalibTracker/SiStripCommon/interface/SiStripDetInfoFileReader.h"
#include <iostream>
#include <fstream>
 
SiStripSummaryBuilder::SiStripSummaryBuilder(const edm::ParameterSet& iConfig)
    : fp_(iConfig.getUntrackedParameter<edm::FileInPath>(
          "file", edm::FileInPath("CalibTracker/SiStripCommon/data/SiStripDetInfo.dat"))),
      printdebug_(iConfig.getUntrackedParameter<uint32_t>("printDebug", 1)),
      iConfig_(iConfig) {}
 
void SiStripSummaryBuilder::analyze(const edm::Event& evt, const edm::EventSetup& iSetup) {
  unsigned int run = evt.id().run();
  edm::LogInfo("SiStripSummaryBuilder") << "... creating dummy SiStripSummary Data for Run " << run << "\n "
                                        << std::endl;
 
  SiStripSummary* obj = new SiStripSummary();
  obj->setRunNr(run);
 
  //* DISCOVER SET OF HISTOGRAMS & QUANTITIES TO BE UPLOADED*//
 
  std::vector<std::string> userDBContent;
  typedef std::vector<edm::ParameterSet> VParameters;
  VParameters histoList = iConfig_.getParameter<VParameters>("histoList");
  VParameters::iterator ithistoList = histoList.begin();
  VParameters::iterator ithistoListEnd = histoList.end();
 
  for (; ithistoList != ithistoListEnd; ++ithistoList) {
    std::string keyName = ithistoList->getUntrackedParameter<std::string>("keyName");
    std::vector<std::string> Quantities =
        ithistoList->getUntrackedParameter<std::vector<std::string> >("quantitiesToExtract");
    for (size_t i = 0; i < Quantities.size(); ++i) {
      if (Quantities[i] == "landau") {
        userDBContent.push_back(keyName + std::string("@") + std::string("landauPeak"));
        userDBContent.push_back(keyName + std::string("@") + std::string("landauPeakErr"));
        userDBContent.push_back(keyName + std::string("@") + std::string("landauSFWHM"));
        userDBContent.push_back(keyName + std::string("@") + std::string("landauChi2NDF"));
      } else if (Quantities[i] == "gauss") {
        userDBContent.push_back(keyName + std::string("@") + std::string("gaussMean"));
        userDBContent.push_back(keyName + std::string("@") + std::string("gaussSigma"));
        userDBContent.push_back(keyName + std::string("@") + std::string("gaussChi2NDF"));
      } else if (Quantities[i] == "stat") {
        userDBContent.push_back(keyName + std::string("@") + std::string("entries"));
        userDBContent.push_back(keyName + std::string("@") + std::string("mean"));
        userDBContent.push_back(keyName + std::string("@") + std::string("rms"));
      } else {
        edm::LogError("SiStripSummaryBuilder")
            << "Quantity " << Quantities[i] << " cannot be handled\nAllowed quantities are"
            << "\n  'stat'   that includes: entries, mean, rms"
            << "\n  'landau' that includes: landauPeak, landauPeakErr, landauSFWHM, landauChi2NDF"
            << "\n  'gauss'  that includes: gaussMean, gaussSigma, gaussChi2NDF" << std::endl;
      }
    }
  }
  obj->setUserDBContent(userDBContent);
 
  std::stringstream ss1;
  ss1 << "QUANTITIES TO BE INSERTED IN DB :"
      << " \n";
  std::vector<std::string> userDBContentA = obj->getUserDBContent();
  for (size_t i = 0; i < userDBContentA.size(); ++i)
    ss1 << userDBContentA[i] << std::endl;
  edm::LogInfo("SiStripSummaryBuilder") << ss1.str();
 
  //* Loop over detids and create dummy data for each *//
 
  std::stringstream ss2;
  for (uint32_t detid = 0; detid < 430; detid++) {
    SiStripSummary::InputVector values;
    for (unsigned int i = 0; i < userDBContent.size(); i++)
      values.push_back((float)CLHEP::RandGauss::shoot(50., 30.));
 
    ss2 << "\n\tdetid " << detid;
    for (size_t j = 0; j < values.size(); ++j)
      ss2 << "\n\t\t " << userDBContent[j] << " " << values[j];
 
    obj->put(detid, values, userDBContent);
 
    // See CondFormats/SiStripObjects/SiStripSummary.h for detid definitions
 
    if (detid == 4)
      detid = 10;
    if (detid == 14)
      detid = 20;
    if (detid == 26)
      detid = 30;
    if (detid == 32)
      detid = 40;
    if (detid == 42)
      detid = 310;
    if (detid == 313)
      detid = 320;
    if (detid == 323)
      detid = 410;
    if (detid == 419)
      detid = 420;
  }
 
  edm::LogInfo("SiStripSummaryBuilder") << ss2.str();
 
  //* Insert summary informations in the DB *//
 
  edm::Service<cond::service::PoolDBOutputService> mydbservice;
 
  if (mydbservice.isAvailable()) {
    if (mydbservice->isNewTagRequest("SiStripSummaryRcd")) {
      mydbservice->createNewIOV<SiStripSummary>(
          obj, mydbservice->beginOfTime(), mydbservice->endOfTime(), "SiStripSummaryRcd");
    } else {
      mydbservice->appendSinceTime<SiStripSummary>(obj, mydbservice->currentTime(), "SiStripSummaryRcd");
    }
  } else {
    edm::LogError("SiStripSummaryBuilder") << "Service is unavailable" << std::endl;
  }
}