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 () override

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)

Public Attributes
L1Analysis::L1AnalysisRecoMuon2 *	muon

L1Analysis::L1AnalysisRecoMuon2DataFormat *	muon_data

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

void	beginJob (void) override

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

void	endJob () override

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

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

Constructor & Destructor Documentation

◆ L1Muon2RecoTreeProducer()

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

explicit

Definition at line 144 of file L1Muon2RecoTreeProducer.cc.

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

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

◆ ~L1Muon2RecoTreeProducer()

L1Muon2RecoTreeProducer::~L1Muon2RecoTreeProducer ( )

override

Definition at line 182 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

◆ analyze()

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

overrideprivatevirtual

Implements edm::EDAnalyzer.

Definition at line 192 of file L1Muon2RecoTreeProducer.cc.

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

References trigger::TriggerEvent::getObjects(), iEvent, edm::HandleBase::isValid(), hltdqm::passTrig(), hltMonBTagIPClient_cfi::pathName, AlCaHLTBitMon_QueryRunRegistry::string, triggerMatchMonitor_cfi::triggerObjects, triggerResults, trigNames, and pwdgSkimBPark_cfi::vertices.

◆ beginJob()

void L1Muon2RecoTreeProducer::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 306 of file L1Muon2RecoTreeProducer.cc.

306 {}

◆ beginRun()

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

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 355 of file L1Muon2RecoTreeProducer.cc.

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

References cms::cuda::assert(), gather_cfg::cout, createfilelist::int, and writedatasetfile::run.

◆ empty_hlt()

void L1Muon2RecoTreeProducer::empty_hlt ( )

Definition at line 311 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);
 }

◆ endJob()

void L1Muon2RecoTreeProducer::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 309 of file L1Muon2RecoTreeProducer.cc.

309 {}

◆ match_trigger()

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

Definition at line 318 of file L1Muon2RecoTreeProducer.cc.

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