#include <HLTMuonDimuonL2Filter.h>

Inheritance diagram for HLTMuonDimuonL2Filter:

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

	HLTMuonDimuonL2Filter (const edm::ParameterSet &)

	~HLTMuonDimuonL2Filter () 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 Attributes
edm::InputTag	beamspotTag_

edm::EDGetTokenT< reco::BeamSpot >	beamspotToken_

edm::InputTag	candTag_

edm::EDGetTokenT< reco::RecoChargedCandidateCollection >	candToken_

int	chargeOpt_

bool	fast_Accept_

double	max_Acop_

double	max_Angle_

double	max_Dr_

double	max_Dz_

double	max_Eta_

double	max_InvMass_

double	max_PtBalance_

double	min_Acop_

double	min_Angle_

double	min_InvMass_

int	min_Nchambers_

int	min_Nhits_

int	min_Nstations_

double	min_PtBalance_

double	min_PtMax_

double	min_PtMin_

double	min_PtPair_

double	nsigma_Pt_

edm::InputTag	previousCandTag_

edm::EDGetTokenT< trigger::TriggerFilterObjectWithRefs >	previousCandToken_

edm::InputTag	seedMapTag_
	input tag of the map from the L2 seed to the sister L2 seeds of cleaned tracks More...

edm::EDGetTokenT< SeedMap >	seedMapToken_

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, P. Garcia, M.Vander Donckt

This class is an HLTFilter (-> EDFilter) implementing a muon pair filter for HLT muons

Author: J. Alcaraz

See header file for documentation

Author: J. Alcaraz, P. Garcia

Definition at line 23 of file HLTMuonDimuonL2Filter.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 40 of file HLTMuonDimuonL2Filter.cc.

References candTag_, chargeOpt_, edm::InputTag::encode(), fast_Accept_, LogDebug, max_Acop_, max_Angle_, max_Dr_, max_Dz_, max_Eta_, max_InvMass_, max_PtBalance_, min_Acop_, min_Angle_, min_InvMass_, min_Nchambers_, min_Nhits_, min_Nstations_, min_PtBalance_, min_PtMax_, min_PtMin_, min_PtPair_, nsigma_Pt_, and ~HLTMuonDimuonL2Filter().

                                                                           : HLTFilter(iConfig),
    beamspotTag_   (iConfig.getParameter< edm::InputTag > ("BeamSpotTag")),
    beamspotToken_ (consumes<reco::BeamSpot>(beamspotTag_)),
    candTag_     (iConfig.getParameter< edm::InputTag > ("CandTag")),
    candToken_   (consumes<reco::RecoChargedCandidateCollection>(candTag_)),
    previousCandTag_   (iConfig.getParameter<InputTag > ("PreviousCandTag")),
    previousCandToken_ (consumes<trigger::TriggerFilterObjectWithRefs>(previousCandTag_)),
    seedMapTag_( iConfig.getParameter<InputTag >("SeedMapTag") ),
    seedMapToken_(consumes<SeedMap>(seedMapTag_)),
    fast_Accept_ (iConfig.getParameter<bool> ("FastAccept")),
    max_Eta_     (iConfig.getParameter<double> ("MaxEta")),
    min_Nhits_   (iConfig.getParameter<int> ("MinNhits")),
    min_Nstations_(iConfig.getParameter<int> ("MinNstations")),
    min_Nchambers_(iConfig.getParameter<int> ("MinNchambers")),
    max_Dr_      (iConfig.getParameter<double> ("MaxDr")),
    max_Dz_      (iConfig.getParameter<double> ("MaxDz")),
    chargeOpt_   (iConfig.getParameter<int> ("ChargeOpt")),
    min_PtPair_  (iConfig.getParameter<double> ("MinPtPair")),
    min_PtMax_   (iConfig.getParameter<double> ("MinPtMax")),
    min_PtMin_   (iConfig.getParameter<double> ("MinPtMin")),
    min_InvMass_ (iConfig.getParameter<double> ("MinInvMass")),
    max_InvMass_ (iConfig.getParameter<double> ("MaxInvMass")),
    min_Acop_    (iConfig.getParameter<double> ("MinAcop")),
    max_Acop_    (iConfig.getParameter<double> ("MaxAcop")),
    min_Angle_   (iConfig.getParameter<double> ("MinAngle")),
    max_Angle_   (iConfig.getParameter<double> ("MaxAngle")),
    min_PtBalance_ (iConfig.getParameter<double> ("MinPtBalance")),
    max_PtBalance_ (iConfig.getParameter<double> ("MaxPtBalance")),
    nsigma_Pt_   (iConfig.getParameter<double> ("NSigmaPt"))
 {
 
    LogDebug("HLTMuonDimuonL2Filter")
       << " CandTag/MinN/MaxEta/MinNhits/MinNstations/MinNchambers/MaxDr/MaxDz/MinPt1/MinPt2/MinInvMass/MaxInvMass/MinAcop/MaxAcop/MinAngle/MaxAngle/MinPtBalance/MaxPtBalance/NSigmaPt : "
       << candTag_.encode()
       << " " << fast_Accept_
       << " " << max_Eta_
       << " " << min_Nhits_
       << " " << min_Nstations_
       << " " << min_Nchambers_
       << " " << max_Dr_
       << " " << max_Dz_
       << " " << chargeOpt_ << " " << min_PtPair_
       << " " << min_PtMax_ << " " << min_PtMin_
       << " " << min_InvMass_ << " " << max_InvMass_
       << " " << min_Acop_ << " " << max_Acop_
       << " " << min_Angle_ << " " << max_Angle_
       << " " << min_PtBalance_ << " " << max_PtBalance_
       << " " << nsigma_Pt_;
 }

HLTMuonDimuonL2Filter::~HLTMuonDimuonL2Filter ( )

overridedefault

Referenced by HLTMuonDimuonL2Filter().

Member Function Documentation

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

static

Definition at line 98 of file HLTMuonDimuonL2Filter.cc.

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

                                                                                   {
   edm::ParameterSetDescription desc;
   makeHLTFilterDescription(desc);
   desc.add<edm::InputTag>("BeamSpotTag",edm::InputTag("hltOfflineBeamSpot"));
   desc.add<edm::InputTag>("CandTag",edm::InputTag("hltL2MuonCandidates"));
   desc.add<edm::InputTag>("PreviousCandTag",edm::InputTag(""));
   desc.add<edm::InputTag>("SeedMapTag",edm::InputTag("hltL2Muons"));
   desc.add<bool>("FastAccept",false);
   desc.add<double>("MaxEta",2.5);
   desc.add<int>("MinNhits",0);
   desc.add<int>("MinNstations",0);
   desc.add<int>("MinNchambers",2);
   desc.add<double>("MaxDr",100.0);
   desc.add<double>("MaxDz",9999.0);
   desc.add<int>("ChargeOpt",0);
   desc.add<double>("MinPtPair",0.0);
   desc.add<double>("MinPtMax",3.0);
   desc.add<double>("MinPtMin",3.0);
   desc.add<double>("MinInvMass",1.6);
   desc.add<double>("MaxInvMass",5.6);
   desc.add<double>("MinAcop",-1.0);
   desc.add<double>("MaxAcop",3.15);
   desc.add<double>("MinAngle",-999.0);
   desc.add<double>("MaxAngle",2.5);
   desc.add<double>("MinPtBalance",-1.0);
   desc.add<double>("MaxPtBalance",999999.0);
   desc.add<double>("NSigmaPt",0.0);
   descriptions.add("hltMuonDimuonL2Filter", desc);
 }

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

overridevirtual

Implements HLTFilter.

Definition at line 130 of file HLTMuonDimuonL2Filter.cc.

References funct::abs(), accept(), trigger::TriggerFilterObjectWithRefs::addCollectionTag(), trigger::TriggerRefsCollections::addObject(), angle(), ecalDrivenElectronSeedsParameters_cff::beamSpot, beamspotToken_, candTag_, candToken_, chargeOpt_, DEFINE_FWK_MODULE, SoftLeptonByDistance_cfi::distance, fast_Accept_, edm::Ref< C, T, F >::get(), edm::Event::getByToken(), trigger::TriggerRefsCollections::getObjects(), mps_fire::i, HLTMuonL2ToL1Map::isTriggeredByL1(), LogDebug, M_PI, max_Acop_, max_Angle_, max_Dr_, max_Dz_, max_Eta_, max_InvMass_, max_PtBalance_, min_Acop_, min_Angle_, min_InvMass_, min_Nchambers_, min_Nhits_, min_Nstations_, min_PtBalance_, min_PtMax_, min_PtMin_, min_PtPair_, gen::n, nsigma_Pt_, AlCaHLTBitMon_ParallelJobs::p, p1, p2, reco::BeamSpot::position(), previousCandToken_, hiDetachedQuadStep_cff::pt1, hiDetachedQuadStep_cff::pt2, HLTFilter::saveTags(), seedMapToken_, mathSSE::sqrt(), and trigger::TriggerMuon.

 {
 
    double const MuMass = 0.106;
    double const MuMass2 = MuMass*MuMass;
    // 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.
 
    // Ref to Candidate object to be recorded in filter object
    RecoChargedCandidateRef ref1;
    RecoChargedCandidateRef ref2;
 
    // get hold of trks
    Handle<RecoChargedCandidateCollection> mucands;
    iEvent.getByToken(candToken_,mucands);
    if (saveTags()) filterproduct.addCollectionTag(candTag_);
 
    BeamSpot beamSpot;
    Handle<BeamSpot> recoBeamSpotHandle;
    iEvent.getByToken(beamspotToken_,recoBeamSpotHandle);
    beamSpot = *recoBeamSpotHandle;
 
    // get the L2 to L1 map object for this event
    HLTMuonL2ToL1Map mapL2ToL1(previousCandToken_, seedMapToken_, iEvent);
 
    // look at all mucands,  check cuts and add to filter object
    int n = 0;
    double e1,e2;
    Particle::LorentzVector p,p1,p2;
 
    RecoChargedCandidateCollection::const_iterator cand1;
    RecoChargedCandidateCollection::const_iterator cand2;
    for (cand1=mucands->begin(); cand1!=mucands->end(); cand1++) {
       TrackRef tk1 = cand1->get<TrackRef>();
       // eta cut
       LogDebug("HLTMuonDimuonL2Filter") << " 1st muon in loop: q*pt= "
             << tk1->charge()*tk1->pt() << ", eta= " << tk1->eta() << ", hits= " << tk1->numberOfValidHits();
       // find the L1 Particle corresponding to the L2 Track
       if (!mapL2ToL1.isTriggeredByL1(tk1)) continue;
 
       if (fabs(tk1->eta())>max_Eta_) continue;
 
       // cut on number of hits
       if (tk1->numberOfValidHits()<min_Nhits_) continue;
 
       // number of stations
       if (tk1->hitPattern().muonStationsWithAnyHits() < min_Nstations_) continue;
 
       // number of chambers
       if(tk1->hitPattern().dtStationsWithAnyHits() +
      tk1->hitPattern().cscStationsWithAnyHits() < min_Nchambers_) continue;
 
       //dr cut
       //if (fabs(tk1->d0())>max_Dr_) continue;
       if (fabs(tk1->dxy(beamSpot.position()))>max_Dr_) continue;
 
       //dz cut
       if (fabs(tk1->dz())>max_Dz_) continue;
 
       // Pt threshold cut
       double pt1 = tk1->pt();
       double err1 = tk1->error(0);
       double abspar1 = fabs(tk1->parameter(0));
       double ptLx1 = pt1;
       // convert 50% efficiency threshold to 90% efficiency threshold
       if (abspar1>0) ptLx1 += nsigma_Pt_*err1/abspar1*pt1;
       LogDebug("HLTMuonDimuonL2Filter") << " ... 1st muon in loop, pt1= "
             << pt1 << ", ptLx1= " << ptLx1;
 
       cand2 = cand1; cand2++;
       for (; cand2!=mucands->end(); cand2++) {
             TrackRef tk2 = cand2->get<TrackRef>();
 
             // eta cut
             LogDebug("HLTMuonDimuonL2Filter") << " 2nd muon in loop: q*pt= " << tk2->charge()*tk2->pt() << ", eta= " << tk2->eta() << ", hits= " << tk2->numberOfValidHits() << ", d0= " << tk2->d0();
         if (!mapL2ToL1.isTriggeredByL1(tk2)) continue;
 
             if (fabs(tk2->eta())>max_Eta_) continue;
 
             // cut on number of hits
             if (tk2->numberOfValidHits()<min_Nhits_) continue;
 
         // number of stations
         if (tk2->hitPattern().muonStationsWithAnyHits() < min_Nstations_) continue;
 
         // number of chambers
         if(tk2->hitPattern().dtStationsWithAnyHits() +
            tk2->hitPattern().cscStationsWithAnyHits() < min_Nchambers_) continue;
 
             //dr cut
             //if (fabs(tk2->d0())>max_Dr_) continue;
         if (fabs(tk2->dxy(beamSpot.position()))>max_Dr_) continue;
 
             //dz cut
             if (fabs(tk2->dz())>max_Dz_) continue;
 
             // Pt threshold cut
             double pt2 = tk2->pt();
             double err2 = tk2->error(0);
             double abspar2 = fabs(tk2->parameter(0));
             double ptLx2 = pt2;
             // convert 50% efficiency threshold to 90% efficiency threshold
             if (abspar2>0) ptLx2 += nsigma_Pt_*err2/abspar2*pt2;
             LogDebug("HLTMuonDimuonL2Filter") << " ... 2nd muon in loop, pt2= "
                   << pt2 << ", ptLx2= " << ptLx2;
 
             if (ptLx1>ptLx2) {
                   if (ptLx1<min_PtMax_) continue;
                   if (ptLx2<min_PtMin_) continue;
             } else {
                   if (ptLx2<min_PtMax_) continue;
                   if (ptLx1<min_PtMin_) continue;
             }
 
             if (chargeOpt_<0) {
                   if (tk1->charge()*tk2->charge()>0) continue;
             } else if (chargeOpt_>0) {
                   if (tk1->charge()*tk2->charge()<0) continue;
             }
 
             // Acoplanarity
             double acop = fabs(tk1->phi()-tk2->phi());
             if (acop>M_PI) acop = 2*M_PI - acop;
             acop = M_PI - acop;
             LogDebug("HLTMuonDimuonL2Filter") << " ... 1-2 acop= " << acop;
             if (acop<min_Acop_) continue;
             if (acop>max_Acop_) continue;
 
         // 3D angle
         double angle = acos((tk1->px()*tk2->px() + tk1->py()*tk2->py() + tk1->pz()*tk2->pz())/(tk1->p()*tk2->p()));
         LogDebug("HLTMuonDimuonL2Filter") << " ... 1-2 angle= " << angle;
         if (angle < min_Angle_) continue;
         if (angle > max_Angle_) continue;
 
             // Pt balance
             double ptbalance = fabs(tk1->pt()-tk2->pt());
             if (ptbalance<min_PtBalance_) continue;
             if (ptbalance>max_PtBalance_) continue;
 
             // Combined dimuon system
             e1 = sqrt(tk1->momentum().Mag2()+MuMass2);
             e2 = sqrt(tk2->momentum().Mag2()+MuMass2);
             p1 = Particle::LorentzVector(tk1->px(),tk1->py(),tk1->pz(),e1);
             p2 = Particle::LorentzVector(tk2->px(),tk2->py(),tk2->pz(),e2);
             p = p1+p2;
 
             double pt12 = p.pt();
             LogDebug("HLTMuonDimuonL2Filter") << " ... 1-2 pt12= " << pt12;
             if (pt12<min_PtPair_) continue;
 
             double invmass = abs(p.mass());
          // if (invmass>0) invmass = sqrt(invmass); else invmass = 0;
             LogDebug("HLTMuonDimuonL2Filter") << " ... 1-2 invmass= " << invmass;
             if (invmass<min_InvMass_) continue;
             if (invmass>max_InvMass_) continue;
 
             // Add this pair
             n++;
             LogDebug("HLTMuonDimuonL2Filter") << " Track1 passing filter: pt= " << tk1->pt() << ", eta: " << tk1->eta();
             LogDebug("HLTMuonDimuonL2Filter") << " Track2 passing filter: pt= " << tk2->pt() << ", eta: " << tk2->eta();
             LogDebug("HLTMuonDimuonL2Filter") << " Invmass= " << invmass;
 
             bool i1done = false;
             bool i2done = false;
             vector<RecoChargedCandidateRef> vref;
         filterproduct.getObjects(TriggerMuon,vref);
             for (auto & i : vref) {
           RecoChargedCandidateRef candref =  RecoChargedCandidateRef(i);
           TrackRef tktmp = candref->get<TrackRef>();
           if (tktmp==tk1) {
                         i1done = true;
           } else if (tktmp==tk2) {
         i2done = true;
           }
           if (i1done && i2done) break;
             }
     
             if (!i1done) {
           ref1=RecoChargedCandidateRef( Ref<RecoChargedCandidateCollection> (mucands,distance(mucands->begin(), cand1)));
           filterproduct.addObject(TriggerMuon,ref1);
             }
             if (!i2done) {
           ref2=RecoChargedCandidateRef( Ref<RecoChargedCandidateCollection> (mucands,distance(mucands->begin(),cand2 )));
         filterproduct.addObject(TriggerMuon,ref2);
             }
 
             if (fast_Accept_) break;
       }
 
    }
 
    // filter decision
    const bool accept (n >= 1);
 
    LogDebug("HLTMuonDimuonL2Filter") << " >>>>> Result of HLTMuonDimuonL2Filter is "<< accept << ", number of muon pairs passing thresholds= " << n;
 
    return accept;
 }