#include <HLTMuonIsoFilter.h>

Inheritance diagram for HLTMuonIsoFilter:

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

	HLTMuonIsoFilter (const edm::ParameterSet &)

	~HLTMuonIsoFilter ()

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

virtual	~HLTFilter ()

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

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

ModuleDescription const &	moduleDescription () const

virtual	~EDFilterBase ()

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

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

virtual	~ProducerBase ()

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 (const std::string &iProcessName, std::vector< const char * > &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 ()

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)

Static Private Member Functions
static bool	triggerdByPreviousLevel (const reco::RecoChargedCandidateRef &, const std::vector< reco::RecoChargedCandidateRef > &)

Private Attributes
edm::InputTag	candTag_

edm::EDGetTokenT < reco::RecoChargedCandidateCollection >	candToken_

edm::EDGetTokenT < edm::ValueMap< bool > >	decMapToken_

std::vector< edm::InputTag >	depTag_

std::vector< edm::EDGetTokenT < edm::ValueMap < reco::IsoDeposit > > >	depToken_

int	min_N_

edm::InputTag	previousCandTag_

edm::EDGetTokenT < trigger::TriggerFilterObjectWithRefs >	previousCandToken_

const muonisolation::MuIsoBaseIsolator *	theDepositIsolator

Additional Inherited Members
Public Types inherited from edm::global::EDFilterBase
typedef EDFilterBase	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

This class is an HLTFilter (-> EDFilter) implementing the isolation filtering for HLT muons

Author: J. Alcaraz

See header file for documentation

Author: J. Alcaraz

Definition at line 24 of file HLTMuonIsoFilter.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 33 of file HLTMuonIsoFilter.cc.

References candTag_, edm::EDConsumerBase::consumes(), edm::EDConsumerBase::consumesCollector(), decMapToken_, depTag_, depToken_, edm::ParameterSet::empty(), edm::InputTag::encode(), reco::get(), edm::ParameterSet::getParameter(), i, LogDebug, min_N_, AlCaHLTBitMon_QueryRunRegistry::string, o2o::tags, and theDepositIsolator.

                                                                  : HLTFilter(iConfig),
    candTag_ (iConfig.getParameter< edm::InputTag > ("CandTag") ),
    candToken_(consumes<reco::RecoChargedCandidateCollection>(candTag_)),
    previousCandTag_ (iConfig.getParameter<edm::InputTag > ("PreviousCandTag")),
    previousCandToken_ (consumes<trigger::TriggerFilterObjectWithRefs>(previousCandTag_)),
    depTag_  (iConfig.getParameter< std::vector< edm::InputTag > >("DepTag" ) ),
    depToken_(0),
    theDepositIsolator(0),
    min_N_   (iConfig.getParameter<int> ("MinN"))
 {
   std::stringstream tags;
   for (unsigned int i=0;i!=depTag_.size();++i) {
     depToken_.push_back(consumes<edm::ValueMap<reco::IsoDeposit> >(depTag_[i]));
     tags<<" IsoTag["<<i<<"] : "<<depTag_[i].encode()<<" \n";
   }
   decMapToken_ = consumes<edm::ValueMap<bool> >(depTag_.front());
 
    LogDebug("HLTMuonIsoFilter") << " candTag : " << candTag_.encode()
                 << "\n" << tags
                 << "  MinN : " << min_N_;
 
    edm::ParameterSet isolatorPSet = iConfig.getParameter<edm::ParameterSet>("IsolatorPSet");
    if (isolatorPSet.empty()) {
      theDepositIsolator=0;
        }else{
      std::string type = isolatorPSet.getParameter<std::string>("ComponentName");
      theDepositIsolator = MuonIsolatorFactory::get()->create(type, isolatorPSet, consumesCollector());
    }
 
    if (theDepositIsolator) produces<edm::ValueMap<bool> >();
 }

HLTMuonIsoFilter::~HLTMuonIsoFilter ( )

Definition at line 65 of file HLTMuonIsoFilter.cc.

66 {

67 }

Member Function Documentation

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

static

Definition at line 73 of file HLTMuonIsoFilter.cc.

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

                                                                              {
   edm::ParameterSetDescription desc;
   makeHLTFilterDescription(desc);
   desc.add<edm::InputTag>("CandTag",edm::InputTag("hltL3MuonCandidates"));
   desc.add<edm::InputTag>("PreviousCandTag",edm::InputTag(""));
   desc.add<int>("MinN",1);
   std::vector<edm::InputTag> depTag(1,edm::InputTag("hltL3MuonIsolations"));
   desc.add<std::vector<edm::InputTag> >("DepTag",depTag);
   edm::ParameterSetDescription isolatorPSet;
   desc.add<edm::ParameterSetDescription>("IsolatorPSet",isolatorPSet);
   descriptions.add("hltMuonIsoFilter", desc);
 }

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

overridevirtual

Implements HLTFilter.

Definition at line 88 of file HLTMuonIsoFilter.cc.

References accept(), trigger::TriggerFilterObjectWithRefs::addCollectionTag(), trigger::TriggerRefsCollections::addObject(), candTag_, candToken_, decMapToken_, depTag_, depToken_, edm::helper::Filler< Map >::fill(), edm::Ref< C, T, F >::get(), edm::Event::getByToken(), i, edm::helper::Filler< Map >::insert(), edm::Ref< C, T, F >::isNonnull(), LogDebug, min_N_, previousCandToken_, reco::IsoDeposit::print(), edm::Event::put(), dt_dqm_sourceclient_common_cff::reco, muonisolation::MuIsoBaseIsolator::result(), HLTFilter::saveTags(), corrVsCorr::selection, summarizeEdmComparisonLogfiles::success, theDepositIsolator, triggerdByPreviousLevel(), trigger::TriggerMuon, and muonisolation::MuIsoBaseIsolator::Result::valBool.

 {
    using namespace std;
    using namespace edm;
    using namespace trigger;
    using namespace reco;
 
    // 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.
 
    //the decision map
    std::auto_ptr<edm::ValueMap<bool> >
      isoMap( new edm::ValueMap<bool> ());
 
    // get hold of trks
    Handle<RecoChargedCandidateCollection> mucands;
    if (saveTags()) filterproduct.addCollectionTag(candTag_);
    iEvent.getByToken(candToken_,mucands);
    Handle<TriggerFilterObjectWithRefs> previousLevelCands;
    iEvent.getByToken(previousCandToken_,previousLevelCands);
    vector<RecoChargedCandidateRef> vcands;
    previousLevelCands->getObjects(TriggerMuon,vcands);
 
    //get hold of energy deposition
    unsigned int nDep=depTag_.size();
    std::vector< Handle<edm::ValueMap<reco::IsoDeposit> > > depMap(nDep);
    Handle<edm::ValueMap<bool> > decisionMap;
    muonisolation::MuIsoBaseIsolator::DepositContainer isoContainer(nDep);
 
    if (theDepositIsolator){
      for (unsigned int i=0;i!=nDep;++i) iEvent.getByToken(depToken_[i],depMap[i]);
    }else{
      bool success = iEvent.getByToken(decMapToken_, decisionMap);
      LogDebug("HLTMuonIsoFilter")<<"get decisionMap " << success;
    }
 
    // look at all mucands,  check cuts and add to filter object
    int nIsolatedMu = 0;
    unsigned int nMu=mucands->size();
    std::vector<bool> isos(nMu, false);
    unsigned int iMu=0;
    for (; iMu<nMu; iMu++) {
      RecoChargedCandidateRef candref(mucands,iMu);
      LogDebug("HLTMuonIsoFilter") << "candref isNonnull " << candref.isNonnull();
 
      //did this candidate triggered at previous stage.
      if (!triggerdByPreviousLevel(candref,vcands)) continue;
 
      //reference to the track
      TrackRef tk = candref->get<TrackRef>();
      LogDebug("HLTMuonIsoFilter") << "tk isNonNull " << tk.isNonnull();
      if (theDepositIsolator){
 
        //get the deposits
        for(unsigned int iDep=0;iDep!=nDep;++iDep){
 
      const edm::ValueMap<reco::IsoDeposit> ::value_type & muonDeposit = (*(depMap[iDep]))[candref];
      LogDebug("HLTMuonIsoFilter") << " Muon with q*pt= " << tk->charge()*tk->pt() << " (" << candref->charge()*candref->pt() << ") " << ", eta= " << tk->eta() << " (" << candref->eta() << ") " << "; has deposit["<<iDep<<"]: " << muonDeposit.print();
      isoContainer[iDep] = muonisolation::MuIsoBaseIsolator::DepositAndVetos(&muonDeposit);
 
        }
 
        //get the selection
        muonisolation::MuIsoBaseIsolator::Result selection = theDepositIsolator->result( isoContainer, *tk );
        isos[iMu]=selection.valBool;
 
      }else{
        //get the decision from the event
        isos[iMu]=(*decisionMap)[candref];
      }
      LogDebug("HLTMuonIsoFilter") << " Muon with q*pt= " << tk->charge()*tk->pt() << ", eta= " << tk->eta() << "; "<<(isos[iMu]?"Is an isolated muon.":"Is NOT an isolated muon.");
 
      if (!isos[iMu]) continue;
 
      nIsolatedMu++;
      filterproduct.addObject(TriggerMuon,candref);
    }
 
    // filter decision
    const bool accept (nIsolatedMu >= min_N_);
 
    if (theDepositIsolator){
      //put the decision map
      if (nMu!=0){
        edm::ValueMap<bool> ::Filler isoFiller(*isoMap);
        isoFiller.insert(mucands, isos.begin(), isos.end());
        isoFiller.fill();
      }
      iEvent.put(isoMap);
    }
 
    LogDebug("HLTMuonIsoFilter") << " >>>>> Result of HLTMuonIsoFilter is " << accept << ", number of muons passing isolation cuts= " << nIsolatedMu;
 
    return accept;
 }

bool HLTMuonIsoFilter::triggerdByPreviousLevel	(	const reco::RecoChargedCandidateRef &	candref,
		const std::vector< reco::RecoChargedCandidateRef > &	vcands
	)

staticprivate

Definition at line 185 of file HLTMuonIsoFilter.cc.

References i, and convertSQLiteXML::ok.

Referenced by hltFilter().

                                                                                                                                                {
   bool ok=false;
   unsigned int i=0;
   unsigned int i_max=vcands.size();
   for (;i!=i_max;++i){
     if (candref == vcands[i]) { ok=true; break;}
   }
 
   return ok;
 }