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 const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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 const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

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
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, 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)

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

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

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

Detailed Description

Definition at line 49 of file MuTriggerAnalyzer.cc.

Constructor & Destructor Documentation

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

Definition at line 113 of file MuTriggerAnalyzer.cc.

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

                                                                 :
   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_);
 
 }

Member Function Documentation

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

overrideprivate

Definition at line 144 of file MuTriggerAnalyzer.cc.

References edm::HLTGlobalStatus::accept(), gather_cfg::cout, DEFINE_FWK_MODULE, edm::InputTag::encode(), PVValHelper::eta, reco::LeafCandidate::eta(), etaMuCut_, trigger::TriggerEvent::filterKeys(), trigger::TriggerEvent::filterTag(), edm::Event::getByToken(), trigger::TriggerEvent::getObjects(), hltPath_, mps_fire::i, reco::Muon::isGlobalMuon(), IsMuMatchedToHLTMu(), gen::k, L3FilterName_, maxDeltaR_, maxDPtRel_, RPCpg::mu, extraflags_cff::muons, muonToken_, dataset::name, AlCaHLTBitMon_ParallelJobs::p, EnergyCorrector::pt, reco::LeafCandidate::pt(), ptMuCut_, edm::HLTGlobalStatus::size(), trigger::TriggerEvent::sizeFilters(), AlCaHLTBitMon_QueryRunRegistry::string, trigEvTag_, trigEvToken_, edm::TriggerNames::triggerName(), edm::Event::triggerNames(), triggerResults, trigNames, trigNames_, and trigToken_.

                                                                      {
           //      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_ );
       }
 
 }

void MuTriggerAnalyzer::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 135 of file MuTriggerAnalyzer.cc.

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

                                {
   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;
 
 
   }
 }

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

private

Definition at line 80 of file MuTriggerAnalyzer.cc.

References gather_cfg::cout, boostedElectronIsolation_cff::deltaR, deltaR_, reco::LeafCandidate::eta(), HcalObjRepresent::Fill(), gen::k, nJettinessAdder_cfi::nPass, EnergyCorrector::pt, and reco::LeafCandidate::pt().

Referenced by analyze().

                                                                                                                             {
   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);
 }

Member Data Documentation

TH1D* MuTriggerAnalyzer::deltaPtOverPt_

private

Definition at line 75 of file MuTriggerAnalyzer.cc.

Referenced by MuTriggerAnalyzer().

TH1D* MuTriggerAnalyzer::deltaR_

private

Definition at line 74 of file MuTriggerAnalyzer.cc.

Referenced by MuTriggerAnalyzer().

double MuTriggerAnalyzer::etaMuCut_

private

Definition at line 62 of file MuTriggerAnalyzer.cc.

Referenced by analyze().

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

private

Definition at line 67 of file MuTriggerAnalyzer.cc.

std::string MuTriggerAnalyzer::hltPath_

private

Definition at line 63 of file MuTriggerAnalyzer.cc.

Referenced by analyze().

TH1D* MuTriggerAnalyzer::hTrigMuonPtDenS_

private

Definition at line 73 of file MuTriggerAnalyzer.cc.

Referenced by endJob(), and MuTriggerAnalyzer().

TH1D* MuTriggerAnalyzer::hTrigMuonPtNumS_

private

Definition at line 72 of file MuTriggerAnalyzer.cc.

Referenced by endJob(), and MuTriggerAnalyzer().

std::string MuTriggerAnalyzer::L3FilterName_

private

Definition at line 64 of file MuTriggerAnalyzer.cc.

Referenced by analyze().

double MuTriggerAnalyzer::maxDeltaR_

private

Definition at line 69 of file MuTriggerAnalyzer.cc.

Referenced by analyze(), and MuTriggerAnalyzer().

double MuTriggerAnalyzer::maxDPtRel_

private

Definition at line 69 of file MuTriggerAnalyzer.cc.

Referenced by analyze(), and MuTriggerAnalyzer().

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

private

Definition at line 61 of file MuTriggerAnalyzer.cc.

Referenced by analyze().

const int MuTriggerAnalyzer::nbins_

private

Definition at line 70 of file MuTriggerAnalyzer.cc.

Referenced by endJob(), and MuTriggerAnalyzer().

const double MuTriggerAnalyzer::ptMax_

private

Definition at line 71 of file MuTriggerAnalyzer.cc.

Referenced by MuTriggerAnalyzer().

double MuTriggerAnalyzer::ptMuCut_

private

Definition at line 62 of file MuTriggerAnalyzer.cc.

Referenced by analyze().

edm::InputTag MuTriggerAnalyzer::trigEvTag_

private

Definition at line 59 of file MuTriggerAnalyzer.cc.

Referenced by analyze().

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

private

Definition at line 60 of file MuTriggerAnalyzer.cc.

Referenced by analyze().

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

private

Definition at line 65 of file MuTriggerAnalyzer.cc.

edm::TriggerNames const* MuTriggerAnalyzer::trigNames_

private

Definition at line 66 of file MuTriggerAnalyzer.cc.

Referenced by analyze().

edm::EDGetTokenT<TriggerResults> MuTriggerAnalyzer::trigToken_

private

Definition at line 58 of file MuTriggerAnalyzer.cc.

Referenced by analyze().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation