DQMDcsInfoClient.cc

Go to the documentation of this file.

/*
 * \file DQMDcsInfoClient.cc
 * \author Andreas Meyer
 * Last Update: 2019-03-26 Migrate to DQMEDHarvester
 *
*/

#include "DQMDcsInfoClient.h"
#include "FWCore/Framework/interface/LuminosityBlock.h"

// -----------------------------
//  constructors and destructor
// -----------------------------

DQMDcsInfoClient::DQMDcsInfoClient( const edm::ParameterSet& ps ) {

  parameters_ = ps;

  subsystemname_ = parameters_.getUntrackedParameter<std::string>("subSystemFolder", "Info") ;
  dcsinfofolder_ = parameters_.getUntrackedParameter<std::string>("dcsInfoFolder", "DcsInfo") ;

}


DQMDcsInfoClient::~DQMDcsInfoClient() = default;

void 
DQMDcsInfoClient::beginRun(const edm::Run& r, const edm::EventSetup& c) 
{
  DCS.clear();
  DCS.resize(10);  // start with 10 LS, resize later
  processedLS_.clear();
}

void 
DQMDcsInfoClient::dqmEndLuminosityBlock(DQMStore::IBooker& ibooker, DQMStore::IGetter& igetter, const edm::LuminosityBlock& l, const edm::EventSetup& c)
{
  unsigned int nlumi = l.id().luminosityBlock() ;
  processedLS_.insert(nlumi);
  // cout << " in lumi section " << nlumi << endl;

  if (nlumi+1 > DCS.size()) 
     DCS.resize(nlumi+1);
  // cout << "DCS size: " << DCS.size() << endl;     

  MonitorElement* DCSbyLS_ = 
            igetter.get(subsystemname_ + "/" + dcsinfofolder_ + "/DCSbyLS" ); 

  if ( DCSbyLS_ ) 
  {
    if ( TH1F * h1 = DCSbyLS_->getTH1F()) 
    {
      int word = 0;
      for (int i = 0; i < 25 ; i++) 
      {

    if ( h1->GetBinContent(i+1) != 0 ) 
           word |= (0x0 << i); // set to 0 because HV was off (!)
    else 
           word |= (0x1 << i); // set to 1 because HV was on (!)
      }
      DCS[nlumi] = word;
    }
  }
  return; 
}

void 
DQMDcsInfoClient::dqmEndJob(DQMStore::IBooker& ibooker, DQMStore::IGetter& igetter)
{
  // Fetch GlobalTag information and fill the string/ME.
  ibooker.cd();  
  ibooker.setCurrentFolder(subsystemname_ +"/CMSSWInfo/");
  const edm::ParameterSet &globalTagPSet =
    edm::getProcessParameterSetContainingModule(moduleDescription())
    .getParameterSet("PoolDBESSource@GlobalTag");

  ibooker.bookString("globalTag_Harvesting", globalTagPSet.getParameter<std::string>("globaltag"));

  ibooker.cd();  
  ibooker.setCurrentFolder(subsystemname_ +"/EventInfo/");

  unsigned int nlsmax = DCS.size();
  reportSummary_=ibooker.bookFloat("reportSummary");
  reportSummary_->Fill(1.);
  
  reportSummaryMap_ = igetter.get(subsystemname_ +"/EventInfo/reportSummaryMap");
  assert(!reportSummaryMap_); // this would be a configuration problem

  reportSummaryMap_ = ibooker.book2D("reportSummaryMap",
                     "HV and GT vs Lumi", nlsmax, 0., nlsmax, 25, 0., 25.);
  if (processedLS_.empty()) return;
  unsigned int lastProcessedLS = *(--processedLS_.end());
  meProcessedLS_ = ibooker.book1D("ProcessedLS",
                                  "Processed Lumisections",
                                  lastProcessedLS+1,
                                  0.,lastProcessedLS+1);
  reportSummaryMap_->setBinLabel( 1, "CSC+"    , 2);
  reportSummaryMap_->setBinLabel( 2, "CSC-"    , 2);
  reportSummaryMap_->setBinLabel( 3, "DT0"     , 2);
  reportSummaryMap_->setBinLabel( 4, "DT+"     , 2);
  reportSummaryMap_->setBinLabel( 5, "DT-"     , 2);
  reportSummaryMap_->setBinLabel( 6, "EB+"     , 2);
  reportSummaryMap_->setBinLabel( 7, "EB-"     , 2);
  reportSummaryMap_->setBinLabel( 8, "EE+"     , 2);
  reportSummaryMap_->setBinLabel( 9, "EE-"     , 2);
  reportSummaryMap_->setBinLabel(10, "ES+"     , 2);
  reportSummaryMap_->setBinLabel(11, "ES-"     , 2);
  reportSummaryMap_->setBinLabel(12, "HBHEa"   , 2);
  reportSummaryMap_->setBinLabel(13, "HBHEb"   , 2);
  reportSummaryMap_->setBinLabel(14, "HBHEc"   , 2);
  reportSummaryMap_->setBinLabel(15, "HF"      , 2);
  reportSummaryMap_->setBinLabel(16, "HO"      , 2);
  reportSummaryMap_->setBinLabel(17, "BPIX"    , 2);
  reportSummaryMap_->setBinLabel(18, "FPIX"    , 2);
  reportSummaryMap_->setBinLabel(19, "RPC"     , 2);
  reportSummaryMap_->setBinLabel(20, "TIBTID"  , 2);
  reportSummaryMap_->setBinLabel(21, "TOB"     , 2);
  reportSummaryMap_->setBinLabel(22, "TECp"    , 2);
  reportSummaryMap_->setBinLabel(23, "TECm"    , 2);
  reportSummaryMap_->setBinLabel(24, "CASTOR"  , 2);
  reportSummaryMap_->setBinLabel(25, "PhysDecl", 2);
  reportSummaryMap_->setAxisTitle("Luminosity Section");

  // fill
  for (unsigned int i = 0 ; i < DCS.size() ; i++ )
  {
    for ( int j = 0 ; j < 25 ; j++ ) 
    {
      if (DCS[i] & (0x1 << j))
        reportSummaryMap_->setBinContent(i,j+1,1.);
      else
        reportSummaryMap_->setBinContent(i,j+1,0.);
    }
  }

  unsigned int lastAccessed = 0;
  
  for (auto ls : processedLS_)
  {
    while (lastAccessed < ls)
    {
      //      std::cout << "Filling " << lastAccessed << " with -1" << std::endl; 
      meProcessedLS_->Fill(lastAccessed, -1.);
      lastAccessed++;
    }
    //    std::cout << "Filling " << *it << " with 1" << std::endl; 
    meProcessedLS_->Fill(ls);
    lastAccessed = ls+1;
  }
}