#include <L2MuonIsolationProducer.h>

Inheritance diagram for L2MuonIsolationProducer:

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

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

	~L2MuonIsolationProducer () override
	destructor More...

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

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

Private Attributes
bool	optOutputDecision

bool	optOutputIsolatorFloat

std::unique_ptr< muonisolation::MuIsoBaseIsolator >	theDepositIsolator

std::unique_ptr< 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

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::L2MuonIsolationProducer ( const edm::ParameterSet & par )

constructor with config

Definition at line 36 of file L2MuonIsolationProducer.cc.

     : theSACollectionLabel(par.getParameter<edm::InputTag>("StandAloneCollectionLabel")) {
   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 = std::unique_ptr<reco::isodeposit::IsoDepositExtractor>{
       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 = std::unique_ptr<muonisolation::MuIsoBaseIsolator>{
         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>>();
   }
 }

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

◆ ~L2MuonIsolationProducer()

L2MuonIsolationProducer::~L2MuonIsolationProducer ( )

override

destructor

Definition at line 69 of file L2MuonIsolationProducer.cc.

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

References LogDebug.

Member Function Documentation

◆ fillDescriptions()

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

static

ParameterSet descriptions.

Definition at line 74 of file L2MuonIsolationProducer.cc.

                                                                                          {
   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);
 }

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

◆ produce()

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

override

Produce isolation maps.

build deposits

do the business of filling iso map

annoying – I will forget it at some point

Definition at line 122 of file L2MuonIsolationProducer.cc.

                                                                                 {
   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";
   auto depMap = std::make_unique<reco::IsoDepositMap>();
   auto isoMap = std::make_unique<edm::ValueMap<bool>>();
   auto isoFloatMap = std::make_unique<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(std::move(depMap));
  
   if (optOutputDecision) {
     edm::ValueMap<bool>::Filler isoFiller(*isoMap);
     isoFiller.insert(muons, isos.begin(), isos.end());
     isoFiller.fill();  
     event.put(std::move(isoMap));
  
     if (optOutputIsolatorFloat) {
       edm::ValueMap<float>::Filler isoFloatFiller(*isoFloatMap);
       isoFloatFiller.insert(muons, isoFloats.begin(), isoFloats.end());
       isoFloatFiller.fill();  
       event.put(std::move(isoFloatMap));
     }
   }
  
   LogDebug(metname) << " Event loaded"
                     << "================================";
 }