#include <EwkMuLumiMonitorDQM.h>

Inheritance diagram for EwkMuLumiMonitorDQM:

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

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

void	dqmBeginRun (const edm::Run &, const edm::EventSetup &) override

	EwkMuLumiMonitorDQM (const edm::ParameterSet &)

void	init_histograms ()

bool	IsMuMatchedToHLTMu (const reco::Muon &, const std::vector< reco::Particle > &, double, double)

double	muIso (const reco::Muon &)

double	tkIso (const reco::Track &, edm::Handle< reco::TrackCollection >, edm::Handle< CaloTowerCollection >)

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

Private Attributes
double	acopCut_

edm::EDGetTokenT< reco::BeamSpot >	beamSpotToken_

edm::EDGetTokenT< CaloTowerCollection >	caloTowerToken_

double	deltaRTrk_

double	deltaRVetoTrk_

double	dxyCut_

double	etaMuCut_

MonitorElement *	highMass1HLT_

MonitorElement *	highMass2HLT_

MonitorElement *	highMassGlbSta_

MonitorElement *	highMassGlbTrk_

MonitorElement *	highMassIsBothGlbTrkThanW_

MonitorElement *	highMassNotIso_

HLTConfigProvider	hltConfigProvider_

bool	isCombinedIso_

double	isoCut03_

bool	isRelativeIso_

bool	isValidHltConfig_

bool	isW_

bool	isZGlbSta_

bool	isZGlbTrk_

bool	isZGolden1HLT_

bool	isZGolden2HLT_

bool	isZGoldenNoIso_

MonitorElement *	mass1HLT_

MonitorElement *	mass2HLT_

MonitorElement *	massGlbSta_

MonitorElement *	massGlbTrk_

MonitorElement *	massIsBothGlbTrkThanW_

MonitorElement *	massNotIso_

double	maxDeltaR_

double	maxDPtRel_

bool	metIncludesMuons_

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

double	mtMax_

double	mtMin_

edm::EDGetTokenT< edm::View< reco::Muon > >	muonToken_

unsigned int	n1hlt

unsigned int	n2hlt

unsigned int	nall

unsigned int	nEvWithHighPtMu

unsigned int	nGlbSta

unsigned int	nGlbTrk

unsigned int	nhlt

unsigned int	nInKinRange

unsigned int	niso

unsigned int	nNotIso

unsigned int	nsel

unsigned int	nTMass

unsigned int	nW

double	ptMuCut_

double	ptThreshold_

MonitorElement *	TMass_

edm::EDGetTokenT< reco::TrackCollection >	trackToken_

edm::EDGetTokenT< trigger::TriggerEvent >	trigEvToken_

edm::InputTag	trigTag_

edm::EDGetTokenT< edm::TriggerResults >	trigToken_

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 Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

Detailed Description

DQM offline for EWK MuLumiMonitor: intended for luminosity purposes using Z/W

Authors: : Michele de Gruttola, INFN Naples - Maria Cepeda, CIEMAT on behalf EWK-Muon group

Definition at line 36 of file EwkMuLumiMonitorDQM.h.

Constructor & Destructor Documentation

◆ EwkMuLumiMonitorDQM()

EwkMuLumiMonitorDQM::EwkMuLumiMonitorDQM ( const edm::ParameterSet & cfg )

Definition at line 35 of file EwkMuLumiMonitorDQM.cc.

References isValidHltConfig_.

     :  // Input collections
       trigTag_(cfg.getUntrackedParameter<edm::InputTag>("TrigTag", edm::InputTag("TriggerResults::HLT"))),
       trigToken_(consumes<edm::TriggerResults>(trigTag_)),
       trigEvToken_(consumes<trigger::TriggerEvent>(cfg.getUntrackedParameter<edm::InputTag>("triggerEvent"))),
       beamSpotToken_(consumes<reco::BeamSpot>(
           cfg.getUntrackedParameter<edm::InputTag>("offlineBeamSpot", edm::InputTag("offlineBeamSpot")))),
       muonToken_(consumes<edm::View<reco::Muon> >(cfg.getUntrackedParameter<edm::InputTag>("muons"))),
       trackToken_(consumes<reco::TrackCollection>(cfg.getUntrackedParameter<edm::InputTag>("tracks"))),
       caloTowerToken_(consumes<CaloTowerCollection>(cfg.getUntrackedParameter<edm::InputTag>("calotower"))),
       metToken_(consumes<edm::View<reco::MET> >(cfg.getUntrackedParameter<edm::InputTag>("metTag"))),
       metIncludesMuons_(cfg.getUntrackedParameter<bool>("METIncludesMuons")),
       // Main cuts
       // massMin_(cfg.getUntrackedParameter<double>("MtMin", 20.)),
       //   massMax_(cfg.getUntrackedParameter<double>("MtMax", 2000.))
       //  hltPath_(cfg.getUntrackedParameter<std::string> ("hltPath")) ,
       //  L3FilterName_(cfg.getUntrackedParameter<std::string>
       // ("L3FilterName")),
       ptMuCut_(cfg.getUntrackedParameter<double>("ptMuCut")),
       etaMuCut_(cfg.getUntrackedParameter<double>("etaMuCut")),
       isRelativeIso_(cfg.getUntrackedParameter<bool>("IsRelativeIso")),
       isCombinedIso_(cfg.getUntrackedParameter<bool>("IsCombinedIso")),
       isoCut03_(cfg.getUntrackedParameter<double>("IsoCut03")),
       //  deltaRTrk_(cfg.getUntrackedParameter<double>("deltaRTrk")),
       ptThreshold_(cfg.getUntrackedParameter<double>("ptThreshold")),
       // deltaRVetoTrk_(cfg.getUntrackedParameter<double>("deltaRVetoTrk")),
       maxDPtRel_(cfg.getUntrackedParameter<double>("maxDPtRel")),
       maxDeltaR_(cfg.getUntrackedParameter<double>("maxDeltaR")),
       mtMin_(cfg.getUntrackedParameter<double>("mtMin")),
       mtMax_(cfg.getUntrackedParameter<double>("mtMax")),
       acopCut_(cfg.getUntrackedParameter<double>("acopCut")),
       dxyCut_(cfg.getUntrackedParameter<double>("DxyCut")) {
   // just to initialize
   isValidHltConfig_ = false;
 }

Member Function Documentation

◆ analyze()

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

overridevirtual

number of prescale sets available

Reimplemented from DQMEDAnalyzer.

Definition at line 185 of file EwkMuLumiMonitorDQM.cc.

References acopCut_, beamSpotToken_, ewkMuLumiMonitorDQM_cfi::calotower, caloTowerToken_, reco::LeafCandidate::charge(), reco::TrackBase::charge(), PVValHelper::dxy, reco::TrackBase::dxy(), dxyCut_, edm::InputTag::encode(), PVValHelper::eta, reco::TrackBase::eta(), etaMuCut_, makeMEIFBenchmarkPlots::ev, f, dqm::impl::MonitorElement::Fill(), trigger::TriggerEvent::filterKeys(), trigger::TriggerEvent::filterTag(), trigger::TriggerEvent::getObjects(), highMass1HLT_, highMass2HLT_, highMassGlbSta_, highMassGlbTrk_, highMassIsBothGlbTrkThanW_, highMassNotIso_, hltConfigProvider_, mps_fire::i, IsMuMatchedToHLTMu(), isoCut03_, isW_, isZGlbSta_, isZGlbTrk_, isZGolden1HLT_, isZGolden2HLT_, isZGoldenNoIso_, dqmiolumiharvest::j, dqmdumpme::k, LogTrace, M_PI, EgHLTOffHistBins_cfi::mass, mass1HLT_, mass2HLT_, massGlbSta_, massGlbTrk_, massIsBothGlbTrkThanW_, massNotIso_, maxDeltaR_, maxDPtRel_, BTaggingMonitor_cfi::met, susyDQM_cfi::metCollection, metIncludesMuons_, metToken_, HLTConfigProvider::moduleLabels(), mtMax_, amptDefaultParameters_cff::mu, muIso(), PDWG_BPHSkim_cff::muons, muonToken_, dqmiodumpmetadata::n, n1hlt, n2hlt, nall, Skims_PA_cff::name, nEvWithHighPtMu, nGlbSta, nGlbTrk, nhlt, niso, nNotIso, nsel, nTMass, nW, MillePedeFileConverter_cfg::out, AlCaHLTBitMon_ParallelJobs::p, reco::TrackBase::p(), reco::LeafCandidate::p4(), reco::LeafCandidate::phi(), reco::BeamSpot::position(), HLTConfigProvider::prescaleSize(), HLTConfigProvider::prescaleValue(), DiDispStaMuonMonitor_cfi::pt, reco::LeafCandidate::pt(), reco::TrackBase::pt(), HLT_2022v14_cff::pt2, ptMuCut_, reco::LeafCandidate::px(), reco::TrackBase::px(), reco::LeafCandidate::py(), reco::TrackBase::py(), reco::TrackBase::pz(), alignCSCRings::s, findQualityFiles::size, trigger::TriggerEvent::sizeFilters(), mathSSE::sqrt(), contentValuesCheck::ss, AlCaHLTBitMon_QueryRunRegistry::string, tkIso(), TMass_, tracks, trackToken_, ClusterTask_cfi::trig, trigEvToken_, HLTConfigProvider::triggerIndex(), L1TEGammaOffline_cfi::triggerNames, HLTConfigProvider::triggerNames(), triggerResults, trigToken_, and findQualityFiles::v.

                                                                     {
   nall++;
   bool hlt_sel = false;
   double iso1 = -1;
   double iso2 = -1;
   bool isMu1Iso = false;
   bool isMu2Iso = false;
   bool singleTrigFlag1 = false;
   bool singleTrigFlag2 = false;
   isZGolden1HLT_ = false;
   isZGolden2HLT_ = false;
   isZGoldenNoIso_ = false;
   isZGlbSta_ = false;
   isZGlbTrk_ = false;
   isW_ = false;
   // Trigger
   bool trigger_fired = false;
 
   Handle<TriggerResults> triggerResults;
   if (!ev.getByToken(trigToken_, triggerResults)) {
     // LogWarning("") << ">>> TRIGGER collection does not exist !!!";
     return;
   }
 
   ev.getByToken(trigToken_, triggerResults);
   /*
     const edm::TriggerNames & trigNames = ev.triggerNames(*triggerResults);
 
 
   for (size_t i=0; i<triggerResults->size(); i++) {
   std::string trigName = trigNames.triggerName(i);
   //std::cout << " trigName == " << trigName << std::endl;
     if ( trigName == hltPath_ && triggerResults->accept(i)) {
     trigger_fired = true;
     hlt_sel=true;
     nhlt++;
     }
     }
   */
   // see the trigger single muon which are present
   string lowestMuonUnprescaledTrig = "";
   bool lowestMuonUnprescaledTrigFound = false;
   const std::vector<std::string>& triggerNames = hltConfigProvider_.triggerNames();
   for (size_t ts = 0; ts < triggerNames.size(); ts++) {
     string trig = triggerNames[ts];
     size_t f = trig.find("HLT_Mu");
     if ((f != std::string::npos)) {
       // std::cout << "single muon trigger present: " << trig << std::endl;
       // See if the trigger is prescaled;
       bool prescaled = false;
       const unsigned int prescaleSize = hltConfigProvider_.prescaleSize();
       for (unsigned int ps = 0; ps < prescaleSize; ps++) {
         if (hltConfigProvider_.prescaleValue<double>(ps, trig) != 1)
           prescaled = true;
       }
       if (!prescaled) {
         // looking now for the lowest hlt path not prescaled, with name of the
         // form HLT_MuX or HLTMuX_vY
         for (unsigned int n = 9; n < 100; n++) {
           string lowestTrig = "HLT_Mu";
           string lowestTrigv0 = "copy";
           std::stringstream out;
           out << n;
           std::string s = out.str();
           lowestTrig.append(s);
           lowestTrigv0 = lowestTrig;
           for (unsigned int v = 1; v < 10; v++) {
             lowestTrig.append("_v");
             std::stringstream oout;
             oout << v;
             std::string ss = oout.str();
             lowestTrig.append(ss);
             if (trig == lowestTrig)
               lowestMuonUnprescaledTrig = trig;
             if (trig == lowestTrig)
               lowestMuonUnprescaledTrigFound = true;
             if (trig == lowestTrig)
               break;
           }
           if (lowestMuonUnprescaledTrigFound)
             break;
 
           lowestTrig = lowestTrigv0;
           if (trig == lowestTrig)
             lowestMuonUnprescaledTrig = trig;
           //      if (trig==lowestTrig) {std::cout << " before break, lowestTrig
           // lowest single muon trigger present unprescaled: " << lowestTrig <<
           // std::endl; }
           if (trig == lowestTrig)
             lowestMuonUnprescaledTrigFound = true;
           if (trig == lowestTrig)
             break;
         }
         if (lowestMuonUnprescaledTrigFound)
           break;
       }
     }
   }
   //  std::cout << "after break, lowest single muon trigger present unprescaled:
   // " << lowestMuonUnprescaledTrig << std::endl;
   unsigned int triggerIndex;  // index of trigger path
 
   // See if event passed trigger paths
   std::string hltPath_ = lowestMuonUnprescaledTrig;
 
   triggerIndex = hltConfigProvider_.triggerIndex(hltPath_);
   if (triggerIndex < triggerResults->size())
     trigger_fired = triggerResults->accept(triggerIndex);
   std::string L3FilterName_ = "";
   if (trigger_fired) {
     const std::vector<std::string>& moduleLabs = hltConfigProvider_.moduleLabels(hltPath_);
     /*for (size_t k =0; k < moduleLabs.size()-1 ; k++){
       std::cout << "moduleLabs[" << k << "] == " << moduleLabs[k] << std::endl;
     }
     */
     // the l3 filter name is just the last module....
     size_t moduleLabsSizeMinus2 = moduleLabs.size() - 2;
     //  std::cout<<"moduleLabs[" << moduleLabsSizeMinus2 << "]== "<<
     // moduleLabs[moduleLabsSizeMinus2] << std::endl;
 
     L3FilterName_ = moduleLabs[moduleLabsSizeMinus2];
   }
 
   edm::Handle<trigger::TriggerEvent> handleTriggerEvent;
   LogTrace("") << ">>> Trigger bit: " << trigger_fired << " (" << hltPath_ << ")";
   if (!ev.getByToken(trigEvToken_, handleTriggerEvent)) {
     // LogWarning( "errorTriggerEventValid" ) << "trigger::TriggerEvent product
     // with InputTag " << trigEv_.encode() << " not in event";
     return;
   }
   ev.getByToken(trigEvToken_, handleTriggerEvent);
   const trigger::TriggerObjectCollection& toc(handleTriggerEvent->getObjects());
   size_t nMuHLT = 0;
   std::vector<reco::Particle> HLTMuMatched;
   for (size_t ia = 0; ia < handleTriggerEvent->sizeFilters(); ++ia) {
     std::string fullname = handleTriggerEvent->filterTag(ia).encode();
     // std::cout<< "fullname::== " << fullname<< std::endl;
     std::string name;
     size_t p = fullname.find_first_of(':');
     if (p != std::string::npos) {
       name = fullname.substr(0, p);
     } else {
       name = fullname;
     }
     if (!toc.empty()) {
       const trigger::Keys& k = handleTriggerEvent->filterKeys(ia);
       for (trigger::Keys::const_iterator ki = k.begin(); ki != k.end(); ++ki) {
         // looking at all the single muon l3 trigger present, for example
         // hltSingleMu15L3Filtered15.....
         if (name == L3FilterName_) {
           // trigger_fired = true;
           hlt_sel = true;
           nhlt++;
           HLTMuMatched.push_back(toc[*ki].particle());
           nMuHLT++;
         }
       }
     }
   }
 
   // Beam spot
   Handle<reco::BeamSpot> beamSpotHandle;
   if (!ev.getByToken(beamSpotToken_, beamSpotHandle)) {
     // LogWarning("") << ">>> No beam spot found !!!";
     return;
   }
 
   //  looping on muon....
   Handle<View<Muon> > muons;
   if (!ev.getByToken(muonToken_, muons)) {
     // LogError("") << ">>> muon collection does not exist !!!";
     return;
   }
 
   ev.getByToken(muonToken_, muons);
   // saving only muons with pt> ptMuCut and eta<etaMuCut, and dxy<dxyCut
   std::vector<reco::Muon> highPtGlbMuons;
   std::vector<reco::Muon> highPtStaMuons;
 
   for (size_t i = 0; i < muons->size(); i++) {
     const reco::Muon& mu = muons->at(i);
     double pt = mu.pt();
     double eta = mu.eta();
     if (pt > ptMuCut_ && fabs(eta) < etaMuCut_) {
       if (mu.isGlobalMuon()) {
         // check the dxy....
         double dxy = mu.innerTrack()->dxy(beamSpotHandle->position());
         if (fabs(dxy) > dxyCut_)
           continue;
         highPtGlbMuons.push_back(mu);
       }
       if (mu.isGlobalMuon())
         continue;
       // if is not, look if it is a standalone....
       if (mu.isStandAloneMuon())
         highPtStaMuons.push_back(mu);
       nEvWithHighPtMu++;
     }
   }
   size_t nHighPtGlbMu = highPtGlbMuons.size();
   size_t nHighPtStaMu = highPtStaMuons.size();
   if (hlt_sel && (nHighPtGlbMu > 0)) {
     // loop on high pt muons if there's at least two with opposite charge build
     // a Z, more then one z candidate is foreseen.........
     // stop the loop after 10 cicles....
     (nHighPtGlbMu > 10) ? nHighPtGlbMu = 10 : 1;
     // Z selection
     if (nHighPtGlbMu > 1) {
       for (unsigned int i = 0; i < nHighPtGlbMu; i++) {
         reco::Muon muon1 = highPtGlbMuons[i];
         const math::XYZTLorentzVector& mu1(muon1.p4());
         // double pt1= muon1.pt();
         for (unsigned int j = i + 1; j < nHighPtGlbMu; ++j) {
           reco::Muon muon2 = highPtGlbMuons[j];
           const math::XYZTLorentzVector& mu2(muon2.p4());
           // double pt2= muon2.pt();
           if (muon1.charge() == muon2.charge())
             continue;
           math::XYZTLorentzVector pair = mu1 + mu2;
           double mass = pair.M();
           // checking isolation and hlt maching
           iso1 = muIso(muon1);
           iso2 = muIso(muon2);
           isMu1Iso = (iso1 < isoCut03_);
           isMu2Iso = (iso2 < isoCut03_);
           singleTrigFlag1 = IsMuMatchedToHLTMu(muon1, HLTMuMatched, maxDeltaR_, maxDPtRel_);
           singleTrigFlag2 = IsMuMatchedToHLTMu(muon2, HLTMuMatched, maxDeltaR_, maxDPtRel_);
           if (singleTrigFlag1 && singleTrigFlag2)
             isZGolden2HLT_ = true;
           if ((singleTrigFlag1 && !singleTrigFlag2) || (!singleTrigFlag1 && singleTrigFlag2))
             isZGolden1HLT_ = true;
           // Z Golden passing all criteria, with 2 or 1 muon matched to an HLT
           // muon. Using the two cathegories as a control sample for the HLT
           // matching efficiency
           if (isZGolden2HLT_ || isZGolden1HLT_) {
             // a Z golden has been found, let's remove the two muons from the
             // list, dome for avoiding resolution effect enter muons in the
             // standalone and tracker collections.........
             highPtGlbMuons.erase(highPtGlbMuons.begin() + i);
             highPtGlbMuons.erase(highPtGlbMuons.begin() + j);
             // and updating the number of high pt muons....
             nHighPtGlbMu = highPtGlbMuons.size();
 
             if (isMu1Iso && isMu2Iso) {
               niso++;
               if (isZGolden2HLT_) {
                 n2hlt++;
                 mass2HLT_->Fill(mass);
                 highMass2HLT_->Fill(mass);
                 /*
                   if (pt1 > pt2) {
                   highest_mupt2HLT_ -> Fill (pt1);
                   lowest_mupt2HLT_ -> Fill (pt2);
                   } else {
                   highest_mupt2HLT_ -> Fill (pt2);
                   lowest_mupt2HLT_ -> Fill (pt1);
                   }
                 */
               }
               if (isZGolden1HLT_) {
                 n1hlt++;
                 mass1HLT_->Fill(mass);
                 highMass1HLT_->Fill(mass);
                 /*
                   if (pt1 >pt2) {
                   highest_mupt1HLT_ -> Fill (pt1);
                   lowest_mupt1HLT_ -> Fill (pt2);
                 } else {
                   highest_mupt1HLT_ -> Fill (pt2);
                   lowest_mupt1HLT_ -> Fill (pt1);
                 }
                 */
               }
             } else {
               // ZGlbGlb when at least one of the two muon is not isolated and
               // at least one HLT matched, used as control sample for the
               // isolation efficiency
               // filling events with one muon not isolated and both hlt mathced
               if (((isMu1Iso && !isMu2Iso) || (!isMu1Iso && isMu2Iso)) && (isZGolden2HLT_ && isZGolden1HLT_)) {
                 isZGoldenNoIso_ = true;
                 nNotIso++;
                 massNotIso_->Fill(mass);
                 highMassNotIso_->Fill(mass);
               }
               /*
                 if (pt1 > pt2) {
                 highest_muptNotIso_ -> Fill (pt1);
                 lowest_muptNotIso_ -> Fill (pt2);
                 } else {
                 highest_muptNotIso_ -> Fill (pt2);
                 lowest_muptNotIso_ -> Fill (pt1);
                 }
               */
             }
           }
         }
       }
     }
     // W->MuNu selection (if at least one muon has been selected)
     // looking for a W if a Z is not found.... let's think if we prefer to
     // exclude zMuMuNotIso or zMuSta....
     if (!(isZGolden2HLT_ || isZGolden1HLT_)) {
       Handle<View<MET> > metCollection;
       if (!ev.getByToken(metToken_, metCollection)) {
         // LogError("") << ">>> MET collection does not exist !!!";
         return;
       }
       const MET& met = metCollection->at(0);
       nW = 0;
       for (unsigned int i = 0; i < nHighPtGlbMu; i++) {
         reco::Muon muon1 = highPtGlbMuons[i];
         /*
                quality cut not implemented
             Quality Cuts           double dxy =
 gm->dxy(beamSpotHandle->position());
             Cut in 0.2           double trackerHits =
 gm->hitPattern().numberOfValidTrackerHits();
             Cut in 11           bool quality = fabs(dxy)<dxyCut_  &&
 muon::isGoodMuon(mu,muon::GlobalMuonPromptTight) && trackerHits>=trackerHitsCut_
 && mu.isTrackerMuon() ;
 if (!quality) continue;
         */
         // isolation cut and hlt maching
         iso1 = muIso(muon1);
         isMu1Iso = (iso1 < isoCut03_);
         if (!isMu1Iso)
           continue;
         // checking if muon is matched to any HLT muon....
         singleTrigFlag1 = IsMuMatchedToHLTMu(muon1, HLTMuMatched, maxDeltaR_, maxDPtRel_);
         if (!singleTrigFlag1)
           continue;
         // std::cout << " is GlobMu macthecd to HLT: " << singleTrigFlag1 <<
         // std::endl;
         // MT cuts
         double met_px = met.px();
         double met_py = met.py();
 
         if (!metIncludesMuons_) {
           for (unsigned int i = 0; i < nHighPtGlbMu; i++) {
             reco::Muon muon1 = highPtGlbMuons[i];
             met_px -= muon1.px();
             met_py -= muon1.py();
           }
         }
         double met_et = met.pt();  // sqrt(met_px*met_px+met_py*met_py);
         LogTrace("") << ">>> MET, MET_px, MET_py: " << met_et << ", " << met_px << ", " << met_py << " [GeV]";
         double w_et = met_et + muon1.pt();
         double w_px = met_px + muon1.px();
         double w_py = met_py + muon1.py();
         double massT = w_et * w_et - w_px * w_px - w_py * w_py;
         massT = (massT > 0) ? sqrt(massT) : 0;
         // Acoplanarity cuts
         Geom::Phi<double> deltaphi(muon1.phi() - atan2(met_py, met_px));
         double acop = M_PI - fabs(deltaphi.value());
         // std::cout << " is acp of W candidate less then cut: " << (acop<
         // acopCut_) << std::endl;
         if (acop > acopCut_)
           continue;  // Cut in 2.0
         TMass_->Fill(massT);
         if (massT < mtMin_ || massT > mtMax_)
           continue;  // Cut in (50,200) GeV
         // we give the number of W only in the tMass selected but we have a look
         // at mass tails to check the QCD background
         isW_ = true;
         nW++;
         nTMass++;
       }
     }
     // if a zGlobal is not selected, look at the dimuonGlobalOneStandAlone and
     // dimuonGlobalOneTrack...., used as a control sample for the tracking
     // efficiency  and muon system efficiency
     if (!(isZGolden2HLT_ || isZGolden1HLT_ || isZGoldenNoIso_)) {
       for (unsigned int i = 0; i < nHighPtGlbMu; ++i) {
         reco::Muon glbMuon = highPtGlbMuons[i];
         const math::XYZTLorentzVector& mu1(glbMuon.p4());
         // double pt1= glbMuon.pt();
         // checking that the global muon is hlt matched otherwise skip the event
         singleTrigFlag1 = IsMuMatchedToHLTMu(glbMuon, HLTMuMatched, maxDeltaR_, maxDPtRel_);
         if (!singleTrigFlag1)
           continue;
         // checking that the global muon is isolated matched otherwise skip the
         // event
         iso1 = muIso(glbMuon);
         isMu1Iso = (iso1 < isoCut03_);
         if (!isMu1Iso)
           continue;
         // look at the standalone muon ....
         // stop the loop after 10 cicles....
         (nHighPtStaMu > 10) ? nHighPtStaMu = 10 : 1;
         for (unsigned int j = 0; j < nHighPtStaMu; ++j) {
           reco::Muon staMuon = highPtStaMuons[j];
           const math::XYZTLorentzVector& mu2(staMuon.p4());
           // double pt2= staMuon.pt();
           if (glbMuon.charge() == staMuon.charge())
             continue;
           math::XYZTLorentzVector pair = mu1 + mu2;
           double mass = pair.M();
           iso2 = muIso(staMuon);
           isMu2Iso = (iso2 < isoCut03_);
           LogTrace("") << "\t... isolation value" << iso1 << ", isolated? " << isMu1Iso;
           LogTrace("") << "\t... isolation value" << iso2 << ", isolated? " << isMu2Iso;
           // requiring theat the global mu is mathed to the HLT  and both are
           // isolated
           if (isMu2Iso) {
             isZGlbSta_ = true;
             nGlbSta++;
             massGlbSta_->Fill(mass);
             highMassGlbSta_->Fill(mass);
             /*
               if (pt1 > pt2) {
               highest_muptGlbSta_ -> Fill (pt1);
               lowest_muptGlbSta_ -> Fill (pt2);
             } else {
               highest_muptGlbSta_ -> Fill (pt2);
               lowest_muptGlbSta_ -> Fill (pt1);
             }
             */
           }
         }
         // look at the tracks....
         Handle<TrackCollection> tracks;
         if (!ev.getByToken(trackToken_, tracks)) {
           // LogError("") << ">>> track collection does not exist !!!";
           return;
         }
         ev.getByToken(trackToken_, tracks);
         Handle<CaloTowerCollection> calotower;
         if (!ev.getByToken(caloTowerToken_, calotower)) {
           // LogError("") << ">>> calotower collection does not exist !!!";
           return;
         }
         ev.getByToken(caloTowerToken_, calotower);
         // avoid to loop on more than 5000 trks
         size_t nTrk = tracks->size();
         (nTrk > 5000) ? nTrk = 5000 : 1;
         for (unsigned int j = 0; j < nTrk; j++) {
           const reco::Track& tk = tracks->at(j);
           if (glbMuon.charge() == tk.charge())
             continue;
           double pt2 = tk.pt();
           double eta = tk.eta();
           double dxy = tk.dxy(beamSpotHandle->position());
           if (pt2 < ptMuCut_ || fabs(eta) > etaMuCut_ || fabs(dxy) > dxyCut_)
             continue;
           // assuming that the track is a mu....
           math::XYZTLorentzVector mu2(tk.px(), tk.py(), tk.pz(), sqrt((tk.p() * tk.p()) + (0.10566 * 0.10566)));
           math::XYZTLorentzVector pair = mu1 + mu2;
           double mass = pair.M();
           // now requiring that the tracks is isolated.......
           iso2 = tkIso(tk, tracks, calotower);
           isMu2Iso = (iso2 < isoCut03_);
           //          std::cout << "found a track with rel/comb iso: " << iso2
           //<< std::endl;
           if (isMu2Iso) {
             isZGlbTrk_ = true;
             nGlbTrk++;
             massGlbTrk_->Fill(mass);
             highMassGlbTrk_->Fill(mass);
             if (!isW_)
               continue;
             massIsBothGlbTrkThanW_->Fill(mass);
             highMassIsBothGlbTrkThanW_->Fill(mass);
             /*
               if (pt1 > pt2) {
               highest_muptGlbTrk_ -> Fill (pt1);
               lowest_muptGlbTrk_ -> Fill (pt2);
               } else {
               highest_muptGlbTrk_ -> Fill (pt2);
               lowest_muptGlbTrk_ -> Fill (pt1);
             }
             */
           }
         }
       }
     }
   }
   if ((hlt_sel || isZGolden2HLT_ || isZGolden1HLT_ || isZGoldenNoIso_ || isZGlbSta_ || isZGlbTrk_ || isW_)) {
     nsel++;
     LogTrace("") << ">>>> Event ACCEPTED";
   } else {
     LogTrace("") << ">>>> Event REJECTED";
   }
   return;
 }

◆ bookHistograms()

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 95 of file EwkMuLumiMonitorDQM.cc.

References dqm::implementation::IBooker::book1D(), highMass1HLT_, highMass2HLT_, highMassGlbSta_, highMassGlbTrk_, highMassIsBothGlbTrkThanW_, highMassNotIso_, mass1HLT_, mass2HLT_, massGlbSta_, massGlbTrk_, massIsBothGlbTrkThanW_, massNotIso_, dqm::implementation::NavigatorBase::setCurrentFolder(), and TMass_.

                                                                                                       {
   ibooker.setCurrentFolder("Physics/EwkMuLumiMonitorDQM");
 
   mass2HLT_ = ibooker.book1D("Z_2HLT_MASS", "Z mass [GeV/c^{2}]", 200, 0., 200.);
   mass1HLT_ = ibooker.book1D("Z_1HLT_MASS", "Z mass [GeV/c^{2}]", 200, 0., 200.);
   massNotIso_ = ibooker.book1D("Z_NOTISO_MASS", "Z mass [GeV/c^{2}]", 200, 0., 200.);
   massGlbSta_ = ibooker.book1D("Z_GLBSTA_MASS", "Z mass [GeV/c^{2}]", 200, 0., 200.);
   massGlbTrk_ = ibooker.book1D("Z_GLBTRK_MASS", "Z mass [GeV/c^{2}]", 200, 0., 200.);
   massIsBothGlbTrkThanW_ = ibooker.book1D("Z_ISBOTHGLBTRKTHANW_MASS", "Z mass [GeV/c^{2}]", 200, 0., 200.);
 
   highMass2HLT_ = ibooker.book1D("Z_2HLT_HIGHMASS", "Z high mass [GeV/c^{2}]", 2000, 0., 2000.);
   highMass1HLT_ = ibooker.book1D("Z_1HLT_HIGHMASS", "Z high mass [GeV/c^{2}]", 2000, 0., 2000.);
   highMassNotIso_ = ibooker.book1D("Z_NOTISO_HIGHMASS", "Z high mass [GeV/c^{2}]", 2000, 0., 2000.);
   highMassGlbSta_ = ibooker.book1D("Z_GLBSTA_HIGHMASS", "Z high mass [GeV/c^{2}]", 2000, 0., 2000.);
   highMassGlbTrk_ = ibooker.book1D("Z_GLBTRK_HIGHMASS", "Z high mass [GeV/c^{2}]", 2000, 0., 2000.);
   highMassIsBothGlbTrkThanW_ =
       ibooker.book1D("Z_ISBOTHGLBTRKTHANW_HIGHMASS", "Z high mass [GeV/c^{2}]", 2000, 0., 2000.);
 
   TMass_ = ibooker.book1D("TMASS", "Transverse mass [GeV]", 300, 0., 300.);
 }

◆ dqmBeginRun()

void EwkMuLumiMonitorDQM::dqmBeginRun	(	const edm::Run &	r,
		const edm::EventSetup &	iSetup
	)

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 71 of file EwkMuLumiMonitorDQM.cc.

References hltConfigProvider_, HLTConfigProvider::init(), isValidHltConfig_, n1hlt, n2hlt, nall, nEvWithHighPtMu, nGlbSta, nGlbTrk, nhlt, nInKinRange, niso, nNotIso, nsel, nTMass, nW, edm::InputTag::process(), alignCSCRings::r, and trigTag_.

                                                                             {
   nall = 0;
   nEvWithHighPtMu = 0;
   nInKinRange = 0;
   nsel = 0;
   niso = 0;
   nhlt = 0;
   n1hlt = 0;
   n2hlt = 0;
   nNotIso = 0;
   nGlbSta = 0;
   nGlbTrk = 0;
   nTMass = 0;
   nW = 0;
 
   // passed as parameter to HLTConfigProvider::init(), not yet used
   bool isConfigChanged = false;
 
   // isValidHltConfig_ used to short-circuit analyze() in case of problems
   isValidHltConfig_ = hltConfigProvider_.init(r, iSetup, trigTag_.process(), isConfigChanged);
   // std::cout << "hlt config trigger is valid??" << isValidHltConfig_ <<
   // std::endl;
 }

◆ init_histograms()

void EwkMuLumiMonitorDQM::init_histograms ( )

◆ IsMuMatchedToHLTMu()

bool EwkMuLumiMonitorDQM::IsMuMatchedToHLTMu	(	const reco::Muon &	mu,
		const std::vector< reco::Particle > &	HLTMu,
		double	DR,
		double	DPtRel
	)

Definition at line 169 of file EwkMuLumiMonitorDQM.cc.

References PbPb_ZMuSkimMuonDPG_cff::deltaR, dqmdumpme::k, amptDefaultParameters_cff::mu, nJettinessAdder_cfi::nPass, and DiDispStaMuonMonitor_cfi::pt.

Referenced by analyze().

                                                             {
   size_t dim = HLTMu.size();
   size_t nPass = 0;
   if (dim == 0)
     return false;
   for (size_t k = 0; k < dim; k++) {
     if ((deltaR(HLTMu[k], mu) < DR) && (fabs(HLTMu[k].pt() - mu.pt()) / HLTMu[k].pt() < DPtRel)) {
       nPass++;
     }
   }
   return (nPass > 0);
 }

◆ muIso()

double EwkMuLumiMonitorDQM::muIso ( const reco::Muon & mu )

Definition at line 116 of file EwkMuLumiMonitorDQM.cc.

References isCombinedIso_, isRelativeIso_, and amptDefaultParameters_cff::mu.

Referenced by analyze().

                                                     {
   double isovar = mu.isolationR03().sumPt;
   if (isCombinedIso_) {
     isovar += mu.isolationR03().emEt;
     isovar += mu.isolationR03().hadEt;
   }
   if (isRelativeIso_)
     isovar /= mu.pt();
   return isovar;
 }

◆ tkIso()

double EwkMuLumiMonitorDQM::tkIso	(	const reco::Track &	tk,
		edm::Handle< reco::TrackCollection >	tracks,
		edm::Handle< CaloTowerCollection >	calotower
	)

Definition at line 127 of file EwkMuLumiMonitorDQM.cc.

References ewkMuLumiMonitorDQM_cfi::calotower, PbPb_ZMuSkimMuonDPG_cff::deltaR, HGC3DClusterGenMatchSelector_cfi::dR, PVValHelper::dz, reco::TrackBase::eta(), mps_fire::i, isCombinedIso_, isRelativeIso_, reco::Track::outerEta(), reco::Track::outerPhi(), reco::TrackBase::phi(), reco::TrackBase::pt(), ptThreshold_, mathSSE::sqrt(), tracks, reco::TrackBase::vx(), reco::TrackBase::vy(), and reco::TrackBase::vz().

Referenced by analyze().

                                                                          {
   double ptSum = 0;
   for (size_t i = 0; i < tracks->size(); ++i) {
     const reco::Track& elem = tracks->at(i);
     double elemPt = elem.pt();
     // same parameter used for muIsolation: dR [0.01, IsoCut03_], |dZ|<0.2,
     // |d_r(xy)|<0.1
     double elemVx = elem.vx();
     double elemVy = elem.vy();
     double elemD0 = sqrt(elemVx * elemVx + elemVy * elemVy);
     if (elemD0 > 0.2)
       continue;
     double dz = fabs(elem.vz() - tk.vz());
     if (dz > 0.1)
       continue;
     // evaluate only for tracks with pt>ptTreshold
     if (elemPt < ptThreshold_)
       continue;
     double dR = deltaR(elem.eta(), elem.phi(), tk.eta(), tk.phi());
     // isolation in a cone with dR=0.3, and vetoing the track itself
     if ((dR < 0.01) || (dR > 0.3))
       continue;
     ptSum += elemPt;
   }
   if (isCombinedIso_) {
     // loop on clusters....
     for (CaloTowerCollection::const_iterator it = calotower->begin(); it != calotower->end(); ++it) {
       double dR = deltaR(it->eta(), it->phi(), tk.outerEta(), tk.outerPhi());
       // veto value is 0.1 for towers....
       if ((dR < 0.1) || (dR > 0.3))
         continue;
       ptSum += it->emEnergy();
       ptSum += it->hadEnergy();
     }
   }
   if (isRelativeIso_)
     ptSum /= tk.pt();
   return ptSum;
 }

Member Data Documentation

◆ acopCut_

double EwkMuLumiMonitorDQM::acopCut_

private

Definition at line 77 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().

◆ beamSpotToken_

edm::EDGetTokenT<reco::BeamSpot> EwkMuLumiMonitorDQM::beamSpotToken_

private

Definition at line 53 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().

◆ caloTowerToken_

edm::EDGetTokenT<CaloTowerCollection> EwkMuLumiMonitorDQM::caloTowerToken_

private

Definition at line 56 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().

◆ deltaRTrk_

double EwkMuLumiMonitorDQM::deltaRTrk_

private

Definition at line 70 of file EwkMuLumiMonitorDQM.h.

◆ deltaRVetoTrk_

double EwkMuLumiMonitorDQM::deltaRVetoTrk_

private

Definition at line 72 of file EwkMuLumiMonitorDQM.h.

◆ dxyCut_

double EwkMuLumiMonitorDQM::dxyCut_

private

Definition at line 78 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().

◆ etaMuCut_

double EwkMuLumiMonitorDQM::etaMuCut_

private

Definition at line 64 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().

◆ highMass1HLT_

MonitorElement* EwkMuLumiMonitorDQM::highMass1HLT_

private

Definition at line 86 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and bookHistograms().

◆ highMass2HLT_

MonitorElement* EwkMuLumiMonitorDQM::highMass2HLT_

private

Definition at line 81 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and bookHistograms().

◆ highMassGlbSta_

MonitorElement* EwkMuLumiMonitorDQM::highMassGlbSta_

private

Definition at line 96 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and bookHistograms().

◆ highMassGlbTrk_

MonitorElement* EwkMuLumiMonitorDQM::highMassGlbTrk_

private

Definition at line 101 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and bookHistograms().

◆ highMassIsBothGlbTrkThanW_

MonitorElement* EwkMuLumiMonitorDQM::highMassIsBothGlbTrkThanW_

private

Definition at line 108 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and bookHistograms().

◆ highMassNotIso_

MonitorElement* EwkMuLumiMonitorDQM::highMassNotIso_

private

Definition at line 91 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and bookHistograms().

◆ hltConfigProvider_

HLTConfigProvider EwkMuLumiMonitorDQM::hltConfigProvider_

private

Definition at line 132 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and dqmBeginRun().

◆ isCombinedIso_

bool EwkMuLumiMonitorDQM::isCombinedIso_

private

Definition at line 67 of file EwkMuLumiMonitorDQM.h.

Referenced by muIso(), and tkIso().

◆ isoCut03_

double EwkMuLumiMonitorDQM::isoCut03_

private

Definition at line 68 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().

◆ isRelativeIso_

bool EwkMuLumiMonitorDQM::isRelativeIso_

private

Definition at line 66 of file EwkMuLumiMonitorDQM.h.

Referenced by muIso(), and tkIso().

◆ isValidHltConfig_

bool EwkMuLumiMonitorDQM::isValidHltConfig_

private

Definition at line 131 of file EwkMuLumiMonitorDQM.h.

Referenced by dqmBeginRun(), and EwkMuLumiMonitorDQM().

◆ isW_

bool EwkMuLumiMonitorDQM::isW_

private

Definition at line 129 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().

◆ isZGlbSta_

bool EwkMuLumiMonitorDQM::isZGlbSta_

private

Definition at line 127 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().

◆ isZGlbTrk_

bool EwkMuLumiMonitorDQM::isZGlbTrk_

private

Definition at line 128 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().

◆ isZGolden1HLT_

bool EwkMuLumiMonitorDQM::isZGolden1HLT_

private

Definition at line 124 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().

◆ isZGolden2HLT_

bool EwkMuLumiMonitorDQM::isZGolden2HLT_

private

Definition at line 125 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().

◆ isZGoldenNoIso_

bool EwkMuLumiMonitorDQM::isZGoldenNoIso_

private

Definition at line 126 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().

◆ mass1HLT_

MonitorElement* EwkMuLumiMonitorDQM::mass1HLT_

private

Definition at line 85 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and bookHistograms().

◆ mass2HLT_

MonitorElement* EwkMuLumiMonitorDQM::mass2HLT_

private

Definition at line 80 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and bookHistograms().

◆ massGlbSta_

MonitorElement* EwkMuLumiMonitorDQM::massGlbSta_

private

Definition at line 95 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and bookHistograms().

◆ massGlbTrk_

MonitorElement* EwkMuLumiMonitorDQM::massGlbTrk_

private

Definition at line 100 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and bookHistograms().

◆ massIsBothGlbTrkThanW_

MonitorElement* EwkMuLumiMonitorDQM::massIsBothGlbTrkThanW_

private

Definition at line 107 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and bookHistograms().

◆ massNotIso_

MonitorElement* EwkMuLumiMonitorDQM::massNotIso_

private

Definition at line 90 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and bookHistograms().

◆ maxDeltaR_

double EwkMuLumiMonitorDQM::maxDeltaR_

private

Definition at line 74 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().

◆ maxDPtRel_

double EwkMuLumiMonitorDQM::maxDPtRel_

private

Definition at line 73 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().

◆ metIncludesMuons_

bool EwkMuLumiMonitorDQM::metIncludesMuons_

private

Definition at line 58 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().

◆ metToken_

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

private

Definition at line 57 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().

◆ mtMax_

double EwkMuLumiMonitorDQM::mtMax_

private

Definition at line 76 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().

◆ mtMin_

double EwkMuLumiMonitorDQM::mtMin_

private

Definition at line 75 of file EwkMuLumiMonitorDQM.h.

◆ muonToken_

edm::EDGetTokenT<edm::View<reco::Muon> > EwkMuLumiMonitorDQM::muonToken_

private

Definition at line 54 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().

◆ n1hlt

unsigned int EwkMuLumiMonitorDQM::n1hlt

private

Definition at line 116 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and dqmBeginRun().

◆ n2hlt

unsigned int EwkMuLumiMonitorDQM::n2hlt

private

Definition at line 117 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and dqmBeginRun().

◆ nall

unsigned int EwkMuLumiMonitorDQM::nall

private

Definition at line 110 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and dqmBeginRun().

◆ nEvWithHighPtMu

unsigned int EwkMuLumiMonitorDQM::nEvWithHighPtMu

private

Definition at line 111 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and dqmBeginRun().

◆ nGlbSta

unsigned int EwkMuLumiMonitorDQM::nGlbSta

private

Definition at line 119 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and dqmBeginRun().

◆ nGlbTrk

unsigned int EwkMuLumiMonitorDQM::nGlbTrk

private

Definition at line 120 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and dqmBeginRun().

◆ nhlt

unsigned int EwkMuLumiMonitorDQM::nhlt

private

Definition at line 115 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and dqmBeginRun().

◆ nInKinRange

unsigned int EwkMuLumiMonitorDQM::nInKinRange

private

Definition at line 112 of file EwkMuLumiMonitorDQM.h.

Referenced by dqmBeginRun().

◆ niso

unsigned int EwkMuLumiMonitorDQM::niso

private

Definition at line 114 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and dqmBeginRun().

◆ nNotIso

unsigned int EwkMuLumiMonitorDQM::nNotIso

private

Definition at line 118 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and dqmBeginRun().

◆ nsel

unsigned int EwkMuLumiMonitorDQM::nsel

private

Definition at line 113 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and dqmBeginRun().

◆ nTMass

unsigned int EwkMuLumiMonitorDQM::nTMass

private

Definition at line 121 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and dqmBeginRun().

◆ nW

unsigned int EwkMuLumiMonitorDQM::nW

private

Definition at line 122 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and dqmBeginRun().

◆ ptMuCut_

double EwkMuLumiMonitorDQM::ptMuCut_

private

Definition at line 63 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().

◆ ptThreshold_

double EwkMuLumiMonitorDQM::ptThreshold_

private

Definition at line 71 of file EwkMuLumiMonitorDQM.h.

Referenced by tkIso().

◆ TMass_

MonitorElement* EwkMuLumiMonitorDQM::TMass_

private

Definition at line 105 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze(), and bookHistograms().

◆ trackToken_

edm::EDGetTokenT<reco::TrackCollection> EwkMuLumiMonitorDQM::trackToken_

private

Definition at line 55 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().

◆ trigEvToken_

edm::EDGetTokenT<trigger::TriggerEvent> EwkMuLumiMonitorDQM::trigEvToken_

private

Definition at line 52 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().

◆ trigTag_

edm::InputTag EwkMuLumiMonitorDQM::trigTag_

private

Definition at line 50 of file EwkMuLumiMonitorDQM.h.

Referenced by dqmBeginRun().

◆ trigToken_

edm::EDGetTokenT<edm::TriggerResults> EwkMuLumiMonitorDQM::trigToken_

private

Definition at line 51 of file EwkMuLumiMonitorDQM.h.

Referenced by analyze().

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ EwkMuLumiMonitorDQM()

Member Function Documentation

◆ analyze()

◆ bookHistograms()

◆ dqmBeginRun()

◆ init_histograms()

◆ IsMuMatchedToHLTMu()

◆ muIso()

◆ tkIso()

Member Data Documentation

◆ acopCut_

◆ beamSpotToken_

◆ caloTowerToken_

◆ deltaRTrk_

◆ deltaRVetoTrk_

◆ dxyCut_

◆ etaMuCut_

◆ highMass1HLT_

◆ highMass2HLT_

◆ highMassGlbSta_

◆ highMassGlbTrk_

◆ highMassIsBothGlbTrkThanW_

◆ highMassNotIso_

◆ hltConfigProvider_

◆ isCombinedIso_

◆ isoCut03_

◆ isRelativeIso_

◆ isValidHltConfig_

◆ isW_

◆ isZGlbSta_

◆ isZGlbTrk_

◆ isZGolden1HLT_

◆ isZGolden2HLT_

◆ isZGoldenNoIso_

◆ mass1HLT_

◆ mass2HLT_

◆ massGlbSta_

◆ massGlbTrk_

◆ massIsBothGlbTrkThanW_

◆ massNotIso_

◆ maxDeltaR_

◆ maxDPtRel_

◆ metIncludesMuons_

◆ metToken_

◆ mtMax_

◆ mtMin_

◆ muonToken_

◆ n1hlt

◆ n2hlt

◆ nall

◆ nEvWithHighPtMu

◆ nGlbSta

◆ nGlbTrk

◆ nhlt

◆ nInKinRange

◆ niso

◆ nNotIso

◆ nsel

◆ nTMass

◆ nW

◆ ptMuCut_

◆ ptThreshold_

◆ TMass_

◆ trackToken_

◆ trigEvToken_

◆ trigTag_

◆ trigToken_