#include <HLTMuonL1Filter.h>

Inheritance diagram for HLTMuonL1Filter:

Public Member Functions
bool	hltFilter (edm::Event &, const edm::EventSetup &, trigger::TriggerFilterObjectWithRefs &filterproduct) const override

	HLTMuonL1Filter (const edm::ParameterSet &)

	~HLTMuonL1Filter () override

Public Member Functions inherited from HLTFilter
	HLTFilter (const edm::ParameterSet &config)

int	module (edm::Event const &) const

const std::string *	moduleLabel () const

int	path (edm::Event const &) const

const std::string *	pathName (edm::Event const &) const

std::pair< int, int >	pmid (edm::Event const &) const

bool	saveTags () const

	~HLTFilter () override

Public Member Functions inherited from edm::global::EDFilter<>
	EDFilter ()=default

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsStreamLuminosityBlocks () const final

bool	wantsStreamRuns () const final

Public Member Functions inherited from edm::global::EDFilterBase
	EDFilterBase ()

ModuleDescription const &	moduleDescription () const

	~EDFilterBase () override

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

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

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, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) 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)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Static Public Member Functions inherited from HLTFilter
static void	makeHLTFilterDescription (edm::ParameterSetDescription &desc)

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
bool	isSingleSegmentCSC (const l1extra::L1MuonParticleRef &muon, L1CSCTrackCollection const &csctfTracks, L1MuTriggerScales const &scales) const
	checks if the passed L1MuExtraParticle is a single segment CSC More...

Private Attributes
edm::InputTag	candTag_
	input tag identifying the product containing muons More...

edm::EDGetTokenT< l1extra::L1MuonParticleCollection >	candToken_

edm::InputTag	csctfTag_
	input tag identifying the product containing CSCTF tracks More...

edm::EDGetTokenT< L1CSCTrackCollection >	csctfToken_

bool	excludeSingleSegmentCSC_
	should we exclude single-segment CSC trigger objects from our counting? More...

double	maxEta_
	max Eta cut More...

int	minN_
	required number of passing candidates to pass the filter More...

double	minPt_
	pT threshold More...

edm::InputTag	previousCandTag_
	input tag identifying the product containing refs to muons passing the previous level More...

edm::EDGetTokenT< trigger::TriggerFilterObjectWithRefs >	previousCandToken_

int	qualityBitMask_

Additional Inherited Members
Public Types inherited from edm::global::EDFilterBase
typedef EDFilterBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

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

See header file for documentation

Author: J. Alcaraz

This class is an HLTFilter (-> EDFilter) implementing a filter on L1 GMT input

Author: J. Alcaraz

Definition at line 29 of file HLTMuonL1Filter.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 32 of file HLTMuonL1Filter.cc.

References candTag_, edm::errors::Configuration, csctfTag_, edm::InputTag::encode(), Exception, excludeSingleSegmentCSC_, mps_fire::i, edm::isDebugEnabled(), LogDebug, maxEta_, minN_, minPt_, previousCandTag_, qualityBitMask_, HLTFilter::saveTags(), and ~HLTMuonL1Filter().

                                                                : HLTFilter(iConfig),
   candTag_( iConfig.getParameter<edm::InputTag>("CandTag") ),
   candToken_( consumes<l1extra::L1MuonParticleCollection>(candTag_)),
   previousCandTag_( iConfig.getParameter<edm::InputTag>("PreviousCandTag") ),
   previousCandToken_( consumes<trigger::TriggerFilterObjectWithRefs>(previousCandTag_)),
   maxEta_( iConfig.getParameter<double>("MaxEta") ),
   minPt_( iConfig.getParameter<double>("MinPt") ),
   minN_( iConfig.getParameter<int>("MinN") ),
   excludeSingleSegmentCSC_( iConfig.getParameter<bool>("ExcludeSingleSegmentCSC") ),
   csctfTag_( iConfig.getParameter<edm::InputTag>("CSCTFtag") ),
   csctfToken_(excludeSingleSegmentCSC_ ? consumes<L1CSCTrackCollection>(csctfTag_) : edm::EDGetTokenT<L1CSCTrackCollection>{})
 {
   using namespace std;
 
   //set the quality bit mask
   qualityBitMask_ = 0;
   vector<int> selectQualities = iConfig.getParameter<vector<int> >("SelectQualities");
   for(int selectQualitie : selectQualities){
     if(selectQualitie > 7){
       throw edm::Exception(edm::errors::Configuration) << "QualityBits must be smaller than 8!";
     }
     qualityBitMask_ |= 1<<selectQualitie;
   }
 
   // dump parameters for debugging
   if(edm::isDebugEnabled()){
     ostringstream ss;
     ss<<"Constructed with parameters:"<<endl;
     ss<<"    CandTag = "<<candTag_.encode()<<endl;
     ss<<"    PreviousCandTag = "<<previousCandTag_.encode()<<endl;
     ss<<"    MaxEta = "<<maxEta_<<endl;
     ss<<"    MinPt = "<<minPt_<<endl;
     ss<<"    SelectQualities =";
     for(size_t i=0; i<8; i++){
       if((qualityBitMask_>>i) % 2) ss<<" "<<i;
     }
     ss<<endl;
     ss<<"    MinN = "<<minN_<<endl;
     ss<<"    ExcludeSingleSegmentCSC = "<<excludeSingleSegmentCSC_<<endl;
     ss<<"    CSCTFtag = "<<csctfTag_.encode()<<endl;
     ss<<"    saveTags= "<<saveTags();
     LogDebug("HLTMuonL1Filter")<<ss.str();
   }
 }

HLTMuonL1Filter::~HLTMuonL1Filter ( )

overridedefault

Referenced by HLTMuonL1Filter().

Member Function Documentation

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

static

Definition at line 80 of file HLTMuonL1Filter.cc.

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

                                                                             {
   edm::ParameterSetDescription desc;
   makeHLTFilterDescription(desc);
   desc.add<edm::InputTag>("CandTag",edm::InputTag("hltL1extraParticles"));
   //  desc.add<edm::InputTag>("PreviousCandTag",edm::InputTag("hltL1sL1DoubleMuOpen"));
   desc.add<edm::InputTag>("PreviousCandTag",edm::InputTag(""));
   desc.add<double>("MaxEta",2.5);
   desc.add<double>("MinPt",0.0);
   desc.add<int>("MinN",1);
   desc.add<bool>("ExcludeSingleSegmentCSC",false);
   //  desc.add<edm::InputTag>("CSCTFtag",edm::InputTag("unused"));
   desc.add<edm::InputTag>("CSCTFtag",edm::InputTag("csctfDigis"));
   {
     std::vector<int> temp1;
     temp1.reserve(0);
     desc.add<std::vector<int> >("SelectQualities",temp1);
   }
   descriptions.add("hltMuonL1Filter",desc);
 }

bool HLTMuonL1Filter::hltFilter	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup,
		trigger::TriggerFilterObjectWithRefs &	filterproduct
	)		const

overridevirtual

handle for CSCTFtracks

Implements HLTFilter.

Definition at line 105 of file HLTMuonL1Filter.cc.

References accept(), trigger::TriggerFilterObjectWithRefs::addCollectionTag(), trigger::TriggerRefsCollections::addObject(), allMuons_cfi::allMuons, candTag_, candToken_, csctfToken_, excludeSingleSegmentCSC_, spr::find(), alignBH_cfg::fixed, edm::EventSetup::get(), edm::Event::getByToken(), trigger::TriggerRefsCollections::getObjects(), mps_fire::i, edm::isDebugEnabled(), isSingleSegmentCSC(), trigger::TriggerRefsCollections::l1muonSize(), LogDebug, maxEta_, minN_, minPt_, RPCpg::mu, metsig::muon, gen::n, previousCandToken_, edm::Handle< T >::product(), edm::ESHandle< T >::product(), jets_cff::quality, qualityBitMask_, HLTFilter::saveTags(), and trigger::TriggerL1Mu.

                                                                                                                                        {
   using namespace std;
   using namespace edm;
   using namespace trigger;
   using namespace l1extra;
 
   // All HLT filters must create and fill an HLT filter object,
   // recording any reconstructed physics objects satisfying (or not)
   // this HLT filter, and place it in the Event.
 
   // get hold of all muons
   Handle<L1MuonParticleCollection> allMuons;
   iEvent.getByToken(candToken_, allMuons);
 
   const L1CSCTrackCollection * csctfTracks = nullptr;
   const L1MuTriggerScales * l1MuTriggerScales = nullptr;
 
   // get hold of CSCTF raw tracks
   if( excludeSingleSegmentCSC_ ) {
     edm::Handle<L1CSCTrackCollection> csctfTracksHandle;
     iEvent.getByToken(csctfToken_, csctfTracksHandle);
     csctfTracks = csctfTracksHandle.product();
 
     // read scales for every event (fast, no need to cache this)
     ESHandle<L1MuTriggerScales> scales;
     iSetup.get<L1MuTriggerScalesRcd>().get(scales);
     l1MuTriggerScales = scales.product();
   }
 
   // get hold of muons that fired the previous level
   Handle<TriggerFilterObjectWithRefs> previousLevelCands;
   iEvent.getByToken(previousCandToken_, previousLevelCands);
   vector<L1MuonParticleRef> prevMuons;
   previousLevelCands->getObjects(TriggerL1Mu, prevMuons);
 
   // look at all muon candidates, check cuts and add to filter object
   int n = 0;
   for(size_t i = 0; i < allMuons->size(); i++){
     L1MuonParticleRef muon(allMuons, i);
 
     //check if triggered by the previous level
     if(find(prevMuons.begin(), prevMuons.end(), muon) == prevMuons.end()) continue;
 
     //check maxEta cut
     if(fabs(muon->eta()) > maxEta_) continue;
 
     //check pT cut
     if(muon->pt() < minPt_) continue;
 
     //check quality cut
     if(qualityBitMask_){
       int quality = muon->gmtMuonCand().empty() ? 0 : (1 << muon->gmtMuonCand().quality());
       if((quality & qualityBitMask_) == 0) continue;
     }
 
     // reject single-segment CSC objects if necessary
     if (excludeSingleSegmentCSC_ and isSingleSegmentCSC(muon, * csctfTracks, * l1MuTriggerScales)) continue;
 
     //we have a good candidate
     n++;
     filterproduct.addObject(TriggerL1Mu,muon);
   }
 
   if (saveTags()) filterproduct.addCollectionTag(candTag_);
 
   // filter decision
   const bool accept(n >= minN_);
 
   // dump event for debugging
   if(edm::isDebugEnabled()){
     ostringstream ss;
     ss.precision(2);
     ss<<"L1mu#"<<'\t'<<"q*pt"<<'\t'<<'\t'<<"eta"<<'\t'<<"phi"<<'\t'<<"quality"<<'\t'<<"isPrev"<<'\t'<<"isFired"<<'\t'<<"isSingleCSC"<<endl;
     ss<<"--------------------------------------------------------------------------"<<endl;
 
     vector<L1MuonParticleRef> firedMuons;
     filterproduct.getObjects(TriggerL1Mu, firedMuons);
     for(size_t i=0; i<allMuons->size(); i++){
       L1MuonParticleRef mu(allMuons, i);
       int quality = mu->gmtMuonCand().empty() ? 0 : mu->gmtMuonCand().quality();
       bool isPrev = find(prevMuons.begin(), prevMuons.end(), mu) != prevMuons.end();
       bool isFired = find(firedMuons.begin(), firedMuons.end(), mu) != firedMuons.end();
       bool isSingleCSC = excludeSingleSegmentCSC_ and isSingleSegmentCSC(mu, * csctfTracks, * l1MuTriggerScales);
       ss<<i<<'\t'<<scientific<<mu->charge()*mu->pt()<<'\t'<<fixed<<mu->eta()<<'\t'<<mu->phi()<<'\t'<<quality<<'\t'<<isPrev<<'\t'<<isFired<<'\t'<<isSingleCSC<<endl;
     }
     ss<<"--------------------------------------------------------------------------"<<endl;
     LogDebug("HLTMuonL1Filter")<<ss.str()<<"Decision of filter is "<<accept<<", number of muons passing = "<<filterproduct.l1muonSize();
   }
 
   return accept;
 }

bool HLTMuonL1Filter::isSingleSegmentCSC	(	const l1extra::L1MuonParticleRef &	muon,
		L1CSCTrackCollection const &	csctfTracks,
		L1MuTriggerScales const &	scales
	)		const

private

checks if the passed L1MuExtraParticle is a single segment CSC

Definition at line 198 of file HLTMuonL1Filter.cc.

References DEFINE_FWK_MODULE, L1MuScale::getCenter(), L1MuScale::getLowEdge(), L1MuTriggerScales::getPhiScale(), L1MuTriggerScales::getRegionalEtaScale(), and Pi.

Referenced by hltFilter().

                                                                                                                                                                            {
   // is the muon matching a csctf track?
   //bool matched   = false;     // unused
   // which csctf track mode?
   // -999: no matching
   //  1: bad phi road. Not good extrapolation, but still triggering
   // 11: singles
   // 15: halo
   // 2->10 and 12->14: coincidence trigger with good extrapolation
   int csctfMode = -999;
 
   // loop over the CSCTF tracks
   for(auto trk = csctfTracks.begin(); trk < csctfTracks.end(); trk++){
 
     int trEndcap = (trk->first.endcap()==2 ? trk->first.endcap()-3 : trk->first.endcap());
     int trSector = 6*(trk->first.endcap()-1)+trk->first.sector();
 
     //... in radians
     // Type 2 is CSC
     float trEtaScale = l1MuTriggerScales.getRegionalEtaScale(2)->getCenter(trk->first.eta_packed());
     float trPhiScale = l1MuTriggerScales.getPhiScale()->getLowEdge(trk->first.localPhi());
 
     double trEta = trEtaScale * trEndcap;
     // there is no L1ExtraParticle below -2.375
     if(trEta<-2.4) trEta=-2.375;
 
     // CSCTF has 6 sectors
     // sector 1 starts at 15 degrees
     // trPhiScale is defined inside a sector
     float trPhi02PI = fmod(trPhiScale +
                            ((trSector-1)*TMath::Pi()/3) +
                            (TMath::Pi()/12) , 2*TMath::Pi());
 
     // L1 information are given from [-Pi,Pi]
     double trPhi = ( trPhi02PI<TMath::Pi()? trPhi02PI : trPhi02PI - 2*TMath::Pi() );
     /*
     std::cout << "\ntrEndcap="               << trEndcap                << std::endl;
     std::cout << "trSector="                 << trSector                << std::endl;
     std::cout << "trk->first.eta_packed()="  << trk->first.eta_packed() << std::endl;
     std::cout << "trk->first.localPhi()="    << trk->first.localPhi()   << std::endl;
     std::cout << "trEtaScale=" << trEtaScale << std::endl;
     std::cout << "trPhiScale=" << trPhiScale << std::endl;
     std::cout << "trEta="      << trEta      << std::endl;
     std::cout << "trPhi="      << trPhi      << std::endl;
     */
     if ( fabs (trEta-muon->eta()) < 0.03 &&
          fabs (trPhi-muon->phi()) < 0.001  ) {
 
       //matched = true;
       ptadd thePtAddress(trk->first.ptLUTAddress());
       csctfMode = thePtAddress.track_mode;
       //std::cout << "is matched -> trMode=" << csctfMode << std::endl;
     }
   }
 
   /*
   std::cout << " ===================================== " << std::endl;
   std::cout << " is matched? " << matched                << std::endl;
   std::cout << " is singles? " << (csctfMode==11 ? 1 :0) << std::endl;
   std::cout << " ===================================== " << std::endl;
   */
 
   // singles are mode 11 "CSCTF tracks"
   return csctfMode==11;
 }

Member Data Documentation

edm::InputTag HLTMuonL1Filter::candTag_

private

input tag identifying the product containing muons

Definition at line 39 of file HLTMuonL1Filter.h.

Referenced by hltFilter(), and HLTMuonL1Filter().

edm::EDGetTokenT<l1extra::L1MuonParticleCollection> HLTMuonL1Filter::candToken_

private

Definition at line 40 of file HLTMuonL1Filter.h.

Referenced by hltFilter().

edm::InputTag HLTMuonL1Filter::csctfTag_

private

input tag identifying the product containing CSCTF tracks

Definition at line 82 of file HLTMuonL1Filter.h.

Referenced by HLTMuonL1Filter().

edm::EDGetTokenT<L1CSCTrackCollection> HLTMuonL1Filter::csctfToken_

private

Definition at line 83 of file HLTMuonL1Filter.h.

Referenced by hltFilter().

bool HLTMuonL1Filter::excludeSingleSegmentCSC_

private

should we exclude single-segment CSC trigger objects from our counting?

Definition at line 76 of file HLTMuonL1Filter.h.

Referenced by hltFilter(), and HLTMuonL1Filter().

double HLTMuonL1Filter::maxEta_

private

max Eta cut

Definition at line 47 of file HLTMuonL1Filter.h.

Referenced by hltFilter(), and HLTMuonL1Filter().

int HLTMuonL1Filter::minN_

private

required number of passing candidates to pass the filter

Definition at line 73 of file HLTMuonL1Filter.h.

Referenced by hltFilter(), and HLTMuonL1Filter().

double HLTMuonL1Filter::minPt_

private

pT threshold

Definition at line 50 of file HLTMuonL1Filter.h.

Referenced by hltFilter(), and HLTMuonL1Filter().

edm::InputTag HLTMuonL1Filter::previousCandTag_

private

input tag identifying the product containing refs to muons passing the previous level

Definition at line 43 of file HLTMuonL1Filter.h.

Referenced by HLTMuonL1Filter().

edm::EDGetTokenT<trigger::TriggerFilterObjectWithRefs> HLTMuonL1Filter::previousCandToken_

private

Definition at line 44 of file HLTMuonL1Filter.h.

Referenced by hltFilter().

int HLTMuonL1Filter::qualityBitMask_

private

Quality codes:

0 .. no muon 1 .. beam halo muon (CSC) 2 .. very low quality level 1 (e.g. ignore in single and di-muon trigger) 3 .. very low quality level 2 (e.g. ignore in single muon trigger use in di-muon trigger) 4 .. very low quality level 3 (e.g. ignore in di-muon trigger, use in single-muon trigger) 5 .. unmatched RPC 6 .. unmatched DT or CSC 7 .. matched DT-RPC or CSC-RPC

attention: try not to rely on quality codes in analysis: they may change again

Quality bit mask:

the eight lowest order or least significant bits correspond to the qulity codes above; if a bit is 1, that code is accepted, otherwise not; example: 11101000 accepts qualities 3, 5, 6, 7

Definition at line 70 of file HLTMuonL1Filter.h.

Referenced by hltFilter(), and HLTMuonL1Filter().

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation