RPCRechitProbabilityClient.cc

Go to the documentation of this file.

// Original Author:  Anna Cimmino

#include "DQM/RPCMonitorClient/interface/RPCRecHitProbabilityClient.h"
//Framework
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
//DQMServices
#include "DQMServices/Core/interface/DQMStore.h"

#include <string>

RPCRecHitProbabilityClient::RPCRecHitProbabilityClient(const edm::ParameterSet &iConfig) {
  edm::LogVerbatim("rpcdqmclient") << "[RPCRecHitProbabilityClient]: Constructor";

  std::string subsystemFolder = iConfig.getUntrackedParameter<std::string>("RPCFolder", "RPC");
  std::string recHitTypeFolder = iConfig.getUntrackedParameter<std::string>("MuonFolder", "Muon");

  std::string summaryFolder = iConfig.getUntrackedParameter<std::string>("GlobalFolder", "SummaryHistograms");

  globalFolder_ = subsystemFolder + "/" + recHitTypeFolder + "/" + summaryFolder;
}

RPCRecHitProbabilityClient::~RPCRecHitProbabilityClient() {}

void RPCRecHitProbabilityClient::beginJob() {
  edm::LogVerbatim("rpcrechitprobabilityclient") << "[RPCRecHitProbabilityClient]: Begin Job";
}

void RPCRecHitProbabilityClient::dqmEndLuminosityBlock(DQMStore::IBooker &,
                                                       DQMStore::IGetter &,
                                                       edm::LuminosityBlock const &,
                                                       edm::EventSetup const &) {}

void RPCRecHitProbabilityClient::dqmEndJob(DQMStore::IBooker &ibooker, DQMStore::IGetter &igetter) {
  edm::LogVerbatim("rpcrechitprobabilityclient") << "[RPCRecHitProbabilityClient]: End Run";

  MonitorElement *NumberOfMuonEta = igetter.get(globalFolder_ + "/NumberOfMuonEta");
  MonitorElement *NumberOfMuonPt_B = igetter.get(globalFolder_ + "/NumberOfMuonPt_Barrel");
  MonitorElement *NumberOfMuonPt_EP = igetter.get(globalFolder_ + "/NumberOfMuonPt_EndcapP");
  MonitorElement *NumberOfMuonPt_EM = igetter.get(globalFolder_ + "/NumberOfMuonPt_EndcapM");
  MonitorElement *NumberOfMuonPhi_B = igetter.get(globalFolder_ + "/NumberOfMuonPhi_Barrel");
  MonitorElement *NumberOfMuonPhi_EP = igetter.get(globalFolder_ + "/NumberOfMuonPhi_EndcapP");
  MonitorElement *NumberOfMuonPhi_EM = igetter.get(globalFolder_ + "/NumberOfMuonPhi_EndcapM");

  if (NumberOfMuonEta == nullptr || NumberOfMuonPt_B == nullptr || NumberOfMuonPt_EP == nullptr ||
      NumberOfMuonPt_EM == nullptr || NumberOfMuonPhi_B == nullptr || NumberOfMuonPhi_EP == nullptr ||
      NumberOfMuonPhi_EM == nullptr)
    return;

  TH1F *NumberOfMuonEtaTH1F = NumberOfMuonEta->getTH1F();
  TH1F *NumberOfMuonPtBTH1F = NumberOfMuonPt_B->getTH1F();
  TH1F *NumberOfMuonPtEPTH1F = NumberOfMuonPt_EP->getTH1F();
  TH1F *NumberOfMuonPtEMTH1F = NumberOfMuonPt_EM->getTH1F();
  TH1F *NumberOfMuonPhiBTH1F = NumberOfMuonPhi_B->getTH1F();
  TH1F *NumberOfMuonPhiEPTH1F = NumberOfMuonPhi_EP->getTH1F();
  TH1F *NumberOfMuonPhiEMTH1F = NumberOfMuonPhi_EM->getTH1F();

  MonitorElement *recHit;
  TH1F *recHitTH1F;
  std::stringstream name;

  for (int i = 1; i <= 6; i++) {
    recHit = nullptr;
    recHitTH1F = nullptr;

    name.str("");
    name << globalFolder_ << "/" << i << "RecHitMuonEta";
    recHit = igetter.get(name.str());

    if (recHit) {
      recHitTH1F = recHit->getTH1F();
      recHitTH1F->Divide(NumberOfMuonEtaTH1F);
    }

    recHit = nullptr;
    recHitTH1F = nullptr;

    name.str("");
    name << globalFolder_ << "/" << i << "RecHitMuonPtB";
    recHit = igetter.get(name.str());

    if (recHit) {
      recHitTH1F = recHit->getTH1F();
      recHitTH1F->Divide(NumberOfMuonPtBTH1F);
    }

    recHit = nullptr;
    recHitTH1F = nullptr;

    name.str("");
    name << globalFolder_ << "/" << i << "RecHitMuonPhiB";
    recHit = igetter.get(name.str());

    if (recHit) {
      recHitTH1F = recHit->getTH1F();
      recHitTH1F->Divide(NumberOfMuonPhiBTH1F);
    }

    recHit = nullptr;
    recHitTH1F = nullptr;

    name.str("");
    name << globalFolder_ << "/" << i << "RecHitMuonPtEP";
    recHit = igetter.get(name.str());

    if (recHit) {
      recHitTH1F = recHit->getTH1F();
      recHitTH1F->Divide(NumberOfMuonPtEPTH1F);
    }

    recHit = nullptr;
    recHitTH1F = nullptr;

    name.str("");
    name << globalFolder_ << "/" << i << "RecHitMuonPhiEP";
    recHit = igetter.get(name.str());

    if (recHit) {
      recHitTH1F = recHit->getTH1F();
      recHitTH1F->Divide(NumberOfMuonPhiEPTH1F);
    }

    recHit = nullptr;
    recHitTH1F = nullptr;

    name.str("");
    name << globalFolder_ << "/" << i << "RecHitMuonPtEM";
    recHit = igetter.get(name.str());

    if (recHit) {
      recHitTH1F = recHit->getTH1F();
      recHitTH1F->Divide(NumberOfMuonPtEMTH1F);
    }

    recHit = nullptr;
    recHitTH1F = nullptr;

    name.str("");
    name << globalFolder_ << "/" << i << "RecHitMuonPhiEM";
    recHit = igetter.get(name.str());

    if (recHit) {
      recHitTH1F = recHit->getTH1F();
      recHitTH1F->Divide(NumberOfMuonPhiEMTH1F);
    }
  }
}