CMS 3D CMS Logo

List of all members | Public Member Functions | Private Attributes
RecoTrackSelectorBase Class Reference

#include <RecoTrackSelectorBase.h>

Inheritance diagram for RecoTrackSelectorBase:
RecoTrackRefSelector RecoTrackSelector RecoTrackViewRefSelector

Public Member Functions

void init (const edm::Event &event, const edm::EventSetup &es)
 
bool operator() (const reco::TrackRef &tref) const
 
bool operator() (const reco::Track &t) const
 
bool operator() (const reco::Track &t, const reco::Track::Point &vertex) const
 
 RecoTrackSelectorBase ()
 
 RecoTrackSelectorBase (const edm::ParameterSet &cfg)
 
 RecoTrackSelectorBase (const edm::ParameterSet &cfg, edm::ConsumesCollector &iC)
 

Private Attributes

std::vector< reco::TrackBase::TrackAlgorithmalgorithm_
 
std::vector< reco::TrackBase::TrackAlgorithmalgorithmMask_
 
edm::EDGetTokenT< reco::BeamSpotbsSrcToken_
 
bool invertRapidityCut_
 
double lip_
 
double maxChi2_
 
double maxRapidity_
 
double meanPhi_
 
int min3DLayer_
 
int minHit_
 
int minLayer_
 
int minPixelHit_
 
double minRapidity_
 
std::vector< reco::TrackBase::TrackAlgorithmoriginalAlgorithm_
 
double ptMin_
 
std::vector< reco::TrackBase::TrackQualityquality_
 
double rangePhi_
 
double tip_
 
bool usePV_
 
reco::Track::Point vertex_
 
edm::EDGetTokenT< reco::VertexCollectionvertexToken_
 

Detailed Description

Definition at line 16 of file RecoTrackSelectorBase.h.

Constructor & Destructor Documentation

RecoTrackSelectorBase::RecoTrackSelectorBase ( )
inline

Definition at line 18 of file RecoTrackSelectorBase.h.

18 {}
RecoTrackSelectorBase::RecoTrackSelectorBase ( const edm::ParameterSet cfg)
inline

Definition at line 19 of file RecoTrackSelectorBase.h.

References reco::TrackBase::algoByName(), qcdUeDQM_cfi::algorithm, algorithm_, algorithmMask_, Exception, edm::ParameterSet::getParameter(), M_PI, HLT_2018_cff::maxPhi, HLT_2018_cff::minPhi, originalAlgorithm_, qcdUeDQM_cfi::quality, quality_, reco::TrackBase::qualityByName(), and AlCaHLTBitMon_QueryRunRegistry::string.

20  : ptMin_(cfg.getParameter<double>("ptMin")),
21  minRapidity_(cfg.getParameter<double>("minRapidity")),
22  maxRapidity_(cfg.getParameter<double>("maxRapidity")),
23  meanPhi_((cfg.getParameter<double>("minPhi") + cfg.getParameter<double>("maxPhi")) / 2.),
24  rangePhi_((cfg.getParameter<double>("maxPhi") - cfg.getParameter<double>("minPhi")) / 2.),
25  tip_(cfg.getParameter<double>("tip")),
26  lip_(cfg.getParameter<double>("lip")),
27  maxChi2_(cfg.getParameter<double>("maxChi2")),
28  minHit_(cfg.getParameter<int>("minHit")),
29  minPixelHit_(cfg.getParameter<int>("minPixelHit")),
30  minLayer_(cfg.getParameter<int>("minLayer")),
31  min3DLayer_(cfg.getParameter<int>("min3DLayer")),
32  usePV_(false),
33  invertRapidityCut_(cfg.getParameter<bool>("invertRapidityCut")) {
34  const auto minPhi = cfg.getParameter<double>("minPhi");
35  const auto maxPhi = cfg.getParameter<double>("maxPhi");
36  if (minPhi >= maxPhi) {
37  throw cms::Exception("Configuration")
38  << "RecoTrackSelectorPhase: minPhi (" << minPhi << ") must be smaller than maxPhi (" << maxPhi
39  << "). The range is constructed from minPhi to maxPhi around their average.";
40  }
41  if (minPhi >= M_PI) {
42  throw cms::Exception("Configuration")
43  << "RecoTrackSelectorPhase: minPhi (" << minPhi
44  << ") must be smaller than PI. The range is constructed from minPhi to maxPhi around their average.";
45  }
46  if (maxPhi <= -M_PI) {
47  throw cms::Exception("Configuration")
48  << "RecoTrackSelectorPhase: maxPhi (" << maxPhi
49  << ") must be larger than -PI. The range is constructed from minPhi to maxPhi around their average.";
50  }
51 
52  for (const std::string& quality : cfg.getParameter<std::vector<std::string> >("quality"))
54  for (const std::string& algorithm : cfg.getParameter<std::vector<std::string> >("algorithm"))
56  for (const std::string& algorithm : cfg.getParameter<std::vector<std::string> >("originalAlgorithm"))
58  for (const std::string& algorithm : cfg.getParameter<std::vector<std::string> >("algorithmMaskContains"))
60  }
T getParameter(std::string const &) const
std::vector< reco::TrackBase::TrackAlgorithm > algorithmMask_
std::vector< reco::TrackBase::TrackQuality > quality_
#define M_PI
static TrackQuality qualityByName(const std::string &name)
Definition: TrackBase.cc:126
std::vector< reco::TrackBase::TrackAlgorithm > algorithm_
static TrackAlgorithm algoByName(const std::string &name)
Definition: TrackBase.cc:137
std::vector< reco::TrackBase::TrackAlgorithm > originalAlgorithm_
RecoTrackSelectorBase::RecoTrackSelectorBase ( const edm::ParameterSet cfg,
edm::ConsumesCollector iC 
)
inline

Definition at line 62 of file RecoTrackSelectorBase.h.

References bsSrcToken_, edm::ConsumesCollector::consumes(), edm::ParameterSet::getParameter(), usePV_, and vertexToken_.

62  : RecoTrackSelectorBase(cfg) {
63  usePV_ = cfg.getParameter<bool>("usePV");
65  if (usePV_)
67  }
EDGetTokenT< ProductType > consumes(edm::InputTag const &tag)
T getParameter(std::string const &) const
std::vector< Vertex > VertexCollection
collection of Vertex objects
Definition: VertexFwd.h:9
edm::EDGetTokenT< reco::VertexCollection > vertexToken_
edm::EDGetTokenT< reco::BeamSpot > bsSrcToken_

Member Function Documentation

void RecoTrackSelectorBase::init ( const edm::Event event,
const edm::EventSetup es 
)
inline

Definition at line 69 of file RecoTrackSelectorBase.h.

References pwdgSkimBPark_cfi::beamSpot, bsSrcToken_, reco::BeamSpot::position(), usePV_, vertex_, and vertexToken_.

Referenced by RecoTrackViewRefSelector::select(), RecoTrackRefSelector::select(), and RecoTrackSelector::select().

69  {
71  event.getByToken(bsSrcToken_, beamSpot);
72  vertex_ = beamSpot->position();
73  if (!usePV_)
74  return;
76  event.getByToken(vertexToken_, hVtx);
77  if (hVtx->empty())
78  return;
79  vertex_ = (*hVtx)[0].position();
80  }
reco::Track::Point vertex_
edm::EDGetTokenT< reco::VertexCollection > vertexToken_
edm::EDGetTokenT< reco::BeamSpot > bsSrcToken_
const Point & position() const
position
Definition: BeamSpot.h:59
bool RecoTrackSelectorBase::operator() ( const reco::TrackRef tref) const
inline

Definition at line 82 of file RecoTrackSelectorBase.h.

82 { return (*this)(*tref); }
bool RecoTrackSelectorBase::operator() ( const reco::Track t) const
inline

Definition at line 84 of file RecoTrackSelectorBase.h.

References OrderedSet::t, and vertex_.

84 { return (*this)(t, vertex_); }
reco::Track::Point vertex_
bool RecoTrackSelectorBase::operator() ( const reco::Track t,
const reco::Track::Point vertex 
) const
inline

Definition at line 86 of file RecoTrackSelectorBase.h.

References ecalcalib_dqm_sourceclient-live_cfg::algo, reco::TrackBase::algo(), reco::TrackBase::algoMask(), algorithm_, algorithmMask_, SiPixelRawToDigiRegional_cfi::deltaPhi, reco::TrackBase::dsz(), reco::TrackBase::dxy(), PVValHelper::eta, spr::find(), reco::TrackBase::hitPattern(), mps_fire::i, invertRapidityCut_, lip_, maxChi2_, maxRapidity_, meanPhi_, min3DLayer_, minHit_, minLayer_, minPixelHit_, minRapidity_, reco::TrackBase::normalizedChi2(), reco::HitPattern::numberOfValidHits(), reco::HitPattern::numberOfValidPixelHits(), reco::HitPattern::numberOfValidStripLayersWithMonoAndStereo(), reco::TrackBase::originalAlgo(), originalAlgorithm_, AlCaHLTBitMon_ParallelJobs::p, reco::TrackBase::phi(), reco::HitPattern::pixelLayersWithMeasurement(), reco::TrackBase::pt(), ptMin_, reco::TrackBase::quality(), quality_, rangePhi_, tip_, and reco::HitPattern::trackerLayersWithMeasurement().

86  {
87  bool quality_ok = true;
88  if (!quality_.empty()) {
89  quality_ok = false;
90  for (unsigned int i = 0; i < quality_.size(); ++i) {
91  if (t.quality(quality_[i])) {
92  quality_ok = true;
93  break;
94  }
95  }
96  }
97 
98  bool algo_ok = true;
99  if (!algorithm_.empty()) {
100  if (std::find(algorithm_.begin(), algorithm_.end(), t.algo()) == algorithm_.end())
101  algo_ok = false;
102  }
103  if (!originalAlgorithm_.empty() && algo_ok) {
105  algo_ok = false;
106  }
107  if (!algorithmMask_.empty() && algo_ok) {
108  if (std::find_if(algorithmMask_.begin(),
109  algorithmMask_.end(),
110  // for some reason I have to either explicitly give the return type, or use static_cast<bool>()
111  [&](reco::TrackBase::TrackAlgorithm algo) -> bool { return t.algoMask()[algo]; }) ==
112  algorithmMask_.end())
113  algo_ok = false;
114  }
115 
116  const auto dphi = deltaPhi(t.phi(), meanPhi_);
117 
118  auto etaOk = [&](const reco::Track& p) -> bool {
119  float eta = p.eta();
120  if (!invertRapidityCut_)
121  return (eta >= minRapidity_) && (eta <= maxRapidity_);
122  else
123  return (eta < minRapidity_ || eta > maxRapidity_);
124  };
125 
126  return ((algo_ok & quality_ok) && t.hitPattern().numberOfValidHits() >= minHit_ &&
130  min3DLayer_ &&
131  fabs(t.pt()) >= ptMin_ && etaOk(t) && dphi >= -rangePhi_ && dphi <= rangePhi_ &&
132  fabs(t.dxy(vertex)) <= tip_ && fabs(t.dsz(vertex)) <= lip_ && t.normalizedChi2() <= maxChi2_);
133  }
double normalizedChi2() const
chi-squared divided by n.d.o.f. (or chi-squared * 1e6 if n.d.o.f. is zero)
Definition: TrackBase.h:572
int numberOfValidHits() const
Definition: HitPattern.h:787
double phi() const
azimuthal angle of momentum vector
Definition: TrackBase.h:614
int pixelLayersWithMeasurement() const
Definition: HitPattern.cc:492
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:19
double dsz() const
dsz parameter (THIS IS NOT the SZ impact parameter to (0,0,0) if refPoint is far from (0...
Definition: TrackBase.h:593
int trackerLayersWithMeasurement() const
Definition: HitPattern.cc:513
TrackAlgorithm
track algorithm
Definition: TrackBase.h:89
TrackAlgorithm algo() const
Definition: TrackBase.h:526
int numberOfValidStripLayersWithMonoAndStereo(uint16_t stripdet, uint16_t layer) const
Definition: HitPattern.cc:348
double pt() const
track transverse momentum
Definition: TrackBase.h:602
std::vector< reco::TrackBase::TrackAlgorithm > algorithmMask_
std::vector< reco::TrackBase::TrackQuality > quality_
TrackAlgorithm originalAlgo() const
Definition: TrackBase.h:527
AlgoMask algoMask() const
Definition: TrackBase.h:394
const HitPattern & hitPattern() const
Access the hit pattern, indicating in which Tracker layers the track has hits.
Definition: TrackBase.h:483
bool quality(const TrackQuality) const
Track quality.
Definition: TrackBase.h:531
std::vector< reco::TrackBase::TrackAlgorithm > algorithm_
int numberOfValidPixelHits() const
Definition: HitPattern.h:801
double dxy() const
dxy parameter. (This is the transverse impact parameter w.r.t. to (0,0,0) ONLY if refPoint is close t...
Definition: TrackBase.h:587
std::vector< reco::TrackBase::TrackAlgorithm > originalAlgorithm_

Member Data Documentation

std::vector<reco::TrackBase::TrackAlgorithm> RecoTrackSelectorBase::algorithm_
private

Definition at line 155 of file RecoTrackSelectorBase.h.

Referenced by operator()(), and RecoTrackSelectorBase().

std::vector<reco::TrackBase::TrackAlgorithm> RecoTrackSelectorBase::algorithmMask_
private

Definition at line 157 of file RecoTrackSelectorBase.h.

Referenced by operator()(), and RecoTrackSelectorBase().

edm::EDGetTokenT<reco::BeamSpot> RecoTrackSelectorBase::bsSrcToken_
private

Definition at line 151 of file RecoTrackSelectorBase.h.

Referenced by init(), and RecoTrackSelectorBase().

bool RecoTrackSelectorBase::invertRapidityCut_
private

Definition at line 149 of file RecoTrackSelectorBase.h.

Referenced by operator()().

double RecoTrackSelectorBase::lip_
private

Definition at line 142 of file RecoTrackSelectorBase.h.

Referenced by operator()().

double RecoTrackSelectorBase::maxChi2_
private

Definition at line 143 of file RecoTrackSelectorBase.h.

Referenced by operator()().

double RecoTrackSelectorBase::maxRapidity_
private

Definition at line 138 of file RecoTrackSelectorBase.h.

Referenced by operator()().

double RecoTrackSelectorBase::meanPhi_
private

Definition at line 139 of file RecoTrackSelectorBase.h.

Referenced by operator()().

int RecoTrackSelectorBase::min3DLayer_
private

Definition at line 147 of file RecoTrackSelectorBase.h.

Referenced by operator()().

int RecoTrackSelectorBase::minHit_
private

Definition at line 144 of file RecoTrackSelectorBase.h.

Referenced by operator()().

int RecoTrackSelectorBase::minLayer_
private

Definition at line 146 of file RecoTrackSelectorBase.h.

Referenced by operator()().

int RecoTrackSelectorBase::minPixelHit_
private

Definition at line 145 of file RecoTrackSelectorBase.h.

Referenced by operator()().

double RecoTrackSelectorBase::minRapidity_
private

Definition at line 137 of file RecoTrackSelectorBase.h.

Referenced by operator()().

std::vector<reco::TrackBase::TrackAlgorithm> RecoTrackSelectorBase::originalAlgorithm_
private

Definition at line 156 of file RecoTrackSelectorBase.h.

Referenced by operator()(), and RecoTrackSelectorBase().

double RecoTrackSelectorBase::ptMin_
private

Definition at line 136 of file RecoTrackSelectorBase.h.

Referenced by operator()().

std::vector<reco::TrackBase::TrackQuality> RecoTrackSelectorBase::quality_
private

Definition at line 154 of file RecoTrackSelectorBase.h.

Referenced by operator()(), and RecoTrackSelectorBase().

double RecoTrackSelectorBase::rangePhi_
private

Definition at line 140 of file RecoTrackSelectorBase.h.

Referenced by operator()().

double RecoTrackSelectorBase::tip_
private

Definition at line 141 of file RecoTrackSelectorBase.h.

Referenced by operator()().

bool RecoTrackSelectorBase::usePV_
private

Definition at line 148 of file RecoTrackSelectorBase.h.

Referenced by init(), and RecoTrackSelectorBase().

reco::Track::Point RecoTrackSelectorBase::vertex_
private

Definition at line 159 of file RecoTrackSelectorBase.h.

Referenced by init(), and operator()().

edm::EDGetTokenT<reco::VertexCollection> RecoTrackSelectorBase::vertexToken_
private

Definition at line 152 of file RecoTrackSelectorBase.h.

Referenced by init(), and RecoTrackSelectorBase().