TrackingParticleSelector Class Reference

#include <TrackingParticleSelector.h>

Public Member Functions
bool	operator() (const TrackingParticle &tp) const
	Operator() performs the selection: e.g. if (tPSelector(tp)) {...}. More...
	TrackingParticleSelector ()
	TrackingParticleSelector (double ptMin, double minRapidity, double maxRapidity, double tip, double lip, int minHit, bool signalOnly, bool chargedOnly, bool stableOnly, const std::vector< int > &pdgId=std::vector< int >())

Private Attributes
bool	chargedOnly_
double	lip_
float	maxRapidity_
int	minHit_
float	minRapidity_
std::vector< int >	pdgId_
double	ptMin_
bool	signalOnly_
bool	stableOnly_
double	tip_

Detailed Description

Definition at line 15 of file TrackingParticleSelector.h.

Constructor & Destructor Documentation

TrackingParticleSelector::TrackingParticleSelector ( )

inline

Definition at line 18 of file TrackingParticleSelector.h.

{}

TrackingParticleSelector::TrackingParticleSelector ( double  ptMin, double  minRapidity, double  maxRapidity, double  tip, double  lip, int  minHit, bool  signalOnly, bool  chargedOnly, bool  stableOnly, const std::vector< int > &  pdgId = std::vector<int>()  )

inline

Definition at line 19 of file TrackingParticleSelector.h.

                                                                 :
    ptMin_( ptMin ), minRapidity_( minRapidity ), maxRapidity_( maxRapidity ),
    tip_( tip ), lip_( lip ), minHit_( minHit ), signalOnly_(signalOnly), chargedOnly_(chargedOnly), stableOnly_(stableOnly), pdgId_( pdgId ) { }

Member Function Documentation

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

inline

Operator() performs the selection: e.g. if (tPSelector(tp)) {...}.

Definition at line 26 of file TrackingParticleSelector.h.

References funct::abs(), EncodedEventId::bunchCrossing(), TrackingParticle::charge(), chargedOnly_, eta(), EncodedEventId::event(), TrackingParticle::eventId(), TrackingParticle::genParticle_begin(), TrackingParticle::genParticle_end(), j, lip_, maxRapidity_, minHit_, minRapidity_, TrackingParticle::momentum(), TrackingParticle::numberOfTrackerLayers(), AlCaHLTBitMon_ParallelJobs::p, TrackingParticle::pdgId(), pdgId_, ptMin_, signalOnly_, TopDecayID::stable, stableOnly_, TrackingParticle::status(), tip_, and TrackingParticle::vertex().

                                                       {
    auto pdgid = tp.pdgId();
    if (chargedOnly_ && tp.charge()==0) return false;//select only if charge!=0

    bool testId = false;
    unsigned int idSize = pdgId_.size();
    if (idSize==0) testId = true;
    else for (unsigned int it=0;it!=idSize;++it){
      if (pdgid==pdgId_[it]) { testId = true; break;}
    }

    bool signal = true;
    if (signalOnly_) signal = (tp.eventId().bunchCrossing()== 0 && tp.eventId().event() == 0); // signal only means no PU particles

    // select only stable particles
    bool stable = true;
    if (stableOnly_) {
      if (!signal) {
    stable = false; // we are not interested into PU particles among the stable ones
      } else {
    for( TrackingParticle::genp_iterator j = tp.genParticle_begin(); j != tp.genParticle_end(); ++ j ) {
      if (j->get()==0 || j->get()->status() != 1) {
        stable = false; break;
      }
    }
       // test for remaining unstabled due to lack of genparticle pointer
       if( (stable  & (tp.status() == -99) ) &&
          (std::abs(pdgid) != 11 && std::abs(pdgid) != 13 && std::abs(pdgid) != 211 &&
           std::abs(pdgid) != 321 && std::abs(pdgid) != 2212 && std::abs(pdgid) != 3112 &&
           std::abs(pdgid) != 3222 && std::abs(pdgid) != 3312 && std::abs(pdgid) != 3334)) stable = false;
      }
    }

    auto etaOk = [&](const TrackingParticle::Vector& p)->bool{ float eta= etaFromXYZ(p.x(),p.y(),p.z()); return (eta>= minRapidity_) & (eta<=maxRapidity_);};
    return (
            (testId & signal & stable) &&
        tp.numberOfTrackerLayers() >= minHit_ &&
        tp.momentum().perp2() >= ptMin_*ptMin_ &&
            etaOk(tp.momentum()) &&
            std::abs(tp.vertex().z()) <= lip_ &&   // vertex last to avoid to load it if not striclty necessary...
        tp.vertex().perp2() <= tip_*tip_ 
        );
  }

Member Data Documentation

bool TrackingParticleSelector::chargedOnly_

private

Definition at line 78 of file TrackingParticleSelector.h.

Referenced by operator()().

double TrackingParticleSelector::lip_

private

Definition at line 75 of file TrackingParticleSelector.h.

Referenced by operator()().

float TrackingParticleSelector::maxRapidity_

private

Definition at line 73 of file TrackingParticleSelector.h.

Referenced by operator()().

int TrackingParticleSelector::minHit_

private

Definition at line 76 of file TrackingParticleSelector.h.

Referenced by operator()().

float TrackingParticleSelector::minRapidity_

private

Definition at line 72 of file TrackingParticleSelector.h.

Referenced by operator()().

std::vector<int> TrackingParticleSelector::pdgId_

private

Definition at line 80 of file TrackingParticleSelector.h.

Referenced by operator()().

double TrackingParticleSelector::ptMin_

private

Definition at line 71 of file TrackingParticleSelector.h.

Referenced by operator()().

bool TrackingParticleSelector::signalOnly_

private

Definition at line 77 of file TrackingParticleSelector.h.

Referenced by operator()().

bool TrackingParticleSelector::stableOnly_

private

Definition at line 79 of file TrackingParticleSelector.h.

Referenced by operator()().

double TrackingParticleSelector::tip_

private

Definition at line 74 of file TrackingParticleSelector.h.

Referenced by operator()().