Inheritance diagram for TrackingParticleRefMuonProducer:

Public Member Functions
void	produce (edm::Event &iEvent, const edm::EventSetup &iSetup) override

	TrackingParticleRefMuonProducer (const edm::ParameterSet &iConfig)

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	pmin_

double	ptmin_

std::string	skim_

edm::EDGetTokenT < TrackingParticleCollection >	tpToken_

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 11 of file TrackingParticleRefMuonProducer.cc.

Constructor & Destructor Documentation

TrackingParticleRefMuonProducer::TrackingParticleRefMuonProducer ( const edm::ParameterSet & iConfig )

Definition at line 24 of file TrackingParticleRefMuonProducer.cc.

References pmin_, ptmin_, and skim_.

     : tpToken_(consumes<TrackingParticleCollection>(iConfig.getParameter<edm::InputTag>("src"))),
       skim_(iConfig.getParameter<std::string>("skim")),
       ptmin_(iConfig.getParameter<double>("ptmin")),
       pmin_(iConfig.getParameter<double>("pmin")) {
   edm::LogVerbatim("TrackingParticleRefMuonProducer")
       << "\n constructing TrackingParticleRefMuonProducer: skim = " << skim_;
   if (skim_ == "mu")
     edm::LogVerbatim("TrackingParticleRefMuonProducer") << "\t ptmin = " << ptmin_ << ", pmin = " << pmin_ << "\n";
   else if (skim_ == "track")
     edm::LogVerbatim("TrackingParticleRefMuonProducer") << "\t ptmin = " << ptmin_ << "\n";
   else if (skim_ == "pf")
     edm::LogVerbatim("TrackingParticleRefMuonProducer") << "\t ptmin = " << ptmin_ << ", pmin = " << pmin_ << "\n";
   else
     edm::LogError("TrackingParticleRefMuonProducer") << "\t undefined skim = " << skim_ << "\n";
 
   produces<TrackingParticleRefVector>();
 }

Member Function Documentation

void TrackingParticleRefMuonProducer::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

override

Definition at line 43 of file TrackingParticleRefMuonProducer.cc.

References funct::abs(), dataset::end, edm::Event::getByToken(), mps_fire::i, eostools::move(), pmin_, ptmin_, edm::Event::put(), skim_, cmsswSequenceInfo::tp, and tpToken_.

                                                                                            {
   edm::Handle<TrackingParticleCollection> tpH;
   iEvent.getByToken(tpToken_, tpH);
 
   auto tpskim = std::make_unique<TrackingParticleRefVector>();
 
   if (skim_ == "mu") {
     for (size_t i = 0, end = tpH->size(); i < end; ++i) {
       auto tp = TrackingParticleRef(tpH, i);
 
       // test if the TP is a muon with pt and p above minimum thresholds
       bool isMu = (std::abs(tp->pdgId()) == 13);
       bool ptpOk = (tp->pt() > ptmin_) && (tp->p() > pmin_);
       if (isMu && ptpOk)
         tpskim->push_back(tp);
       else {
         // test if the TP has muon hits
         int n_muon_hits = tp->numberOfHits() - tp->numberOfTrackerHits();
         if (n_muon_hits > 0)
           tpskim->push_back(tp);
       }
     }
   } else if (skim_ == "track") {
     for (size_t i = 0, end = tpH->size(); i < end; ++i) {
       auto tp = TrackingParticleRef(tpH, i);
 
       // test if the pt is above a minimum cut
       if (tp->pt() > ptmin_)
         tpskim->push_back(tp);
     }
   } else if (skim_ == "pf") {
     for (size_t i = 0, end = tpH->size(); i < end; ++i) {
       auto tp = TrackingParticleRef(tpH, i);
 
       // test if p and pt are above minimum cuts
       if ((tp->pt() > ptmin_) && (tp->p() > pmin_))
         tpskim->push_back(tp);
     }
   }
 
   iEvent.put(std::move(tpskim));
 }