#include <L3MuonIsolationProducer.h>

Inheritance diagram for L3MuonIsolationProducer:

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

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

	~L3MuonIsolationProducer () 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

Private Attributes
double	optOutputIsoDeposits

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

edm::ParameterSet	theConfig

muonisolation::Cuts	theCuts

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

int	theMaxNTracks

edm::InputTag	theMuonCollectionLabel

edm::EDGetTokenT< reco::RecoChargedCandidateCollection >	theMuonCollectionToken

double	theTrackPt_Min

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 24 of file L3MuonIsolationProducer.h.

Constructor & Destructor Documentation

◆ L3MuonIsolationProducer()

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

constructor with config

get min pt for the track to go into sumPt

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

Definition at line 36 of file L3MuonIsolationProducer.cc.

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

     : theConfig(par),
       theMuonCollectionLabel(par.getParameter<InputTag>("inputMuonCollection")),
       optOutputIsoDeposits(par.getParameter<bool>("OutputMuIsoDeposits")),
       theTrackPt_Min(-1) {
   LogDebug("RecoMuon|L3MuonIsolationProducer") << " L3MuonIsolationProducer CTOR";
 
   theMuonCollectionToken = consumes<RecoChargedCandidateCollection>(theMuonCollectionLabel);
 
   if (optOutputIsoDeposits)
     produces<reco::IsoDepositMap>();
   produces<edm::ValueMap<bool>>();
 
   //
   // Extractor
   //
   edm::ParameterSet extractorPSet = theConfig.getParameter<edm::ParameterSet>("ExtractorPSet");
   theTrackPt_Min = theConfig.getParameter<double>("TrackPt_Min");
   std::string extractorName = extractorPSet.getParameter<std::string>("ComponentName");
   theExtractor = std::unique_ptr<reco::isodeposit::IsoDepositExtractor>{
       IsoDepositExtractorFactory::get()->create(extractorName, extractorPSet, consumesCollector())};
   std::string depositType = extractorPSet.getUntrackedParameter<std::string>("DepositLabel");
 
   //
   // Cuts
   //
   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("L3MuonIsolationProducer::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");
 }

◆ ~L3MuonIsolationProducer()

L3MuonIsolationProducer::~L3MuonIsolationProducer ( )

override

destructor

Definition at line 85 of file L3MuonIsolationProducer.cc.

References LogDebug.

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

Member Function Documentation

◆ produce()

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

override

Produce isolation maps.

Definition at line 89 of file L3MuonIsolationProducer.cc.

References submitPVResolutionJobs::count, PA_MinBiasSkim_cff::cut, reco::IsoDeposit::depositAndCountWithin(), symbols::deps, options_cfi::eventSetup, edm::helper::Filler< Map >::fill(), mps_fire::i, edm::helper::Filler< Map >::insert(), LogDebug, LogTrace, metname, eostools::move(), amptDefaultParameters_cff::mu, PDWG_BPHSkim_cff::muons, custom_jme_cff::nMuons, optOutputIsoDeposits, reco::IsoDeposit::print(), mps_fire::result, AlCaHLTBitMon_QueryRunRegistry::string, theApplyCutsORmaxNTracks, theCuts, theExtractor, theMaxNTracks, theMuonCollectionLabel, theMuonCollectionToken, theTrackPt_Min, and boostedElectronIsolation_cff::vetos.

                                                                                 {
   std::string metname = "RecoMuon|L3MuonIsolationProducer";
 
   LogDebug(metname) << " L3 Muon Isolation producing..."
                     << " BEGINING OF EVENT "
                     << "================================";
 
   // Take the SA container
   LogTrace(metname) << " Taking the muons: " << theMuonCollectionLabel;
   Handle<TrackCollection> muons;
   event.getByToken(theMuonCollectionToken, muons);
 
   auto depMap = std::make_unique<reco::IsoDepositMap>();
   auto isoMap = std::make_unique<edm::ValueMap<bool>>();
 
   //
   // get Vetos and deposits
   //
   unsigned int nMuons = muons->size();
 
   IsoDeposit::Vetos vetos(nMuons);
 
   std::vector<IsoDeposit> deps(nMuons);
   std::vector<bool> isos(nMuons, false);
 
   for (unsigned int i = 0; i < nMuons; i++) {
     TrackRef mu(muons, i);
     deps[i] = theExtractor->deposit(event, eventSetup, *mu);
     vetos[i] = deps[i].veto();
   }
 
   //
   // add here additional vetos
   //
   //.....
 
   //
   // actual cut step
   //
   for (unsigned int iMu = 0; iMu < nMuons; ++iMu) {
     const reco::Track* mu = &(*muons)[iMu];
 
     const IsoDeposit& deposit = deps[iMu];
     LogTrace(metname) << deposit.print();
 
     const Cuts::CutSpec& cut = theCuts(mu->eta());
     std::pair<double, int> sumAndCount = deposit.depositAndCountWithin(cut.conesize, vetos, theTrackPt_Min);
 
     double value = sumAndCount.first;
     int count = sumAndCount.second;
 
     bool result = (value < cut.threshold);
     if (theApplyCutsORmaxNTracks)
       result |= count <= theMaxNTracks;
     LogTrace(metname) << "deposit in cone: " << value << "with count " << count << " is isolated: " << result;
 
     isos[iMu] = result;
   }
 
   //
   // store
   //
   if (optOutputIsoDeposits) {
     reco::IsoDepositMap::Filler depFiller(*depMap);
     depFiller.insert(muons, deps.begin(), deps.end());
     depFiller.fill();
     event.put(std::move(depMap));
   }
   edm::ValueMap<bool>::Filler isoFiller(*isoMap);
   isoFiller.insert(muons, isos.begin(), isos.end());
   isoFiller.fill();
   event.put(std::move(isoMap));
 
   LogTrace(metname) << " END OF EVENT "
                     << "================================";
 }