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);
  }
}