#include <HLTMuonL2SelectorForL3IO.h>

Inheritance diagram for HLTMuonL2SelectorForL3IO:

Public Member Functions
	HLTMuonL2SelectorForL3IO (const edm::ParameterSet &)
	constructor with config More...

virtual void	produce (edm::Event &, const edm::EventSetup &)
	select muons More...

virtual	~HLTMuonL2SelectorForL3IO ()
	destructor More...

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

Public Member Functions inherited from edm::stream::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducerBase ()

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

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, std::unordered_multimap< std::string, edm::ProductResolverIndex > const &iIndicies, std::string const &moduleLabel)

virtual	~ProducerBase () noexcept(false)

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

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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

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

virtual	~EDConsumerBase () noexcept(false)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)
	default values More...

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Attributes
const bool	applyL3Filters_

const edm::EDGetTokenT< reco::TrackCollection >	l2Src_

const edm::EDGetTokenT< reco::MuonTrackLinksCollection >	l3linkToken_

const edm::EDGetTokenT< reco::RecoChargedCandidateCollection >	l3OISrc_

const double	max_NormalizedChi2_

const double	max_PtDifference_

const int	min_Nhits_

const int	min_NmuonHits_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Public Types inherited from edm::stream::EDProducerBase
typedef EDProducerAdaptorBase	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

L2 muon selector for L3 IO: finds L2 muons not previous converted into (good) L3 muons

Author: Benjamin Radburn-Smith, Santiago Folgueras - Purdue University

Definition at line 23 of file HLTMuonL2SelectorForL3IO.h.

Constructor & Destructor Documentation

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

constructor with config

Definition at line 17 of file HLTMuonL2SelectorForL3IO.cc.

References LogTrace.

                                                                                 :  
   l2Src_(consumes<reco::TrackCollection>(iConfig.getParameter<edm::InputTag>("l2Src"))),
   l3OISrc_(consumes<reco::RecoChargedCandidateCollection>(iConfig.getParameter<edm::InputTag>("l3OISrc"))),
   l3linkToken_(consumes<reco::MuonTrackLinksCollection>(iConfig.getParameter<edm::InputTag>("InputLinks"))),
   applyL3Filters_(iConfig.getParameter<bool>("applyL3Filters")),
   max_NormalizedChi2_ (iConfig.getParameter<double> ("MaxNormalizedChi2")),
   max_PtDifference_ (iConfig.getParameter<double> ("MaxPtDifference")),
   min_Nhits_ (iConfig.getParameter<int> ("MinNhits")),
   min_NmuonHits_ (iConfig.getParameter<int> ("MinNmuonHits"))  
 {
    LogTrace("Muon|RecoMuon|HLTMuonL2SelectorForL3IO")<<"constructor called";
    produces<reco::TrackCollection>();
 }

HLTMuonL2SelectorForL3IO::~HLTMuonL2SelectorForL3IO ( )

virtual

destructor

Definition at line 32 of file HLTMuonL2SelectorForL3IO.cc.

32 {}

Member Function Documentation

void HLTMuonL2SelectorForL3IO::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

default values

Definition at line 94 of file HLTMuonL2SelectorForL3IO.cc.

References edm::ConfigurationDescriptions::add(), and edm::ParameterSetDescription::add().

                                                                                           {
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("l2Src",edm::InputTag("hltL2Muons","UpdatedAtVtx"));
   desc.add<edm::InputTag>("l3OISrc",edm::InputTag("hltNewOIL3MuonCandidates"));
   desc.add<edm::InputTag>("InputLinks",edm::InputTag("hltNewOIL3MuonsLinksCombination"));
   desc.add<bool>("applyL3Filters",true); 
   desc.add<int>("MinNhits",1);
   desc.add<double>("MaxNormalizedChi2",20.0);
   desc.add<int>("MinNmuonHits",0);
   desc.add<double>("MaxPtDifference",999.0); //relative difference
   descriptions.add("HLTMuonL2SelectorForL3IO",desc);
 }

void HLTMuonL2SelectorForL3IO::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

virtual

select muons

create collection of L2 muons not already reconstructed as L3 muons

Definition at line 35 of file HLTMuonL2SelectorForL3IO.cc.

References funct::abs(), applyL3Filters_, reco::deltaR2(), reco::TrackBase::eta(), edm::Event::getByToken(), reco::MuonTrackLinks::globalTrack(), cutBasedMuonId_MuonPOG_V0_cff::globalTrack, reco::TrackBase::hitPattern(), checklumidiff::l, l2Src_, l3linkToken_, l3OISrc_, max_NormalizedChi2_, max_PtDifference_, metname, min_Nhits_, min_NmuonHits_, eostools::move(), reco::TrackBase::normalizedChi2(), reco::TrackBase::numberOfValidHits(), reco::HitPattern::numberOfValidMuonHits(), reco::TrackBase::phi(), reco::TrackBase::pt(), edm::Event::put(), mps_fire::result, reco::MuonTrackLinks::standAloneTrack(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                    {
   const std::string metname = "Muon|RecoMuon|HLTMuonL2SelectorForL3IO";
     
   //  IN:
   edm::Handle<reco::TrackCollection> l2muonH;
   iEvent.getByToken(l2Src_, l2muonH);
 
   edm::Handle<reco::RecoChargedCandidateCollection> l3muonH;
   iEvent.getByToken(l3OISrc_,l3muonH);
 
   // Read Links collection:
   edm::Handle<reco::MuonTrackLinksCollection> links;
   iEvent.getByToken(l3linkToken_, links);  
   
   //    OUT:
   std::unique_ptr<reco::TrackCollection > result(new reco::TrackCollection());
     
   for (unsigned int il2=0; il2 != l2muonH->size(); ++il2){
     reco::TrackRef l2muRef(l2muonH, il2);
     bool re_do_this_L2=true;
     
     for (unsigned int il3=0; il3 != l3muonH->size(); ++il3){
       reco::RecoChargedCandidateRef cand(l3muonH,il3);
       reco::TrackRef tk = cand->track();
       
       bool useThisLink = false;
       for (unsigned int l(0); l<links->size() && !useThisLink; ++l){
     const reco::MuonTrackLinks* link = &links->at(l);
     
     // Check if the L3 link matches the L3 candidate
     const reco::Track& globalTrack = *link->globalTrack();
     float dR2 = deltaR2(tk->eta(),tk->phi(),globalTrack.eta(),globalTrack.phi());
     if (dR2 < 0.02*0.02 and std::abs(tk->pt() - globalTrack.pt()) < 0.001*tk->pt()) {
       useThisLink=true;
     }
     
     if (!useThisLink) continue;
     
     // Check whether the stand-alone track matches a L2, if not, we will re-use this L2
     const reco::TrackRef staTrack = link->standAloneTrack();    
     if (l2muRef==staTrack) re_do_this_L2=false; 
       
     // Check the quality of the reconstructed L3, if poor quality, we will re-use this L2
     if (staTrack==l2muRef && applyL3Filters_) {
       re_do_this_L2 = true;
       const reco::Track& globalTrack = *link->globalTrack();
       if (globalTrack.numberOfValidHits()< min_Nhits_) continue;    // cut on number of hits
       if (globalTrack.normalizedChi2() > max_NormalizedChi2_ ) continue;        //normalizedChi2 cut
       if (globalTrack.hitPattern().numberOfValidMuonHits() < min_NmuonHits_ ) continue;     //min muon hits cut
       if (std::abs(globalTrack.pt()-l2muRef->pt()) > max_PtDifference_*globalTrack.pt()) continue; // pt difference 
       re_do_this_L2 = false;
     }
       }
     }
     if (re_do_this_L2) result->push_back(*l2muRef);  // used the L2 if no L3 if matched or if the matched L3 has poor quality cuts. 
   }
   iEvent.put(std::move(result));
 }

Member Data Documentation

const bool HLTMuonL2SelectorForL3IO::applyL3Filters_

private

Definition at line 41 of file HLTMuonL2SelectorForL3IO.h.

Referenced by produce().

const edm::EDGetTokenT<reco::TrackCollection> HLTMuonL2SelectorForL3IO::l2Src_

private

Definition at line 38 of file HLTMuonL2SelectorForL3IO.h.

Referenced by produce().

const edm::EDGetTokenT<reco::MuonTrackLinksCollection> HLTMuonL2SelectorForL3IO::l3linkToken_

private

Definition at line 40 of file HLTMuonL2SelectorForL3IO.h.

Referenced by produce().

const edm::EDGetTokenT<reco::RecoChargedCandidateCollection> HLTMuonL2SelectorForL3IO::l3OISrc_

private

Definition at line 39 of file HLTMuonL2SelectorForL3IO.h.

Referenced by produce().

const double HLTMuonL2SelectorForL3IO::max_NormalizedChi2_

private

Definition at line 42 of file HLTMuonL2SelectorForL3IO.h.

Referenced by produce().

const double HLTMuonL2SelectorForL3IO::max_PtDifference_

private

Definition at line 42 of file HLTMuonL2SelectorForL3IO.h.

Referenced by produce().

const int HLTMuonL2SelectorForL3IO::min_Nhits_

private

Definition at line 43 of file HLTMuonL2SelectorForL3IO.h.

Referenced by produce().

const int HLTMuonL2SelectorForL3IO::min_NmuonHits_

private

Definition at line 43 of file HLTMuonL2SelectorForL3IO.h.

Referenced by produce().

Public Member Functions

Static Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation