ZCounting Class Reference

#include <ZCounting.h>

Inheritance diagram for ZCounting:

Public Types
enum	MuonIDTypes {   NoneID, LooseID, MediumID, TightID,   CustomTightID }
enum	MuonIsoTypes { NoneIso, TrackerIso, PFIso }
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 >
using	CacheTypes = CacheContexts< T... >
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T... >
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Public Member Functions
	ZCounting (const edm::ParameterSet &ps)
	~ZCounting () 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
	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
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

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

Private Member Functions
bool	isCustomTightMuon (const reco::Muon &muon)
bool	passMuonID (const reco::Muon &muon, const reco::Vertex *vtx)
bool	passMuonIso (const reco::Muon &muon)

Private Attributes
const double	DRMAX = 0.1
const double	EtaCutL1_
const double	EtaCutL2_
std::vector< std::string >	fMuonHLTNames
std::vector< std::string >	fMuonHLTObjectNames
edm::EDGetTokenT< reco::MuonCollection >	fMuonName_token
edm::EDGetTokenT< reco::VertexCollection >	fPVName_token
edm::EDGetTokenT< reco::TrackCollection >	fTrackName_token
MonitorElement *	h_mass_1HLT_BB
MonitorElement *	h_mass_1HLT_BE
MonitorElement *	h_mass_1HLT_EE
MonitorElement *	h_mass_2HLT_BB
MonitorElement *	h_mass_2HLT_BE
MonitorElement *	h_mass_2HLT_EE
MonitorElement *	h_mass_Glo_fail_BB
MonitorElement *	h_mass_Glo_fail_BE
MonitorElement *	h_mass_Glo_fail_EE
MonitorElement *	h_mass_SIT_fail_BB
MonitorElement *	h_mass_SIT_fail_BE
MonitorElement *	h_mass_SIT_fail_EE
MonitorElement *	h_npv
HLTConfigProvider	hltConfigProvider_
MuonIDTypes	IDType_ {NoneID}
const std::string	IDTypestr_
const double	IsoCut_
MuonIsoTypes	IsoType_ {NoneIso}
const std::string	IsoTypestr_
const int	LumiBin_
const double	LumiMax_
const double	LumiMin_
const int	MassBin_
const double	MassMax_
const double	MassMin_
const double	MUON_BOUND = 0.9
const double	MUON_MASS = 0.105658369
const double	PtCutL1_
const double	PtCutL2_
const int	PVBin_
const double	PVMax_
const double	PVMin_
const edm::InputTag	triggerResultsInputTag_
TriggerTools *	triggers
const double	VtxAbsZCut_
const double	VtxNdofCut_
const double	VtxNTracksFitCut_
const double	VtxRhoCut_

Additional Inherited Members
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 26 of file ZCounting.h.

Member Enumeration Documentation

MuonIDTypes

enum ZCounting::MuonIDTypes

Enumerator
NoneID
LooseID
MediumID
TightID
CustomTightID

Definition at line 31 of file ZCounting.h.

{ NoneID, LooseID, MediumID, TightID, CustomTightID };

MuonIsoTypes

enum ZCounting::MuonIsoTypes

Enumerator
NoneIso
TrackerIso
PFIso

Definition at line 32 of file ZCounting.h.

{ NoneIso, TrackerIso, PFIso };

Constructor & Destructor Documentation

ZCounting()

ZCounting::ZCounting

(

const edm::ParameterSet &

ps

)

Definition at line 20 of file ZCounting.cc.

References TriggerTools::addTriggerRecord(), CustomTightID, DRMAX, edm::ParameterSet::getParameter(), IDType_, IDTypestr_, IsoType_, IsoTypestr_, LooseID, MediumID, NoneID, NoneIso, PFIso, TriggerTools::setDRMAX(), TriggerTools::setTriggerEventToken(), TriggerTools::setTriggerResultsToken(), AlCaHLTBitMon_QueryRunRegistry::string, TightID, TrackerIso, triggerResultsInputTag_, and triggers.

    : triggerResultsInputTag_(iConfig.getParameter<edm::InputTag>("TriggerResults")),
      fPVName_token(consumes<reco::VertexCollection>(
          iConfig.getUntrackedParameter<std::string>("edmPVName", "offlinePrimaryVertices"))),
      fMuonName_token(consumes<reco::MuonCollection>(iConfig.getUntrackedParameter<std::string>("edmName", "muons"))),
      fTrackName_token(
          consumes<reco::TrackCollection>(iConfig.getUntrackedParameter<std::string>("edmTrackName", "generalTracks"))),

      PtCutL1_(iConfig.getUntrackedParameter<double>("PtCutL1")),
      PtCutL2_(iConfig.getUntrackedParameter<double>("PtCutL2")),
      EtaCutL1_(iConfig.getUntrackedParameter<double>("EtaCutL1")),
      EtaCutL2_(iConfig.getUntrackedParameter<double>("EtaCutL2")),

      MassBin_(iConfig.getUntrackedParameter<int>("MassBin")),
      MassMin_(iConfig.getUntrackedParameter<double>("MassMin")),
      MassMax_(iConfig.getUntrackedParameter<double>("MassMax")),

      LumiBin_(iConfig.getUntrackedParameter<int>("LumiBin")),
      LumiMin_(iConfig.getUntrackedParameter<double>("LumiMin")),
      LumiMax_(iConfig.getUntrackedParameter<double>("LumiMax")),

      PVBin_(iConfig.getUntrackedParameter<int>("PVBin")),
      PVMin_(iConfig.getUntrackedParameter<double>("PVMin")),
      PVMax_(iConfig.getUntrackedParameter<double>("PVMax")),

      VtxNTracksFitCut_(iConfig.getUntrackedParameter<double>("VtxNTracksFitMin")),
      VtxNdofCut_(iConfig.getUntrackedParameter<double>("VtxNdofMin")),
      VtxAbsZCut_(iConfig.getUntrackedParameter<double>("VtxAbsZMax")),
      VtxRhoCut_(iConfig.getUntrackedParameter<double>("VtxRhoMax")),

      IDTypestr_(iConfig.getUntrackedParameter<std::string>("IDType")),
      IsoTypestr_(iConfig.getUntrackedParameter<std::string>("IsoType")),
      IsoCut_(iConfig.getUntrackedParameter<double>("IsoCut")) {
  edm::LogInfo("ZCounting") << "Constructor  ZCounting::ZCounting " << std::endl;

  // Trigger settings
  triggers = new TriggerTools();
  triggers->setTriggerResultsToken(consumes<edm::TriggerResults>(triggerResultsInputTag_));
  triggers->setTriggerEventToken(consumes<trigger::TriggerEvent>(iConfig.getParameter<edm::InputTag>("TriggerEvent")));
  triggers->setDRMAX(DRMAX);

  edm::LogVerbatim("ZCounting") << "ZCounting::ZCounting set trigger names";
  const std::vector<std::string> patterns_ = iConfig.getParameter<std::vector<std::string>>("MuonTriggerNames");
  for (const std::string& pattern_ : patterns_) {
    triggers->addTriggerRecord(pattern_);
  }

  if (IDTypestr_ == "Loose")
    IDType_ = LooseID;
  else if (IDTypestr_ == "Medium")
    IDType_ = MediumID;
  else if (IDTypestr_ == "Tight")
    IDType_ = TightID;
  else if (IDTypestr_ == "CustomTight")
    IDType_ = CustomTightID;
  else
    IDType_ = NoneID;

  if (IsoTypestr_ == "Tracker-based")
    IsoType_ = TrackerIso;
  else if (IsoTypestr_ == "PF-based")
    IsoType_ = PFIso;
  else
    IsoType_ = NoneIso;
}

~ZCounting()

ZCounting::~ZCounting ( )

override

Definition at line 89 of file ZCounting.cc.

{ edm::LogInfo("ZCounting") << "Destructor ZCounting::~ZCounting " << std::endl; }

Member Function Documentation

analyze()

void ZCounting::analyze ( edm::Event const &  e, edm::EventSetup const &  eSetup  )

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 263 of file ZCounting.cc.

References HLT_2022v12_cff::eta1, HLT_2022v12_cff::eta2, EtaCutL1_, EtaCutL2_, dqm::impl::MonitorElement::Fill(), fMuonName_token, fPVName_token, fTrackName_token, h_mass_1HLT_BB, h_mass_1HLT_BE, h_mass_1HLT_EE, h_mass_2HLT_BB, h_mass_2HLT_BE, h_mass_2HLT_EE, h_mass_Glo_fail_BB, h_mass_Glo_fail_BE, h_mass_Glo_fail_EE, h_mass_SIT_fail_BB, h_mass_SIT_fail_BE, h_mass_SIT_fail_EE, h_npv, iEvent, reco::isMuon(), edm::HandleBase::isValid(), MassMax_, MassMin_, MUON_BOUND, MUON_MASS, TriggerTools::pass(), passMuonID(), passMuonIso(), TriggerTools::passObj(), HLT_2022v12_cff::pt1, HLT_2022v12_cff::pt2, PtCutL1_, PtCutL2_, MetAnalyzer::pv(), TriggerTools::readEvent(), triggers, VtxAbsZCut_, VtxNdofCut_, VtxNTracksFitCut_, and VtxRhoCut_.

                                                                           {  // Fill event tree on the fly
  edm::LogInfo("ZCounting") << "ZCounting::analyze" << std::endl;

  //-------------------------------
  //--- Vertex
  //-------------------------------
  edm::Handle<reco::VertexCollection> hVertexProduct;
  iEvent.getByToken(fPVName_token, hVertexProduct);
  if (!hVertexProduct.isValid()) {
    edm::LogWarning("ZCounting") << "ZCounting::analyze - no valid primary vertex product found" << std::endl;
    return;
  }

  const reco::Vertex* pv = nullptr;
  int nvtx = 0;

  for (auto const& itVtx : *hVertexProduct) {
    if (itVtx.isFake())
      continue;
    if (itVtx.tracksSize() < VtxNTracksFitCut_)
      continue;
    if (itVtx.ndof() < VtxNdofCut_)
      continue;
    if (fabs(itVtx.z()) > VtxAbsZCut_)
      continue;
    if (itVtx.position().Rho() > VtxRhoCut_)
      continue;

    if (nvtx == 0) {
      pv = &itVtx;
    }
    nvtx++;
  }

  h_npv->Fill(iEvent.luminosityBlock(), nvtx);

  //-------------------------------
  //--- Trigger
  //-------------------------------
  triggers->readEvent(iEvent);

  // Trigger requirement
  if (!triggers->pass())
    return;

  //-------------------------------
  //--- Muons and Tracks
  //-------------------------------
  edm::Handle<reco::MuonCollection> hMuonProduct;
  iEvent.getByToken(fMuonName_token, hMuonProduct);
  if (!hMuonProduct.isValid())
    return;

  edm::Handle<reco::TrackCollection> hTrackProduct;
  iEvent.getByToken(fTrackName_token, hTrackProduct);
  if (!hTrackProduct.isValid())
    return;

  TLorentzVector vTag(0., 0., 0., 0.);
  TLorentzVector vProbe(0., 0., 0., 0.);
  TLorentzVector vTrack(0., 0., 0., 0.);

  // Tag loop
  for (auto const& itMu1 : *hMuonProduct) {
    const float pt1 = itMu1.muonBestTrack()->pt();
    const float eta1 = itMu1.muonBestTrack()->eta();
    const float phi1 = itMu1.muonBestTrack()->phi();
    const float q1 = itMu1.muonBestTrack()->charge();

    // Tag selection: kinematic cuts, lepton selection and trigger matching
    if (pt1 < PtCutL1_)
      continue;
    if (fabs(eta1) > EtaCutL1_)
      continue;
    if (!(passMuonID(itMu1, pv) && passMuonIso(itMu1)))
      continue;
    if (!triggers->passObj(eta1, phi1))
      continue;

    vTag.SetPtEtaPhiM(pt1, eta1, phi1, MUON_MASS);

    bool isTagCentral = false;
    if (fabs(eta1) < MUON_BOUND)
      isTagCentral = true;

    // Probe loop over muons
    for (auto const& itMu2 : *hMuonProduct) {
      if (&itMu2 == &itMu1)
        continue;

      const float pt2 = itMu2.muonBestTrack()->pt();
      const float eta2 = itMu2.muonBestTrack()->eta();
      const float phi2 = itMu2.muonBestTrack()->phi();
      const float q2 = itMu2.muonBestTrack()->charge();

      // Probe selection: kinematic cuts and opposite charge requirement
      if (pt2 < PtCutL2_)
        continue;
      if (fabs(eta2) > EtaCutL2_)
        continue;
      if (q1 == q2)
        continue;

      vProbe.SetPtEtaPhiM(pt2, eta2, phi2, MUON_MASS);

      // Mass window
      TLorentzVector vDilep = vTag + vProbe;
      float dilepMass = vDilep.M();
      if ((dilepMass < MassMin_) || (dilepMass > MassMax_))
        continue;

      bool isProbeCentral = false;
      if (fabs(eta2) < MUON_BOUND)
        isProbeCentral = true;

      // Determine event category for efficiency calculation
      if (passMuonID(itMu2, pv) && passMuonIso(itMu2)) {
        if (triggers->passObj(eta2, phi2)) {
          // category 2HLT: both muons passing trigger requirements
          if (&itMu1 > &itMu2)
            continue;  // make sure we don't double count MuMu2HLT category

          if (isTagCentral && isProbeCentral) {
            h_mass_2HLT_BB->Fill(iEvent.luminosityBlock(), dilepMass);
          } else if (!isTagCentral && !isProbeCentral) {
            h_mass_2HLT_EE->Fill(iEvent.luminosityBlock(), dilepMass);
          } else {
            h_mass_2HLT_BE->Fill(iEvent.luminosityBlock(), dilepMass);
          }
        } else {
          // category 1HLT: only one muon passes trigger
          if (isTagCentral && isProbeCentral) {
            h_mass_1HLT_BB->Fill(iEvent.luminosityBlock(), dilepMass);
          } else if (!isTagCentral && !isProbeCentral) {
            h_mass_1HLT_EE->Fill(iEvent.luminosityBlock(), dilepMass);
          } else {
            h_mass_1HLT_BE->Fill(iEvent.luminosityBlock(), dilepMass);
          }
        }
      } else if (itMu2.isGlobalMuon()) {
        // category Glo: probe is a Global muon but failing selection
        if (isTagCentral && isProbeCentral) {
          h_mass_SIT_fail_BB->Fill(iEvent.luminosityBlock(), dilepMass);
        } else if (!isTagCentral && !isProbeCentral) {
          h_mass_SIT_fail_EE->Fill(iEvent.luminosityBlock(), dilepMass);
        } else {
          h_mass_SIT_fail_BE->Fill(iEvent.luminosityBlock(), dilepMass);
        }
      } else if (itMu2.isStandAloneMuon()) {
        // category Sta: probe is a Standalone muon
        if (isTagCentral && isProbeCentral) {
          h_mass_Glo_fail_BB->Fill(iEvent.luminosityBlock(), dilepMass);
        } else if (!isTagCentral && !isProbeCentral) {
          h_mass_Glo_fail_EE->Fill(iEvent.luminosityBlock(), dilepMass);
        } else {
          h_mass_Glo_fail_BE->Fill(iEvent.luminosityBlock(), dilepMass);
        }
      } else if (itMu2.innerTrack()->hitPattern().trackerLayersWithMeasurement() >= 6 &&
                 itMu2.innerTrack()->hitPattern().numberOfValidPixelHits() >= 1) {
        // cateogry Trk: probe is a tracker track
        if (isTagCentral && isProbeCentral) {
          h_mass_Glo_fail_BB->Fill(iEvent.luminosityBlock(), dilepMass);
        } else if (!isTagCentral && !isProbeCentral) {
          h_mass_Glo_fail_EE->Fill(iEvent.luminosityBlock(), dilepMass);
        } else {
          h_mass_Glo_fail_BE->Fill(iEvent.luminosityBlock(), dilepMass);
        }
      }
    }  // End of probe loop over muons

    // Probe loop over tracks, only for standalone efficiency calculation
    for (auto const& itTrk : *hTrackProduct) {
      // Check track is not a muon
      bool isMuon = false;
      for (auto const& itMu : *hMuonProduct) {
        if (itMu.innerTrack().isNonnull() && itMu.innerTrack().get() == &itTrk) {
          isMuon = true;
          break;
        }
      }
      if (isMuon)
        continue;

      const float pt2 = itTrk.pt();
      const float eta2 = itTrk.eta();
      const float phi2 = itTrk.phi();
      const float q2 = itTrk.charge();

      // Probe selection:  kinematic cuts and opposite charge requirement
      if (pt2 < PtCutL2_)
        continue;
      if (fabs(eta2) > EtaCutL2_)
        continue;
      if (q1 == q2)
        continue;

      vTrack.SetPtEtaPhiM(pt2, eta2, phi2, MUON_MASS);

      TLorentzVector vDilep = vTag + vTrack;
      float dilepMass = vDilep.M();
      if ((dilepMass < MassMin_) || (dilepMass > MassMax_))
        continue;

      bool isTrackCentral = false;
      if (fabs(eta2) < MUON_BOUND)
        isTrackCentral = true;

      if (itTrk.hitPattern().trackerLayersWithMeasurement() >= 6 && itTrk.hitPattern().numberOfValidPixelHits() >= 1) {
        if (isTagCentral && isTrackCentral) {
          h_mass_Glo_fail_BB->Fill(iEvent.luminosityBlock(), dilepMass);
        } else if (!isTagCentral && !isTrackCentral) {
          h_mass_Glo_fail_EE->Fill(iEvent.luminosityBlock(), dilepMass);
        } else {
          h_mass_Glo_fail_BE->Fill(iEvent.luminosityBlock(), dilepMass);
        }
      }
    }  //End of probe loop over tracks
  }    //End of tag loop
}

bookHistograms()

void ZCounting::bookHistograms ( DQMStore::IBooker &  ibooker_, edm::Run const &  , edm::EventSetup const &    )

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 120 of file ZCounting.cc.

References dqm::implementation::IBooker::book2D(), dqm::implementation::NavigatorBase::cd(), h_mass_1HLT_BB, h_mass_1HLT_BE, h_mass_1HLT_EE, h_mass_2HLT_BB, h_mass_2HLT_BE, h_mass_2HLT_EE, h_mass_Glo_fail_BB, h_mass_Glo_fail_BE, h_mass_Glo_fail_EE, h_mass_SIT_fail_BB, h_mass_SIT_fail_BE, h_mass_SIT_fail_EE, h_npv, LumiBin_, LumiMax_, LumiMin_, MassBin_, MassMax_, MassMin_, PVBin_, PVMax_, PVMin_, dqm::impl::MonitorElement::setAxisTitle(), and dqm::implementation::NavigatorBase::setCurrentFolder().

                                                                                             {
  edm::LogInfo("ZCounting") << "ZCounting::bookHistograms" << std::endl;
  ibooker_.cd();
  ibooker_.setCurrentFolder("ZCounting/Histograms");

  // Muon histograms
  h_mass_2HLT_BB = ibooker_.book2D("h_mass_2HLT_BB",
                                   "Both muon pass HLT in barrel-barrel",
                                   LumiBin_,
                                   LumiMin_,
                                   LumiMax_,
                                   MassBin_,
                                   MassMin_,
                                   MassMax_);
  h_mass_2HLT_BE = ibooker_.book2D("h_mass_2HLT_BE",
                                   "Both muon pass HLT passing in barrel-endcap",
                                   LumiBin_,
                                   LumiMin_,
                                   LumiMax_,
                                   MassBin_,
                                   MassMin_,
                                   MassMax_);
  h_mass_2HLT_EE = ibooker_.book2D("h_mass_2HLT_EE",
                                   "Both muon pass HLT passing in endcap-endcap",
                                   LumiBin_,
                                   LumiMin_,
                                   LumiMax_,
                                   MassBin_,
                                   MassMin_,
                                   MassMax_);
  h_mass_1HLT_BB = ibooker_.book2D("h_mass_1HLT_BB",
                                   "One muon pass HLT in barrel-barrel",
                                   LumiBin_,
                                   LumiMin_,
                                   LumiMax_,
                                   MassBin_,
                                   MassMin_,
                                   MassMax_);
  h_mass_1HLT_BE = ibooker_.book2D("h_mass_1HLT_BE",
                                   "One muon pass HLT passing in barrel-endcap",
                                   LumiBin_,
                                   LumiMin_,
                                   LumiMax_,
                                   MassBin_,
                                   MassMin_,
                                   MassMax_);
  h_mass_1HLT_EE = ibooker_.book2D("h_mass_1HLT_EE",
                                   "One muon pass HLT passing in endcap-endcap",
                                   LumiBin_,
                                   LumiMin_,
                                   LumiMax_,
                                   MassBin_,
                                   MassMin_,
                                   MassMax_);

  h_mass_SIT_fail_BB = ibooker_.book2D("h_mass_SIT_fail_BB",
                                       "Muon SIT failing barrel-barrel",
                                       LumiBin_,
                                       LumiMin_,
                                       LumiMax_,
                                       MassBin_,
                                       MassMin_,
                                       MassMax_);
  h_mass_SIT_fail_BE = ibooker_.book2D("h_mass_SIT_fail_BE",
                                       "Muon SIT failing barrel-endcap",
                                       LumiBin_,
                                       LumiMin_,
                                       LumiMax_,
                                       MassBin_,
                                       MassMin_,
                                       MassMax_);

  h_mass_SIT_fail_EE = ibooker_.book2D("h_mass_SIT_fail_EE",
                                       "Muon SIT failing endcap-endcap",
                                       LumiBin_,
                                       LumiMin_,
                                       LumiMax_,
                                       MassBin_,
                                       MassMin_,
                                       MassMax_);

  h_mass_Glo_fail_BB = ibooker_.book2D("h_mass_Glo_fail_BB",
                                       "Muon Glo failing barrel-barrel",
                                       LumiBin_,
                                       LumiMin_,
                                       LumiMax_,
                                       MassBin_,
                                       MassMin_,
                                       MassMax_);
  h_mass_Glo_fail_BE = ibooker_.book2D("h_mass_Glo_fail_BE",
                                       "Muon Glo failing barrel-endcap",
                                       LumiBin_,
                                       LumiMin_,
                                       LumiMax_,
                                       MassBin_,
                                       MassMin_,
                                       MassMax_);

  h_mass_Glo_fail_EE = ibooker_.book2D("h_mass_Glo_fail_EE",
                                       "Muon Glo failing endcap-endcap",
                                       LumiBin_,
                                       LumiMin_,
                                       LumiMax_,
                                       MassBin_,
                                       MassMin_,
                                       MassMax_);

  h_npv = ibooker_.book2D(
      "h_npv", "Events with valid primary vertex", LumiBin_, LumiMin_, LumiMax_, PVBin_, PVMin_, PVMax_);

  // Axis titles
  h_mass_2HLT_BB->setAxisTitle("luminosity section", 1);
  h_mass_2HLT_BE->setAxisTitle("luminosity section", 1);
  h_mass_2HLT_EE->setAxisTitle("luminosity section", 1);
  h_mass_1HLT_BB->setAxisTitle("luminosity section", 1);
  h_mass_1HLT_BE->setAxisTitle("luminosity section", 1);
  h_mass_1HLT_EE->setAxisTitle("luminosity section", 1);
  h_mass_SIT_fail_BB->setAxisTitle("luminosity section", 1);
  h_mass_SIT_fail_BE->setAxisTitle("luminosity section", 1);
  h_mass_SIT_fail_EE->setAxisTitle("luminosity section", 1);
  h_mass_Glo_fail_BB->setAxisTitle("luminosity section", 1);
  h_mass_Glo_fail_BE->setAxisTitle("luminosity section", 1);
  h_mass_Glo_fail_EE->setAxisTitle("luminosity section", 1);
  h_mass_2HLT_BB->setAxisTitle("tag and probe mass", 2);
  h_mass_2HLT_BE->setAxisTitle("tag and probe mass", 2);
  h_mass_2HLT_EE->setAxisTitle("tag and probe mass", 2);
  h_mass_1HLT_BB->setAxisTitle("tag and probe mass", 2);
  h_mass_1HLT_BE->setAxisTitle("tag and probe mass", 2);
  h_mass_1HLT_EE->setAxisTitle("tag and probe mass", 2);
  h_mass_SIT_fail_BB->setAxisTitle("tag and probe mass", 2);
  h_mass_SIT_fail_BE->setAxisTitle("tag and probe mass", 2);
  h_mass_SIT_fail_EE->setAxisTitle("tag and probe mass", 2);
  h_mass_Glo_fail_BB->setAxisTitle("tag and probe mass", 2);
  h_mass_Glo_fail_BE->setAxisTitle("tag and probe mass", 2);
  h_mass_Glo_fail_EE->setAxisTitle("tag and probe mass", 2);
  h_npv->setAxisTitle("luminosity section", 1);
  h_npv->setAxisTitle("number of primary vertices", 2);
}

dqmBeginRun()

void ZCounting::dqmBeginRun ( edm::Run const &  iRun, edm::EventSetup const &  iSetup  )

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 94 of file ZCounting.cc.

References hltConfigProvider_, edm::RunBase::id(), HLTConfigProvider::init(), TriggerTools::initHLTObjects(), edm::InputTag::process(), HLTConfigProvider::processName(), HLTConfigProvider::size(), HLTConfigProvider::tableName(), triggerResultsInputTag_, and triggers.

                                                                           {
  edm::LogInfo("ZCounting") << "ZCounting::beginRun" << std::endl;

  // initialize triggers

  edm::LogVerbatim("ZCounting") << "ZCounting::dqmBeginRun now at " << iRun.id();
  bool hltChanged_ = true;
  if (hltConfigProvider_.init(iRun, iSetup, triggerResultsInputTag_.process(), hltChanged_)) {
    edm::LogVerbatim("ZCounting") << "ZCounting::dqmBeginRun [TriggerObjMatchValueMapsProducer::beginRun] "
                                     "HLTConfigProvider initialized [processName() = \""
                                  << hltConfigProvider_.processName() << "\", tableName() = \""
                                  << hltConfigProvider_.tableName() << "\", size() = " << hltConfigProvider_.size()
                                  << "]";
  } else {
    edm::LogError("ZCounting") << "ZCounting::dqmBeginRun Initialization of HLTConfigProvider failed for Run="
                               << iRun.id() << " (process=\"" << triggerResultsInputTag_.process()
                               << "\") -> plugin will not produce outputs for this Run";
    return;
  }

  triggers->initHLTObjects(hltConfigProvider_);
}

isCustomTightMuon()

bool ZCounting::isCustomTightMuon ( const reco::Muon &  muon )

private

Definition at line 489 of file ZCounting.cc.

References muon::GlobalMuonPromptTight, and muon::isGoodMuon().

Referenced by passMuonID().

                                                      {
  if (!muon.isPFMuon() || !muon.isGlobalMuon())
    return false;

  bool muID = isGoodMuon(muon, muon::GlobalMuonPromptTight) && (muon.numberOfMatchedStations() > 1);

  bool muIdAndHits = muID && muon.innerTrack()->hitPattern().trackerLayersWithMeasurement() > 5 &&
                     muon.innerTrack()->hitPattern().numberOfValidPixelHits() > 0;

  return muIdAndHits;
}

passMuonID()

bool ZCounting::passMuonID ( const reco::Muon &  muon, const reco::Vertex *  vtx  )

private

Definition at line 502 of file ZCounting.cc.

References CustomTightID, IDType_, isCustomTightMuon(), muon::isLooseMuon(), muon::isMediumMuon(), muon::isTightMuon(), LooseID, MediumID, NoneID, TightID, and extraflags_cff::vtx.

Referenced by analyze().

                                                                      {
  // Muon ID selection, using internal function "DataFormats/MuonReco/src/MuonSelectors.cc
  switch (IDType_) {
    case LooseID:
      return muon::isLooseMuon(muon);
    case MediumID:
      return muon::isMediumMuon(muon);
    case CustomTightID:
      return isCustomTightMuon(muon);
    case TightID:
      return vtx != nullptr && muon::isTightMuon(muon, *vtx);
    case NoneID:
      return true;
  }
  return false;
}

passMuonIso()

bool ZCounting::passMuonIso ( const reco::Muon &  muon )

private

Definition at line 520 of file ZCounting.cc.

References IsoCut_, IsoType_, SiStripPI::max, NoneIso, PFIso, and TrackerIso.

Referenced by analyze().

                                                {
  //Muon isolation selection, up-to-date with MUO POG recommendation
  switch (IsoType_) {
    case TrackerIso:
      return muon.isolationR03().sumPt < IsoCut_;
    case PFIso:
      return muon.pfIsolationR04().sumChargedHadronPt +
                 std::max(0.,
                          muon.pfIsolationR04().sumNeutralHadronEt + muon.pfIsolationR04().sumPhotonEt -
                              0.5 * muon.pfIsolationR04().sumPUPt) <
             IsoCut_;
    case NoneIso:
      return true;
  }

  return false;
}

Member Data Documentation

DRMAX

const double ZCounting::DRMAX = 0.1

private

Definition at line 93 of file ZCounting.h.

Referenced by ZCounting().

EtaCutL1_

const double ZCounting::EtaCutL1_

private

Definition at line 60 of file ZCounting.h.

Referenced by analyze().

EtaCutL2_

const double ZCounting::EtaCutL2_

private

Definition at line 61 of file ZCounting.h.

Referenced by analyze().

fMuonHLTNames

std::vector<std::string> ZCounting::fMuonHLTNames

private

Definition at line 51 of file ZCounting.h.

fMuonHLTObjectNames

std::vector<std::string> ZCounting::fMuonHLTObjectNames

private

Definition at line 52 of file ZCounting.h.

fMuonName_token

edm::EDGetTokenT<reco::MuonCollection> ZCounting::fMuonName_token

private

Definition at line 50 of file ZCounting.h.

Referenced by analyze().

fPVName_token

edm::EDGetTokenT<reco::VertexCollection> ZCounting::fPVName_token

private

Definition at line 47 of file ZCounting.h.

Referenced by analyze().

fTrackName_token

edm::EDGetTokenT<reco::TrackCollection> ZCounting::fTrackName_token

private

Definition at line 55 of file ZCounting.h.

Referenced by analyze().

h_mass_1HLT_BB

MonitorElement* ZCounting::h_mass_1HLT_BB

private

Definition at line 106 of file ZCounting.h.

Referenced by analyze(), and bookHistograms().

h_mass_1HLT_BE

MonitorElement* ZCounting::h_mass_1HLT_BE

private

Definition at line 107 of file ZCounting.h.

Referenced by analyze(), and bookHistograms().

h_mass_1HLT_EE

MonitorElement* ZCounting::h_mass_1HLT_EE

private

Definition at line 108 of file ZCounting.h.

Referenced by analyze(), and bookHistograms().

h_mass_2HLT_BB

MonitorElement* ZCounting::h_mass_2HLT_BB

private

Definition at line 102 of file ZCounting.h.

Referenced by analyze(), and bookHistograms().

h_mass_2HLT_BE

MonitorElement* ZCounting::h_mass_2HLT_BE

private

Definition at line 103 of file ZCounting.h.

Referenced by analyze(), and bookHistograms().

h_mass_2HLT_EE

MonitorElement* ZCounting::h_mass_2HLT_EE

private

Definition at line 104 of file ZCounting.h.

Referenced by analyze(), and bookHistograms().

h_mass_Glo_fail_BB

MonitorElement* ZCounting::h_mass_Glo_fail_BB

private

Definition at line 114 of file ZCounting.h.

Referenced by analyze(), and bookHistograms().

h_mass_Glo_fail_BE

MonitorElement* ZCounting::h_mass_Glo_fail_BE

private

Definition at line 115 of file ZCounting.h.

Referenced by analyze(), and bookHistograms().

h_mass_Glo_fail_EE

MonitorElement* ZCounting::h_mass_Glo_fail_EE

private

Definition at line 116 of file ZCounting.h.

Referenced by analyze(), and bookHistograms().

h_mass_SIT_fail_BB

MonitorElement* ZCounting::h_mass_SIT_fail_BB

private

Definition at line 110 of file ZCounting.h.

Referenced by analyze(), and bookHistograms().

h_mass_SIT_fail_BE

MonitorElement* ZCounting::h_mass_SIT_fail_BE

private

Definition at line 111 of file ZCounting.h.

Referenced by analyze(), and bookHistograms().

h_mass_SIT_fail_EE

MonitorElement* ZCounting::h_mass_SIT_fail_EE

private

Definition at line 112 of file ZCounting.h.

Referenced by analyze(), and bookHistograms().

h_npv

MonitorElement* ZCounting::h_npv

private

Definition at line 99 of file ZCounting.h.

Referenced by analyze(), and bookHistograms().

hltConfigProvider_

HLTConfigProvider ZCounting::hltConfigProvider_

private

Definition at line 89 of file ZCounting.h.

Referenced by dqmBeginRun().

IDType_

MuonIDTypes ZCounting::IDType_ {NoneID}

private

Definition at line 85 of file ZCounting.h.

Referenced by passMuonID(), and ZCounting().

IDTypestr_

const std::string ZCounting::IDTypestr_

private

Definition at line 80 of file ZCounting.h.

Referenced by ZCounting().

IsoCut_

const double ZCounting::IsoCut_

private

Definition at line 82 of file ZCounting.h.

Referenced by passMuonIso().

IsoType_

MuonIsoTypes ZCounting::IsoType_ {NoneIso}

private

Definition at line 86 of file ZCounting.h.

Referenced by passMuonIso(), and ZCounting().

IsoTypestr_

const std::string ZCounting::IsoTypestr_

private

Definition at line 81 of file ZCounting.h.

Referenced by ZCounting().

LumiBin_

const int ZCounting::LumiBin_

private

Definition at line 67 of file ZCounting.h.

Referenced by bookHistograms().

LumiMax_

const double ZCounting::LumiMax_

private

Definition at line 69 of file ZCounting.h.

Referenced by bookHistograms().

LumiMin_

const double ZCounting::LumiMin_

private

Definition at line 68 of file ZCounting.h.

Referenced by bookHistograms().

MassBin_

const int ZCounting::MassBin_

private

Definition at line 63 of file ZCounting.h.

Referenced by bookHistograms().

MassMax_

const double ZCounting::MassMax_

private

Definition at line 65 of file ZCounting.h.

Referenced by analyze(), and bookHistograms().

MassMin_

const double ZCounting::MassMin_

private

Definition at line 64 of file ZCounting.h.

Referenced by analyze(), and bookHistograms().

MUON_BOUND

const double ZCounting::MUON_BOUND = 0.9

private

Definition at line 96 of file ZCounting.h.

Referenced by analyze().

MUON_MASS

const double ZCounting::MUON_MASS = 0.105658369

private

Definition at line 95 of file ZCounting.h.

Referenced by analyze().

PtCutL1_

const double ZCounting::PtCutL1_

private

Definition at line 58 of file ZCounting.h.

Referenced by analyze().

PtCutL2_

const double ZCounting::PtCutL2_

private

Definition at line 59 of file ZCounting.h.

Referenced by analyze().

PVBin_

const int ZCounting::PVBin_

private

Definition at line 71 of file ZCounting.h.

Referenced by bookHistograms().

PVMax_

const double ZCounting::PVMax_

private

Definition at line 73 of file ZCounting.h.

Referenced by bookHistograms().

PVMin_

const double ZCounting::PVMin_

private

Definition at line 72 of file ZCounting.h.

Referenced by bookHistograms().

triggerResultsInputTag_

const edm::InputTag ZCounting::triggerResultsInputTag_

private

Definition at line 46 of file ZCounting.h.

Referenced by dqmBeginRun(), and ZCounting().

triggers

TriggerTools* ZCounting::triggers

private

Definition at line 90 of file ZCounting.h.

Referenced by analyze(), dqmBeginRun(), and ZCounting().

VtxAbsZCut_

const double ZCounting::VtxAbsZCut_

private

Definition at line 77 of file ZCounting.h.

Referenced by analyze().

VtxNdofCut_

const double ZCounting::VtxNdofCut_

private

Definition at line 76 of file ZCounting.h.

Referenced by analyze().

VtxNTracksFitCut_

const double ZCounting::VtxNTracksFitCut_

private

Definition at line 75 of file ZCounting.h.

Referenced by analyze().

VtxRhoCut_

const double ZCounting::VtxRhoCut_

private

Definition at line 78 of file ZCounting.h.

Referenced by analyze().