#include <L3MuonCombinedRelativeIsolationProducer.h>

Inheritance diagram for L3MuonCombinedRelativeIsolationProducer:

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

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

	~L3MuonCombinedRelativeIsolationProducer () override
	destructor More...

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

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

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

Private Attributes
std::unique_ptr< reco::isodeposit::IsoDepositExtractor >	caloExtractor

bool	optOutputIsoDeposits

bool	printDebug

bool	theApplyCutsORmaxNTracks
	apply or not the maxN cut on top of the sumPt (or nominall eff) < cuts More...

edm::InputTag	theCaloDepsLabel

edm::EDGetTokenT< edm::ValueMap< float > >	theCaloDepsToken

edm::ParameterSet	theConfig

muonisolation::Cuts	theCuts

int	theMaxNTracks

edm::InputTag	theMuonCollectionLabel

edm::EDGetTokenT< reco::RecoChargedCandidateCollection >	theMuonCollectionToken

double	theTrackPt_Min

std::unique_ptr< reco::isodeposit::IsoDepositExtractor >	trkExtractor

bool	useCaloIso
	flag to include or exclude calo iso from calculation More...

bool	useRhoCorrectedCaloDeps

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
using	CacheTypes = CacheContexts< T... >

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T... >

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

Definition at line 26 of file L3MuonCombinedRelativeIsolationProducer.h.

Constructor & Destructor Documentation

◆ L3MuonCombinedRelativeIsolationProducer()

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

constructor with config

test cutsName only. The depositType is informational only (has not been used so far) [VK]

Definition at line 38 of file L3MuonCombinedRelativeIsolationProducer.cc.

References caloExtractor, L3NominalEfficiencyConfigurator::cuts(), get, edm::ParameterSet::getParameter(), LogDebug, LogTrace, optOutputIsoDeposits, muonisolation::Cuts::print(), AlCaHLTBitMon_QueryRunRegistry::string, theApplyCutsORmaxNTracks, theCaloDepsLabel, theCaloDepsToken, theConfig, theCuts, theMaxNTracks, theMuonCollectionLabel, theMuonCollectionToken, theTrackPt_Min, trkExtractor, useCaloIso, and useRhoCorrectedCaloDeps.

     : theConfig(par),
       theMuonCollectionLabel(par.getParameter<InputTag>("inputMuonCollection")),
       optOutputIsoDeposits(par.getParameter<bool>("OutputMuIsoDeposits")),
       useCaloIso(par.existsAs<bool>("UseCaloIso") ? par.getParameter<bool>("UseCaloIso") : true),
       useRhoCorrectedCaloDeps(par.existsAs<bool>("UseRhoCorrectedCaloDeposits")
                                   ? par.getParameter<bool>("UseRhoCorrectedCaloDeposits")
                                   : false),
       theCaloDepsLabel(par.existsAs<InputTag>("CaloDepositsLabel") ? par.getParameter<InputTag>("CaloDepositsLabel")
                                                                    : InputTag("hltL3CaloMuonCorrectedIsolations")),
       theTrackPt_Min(-1),
       printDebug(par.getParameter<bool>("printDebug")) {
   LogDebug("RecoMuon|L3MuonCombinedRelativeIsolationProducer") << " L3MuonCombinedRelativeIsolationProducer CTOR";
 
   theMuonCollectionToken = consumes<RecoChargedCandidateCollection>(theMuonCollectionLabel);
   theCaloDepsToken = consumes<edm::ValueMap<float>>(theCaloDepsLabel);
 
   if (optOutputIsoDeposits) {
     produces<reco::IsoDepositMap>("trkIsoDeposits");
     if (useRhoCorrectedCaloDeps == false)  // otherwise, calo deposits have been previously computed
       produces<reco::IsoDepositMap>("caloIsoDeposits");
     //produces<std::vector<double> >("combinedRelativeIsoDeposits");
     produces<edm::ValueMap<double>>("combinedRelativeIsoDeposits");
   }
   produces<edm::ValueMap<bool>>();
 
   //
   // Extractor
   //
   // Calorimeters (ONLY if not previously computed)
   //
   if (useCaloIso && (useRhoCorrectedCaloDeps == false)) {
     edm::ParameterSet caloExtractorPSet = theConfig.getParameter<edm::ParameterSet>("CaloExtractorPSet");
 
     std::string caloExtractorName = caloExtractorPSet.getParameter<std::string>("ComponentName");
     caloExtractor = std::unique_ptr<reco::isodeposit::IsoDepositExtractor>{
         IsoDepositExtractorFactory::get()->create(caloExtractorName, caloExtractorPSet, consumesCollector())};
     //std::string caloDepositType = caloExtractorPSet.getUntrackedParameter<std::string>("DepositLabel"); // N.B. Not used in the following!
   }
 
   // Tracker
   //
   edm::ParameterSet trkExtractorPSet = theConfig.getParameter<edm::ParameterSet>("TrkExtractorPSet");
 
   theTrackPt_Min = theConfig.getParameter<double>("TrackPt_Min");
   std::string trkExtractorName = trkExtractorPSet.getParameter<std::string>("ComponentName");
   trkExtractor = std::unique_ptr<reco::isodeposit::IsoDepositExtractor>{
       IsoDepositExtractorFactory::get()->create(trkExtractorName, trkExtractorPSet, consumesCollector())};
   //std::string trkDepositType = trkExtractorPSet.getUntrackedParameter<std::string>("DepositLabel"); // N.B. Not used in the following!
 
   //
   // Cuts for track isolation
   //
   edm::ParameterSet cutsPSet = theConfig.getParameter<edm::ParameterSet>("CutsPSet");
   std::string cutsName = cutsPSet.getParameter<std::string>("ComponentName");
   if (cutsName == "SimpleCuts") {
     theCuts = Cuts(cutsPSet);
   } else if (
       //        (cutsName== "L3NominalEfficiencyCuts_PXLS" && depositType=="PXLS")
       //     || (cutsName== "L3NominalEfficiencyCuts_TRKS" && depositType=="TRKS")
       (cutsName == "L3NominalEfficiencyCuts_PXLS") || (cutsName == "L3NominalEfficiencyCuts_TRKS")) {
     theCuts = L3NominalEfficiencyConfigurator(cutsPSet).cuts();
   } else {
     LogError("L3MuonCombinedRelativeIsolationProducer::beginJob") << "cutsName: " << cutsPSet << " is not recognized:"
                                                                   << " theCuts not set!";
   }
   LogTrace("") << theCuts.print();
 
   // (kludge) additional cut on the number of tracks
   theMaxNTracks = cutsPSet.getParameter<int>("maxNTracks");
   theApplyCutsORmaxNTracks = cutsPSet.getParameter<bool>("applyCutsORmaxNTracks");
 }

◆ ~L3MuonCombinedRelativeIsolationProducer()

L3MuonCombinedRelativeIsolationProducer::~L3MuonCombinedRelativeIsolationProducer ( )

override

destructor

Definition at line 113 of file L3MuonCombinedRelativeIsolationProducer.cc.

References LogDebug.

                                                                                   {
   LogDebug("RecoMuon|L3MuonCombinedRelativeIsolationProducer") << " L3MuonCombinedRelativeIsolationProducer DTOR";
 }

Member Function Documentation

◆ fillDescriptions()

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

static

ParameterSet descriptions.

Definition at line 118 of file L3MuonCombinedRelativeIsolationProducer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), submitPVResolutionJobs::desc, HLT_2022v11_cff::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                                          {
   edm::ParameterSetDescription desc;
   desc.add<bool>("UseRhoCorrectedCaloDeposits", false);
   desc.add<bool>("UseCaloIso", true);
   desc.add<edm::InputTag>("CaloDepositsLabel", edm::InputTag("hltL3CaloMuonCorrectedIsolations"));
   desc.add<edm::InputTag>("inputMuonCollection", edm::InputTag("hltL3MuonCandidates"));
   desc.add<bool>("OutputMuIsoDeposits", true);
   desc.add<double>("TrackPt_Min", -1.0);
   desc.add<bool>("printDebug", false);
   edm::ParameterSetDescription cutsPSet;
   {
     cutsPSet.add<std::vector<double>>("ConeSizes", std::vector<double>(1, 0.24));
     cutsPSet.add<std::string>("ComponentName", "SimpleCuts");
     cutsPSet.add<std::vector<double>>("Thresholds", std::vector<double>(1, 0.1));
     cutsPSet.add<int>("maxNTracks", -1);
     cutsPSet.add<std::vector<double>>("EtaBounds", std::vector<double>(1, 2.411));
     cutsPSet.add<bool>("applyCutsORmaxNTracks", false);
   }
   desc.add<edm::ParameterSetDescription>("CutsPSet", cutsPSet);
   edm::ParameterSetDescription trkExtractorPSet;
   {
     trkExtractorPSet.add<double>("Chi2Prob_Min", -1.0);
     trkExtractorPSet.add<double>("Chi2Ndof_Max", 1.0E64);
     trkExtractorPSet.add<double>("Diff_z", 0.2);
     trkExtractorPSet.add<edm::InputTag>("inputTrackCollection", edm::InputTag("hltPixelTracks"));
     trkExtractorPSet.add<double>("ReferenceRadius", 6.0);
     trkExtractorPSet.add<edm::InputTag>("BeamSpotLabel", edm::InputTag("hltOnlineBeamSpot"));
     trkExtractorPSet.add<std::string>("ComponentName", "PixelTrackExtractor");
     trkExtractorPSet.add<double>("DR_Max", 0.24);
     trkExtractorPSet.add<double>("Diff_r", 0.1);
     trkExtractorPSet.add<bool>("VetoLeadingTrack", true);
     trkExtractorPSet.add<double>("DR_VetoPt", 0.025);
     trkExtractorPSet.add<double>("DR_Veto", 0.01);
     trkExtractorPSet.add<unsigned int>("NHits_Min", 0);
     trkExtractorPSet.add<double>("Pt_Min", -1.0);
     trkExtractorPSet.addUntracked<std::string>("DepositLabel", "PXLS");
     trkExtractorPSet.add<std::string>("BeamlineOption", "BeamSpotFromEvent");
     trkExtractorPSet.add<bool>("PropagateTracksToRadius", true);
     trkExtractorPSet.add<double>("PtVeto_Min", 2.0);
   }
   desc.add<edm::ParameterSetDescription>("TrkExtractorPSet", trkExtractorPSet);
   edm::ParameterSetDescription caloExtractorPSet;
   {
     caloExtractorPSet.add<double>("DR_Veto_H", 0.1);
     caloExtractorPSet.add<bool>("Vertex_Constraint_Z", false);
     caloExtractorPSet.add<double>("Threshold_H", 0.5);
     caloExtractorPSet.add<std::string>("ComponentName", "CaloExtractor");
     caloExtractorPSet.add<double>("Threshold_E", 0.2);
     caloExtractorPSet.add<double>("DR_Max", 0.24);
     caloExtractorPSet.add<double>("DR_Veto_E", 0.07);
     caloExtractorPSet.add<double>("Weight_E", 1.5);
     caloExtractorPSet.add<bool>("Vertex_Constraint_XY", false);
     caloExtractorPSet.addUntracked<std::string>("DepositLabel", "EcalPlusHcal");
     caloExtractorPSet.add<edm::InputTag>("CaloTowerCollectionLabel", edm::InputTag("hltTowerMakerForMuons"));
     caloExtractorPSet.add<double>("Weight_H", 1.0);
   }
   desc.add<edm::ParameterSetDescription>("CaloExtractorPSet", caloExtractorPSet);
   descriptions.add("hltL3MuonIsolations", desc);
 }

◆ produce()

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

override

Produce isolation maps.

Definition at line 178 of file L3MuonCombinedRelativeIsolationProducer.cc.

References caloExtractor, submitPVResolutionJobs::count, gather_cfg::cout, PA_MinBiasSkim_cff::cut, reco::IsoDeposit::depositAndCountWithin(), reco::IsoDeposit::depositWithin(), options_cfi::eventSetup, edm::helper::Filler< Map >::fill(), mps_fire::i, edm::helper::Filler< Map >::insert(), metname, eostools::move(), amptDefaultParameters_cff::mu, PDWG_BPHSkim_cff::muons, custom_jme_cff::nMuons, optOutputIsoDeposits, reco::IsoDeposit::print(), printDebug, mps_fire::result, AlCaHLTBitMon_QueryRunRegistry::string, theApplyCutsORmaxNTracks, theCaloDepsToken, theCuts, theMaxNTracks, theMuonCollectionLabel, theMuonCollectionToken, theTrackPt_Min, trkExtractor, useCaloIso, and useRhoCorrectedCaloDeps.

                                                                                                 {
   std::string metname = "RecoMuon|L3MuonCombinedRelativeIsolationProducer";
 
   if (printDebug)
     std::cout << " L3 Muon Isolation producing..."
               << " BEGINING OF EVENT "
               << "================================" << std::endl;
 
   // Take the SA container
   if (printDebug)
     std::cout << " Taking the muons: " << theMuonCollectionLabel << std::endl;
   Handle<RecoChargedCandidateCollection> muons;
   event.getByToken(theMuonCollectionToken, muons);
 
   // Take calo deposits with rho corrections (ONLY if previously computed)
   Handle<edm::ValueMap<float>> caloDepWithCorrMap;
   if (useRhoCorrectedCaloDeps)
     event.getByToken(theCaloDepsToken, caloDepWithCorrMap);
 
   auto caloDepMap = std::make_unique<reco::IsoDepositMap>();
   auto trkDepMap = std::make_unique<reco::IsoDepositMap>();
 
   auto comboIsoDepMap = std::make_unique<edm::ValueMap<bool>>();
 
   //auto combinedRelativeDeps = std::make_unique<std::vector<double>>();
   auto combinedRelativeDepMap = std::make_unique<edm::ValueMap<double>>();
 
   //
   // get Vetos and deposits
   //
   unsigned int nMuons = muons->size();
 
   IsoDeposit::Vetos trkVetos(nMuons);
   std::vector<IsoDeposit> trkDeps(nMuons);
 
   // IsoDeposit::Vetos caloVetos(nMuons);
   // std::vector<IsoDeposit> caloDeps(nMuons);
   // std::vector<float> caloCorrDeps(nMuons, 0.);  // if calo deposits with corrections available
 
   IsoDeposit::Vetos caloVetos;
   std::vector<IsoDeposit> caloDeps;
   std::vector<float> caloCorrDeps;  // if calo deposits with corrections available
 
   if (useCaloIso && useRhoCorrectedCaloDeps) {
     caloCorrDeps.resize(nMuons, 0.);
   } else if (useCaloIso) {
     caloVetos.resize(nMuons);
     caloDeps.resize(nMuons);
   }
 
   std::vector<double> combinedRelativeDeps(nMuons, 0.);
   std::vector<bool> combinedRelativeIsos(nMuons, false);
 
   for (unsigned int i = 0; i < nMuons; i++) {
     RecoChargedCandidateRef candref(muons, i);
     TrackRef mu = candref->track();
 
     trkDeps[i] = trkExtractor->deposit(event, eventSetup, *mu);
     trkVetos[i] = trkDeps[i].veto();
 
     if (useCaloIso && useRhoCorrectedCaloDeps) {
       caloCorrDeps[i] = (*caloDepWithCorrMap)[candref];
     } else if (useCaloIso) {
       caloDeps[i] = caloExtractor->deposit(event, eventSetup, *mu);
       caloVetos[i] = caloDeps[i].veto();
     }
   }
 
   //
   // add here additional vetos
   //
   //.....
 
   //
   // actual cut step
   //
 
   if (printDebug)
     std::cout << "Looping over deposits...." << std::endl;
   for (unsigned int iMu = 0; iMu < nMuons; ++iMu) {
     if (printDebug)
       std::cout << "Muon number = " << iMu << std::endl;
     TrackRef mu = (*muons)[iMu].track();
 
     // cuts
     const Cuts::CutSpec& cut = theCuts(mu->eta());
 
     if (printDebug)
       std::cout << "CUTDEBUG: Muon eta = " << mu->eta() << std::endl
                 << "CUTDEBUG: Muon pt  = " << mu->pt() << std::endl
                 << "CUTDEBUG: minEta   = " << cut.etaRange.min() << std::endl
                 << "CUTDEBUG: maxEta   = " << cut.etaRange.max() << std::endl
                 << "CUTDEBUG: consize  = " << cut.conesize << std::endl
                 << "CUTDEBUG: thresho  = " << cut.threshold << std::endl;
 
     const IsoDeposit& trkDeposit = trkDeps[iMu];
     if (printDebug)
       std::cout << trkDeposit.print();
     std::pair<double, int> trkIsoSumAndCount = trkDeposit.depositAndCountWithin(cut.conesize, trkVetos, theTrackPt_Min);
 
     double caloIsoSum = 0.;
     if (useCaloIso && useRhoCorrectedCaloDeps) {
       caloIsoSum = caloCorrDeps[iMu];
       if (caloIsoSum < 0.)
         caloIsoSum = 0.;
       if (printDebug)
         std::cout << "Rho-corrected calo deposit (min. 0) = " << caloIsoSum << std::endl;
     } else if (useCaloIso) {
       const IsoDeposit& caloDeposit = caloDeps[iMu];
       if (printDebug)
         std::cout << caloDeposit.print();
       caloIsoSum = caloDeposit.depositWithin(cut.conesize, caloVetos);
     }
 
     double trkIsoSum = trkIsoSumAndCount.first;
     int count = trkIsoSumAndCount.second;
 
     double muPt = mu->pt();
     if (muPt < 1.)
       muPt = 1.;
     double combinedRelativeDeposit = ((trkIsoSum + caloIsoSum) / muPt);
     bool result = (combinedRelativeDeposit < cut.threshold);
     if (theApplyCutsORmaxNTracks)
       result |= count <= theMaxNTracks;
     if (printDebug)
       std::cout << "  trk dep in cone:  " << trkIsoSum << "  with count " << count << std::endl
                 << "  calo dep in cone: " << caloIsoSum << std::endl
                 << "  muPt: " << muPt << std::endl
                 << "  relIso:  " << combinedRelativeDeposit << std::endl
                 << "  is isolated: " << result << std::endl;
 
     combinedRelativeIsos[iMu] = result;
     //combinedRelativeDeps->push_back(combinedRelativeDeposit);
     combinedRelativeDeps[iMu] = combinedRelativeDeposit;
   }
 
   //
   // store
   //
   if (optOutputIsoDeposits) {
     reco::IsoDepositMap::Filler depFillerTrk(*trkDepMap);
     depFillerTrk.insert(muons, trkDeps.begin(), trkDeps.end());
     depFillerTrk.fill();
     event.put(std::move(trkDepMap), "trkIsoDeposits");
 
     if (useCaloIso && (useRhoCorrectedCaloDeps == false)) {
       reco::IsoDepositMap::Filler depFillerCalo(*caloDepMap);
       depFillerCalo.insert(muons, caloDeps.begin(), caloDeps.end());
       depFillerCalo.fill();
       event.put(std::move(caloDepMap), "caloIsoDeposits");
     }
 
     //event.put(std::move(combinedRelativeDeps, "combinedRelativeIsoDeposits");
     edm::ValueMap<double>::Filler depFillerCombRel(*combinedRelativeDepMap);
     depFillerCombRel.insert(muons, combinedRelativeDeps.begin(), combinedRelativeDeps.end());
     depFillerCombRel.fill();
     event.put(std::move(combinedRelativeDepMap), "combinedRelativeIsoDeposits");
   }
   edm::ValueMap<bool>::Filler isoFiller(*comboIsoDepMap);
   isoFiller.insert(muons, combinedRelativeIsos.begin(), combinedRelativeIsos.end());
   isoFiller.fill();
   event.put(std::move(comboIsoDepMap));
 
   if (printDebug)
     std::cout << " END OF EVENT "
               << "================================";
 }