#include <L2MuonIsolationProducer.h>

Inheritance diagram for L2MuonIsolationProducer:

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

virtual void	produce (edm::Event &, const edm::EventSetup &) override
	Produce isolation maps. More...

virtual	~L2MuonIsolationProducer ()
	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
	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)
	ParameterSet descriptions. 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
bool	optOutputDecision

bool	optOutputIsolatorFloat

muonisolation::MuIsoBaseIsolator *	theDepositIsolator

reco::isodeposit::IsoDepositExtractor *	theExtractor

edm::InputTag	theSACollectionLabel

edm::EDGetTokenT < reco::RecoChargedCandidateCollection >	theSACollectionToken

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

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

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

Author: J. Alcaraz

L2 HLT muon isolation producer

Author: J.Alcaraz

Definition at line 22 of file L2MuonIsolationProducer.h.

Constructor & Destructor Documentation

L2MuonIsolationProducer::L2MuonIsolationProducer ( const edm::ParameterSet & par )

constructor with config

Definition at line 36 of file L2MuonIsolationProducer.cc.

References edm::EDConsumerBase::consumesCollector(), edm::ParameterSet::empty(), reco::get(), edm::ParameterSet::getParameter(), LogDebug, optOutputDecision, optOutputIsolatorFloat, AlCaHLTBitMon_QueryRunRegistry::string, theDepositIsolator, theExtractor, theSACollectionLabel, and theSACollectionToken.

                                                                         :
   theSACollectionLabel(par.getParameter<edm::InputTag>("StandAloneCollectionLabel")),
   theExtractor(0),
   theDepositIsolator(0)
 {
   LogDebug("Muon|RecoMuon|L2MuonIsolationProducer")<<" L2MuonIsolationProducer constructor called";
 
   theSACollectionToken = consumes<RecoChargedCandidateCollection>(theSACollectionLabel);
 
   //
   // Extractor
   //
   edm::ParameterSet extractorPSet = par.getParameter<edm::ParameterSet>("ExtractorPSet");
   std::string extractorName = extractorPSet.getParameter<std::string>("ComponentName");
   theExtractor = IsoDepositExtractorFactory::get()->create( extractorName, extractorPSet, consumesCollector());
 
 
   edm::ParameterSet isolatorPSet = par.getParameter<edm::ParameterSet>("IsolatorPSet");
   bool haveIsolator = !isolatorPSet.empty();
   optOutputDecision = haveIsolator;
   if (optOutputDecision){
     std::string type = isolatorPSet.getParameter<std::string>("ComponentName");
     theDepositIsolator = MuonIsolatorFactory::get()->create(type, isolatorPSet, consumesCollector());
   }
   if (optOutputDecision) produces<edm::ValueMap<bool> >();
   produces<reco::IsoDepositMap>();
 
   optOutputIsolatorFloat = par.getParameter<bool>("WriteIsolatorFloat");
   if (optOutputIsolatorFloat && haveIsolator){
     produces<edm::ValueMap<float> >();
   }
 }

L2MuonIsolationProducer::~L2MuonIsolationProducer ( )

virtual

destructor

Definition at line 70 of file L2MuonIsolationProducer.cc.

References LogDebug, and theExtractor.

                                                  {
   LogDebug("Muon|RecoMuon|L2MuonIsolationProducer")<<" L2MuonIsolationProducer destructor called";
   if (theExtractor) delete theExtractor;
 }

Member Function Documentation

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

static

ParameterSet descriptions.

Definition at line 76 of file L2MuonIsolationProducer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), HLT_25ns14e33_v1_cff::InputTag, AlCaHLTBitMon_QueryRunRegistry::string, and groupFilesInBlocks::temp.

                                                                                          {
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("StandAloneCollectionLabel",edm::InputTag("hltL2MuonCandidates"));
   edm::ParameterSetDescription extractorPSet;
   {
     extractorPSet.add<double>("DR_Veto_H",0.1);
     extractorPSet.add<bool>("Vertex_Constraint_Z",false);
     extractorPSet.add<double>("Threshold_H",0.5);
     extractorPSet.add<std::string>("ComponentName","CaloExtractor");
     extractorPSet.add<double>("Threshold_E",0.2);
     extractorPSet.add<double>("DR_Max",1.0);
     extractorPSet.add<double>("DR_Veto_E",0.07);
     extractorPSet.add<double>("Weight_E",1.5);
     extractorPSet.add<bool>("Vertex_Constraint_XY",false);
     extractorPSet.addUntracked<std::string>("DepositLabel","EcalPlusHcal");
     extractorPSet.add<edm::InputTag>("CaloTowerCollectionLabel",edm::InputTag("towerMaker"));
     extractorPSet.add<double>("Weight_H",1.0);
   }
   desc.add<edm::ParameterSetDescription>("ExtractorPSet",extractorPSet);
   edm::ParameterSetDescription isolatorPSet;
   {
     std::vector<double> temp;
     isolatorPSet.add<std::vector<double> >("ConeSizesRel",std::vector<double>(1, 0.3));
     isolatorPSet.add<double>("EffAreaSFEndcap",1.0);
     isolatorPSet.add<bool>("CutAbsoluteIso",true);
     isolatorPSet.add<bool>("AndOrCuts",true);
     isolatorPSet.add<edm::InputTag>("RhoSrc",edm::InputTag("hltKT6CaloJetsForMuons","rho"));
     isolatorPSet.add<std::vector<double> >("ConeSizes",std::vector<double>(1, 0.3));
     isolatorPSet.add<std::string>("ComponentName","CutsIsolatorWithCorrection");
     isolatorPSet.add<bool>("ReturnRelativeSum",false);
     isolatorPSet.add<double>("RhoScaleBarrel",1.0);
     isolatorPSet.add<double>("EffAreaSFBarrel",1.0);
     isolatorPSet.add<bool>("CutRelativeIso",false);
     isolatorPSet.add<std::vector<double> >("EtaBounds",std::vector<double>(1, 2.411));
     isolatorPSet.add<std::vector<double> >("Thresholds",std::vector<double>(1, 9.9999999E7));
     isolatorPSet.add<bool>("ReturnAbsoluteSum",true);
     isolatorPSet.add<std::vector<double> >("EtaBoundsRel",std::vector<double>(1, 2.411));
     isolatorPSet.add<std::vector<double> >("ThresholdsRel",std::vector<double>(1, 9.9999999E7));
     isolatorPSet.add<double>("RhoScaleEndcap",1.0);
     isolatorPSet.add<double>("RhoMax",9.9999999E7);
     isolatorPSet.add<bool>("UseRhoCorrection",true);
   }
   desc.add<edm::ParameterSetDescription>("IsolatorPSet",isolatorPSet);
   desc.add<bool>("WriteIsolatorFloat",false);
   descriptions.add("hltL2MuonIsolations", desc);
 }

void L2MuonIsolationProducer::produce	(	edm::Event &	event,
		const edm::EventSetup &	eventSetup
	)

overridevirtual

Produce isolation maps.

build deposits

do the business of filling iso map

annoying – I will forget it at some point

Implements edm::stream::EDProducerBase.

Definition at line 124 of file L2MuonIsolationProducer.cc.

References reco::isodeposit::IsoDepositExtractor::deposit(), edm::helper::Filler< Map >::fill(), reco::isodeposit::IsoDepositExtractor::fillVetos(), i, edm::helper::Filler< Map >::insert(), LogDebug, metname, patZpeak::muons, optOutputDecision, optOutputIsolatorFloat, muonisolation::MuIsoBaseIsolator::result(), AlCaHLTBitMon_QueryRunRegistry::string, theDepositIsolator, theExtractor, theSACollectionLabel, theSACollectionToken, muonisolation::MuIsoBaseIsolator::Result::valBool, and muonisolation::MuIsoBaseIsolator::Result::valFloat.

                                                                                {
   std::string metname = "Muon|RecoMuon|L2MuonIsolationProducer";
 
   LogDebug(metname)<<" L2 Muon Isolation producing...";
 
   // Take the SA container
   LogDebug(metname)<<" Taking the StandAlone muons: "<<theSACollectionLabel;
   Handle<RecoChargedCandidateCollection> muons;
   event.getByToken(theSACollectionToken,muons);
 
   // Find deposits and load into event
   LogDebug(metname)<<" Get energy around";
   std::auto_ptr<reco::IsoDepositMap> depMap( new reco::IsoDepositMap());
   std::auto_ptr<edm::ValueMap<bool> > isoMap( new edm::ValueMap<bool> ());
   std::auto_ptr<edm::ValueMap<float> > isoFloatMap( new edm::ValueMap<float> ());
 
   unsigned int nMuons = muons->size();
   std::vector<IsoDeposit> deps(nMuons);
   std::vector<bool> isos(nMuons, false);
   std::vector<float> isoFloats(nMuons, 0);
 
   // fill track collection to use for vetos calculation
   TrackCollection muonTracks;
   for (unsigned int i=0; i<nMuons; i++) {
     TrackRef tk = (*muons)[i].track();
     muonTracks.push_back(*tk);
   }
 
   theExtractor->fillVetos(event,eventSetup,muonTracks);
 
   for (unsigned int i=0; i<nMuons; i++) {
     TrackRef tk = (*muons)[i].track();
 
       deps[i] = theExtractor->deposit(event, eventSetup, *tk);
 
       if (optOutputDecision){
     muonisolation::MuIsoBaseIsolator::DepositContainer isoContainer(1,muonisolation::MuIsoBaseIsolator::DepositAndVetos(&deps[i]));
     muonisolation::MuIsoBaseIsolator::Result isoResult = theDepositIsolator->result( isoContainer, *tk, &event );
     isos[i] = isoResult.valBool;
     isoFloats[i] = isoResult.valFloat;
       }
   }
 
 
 
   reco::IsoDepositMap::Filler depFiller(*depMap);
   depFiller.insert(muons, deps.begin(), deps.end());
   depFiller.fill();
   event.put(depMap);
 
   if (optOutputDecision){
     edm::ValueMap<bool> ::Filler isoFiller(*isoMap);
     isoFiller.insert(muons, isos.begin(), isos.end());
     isoFiller.fill();
     event.put(isoMap);
 
     if (optOutputIsolatorFloat){
       edm::ValueMap<float> ::Filler isoFloatFiller(*isoFloatMap);
       isoFloatFiller.insert(muons, isoFloats.begin(), isoFloats.end());
       isoFloatFiller.fill();
       event.put(isoFloatMap);
     }
   }
 
   LogDebug(metname) <<" Event loaded"
             <<"================================";
 }