#include <TrackingParticleSelector.h>

Public Member Functions
bool	isCharged (const TrackingParticle *tp) const

bool	isInTime (const TrackingParticle *tp) const

bool	isSignal (const TrackingParticle *tp) const

bool	isStable (const TrackingParticle *tp) const

bool	operator() (const TrackingParticle &tp) const

bool	operator() (const TrackingParticle *tp) const

bool	select (const TrackingParticle *tp) const

bool	selectKinematics (const TrackingParticle *tp) const

bool	selectParticleType (const TrackingParticle *tp) const

	TrackingParticleSelector ()

	TrackingParticleSelector (double ptMin, double ptMax, double minRapidity, double maxRapidity, double tip, double lip, int minHit, bool signalOnly, bool intimeOnly, bool chargedOnly, bool stableOnly, const std::vector< int > &pdgId=std::vector< int >(), bool invertRapidityCut=false, double minPhi=-3.2, double maxPhi=3.2)

Private Attributes
bool	chargedOnly_

bool	intimeOnly_

bool	invertRapidityCut_

double	lip_

float	maxRapidity_

float	meanPhi_

int	minHit_

float	minRapidity_

std::vector< int >	pdgId_

double	ptMax2_

double	ptMin2_

float	rangePhi_

bool	signalOnly_

bool	stableOnly_

double	tip2_

Detailed Description

Definition at line 16 of file TrackingParticleSelector.h.

Constructor & Destructor Documentation

TrackingParticleSelector::TrackingParticleSelector ( )

inline

Definition at line 18 of file TrackingParticleSelector.h.

18 {}

TrackingParticleSelector::TrackingParticleSelector	(	double	ptMin,
		double	ptMax,
		double	minRapidity,
		double	maxRapidity,
		double	tip,
		double	lip,
		int	minHit,
		bool	signalOnly,
		bool	intimeOnly,
		bool	chargedOnly,
		bool	stableOnly,
		const std::vector< int > &	pdgId = `std::vector<int>()`,
		bool	invertRapidityCut = `false`,
		double	minPhi = `-3.2`,
		double	maxPhi = `3.2`
	)

inline

Definition at line 19 of file TrackingParticleSelector.h.

References Exception, M_PI, HLT_FULL_cff::maxPhi, and HLT_FULL_cff::minPhi.

       : ptMin2_(ptMin * ptMin),
         ptMax2_(ptMax * ptMax),
         minRapidity_(minRapidity),
         maxRapidity_(maxRapidity),
         meanPhi_((minPhi + maxPhi) / 2.),
         rangePhi_((maxPhi - minPhi) / 2.),
         tip2_(tip * tip),
         lip_(lip),
         minHit_(minHit),
         signalOnly_(signalOnly),
         intimeOnly_(intimeOnly),
         chargedOnly_(chargedOnly),
         stableOnly_(stableOnly),
         pdgId_(pdgId),
         invertRapidityCut_(invertRapidityCut) {
     if (minPhi >= maxPhi) {
       throw cms::Exception("Configuration")
           << "TrackingParticleSelector: minPhi (" << minPhi << ") must be smaller than maxPhi (" << maxPhi
           << "). The range is constructed from minPhi to maxPhi around their "
              "average.";
     }
     if (minPhi >= M_PI) {
       throw cms::Exception("Configuration") << "TrackingParticleSelector: minPhi (" << minPhi
                                             << ") must be smaller than PI. The range is constructed from minPhi "
                                                "to maxPhi around their average.";
     }
     if (maxPhi <= -M_PI) {
       throw cms::Exception("Configuration") << "TrackingParticleSelector: maxPhi (" << maxPhi
                                             << ") must be larger than -PI. The range is constructed from minPhi "
                                                "to maxPhi around their average.";
     }
   }

Member Function Documentation

bool TrackingParticleSelector::isCharged ( const TrackingParticle * tp ) const

inline

Definition at line 67 of file TrackingParticleSelector.h.

References TrackingParticle::charge().

Referenced by select().

67 { return (tp->charge() == 0 ? false : true); }

TrackingParticle::charge

float charge() const

Electric charge. Note this is taken from the first SimTrack only.

Definition: TrackingParticle.h:98

bool TrackingParticleSelector::isInTime ( const TrackingParticle * tp ) const

inline

Definition at line 69 of file TrackingParticleSelector.h.

References EncodedEventId::bunchCrossing(), and TrackingParticle::eventId().

Referenced by select().

69 { return (tp->eventId().bunchCrossing() == 0); }

EncodedEventId::bunchCrossing

int bunchCrossing() const

get the detector field from this detid

Definition: EncodedEventId.h:25

TrackingParticle::eventId

EncodedEventId eventId() const

Signal source, crossing number.

Definition: TrackingParticle.h:72

bool TrackingParticleSelector::isSignal ( const TrackingParticle * tp ) const

inline

Definition at line 71 of file TrackingParticleSelector.h.

References EncodedEventId::bunchCrossing(), EncodedEventId::event(), and TrackingParticle::eventId().

Referenced by select().

                                                   {
     return (tp->eventId().bunchCrossing() == 0 && tp->eventId().event() == 0);
   }

bool TrackingParticleSelector::isStable ( const TrackingParticle * tp ) const

inline

Definition at line 75 of file TrackingParticleSelector.h.

References funct::abs(), TrackingParticle::genParticle_begin(), TrackingParticle::genParticle_end(), dqmiolumiharvest::j, TrackingParticle::pdgId(), and TrackingParticle::status().

Referenced by select().

                                                   {
     for (TrackingParticle::genp_iterator j = tp->genParticle_begin(); j != tp->genParticle_end(); ++j) {
       if (j->get() == nullptr || j->get()->status() != 1) {
         return false;
       }
     }
     // test for remaining unstabled due to lack of genparticle pointer
     int pdgid = std::abs(tp->pdgId());
     if (tp->status() == -99 && (pdgid != 11 && pdgid != 13 && pdgid != 211 && pdgid != 321 && pdgid != 2212 &&
                                 pdgid != 3112 && pdgid != 3222 && pdgid != 3312 && pdgid != 3334)) {
       return false;
     }
     return true;
   }

bool TrackingParticleSelector::operator() ( const TrackingParticle & tp ) const

inline

Operator() performs the selection: e.g. if (tPSelector(tp)) {...} https://stackoverflow.com/questions/14466620/c-template-specialization-calling-methods-on-types-that-could-be-pointers-or/14466705

Definition at line 92 of file TrackingParticleSelector.h.

References select().

92 { return select(&tp); }

TrackingParticleSelector::select

bool select(const TrackingParticle *tp) const

Definition: TrackingParticleSelector.h:95

bool TrackingParticleSelector::operator() ( const TrackingParticle * tp ) const

inline

Definition at line 93 of file TrackingParticleSelector.h.

References select().

93 { return select(tp); }

TrackingParticleSelector::select

bool select(const TrackingParticle *tp) const

Definition: TrackingParticleSelector.h:95

bool TrackingParticleSelector::select ( const TrackingParticle * tp ) const

inline

Definition at line 95 of file TrackingParticleSelector.h.

References chargedOnly_, intimeOnly_, isCharged(), isInTime(), isSignal(), isStable(), selectKinematics(), selectParticleType(), signalOnly_, and stableOnly_.

Referenced by operator()().

                                                 {
     // signal only means no PU particles
     if (signalOnly_ && !isSignal(tp))
       return false;
     // intime only means no OOT PU particles
     if (intimeOnly_ && !isInTime(tp))
       return false;
 
     // select only if charge!=0
     if (chargedOnly_ && !isCharged(tp))
       return false;
 
     // select for particle type
     if (!selectParticleType(tp)) {
       return false;
     }
 
     // select only stable particles
     if (stableOnly_ && !isStable(tp)) {
       return false;
     }
 
     return selectKinematics(tp);
   }

bool TrackingParticleSelector::selectKinematics ( const TrackingParticle * tp ) const

inline

Definition at line 120 of file TrackingParticleSelector.h.

References funct::abs(), srCondWrite_cfg::deltaPhi, PVValHelper::eta, invertRapidityCut_, lip_, maxRapidity_, meanPhi_, minHit_, minRapidity_, TrackingParticle::numberOfTrackerLayers(), AlCaHLTBitMon_ParallelJobs::p, TrackingParticle::p4(), ptMax2_, ptMin2_, rangePhi_, tip2_, and TrackingParticle::vertex().

Referenced by select().

                                                           {
     auto etaOk = [&](const TrackingParticle *p) -> bool {
       float eta = etaFromXYZ(p->px(), p->py(), p->pz());
       if (!invertRapidityCut_)
         return (eta >= minRapidity_) && (eta <= maxRapidity_);
       else
         return (eta < minRapidity_ || eta > maxRapidity_);
     };
     auto phiOk = [&](const TrackingParticle *p) {
       float dphi = deltaPhi(atan2f(p->py(), p->px()), meanPhi_);
       return dphi >= -rangePhi_ && dphi <= rangePhi_;
     };
     auto ptOk = [&](const TrackingParticle *p) {
       double pt2 = tp->p4().perp2();
       return pt2 >= ptMin2_ && pt2 <= ptMax2_;
     };
     return (tp->numberOfTrackerLayers() >= minHit_ && ptOk(tp) && etaOk(tp) && phiOk(tp) &&
             std::abs(tp->vertex().z()) <= lip_ &&  // vertex last to avoid to load it if not striclty
                                                    // necessary...
             tp->vertex().perp2() <= tip2_);
   }

bool TrackingParticleSelector::selectParticleType ( const TrackingParticle * tp ) const

inline

Definition at line 142 of file TrackingParticleSelector.h.

References TrackingParticle::pdgId(), and pdgId_.

Referenced by select().

                                                             {
     auto pdgid = tp->pdgId();
     if (!pdgId_.empty()) {
       for (auto id : pdgId_) {
         if (id == pdgid) {
           return true;
         }
       }
     } else {
       return true;
     }
     return false;
   }