Inheritance diagram for L1Muon2RecoTreeProducer:

Public Member Functions
void	empty_hlt ()

	L1Muon2RecoTreeProducer (const edm::ParameterSet &)

double	match_trigger (std::vector< int > &trigIndices, const trigger::TriggerObjectCollection &trigObjs, edm::Handle< trigger::TriggerEvent > &triggerEvent, const reco::Muon &mu)

	~L1Muon2RecoTreeProducer ()

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

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

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

	EDConsumerBase ()

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

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

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

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const

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

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Public Attributes
L1Analysis::L1AnalysisRecoMuon2 *	muon

L1Analysis::L1AnalysisRecoMuon2DataFormat *	muon_data

Private Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)

virtual void	beginJob (void)

virtual void	beginRun (const edm::Run &, const edm::EventSetup &)

virtual void	endJob ()

Private Attributes
edm::Service< TFileService >	fs_

HLTConfigProvider	hltConfig_

std::vector< int >	isoTriggerIndices_

std::vector< std::string >	isoTriggerNames_

edm::Handle< std::vector < std::string > >	isoTriggerNamesToken_

edm::Handle< edm::TriggerResults >	isoTriggerToken_

unsigned int	maxMuon_

edm::EDGetTokenT < reco::PFMETCollection >	metToken_

edm::EDGetTokenT < reco::MuonCollection >	MuonToken_

TTree *	tree_

std::vector< int >	triggerIndices_

bool	triggerMatching_

double	triggerMaxDeltaR_

std::vector< std::string >	triggerNames_

std::string	triggerProcessLabel_

edm::EDGetTokenT < edm::TriggerResults >	TriggerResultsToken_

edm::EDGetTokenT < trigger::TriggerEvent >	triggerSummaryLabelToken_

edm::EDGetTokenT < reco::VertexCollection >	VtxToken_

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

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 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 87 of file L1Muon2RecoTreeProducer.cc.

Constructor & Destructor Documentation

L1Muon2RecoTreeProducer::L1Muon2RecoTreeProducer ( const edm::ParameterSet & iConfig )

explicit

Definition at line 149 of file L1Muon2RecoTreeProducer.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), metsig::muon, and AlCaHLTBitMon_QueryRunRegistry::string.

 {
 
   maxMuon_         = iConfig.getParameter<unsigned int>("maxMuon");
   isoTriggerNames_         = iConfig.getParameter<std::vector<std::string>>("isoTriggerNames");
   // isoTriggerToken_         = iConfig.getParameter<std::vector<std::string>>("isoTriggerNames");
   // TriggerToken_         = iConfig.getParameter<std::vector<std::string>>("triggerNames");
   MuonToken_ = consumes<reco::MuonCollection>(iConfig.getUntrackedParameter("MuonToken",edm::InputTag("muons")));
   VtxToken_  = consumes<reco::VertexCollection>(iConfig.getUntrackedParameter("VertexToken",edm::InputTag("offlinePrimaryVertices"))); 
     
   TriggerResultsToken_      = consumes<edm::TriggerResults>(iConfig.getUntrackedParameter("TriggerResultsToken",edm::InputTag("TriggerResults","","HLT")));
   triggerSummaryLabelToken_ = consumes<trigger::TriggerEvent>(iConfig.getUntrackedParameter("triggerSummaryLabelToken",edm::InputTag("hltTriggerSummaryAOD","","HLT")));
   // TriggerResultsToken_      = consumes<edm::TriggerResults>(iConfig.getUntrackedParameter("TriggerResultsToken",edm::InputTag("triggerSummaryLabel")));
   // triggerSummaryLabelToken_ = consumes<trigger::TriggerEvent>(iConfig.getUntrackedParameter("triggerSummaryLabelToken",edm::InputTag("triggerSummaryLabel")));
     //iConfig.getParameter<edm::InputTag>("triggerSummaryLabel");
   metToken_               = consumes<reco::PFMETCollection>(iConfig.getUntrackedParameter("metToken",edm::InputTag("pfMet")));
 
   muon           = new L1Analysis::L1AnalysisRecoMuon2(iConfig);
   muon_data           = muon->getData();
 
   tree_=fs_->make<TTree>("Muon2RecoTree", "Muon2RecoTree");
   tree_->Branch("Muon",           "L1Analysis::L1AnalysisRecoMuon2DataFormat",         &muon_data,                32000, 3);
 
   triggerMaxDeltaR_    = iConfig.getParameter<double>("triggerMaxDeltaR");
   triggerMatching_     = iConfig.getUntrackedParameter<bool>("triggerMatching");
   triggerProcessLabel_ = iConfig.getUntrackedParameter<std::string>("triggerProcessLabel");
   isoTriggerNames_ = iConfig.getParameter<std::vector<std::string>>("isoTriggerNames");
   triggerNames_ = iConfig.getParameter<std::vector<std::string>>("triggerNames");
 
   // triggerMatching_     = iConfig.getUntrackedParameter<bool>("triggerMatching");
   // _ = iConfig.getParameter<edm::InputTag>("triggerSummaryLabel");
   // triggerProcessLabel_ = iConfig.getUntrackedParameter<std::string>("triggerProcessLabel");
   // triggerNames_        = iConfig.getParameter<std::vector<std::string> > ("triggerNames");
   // isoTriggerNames_     = iConfig.getParameter<std::vector<std::string> > ("isoTriggerNames");
   // triggerMaxDeltaR_    = iConfig.getParameter<double> ("triggerMaxDeltaR");
   
 }

L1Muon2RecoTreeProducer::~L1Muon2RecoTreeProducer ( )

Definition at line 188 of file L1Muon2RecoTreeProducer.cc.

 {
 
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)
 
 }

Member Function Documentation

void L1Muon2RecoTreeProducer::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 202 of file L1Muon2RecoTreeProducer.cc.

References edm::Event::getByToken(), edm::HandleBase::isValid(), metsig::muon, EgammaValidation_cff::pathName, edm::TriggerNames::size(), AlCaHLTBitMon_QueryRunRegistry::string, edm::TriggerNames::triggerIndex(), edm::Event::triggerNames(), triggerResults, trigNames, and HLT_25ns10e33_v2_cff::vertices.

 {
 
   muon->Reset();
   edm::Handle<reco::MuonCollection> recoMuons;
   iEvent.getByToken(MuonToken_, recoMuons);
 
   edm::Handle<reco::VertexCollection> vertices;
   iEvent.getByToken(VtxToken_, vertices);
 
   edm::Handle<edm::TriggerResults> triggerResults;
   iEvent.getByToken(TriggerResultsToken_, triggerResults);
 
   edm::Handle<trigger::TriggerEvent> triggerSummaryLabel_;
   iEvent.getByToken(triggerSummaryLabelToken_, triggerSummaryLabel_);
 
   edm::Handle<reco::PFMETCollection> metLabel_;
   iEvent.getByToken(metToken_, metLabel_);
 
   int counter_met = 0;
 
   double METx = 0.;
   double METy = 0.;
 
   for(reco::PFMETCollection::const_iterator imet = metLabel_->begin(); 
       imet != metLabel_->end() && (unsigned) counter_met < 1; imet++) {
     
     METx = imet->px();
     METy = imet->py();
     
   }
 
   if (recoMuons.isValid()) {
     muon->SetMuon(iEvent, iSetup, recoMuons, vertices, METx, METy, maxMuon_);
   }
   else {
 
   }
 
 
   int counter_mu = 0;
   for(reco::MuonCollection::const_iterator imu = recoMuons->begin(); 
       imu != recoMuons->end() && (unsigned) counter_mu < maxMuon_; imu++) {
 
     //---------------------------------------------------------------------
     // TRIGGER MATCHING:
     // if specified the reconstructed muons are matched to a trigger
     //---------------------------------------------------------------------
     if (triggerMatching_) {
       double isoMatchDeltaR = 9999.;
       double matchDeltaR = 9999.;
       int hasIsoTriggered = 0;
       int hasTriggered = 0;
 
       int passesSingleMuonFlag = 0;
       
       // first check if the trigger results are valid:
       if (triggerResults.isValid()) {   
 
     if (triggerSummaryLabel_.isValid()) {   
 
     const edm::TriggerNames& trigNames = iEvent.triggerNames(*triggerResults); 
     // for(UInt_t iPath = 0 ; iPath < trigNames.size() ; ++iPath)
     //   {
     //     cout<<iPath<<": "<<trigNames.triggerName(iPath)<<endl;
     //   }
     
     for(UInt_t iPath = 0 ; iPath < isoTriggerNames_.size() ; ++iPath)
       {
         if(passesSingleMuonFlag==1) continue;
         std::string pathName=isoTriggerNames_.at(iPath);
 
         bool passTrig=false;
         //cout<<"testing pathName                         = "<<pathName<<endl;
         //cout<<"trigNames.triggerIndex(pathName) = "<<trigNames.triggerIndex(pathName)<<endl;
         // cout<<"trigNames.triggerIndex(pathName)<trigNames.size()= "<<(trigNames.triggerIndex(pathName)<trigNames.size())<<endl;
         
         if(trigNames.triggerIndex(pathName)<trigNames.size()) passTrig=triggerResults->accept(trigNames.triggerIndex(pathName));
         // cout<<"pass = "<<passTrig<<endl;
         if(passTrig) passesSingleMuonFlag=1;
       }   
     
     if (triggerSummaryLabel_.isValid()) {
       // get trigger objects:
       const trigger::TriggerObjectCollection triggerObjects = triggerSummaryLabel_->getObjects();
 
       matchDeltaR = match_trigger(triggerIndices_, triggerObjects, triggerSummaryLabel_, (*imu));
       if (matchDeltaR < triggerMaxDeltaR_)
         hasTriggered = 1;
 
       // cout<<"isoTriggerIndices_.size() = "<<isoTriggerIndices_.size()<<endl;
       // cout<<"triggerObjects.size() = "<<triggerObjects.size()<<endl;
       // cout<<"imu->eta() = "<<imu->eta()<<endl;
 
       isoMatchDeltaR = match_trigger(isoTriggerIndices_, triggerObjects, triggerSummaryLabel_, (*imu));
       if (isoMatchDeltaR < triggerMaxDeltaR_) 
         hasIsoTriggered = 1;
       
     } // end if (triggerEvent.isValid())
     
       } // end if (triggerResults.isValid())
 
       } // end if (triggerResults.isValid())
 
         // fill trigger matching variables:
       muon_data->hlt_isomu.push_back(hasIsoTriggered);
       muon_data->hlt_mu.push_back(hasTriggered);
       muon_data->hlt_isoDeltaR.push_back(isoMatchDeltaR);
       muon_data->hlt_deltaR.push_back(matchDeltaR);
       muon_data->passesSingleMuon.push_back(passesSingleMuonFlag);
       
     } else {
       empty_hlt();
     } // end if (triggerMatching_)
   }
   
 
   tree_->Fill();
 
 }

void L1Muon2RecoTreeProducer::beginJob ( void )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 325 of file L1Muon2RecoTreeProducer.cc.

326 {

327 }

void L1Muon2RecoTreeProducer::beginRun	(	const edm::Run &	run,
		const edm::EventSetup &	eventSetup
	)

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 380 of file L1Muon2RecoTreeProducer.cc.

References assert(), gather_cfg::cout, and metsig::muon.

                                                                                          {
   // Prepare for trigger matching for each new run: 
   // Look up triggetIndices in the HLT config for the different paths
   if (triggerMatching_) {
     bool changed = true;
     if (!hltConfig_.init(run, eventSetup, triggerProcessLabel_, changed)) {
       // if you can't initialize hlt configuration, crash!
       std::cout << "Error: didn't find process" << triggerProcessLabel_ << std::endl;
       assert(false);
     }
 
     bool enableWildcard = true;
     for (size_t iTrig = 0; iTrig < triggerNames_.size(); ++iTrig) { 
       // prepare for regular expression (with wildcards) functionality:
       TString tNameTmp = TString(triggerNames_[iTrig]);
       TRegexp tNamePattern = TRegexp(tNameTmp, enableWildcard);
       int tIndex = -1;
       // find the trigger index:
       for (unsigned ipath = 0; ipath < hltConfig_.size(); ++ipath) {
         // use TString since it provides reg exp functionality:
         TString tmpName = TString(hltConfig_.triggerName(ipath));
         if (tmpName.Contains(tNamePattern)) {
           tIndex = int(ipath);
           triggerIndices_.push_back(tIndex);
         }
       }
       if (tIndex < 0) { // if can't find trigger path at all, give warning:
         std::cout << "Warning: Could not find trigger" << triggerNames_[iTrig] << std::endl;
         //assert(false);
       }
     } // end for triggerNames
     for (size_t iTrig = 0; iTrig < isoTriggerNames_.size(); ++iTrig) { 
 
       // prepare for regular expression functionality:
       TString tNameTmp = TString(isoTriggerNames_[iTrig]);
       TRegexp tNamePattern = TRegexp(tNameTmp, enableWildcard);
       int tIndex = -1;
       // find the trigger index:
       for (unsigned ipath = 0; ipath < hltConfig_.size(); ++ipath) {
         // use TString since it provides reg exp functionality:
         TString tmpName = TString(hltConfig_.triggerName(ipath));
         if (tmpName.Contains(tNamePattern)) {
           tIndex = int(ipath);
           isoTriggerIndices_.push_back(tIndex);
         }
       }
       if (tIndex < 0) { // if can't find trigger path at all, give warning:
         std::cout << "Warning: Could not find trigger" << isoTriggerNames_[iTrig] << std::endl;
         //assert(false);
       }
     } // end for isoTriggerNames
   } // end if (triggerMatching_)
 
   muon->init(eventSetup);
 }

void L1Muon2RecoTreeProducer::empty_hlt ( )

Definition at line 334 of file L1Muon2RecoTreeProducer.cc.

                                        {
 
   muon_data->hlt_isomu.push_back(-9999);
   muon_data->hlt_mu.push_back(-9999);
   muon_data->hlt_isoDeltaR.push_back(-9999);
   muon_data->hlt_deltaR.push_back(-9999);
   
 }

void L1Muon2RecoTreeProducer::endJob ( void )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 331 of file L1Muon2RecoTreeProducer.cc.

331 {

332 }

double L1Muon2RecoTreeProducer::match_trigger	(	std::vector< int > &	trigIndices,
		const trigger::TriggerObjectCollection &	trigObjs,
		edm::Handle< trigger::TriggerEvent > &	triggerEvent,
		const reco::Muon &	mu
	)

Definition at line 343 of file L1Muon2RecoTreeProducer.cc.

References deltaR().

 {
   double matchDeltaR = 9999;
 
   for(size_t iTrigIndex = 0; iTrigIndex < trigIndices.size(); ++iTrigIndex) {
     int triggerIndex = trigIndices[iTrigIndex];
     const std::vector<std::string> moduleLabels(hltConfig_.moduleLabels(triggerIndex));
     // find index of the last module:
     const unsigned moduleIndex = hltConfig_.size(triggerIndex)-2;
     // find index of HLT trigger name:
     const unsigned hltFilterIndex = triggerEvent->filterIndex( edm::InputTag ( moduleLabels[moduleIndex], "", triggerProcessLabel_ ) );
 
     if (hltFilterIndex < triggerEvent->sizeFilters()) {
       const trigger::Keys triggerKeys(triggerEvent->filterKeys(hltFilterIndex));
       const trigger::Vids triggerVids(triggerEvent->filterIds(hltFilterIndex));
 
       const unsigned nTriggers = triggerVids.size();
       for (size_t iTrig = 0; iTrig < nTriggers; ++iTrig) {
         // loop over all trigger objects:
         const trigger::TriggerObject trigObject = trigObjs[triggerKeys[iTrig]];
 
         double dRtmp = deltaR( mu, trigObject );
 
         if ( dRtmp < matchDeltaR ) {
           matchDeltaR = dRtmp;
         }
 
       } // loop over different trigger objects
     } // if trigger is in event (should apply hltFilter with used trigger...)
   } // loop over muon candidates
 
   return matchDeltaR;
 }