EwkMuDQM Class Reference

#include <EwkMuDQM.h>

Inheritance diagram for EwkMuDQM:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
	EwkMuDQM (const edm::ParameterSet &)
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	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
void	dqmBeginRun (const edm::Run &, const edm::EventSetup &) override
void	init_histograms ()
Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Private Attributes
MonitorElement *	acop_after_
MonitorElement *	acop_before_
double	acopCut_
edm::EDGetTokenT< reco::BeamSpot >	beamSpotTag_
MonitorElement *	dimuonmass_afterZ_
MonitorElement *	dimuonmass_before_
double	dimuonMassMax_
double	dimuonMassMin_
MonitorElement *	dxy1_afterZ_
MonitorElement *	dxy2_afterZ_
MonitorElement *	dxy_after_
MonitorElement *	dxy_before_
double	dxyCut_
double	eJetMin_
MonitorElement *	eta1_afterZ_
MonitorElement *	eta2_afterZ_
MonitorElement *	eta_after_
MonitorElement *	eta_before_
double	etaCut_
MonitorElement *	goodewkmuon1_afterZ_
MonitorElement *	goodewkmuon2_afterZ_
MonitorElement *	goodewkmuon_after_
MonitorElement *	goodewkmuon_before_
HLTPrescaleProvider	hltPrescaleProvider_
bool	isAlsoTrackerMuon_
bool	isCombinedIso_
MonitorElement *	iso1_afterZ_
MonitorElement *	iso2_afterZ_
MonitorElement *	iso_after_
MonitorElement *	iso_before_
double	isoCut03_
bool	isRelativeIso_
bool	isValidHltConfig_
edm::InputTag	jetTag_
edm::EDGetTokenT< edm::View< reco::Jet > >	jetToken_
MonitorElement *	leadingjet_eta_after_
MonitorElement *	leadingjet_eta_afterZ_
MonitorElement *	leadingjet_eta_before_
MonitorElement *	leadingjet_pt_after_
MonitorElement *	leadingjet_pt_afterZ_
MonitorElement *	leadingjet_pt_before_
MonitorElement *	met_after_
MonitorElement *	met_afterZ_
MonitorElement *	met_before_
double	metMax_
double	metMin_
edm::InputTag	metTag_
edm::EDGetTokenT< edm::View< reco::MET > >	metToken_
MonitorElement *	mt_after_
MonitorElement *	mt_before_
double	mtMax_
double	mtMin_
MonitorElement *	muoncharge_after_
MonitorElement *	muoncharge_afterZ_
MonitorElement *	muoncharge_before_
int	muonHitsCut_
edm::EDGetTokenT< edm::View< reco::Muon > >	muonTag_
MonitorElement *	n_zselPt1thr_
MonitorElement *	n_zselPt2thr_
unsigned int	nall
MonitorElement *	ngoodmuons_
unsigned int	nhlt
unsigned int	niso
int	nJetMax_
MonitorElement *	njets_after_
MonitorElement *	njets_afterZ_
MonitorElement *	njets_before_
int	nMatchesCut_
unsigned int	nmet
MonitorElement *	nmuons_
double	normalizedChi2Cut_
MonitorElement *	npfph_
MonitorElement *	nph_
int	nPhoMax_
MonitorElement *	npvs_after_
MonitorElement *	npvs_afterZ_
MonitorElement *	npvs_before_
unsigned int	nrec
unsigned int	nsel
unsigned int	nz
MonitorElement *	pfphEta_
MonitorElement *	pfphPt_
MonitorElement *	phEta_
edm::EDGetTokenT< edm::View< reco::Photon > >	phoTag_
MonitorElement *	phPt_
int	pixelHitsCut_
MonitorElement *	pt1_afterZ_
MonitorElement *	pt2_afterZ_
MonitorElement *	pt_after_
MonitorElement *	pt_before_
double	ptCut_
MonitorElement *	ptDiffPM_afterZ_
MonitorElement *	ptDiffPM_before_
double	ptThrForPhoton_
double	ptThrForZ1_
double	ptThrForZ2_
int	trackerHitsCut_
MonitorElement *	trig_after_
MonitorElement *	trig_before_
std::vector< std::string >	trigPathNames_
edm::EDGetTokenT< edm::TriggerResults >	trigTag_
edm::EDGetTokenT< edm::View< reco::Vertex > >	vertexTag_
MonitorElement *	ztrig_afterZ_

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 >
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
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

DQM offline for EWKMu

Definition at line 27 of file EwkMuDQM.h.

Constructor & Destructor Documentation

EwkMuDQM()

EwkMuDQM::EwkMuDQM

(

const edm::ParameterSet &

cfg

)

Definition at line 36 of file EwkMuDQM.cc.

References isValidHltConfig_.

    :  // Input collections
      metTag_(cfg.getUntrackedParameter<edm::InputTag>("METTag", edm::InputTag("pfmet"))),
      jetTag_(cfg.getUntrackedParameter<edm::InputTag>("JetTag", edm::InputTag("ak4PFJets"))),
      //trigTag_(consumes<edm::TriggerResults>(
      // cfg.getUntrackedParameter<edm::InputTag>(
      // "TrigTag", edm::InputTag("TriggerResults::HLT")))),
      muonTag_(consumes<edm::View<reco::Muon> >(
          cfg.getUntrackedParameter<edm::InputTag>("MuonTag", edm::InputTag("muons")))),
      metToken_(
          consumes<edm::View<reco::MET> >(cfg.getUntrackedParameter<edm::InputTag>("METTag", edm::InputTag("pfmet")))),
      jetToken_(consumes<edm::View<reco::Jet> >(
          cfg.getUntrackedParameter<edm::InputTag>("JetTag", edm::InputTag("ak4PFJets")))),
      phoTag_(consumes<edm::View<reco::Photon> >(
          cfg.getUntrackedParameter<edm::InputTag>("phoTag", edm::InputTag("photons")))),
      vertexTag_(consumes<edm::View<reco::Vertex> >(
          cfg.getUntrackedParameter<edm::InputTag>("VertexTag", edm::InputTag("offlinePrimaryVertices")))),
      beamSpotTag_(consumes<reco::BeamSpot>(
          cfg.getUntrackedParameter<edm::InputTag>("beamSpotTag", edm::InputTag("offlineBeamSpot")))),
      // trigPathNames_(cfg.getUntrackedParameter<std::vector<std::string> >(
      //   "TrigPathNames")),

      // Muon quality cuts
      isAlsoTrackerMuon_(cfg.getUntrackedParameter<bool>("IsAlsoTrackerMuon", true)),  // Glb muon also tracker muon
      dxyCut_(cfg.getUntrackedParameter<double>("DxyCut", 0.2)),                       // dxy < 0.2 cm
      normalizedChi2Cut_(
          cfg.getUntrackedParameter<double>("NormalizedChi2Cut", 10.)),  // chi2/ndof (of global fit) <10.0
      trackerHitsCut_(cfg.getUntrackedParameter<int>("TrackerHitsCut",
                                                     11)),               // Tracker Hits >10
      pixelHitsCut_(cfg.getUntrackedParameter<int>("PixelHitsCut", 1)),  // Pixel Hits >0
      muonHitsCut_(cfg.getUntrackedParameter<int>("MuonHitsCut",
                                                  1)),                 // Valid Muon Hits >0
      nMatchesCut_(cfg.getUntrackedParameter<int>("NMatchesCut", 2)),  // At least 2 Chambers with matches

      // W-boson cuts
      isRelativeIso_(cfg.getUntrackedParameter<bool>("IsRelativeIso", true)),
      isCombinedIso_(cfg.getUntrackedParameter<bool>("IsCombinedIso", false)),
      isoCut03_(cfg.getUntrackedParameter<double>("IsoCut03", 0.1)),
      acopCut_(cfg.getUntrackedParameter<double>("AcopCut", 999.)),
      metMin_(cfg.getUntrackedParameter<double>("MetMin", -999999.)),
      metMax_(cfg.getUntrackedParameter<double>("MetMax", 999999.)),
      mtMin_(cfg.getUntrackedParameter<double>("MtMin", 50.)),
      mtMax_(cfg.getUntrackedParameter<double>("MtMax", 200.)),
      ptCut_(cfg.getUntrackedParameter<double>("PtCut", 20.)),
      etaCut_(cfg.getUntrackedParameter<double>("EtaCut", 2.4)),

      // Z rejection
      ptThrForZ1_(cfg.getUntrackedParameter<double>("PtThrForZ1", 20.)),
      ptThrForZ2_(cfg.getUntrackedParameter<double>("PtThrForZ2", 10.)),

      // Z selection
      dimuonMassMin_(cfg.getUntrackedParameter<double>("dimuonMassMin", 80.)),
      dimuonMassMax_(cfg.getUntrackedParameter<double>("dimuonMassMax", 120.)),

      // Top rejection
      eJetMin_(cfg.getUntrackedParameter<double>("EJetMin", 999999.)),
      nJetMax_(cfg.getUntrackedParameter<int>("NJetMax", 999999)),

      // Photon cuts
      ptThrForPhoton_(cfg.getUntrackedParameter<double>("ptThrForPhoton", 5.)),
      nPhoMax_(cfg.getUntrackedParameter<int>("nPhoMax", 999999)),
      hltPrescaleProvider_(cfg, consumesCollector(), *this) {
  isValidHltConfig_ = false;
}

Member Function Documentation

analyze()

void EwkMuDQM::analyze ( const edm::Event &  ev, const edm::EventSetup &  iSet  )

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 243 of file EwkMuDQM.cc.

References acop_after_, acop_before_, acopCut_, beamSpotTag_, ALCARECOTkAlJpsiMuMu_cff::charge, reco::LeafCandidate::charge(), EcalCondDBWriter_cfi::cosmic, dimuonmass_afterZ_, dimuonmass_before_, dimuonMassMax_, dimuonMassMin_, HLT_2024v11_cff::distance, PVValHelper::dxy, dxy1_afterZ_, dxy2_afterZ_, dxy_after_, dxy_before_, dxyCut_, eJetMin_, reco::MuonIsolation::emEt, PVValHelper::eta, reco::LeafCandidate::eta(), eta1_afterZ_, HLT_2024v11_cff::eta2, eta2_afterZ_, eta_after_, eta_before_, etaCut_, makeMEIFBenchmarkPlots::ev, dqm::impl::MonitorElement::Fill(), reco::Muon::globalTrack(), goodewkmuon1_afterZ_, goodewkmuon2_afterZ_, goodewkmuon_after_, goodewkmuon_before_, reco::MuonIsolation::hadEt, hltPrescaleProvider_, mps_fire::i, reco::Muon::innerTrack(), isCombinedIso_, reco::Muon::isGlobalMuon(), iso1_afterZ_, iso2_afterZ_, iso_after_, iso_before_, isoCut03_, reco::Muon::isolationR03(), isRelativeIso_, reco::Muon::isTrackerMuon(), dqmiolumiharvest::j, metsig::jet, jetfilter_cfi::jetCollection, jetToken_, findQualityFiles::jj, leadingjet_eta_after_, leadingjet_eta_afterZ_, leadingjet_eta_before_, leadingjet_pt_after_, leadingjet_pt_afterZ_, leadingjet_pt_before_, LogTrace, M_PI, BTaggingMonitor_cfi::met, met_after_, met_afterZ_, met_before_, susyDQM_cfi::metCollection, metMax_, metMin_, metToken_, mt_after_, mt_before_, mtMax_, mtMin_, amptDefaultParameters_cff::mu, muoncharge_after_, muoncharge_afterZ_, muoncharge_before_, pdwgLeptonRecoSkim_cfi::muonCollection, muonHitsCut_, muonTag_, n_zselPt1thr_, n_zselPt2thr_, nall, ngoodmuons_, nJetMax_, BTaggingMonitoring_cff::njets, njets_after_, njets_afterZ_, njets_before_, nMatchesCut_, nmuons_, TrackCollections2monitor_cff::normalizedChi2, normalizedChi2Cut_, nph_, npvs_after_, npvs_afterZ_, npvs_before_, reco::Muon::numberOfMatches(), reco::LeafCandidate::p(), phEta_, phoTag_, ExoticaDQM_cfi::photonCollection, HLTObjectMonitor_cfi::photonPt, phPt_, InitialStepPreSplitting_cff::pixelHits, pixelHitsCut_, reco::BeamSpot::position(), HLTPrescaleProvider::prescaleSet(), DiDispStaMuonMonitor_cfi::pt, reco::LeafCandidate::pt(), pt1_afterZ_, HLT_2024v11_cff::pt2, pt2_afterZ_, pt_after_, pt_before_, ptCut_, ptDiffPM_afterZ_, ptDiffPM_before_, ptThrForPhoton_, ptThrForZ1_, ptThrForZ2_, reco::LeafCandidate::px(), reco::LeafCandidate::py(), reco::LeafCandidate::pz(), quality, mathSSE::sqrt(), reco::MuonIsolation::sumPt, trackerHitsCut_, Geom::Phi< T1, Range >::value(), bphysicsOniaDQM_cfi::vertex, spclusmultinvestigator_cfi::vertexCollection, and vertexTag_.

                                                              {
  // Muon collection
  Handle<View<Muon> > muonCollection;
  if (!ev.getByToken(muonTag_, muonCollection)) {
    // LogWarning("") << ">>> Muon collection does not exist !!!";
    return;
  }
  unsigned int muonCollectionSize = muonCollection->size();

  // Beam spot
  Handle<reco::BeamSpot> beamSpotHandle;
  if (!ev.getByToken(beamSpotTag_, beamSpotHandle)) {
    // LogWarning("") << ">>> No beam spot found !!!";
    return;
  }

  // Loop to reject/control Z->mumu is done separately
  unsigned int nmuonsForZ1 = 0;
  unsigned int nmuonsForZ2 = 0;
  bool cosmic = false;
  for (unsigned int i = 0; i < muonCollectionSize; i++) {
    const Muon& mu = muonCollection->at(i);
    if (!mu.isGlobalMuon())
      continue;
    double pt = mu.pt();
    double dxy = mu.innerTrack()->dxy(beamSpotHandle->position());

    if (fabs(dxy) > 1) {
      cosmic = true;
      break;
    }

    if (pt > ptThrForZ1_)
      nmuonsForZ1++;
    if (pt > ptThrForZ2_)
      nmuonsForZ2++;

    for (unsigned int j = i + 1; j < muonCollectionSize; j++) {
      const Muon& mu2 = muonCollection->at(j);
      if (mu2.isGlobalMuon() && (mu.charge() * mu2.charge() == -1)) {
        const math::XYZTLorentzVector ZRecoGlb(
            mu.px() + mu2.px(), mu.py() + mu2.py(), mu.pz() + mu2.pz(), mu.p() + mu2.p());
        dimuonmass_before_->Fill(ZRecoGlb.mass());
        if (mu.charge() > 0) {
          ptDiffPM_before_->Fill(mu.pt() - mu2.pt());
        } else {
          ptDiffPM_before_->Fill(mu2.pt() - mu.pt());
        }
      }
    }
  }
  if (cosmic)
    return;

  LogTrace("") << "> Z rejection: muons above " << ptThrForZ1_ << " [GeV]: " << nmuonsForZ1;
  LogTrace("") << "> Z rejection: muons above " << ptThrForZ2_ << " [GeV]: " << nmuonsForZ2;

  // MET
  Handle<View<MET> > metCollection;
  if (!ev.getByToken(metToken_, metCollection)) {
    // LogWarning("") << ">>> MET collection does not exist !!!";
    return;
  }
  const MET& met = metCollection->at(0);
  double met_et = met.pt();
  LogTrace("") << ">>> MET, MET_px, MET_py: " << met_et << ", " << met.px() << ", " << met.py() << " [GeV]";
  met_before_->Fill(met_et);

  // Vertices in the event
  Handle<View<reco::Vertex> > vertexCollection;
  if (!ev.getByToken(vertexTag_, vertexCollection)) {
    LogError("") << ">>> Vertex collection does not exist !!!";
    return;
  }
  unsigned int vertexCollectionSize = vertexCollection->size();

  int nvvertex = 0;
  for (unsigned int i = 0; i < vertexCollectionSize; i++) {
    const Vertex& vertex = vertexCollection->at(i);
    if (vertex.isValid())
      nvvertex++;
  }

  npvs_before_->Fill(nvvertex);

  // bool trigger_fired = false;
  //Handle<TriggerResults> triggerResults;
  // if (!ev.getByToken(trigTag_, triggerResults)) {
  // LogWarning("") << ">>> TRIGGER collection does not exist !!!";
  // return;
  // }
  // const edm::TriggerNames& trigNames = ev.triggerNames(*triggerResults);
  //  LogWarning("")<<"Loop over triggers";

  //HLTConfigProvider const&  hltConfigProvider = hltPrescaleProvider_.hltConfigProvider();

  /*  change faulty logic of triggering
  for (unsigned int i=0; i<triggerResults->size(); i++)
  {
          const std::string trigName = trigNames.triggerName(i);

          bool found=false;
          for(unsigned int index=0; index<trigPathNames_.size() && found==false;
  index++) {
               size_t trigPath = trigName.find(trigPathNames_[index]); // 0 if
  found, pos if not
               if (trigPath==0) found=true;
          }
          if(!found) {continue;}

          bool prescaled=false;
          for (unsigned int ps= 0; ps<  hltConfigProvider.prescaleSize();
  ps++){
              const unsigned int prescaleValue =
  hltConfigProvider.prescaleValue(ps, trigName) ;
              if (prescaleValue != 1) prescaled =true;
          }

          if( triggerResults->accept(i) && !prescaled){   trigger_fired=true;}
                    // LogWarning("")<<"TrigNo: "<<i<<"  "<<found<<"
  "<<trigName<<" ---> FIRED";}
  }
  */

  // get the prescale set for this event
  const int prescaleSet = hltPrescaleProvider_.prescaleSet(ev, iSet);
  if (prescaleSet == -1) {
    LogTrace("") << "Failed to determine prescaleSet\n";
    // std::cout << "Failed to determine prescaleSet. Check the GlobalTag in
    // cfg\n";
    return;
  }

  // for (unsigned int i = 0;
  //    (i < triggerResults->size()) && (trigger_fired == false); i++) {
  // skip trigger, if it did not fire
  //if (!triggerResults->accept(i)) continue;

  // skip trigger, if it is not on our list
  //bool found = false;
  //const std::string trigName = trigNames.triggerName(i);
  //for (unsigned int index = 0;
  //   index < trigPathNames_.size() && found == false; index++) {
  // if (trigName.find(trigPathNames_.at(index)) == 0) found = true;
  // }
  // if (!found) continue;

  // skip trigger, if it is prescaled
  /* if (prescaleSet != -1) {
      if (hltConfigProvider.prescaleValue(prescaleSet, trigName) != 1)
        continue;
    } else {
      // prescaleSet is not known.
      // This branch is not needed, if prescaleSet=-1 forces to skip event
      int prescaled = 0;
      for (unsigned int ps = 0;
           !prescaled && (ps < hltConfigProvider.prescaleSize()); ++ps) {
        if (hltConfigProvider.prescaleValue(ps, trigName) != 1) {
          prescaled = 1;
        }
      }
      if (prescaled) {
        // std::cout << "trigger prescaled\n";
        continue;
      }
    }
*/
  // std::cout << "found unprescaled trigger that fired: " << trigName <<
  // "\n";
  // trigger_fired = true;
  // }
  // if (trigger_fired) std::cout << "\n\tGot Trigger\n";

  // trig_before_->Fill(trigger_fired);

  // Jet collection
  Handle<View<Jet> > jetCollection;
  if (!ev.getByToken(jetToken_, jetCollection)) {
    // LogError("") << ">>> JET collection does not exist !!!";
    return;
  }
  unsigned int jetCollectionSize = jetCollection->size();
  int njets = 0;
  int LEADJET = -1;
  double max_pt = 0;
  for (unsigned int i = 0; i < jetCollectionSize; i++) {
    const Jet& jet = jetCollection->at(i);
    double minDistance = 99999;  // This is in order to use PFJets
    for (unsigned int j = 0; j < muonCollectionSize; j++) {
      const Muon& mu = muonCollection->at(j);
      double distance =
          sqrt((mu.eta() - jet.eta()) * (mu.eta() - jet.eta()) + (mu.phi() - jet.phi()) * (mu.phi() - jet.phi()));
      if (minDistance > distance)
        minDistance = distance;
    }
    if (minDistance < 0.3)
      continue;  // 0.3 is the isolation cone around the muon
    if (jet.et() > max_pt) {
      LEADJET = i;
      max_pt = jet.et();
    }
    if (jet.et() > eJetMin_) {
      njets++;
    }
  }

  LogTrace("") << ">>> Total number of jets: " << jetCollectionSize;
  LogTrace("") << ">>> Number of jets above " << eJetMin_ << " [GeV]: " << njets;
  njets_before_->Fill(njets);
  double lead_jet_pt = -1;
  double lead_jet_eta = -100;
  if (LEADJET != -1) {
    const Jet& leadJet = jetCollection->at(LEADJET);
    leadingjet_pt_before_->Fill(leadJet.pt());
    leadingjet_eta_before_->Fill(leadJet.eta());
    lead_jet_pt = leadJet.pt();
    lead_jet_eta = leadJet.eta();
  }
  // Photon Collection
  Handle<View<Photon> > photonCollection;
  if (!ev.getByToken(phoTag_, photonCollection)) {
    // LogError("")
    return;
  }
  unsigned int ngam = 0;

  for (unsigned int i = 0; i < photonCollection->size(); i++) {
    const Photon& ph = photonCollection->at(i);
    double photonPt = ph.pt();
    if (photonPt > ptThrForPhoton_) {
      ngam++;
      phEta_->Fill(ph.eta());
    }
    phPt_->Fill(photonPt);
  }
  nph_->Fill(ngam);
  LogTrace("") << " >>> N photons " << ngam << std::endl;

  nmuons_->Fill(muonCollectionSize);

  // Start counting
  nall++;

  // Histograms per event should be done only once, so keep track of them
  // bool hlt_hist_done = false;
  // bool zhlt_hist_done = false;
  bool zjets_hist_done = false;
  bool zfullsel_hist_done = false;
  bool met_hist_done = false;
  bool njets_hist_done = false;
  bool wfullsel_hist_done = false;

  // Central W->mu nu selection criteria
  const int NFLAGS = 10;
  bool muon_sel[NFLAGS];
  const int NFLAGSZ = 12;
  bool zmuon_sel[NFLAGSZ];
  bool muon4Z = false;

  double number_of_goodMuons = 0;

  for (unsigned int i = 0; i < muonCollectionSize; i++) {
    for (int j = 0; j < NFLAGS; ++j) {
      muon_sel[j] = false;
    }

    const Muon& mu = muonCollection->at(i);
    if (!mu.isGlobalMuon())
      continue;
    if (mu.globalTrack().isNull())
      continue;
    if (mu.innerTrack().isNull())
      continue;

    LogTrace("") << "> Wsel: processing muon number " << i << "...";
    reco::TrackRef gm = mu.globalTrack();
    reco::TrackRef tk = mu.innerTrack();

    // Pt,eta cuts
    double pt = mu.pt();
    double eta = mu.eta();
    LogTrace("") << "\t... pt, eta: " << pt << " [GeV], " << eta;
    ;
    if (pt > ptCut_)
      muon_sel[0] = true;
    if (fabs(eta) < etaCut_)
      muon_sel[1] = true;

    double charge = mu.charge();

    // d0, chi2, nhits quality cuts
    double dxy = gm->dxy(beamSpotHandle->position());
    double normalizedChi2 = gm->normalizedChi2();
    double trackerHits = tk->hitPattern().numberOfValidTrackerHits();
    int pixelHits = tk->hitPattern().numberOfValidPixelHits();
    int muonHits = gm->hitPattern().numberOfValidMuonHits();
    int nMatches = mu.numberOfMatches();

    LogTrace("") << "\t... dxy, normalizedChi2, trackerHits, isTrackerMuon?: " << dxy << " [cm], " << normalizedChi2
                 << ", " << trackerHits << ", " << mu.isTrackerMuon();
    if (fabs(dxy) < dxyCut_)
      muon_sel[2] = true;

    bool quality = true;

    if (normalizedChi2 > normalizedChi2Cut_)
      quality = false;
    if (trackerHits < trackerHitsCut_)
      quality = false;
    if (pixelHits < pixelHitsCut_)
      quality = false;
    if (muonHits < muonHitsCut_)
      quality = false;
    ;
    if (!mu.isTrackerMuon())
      quality = false;
    if (nMatches < nMatchesCut_)
      quality = false;
    muon_sel[3] = quality;
    if (quality)
      number_of_goodMuons++;

    pt_before_->Fill(pt);
    eta_before_->Fill(eta);
    dxy_before_->Fill(dxy);
    muoncharge_before_->Fill(charge);
    goodewkmuon_before_->Fill(quality);

    // Charge asymmetry
    // if (quality) {
    //  if (charge>0) ptPlus_before_->Fill(pt);
    //  if (charge<0) ptMinus_before_->Fill(pt);
    //}

    // Isolation cuts
    double isovar = mu.isolationR03().sumPt;
    if (isCombinedIso_) {
      isovar += mu.isolationR03().emEt;
      isovar += mu.isolationR03().hadEt;
    }
    if (isRelativeIso_)
      isovar /= pt;
    if (isovar < isoCut03_)
      muon_sel[4] = true;

    LogTrace("") << "\t... isolation value" << isovar << ", isolated? " << muon_sel[6];
    iso_before_->Fill(isovar);

    // HLT (not mtched to muon for the time being)
    // if (trigger_fired) muon_sel[5] = true;

    // For Z:
    if (pt > ptThrForZ1_ && fabs(eta) < etaCut_ && fabs(dxy) < dxyCut_ && quality && isovar < isoCut03_) {
      muon4Z = true;
    }

    // MET/MT cuts
    double w_et = met_et + mu.pt();
    double w_px = met.px() + mu.px();
    double w_py = met.py() + mu.py();

    double massT = w_et * w_et - w_px * w_px - w_py * w_py;
    massT = (massT > 0) ? sqrt(massT) : 0;

    LogTrace("") << "\t... W mass, W_et, W_px, W_py: " << massT << ", " << w_et << ", " << w_px << ", " << w_py
                 << " [GeV]";
    if (massT > mtMin_ && massT < mtMax_)
      muon_sel[5] = true;
    mt_before_->Fill(massT);
    if (met_et > metMin_ && met_et < metMax_)
      muon_sel[6] = true;

    // Acoplanarity cuts
    Geom::Phi<double> deltaphi(mu.phi() - atan2(met.py(), met.px()));
    double acop = deltaphi.value();
    if (acop < 0)
      acop = -acop;
    acop = M_PI - acop;
    LogTrace("") << "\t... acoplanarity: " << acop;
    if (acop < acopCut_)
      muon_sel[7] = true;
    acop_before_->Fill(acop);

    // Remaining flags (from global event information)
    if (nmuonsForZ1 < 1 || nmuonsForZ2 < 2)
      muon_sel[8] = true;
    if (njets <= nJetMax_)
      muon_sel[9] = true;

    // Collect necessary flags "per muon"
    int flags_passed = 0;
    for (int j = 0; j < NFLAGS; ++j) {
      if (muon_sel[j])
        flags_passed += 1;
    }

    // Do N-1 histograms now (and only once for global event quantities)
    if (flags_passed >= (NFLAGS - 1)) {
      if (!muon_sel[0] || flags_passed == NFLAGS)
        pt_after_->Fill(pt);
      if (!muon_sel[1] || flags_passed == NFLAGS)
        eta_after_->Fill(eta);
      if (!muon_sel[2] || flags_passed == NFLAGS)
        dxy_after_->Fill(dxy);
      if (!muon_sel[3] || flags_passed == NFLAGS)
        goodewkmuon_after_->Fill(quality);
      if (!muon_sel[4] || flags_passed == NFLAGS)
        iso_after_->Fill(isovar);
      // if (!muon_sel[5] || flags_passed == NFLAGS)
      //   if (!hlt_hist_done) trig_after_->Fill(trigger_fired);
      //hlt_hist_done = true;
      if (!muon_sel[5] || flags_passed == NFLAGS)
        mt_after_->Fill(massT);
      if (!muon_sel[6] || flags_passed == NFLAGS)
        if (!met_hist_done)
          met_after_->Fill(met_et);
      met_hist_done = true;
      if (!muon_sel[7] || flags_passed == NFLAGS)
        acop_after_->Fill(acop);
      // no action here for muon_sel[8]
      if (!muon_sel[9] || flags_passed == NFLAGS) {
        if (!njets_hist_done) {
          njets_after_->Fill(njets);
          leadingjet_pt_after_->Fill(lead_jet_pt);
          leadingjet_eta_after_->Fill(lead_jet_eta);
        }
        njets_hist_done = true;
      }
      if (flags_passed == NFLAGS) {
        if (!wfullsel_hist_done) {
          npvs_after_->Fill(nvvertex);
          muoncharge_after_->Fill(charge);
          // if (charge>0) ptPlus_afterW_->Fill(pt);
          // if (charge<0) ptMinus_afterW_->Fill(pt);
        }
        wfullsel_hist_done = true;
      }
    }

    // The cases in which the event is rejected as a Z are considered
    // independently:
    if (muon4Z && !muon_sel[8]) {
      // Plots for 2 muons
      for (unsigned int j = i + 1; j < muonCollectionSize; j++) {
        for (int ij = 0; ij < NFLAGSZ; ++ij) {
          zmuon_sel[ij] = false;
        }

        for (int ji = 0; ji < 5; ++ji) {
          zmuon_sel[ji] = muon_sel[ji];
        }

        const Muon& mu2 = muonCollection->at(j);
        if (!mu2.isGlobalMuon())
          continue;
        if (mu2.charge() * charge != -1)
          continue;
        reco::TrackRef gm2 = mu2.globalTrack();
        reco::TrackRef tk2 = mu2.innerTrack();
        double pt2 = mu2.pt();
        if (pt2 > ptThrForZ2_)
          zmuon_sel[5] = true;
        double eta2 = mu2.eta();
        if (fabs(eta2) < etaCut_)
          zmuon_sel[6] = true;
        double dxy2 = gm2->dxy(beamSpotHandle->position());
        if (fabs(dxy2) < dxyCut_)
          zmuon_sel[7] = true;
        double normalizedChi22 = gm2->normalizedChi2();
        double trackerHits2 = tk2->hitPattern().numberOfValidTrackerHits();
        int pixelHits2 = tk2->hitPattern().numberOfValidPixelHits();
        int muonHits2 = gm2->hitPattern().numberOfValidMuonHits();
        int nMatches2 = mu2.numberOfMatches();
        bool quality2 = true;
        if (normalizedChi22 > normalizedChi2Cut_)
          quality2 = false;
        if (trackerHits2 < trackerHitsCut_)
          quality2 = false;
        if (pixelHits2 < pixelHitsCut_)
          quality2 = false;
        if (muonHits2 < muonHitsCut_)
          quality2 = false;
        if (!mu2.isTrackerMuon())
          quality2 = false;
        if (nMatches2 < nMatchesCut_)
          quality2 = false;
        zmuon_sel[8] = quality2;
        double isovar2 = mu2.isolationR03().sumPt;
        if (isCombinedIso_) {
          isovar2 += mu2.isolationR03().emEt;
          isovar2 += mu2.isolationR03().hadEt;
        }
        if (isRelativeIso_)
          isovar2 /= pt2;
        if (isovar2 < isoCut03_)
          zmuon_sel[9] = true;
        //  if (trigger_fired) zmuon_sel[10] = true;
        const math::XYZTLorentzVector ZRecoGlb(
            mu.px() + mu2.px(), mu.py() + mu2.py(), mu.pz() + mu2.pz(), mu.p() + mu2.p());
        if (ZRecoGlb.mass() > dimuonMassMin_ && ZRecoGlb.mass() < dimuonMassMax_)
          zmuon_sel[10] = true;

        // jet flag
        if (njets <= nJetMax_)
          zmuon_sel[11] = true;

        // start filling histos: N-1 plots
        int flags_passed_z = 0;

        for (int jj = 0; jj < NFLAGSZ; ++jj) {
          if (zmuon_sel[jj])
            ++flags_passed_z;
        }

        if (flags_passed_z >= (NFLAGSZ - 1)) {
          if (!zmuon_sel[0] || flags_passed_z == NFLAGSZ) {
            pt1_afterZ_->Fill(pt);
          }
          if (!zmuon_sel[1] || flags_passed_z == NFLAGSZ) {
            eta1_afterZ_->Fill(eta);
          }
          if (!zmuon_sel[2] || flags_passed_z == NFLAGSZ) {
            dxy1_afterZ_->Fill(dxy);
          }
          if (!zmuon_sel[3] || flags_passed_z == NFLAGSZ) {
            goodewkmuon1_afterZ_->Fill(quality);
          }
          if (!zmuon_sel[4] || flags_passed_z == NFLAGSZ) {
            iso1_afterZ_->Fill(isovar);
          }
          if (!zmuon_sel[5] || flags_passed_z == NFLAGSZ) {
            pt2_afterZ_->Fill(pt2);
          }
          if (!zmuon_sel[6] || flags_passed_z == NFLAGSZ) {
            eta2_afterZ_->Fill(eta2);
          }
          if (!zmuon_sel[7] || flags_passed_z == NFLAGSZ) {
            dxy2_afterZ_->Fill(dxy2);
          }
          if (!zmuon_sel[8] || flags_passed_z == NFLAGSZ) {
            goodewkmuon2_afterZ_->Fill(quality2);
          }
          if (!zmuon_sel[9] || flags_passed_z == NFLAGSZ) {
            iso2_afterZ_->Fill(isovar2);
          }
          //     if (!zmuon_sel[10] || flags_passed_z == NFLAGSZ) {
          //     if (!zhlt_hist_done) ztrig_afterZ_->Fill(trigger_fired);
          //   zhlt_hist_done = true;
          //  }
          if (!zmuon_sel[10] || flags_passed_z == NFLAGSZ) {
            dimuonmass_afterZ_->Fill(ZRecoGlb.mass());
          }
          if (!zmuon_sel[11] || flags_passed_z == NFLAGSZ) {
            if (!zjets_hist_done) {
              njets_afterZ_->Fill(njets);
              leadingjet_pt_afterZ_->Fill(lead_jet_pt);
              leadingjet_eta_afterZ_->Fill(lead_jet_eta);
            }
            zjets_hist_done = true;
          }
          if (flags_passed_z == NFLAGSZ) {
            met_afterZ_->Fill(met_et);
            if (!zfullsel_hist_done) {
              npvs_afterZ_->Fill(nvvertex);
              muoncharge_afterZ_->Fill(charge);
              if (charge > 0) {
                // ptPlus_afterZ_->Fill(mu.pt());
                // ptMinus_afterZ_->Fill(mu2.pt());
                ptDiffPM_afterZ_->Fill(mu.pt() - mu2.pt());
              } else {
                // ptPlus_afterZ_->Fill(mu2.pt());
                // ptMinus_afterZ_->Fill(mu.pt());
                ptDiffPM_afterZ_->Fill(mu2.pt() - mu.pt());
              }
            }
            zfullsel_hist_done = true;
          }
        }
      }
    }
  }

  if (zfullsel_hist_done) {
    // here was a Z candidate
    n_zselPt1thr_->Fill(nmuonsForZ1);
    n_zselPt2thr_->Fill(nmuonsForZ2);
  }

  // nmuons_->Fill(muonCollectionSize);
  ngoodmuons_->Fill(number_of_goodMuons);

  return;
}

bookHistograms()

void EwkMuDQM::bookHistograms ( DQMStore::IBooker &  ibooker, edm::Run const &  , edm::EventSetup const &    )

overrideprotectedvirtual

\ For Z-boson events

Implements DQMEDAnalyzer.

Definition at line 117 of file EwkMuDQM.cc.

References acop_after_, acop_before_, dqm::implementation::IBooker::book1D(), dimuonmass_afterZ_, dimuonmass_before_, dxy1_afterZ_, dxy2_afterZ_, dxy_after_, dxy_before_, eJetMin_, eta1_afterZ_, eta2_afterZ_, eta_after_, eta_before_, goodewkmuon1_afterZ_, goodewkmuon2_afterZ_, goodewkmuon_after_, goodewkmuon_before_, isCombinedIso_, iso1_afterZ_, iso2_afterZ_, iso_after_, iso_before_, isRelativeIso_, jetTag_, edm::InputTag::label(), leadingjet_eta_after_, leadingjet_eta_afterZ_, leadingjet_eta_before_, leadingjet_pt_after_, leadingjet_pt_afterZ_, leadingjet_pt_before_, M_PI, met_after_, met_afterZ_, met_before_, metTag_, mt_after_, mt_before_, muoncharge_after_, muoncharge_afterZ_, muoncharge_before_, n_zselPt1thr_, n_zselPt2thr_, ngoodmuons_, njets_after_, njets_afterZ_, njets_before_, nmuons_, nph_, npvs_after_, npvs_afterZ_, npvs_before_, phEta_, phPt_, pt1_afterZ_, pt2_afterZ_, pt_after_, pt_before_, ptDiffPM_afterZ_, ptDiffPM_before_, ptThrForPhoton_, ptThrForZ1_, ptThrForZ2_, and dqm::implementation::NavigatorBase::setCurrentFolder().

                                                                                           {
  ibooker.setCurrentFolder("Physics/EwkMuDQM");

  char chtitle[256] = "";

  pt_before_ = ibooker.book1D("PT_BEFORECUTS", "Muon transverse momentum (global muon) [GeV]", 100, 0., 100.);
  pt_after_ = ibooker.book1D("PT_AFTERWCUTS", "Muon transverse momentum (global muon) [GeV]", 100, 0., 100.);

  eta_before_ = ibooker.book1D("ETA_BEFORECUTS", "Muon pseudo-rapidity", 50, -2.5, 2.5);
  eta_after_ = ibooker.book1D("ETA_AFTERWCUTS", "Muon pseudo-rapidity", 50, -2.5, 2.5);

  dxy_before_ = ibooker.book1D("DXY_BEFORECUTS", "Muon transverse distance to beam spot [cm]", 100, -0.5, 0.5);
  dxy_after_ = ibooker.book1D("DXY_AFTERWCUTS", "Muon transverse distance to beam spot [cm]", 100, -0.5, 0.5);

  goodewkmuon_before_ = ibooker.book1D("GOODEWKMUON_BEFORECUTS", "Quality-muon flag", 2, -0.5, 1.5);
  goodewkmuon_after_ = ibooker.book1D("GOODEWKMUON_AFTERWCUTS", "Quality-muon flag", 2, -0.5, 1.5);

  if (isRelativeIso_) {
    if (isCombinedIso_) {
      iso_before_ = ibooker.book1D("ISO_BEFORECUTS", "Relative (combined) isolation variable", 100, 0., 1.);
      iso_after_ = ibooker.book1D("ISO_AFTERWCUTS", "Relative (combined) isolation variable", 100, 0., 1.);
    } else {
      iso_before_ = ibooker.book1D("ISO_BEFORECUTS", "Relative (tracker) isolation variable", 100, 0., 1.);
      iso_after_ = ibooker.book1D("ISO_AFTERWCUTS", "Relative (tracker) isolation variable", 100, 0., 1.);
    }
  } else {
    if (isCombinedIso_) {
      iso_before_ = ibooker.book1D("ISO_BEFORECUTS", "Absolute (combined) isolation variable [GeV]", 100, 0., 20.);
      iso_after_ = ibooker.book1D("ISO_AFTERWCUTS", "Absolute (combined) isolation variable [GeV]", 100, 0., 20.);
    } else {
      iso_before_ = ibooker.book1D("ISO_BEFORECUTS", "Absolute (tracker) isolation variable [GeV]", 100, 0., 20.);
      iso_after_ = ibooker.book1D("ISO_AFTERWCUTS", "Absolute (tracker) isolation variable [GeV]", 100, 0., 20.);
    }
  }

  /*  trig_before_ = ibooker.book1D("TRIG_BEFORECUTS",
      "Trigger response (boolean of muon triggers)", 2, -0.5, 1.5);
  trig_after_ = ibooker.book1D("TRIG_AFTERWCUTS",
      "Trigger response (boolean of muon triggers)", 2, -0.5, 1.5);
*/
  snprintf(chtitle, 255, "Transverse mass (%s) [GeV]", metTag_.label().data());
  mt_before_ = ibooker.book1D("MT_BEFORECUTS", chtitle, 150, 0., 300.);
  mt_after_ = ibooker.book1D("MT_AFTERWCUTS", chtitle, 150, 0., 300.);

  snprintf(chtitle, 255, "Missing transverse energy (%s) [GeV]", metTag_.label().data());
  met_before_ = ibooker.book1D("MET_BEFORECUTS", chtitle, 100, 0., 200.);
  met_after_ = ibooker.book1D("MET_AFTERWCUTS", chtitle, 100, 0., 200.);
  met_afterZ_ = ibooker.book1D("MET_AFTERZCUTS", chtitle, 100, 0., 200.);

  snprintf(chtitle, 255, "MU-MET (%s) acoplanarity", metTag_.label().data());
  acop_before_ = ibooker.book1D("ACOP_BEFORECUTS", chtitle, 50, 0., M_PI);
  acop_after_ = ibooker.book1D("ACOP_AFTERWCUTS", chtitle, 50, 0., M_PI);

  snprintf(chtitle, 255, "Z selection: muons above %.2f GeV", ptThrForZ1_);
  n_zselPt1thr_ = ibooker.book1D("NZSELPT1THR", chtitle, 10, -0.5, 9.5);
  snprintf(chtitle, 255, "Z selection: muons above %.2f GeV", ptThrForZ2_);
  n_zselPt2thr_ = ibooker.book1D("NZSELPT2THR", chtitle, 10, -0.5, 9.5);

  snprintf(chtitle, 255, "Number of jets (%s) above %.2f GeV", jetTag_.label().data(), eJetMin_);
  njets_before_ = ibooker.book1D("NJETS_BEFORECUTS", chtitle, 16, -0.5, 15.5);
  njets_after_ = ibooker.book1D("NJETS_AFTERWCUTS", chtitle, 16, -0.5, 15.5);
  njets_afterZ_ = ibooker.book1D("NJETS_AFTERZCUTS", chtitle, 16, -0.5, 15.5);

  leadingjet_pt_before_ = ibooker.book1D("LEADINGJET_PT_BEFORECUTS", "Leading Jet transverse momentum", 300, 0., 300.);
  leadingjet_pt_after_ = ibooker.book1D("LEADINGJET_PT_AFTERWCUTS", "Leading Jet transverse momentum", 300, 0., 300.);
  leadingjet_pt_afterZ_ = ibooker.book1D("LEADINGJET_PT_AFTERZCUTS", "Leading Jet transverse momentum", 300, 0., 300.);

  leadingjet_eta_before_ = ibooker.book1D("LEADINGJET_ETA_BEFORECUTS", "Leading Jet pseudo-rapidity", 50, -2.5, 2.5);
  leadingjet_eta_after_ = ibooker.book1D("LEADINGJET_ETA_AFTERWCUTS", "Leading Jet pseudo-rapidity", 50, -2.5, 2.5);
  leadingjet_eta_afterZ_ = ibooker.book1D("LEADINGJET_ETA_AFTERZCUTS", "Leading Jet pseudo-rapidity", 50, -2.5, 2.5);

  ptDiffPM_before_ = ibooker.book1D("PTDIFFPM_BEFORE_CUTS", "pt(Muon+)-pt(Muon-) after Z cuts [GeV]", 200, -100., 100.);
  ptDiffPM_afterZ_ = ibooker.book1D("PTDIFFPM_AFTERZ_CUTS", "pt(Muon+)-pt(Muon-) after Z cuts [GeV]", 200, -100., 100.);

  pt1_afterZ_ = ibooker.book1D("PT1_AFTERZCUTS", "Muon transverse momentum (global muon) [GeV]", 100, 0., 100.);
  eta1_afterZ_ = ibooker.book1D("ETA1_AFTERZCUTS", "Muon pseudo-rapidity", 50, -2.5, 2.5);
  dxy1_afterZ_ = ibooker.book1D("DXY1_AFTERZCUTS", "Muon transverse distance to beam spot [cm]", 100, -0.5, 0.5);
  goodewkmuon1_afterZ_ = ibooker.book1D("GOODEWKMUON1_AFTERZCUTS", "Quality-muon flag", 2, -0.5, 1.5);

  if (isRelativeIso_) {
    if (isCombinedIso_) {
      iso1_afterZ_ = ibooker.book1D("ISO1_AFTERZCUTS", "Relative (combined) isolation variable", 100, 0., 1.);
      iso2_afterZ_ = ibooker.book1D("ISO2_AFTERZCUTS", "Relative (combined) isolation variable", 100, 0., 1.);
    } else {
      iso1_afterZ_ = ibooker.book1D("ISO1_AFTERZCUTS", "Relative (tracker) isolation variable", 100, 0., 1.);
      iso2_afterZ_ = ibooker.book1D("ISO2_AFTERZCUTS", "Relative (tracker) isolation variable", 100, 0., 1.);
    }
  } else {
    if (isCombinedIso_) {
      iso1_afterZ_ = ibooker.book1D("ISO1_AFTERZCUTS", "Absolute (combined) isolation variable [GeV]", 100, 0., 20.);
      iso2_afterZ_ = ibooker.book1D("ISO2_AFTERZCUTS", "Absolute (combined) isolation variable [GeV]", 100, 0., 20.);
    } else {
      iso1_afterZ_ = ibooker.book1D("ISO1_AFTERZCUTS", "Absolute (tracker) isolation variable [GeV]", 100, 0., 20.);
      iso2_afterZ_ = ibooker.book1D("ISO2_AFTERZCUTS", "Absolute (tracker) isolation variable [GeV]", 100, 0., 20.);
    }
  }

  pt2_afterZ_ = ibooker.book1D("PT2_AFTERZCUTS", "Muon transverse momentum (global muon) [GeV]", 100, 0., 100.);
  eta2_afterZ_ = ibooker.book1D("ETA2_AFTERZCUTS", "Muon pseudo-rapidity", 50, -2.5, 2.5);
  dxy2_afterZ_ = ibooker.book1D("DXY2_AFTERZCUTS", "Muon transverse distance to beam spot [cm]", 100, -0.5, 0.5);
  goodewkmuon2_afterZ_ = ibooker.book1D("GOODEWKMUON2_AFTERZCUTS", "Quality-muon flag", 2, -0.5, 1.5);
  /*  ztrig_afterZ_ = ibooker.book1D("ZTRIG_AFTERZCUTS",
      "Trigger response (boolean of muon triggers)", 2, -0.5, 1.5);
 */
  dimuonmass_before_ = ibooker.book1D("DIMUONMASS_BEFORECUTS", "DiMuonMass (2 globals)", 100, 0, 200);
  dimuonmass_afterZ_ = ibooker.book1D("DIMUONMASS_AFTERZCUTS", "DiMuonMass (2 globals)", 100, 0, 200);
  npvs_before_ = ibooker.book1D("NPVs_BEFORECUTS", "Number of Valid Primary Vertices", 51, -0.5, 50.5);
  npvs_after_ = ibooker.book1D("NPVs_AFTERWCUTS", "Number of Valid Primary Vertices", 51, -0.5, 50.5);
  npvs_afterZ_ = ibooker.book1D("NPVs_AFTERZCUTS", "Number of Valid Primary Vertices", 51, -0.5, 50.5);
  muoncharge_before_ = ibooker.book1D("MUONCHARGE_BEFORECUTS", "Muon Charge", 3, -1.5, 1.5);
  muoncharge_after_ = ibooker.book1D("MUONCHARGE_AFTERWCUTS", "Muon Charge", 3, -1.5, 1.5);
  muoncharge_afterZ_ = ibooker.book1D("MUONCHARGE_AFTERZCUTS", "Muon Charge", 3, -1.5, 1.5);

  // Adding these to replace the NZ ones (more useful, since they are more
  // general?)
  nmuons_ = ibooker.book1D("NMuons", "Number of muons in the event", 10, -0.5, 9.5);
  ngoodmuons_ = ibooker.book1D("NGoodMuons", "Number of muons passing the quality criteria", 10, -0.5, 9.5);

  nph_ = ibooker.book1D("nph", "Number of photons in the event", 20, 0., 20.);
  phPt_ = ibooker.book1D("phPt", "Photon transverse momentum [GeV]", 100, 0., 1000.);
  snprintf(chtitle, 255, "Photon pseudorapidity (pT>%4.1f)", ptThrForPhoton_);
  phEta_ = ibooker.book1D("phEta", chtitle, 100, -2.5, 2.5);
}

dqmBeginRun()

void EwkMuDQM::dqmBeginRun ( const edm::Run &  iRun, const edm::EventSetup &  iSet  )

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 101 of file EwkMuDQM.cc.

References hltPrescaleProvider_, HLTPrescaleProvider::init(), isValidHltConfig_, nall, nhlt, niso, nmet, nrec, nsel, and nz.

                                                                  {
  nall = 0;
  nsel = 0;
  nz = 0;

  nrec = 0;
  niso = 0;
  nhlt = 0;
  nmet = 0;

  // passed as parameter to HLTConfigProvider::init(), not yet used
  bool isConfigChanged = false;
  // isValidHltConfig_ used to short-circuit analyze() in case of problems
  isValidHltConfig_ = hltPrescaleProvider_.init(iRun, iSet, "HLT", isConfigChanged);
}

init_histograms()

void EwkMuDQM::init_histograms ( )

protected

Member Data Documentation

acop_after_

MonitorElement* EwkMuDQM::acop_after_

private

Definition at line 112 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

acop_before_

MonitorElement* EwkMuDQM::acop_before_

private

Definition at line 111 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

acopCut_

double EwkMuDQM::acopCut_

private

Definition at line 63 of file EwkMuDQM.h.

Referenced by analyze().

beamSpotTag_

edm::EDGetTokenT<reco::BeamSpot> EwkMuDQM::beamSpotTag_

private

Definition at line 48 of file EwkMuDQM.h.

Referenced by analyze().

dimuonmass_afterZ_

MonitorElement* EwkMuDQM::dimuonmass_afterZ_

private

Definition at line 151 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

dimuonmass_before_

MonitorElement* EwkMuDQM::dimuonmass_before_

private

Definition at line 150 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

dimuonMassMax_

double EwkMuDQM::dimuonMassMax_

private

Definition at line 76 of file EwkMuDQM.h.

Referenced by analyze().

dimuonMassMin_

double EwkMuDQM::dimuonMassMin_

private

Definition at line 75 of file EwkMuDQM.h.

Referenced by analyze().

dxy1_afterZ_

MonitorElement* EwkMuDQM::dxy1_afterZ_

private

Definition at line 136 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

dxy2_afterZ_

MonitorElement* EwkMuDQM::dxy2_afterZ_

private

Definition at line 141 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

dxy_after_

MonitorElement* EwkMuDQM::dxy_after_

private

Definition at line 100 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

dxy_before_

MonitorElement* EwkMuDQM::dxy_before_

private

Definition at line 99 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

dxyCut_

double EwkMuDQM::dxyCut_

private

Definition at line 52 of file EwkMuDQM.h.

Referenced by analyze().

eJetMin_

double EwkMuDQM::eJetMin_

private

Definition at line 78 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

eta1_afterZ_

MonitorElement* EwkMuDQM::eta1_afterZ_

private

Definition at line 135 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

eta2_afterZ_

MonitorElement* EwkMuDQM::eta2_afterZ_

private

Definition at line 140 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

eta_after_

MonitorElement* EwkMuDQM::eta_after_

private

Definition at line 98 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

eta_before_

MonitorElement* EwkMuDQM::eta_before_

private

Definition at line 97 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

etaCut_

double EwkMuDQM::etaCut_

private

Definition at line 70 of file EwkMuDQM.h.

Referenced by analyze().

goodewkmuon1_afterZ_

MonitorElement* EwkMuDQM::goodewkmuon1_afterZ_

private

Definition at line 137 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

goodewkmuon2_afterZ_

MonitorElement* EwkMuDQM::goodewkmuon2_afterZ_

private

Definition at line 142 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

goodewkmuon_after_

MonitorElement* EwkMuDQM::goodewkmuon_after_

private

Definition at line 102 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

goodewkmuon_before_

MonitorElement* EwkMuDQM::goodewkmuon_before_

private

Definition at line 101 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

hltPrescaleProvider_

HLTPrescaleProvider EwkMuDQM::hltPrescaleProvider_

private

Definition at line 85 of file EwkMuDQM.h.

Referenced by analyze(), and dqmBeginRun().

isAlsoTrackerMuon_

bool EwkMuDQM::isAlsoTrackerMuon_

private

Definition at line 51 of file EwkMuDQM.h.

isCombinedIso_

bool EwkMuDQM::isCombinedIso_

private

Definition at line 60 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

iso1_afterZ_

MonitorElement* EwkMuDQM::iso1_afterZ_

private

Definition at line 138 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

iso2_afterZ_

MonitorElement* EwkMuDQM::iso2_afterZ_

private

Definition at line 143 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

iso_after_

MonitorElement* EwkMuDQM::iso_after_

private

Definition at line 104 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

iso_before_

MonitorElement* EwkMuDQM::iso_before_

private

Definition at line 103 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

isoCut03_

double EwkMuDQM::isoCut03_

private

Definition at line 61 of file EwkMuDQM.h.

Referenced by analyze().

isRelativeIso_

bool EwkMuDQM::isRelativeIso_

private

Definition at line 59 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

isValidHltConfig_

bool EwkMuDQM::isValidHltConfig_

private

Definition at line 84 of file EwkMuDQM.h.

Referenced by dqmBeginRun(), and EwkMuDQM().

jetTag_

edm::InputTag EwkMuDQM::jetTag_

private

Definition at line 41 of file EwkMuDQM.h.

Referenced by bookHistograms().

jetToken_

edm::EDGetTokenT<edm::View<reco::Jet> > EwkMuDQM::jetToken_

private

Definition at line 45 of file EwkMuDQM.h.

Referenced by analyze().

leadingjet_eta_after_

MonitorElement* EwkMuDQM::leadingjet_eta_after_

private

Definition at line 121 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

leadingjet_eta_afterZ_

MonitorElement* EwkMuDQM::leadingjet_eta_afterZ_

private

Definition at line 122 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

leadingjet_eta_before_

MonitorElement* EwkMuDQM::leadingjet_eta_before_

private

Definition at line 120 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

leadingjet_pt_after_

MonitorElement* EwkMuDQM::leadingjet_pt_after_

private

Definition at line 118 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

leadingjet_pt_afterZ_

MonitorElement* EwkMuDQM::leadingjet_pt_afterZ_

private

Definition at line 119 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

leadingjet_pt_before_

MonitorElement* EwkMuDQM::leadingjet_pt_before_

private

Definition at line 117 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

met_after_

MonitorElement* EwkMuDQM::met_after_

private

Definition at line 110 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

met_afterZ_

MonitorElement* EwkMuDQM::met_afterZ_

private

Definition at line 133 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

met_before_

MonitorElement* EwkMuDQM::met_before_

private

Definition at line 109 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

metMax_

double EwkMuDQM::metMax_

private

Definition at line 65 of file EwkMuDQM.h.

Referenced by analyze().

metMin_

double EwkMuDQM::metMin_

private

Definition at line 64 of file EwkMuDQM.h.

Referenced by analyze().

metTag_

edm::InputTag EwkMuDQM::metTag_

private

Definition at line 40 of file EwkMuDQM.h.

Referenced by bookHistograms().

metToken_

edm::EDGetTokenT<edm::View<reco::MET> > EwkMuDQM::metToken_

private

Definition at line 44 of file EwkMuDQM.h.

Referenced by analyze().

mt_after_

MonitorElement* EwkMuDQM::mt_after_

private

Definition at line 108 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

mt_before_

MonitorElement* EwkMuDQM::mt_before_

private

Definition at line 107 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

mtMax_

double EwkMuDQM::mtMax_

private

Definition at line 67 of file EwkMuDQM.h.

Referenced by analyze().

mtMin_

double EwkMuDQM::mtMin_

private

Definition at line 66 of file EwkMuDQM.h.

Referenced by analyze().

muoncharge_after_

MonitorElement* EwkMuDQM::muoncharge_after_

private

Definition at line 158 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

muoncharge_afterZ_

MonitorElement* EwkMuDQM::muoncharge_afterZ_

private

Definition at line 159 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

muoncharge_before_

MonitorElement* EwkMuDQM::muoncharge_before_

private

Definition at line 157 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

muonHitsCut_

int EwkMuDQM::muonHitsCut_

private

Definition at line 56 of file EwkMuDQM.h.

Referenced by analyze().

muonTag_

edm::EDGetTokenT<edm::View<reco::Muon> > EwkMuDQM::muonTag_

private

Definition at line 43 of file EwkMuDQM.h.

Referenced by analyze().

n_zselPt1thr_

MonitorElement* EwkMuDQM::n_zselPt1thr_

private

Definition at line 146 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

n_zselPt2thr_

MonitorElement* EwkMuDQM::n_zselPt2thr_

private

Definition at line 147 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

nall

unsigned int EwkMuDQM::nall

private

Definition at line 87 of file EwkMuDQM.h.

Referenced by analyze(), and dqmBeginRun().

ngoodmuons_

MonitorElement* EwkMuDQM::ngoodmuons_

private

Definition at line 162 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

nhlt

unsigned int EwkMuDQM::nhlt

private

Definition at line 90 of file EwkMuDQM.h.

Referenced by dqmBeginRun().

niso

unsigned int EwkMuDQM::niso

private

Definition at line 89 of file EwkMuDQM.h.

Referenced by dqmBeginRun().

nJetMax_

int EwkMuDQM::nJetMax_

private

Definition at line 79 of file EwkMuDQM.h.

Referenced by analyze().

njets_after_

MonitorElement* EwkMuDQM::njets_after_

private

Definition at line 115 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

njets_afterZ_

MonitorElement* EwkMuDQM::njets_afterZ_

private

Definition at line 116 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

njets_before_

MonitorElement* EwkMuDQM::njets_before_

private

Definition at line 114 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

nMatchesCut_

int EwkMuDQM::nMatchesCut_

private

Definition at line 57 of file EwkMuDQM.h.

Referenced by analyze().

nmet

unsigned int EwkMuDQM::nmet

private

Definition at line 91 of file EwkMuDQM.h.

Referenced by dqmBeginRun().

nmuons_

MonitorElement* EwkMuDQM::nmuons_

private

Definition at line 161 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

normalizedChi2Cut_

double EwkMuDQM::normalizedChi2Cut_

private

Definition at line 53 of file EwkMuDQM.h.

Referenced by analyze().

npfph_

MonitorElement* EwkMuDQM::npfph_

private

Definition at line 164 of file EwkMuDQM.h.

nph_

MonitorElement* EwkMuDQM::nph_

private

Definition at line 165 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

nPhoMax_

int EwkMuDQM::nPhoMax_

private

Definition at line 82 of file EwkMuDQM.h.

npvs_after_

MonitorElement* EwkMuDQM::npvs_after_

private

Definition at line 154 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

npvs_afterZ_

MonitorElement* EwkMuDQM::npvs_afterZ_

private

Definition at line 155 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

npvs_before_

MonitorElement* EwkMuDQM::npvs_before_

private

Definition at line 153 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

nrec

unsigned int EwkMuDQM::nrec

private

Definition at line 88 of file EwkMuDQM.h.

Referenced by dqmBeginRun().

nsel

unsigned int EwkMuDQM::nsel

private

Definition at line 92 of file EwkMuDQM.h.

Referenced by dqmBeginRun().

nz

unsigned int EwkMuDQM::nz

private

Definition at line 93 of file EwkMuDQM.h.

Referenced by dqmBeginRun().

pfphEta_

MonitorElement* EwkMuDQM::pfphEta_

private

Definition at line 168 of file EwkMuDQM.h.

pfphPt_

MonitorElement* EwkMuDQM::pfphPt_

private

Definition at line 166 of file EwkMuDQM.h.

phEta_

MonitorElement* EwkMuDQM::phEta_

private

Definition at line 169 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

phoTag_

edm::EDGetTokenT<edm::View<reco::Photon> > EwkMuDQM::phoTag_

private

Definition at line 46 of file EwkMuDQM.h.

Referenced by analyze().

phPt_

MonitorElement* EwkMuDQM::phPt_

private

Definition at line 167 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

pixelHitsCut_

int EwkMuDQM::pixelHitsCut_

private

Definition at line 55 of file EwkMuDQM.h.

Referenced by analyze().

pt1_afterZ_

MonitorElement* EwkMuDQM::pt1_afterZ_

private

Definition at line 134 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

pt2_afterZ_

MonitorElement* EwkMuDQM::pt2_afterZ_

private

Definition at line 139 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

pt_after_

MonitorElement* EwkMuDQM::pt_after_

private

Definition at line 96 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

pt_before_

MonitorElement* EwkMuDQM::pt_before_

private

Definition at line 95 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

ptCut_

double EwkMuDQM::ptCut_

private

Definition at line 69 of file EwkMuDQM.h.

Referenced by analyze(), and badGlobalMuonTaggerFWLite.BadGlobalMuonTagger::badMuons().

ptDiffPM_afterZ_

MonitorElement* EwkMuDQM::ptDiffPM_afterZ_

private

Definition at line 131 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

ptDiffPM_before_

MonitorElement* EwkMuDQM::ptDiffPM_before_

private

Definition at line 126 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

ptThrForPhoton_

double EwkMuDQM::ptThrForPhoton_

private

Definition at line 81 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

ptThrForZ1_

double EwkMuDQM::ptThrForZ1_

private

Definition at line 72 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

ptThrForZ2_

double EwkMuDQM::ptThrForZ2_

private

Definition at line 73 of file EwkMuDQM.h.

Referenced by analyze(), and bookHistograms().

trackerHitsCut_

int EwkMuDQM::trackerHitsCut_

private

Definition at line 54 of file EwkMuDQM.h.

Referenced by analyze().

trig_after_

MonitorElement* EwkMuDQM::trig_after_

private

Definition at line 106 of file EwkMuDQM.h.

trig_before_

MonitorElement* EwkMuDQM::trig_before_

private

Definition at line 105 of file EwkMuDQM.h.

trigPathNames_

std::vector<std::string> EwkMuDQM::trigPathNames_

private

Definition at line 49 of file EwkMuDQM.h.

trigTag_

edm::EDGetTokenT<edm::TriggerResults> EwkMuDQM::trigTag_

private

Definition at line 42 of file EwkMuDQM.h.

vertexTag_

edm::EDGetTokenT<edm::View<reco::Vertex> > EwkMuDQM::vertexTag_

private

Definition at line 47 of file EwkMuDQM.h.

Referenced by analyze().

ztrig_afterZ_

MonitorElement* EwkMuDQM::ztrig_afterZ_

private

Definition at line 149 of file EwkMuDQM.h.