TrackingParticleSelector Class Reference

#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() [1/2]

TrackingParticleSelector::TrackingParticleSelector ( )

inline

Definition at line 18 of file TrackingParticleSelector.h.

{}

TrackingParticleSelector() [2/2]

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_2022v15_cff::maxPhi, and HLT_2022v15_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

isCharged()

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

inline

Definition at line 67 of file TrackingParticleSelector.h.

References funct::false, and cmsswSequenceInfo::tp.

Referenced by select().

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

isInTime()

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

inline

Definition at line 69 of file TrackingParticleSelector.h.

References cmsswSequenceInfo::tp.

Referenced by select().

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

isSignal()

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

inline

Definition at line 71 of file TrackingParticleSelector.h.

References cmsswSequenceInfo::tp.

Referenced by select().

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

isStable()

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

inline

Definition at line 75 of file TrackingParticleSelector.h.

References funct::abs(), dqmiolumiharvest::j, EgammaValidation_cff::pdgid, and cmsswSequenceInfo::tp.

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

operator()() [1/2]

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(), and cmsswSequenceInfo::tp.

{ return select(&tp); }

operator()() [2/2]

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

inline

Definition at line 93 of file TrackingParticleSelector.h.

References select(), and cmsswSequenceInfo::tp.

{ return select(tp); }

select()

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_, stableOnly_, and cmsswSequenceInfo::tp.

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

selectKinematics()

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

inline

Definition at line 120 of file TrackingParticleSelector.h.

References funct::abs(), SiPixelRawToDigiRegional_cfi::deltaPhi, PVValHelper::eta, invertRapidityCut_, lip_, maxRapidity_, meanPhi_, minHit_, minRapidity_, AlCaHLTBitMon_ParallelJobs::p, HLT_2022v15_cff::pt2, ptMax2_, ptMin2_, rangePhi_, tip2_, and cmsswSequenceInfo::tp.

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_);
  }

selectParticleType()

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

inline

Definition at line 142 of file TrackingParticleSelector.h.

References EgammaValidation_cff::pdgid, pdgId_, and cmsswSequenceInfo::tp.

Referenced by select().

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

Member Data Documentation

chargedOnly_

bool TrackingParticleSelector::chargedOnly_

private

Definition at line 168 of file TrackingParticleSelector.h.

Referenced by select().

intimeOnly_

bool TrackingParticleSelector::intimeOnly_

private

Definition at line 167 of file TrackingParticleSelector.h.

Referenced by select().

invertRapidityCut_

bool TrackingParticleSelector::invertRapidityCut_

private

Definition at line 171 of file TrackingParticleSelector.h.

Referenced by selectKinematics().

lip_

double TrackingParticleSelector::lip_

private

Definition at line 164 of file TrackingParticleSelector.h.

Referenced by selectKinematics().

maxRapidity_

float TrackingParticleSelector::maxRapidity_

private

Definition at line 160 of file TrackingParticleSelector.h.

Referenced by selectKinematics().

meanPhi_

float TrackingParticleSelector::meanPhi_

private

Definition at line 161 of file TrackingParticleSelector.h.

Referenced by selectKinematics().

minHit_

int TrackingParticleSelector::minHit_

private

Definition at line 165 of file TrackingParticleSelector.h.

Referenced by selectKinematics().

minRapidity_

float TrackingParticleSelector::minRapidity_

private

Definition at line 159 of file TrackingParticleSelector.h.

Referenced by selectKinematics().

pdgId_

std::vector<int> TrackingParticleSelector::pdgId_

private

Definition at line 170 of file TrackingParticleSelector.h.

Referenced by selectParticleType().

ptMax2_

double TrackingParticleSelector::ptMax2_

private

Definition at line 158 of file TrackingParticleSelector.h.

Referenced by selectKinematics().

ptMin2_

double TrackingParticleSelector::ptMin2_

private

Definition at line 157 of file TrackingParticleSelector.h.

Referenced by selectKinematics().

rangePhi_

float TrackingParticleSelector::rangePhi_

private

Definition at line 162 of file TrackingParticleSelector.h.

Referenced by selectKinematics().

signalOnly_

bool TrackingParticleSelector::signalOnly_

private

Definition at line 166 of file TrackingParticleSelector.h.

Referenced by select().

stableOnly_

bool TrackingParticleSelector::stableOnly_

private

Definition at line 169 of file TrackingParticleSelector.h.

Referenced by select().

tip2_

double TrackingParticleSelector::tip2_

private

Definition at line 163 of file TrackingParticleSelector.h.

Referenced by selectKinematics().