#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

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

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

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<>
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

Definition at line 20 of file L3MuonIsolationProducer.h.

Constructor & Destructor Documentation

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(), reco::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")),
   theExtractor(nullptr),
   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 = 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 ( )

override

destructor

Definition at line 89 of file L3MuonIsolationProducer.cc.

References LogDebug, and theExtractor.

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

Member Function Documentation

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

override

Produce isolation maps.

Definition at line 94 of file L3MuonIsolationProducer.cc.

References KineDebug3::count(), TkAlMuonSelectors_cfi::cut, reco::isodeposit::IsoDepositExtractor::deposit(), reco::IsoDeposit::depositAndCountWithin(), symbols::deps, reco::TrackBase::eta(), edm::helper::Filler< Map >::fill(), mps_fire::i, edm::helper::Filler< Map >::insert(), LogDebug, LogTrace, metname, eostools::move(), RPCpg::mu, extraflags_cff::muons, jets_cff::nMuons, optOutputIsoDeposits, reco::IsoDeposit::print(), mps_fire::result, AlCaHLTBitMon_QueryRunRegistry::string, theApplyCutsORmaxNTracks, theCuts, theExtractor, theMaxNTracks, theMuonCollectionLabel, theMuonCollectionToken, theTrackPt_Min, reco::IsoDeposit::veto(), 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 " <<"================================";
 }