L1TDTTPGClient.cc

Go to the documentation of this file.

#include "DQM/L1TMonitorClient/interface/L1TDTTPGClient.h"
 
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "TRandom.h"
 
#include <TF1.h>
#include <cstdio>
#include <sstream>
#include <cmath>
#include <TProfile.h>
#include <TProfile2D.h>
 
using namespace edm;
using namespace std;
 
L1TDTTPGClient::L1TDTTPGClient(const edm::ParameterSet &ps) {
  parameters_ = ps;
  initialize();
}
 
L1TDTTPGClient::~L1TDTTPGClient() { LogInfo("TriggerDQM") << "[TriggerDQM]: ending... "; }
 
//--------------------------------------------------------
void L1TDTTPGClient::initialize() {
  counterLS_ = 0;
  counterEvt_ = 0;
 
  // get back-end interface
  //dbe_ = Service<DQMStore>().operator->();
 
  // base folder for the contents of this job
  monitorName_ = parameters_.getUntrackedParameter<string>("monitorName", "");
  //  cout << "Monitor name = " << monitorName_ << endl;
  prescaleLS_ = parameters_.getUntrackedParameter<int>("prescaleLS", -1);
  //  cout << "DQM lumi section prescale = " << prescaleLS_ << " lumi section(s)"<< endl;
  prescaleEvt_ = parameters_.getUntrackedParameter<int>("prescaleEvt", -1);
  //  cout << "DQM event prescale = " << prescaleEvt_ << " events(s)"<< endl;
  output_dir_ = parameters_.getUntrackedParameter<string>("output_dir", "");
  //  cout << "DQM output dir = " << output_dir_ << endl;
  input_dir_ = parameters_.getUntrackedParameter<string>("input_dir", "");
  //  cout << "DQM input dir = " << input_dir_ << endl;
 
  LogInfo("TriggerDQM");
}
 
//--------------------------------------------------------
void L1TDTTPGClient::dqmEndJob(DQMStore::IBooker &ibooker, DQMStore::IGetter &igetter) {
  ibooker.setCurrentFolder(output_dir_);
 
  // booking
 
  dttpgphmapcorrf = ibooker.book2D(
      "DT_TPG_phi_map_corr_frac", "Fraction of correlated best triggers per station", 20, 1, 21, 12, 0, 12);
  dttpgphmap2ndf =
      ibooker.book2D("DT_TPG_phi_map_2nd_frac", "Fraction of second tracks per station", 20, 1, 21, 12, 0, 12);
  dttpgphmapbxf[0] =
      ibooker.book2D("DT_TPG_phi_map_bx-1_frac", "Fraction of triggers per station (BX=-1)", 20, 1, 21, 12, 0, 12);
  dttpgphmapbxf[1] =
      ibooker.book2D("DT_TPG_phi_map_bx0_frac", "Fraction of triggers per station (BX=0)", 20, 1, 21, 12, 0, 12);
  dttpgphmapbxf[2] =
      ibooker.book2D("DT_TPG_phi_map_bx+1_frac", "Fraction of triggers per station (BX=1)", 20, 1, 21, 12, 0, 12);
  setMapPhLabel(dttpgphmapcorrf);
  setMapPhLabel(dttpgphmap2ndf);
  setMapPhLabel(dttpgphmapbxf[0]);
  setMapPhLabel(dttpgphmapbxf[1]);
  setMapPhLabel(dttpgphmapbxf[2]);
 
  dttpgthmaphf = ibooker.book2D(
      "DT_TPG_theta_map_corr_frac", "Fraction of H quality best triggers per station", 15, 1, 16, 12, 0, 12);
  dttpgthmapbxf[0] =
      ibooker.book2D("DT_TPG_theta_map_bx-1_frac", "Fraction of triggers per station (BX=-1)", 15, 1, 16, 12, 0, 12);
  dttpgthmapbxf[1] =
      ibooker.book2D("DT_TPG_theta_map_bx0_frac", "Fraction of triggers per station (BX=0)", 15, 1, 16, 12, 0, 12);
  dttpgthmapbxf[2] =
      ibooker.book2D("DT_TPG_theta_map_bx+1_frac", "Fraction of triggers per station (BX=1)", 15, 1, 16, 12, 0, 12);
  setMapThLabel(dttpgthmaphf);
  setMapThLabel(dttpgthmapbxf[0]);
  setMapThLabel(dttpgthmapbxf[1]);
  setMapThLabel(dttpgthmapbxf[2]);
 
  //   cout << "L1TDTTPGClient::analyze" << endl;
  counterEvt_++;
  if (prescaleEvt_ < 1)
    return;
  if (prescaleEvt_ > 0 && counterEvt_ % prescaleEvt_ != 0)
    return;
 
  string nName = "DT_TPG_phi_best_map_corr";
  string dName = "DT_TPG_phi_best_map";
  makeRatioHisto(igetter, dttpgphmapcorrf, nName, dName);
  dName = "DT_TPG_phi_map";
  nName = "DT_TPG_phi_map_2nd";
  makeRatioHisto(igetter, dttpgphmap2ndf, nName, dName);
  nName = "DT_TPG_phi_map_bx-1";
  makeRatioHisto(igetter, dttpgphmapbxf[0], nName, dName);
  nName = "DT_TPG_phi_map_bx0";
  makeRatioHisto(igetter, dttpgphmapbxf[1], nName, dName);
  nName = "DT_TPG_phi_map_bx+1";
  makeRatioHisto(igetter, dttpgphmapbxf[2], nName, dName);
 
  nName = "DT_TPG_theta_best_map_h";
  dName = "DT_TPG_theta_best_map";
  makeRatioHisto(igetter, dttpgthmaphf, nName, dName);
  dName = "DT_TPG_theta_map";
  nName = "DT_TPG_theta_map_bx-1";
  makeRatioHisto(igetter, dttpgthmapbxf[0], nName, dName);
  nName = "DT_TPG_theta_map_bx0";
  makeRatioHisto(igetter, dttpgthmapbxf[1], nName, dName);
  nName = "DT_TPG_theta_map_bx+1";
  makeRatioHisto(igetter, dttpgthmapbxf[2], nName, dName);
}
 
void L1TDTTPGClient::makeRatioHisto(DQMStore::IGetter &igetter, MonitorElement *ratioME, string &nName, string &dName) {
  TH2F *numerator;
  TH2F *denominator;
 
  denominator = this->get2DHisto(input_dir_ + "/" + dName, igetter);
  numerator = this->get2DHisto(input_dir_ + "/" + nName, igetter);
 
  if (numerator && denominator) {
    TH2F *ratio = ratioME->getTH2F();
    if (ratio) {
      ratio->Divide(numerator, denominator);
    } else {
      LogInfo("TriggerDQM") << "[TriggerDQM]: ratio histo named \"" << ratioME->getName() << "\" not found!" << endl;
    }
  } else {
    if (!numerator)
      LogInfo("TriggerDQM") << "[TriggerDQM]: numerator histo \"" << nName << "\" not found!" << endl;
    if (!denominator)
      LogInfo("TriggerDQM") << "[TriggerDQM]: denominator histo \"" << dName << "\" not found!" << endl;
  }
}
 
TH1F *L1TDTTPGClient::get1DHisto(string meName, DQMStore::IGetter &igetter) {
  MonitorElement *me_ = igetter.get(meName);
 
  if (!me_) {
    LogInfo("TriggerDQM") << "ME NOT FOUND.";
    return nullptr;
  }
 
  return me_->getTH1F();
}
 
TH2F *L1TDTTPGClient::get2DHisto(string meName, DQMStore::IGetter &igetter) {
  MonitorElement *me_ = igetter.get(meName);
 
  if (!me_) {
    LogInfo("TriggerDQM") << "ME NOT FOUND.";
    return nullptr;
  }
 
  return me_->getTH2F();
}
 
TProfile2D *L1TDTTPGClient::get2DProfile(string meName, DQMStore::IGetter &igetter) {
  MonitorElement *me_ = igetter.get(meName);
 
  if (!me_) {
    LogInfo("TriggerDQM") << "ME NOT FOUND.";
    return nullptr;
  }
 
  return me_->getTProfile2D();
}
 
TProfile *L1TDTTPGClient::get1DProfile(string meName, DQMStore::IGetter &igetter) {
  MonitorElement *me_ = igetter.get(meName);
 
  if (!me_) {
    LogInfo("TriggerDQM") << "ME NOT FOUND.";
    return nullptr;
  }
 
  return me_->getTProfile();
}
 
void L1TDTTPGClient::setMapPhLabel(MonitorElement *me) {
  me->setAxisTitle("DTTF Sector", 2);
  for (int i = 0; i < 5; i++) {
    ostringstream wheel;
    wheel << i - 2;
    me->setBinLabel(1 + i * 4, "Wheel " + wheel.str(), 1);
  }
}
 
void L1TDTTPGClient::setMapThLabel(MonitorElement *me) {
  me->setAxisTitle("DTTF Sector", 2);
  for (int i = 0; i < 5; i++) {
    ostringstream wheel;
    wheel << i - 2;
    me->setBinLabel(1 + i * 3, "Wheel " + wheel.str(), 1);
  }
}