#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	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () 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

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

     : 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");
 }

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.

◆ ~L3MuonIsolationProducer()

L3MuonIsolationProducer::~L3MuonIsolationProducer ( )

override

destructor

Definition at line 85 of file L3MuonIsolationProducer.cc.

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

References LogDebug.

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.

                                                                                 {
   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 "
                     << "================================";
 }