egammaisolation::EgammaTrackSelector Class Reference

#include <EgammaTrackSelector.h>

Classes
struct	Parameters
	config parameters More...

Public Types
enum	{ dz = 0, vz, bs, vtx }
typedef reco::TrackBase::Point	BeamPoint
typedef edm::View< reco::Track >	input_type
typedef std::pair< float, float >	Range
typedef std::list< const reco::Track * >	result_type

Public Member Functions
	EgammaTrackSelector (const Parameters &pars)
result_type	operator() (const input_type &tracks) const

Private Attributes
Parameters	thePars

Detailed Description

Definition at line 13 of file EgammaTrackSelector.h.

Member Typedef Documentation

BeamPoint

typedef reco::TrackBase::Point egammaisolation::EgammaTrackSelector::BeamPoint

Definition at line 18 of file EgammaTrackSelector.h.

input_type

typedef edm::View<reco::Track> egammaisolation::EgammaTrackSelector::input_type

Definition at line 17 of file EgammaTrackSelector.h.

Range

typedef std::pair<float, float> egammaisolation::EgammaTrackSelector::Range

Definition at line 15 of file EgammaTrackSelector.h.

result_type

typedef std::list<const reco::Track*> egammaisolation::EgammaTrackSelector::result_type

Definition at line 16 of file EgammaTrackSelector.h.

Member Enumeration Documentation

anonymous enum

anonymous enum

Enumerator
dz
vz
bs
vtx

Definition at line 20 of file EgammaTrackSelector.h.

{ dz = 0, vz, bs, vtx };

Constructor & Destructor Documentation

EgammaTrackSelector()

egammaisolation::EgammaTrackSelector::EgammaTrackSelector ( const Parameters &  pars )

inline

Definition at line 60 of file EgammaTrackSelector.h.

: thePars(pars) {}

Member Function Documentation

operator()()

EgammaTrackSelector::result_type EgammaTrackSelector::operator()

(

const input_type &

tracks

)

const

pick/read variables in order to cut down on unnecessary calls someone will have some fun reading the log if Debug is on the biggest reason is the numberOfValidHits call (the rest are not as costly)


access to the remaining vars is slow

skip if min Hits == 0; assumes any track has at least one valid hit

similarly here

Definition at line 15 of file EgammaTrackSelector.cc.

References cms::cuda::bs, ChiSquaredProbability(), PVValHelper::dz, metname, mps_fire::result, AlCaHLTBitMon_QueryRunRegistry::string, and pwdgSkimBPark_cfi::tracks.

                                                                                                                {
  static const std::string metname = "EgammaIsolationAlgos|EgammaTrackSelector";
  result_type result;
  for (input_type::const_iterator it = tracks.begin(); it != tracks.end(); it++) {

    float tZ;
    switch (thePars.dzOption) {
      case dz:
        tZ = it->dz();
        break;
      case vz:
        tZ = it->vz();
        break;
      case bs:
        tZ = it->dz(thePars.beamPoint);
        break;
      default:
        tZ = it->vz();
        break;
    }

    float tPt = it->pt();
    //float tD0 = fabs(it->d0());  //currently not used.
    float tD0Cor = fabs(it->dxy(thePars.beamPoint));
    float tEta = it->eta();
    float tPhi = it->phi();
    float tChi2Ndof = it->normalizedChi2();


    if (!inside(tZ, thePars.zRange))
      continue;
    if (tPt < thePars.ptMin)
      continue;
    if (!inside(tD0Cor, thePars.rRange))
      continue;
    if (thePars.dir.deltaR(reco::isodeposit::Direction(tEta, tPhi)) > thePars.drMax)
      continue;
    if (tChi2Ndof > thePars.chi2NdofMax)
      continue;

    if (thePars.nHitsMin > 0) {
      unsigned int tHits = it->numberOfValidHits();
      if (tHits < thePars.nHitsMin)
        continue;
    }

    if (thePars.chi2ProbMin > 0) {
      float tChi2Prob = ChiSquaredProbability(it->chi2(), it->ndof());
      if (tChi2Prob < thePars.chi2ProbMin)
        continue;
    }
    result.push_back(&*it);
  }
  return result;
}

Member Data Documentation

thePars

Parameters egammaisolation::EgammaTrackSelector::thePars

private

Definition at line 65 of file EgammaTrackSelector.h.