#include <HLTMuonIsoFilter.h>

Inheritance diagram for HLTMuonIsoFilter:

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

	HLTMuonIsoFilter (const edm::ParameterSet &)

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

	~HLTMuonIsoFilter ()

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

int	module () const

const std::string *	moduleLabel () const

int	path () const

const std::string *	pathName () const

std::pair< int, int >	pmid () const

bool	saveTags () const

virtual	~HLTFilter ()

Public Member Functions inherited from edm::EDFilter
	EDFilter ()

virtual	~EDFilter ()

Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

virtual	~ProducerBase ()

Private Attributes
edm::InputTag	candTag_

std::vector< edm::InputTag >	depTag_

int	min_N_

edm::InputTag	previousCandTag_

const muonisolation::MuIsoBaseIsolator *	theDepositIsolator

Additional Inherited Members
Public Types inherited from edm::EDFilter
typedef EDFilter	ModuleType

typedef WorkerT< EDFilter >	WorkerType

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

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

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

Protected Member Functions inherited from edm::EDFilter
CurrentProcessingContext const *	currentContext () const

Protected Member Functions inherited from edm::ProducerBase
template<class TProducer , class TMethod >
void	callWhenNewProductsRegistered (TProducer *iProd, TMethod iMethod)

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 18 of file HLTMuonIsoFilter.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 29 of file HLTMuonIsoFilter.cc.

References candTag_, depTag_, edm::ParameterSet::empty(), edm::InputTag::encode(), alcaDQMUpload::encode(), reco::get(), edm::ParameterSet::getParameter(), i, LogDebug, min_N_, o2o::tags, and theDepositIsolator.

                                                                  : HLTFilter(iConfig),
    candTag_ (iConfig.getParameter< edm::InputTag > ("CandTag") ),
    previousCandTag_ (iConfig.getParameter<edm::InputTag > ("PreviousCandTag")),
    depTag_  (iConfig.getParameter< std::vector< edm::InputTag > >("DepTag" ) ),
    theDepositIsolator(0),
    min_N_   (iConfig.getParameter<int> ("MinN"))
 {
   std::stringstream tags;
   for (unsigned int i=0;i!=depTag_.size();++i)
     tags<<" IsoTag["<<i<<"] : "<<depTag_[i].encode()<<" \n";
    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);
    }
    
    if (theDepositIsolator) produces<edm::ValueMap<bool> >();
 }

HLTMuonIsoFilter::~HLTMuonIsoFilter ( )

Definition at line 54 of file HLTMuonIsoFilter.cc.

55 {

56 }

Member Function Documentation

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

virtual

Implements HLTFilter.

Definition at line 64 of file HLTMuonIsoFilter.cc.

References accept(), trigger::TriggerFilterObjectWithRefs::addCollectionTag(), trigger::TriggerRefsCollections::addObject(), candTag_, depTag_, edm::helper::Filler< Map >::fill(), edm::Event::get(), edm::Ref< C, T, F >::get(), edm::Event::getByLabel(), i, edm::Ref< C, T, F >::id(), edm::helper::Filler< Map >::insert(), edm::Ref< C, T, F >::isNonnull(), LogDebug, min_N_, previousCandTag_, 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.getByLabel (candTag_,mucands);
    Handle<TriggerFilterObjectWithRefs> previousLevelCands;
    iEvent.getByLabel (previousCandTag_,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.getByLabel (depTag_[i],depMap[i]);
    }else{
      bool success = iEvent.getByLabel(depTag_.front(), 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]))[tk];
      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)[tk];
      }
      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);     
        // get a track ref
 
        TrackRef aRef = mucands->front().get<TrackRef>();
        // get the corresponding handle
        edm::Handle<reco::TrackCollection> HandleToTrackRef;
        iEvent.get(aRef.id(), HandleToTrackRef);
        isoFiller.insert(HandleToTrackRef, 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
	)

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