TopDiLeptonDQM Class Reference

#include <TopDiLeptonDQM.h>

Inheritance diagram for TopDiLeptonDQM:

Public Member Functions
	TopDiLeptonDQM (const edm::ParameterSet &)
	~TopDiLeptonDQM () override
Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final
void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final
void	beginStream (edm::StreamID id) final
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	DQMEDAnalyzer ()
void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	endRun (edm::Run const &run, edm::EventSetup const &setup) final
virtual bool	getCanSaveByLumi ()
Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndRuns () const final

Protected Member Functions
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Private Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

Private Attributes
MonitorElement *	D_eta_elecs_
MonitorElement *	D_eta_lepts_
MonitorElement *	D_eta_muons_
MonitorElement *	D_phi_elecs_
MonitorElement *	D_phi_lepts_
MonitorElement *	D_phi_muons_
MonitorElement *	dimassRC_
MonitorElement *	dimassRC_LOG10_
MonitorElement *	dimassRC_LOGX_
MonitorElement *	dimassWC_
MonitorElement *	dimassWC_LOG10_
MonitorElement *	dimassWC_LOGX_
float	Eff [100]
double	elec_emf_cut_
double	elec_eta_cut_
double	elec_iso_cut_
double	elec_pT_cut_
MonitorElement *	ElecIso_cal_
MonitorElement *	ElecIso_CombRelIso_
MonitorElement *	ElecIso_trk_
edm::EDGetTokenT< reco::GsfElectronCollection >	elecs_
MonitorElement *	eta_elecs_
MonitorElement *	eta_muons_
MonitorElement *	Events_
std::vector< std::string >	hltPaths_
std::vector< std::string >	hltPaths_sig_
std::vector< std::string >	hltPaths_trig_
MonitorElement *	HoverE_elecs_
double	MassWindow_down_
double	MassWindow_up_
std::string	moduleName_
MonitorElement *	MuIso_CombRelIso03_
MonitorElement *	MuIso_emEt03_
MonitorElement *	MuIso_hadEt03_
MonitorElement *	MuIso_hoEt03_
MonitorElement *	MuIso_nJets03_
MonitorElement *	MuIso_nTracks03_
MonitorElement *	MuIso_sumPt03_
double	muon_eta_cut_
double	muon_iso_cut_
double	muon_pT_cut_
edm::EDGetTokenT< reco::MuonCollection >	muons_
int	N_elel
int	N_muel
int	N_mumu
int	N_sig [100]
int	N_trig [100]
MonitorElement *	Nelecs_
MonitorElement *	Nelecs_charge_
MonitorElement *	Nelecs_iso_
MonitorElement *	Nmuons_
MonitorElement *	Nmuons_charge_
MonitorElement *	Nmuons_iso_
MonitorElement *	Ntracks_
std::ofstream	outfile
std::string	outputFile_
MonitorElement *	phi_elecs_
MonitorElement *	phi_muons_
MonitorElement *	pT_elecs_
MonitorElement *	pT_muons_
MonitorElement *	TriggerEff_
edm::EDGetTokenT< edm::TriggerResults >	triggerResults_
MonitorElement *	Trigs_
edm::EDGetTokenT< reco::VertexCollection >	vertex_
double	vertex_X_cut_
double	vertex_Y_cut_
double	vertex_Z_cut_
MonitorElement *	VxVy_muons_
MonitorElement *	Vz_muons_

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore
typedef dqm::reco::MonitorElement	MonitorElement
Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
typedef CacheContexts< T... >	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T... >	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache
Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)
static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)
static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)
static std::unique_ptr< DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)
Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_
unsigned int	streamId_

Detailed Description

Definition at line 38 of file TopDiLeptonDQM.h.

Constructor & Destructor Documentation

TopDiLeptonDQM()

TopDiLeptonDQM::TopDiLeptonDQM ( const edm::ParameterSet &  ps )

explicit

Definition at line 16 of file TopDiLeptonDQM.cc.

                                                        {
  moduleName_ = ps.getUntrackedParameter<string>("moduleName");
  outputFile_ = ps.getUntrackedParameter<string>("outputFile");
  triggerResults_ = consumes<TriggerResults>(ps.getParameter<edm::InputTag>("TriggerResults"));
  hltPaths_ = ps.getParameter<vector<string> >("hltPaths");
  hltPaths_sig_ = ps.getParameter<vector<string> >("hltPaths_sig");
  hltPaths_trig_ = ps.getParameter<vector<string> >("hltPaths_trig");
 
  vertex_ = consumes<reco::VertexCollection>(ps.getParameter<edm::InputTag>("vertexCollection"));
  vertex_X_cut_ = ps.getParameter<double>("vertex_X_cut");
  vertex_Y_cut_ = ps.getParameter<double>("vertex_Y_cut");
  vertex_Z_cut_ = ps.getParameter<double>("vertex_Z_cut");
 
  muons_ = consumes<reco::MuonCollection>(ps.getParameter<edm::InputTag>("muonCollection"));
  muon_pT_cut_ = ps.getParameter<double>("muon_pT_cut");
  muon_eta_cut_ = ps.getParameter<double>("muon_eta_cut");
  muon_iso_cut_ = ps.getParameter<double>("muon_iso_cut");
 
  elecs_ = consumes<reco::GsfElectronCollection>(ps.getParameter<edm::InputTag>("elecCollection"));
  elec_pT_cut_ = ps.getParameter<double>("elec_pT_cut");
  elec_eta_cut_ = ps.getParameter<double>("elec_eta_cut");
  elec_iso_cut_ = ps.getParameter<double>("elec_iso_cut");
  elec_emf_cut_ = ps.getParameter<double>("elec_emf_cut");
 
  MassWindow_up_ = ps.getParameter<double>("MassWindow_up");
  MassWindow_down_ = ps.getParameter<double>("MassWindow_down");
 
  for (int i = 0; i < 100; ++i) {
    N_sig[i] = 0;
    N_trig[i] = 0;
    Eff[i] = 0.;
  }
 
  N_mumu = 0;
  N_muel = 0;
  N_elel = 0;
  Events_ = nullptr;
  Trigs_ = nullptr;
  TriggerEff_ = nullptr;
  Ntracks_ = nullptr;
  Nmuons_ = nullptr;
  Nmuons_iso_ = nullptr;
  Nmuons_charge_ = nullptr;
  VxVy_muons_ = nullptr;
  Vz_muons_ = nullptr;
  pT_muons_ = nullptr;
  eta_muons_ = nullptr;
  phi_muons_ = nullptr;
  Nelecs_ = nullptr;
  Nelecs_iso_ = nullptr;
  Nelecs_charge_ = nullptr;
  HoverE_elecs_ = nullptr;
  pT_elecs_ = nullptr;
  eta_elecs_ = nullptr;
  phi_elecs_ = nullptr;
  MuIso_emEt03_ = nullptr;
  MuIso_hadEt03_ = nullptr;
  MuIso_hoEt03_ = nullptr;
  MuIso_nJets03_ = nullptr;
  MuIso_nTracks03_ = nullptr;
  MuIso_sumPt03_ = nullptr;
  MuIso_CombRelIso03_ = nullptr;
  ElecIso_cal_ = nullptr;
  ElecIso_trk_ = nullptr;
  ElecIso_CombRelIso_ = nullptr;
  dimassRC_ = nullptr;
  dimassWC_ = nullptr;
  dimassRC_LOGX_ = nullptr;
  dimassWC_LOGX_ = nullptr;
  dimassRC_LOG10_ = nullptr;
  dimassWC_LOG10_ = nullptr;
  D_eta_muons_ = nullptr;
  D_phi_muons_ = nullptr;
  D_eta_elecs_ = nullptr;
  D_phi_elecs_ = nullptr;
  D_eta_lepts_ = nullptr;
  D_phi_lepts_ = nullptr;
}

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), mps_fire::i, and commonCuts_cff::triggerResults_.

~TopDiLeptonDQM()

TopDiLeptonDQM::~TopDiLeptonDQM ( )

override

Definition at line 95 of file TopDiLeptonDQM.cc.

{}

Member Function Documentation

analyze()

void TopDiLeptonDQM::analyze ( const edm::Event &  evt, const edm::EventSetup &  context  )

overrideprivatevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 168 of file TopDiLeptonDQM.cc.

                                                                              {
  // ------------------------
  //  Global Event Variables
  // ------------------------
 
  const int N_TriggerPaths = hltPaths_.size();
  const int N_SignalPaths = hltPaths_sig_.size();
  const int N_ControlPaths = hltPaths_trig_.size();
 
  bool Fired_Signal_Trigger[100] = {false};
  bool Fired_Control_Trigger[100] = {false};
 
  int N_leptons = 0;
  int N_iso_mu = 0;
  int N_iso_el = 0;
 
  double DilepMass = 0.;
 
  double vertex_X = 100.;
  double vertex_Y = 100.;
  double vertex_Z = 100.;
 
  // ------------------------
  //  Analyze Primary Vertex
  // ------------------------
 
  edm::Handle<reco::VertexCollection> vertexs;
  evt.getByToken(vertex_, vertexs);
 
  if (!vertexs.failedToGet()) {
    reco::Vertex primaryVertex = vertexs->front();
 
    int numberTracks = primaryVertex.tracksSize();
    //    double ndof      = primaryVertex.ndof();
    bool fake = primaryVertex.isFake();
 
    Ntracks_->Fill(numberTracks);
 
    if (!fake && numberTracks > 3) {
      vertex_X = primaryVertex.x();
      vertex_Y = primaryVertex.y();
      vertex_Z = primaryVertex.z();
    }
  }
 
  // -------------------------
  //  Analyze Trigger Results
  // -------------------------
 
  edm::Handle<TriggerResults> trigResults;
  evt.getByToken(triggerResults_, trigResults);
 
  if (!trigResults.failedToGet()) {
    int N_Triggers = trigResults->size();
 
    const edm::TriggerNames& trigName = evt.triggerNames(*trigResults);
 
    for (int i_Trig = 0; i_Trig < N_Triggers; ++i_Trig) {
      if (trigResults.product()->accept(i_Trig)) {
        // Check for all trigger paths
 
        for (int i = 0; i < N_TriggerPaths; i++) {
          if (trigName.triggerName(i_Trig) == hltPaths_[i]) {
            Trigs_->Fill(i);
            Trigs_->setBinLabel(i + 1, hltPaths_[i], 1);
          }
        }
 
        // Check for signal & control trigger paths
 
        for (int j = 0; j < N_SignalPaths; ++j) {
          if (trigName.triggerName(i_Trig) == hltPaths_sig_[j])
            Fired_Signal_Trigger[j] = true;
        }
 
        for (int k = 0; k < N_ControlPaths; ++k) {
          if (trigName.triggerName(i_Trig) == hltPaths_trig_[k])
            Fired_Control_Trigger[k] = true;
        }
      }
    }
  }
 
  // ------------------------
  //  Analyze Muon Isolation
  // ------------------------
 
  edm::Handle<reco::MuonCollection> muons;
  evt.getByToken(muons_, muons);
 
  reco::MuonCollection::const_iterator muon;
 
  if (!muons.failedToGet()) {
    Nmuons_->Fill(muons->size());
 
    N_leptons = N_leptons + muons->size();
 
    for (muon = muons->begin(); muon != muons->end(); ++muon) {
      float N_muons = muons->size();
      float Q_muon = muon->charge();
 
      Nmuons_charge_->Fill(N_muons * Q_muon);
 
      double track_X = 100.;
      double track_Y = 100.;
      double track_Z = 100.;
 
      if (muon->isGlobalMuon()) {
        reco::TrackRef track = muon->globalTrack();
 
        track_X = track->vx();
        track_Y = track->vy();
        track_Z = track->vz();
 
        VxVy_muons_->Fill(track_X, track_Y);
        Vz_muons_->Fill(track_Z);
      }
 
      // Vertex and kinematic cuts
 
      if (track_X > vertex_X_cut_)
        continue;
      if (track_Y > vertex_Y_cut_)
        continue;
      if (track_Z > vertex_Z_cut_)
        continue;
      if (muon->pt() < muon_pT_cut_)
        continue;
      if (abs(muon->eta()) > muon_eta_cut_)
        continue;
 
      reco::MuonIsolation muIso03 = muon->isolationR03();
 
      double muonCombRelIso = 1.;
 
      if (muon->pt() != 0.)
        muonCombRelIso = (muIso03.emEt + muIso03.hadEt + muIso03.hoEt + muIso03.sumPt) / muon->pt();
 
      MuIso_CombRelIso03_->Fill(muonCombRelIso);
 
      MuIso_emEt03_->Fill(muIso03.emEt);
      MuIso_hadEt03_->Fill(muIso03.hadEt);
      MuIso_hoEt03_->Fill(muIso03.hoEt);
      MuIso_nJets03_->Fill(muIso03.nJets);
      MuIso_nTracks03_->Fill(muIso03.nTracks);
      MuIso_sumPt03_->Fill(muIso03.sumPt);
 
      if (muonCombRelIso < muon_iso_cut_)
        ++N_iso_mu;
    }
 
    Nmuons_iso_->Fill(N_iso_mu);
  }
 
  // ----------------------------
  //  Analyze Electron Isolation
  // ----------------------------
 
  edm::Handle<reco::GsfElectronCollection> elecs;
  evt.getByToken(elecs_, elecs);
 
  reco::GsfElectronCollection::const_iterator elec;
 
  if (!elecs.failedToGet()) {
    Nelecs_->Fill(elecs->size());
 
    N_leptons = N_leptons + elecs->size();
 
    for (elec = elecs->begin(); elec != elecs->end(); ++elec) {
      float N_elecs = elecs->size();
      float Q_elec = elec->charge();
      float HoverE = elec->hcalOverEcal();
 
      HoverE_elecs_->Fill(HoverE);
 
      Nelecs_charge_->Fill(N_elecs * Q_elec);
 
      double track_X = 100.;
      double track_Y = 100.;
      double track_Z = 100.;
 
      reco::GsfTrackRef track = elec->gsfTrack();
 
      track_X = track->vx();
      track_Y = track->vy();
      track_Z = track->vz();
 
      // Vertex and kinematic cuts
 
      if (track_X > vertex_X_cut_)
        continue;
      if (track_Y > vertex_Y_cut_)
        continue;
      if (track_Z > vertex_Z_cut_)
        continue;
      if (elec->pt() < elec_pT_cut_)
        continue;
      if (abs(elec->eta()) > elec_eta_cut_)
        continue;
      if (HoverE > elec_emf_cut_)
        continue;
 
      reco::GsfElectron::IsolationVariables elecIso = elec->dr03IsolationVariables();
 
      double elecCombRelIso = 1.;
 
      if (elec->et() != 0.)
        elecCombRelIso = (elecIso.ecalRecHitSumEt + elecIso.hcalDepth1TowerSumEt + elecIso.tkSumPt) / elec->et();
 
      ElecIso_CombRelIso_->Fill(elecCombRelIso);
 
      ElecIso_cal_->Fill(elecIso.ecalRecHitSumEt);
      ElecIso_trk_->Fill(elecIso.tkSumPt);
 
      if (elecCombRelIso < elec_iso_cut_)
        ++N_iso_el;
    }
 
    Nelecs_iso_->Fill(N_iso_el);
  }
 
  // --------------------
  //  TWO Isolated MUONS
  // --------------------
 
  if (N_iso_mu > 1) {
    // Vertex cut
 
    if (vertex_X < vertex_X_cut_ && vertex_Y < vertex_Y_cut_ && vertex_Z < vertex_Z_cut_) {
      ++N_mumu;
 
      Events_->Fill(1.);
 
      reco::MuonCollection::const_reference mu1 = muons->at(0);
      reco::MuonCollection::const_reference mu2 = muons->at(1);
 
      DilepMass = sqrt((mu1.energy() + mu2.energy()) * (mu1.energy() + mu2.energy()) -
                       (mu1.px() + mu2.px()) * (mu1.px() + mu2.px()) - (mu1.py() + mu2.py()) * (mu1.py() + mu2.py()) -
                       (mu1.pz() + mu2.pz()) * (mu1.pz() + mu2.pz()));
 
      // Opposite muon charges -> Right Charge (RC)
 
      if (mu1.charge() * mu2.charge() < 0.) {
        dimassRC_LOG10_->Fill(log10(DilepMass));
        dimassRC_->Fill(DilepMass);
        dimassRC_LOGX_->Fill(DilepMass);
 
        if (DilepMass > MassWindow_down_ && DilepMass < MassWindow_up_) {
          for (muon = muons->begin(); muon != muons->end(); ++muon) {
            pT_muons_->Fill(muon->pt());
            eta_muons_->Fill(muon->eta());
            phi_muons_->Fill(muon->phi());
          }
 
          D_eta_muons_->Fill(mu1.eta() - mu2.eta());
          D_phi_muons_->Fill(mu1.phi() - mu2.phi());
 
          // Determinating trigger efficiencies
 
          for (int k = 0; k < N_SignalPaths; ++k) {
            if (Fired_Signal_Trigger[k] && Fired_Control_Trigger[k])
              ++N_sig[k];
 
            if (Fired_Control_Trigger[k])
              ++N_trig[k];
 
            if (N_trig[k] != 0)
              Eff[k] = N_sig[k] / static_cast<float>(N_trig[k]);
 
            TriggerEff_->setBinContent(k + 1, Eff[k]);
            TriggerEff_->setBinLabel(k + 1, "#frac{[" + hltPaths_sig_[k] + "]}{vs. [" + hltPaths_trig_[k] + "]}", 1);
          }
        }
      }
 
      // Same muon charges -> Wrong Charge (WC)
 
      if (mu1.charge() * mu2.charge() > 0.) {
        dimassWC_LOG10_->Fill(log10(DilepMass));
        dimassWC_->Fill(DilepMass);
        dimassWC_LOGX_->Fill(DilepMass);
      }
    }
  }
 
  // -----------------------------
  //  TWO Isolated LEPTONS (mu/e)
  // -----------------------------
 
  if (N_iso_el > 0 && N_iso_mu > 0) {
    // Vertex cut
 
    if (vertex_X < vertex_X_cut_ && vertex_Y < vertex_Y_cut_ && vertex_Z < vertex_Z_cut_) {
      ++N_muel;
 
      Events_->Fill(2.);
 
      reco::MuonCollection::const_reference mu1 = muons->at(0);
      reco::GsfElectronCollection::const_reference el1 = elecs->at(0);
 
      DilepMass = sqrt((mu1.energy() + el1.energy()) * (mu1.energy() + el1.energy()) -
                       (mu1.px() + el1.px()) * (mu1.px() + el1.px()) - (mu1.py() + el1.py()) * (mu1.py() + el1.py()) -
                       (mu1.pz() + el1.pz()) * (mu1.pz() + el1.pz()));
 
      // Opposite lepton charges -> Right Charge (RC)
 
      if (mu1.charge() * el1.charge() < 0.) {
        dimassRC_LOG10_->Fill(log10(DilepMass));
        dimassRC_->Fill(DilepMass);
        dimassRC_LOGX_->Fill(DilepMass);
 
        if (DilepMass > MassWindow_down_ && DilepMass < MassWindow_up_) {
          for (muon = muons->begin(); muon != muons->end(); ++muon) {
            pT_muons_->Fill(muon->pt());
            eta_muons_->Fill(muon->eta());
            phi_muons_->Fill(muon->phi());
          }
 
          for (elec = elecs->begin(); elec != elecs->end(); ++elec) {
            pT_elecs_->Fill(elec->pt());
            eta_elecs_->Fill(elec->eta());
            phi_elecs_->Fill(elec->phi());
          }
 
          D_eta_lepts_->Fill(mu1.eta() - el1.eta());
          D_phi_lepts_->Fill(mu1.phi() - el1.phi());
 
          // Determinating trigger efficiencies
 
          for (int k = 0; k < N_SignalPaths; ++k) {
            if (Fired_Signal_Trigger[k] && Fired_Control_Trigger[k])
              ++N_sig[k];
 
            if (Fired_Control_Trigger[k])
              ++N_trig[k];
 
            if (N_trig[k] != 0)
              Eff[k] = N_sig[k] / static_cast<float>(N_trig[k]);
 
            TriggerEff_->setBinContent(k + 1, Eff[k]);
            TriggerEff_->setBinLabel(k + 1, "#frac{[" + hltPaths_sig_[k] + "]}{vs. [" + hltPaths_trig_[k] + "]}", 1);
          }
        }
      }
 
      // Same muon charges -> Wrong Charge (WC)
 
      if (mu1.charge() * el1.charge() > 0.) {
        dimassWC_LOG10_->Fill(log10(DilepMass));
        dimassWC_->Fill(DilepMass);
        dimassWC_LOGX_->Fill(DilepMass);
      }
    }
  }
 
  // ------------------------
  //  TWO Isolated ELECTRONS
  // ------------------------
 
  if (N_iso_el > 1) {
    // Vertex cut
 
    if (vertex_X < vertex_X_cut_ && vertex_Y < vertex_Y_cut_ && vertex_Z < vertex_Z_cut_) {
      ++N_elel;
 
      Events_->Fill(3.);
 
      reco::GsfElectronCollection::const_reference el1 = elecs->at(0);
      reco::GsfElectronCollection::const_reference el2 = elecs->at(1);
 
      DilepMass = sqrt((el1.energy() + el2.energy()) * (el1.energy() + el2.energy()) -
                       (el1.px() + el2.px()) * (el1.px() + el2.px()) - (el1.py() + el2.py()) * (el1.py() + el2.py()) -
                       (el1.pz() + el2.pz()) * (el1.pz() + el2.pz()));
 
      // Opposite lepton charges -> Right Charge (RC)
 
      if (el1.charge() * el2.charge() < 0.) {
        dimassRC_LOG10_->Fill(log10(DilepMass));
        dimassRC_->Fill(DilepMass);
        dimassRC_LOGX_->Fill(DilepMass);
 
        if (DilepMass > MassWindow_down_ && DilepMass < MassWindow_up_) {
          for (elec = elecs->begin(); elec != elecs->end(); ++elec) {
            pT_elecs_->Fill(elec->pt());
            eta_elecs_->Fill(elec->eta());
            phi_elecs_->Fill(elec->phi());
          }
 
          D_eta_elecs_->Fill(el1.eta() - el2.eta());
          D_phi_elecs_->Fill(el1.phi() - el2.phi());
 
          // Determinating trigger efficiencies
 
          for (int k = 0; k < N_SignalPaths; ++k) {
            if (Fired_Signal_Trigger[k] && Fired_Control_Trigger[k])
              ++N_sig[k];
 
            if (Fired_Control_Trigger[k])
              ++N_trig[k];
 
            if (N_trig[k] != 0)
              Eff[k] = N_sig[k] / static_cast<float>(N_trig[k]);
 
            TriggerEff_->setBinContent(k + 1, Eff[k]);
            TriggerEff_->setBinLabel(k + 1, "#frac{[" + hltPaths_sig_[k] + "]}{vs. [" + hltPaths_trig_[k] + "]}", 1);
          }
        }
      }
 
      // Same muon charges -> Wrong Charge (WC)
 
      if (el1.charge() * el2.charge() > 0.) {
        dimassWC_LOG10_->Fill(log10(DilepMass));
        dimassWC_->Fill(DilepMass);
        dimassWC_LOGX_->Fill(DilepMass);
      }
    }
  }
}

References funct::abs(), edm::HLTGlobalStatus::accept(), reco::GsfElectron::IsolationVariables::ecalRecHitSumEt, singleTopDQM_cfi::elecs, reco::MuonIsolation::emEt, edm::HandleBase::failedToGet(), edm::Event::getByToken(), reco::MuonIsolation::hadEt, reco::GsfElectron::IsolationVariables::hcalDepth1TowerSumEt, reco::MuonIsolation::hoEt, mps_fire::i, dqmiolumiharvest::j, dqmdumpme::k, HLT_2018_cff::muon, PDWG_BPHSkim_cff::muons, reco::MuonIsolation::nJets, reco::MuonIsolation::nTracks, BeamMonitor_cff::primaryVertex, edm::Handle< T >::product(), edm::HLTGlobalStatus::size(), mathSSE::sqrt(), reco::MuonIsolation::sumPt, reco::GsfElectron::IsolationVariables::tkSumPt, HLT_2018_cff::track, edm::Event::triggerNames(), commonCuts_cff::triggerResults_, and EgHLTOffTrigSelection_cfi::trigName.

bookHistograms()

void TopDiLeptonDQM::bookHistograms ( DQMStore::IBooker &  iBooker, edm::Run const &  , edm::EventSetup const &    )

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 97 of file TopDiLeptonDQM.cc.

                                                                                                 {
  iBooker.setCurrentFolder(moduleName_);
 
  Events_ = iBooker.book1D("00_Events", "Isolated dilepton events", 5, 0., 5.);
  Events_->setBinLabel(2, "#mu #mu", 1);
  Events_->setBinLabel(3, "#mu e", 1);
  Events_->setBinLabel(4, "e e", 1);
 
  Trigs_ = iBooker.book1D("01_Trigs", "Fired muon/electron triggers", 15, 0., 15.);
  TriggerEff_ = iBooker.book1D("02_TriggerEff", "HL Trigger Efficiencies", 10, 0., 10.);
  TriggerEff_->setTitle(
      "HL Trigger Efficiencies #epsilon_{signal} = #frac{[signal] && "
      "[control]}{[control]}");
  Ntracks_ = iBooker.book1D("Ntracks", "Number of tracks", 50, 0., 50.);
 
  Nmuons_ = iBooker.book1D("03_Nmuons", "Number of muons", 20, 0., 10.);
  Nmuons_iso_ = iBooker.book1D("04_Nmuons_iso", "Number of isolated muons", 20, 0., 10.);
  Nmuons_charge_ = iBooker.book1D("Nmuons_charge", "Number of muons * moun charge", 19, -10., 10.);
  VxVy_muons_ = iBooker.book2D("VxVy_muons", "Vertex x-y-positon (global)", 40, -1., 1., 40, -1., 1.);
  Vz_muons_ = iBooker.book1D("Vz_muons", "Vertex z-positon (global)", 40, -20., 20.);
  pT_muons_ = iBooker.book1D("pT_muons", "P_T of muons", 40, 0., 200.);
  eta_muons_ = iBooker.book1D("eta_muons", "Eta of muons", 50, -5., 5.);
  phi_muons_ = iBooker.book1D("phi_muons", "Phi of muons", 40, -4., 4.);
 
  Nelecs_ = iBooker.book1D("05_Nelecs", "Number of electrons", 20, 0., 10.);
  Nelecs_iso_ = iBooker.book1D("06_Nelecs_iso", "Number of isolated electrons", 20, 0., 10.);
  Nelecs_charge_ = iBooker.book1D("Nelecs_charge", "Number of elecs * elec charge", 19, -10., 10.);
  HoverE_elecs_ = iBooker.book1D("HoverE_elecs", "Hadronic over Ecal energy", 50, 0., 1.);
  pT_elecs_ = iBooker.book1D("pT_elecs", "P_T of electrons", 40, 0., 200.);
  eta_elecs_ = iBooker.book1D("eta_elecs", "Eta of electrons", 50, -5., 5.);
  phi_elecs_ = iBooker.book1D("phi_elecs", "Phi of electrons", 40, -4., 4.);
 
  MuIso_emEt03_ = iBooker.book1D("MuIso_emEt03", "Muon emEt03", 20, 0., 20.);
  MuIso_hadEt03_ = iBooker.book1D("MuIso_hadEt03", "Muon hadEt03", 20, 0., 20.);
  MuIso_hoEt03_ = iBooker.book1D("MuIso_hoEt03", "Muon hoEt03", 20, 0., 20.);
  MuIso_nJets03_ = iBooker.book1D("MuIso_nJets03", "Muon nJets03", 10, 0., 10.);
  MuIso_nTracks03_ = iBooker.book1D("MuIso_nTracks03", "Muon nTracks03", 20, 0., 20.);
  MuIso_sumPt03_ = iBooker.book1D("MuIso_sumPt03", "Muon sumPt03", 20, 0., 40.);
  MuIso_CombRelIso03_ = iBooker.book1D("07_MuIso_CombRelIso03", "Muon CombRelIso03", 20, 0., 1.);
 
  ElecIso_cal_ = iBooker.book1D("ElecIso_cal", "Electron Iso_cal", 21, -1., 20.);
  ElecIso_trk_ = iBooker.book1D("ElecIso_trk", "Electron Iso_trk", 21, -2., 40.);
  ElecIso_CombRelIso_ = iBooker.book1D("08_ElecIso_CombRelIso", "Electron CombRelIso", 20, 0., 1.);
 
  const int nbins = 200;
 
  double logmin = 0.;
  double logmax = 3.;  // 10^(3.)=1000
 
  float bins[nbins + 1];
 
  for (int i = 0; i <= nbins; i++) {
    double log = logmin + (logmax - logmin) * i / nbins;
    bins[i] = std::pow(10.0, log);
  }
 
  dimassRC_ = iBooker.book1D("09_dimassRC", "Dilepton mass RC", 50, 0., 200.);
  dimassWC_ = iBooker.book1D("11_dimassWC", "Dilepton mass WC", 50, 0., 200.);
  dimassRC_LOGX_ = iBooker.book1D("10_dimassRC_LOGX", "Dilepton mass RC LOG", nbins, &bins[0]);
  dimassWC_LOGX_ = iBooker.book1D("12_dimassWC_LOGX", "Dilepton mass WC LOG", nbins, &bins[0]);
  dimassRC_LOG10_ = iBooker.book1D("dimassRC_LOG10", "Dilepton mass RC LOG", 50, 0., 2.5);
  dimassWC_LOG10_ = iBooker.book1D("dimassWC_LOG10", "Dilepton mass WC LOG", 50, 0., 2.5);
 
  D_eta_muons_ = iBooker.book1D("13_D_eta_muons", "#Delta eta_muons", 20, -5., 5.);
  D_phi_muons_ = iBooker.book1D("14_D_phi_muons", "#Delta phi_muons", 20, -5., 5.);
  D_eta_elecs_ = iBooker.book1D("D_eta_elecs", "#Delta eta_elecs", 20, -5., 5.);
  D_phi_elecs_ = iBooker.book1D("D_phi_elecs", "#Delta phi_elecs", 20, -5., 5.);
  D_eta_lepts_ = iBooker.book1D("D_eta_lepts", "#Delta eta_lepts", 20, -5., 5.);
  D_phi_lepts_ = iBooker.book1D("D_phi_lepts", "#Delta phi_lepts", 20, -5., 5.);
}

References trigObjTnPSource_cfi::bins, dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), mps_fire::i, dqm-mbProfile::log, LaserClient_cfi::nbins, funct::pow(), dqm::impl::MonitorElement::setBinLabel(), dqm::implementation::NavigatorBase::setCurrentFolder(), and dqm::impl::MonitorElement::setTitle().

Member Data Documentation

D_eta_elecs_

MonitorElement* TopDiLeptonDQM::D_eta_elecs_

private

Definition at line 128 of file TopDiLeptonDQM.h.

D_eta_lepts_

MonitorElement* TopDiLeptonDQM::D_eta_lepts_

private

Definition at line 130 of file TopDiLeptonDQM.h.

D_eta_muons_

MonitorElement* TopDiLeptonDQM::D_eta_muons_

private

Definition at line 126 of file TopDiLeptonDQM.h.

D_phi_elecs_

MonitorElement* TopDiLeptonDQM::D_phi_elecs_

private

Definition at line 129 of file TopDiLeptonDQM.h.

D_phi_lepts_

MonitorElement* TopDiLeptonDQM::D_phi_lepts_

private

Definition at line 131 of file TopDiLeptonDQM.h.

D_phi_muons_

MonitorElement* TopDiLeptonDQM::D_phi_muons_

private

Definition at line 127 of file TopDiLeptonDQM.h.

dimassRC_

MonitorElement* TopDiLeptonDQM::dimassRC_

private

Definition at line 119 of file TopDiLeptonDQM.h.

dimassRC_LOG10_

MonitorElement* TopDiLeptonDQM::dimassRC_LOG10_

private

Definition at line 123 of file TopDiLeptonDQM.h.

dimassRC_LOGX_

MonitorElement* TopDiLeptonDQM::dimassRC_LOGX_

private

Definition at line 121 of file TopDiLeptonDQM.h.

dimassWC_

MonitorElement* TopDiLeptonDQM::dimassWC_

private

Definition at line 120 of file TopDiLeptonDQM.h.

dimassWC_LOG10_

MonitorElement* TopDiLeptonDQM::dimassWC_LOG10_

private

Definition at line 124 of file TopDiLeptonDQM.h.

dimassWC_LOGX_

MonitorElement* TopDiLeptonDQM::dimassWC_LOGX_

private

Definition at line 122 of file TopDiLeptonDQM.h.

Eff

float TopDiLeptonDQM::Eff[100]

private

Definition at line 60 of file TopDiLeptonDQM.h.

elec_emf_cut_

double TopDiLeptonDQM::elec_emf_cut_

private

Definition at line 80 of file TopDiLeptonDQM.h.

elec_eta_cut_

double TopDiLeptonDQM::elec_eta_cut_

private

Definition at line 78 of file TopDiLeptonDQM.h.

elec_iso_cut_

double TopDiLeptonDQM::elec_iso_cut_

private

Definition at line 79 of file TopDiLeptonDQM.h.

elec_pT_cut_

double TopDiLeptonDQM::elec_pT_cut_

private

Definition at line 77 of file TopDiLeptonDQM.h.

ElecIso_cal_

MonitorElement* TopDiLeptonDQM::ElecIso_cal_

private

Definition at line 115 of file TopDiLeptonDQM.h.

ElecIso_CombRelIso_

MonitorElement* TopDiLeptonDQM::ElecIso_CombRelIso_

private

Definition at line 117 of file TopDiLeptonDQM.h.

ElecIso_trk_

MonitorElement* TopDiLeptonDQM::ElecIso_trk_

private

Definition at line 116 of file TopDiLeptonDQM.h.

elecs_

edm::EDGetTokenT<reco::GsfElectronCollection> TopDiLeptonDQM::elecs_

private

Definition at line 76 of file TopDiLeptonDQM.h.

eta_elecs_

MonitorElement* TopDiLeptonDQM::eta_elecs_

private

Definition at line 104 of file TopDiLeptonDQM.h.

eta_muons_

MonitorElement* TopDiLeptonDQM::eta_muons_

private

Definition at line 96 of file TopDiLeptonDQM.h.

Events_

MonitorElement* TopDiLeptonDQM::Events_

private

Definition at line 85 of file TopDiLeptonDQM.h.

hltPaths_

std::vector<std::string> TopDiLeptonDQM::hltPaths_

private

Definition at line 52 of file TopDiLeptonDQM.h.

hltPaths_sig_

std::vector<std::string> TopDiLeptonDQM::hltPaths_sig_

private

Definition at line 53 of file TopDiLeptonDQM.h.

hltPaths_trig_

std::vector<std::string> TopDiLeptonDQM::hltPaths_trig_

private

Definition at line 54 of file TopDiLeptonDQM.h.

HoverE_elecs_

MonitorElement* TopDiLeptonDQM::HoverE_elecs_

private

Definition at line 102 of file TopDiLeptonDQM.h.

MassWindow_down_

double TopDiLeptonDQM::MassWindow_down_

private

Definition at line 83 of file TopDiLeptonDQM.h.

MassWindow_up_

double TopDiLeptonDQM::MassWindow_up_

private

Definition at line 82 of file TopDiLeptonDQM.h.

moduleName_

std::string TopDiLeptonDQM::moduleName_

private

Definition at line 49 of file TopDiLeptonDQM.h.

MuIso_CombRelIso03_

MonitorElement* TopDiLeptonDQM::MuIso_CombRelIso03_

private

Definition at line 113 of file TopDiLeptonDQM.h.

MuIso_emEt03_

MonitorElement* TopDiLeptonDQM::MuIso_emEt03_

private

Definition at line 107 of file TopDiLeptonDQM.h.

MuIso_hadEt03_

MonitorElement* TopDiLeptonDQM::MuIso_hadEt03_

private

Definition at line 108 of file TopDiLeptonDQM.h.

MuIso_hoEt03_

MonitorElement* TopDiLeptonDQM::MuIso_hoEt03_

private

Definition at line 109 of file TopDiLeptonDQM.h.

MuIso_nJets03_

MonitorElement* TopDiLeptonDQM::MuIso_nJets03_

private

Definition at line 110 of file TopDiLeptonDQM.h.

MuIso_nTracks03_

MonitorElement* TopDiLeptonDQM::MuIso_nTracks03_

private

Definition at line 111 of file TopDiLeptonDQM.h.

MuIso_sumPt03_

MonitorElement* TopDiLeptonDQM::MuIso_sumPt03_

private

Definition at line 112 of file TopDiLeptonDQM.h.

muon_eta_cut_

double TopDiLeptonDQM::muon_eta_cut_

private

Definition at line 73 of file TopDiLeptonDQM.h.

muon_iso_cut_

double TopDiLeptonDQM::muon_iso_cut_

private

Definition at line 74 of file TopDiLeptonDQM.h.

muon_pT_cut_

double TopDiLeptonDQM::muon_pT_cut_

private

Definition at line 72 of file TopDiLeptonDQM.h.

muons_

edm::EDGetTokenT<reco::MuonCollection> TopDiLeptonDQM::muons_

private

Definition at line 71 of file TopDiLeptonDQM.h.

N_elel

int TopDiLeptonDQM::N_elel

private

Definition at line 64 of file TopDiLeptonDQM.h.

N_muel

int TopDiLeptonDQM::N_muel

private

Definition at line 63 of file TopDiLeptonDQM.h.

N_mumu

int TopDiLeptonDQM::N_mumu

private

Definition at line 62 of file TopDiLeptonDQM.h.

N_sig

int TopDiLeptonDQM::N_sig[100]

private

Definition at line 58 of file TopDiLeptonDQM.h.

N_trig

int TopDiLeptonDQM::N_trig[100]

private

Definition at line 59 of file TopDiLeptonDQM.h.

Nelecs_

MonitorElement* TopDiLeptonDQM::Nelecs_

private

Definition at line 99 of file TopDiLeptonDQM.h.

Nelecs_charge_

MonitorElement* TopDiLeptonDQM::Nelecs_charge_

private

Definition at line 101 of file TopDiLeptonDQM.h.

Nelecs_iso_

MonitorElement* TopDiLeptonDQM::Nelecs_iso_

private

Definition at line 100 of file TopDiLeptonDQM.h.

Nmuons_

MonitorElement* TopDiLeptonDQM::Nmuons_

private

Definition at line 90 of file TopDiLeptonDQM.h.

Nmuons_charge_

MonitorElement* TopDiLeptonDQM::Nmuons_charge_

private

Definition at line 92 of file TopDiLeptonDQM.h.

Nmuons_iso_

MonitorElement* TopDiLeptonDQM::Nmuons_iso_

private

Definition at line 91 of file TopDiLeptonDQM.h.

Ntracks_

MonitorElement* TopDiLeptonDQM::Ntracks_

private

Definition at line 88 of file TopDiLeptonDQM.h.

outfile

std::ofstream TopDiLeptonDQM::outfile

private

Definition at line 56 of file TopDiLeptonDQM.h.

outputFile_

std::string TopDiLeptonDQM::outputFile_

private

Definition at line 50 of file TopDiLeptonDQM.h.

phi_elecs_

MonitorElement* TopDiLeptonDQM::phi_elecs_

private

Definition at line 105 of file TopDiLeptonDQM.h.

phi_muons_

MonitorElement* TopDiLeptonDQM::phi_muons_

private

Definition at line 97 of file TopDiLeptonDQM.h.

pT_elecs_

MonitorElement* TopDiLeptonDQM::pT_elecs_

private

Definition at line 103 of file TopDiLeptonDQM.h.

pT_muons_

MonitorElement* TopDiLeptonDQM::pT_muons_

private

Definition at line 95 of file TopDiLeptonDQM.h.

TriggerEff_

MonitorElement* TopDiLeptonDQM::TriggerEff_

private

Definition at line 87 of file TopDiLeptonDQM.h.

triggerResults_

edm::EDGetTokenT<edm::TriggerResults> TopDiLeptonDQM::triggerResults_

private

Definition at line 51 of file TopDiLeptonDQM.h.

Trigs_

MonitorElement* TopDiLeptonDQM::Trigs_

private

Definition at line 86 of file TopDiLeptonDQM.h.

vertex_

edm::EDGetTokenT<reco::VertexCollection> TopDiLeptonDQM::vertex_

private

Definition at line 66 of file TopDiLeptonDQM.h.

vertex_X_cut_

double TopDiLeptonDQM::vertex_X_cut_

private

Definition at line 67 of file TopDiLeptonDQM.h.

vertex_Y_cut_

double TopDiLeptonDQM::vertex_Y_cut_

private

Definition at line 68 of file TopDiLeptonDQM.h.

vertex_Z_cut_

double TopDiLeptonDQM::vertex_Z_cut_

private

Definition at line 69 of file TopDiLeptonDQM.h.

VxVy_muons_

MonitorElement* TopDiLeptonDQM::VxVy_muons_

private

Definition at line 93 of file TopDiLeptonDQM.h.

Vz_muons_

MonitorElement* TopDiLeptonDQM::Vz_muons_

private

Definition at line 94 of file TopDiLeptonDQM.h.