SiStripAnalyser.cc

Go to the documentation of this file.

/*
 * \file SiStripAnalyser.cc
 *
 * \author  S. Dutta INFN-Pisa
 *
 */

#include "SiStripAnalyser.h"

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

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

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

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

#include "DQM/SiStripCommon/interface/SiStripFolderOrganizer.h"
#include "DQM/SiStripMonitorClient/interface/SiStripUtility.h"

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

SiStripAnalyser::SiStripAnalyser(edm::ParameterSet const& ps)
  : condDataMon_{ps},
    actionExecutor_{ps},
    tkMapPSet_{ps.getParameter<edm::ParameterSet>("TkmapParameters")},
    summaryFrequency_{
      ps.getUntrackedParameter<int>("SummaryCreationFrequency", 1)},
    staticUpdateFrequency_{
      ps.getUntrackedParameter<int>("StaticUpdateFrequency", 1)},
    globalStatusFilling_{
      ps.getUntrackedParameter<int>("GlobalStatusFilling", 1)},
    shiftReportFrequency_{
      ps.getUntrackedParameter<int>("ShiftReportFrequency", 1)},
    rawDataToken_{consumes<FEDRawDataCollection>(
      ps.getUntrackedParameter<edm::InputTag>("RawDataTag"))},
    printFaultyModuleList_{
      ps.getUntrackedParameter<bool>("PrintFaultyModuleList", true)}
{
  std::string const localPath{"DQM/SiStripMonitorClient/test/loader.html"};
  std::ifstream fin{edm::FileInPath(localPath).fullPath(), std::ios::in};
  if (!fin) {
    std::cerr << "Input File: loader.html"
              << " could not be opened!" << std::endl;
    return;
  }
  edm::LogInfo("SiStripAnalyser")
    << " SiStripAnalyser::Creating SiStripAnalyser ";
}

SiStripAnalyser::~SiStripAnalyser()
{
  edm::LogInfo("SiStripAnalyser")
    << "SiStripAnalyser::Deleting SiStripAnalyser ";
}

void
SiStripAnalyser::beginJob()
{
  // Read the summary configuration file
  if (!actionExecutor_.readConfiguration()) {
    edm::LogInfo("SiStripAnalyser")
      << "SiStripAnalyser:: Error to read configuration file!! Summary will "
         "not be produced!!!";
  }
}

void
SiStripAnalyser::beginRun(edm::Run const& run, edm::EventSetup const& eSetup)
{
  edm::LogInfo("SiStripAnalyser") << "SiStripAnalyser:: Begining of Run";

  // Check latest Fed cabling and create TrackerMapCreator
  unsigned long long const cacheID =
    eSetup.get<SiStripFedCablingRcd>().cacheIdentifier();
  if (m_cacheID_ != cacheID) {
    m_cacheID_ = cacheID;
    edm::LogInfo("SiStripAnalyser")
      << "SiStripAnalyser::beginRun: "
      << " Change in Cabling, recrated TrackerMap";
    if (!actionExecutor_.readTkMapConfiguration(eSetup)) {
      edm::LogInfo("SiStripAnalyser")
        << "SiStripAnalyser:: Error to read configuration file!! TrackerMap "
           "will not be produced!!!";
    }
    eSetup.get<SiStripFedCablingRcd>().get(fedCabling_);
    eSetup.get<SiStripDetCablingRcd>().get(detCabling_);
  }
  condDataMon_.beginRun(run.run(), eSetup);
  if (globalStatusFilling_) {
    auto& dqm_store = *edm::Service<DQMStore>{};
    actionExecutor_.createStatus(dqm_store);
  }
}

void
SiStripAnalyser::beginLuminosityBlock(edm::LuminosityBlock const& lumiSeg,
                                      edm::EventSetup const& eSetup)
{
  edm::LogInfo("SiStripAnalyser") << "SiStripAnalyser:: Begin of LS transition";
}

void
SiStripAnalyser::analyze(edm::Event const& e, edm::EventSetup const& eSetup)
{
  ++nEvents_;
  if (nEvents_ == 1 && globalStatusFilling_ > 0) {
    checkTrackerFEDs(e);
    if (!trackerFEDsFound_) {
      actionExecutor_.fillDummyStatus();
      actionExecutor_.createDummyShiftReport();
    } else {
      auto& dqm_store = *edm::Service<DQMStore>{};
      actionExecutor_.fillStatus(dqm_store, detCabling_, eSetup);
      if (shiftReportFrequency_ != -1)
        actionExecutor_.createShiftReport(dqm_store);
    }
  }
}

void
SiStripAnalyser::endLuminosityBlock(edm::LuminosityBlock const& lumiSeg,
                                    edm::EventSetup const& eSetup)
{
  edm::LogInfo("SiStripAnalyser") << "SiStripAnalyser:: End of LS transition, "
                                     "performing the DQM client operation";
  ++nLumiSecs_;

  if (!trackerFEDsFound_) {
    actionExecutor_.fillDummyStatus();
    return;
  }
  endLumiAnalysisOn_ = true;

  std::cout << "====================================================== "
            << std::endl;
  std::cout << " ===> Iteration # " << nLumiSecs_ << " "
            << lumiSeg.luminosityBlock() << std::endl;
  std::cout << "====================================================== "
            << std::endl;

  auto& dqm_store = *edm::Service<DQMStore>{};
  // Fill Global Status
  if (globalStatusFilling_ > 0) {
    actionExecutor_.fillStatus(dqm_store, detCabling_, eSetup);
  }
  // -- Create summary monitor elements according to the frequency
  if (summaryFrequency_ != -1 && nLumiSecs_ > 0 &&
      nLumiSecs_ % summaryFrequency_ == 0) {
    std::cout << " Creating Summary " << std::endl;
    actionExecutor_.createSummary(dqm_store);
  }
  endLumiAnalysisOn_ = false;
}

void
SiStripAnalyser::endRun(edm::Run const& run, edm::EventSetup const& eSetup)
{
  edm::LogInfo("SiStripAnalyser") << "SiStripAnalyser:: End of Run";
}

void
SiStripAnalyser::endJob()
{
  edm::LogInfo("SiStripAnalyser") << "SiStripAnalyser:: endjob called!";
  if (printFaultyModuleList_) {
    std::ostringstream str_val;
    auto& dqm_store = *edm::Service<DQMStore>{};
    actionExecutor_.printFaultyModuleList(dqm_store, str_val);
    std::cout << str_val.str() << std::endl;
  }
}

void
SiStripAnalyser::checkTrackerFEDs(edm::Event const& e)
{
  edm::Handle<FEDRawDataCollection> rawDataHandle;
  e.getByToken(rawDataToken_, rawDataHandle);
  if (!rawDataHandle.isValid()) return;

  auto const& rawDataCollection = *rawDataHandle;
  constexpr int siStripFedIdMin{FEDNumbering::MINSiStripFEDID};
  constexpr int siStripFedIdMax{FEDNumbering::MAXSiStripFEDID};

  for (int i = siStripFedIdMin; i <= siStripFedIdMax; ++i) {
    auto const& fedData = rawDataCollection.FEDData(i);
    if (fedData.size() && fedData.data()) {
      trackerFEDsFound_ = true;
      return;
    }
  }
}
#include "FWCore/Framework/interface/MakerMacros.h"
DEFINE_FWK_MODULE(SiStripAnalyser);