Inheritance diagram for MuTriggerAnalyzer:

Public Member Functions
	MuTriggerAnalyzer (const edm::ParameterSet &pset)

Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

	~EDAnalyzer () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase &	operator= (EDConsumerBase &&)=default

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

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

void	endJob () override

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

Private Attributes
TH1D *	deltaPtOverPt_

TH1D *	deltaR_

double	etaMuCut_

edm::Handle< trigger::TriggerEvent >	handleTriggerEvent_

std::string	hltPath_

TH1D *	hTrigMuonPtDenS_

TH1D *	hTrigMuonPtNumS_

std::string	L3FilterName_

double	maxDeltaR_

double	maxDPtRel_

edm::EDGetTokenT< View< Muon > >	muonToken_

const int	nbins_

const double	ptMax_

double	ptMuCut_

edm::InputTag	trigEvTag_

edm::EDGetTokenT< trigger::TriggerEvent >	trigEvToken_

edm::Handle< edm::TriggerResults >	triggerResults_

edm::TriggerNames const *	trigNames_

edm::EDGetTokenT< TriggerResults >	trigToken_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Detailed Description

Definition at line 47 of file MuTriggerAnalyzer.cc.

Constructor & Destructor Documentation

◆ MuTriggerAnalyzer()

MuTriggerAnalyzer::MuTriggerAnalyzer ( const edm::ParameterSet & pset )

Definition at line 110 of file MuTriggerAnalyzer.cc.

     : trigToken_(consumes<TriggerResults>(cfg.getParameter<edm::InputTag>("TrigTag"))),
       trigEvTag_(cfg.getParameter<edm::InputTag>("triggerEvent")),
       trigEvToken_(consumes<trigger::TriggerEvent>(trigEvTag_)),
       muonToken_(consumes<View<Muon> >(cfg.getUntrackedParameter<edm::InputTag>("muons"))),
       ptMuCut_(cfg.getUntrackedParameter<double>("ptMuCut")),
       etaMuCut_(cfg.getUntrackedParameter<double>("etaMuCut")),
       hltPath_(cfg.getParameter<std::string>("hltPath")),
       L3FilterName_(cfg.getParameter<std::string>("L3FilterName")),
       maxDPtRel_(cfg.getParameter<double>("maxDPtRel")),
       maxDeltaR_(cfg.getParameter<double>("maxDeltaR")),
       nbins_(cfg.getParameter<double>("ptMax_")),
       ptMax_(cfg.getParameter<double>("ptMax_")) {
   Service<TFileService> fs;
  
   hTrigMuonPtNumS_ = fs->make<TH1D>("hTrigMuonPtNumS", "hTrigMuonPtNumS", nbins_ + 1, 0, ptMax_);
   hTrigMuonPtDenS_ = fs->make<TH1D>("hTrigMuonPtDenS", "hTrigMuonPtDenS", nbins_ + 1, 0, ptMax_);
   deltaR_ = fs->make<TH1D>("deltaR", "deltaR", nbins_ + 1, 0, maxDeltaR_);
   deltaPtOverPt_ = fs->make<TH1D>("deltaPtOverPt", "deltaPtOverPt", nbins_ + 1, 0, maxDPtRel_);
 }

References deltaPtOverPt_, deltaR_, hTrigMuonPtDenS_, hTrigMuonPtNumS_, TFileService::make(), maxDeltaR_, maxDPtRel_, nbins_, and ptMax_.

Member Function Documentation

◆ analyze()

void MuTriggerAnalyzer::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	setup
	)

overrideprivatevirtual

Implements edm::EDAnalyzer.

Definition at line 138 of file MuTriggerAnalyzer.cc.

                                                                   {
   //      bool overlap = false;
  
   // Trigger
   Handle<TriggerResults> triggerResults;
   TriggerNames trigNames;
   if (!ev.getByToken(trigToken_, triggerResults)) {
     LogWarning("") << ">>> TRIGGER collection does not exist !!!";
     return;
   }
   // trigNames.init(*triggerResults);
   trigNames_ = &ev.triggerNames(*triggerResults);
   //bool trigger_fired = false;
   for (unsigned int i = 0; i < triggerResults->size(); i++) {
     std::string trigName = trigNames_->triggerName(i);
     if (trigName == hltPath_ && triggerResults->accept(i)) {
       //trigger_fired = true;
     }
   }
   edm::Handle<trigger::TriggerEvent> handleTriggerEvent;
   //  LogTrace("") << ">>> Trigger bit: " << trigger_fired << " (" << hltPath_ << ")";
   if (!ev.getByToken(trigEvToken_, handleTriggerEvent)) {
     LogWarning("errorTriggerEventValid") << "trigger::TriggerEvent product with InputTag " << trigEvTag_.encode()
                                          << " not in event";
     return;
   }
   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::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) {
         if (name == L3FilterName_) {
           HLTMuMatched.push_back(toc[*ki].particle());
           nMuHLT++;
         }
       }
     }
   }
  
   //  looping on muon....
   Handle<View<Muon> > muons;
   if (!ev.getByToken(muonToken_, muons)) {
     LogError("") << ">>> muon collection does not exist !!!";
     return;
   }
   //saving only muons with pt> ptMuCut and eta<etaMuCut
   std::vector<reco::Muon> highPtGlbMuons;
  
   for (unsigned int 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())
         highPtGlbMuons.push_back(mu);
     }
   }
   unsigned int nHighPtGlbMu = highPtGlbMuons.size();
   std::cout << "I've got " << nHighPtGlbMu << " nHighPtGlbMu" << std::endl;
   // unsigned int nHighPtStaMu = highPtStaMuons.size();
   // stop the loop after 10 cicles....
   if (nHighPtGlbMu > 10)
     nHighPtGlbMu = 10;
  
   for (unsigned int i = 0; i < nHighPtGlbMu; i++) {
     IsMuMatchedToHLTMu(highPtGlbMuons[i], HLTMuMatched, maxDeltaR_, maxDPtRel_);
   }
 }

References gather_cfg::cout, edm::InputTag::encode(), PVValHelper::eta, etaMuCut_, ev, trigger::TriggerEvent::filterKeys(), trigger::TriggerEvent::filterTag(), trigger::TriggerEvent::getObjects(), hltPath_, mps_fire::i, IsMuMatchedToHLTMu(), dqmdumpme::k, L3FilterName_, maxDeltaR_, maxDPtRel_, amptDefaultParameters_cff::mu, PDWG_BPHSkim_cff::muons, muonToken_, Skims_PA_cff::name, AlCaHLTBitMon_ParallelJobs::p, DiDispStaMuonMonitor_cfi::pt, ptMuCut_, trigger::TriggerEvent::sizeFilters(), AlCaHLTBitMon_QueryRunRegistry::string, trigEvTag_, trigEvToken_, edm::TriggerNames::triggerName(), triggerResults, EgHLTOffTrigSelection_cfi::trigName, trigNames, trigNames_, and trigToken_.

◆ endJob()

void MuTriggerAnalyzer::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 131 of file MuTriggerAnalyzer.cc.

                                {
   for (int i = 0; i < nbins_ + 1; ++i) {
     std::cout << "number of reco muon in bin " << i << " = " << hTrigMuonPtDenS_->GetBinContent(i) << std::endl;
     std::cout << "number of hlt muon in bin " << i << " = " << hTrigMuonPtNumS_->GetBinContent(i) << std::endl;
   }
 }

References gather_cfg::cout, hTrigMuonPtDenS_, hTrigMuonPtNumS_, mps_fire::i, and nbins_.

◆ IsMuMatchedToHLTMu()

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

private

Definition at line 76 of file MuTriggerAnalyzer.cc.

                                                           {
   size_t dim = HLTMu.size();
   size_t nPass = 0;
  
   // filling the denumerator;
   double muRecoPt = mu.pt();
   hTrigMuonPtDenS_->Fill(muRecoPt);
  
   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++;
       std::cout << " matching a muon " << std::endl;
       std::cout << "muon reco  pt : " << muRecoPt << std::endl;
       std::cout << "muon reco  eta " << mu.eta() << std::endl;
       std::cout << "muon trigger  pt " << HLTMu[k].pt() << std::endl;
       // filling the numerator, at the same bin as the denum.....
       hTrigMuonPtNumS_->Fill(muRecoPt);
       deltaR_->Fill(deltaR(HLTMu[k], mu));
       deltaPtOverPt_->Fill(fabs(HLTMu[k].pt() - mu.pt()) / HLTMu[k].pt());
  
       std::cout << "muon trigger  eta : " << HLTMu[k].eta() << std::endl;
       std::cout << "deltaR((HLTMu[k], mu)): " << deltaR(HLTMu[k], mu) << std::endl;
       std::cout << "deltaPtOverPt: " << fabs(HLTMu[k].pt() - mu.pt()) / HLTMu[k].pt() << std::endl;
     }
   }
  
   return (nPass > 0);
 }

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

Referenced by analyze().

Member Data Documentation

◆ deltaPtOverPt_

TH1D* MuTriggerAnalyzer::deltaPtOverPt_

private

Definition at line 73 of file MuTriggerAnalyzer.cc.

Referenced by MuTriggerAnalyzer().

◆ deltaR_

TH1D* MuTriggerAnalyzer::deltaR_

private

Definition at line 72 of file MuTriggerAnalyzer.cc.

Referenced by MuTriggerAnalyzer().

◆ etaMuCut_

double MuTriggerAnalyzer::etaMuCut_

private

Definition at line 60 of file MuTriggerAnalyzer.cc.

Referenced by analyze().

◆ handleTriggerEvent_

edm::Handle<trigger::TriggerEvent> MuTriggerAnalyzer::handleTriggerEvent_

private

Definition at line 65 of file MuTriggerAnalyzer.cc.

◆ hltPath_

std::string MuTriggerAnalyzer::hltPath_

private

Definition at line 61 of file MuTriggerAnalyzer.cc.

Referenced by analyze().

◆ hTrigMuonPtDenS_

TH1D* MuTriggerAnalyzer::hTrigMuonPtDenS_

private

Definition at line 71 of file MuTriggerAnalyzer.cc.

Referenced by endJob(), and MuTriggerAnalyzer().

◆ hTrigMuonPtNumS_

TH1D* MuTriggerAnalyzer::hTrigMuonPtNumS_

private

Definition at line 70 of file MuTriggerAnalyzer.cc.

Referenced by endJob(), and MuTriggerAnalyzer().

◆ L3FilterName_

std::string MuTriggerAnalyzer::L3FilterName_

private

Definition at line 62 of file MuTriggerAnalyzer.cc.

Referenced by analyze().

◆ maxDeltaR_

double MuTriggerAnalyzer::maxDeltaR_

private

Definition at line 67 of file MuTriggerAnalyzer.cc.

Referenced by analyze(), and MuTriggerAnalyzer().

◆ maxDPtRel_

double MuTriggerAnalyzer::maxDPtRel_

private

Definition at line 67 of file MuTriggerAnalyzer.cc.

Referenced by analyze(), and MuTriggerAnalyzer().

◆ muonToken_

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

private

Definition at line 59 of file MuTriggerAnalyzer.cc.

Referenced by analyze().

◆ nbins_

const int MuTriggerAnalyzer::nbins_

private

Definition at line 68 of file MuTriggerAnalyzer.cc.

Referenced by endJob(), and MuTriggerAnalyzer().

◆ ptMax_

const double MuTriggerAnalyzer::ptMax_

private

Definition at line 69 of file MuTriggerAnalyzer.cc.

Referenced by MuTriggerAnalyzer().

◆ ptMuCut_

double MuTriggerAnalyzer::ptMuCut_

private

Definition at line 60 of file MuTriggerAnalyzer.cc.

Referenced by analyze().

◆ trigEvTag_

edm::InputTag MuTriggerAnalyzer::trigEvTag_

private

Definition at line 57 of file MuTriggerAnalyzer.cc.

Referenced by analyze().

◆ trigEvToken_

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

private

Definition at line 58 of file MuTriggerAnalyzer.cc.

Referenced by analyze().

◆ triggerResults_

edm::Handle<edm::TriggerResults> MuTriggerAnalyzer::triggerResults_

private

Definition at line 63 of file MuTriggerAnalyzer.cc.

◆ trigNames_

edm::TriggerNames const* MuTriggerAnalyzer::trigNames_

private

Definition at line 64 of file MuTriggerAnalyzer.cc.

Referenced by analyze().

◆ trigToken_

edm::EDGetTokenT<TriggerResults> MuTriggerAnalyzer::trigToken_

private

Definition at line 56 of file MuTriggerAnalyzer.cc.

Referenced by analyze().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ MuTriggerAnalyzer()

Member Function Documentation

◆ analyze()

◆ endJob()

◆ IsMuMatchedToHLTMu()

Member Data Documentation

◆ deltaPtOverPt_

◆ deltaR_

◆ etaMuCut_

◆ handleTriggerEvent_

◆ hltPath_

◆ hTrigMuonPtDenS_

◆ hTrigMuonPtNumS_

◆ L3FilterName_

◆ maxDeltaR_

◆ maxDPtRel_

◆ muonToken_

◆ nbins_

◆ ptMax_

◆ ptMuCut_

◆ trigEvTag_

◆ trigEvToken_

◆ triggerResults_

◆ trigNames_

◆ trigToken_