L1TGCTClient.cc

Go to the documentation of this file.

#include "DQM/L1TMonitorClient/interface/L1TGCTClient.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/QReport.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "DQMServices/Core/interface/MonitorElement.h"

using namespace edm;
using namespace std;

// Define statics for bins etc.
const unsigned int ETABINS = 22;
const float ETAMIN = -0.5;
const float ETAMAX = 21.5;

const unsigned int PHIBINS = 18;
const float PHIMIN = -0.5;
const float PHIMAX = 17.5;

L1TGCTClient::L1TGCTClient(const edm::ParameterSet& ps):
  monitorDir_(ps.getUntrackedParameter<string>("monitorDir","")),
  counterLS_(0), 
  counterEvt_(0), 
  prescaleLS_(ps.getUntrackedParameter<int>("prescaleLS", -1)),
  prescaleEvt_(ps.getUntrackedParameter<int>("prescaleEvt", -1)),
  m_runInEventLoop(ps.getUntrackedParameter<bool>("runInEventLoop", false)),
  m_runInEndLumi(ps.getUntrackedParameter<bool>("runInEndLumi", false)),
  m_runInEndRun(ps.getUntrackedParameter<bool>("runInEndRun", false)),
  m_runInEndJob(ps.getUntrackedParameter<bool>("runInEndJob", false)),
  m_stage1_layer2_(ps.getUntrackedParameter<bool>("stage1_layer2_"))

{
}

L1TGCTClient::~L1TGCTClient(){}

void L1TGCTClient::book(DQMStore::IBooker &ibooker)
{
  // Set to directory with ME in
  ibooker.setCurrentFolder(monitorDir_);
  
  l1GctIsoEmOccEta_ = ibooker.book1D("IsoEmOccEta","ISO EM  #eta OCCUPANCY", ETABINS, ETAMIN, ETAMAX);
  l1GctIsoEmOccPhi_ = ibooker.book1D("IsoEmOccPhi","ISO EM  #phi OCCUPANCY", PHIBINS, PHIMIN, PHIMAX);
  l1GctNonIsoEmOccEta_ = ibooker.book1D("NonIsoEmOccEta","NON-ISO EM  #eta OCCUPANCY", ETABINS, ETAMIN, ETAMAX);
  l1GctNonIsoEmOccPhi_ = ibooker.book1D("NonIsoEmOccPhi","NON-ISO EM  #phi OCCUPANCY", PHIBINS, PHIMIN, PHIMAX);
  l1GctAllJetsOccEta_ = ibooker.book1D("AllJetsOccEta","CENTRAL AND FORWARD JET  #eta OCCUPANCY", ETABINS, ETAMIN, ETAMAX);
  l1GctAllJetsOccPhi_ = ibooker.book1D("AllJetsOccPhi","CENTRAL AND FORWARD JET  #phi OCCUPANCY", PHIBINS, PHIMIN, PHIMAX);
  l1GctCenJetsOccEta_ = ibooker.book1D("CenJetsOccEta","CENTRAL JET  #eta OCCUPANCY", ETABINS, ETAMIN, ETAMAX);
  l1GctCenJetsOccPhi_ = ibooker.book1D("CenJetsOccPhi","CENTRAL JET  #phi OCCUPANCY", PHIBINS, PHIMIN, PHIMAX);
  l1GctForJetsOccEta_ = ibooker.book1D("ForJetsOccEta","FORWARD JET  #eta OCCUPANCY", ETABINS, ETAMIN, ETAMAX);
  l1GctForJetsOccPhi_ = ibooker.book1D("ForJetsOccPhi","FORWARD JET  #phi OCCUPANCY", PHIBINS, PHIMIN, PHIMAX);
  l1GctTauJetsOccEta_ = ibooker.book1D("TauJetsOccEta","TAU JET  #eta OCCUPANCY", ETABINS, ETAMIN, ETAMAX);
  l1GctTauJetsOccPhi_ = ibooker.book1D("TauJetsOccPhi","TAU JET  #phi OCCUPANCY", PHIBINS, PHIMIN, PHIMAX);
  if (m_stage1_layer2_ == true){
    l1GctIsoTauJetsOccEta_ = ibooker.book1D("IsoTauJetsOccEta", "TAU JET #eta OCCUPANCY", ETABINS, ETAMIN, ETAMAX);
    l1GctIsoTauJetsOccPhi_ = ibooker.book1D("IsoTauJetsOccPhi", "TAU JET #phi OCCUPANCY", PHIBINS, PHIMIN, PHIMAX);
  }
    
}

void L1TGCTClient::dqmEndLuminosityBlock(DQMStore::IBooker &ibooker,DQMStore::IGetter &igetter,const edm::LuminosityBlock& lumiSeg, const edm::EventSetup& c) 
{

    if (m_runInEndLumi) {
      book(ibooker);
      processHistograms(igetter);
    }

}

void L1TGCTClient::dqmEndJob(DQMStore::IBooker &ibooker, DQMStore::IGetter &igetter) {

    if (m_runInEndRun) {
      book(ibooker);
      processHistograms(igetter);
    }

}

void L1TGCTClient::processHistograms(DQMStore::IGetter &igetter) {

    MonitorElement* Input;

    if(m_stage1_layer2_ ==false){
      InputDir = "L1T/L1TGCT/";
    } else{
      InputDir = "L1T/L1TStage1Layer2/";
    }

    Input = igetter.get(InputDir + "IsoEmOccEtaPhi");
    if (Input!=nullptr){
      makeXProjection(Input->getTH2F(),l1GctIsoEmOccEta_);
      makeYProjection(Input->getTH2F(),l1GctIsoEmOccPhi_);
    }

    Input = igetter.get(InputDir + "NonIsoEmOccEtaPhi");
    if (Input!=nullptr){
      makeXProjection(Input->getTH2F(),l1GctNonIsoEmOccEta_);
      makeYProjection(Input->getTH2F(),l1GctNonIsoEmOccPhi_);
    }

    Input = igetter.get(InputDir + "AllJetsOccEtaPhi");
    if (Input!=nullptr){
      makeXProjection(Input->getTH2F(),l1GctAllJetsOccEta_);
      makeYProjection(Input->getTH2F(),l1GctAllJetsOccPhi_);
    }

    Input = igetter.get(InputDir + "CenJetsOccEtaPhi");
    if (Input!=nullptr){
      makeXProjection(Input->getTH2F(),l1GctCenJetsOccEta_);
      makeYProjection(Input->getTH2F(),l1GctCenJetsOccPhi_);
    }

    Input = igetter.get(InputDir + "ForJetsOccEtaPhi");
    if (Input!=nullptr){
      makeXProjection(Input->getTH2F(),l1GctForJetsOccEta_);
      makeYProjection(Input->getTH2F(),l1GctForJetsOccPhi_);
    }

    Input = igetter.get(InputDir + "TauJetsOccEtaPhi");
    if (Input!=nullptr){
      makeXProjection(Input->getTH2F(),l1GctTauJetsOccEta_);
      makeYProjection(Input->getTH2F(),l1GctTauJetsOccPhi_);
    }

    if (m_stage1_layer2_ == true){
      Input = igetter.get(InputDir + "IsoTauJetsOccEtaPhi");
      if (Input!=nullptr){
        makeXProjection(Input->getTH2F(),l1GctIsoTauJetsOccEta_);
        makeYProjection(Input->getTH2F(),l1GctIsoTauJetsOccPhi_);
      }      
    }
}


void L1TGCTClient::makeXProjection(TH2F* input, MonitorElement* output)
{
  // Are the provided input and output consistent
  if (input->GetNbinsX() != output->getNbinsX()) return;
  
  // Make the projection
  TH1D* projX = input->ProjectionX();
  
  for (Int_t i=0; i<projX->GetNbinsX(); i++) {
    output->setBinContent(i+1,projX->GetBinContent(i+1));
  }
  delete projX;
}

void L1TGCTClient::makeYProjection(TH2F* input, MonitorElement* output)
{
  // Are the provided input and output consistent
  if (input->GetNbinsY() != output->getNbinsX()) return;
  
  // Make the projection
  TH1D* projY = input->ProjectionY();
  
  for (Int_t i=0; i<projY->GetNbinsX(); i++) {
    output->setBinContent(i+1,projY->GetBinContent(i+1));
  }
  delete projY;
}