TrackingAnalyzer.cc

Go to the documentation of this file.

#include "DQM/TrackingMonitorClient/plugins/TrackingAnalyzer.h"

#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include "FWCore/Framework/interface/LuminosityBlock.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/ParameterSet/interface/FileInPath.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Utilities/interface/InputTag.h"

#include "DQMServices/Core/interface/DQMStore.h"

#include "DataFormats/SiStripDetId/interface/StripSubdetector.h"

#include "DataFormats/FEDRawData/interface/FEDRawDataCollection.h"
#include "DataFormats/FEDRawData/interface/FEDRawData.h"
#include "DataFormats/FEDRawData/interface/FEDNumbering.h"

#include "CondFormats/DataRecord/interface/SiStripCondDataRecords.h"
#include "CondFormats/SiStripObjects/interface/SiStripFedCabling.h"
#include "CalibTracker/Records/interface/SiStripDetCablingRcd.h"
#include "CalibFormats/SiStripObjects/interface/SiStripDetCabling.h"

#include "DQM/SiStripCommon/interface/SiStripFolderOrganizer.h"
#include "DQM/TrackingMonitorClient/interface/TrackingActionExecutor.h"
#include "DQM/TrackingMonitorClient/interface/TrackingUtility.h"

#include <iomanip>
#include <cstdio>
#include <string>
#include <sstream>
#include <cmath>

#define BUF_SIZE 256

//
// -- Constructor
//
TrackingAnalyser::TrackingAnalyser(edm::ParameterSet const& ps)
    : verbose_(ps.getUntrackedParameter<bool>("verbose", false)),
      fedCablingToken_(esConsumes<SiStripFedCabling, SiStripFedCablingRcd, edm::Transition::BeginRun>()),
      detCablingToken_(esConsumes<SiStripDetCabling, SiStripDetCablingRcd, edm::Transition::BeginRun>()) {
  edm::LogInfo("TrackingAnalyser") << "Creating TrackingAnalyser ";

  // Get TkMap ParameterSet
  //  tkMapPSet_ = ps.getParameter<edm::ParameterSet>("TkmapParameters");

  std::string localPath = std::string("DQM/TrackingMonitorClient/test/loader.html");
  std::ifstream fin(edm::FileInPath(localPath).fullPath().c_str(), std::ios::in);
  char buf[BUF_SIZE];

  if (!fin) {
    edm::LogError("TrackingAnalyzer") << "Input File: loader.html"
                                      << " could not be opened!" << std::endl;
    return;
  }

  while (fin.getline(buf, BUF_SIZE, '\n')) {  // pops off the newline character
    html_out_ << buf;
  }
  fin.close();

  staticUpdateFrequency_ = ps.getUntrackedParameter<int>("StaticUpdateFrequency", 1);
  globalStatusFilling_ = ps.getUntrackedParameter<int>("GlobalStatusFilling", 1);
  shiftReportFrequency_ = ps.getUntrackedParameter<int>("ShiftReportFrequency", 1);

  edm::InputTag rawDataTag = ps.getUntrackedParameter<edm::InputTag>("RawDataTag");
  rawDataToken_ = consumes<FEDRawDataCollection>(rawDataTag);

  // instantiate web interface
  actionExecutor_ = new TrackingActionExecutor(ps);
  trackerFEDsFound_ = false;
  trackerFEDsWdataFound_ = false;

  nFEDinfoDir_ = ps.getParameter<std::string>("nFEDinfoDir");
  nFEDinVsLSname_ = ps.getParameter<std::string>("nFEDinVsLSname");
  nFEDinWdataVsLSname_ = ps.getParameter<std::string>("nFEDinWdataVsLSname");
}
//
// -- Destructor
//
TrackingAnalyser::~TrackingAnalyser() { edm::LogInfo("TrackingAnalyser") << "Deleting TrackingAnalyser "; }
//
// -- Begin Job
//
void TrackingAnalyser::beginJob() { nLumiSecs_ = 0; }
//
// -- Begin Run
//
void TrackingAnalyser::beginRun(edm::Run const& run, edm::EventSetup const& eSetup) {
  edm::LogInfo("TrackingAnalyser") << " Begining of Run";

  // Check latest Fed cabling and create TrackerMapCreator
  if (fedCablingWatcher_.check(eSetup)) {  //this should check if cabling record has changed
    edm::LogInfo("TrackingAnalyser") << "beginRun: "
                                     << " Change in Cabling, recrated TrackerMap";
    fedCabling_ = &eSetup.getData(fedCablingToken_);
    detCabling_ = &eSetup.getData(detCablingToken_);
  }
}
//
// -- Begin Luminosity Block
//
void TrackingAnalyser::dqmBeginLuminosityBlock(DQMStore::IBooker& ibooker_,
                                               DQMStore::IGetter& igetter_,
                                               edm::LuminosityBlock const& lumiSeg,
                                               edm::EventSetup const& eSetup) {
  edm::LogInfo("TrackingAnalyser") << " Begin of LS transition";
}

//
// -- End Luminosity Block
//
void TrackingAnalyser::dqmEndLuminosityBlock(DQMStore::IBooker& ibooker_,
                                             DQMStore::IGetter& igetter_,
                                             edm::LuminosityBlock const& lumiSeg,
                                             edm::EventSetup const& eSetup) {
  edm::LogInfo("TrackingAnalyser") << " End of LS transition, performing the DQM client operation";

  nLumiSecs_++;

  if (verbose_)
    edm::LogInfo("TrackingAnalyser") << "[TrackingAnalyser::endLuminosityBlock] globalStatusFilling_ "
                                     << (globalStatusFilling_ ? "YES" : "NOPE") << std::endl;
  if (globalStatusFilling_)
    actionExecutor_->createGlobalStatus(ibooker_, igetter_);

  double iLS = lumiSeg.id().luminosityBlock();
  checkTrackerFEDsInLS(igetter_, iLS);
  checkTrackerFEDsWdataInLS(igetter_, iLS);
  if (verbose_)
    edm::LogInfo("TrackingAnalyser") << "endLuminosityBlock trackerFEDsFound_ " << (trackerFEDsFound_ ? "YES" : "NOPE")
                                     << std::endl;
  if (verbose_)
    edm::LogInfo("TrackingAnalyser") << "endLuminosityBlock trackerFEDsWdataFound_ "
                                     << (trackerFEDsWdataFound_ ? "YES" : "NOPE") << std::endl;

  if (!trackerFEDsFound_) {
    actionExecutor_->fillDummyGlobalStatus();
    actionExecutor_->createDummyShiftReport();
  } else {
    if (trackerFEDsWdataFound_) {
      actionExecutor_->fillGlobalStatus(ibooker_, igetter_);
      if (shiftReportFrequency_ != -1)
        actionExecutor_->createShiftReport(ibooker_, igetter_);
    }
  }

  if (!trackerFEDsFound_) {
    actionExecutor_->fillDummyLSStatus();
    return;
  }

  if (verbose_)
    edm::LogInfo("TrackingAnalyser") << "====================================================== " << std::endl;
  if (verbose_)
    edm::LogInfo("TrackingAnalyser") << " ===> Iteration # " << nLumiSecs_ << " " << lumiSeg.luminosityBlock()
                                     << std::endl;
  if (verbose_)
    edm::LogInfo("TrackingAnalyser") << "====================================================== " << std::endl;
}

//
// -- End Job
//
void TrackingAnalyser::dqmEndJob(DQMStore::IBooker& ibooker_, DQMStore::IGetter& igetter_) {
  edm::LogInfo("TrackingAnalyser") << " endjob called!";

  if (globalStatusFilling_)
    actionExecutor_->createGlobalStatus(ibooker_, igetter_);
  // Fill Global Status
  if (globalStatusFilling_ > 0) {
    actionExecutor_->fillGlobalStatus(ibooker_, igetter_);
  }
}
//
// Check Tracker FEDs
//
void TrackingAnalyser::checkTrackerFEDsInLS(DQMStore::IGetter& igetter, double iLS) {
  double nFEDinLS = 0.;
  MonitorElement* tmpME = igetter.get(nFEDinfoDir_ + "/" + nFEDinVsLSname_);
  if (tmpME) {
    TProfile* tmpP = tmpME->getTProfile();
    size_t ibin = tmpP->GetXaxis()->FindBin(iLS);
    if (verbose_)
      edm::LogInfo("TrackingAnalyser") << "iLS: " << iLS << " ibin: " << ibin;
    nFEDinLS = tmpME->getBinContent(ibin);
    if (verbose_)
      edm::LogInfo("TrackingAnalyser") << " ---> nFEDinLS: " << nFEDinLS;
  }

  trackerFEDsFound_ = (nFEDinLS > 0);
  if (verbose_)
    edm::LogInfo("TrackingAnalyser") << " ---> trackerFEDsFound_: " << trackerFEDsFound_ << std::endl;
}

void TrackingAnalyser::checkTrackerFEDsWdataInLS(DQMStore::IGetter& igetter, double iLS) {
  double nFEDinLS = 0.;
  MonitorElement* tmpME = igetter.get(nFEDinfoDir_ + "/" + nFEDinWdataVsLSname_);
  if (verbose_)
    edm::LogInfo("TrackingAnalyser") << "found " << nFEDinfoDir_ << "/" << nFEDinWdataVsLSname_ << " ? "
                                     << (tmpME ? "YES" : "NOPE") << std::endl;
  if (tmpME) {
    TProfile* tmpP = tmpME->getTProfile();
    size_t ibin = tmpP->GetXaxis()->FindBin(iLS);
    if (verbose_)
      edm::LogInfo("TrackingAnalyser") << "iLS: " << iLS << " ibin: " << ibin;
    nFEDinLS = tmpME->getBinContent(ibin);
    if (verbose_)
      edm::LogInfo("TrackingAnalyser") << " ---> nFEDinLS: " << nFEDinLS;
  }

  trackerFEDsWdataFound_ = (nFEDinLS > 0);
  if (verbose_)
    edm::LogInfo("TrackingAnalyser") << " ---> trackerFEDsWdataFound_: " << trackerFEDsWdataFound_ << std::endl;
}

#include "FWCore/Framework/interface/MakerMacros.h"
DEFINE_FWK_MODULE(TrackingAnalyser);