CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Classes | Public Types | Public Member Functions | Private Member Functions | Private Attributes
AlignmentTrackSelector Class Reference

#include <AlignmentTrackSelector.h>

Classes

struct  ComparePt
 compare two tracks in pt (used by theNHighestPtTracks) More...
 

Public Types

typedef std::vector< const
reco::Track * > 
Tracks
 

Public Member Functions

 AlignmentTrackSelector (const edm::ParameterSet &cfg)
 constructor More...
 
Tracks select (const Tracks &tracks, const edm::Event &evt) const
 select tracks More...
 
bool useThisFilter ()
 returns if any of the Filters is used. More...
 
 ~AlignmentTrackSelector ()
 destructor More...
 

Private Member Functions

Tracks basicCuts (const Tracks &tracks, const edm::Event &evt) const
 apply basic cuts on pt,eta,phi,nhit More...
 
Tracks checkPrescaledHits (const Tracks &tracks, const edm::Event &evt) const
 
bool detailedHitsCheck (const reco::Track *track, const edm::Event &evt) const
 checking hit requirements beyond simple number of valid hits More...
 
bool isHit2D (const TrackingRecHit &hit) const
 
bool isIsolated (const TrackingRecHit *therechit, const edm::Event &evt) const
 
bool isOkCharge (const TrackingRecHit *therechit) const
 if valid, check for minimum charge (currently only in strip), if invalid give true More...
 
bool isOkChargeStripHit (const SiStripRecHit1D &siStripRecHit1D) const
 
bool isOkChargeStripHit (const SiStripRecHit2D &siStripRecHit2D) const
 
bool isOkTrkQuality (const reco::Track *track) const
 
Tracks theNHighestPtTracks (const Tracks &tracks) const
 filter the n highest pt tracks More...
 

Private Attributes

const bool applyBasicCuts_
 
const bool applyIsolation_
 
bool applyIterStepCheck_
 
const bool applyMultiplicityFilter_
 
const bool applyNHighestPt_
 
const bool applyPrescaledHitsFilter_
 
bool applyTrkQualityCheck_
 
const bool chargeCheck_
 
const double chi2nMax_
 
const edm::InputTag clusterValueMapTag_
 
const bool countStereoHitAs2D_
 
const double d0Max_
 
const double d0Min_
 
const double dzMax_
 
const double dzMin_
 
const double etaMax_
 
const double etaMin_
 
const edm::InputTag matchedrecHitsTag_
 
const double maxHitDiffEndcaps_
 
const int maxMultiplicity_
 
const double minHitChargeStrip_
 
const double minHitIsolation_
 
const int minHitsinBPIX_
 
const int minHitsinENDCAP_
 
const int minHitsinENDCAPminus_
 
const int minHitsinENDCAPplus_
 
const int minHitsinFPIX_
 
const int minHitsinFPIXminus_
 
const int minHitsinFPIXplus_
 
const int minHitsinPIX_
 
const int minHitsinTEC_
 
const int minHitsinTECminus_
 
const int minHitsinTECplus_
 
const int minHitsinTIB_
 
const int minHitsinTID_
 
const int minHitsinTIDminus_
 
const int minHitsinTIDplus_
 
const int minHitsinTOB_
 
const int minMultiplicity_
 
const int minPrescaledHits_
 
const bool multiplicityOnInput_
 
const int nHighestPt_
 
const double nHitMax_
 
const unsigned int nHitMin2D_
 
const double nHitMin_
 
const double nLostHitMax_
 
const double phiMax_
 
const double phiMin_
 
const double pMax_
 
const double pMin_
 
ComparePt ptComparator
 
const double ptMax_
 
const double ptMin_
 if true, cut min/maxMultiplicity on input instead of on final result More...
 
const edm::InputTag rphirecHitsTag_
 
const bool seedOnlyFromAbove_
 
const int theCharge_
 
std::vector
< reco::TrackBase::TrackQuality
trkQualities_
 
std::vector
< reco::TrackBase::TrackAlgorithm
trkSteps_
 

Detailed Description

Definition at line 17 of file AlignmentTrackSelector.h.

Member Typedef Documentation

typedef std::vector<const reco::Track*> AlignmentTrackSelector::Tracks

Definition at line 22 of file AlignmentTrackSelector.h.

Constructor & Destructor Documentation

AlignmentTrackSelector::AlignmentTrackSelector ( const edm::ParameterSet cfg)

constructor

Definition at line 32 of file AlignmentTrackSelector.cc.

References reco::TrackBase::algoByName(), applyBasicCuts_, applyIsolation_, applyIterStepCheck_, applyMultiplicityFilter_, applyNHighestPt_, applyPrescaledHitsFilter_, applyTrkQualityCheck_, chargeCheck_, chi2nMax_, clusterValueMapTag_, countStereoHitAs2D_, etaMax_, etaMin_, edm::ParameterSet::getParameter(), i, maxHitDiffEndcaps_, maxMultiplicity_, minHitChargeStrip_, minHitIsolation_, minHitsinBPIX_, minHitsinENDCAP_, minHitsinENDCAPminus_, minHitsinENDCAPplus_, minHitsinFPIX_, minHitsinFPIXminus_, minHitsinFPIXplus_, minHitsinPIX_, minHitsinTEC_, minHitsinTECminus_, minHitsinTECplus_, minHitsinTIB_, minHitsinTID_, minHitsinTIDminus_, minHitsinTIDplus_, minHitsinTOB_, minMultiplicity_, minPrescaledHits_, multiplicityOnInput_, nHighestPt_, nHitMax_, nHitMin2D_, nHitMin_, nLostHitMax_, phiMax_, phiMin_, pMax_, pMin_, ptMax_, ptMin_, reco::TrackBase::qualityByName(), trkQualities_, and trkSteps_.

32  :
33  applyBasicCuts_( cfg.getParameter<bool>( "applyBasicCuts" ) ),
34  applyNHighestPt_( cfg.getParameter<bool>( "applyNHighestPt" ) ),
35  applyMultiplicityFilter_( cfg.getParameter<bool>( "applyMultiplicityFilter" ) ),
36  seedOnlyFromAbove_( cfg.getParameter<int>( "seedOnlyFrom" ) ),
37  applyIsolation_( cfg.getParameter<bool>( "applyIsolationCut" ) ),
38  chargeCheck_( cfg.getParameter<bool>( "applyChargeCheck" ) ),
39  nHighestPt_( cfg.getParameter<int>( "nHighestPt" ) ),
40  minMultiplicity_ ( cfg.getParameter<int>( "minMultiplicity" ) ),
41  maxMultiplicity_ ( cfg.getParameter<int>( "maxMultiplicity" ) ),
42  multiplicityOnInput_ ( cfg.getParameter<bool>( "multiplicityOnInput" ) ),
43  ptMin_( cfg.getParameter<double>( "ptMin" ) ),
44  ptMax_( cfg.getParameter<double>( "ptMax" ) ),
45  pMin_( cfg.getParameter<double>( "pMin" ) ),
46  pMax_( cfg.getParameter<double>( "pMax" ) ),
47  etaMin_( cfg.getParameter<double>( "etaMin" ) ),
48  etaMax_( cfg.getParameter<double>( "etaMax" ) ),
49  phiMin_( cfg.getParameter<double>( "phiMin" ) ),
50  phiMax_( cfg.getParameter<double>( "phiMax" ) ),
51  nHitMin_( cfg.getParameter<double>( "nHitMin" ) ),
52  nHitMax_( cfg.getParameter<double>( "nHitMax" ) ),
53  chi2nMax_( cfg.getParameter<double>( "chi2nMax" ) ),
54  d0Min_( cfg.getParameter<double>( "d0Min" ) ),
55  d0Max_( cfg.getParameter<double>( "d0Max" ) ),
56  dzMin_( cfg.getParameter<double>( "dzMin" ) ),
57  dzMax_( cfg.getParameter<double>( "dzMax" ) ),
58  theCharge_( cfg.getParameter<int>( "theCharge" ) ),
59  minHitChargeStrip_( cfg.getParameter<double>( "minHitChargeStrip" ) ),
60  minHitIsolation_( cfg.getParameter<double>( "minHitIsolation" ) ),
61  rphirecHitsTag_( cfg.getParameter<edm::InputTag>("rphirecHits") ),
62  matchedrecHitsTag_( cfg.getParameter<edm::InputTag>("matchedrecHits") ),
63  countStereoHitAs2D_( cfg.getParameter<bool>( "countStereoHitAs2D" ) ),
64  nHitMin2D_( cfg.getParameter<unsigned int>( "nHitMin2D" ) ),
65  // Ugly to use the same getParameter n times, but this allows const cut variables...
66  minHitsinTIB_(cfg.getParameter<edm::ParameterSet>( "minHitsPerSubDet" ).getParameter<int>( "inTIB" ) ),
67  minHitsinTOB_ (cfg.getParameter<edm::ParameterSet>( "minHitsPerSubDet" ).getParameter<int>( "inTOB" ) ),
68  minHitsinTID_ (cfg.getParameter<edm::ParameterSet>( "minHitsPerSubDet" ).getParameter<int>( "inTID" ) ),
69  minHitsinTEC_ (cfg.getParameter<edm::ParameterSet>( "minHitsPerSubDet" ).getParameter<int>( "inTEC" ) ),
70  minHitsinBPIX_ (cfg.getParameter<edm::ParameterSet>( "minHitsPerSubDet" ).getParameter<int>( "inBPIX" ) ),
71  minHitsinFPIX_ (cfg.getParameter<edm::ParameterSet>( "minHitsPerSubDet" ).getParameter<int>( "inFPIX" ) ),
72  minHitsinPIX_ (cfg.getParameter<edm::ParameterSet>( "minHitsPerSubDet" ).getParameter<int>( "inPIXEL" ) ),
73  minHitsinTIDplus_ (cfg.getParameter<edm::ParameterSet>( "minHitsPerSubDet" ).getParameter<int>( "inTIDplus" ) ),
74  minHitsinTIDminus_ (cfg.getParameter<edm::ParameterSet>( "minHitsPerSubDet" ).getParameter<int>( "inTIDminus" ) ),
75  minHitsinTECplus_ (cfg.getParameter<edm::ParameterSet>( "minHitsPerSubDet" ).getParameter<int>( "inTECplus" ) ),
76  minHitsinTECminus_ (cfg.getParameter<edm::ParameterSet>( "minHitsPerSubDet" ).getParameter<int>( "inTECminus" ) ),
77  minHitsinFPIXplus_ (cfg.getParameter<edm::ParameterSet>( "minHitsPerSubDet" ).getParameter<int>( "inFPIXplus" ) ),
78  minHitsinFPIXminus_ (cfg.getParameter<edm::ParameterSet>( "minHitsPerSubDet" ).getParameter<int>( "inFPIXminus" ) ),
79  minHitsinENDCAP_ (cfg.getParameter<edm::ParameterSet>( "minHitsPerSubDet" ).getParameter<int>( "inENDCAP" ) ),
80  minHitsinENDCAPplus_ (cfg.getParameter<edm::ParameterSet>( "minHitsPerSubDet" ).getParameter<int>( "inENDCAPplus" ) ),
81  minHitsinENDCAPminus_ (cfg.getParameter<edm::ParameterSet>( "minHitsPerSubDet" ).getParameter<int>( "inENDCAPminus" ) ),
82  maxHitDiffEndcaps_( cfg.getParameter<double>( "maxHitDiffEndcaps" ) ),
83  nLostHitMax_( cfg.getParameter<double>( "nLostHitMax" ) ),
84  clusterValueMapTag_(cfg.getParameter<edm::InputTag>("hitPrescaleMapTag")),
85  minPrescaledHits_( cfg.getParameter<int>("minPrescaledHits")),
87 {
88 
89  //convert track quality from string to enum
90  std::vector<std::string> trkQualityStrings(cfg.getParameter<std::vector<std::string> >("trackQualities"));
91  std::string qualities;
92  if(trkQualityStrings.size()>0){
94  for (unsigned int i = 0; i < trkQualityStrings.size(); ++i) {
95  (qualities += trkQualityStrings[i]) += ", ";
96  trkQualities_.push_back(reco::TrackBase::qualityByName(trkQualityStrings[i]));
97  }
98  }
99  else applyTrkQualityCheck_=false;
100 
101  std::vector<std::string> trkIterStrings(cfg.getParameter<std::vector<std::string> >("iterativeTrackingSteps"));
102  if(trkIterStrings.size()>0){
103  applyIterStepCheck_=true;
104  std::string tracksteps;
105  for (unsigned int i = 0; i < trkIterStrings.size(); ++i) {
106  (tracksteps += trkIterStrings[i]) += ", ";
107  trkSteps_.push_back(reco::TrackBase::algoByName(trkIterStrings[i]));
108  }
109  }
110  else applyIterStepCheck_=false;
111 
112  if (applyBasicCuts_){
113  edm::LogInfo("AlignmentTrackSelector")
114  << "applying basic track cuts ..."
115  << "\nptmin,ptmax: " << ptMin_ << "," << ptMax_
116  << "\npmin,pmax: " << pMin_ << "," << pMax_
117  << "\netamin,etamax: " << etaMin_ << "," << etaMax_
118  << "\nphimin,phimax: " << phiMin_ << "," << phiMax_
119  << "\nnhitmin,nhitmax: " << nHitMin_ << "," << nHitMax_
120  << "\nnlosthitmax: " << nLostHitMax_
121  << "\nnhitmin2D: " << nHitMin2D_
122  << (countStereoHitAs2D_ ? "," : ", not") << " counting hits on SiStrip stereo modules as 2D"
123  << "\nchi2nmax: " << chi2nMax_;
124 
125  if (applyIsolation_)
126  edm::LogInfo("AlignmentTrackSelector")
127  << "only retain tracks isolated at least by " << minHitIsolation_ <<
128  " cm from other rechits";
129 
130  if (chargeCheck_)
131  edm::LogInfo("AlignmentTrackSelector")
132  << "only retain hits with at least " << minHitChargeStrip_ <<
133  " ADC counts of total cluster charge";
134 
135  edm::LogInfo("AlignmentTrackSelector")
136  << "Minimum number of hits in TIB/TID/TOB/TEC/BPIX/FPIX/PIXEL = "
137  << minHitsinTIB_ << "/" << minHitsinTID_ << "/" << minHitsinTOB_
138  << "/" << minHitsinTEC_ << "/" << minHitsinBPIX_ << "/" << minHitsinFPIX_ << "/" << minHitsinPIX_;
139 
140  edm::LogInfo("AlignmentTrackSelector")
141  << "Minimum number of hits in TID+/TID-/TEC+/TEC-/FPIX+/FPIX- = "
143  << "/" << minHitsinTECplus_ << "/" << minHitsinTECminus_
144  << "/" << minHitsinFPIXplus_ << "/" << minHitsinFPIXminus_;
145 
146  edm::LogInfo("AlignmentTrackSelector")
147  << "Minimum number of hits in EndCap (TID+TEC)/EndCap+/EndCap- = "
149 
150  edm::LogInfo("AlignmentTrackSelector")
151  << "Max value of |nHitsinENDCAPplus - nHitsinENDCAPminus| = "
153 
154  if (trkQualityStrings.size()) {
155  edm::LogInfo("AlignmentTrackSelector")
156  << "Select tracks with these qualities: " << qualities;
157  }
158  }
159 
160  if (applyNHighestPt_)
161  edm::LogInfo("AlignmentTrackSelector")
162  << "filter N tracks with highest Pt N=" << nHighestPt_;
163 
165  edm::LogInfo("AlignmentTrackSelector")
166  << "apply multiplicity filter N>= " << minMultiplicity_ << "and N<= " << maxMultiplicity_
167  << " on " << (multiplicityOnInput_ ? "input" : "output");
168 
170  edm::LogInfo("AlignmentTrackSelector")
171  << "apply cut on number of prescaled hits N>= " << minPrescaledHits_
172  << " (prescale info from " << clusterValueMapTag_ << ")";
173 
174  }
175 
176 }
T getParameter(std::string const &) const
int i
Definition: DBlmapReader.cc:9
std::vector< reco::TrackBase::TrackAlgorithm > trkSteps_
std::vector< reco::TrackBase::TrackQuality > trkQualities_
const unsigned int nHitMin2D_
std::string encode() const
Definition: InputTag.cc:72
const edm::InputTag matchedrecHitsTag_
const double ptMin_
if true, cut min/maxMultiplicity on input instead of on final result
static TrackQuality qualityByName(const std::string &name)
Definition: TrackBase.cc:46
const edm::InputTag clusterValueMapTag_
const edm::InputTag rphirecHitsTag_
static TrackAlgorithm algoByName(const std::string &name)
Definition: TrackBase.cc:55
AlignmentTrackSelector::~AlignmentTrackSelector ( )

destructor

Definition at line 180 of file AlignmentTrackSelector.cc.

181 {}

Member Function Documentation

AlignmentTrackSelector::Tracks AlignmentTrackSelector::basicCuts ( const Tracks tracks,
const edm::Event evt 
) const
private

apply basic cuts on pt,eta,phi,nhit

Definition at line 230 of file AlignmentTrackSelector.cc.

References applyIterStepCheck_, applyTrkQualityCheck_, reco::TrackBase::charge(), chi2nMax_, reco::TrackBase::d0(), d0Max_, d0Min_, detailedHitsCheck(), reco::TrackBase::dz(), dzMax_, dzMin_, reco::TrackBase::eta(), eta(), etaMax_, etaMin_, funct::false, isOkTrkQuality(), nHitMax_, nHitMin_, nLostHitMax_, reco::TrackBase::normalizedChi2(), reco::TrackBase::numberOfLostHits(), reco::TrackBase::numberOfValidHits(), reco::TrackBase::p(), AlCaHLTBitMon_ParallelJobs::p, phi, reco::TrackBase::phi(), phiMax_, phiMin_, pMax_, pMin_, reco::TrackBase::pt(), ptMax_, ptMin_, lumiQueryAPI::q, query::result, theCharge_, and funct::true.

Referenced by select().

231 {
232  Tracks result;
233 
234  for (Tracks::const_iterator it=tracks.begin(); it != tracks.end(); ++it) {
235  const reco::Track* trackp=*it;
236  float pt=trackp->pt();
237  float p=trackp->p();
238  float eta=trackp->eta();
239  float phi=trackp->phi();
240  int nhit = trackp->numberOfValidHits();
241  int nlosthit = trackp->numberOfLostHits();
242  float chi2n = trackp->normalizedChi2();
243 
244  int q = trackp->charge();
245  bool isChargeOk = false;
246  if(theCharge_==-1 && q<0) isChargeOk = true;
247  else if (theCharge_==1 && q>0) isChargeOk = true;
248  else if (theCharge_==0) isChargeOk = true;
249 
250  float d0 = trackp->d0();
251  float dz = trackp->dz();
252 
253  // edm::LogDebug("AlignmentTrackSelector") << " pt,eta,phi,nhit: "
254  // <<pt<<","<<eta<<","<<phi<<","<<nhit;
255 
256  if (pt>ptMin_ && pt<ptMax_
257  && p>pMin_ && p<pMax_
258  && eta>etaMin_ && eta<etaMax_
259  && phi>phiMin_ && phi<phiMax_
260  && nhit>=nHitMin_ && nhit<=nHitMax_
261  && nlosthit<=nLostHitMax_
262  && chi2n<chi2nMax_
263  && isChargeOk
264  && d0>=d0Min_ && d0<=d0Max_
265  && dz>=dzMin_ && dz<=dzMax_) {
266  bool trkQualityOk=false ;
267  if (!applyTrkQualityCheck_ &&!applyIterStepCheck_)trkQualityOk=true ; // nothing required
268  else trkQualityOk = this->isOkTrkQuality(trackp);
269 
270  bool hitsCheckOk=this->detailedHitsCheck(trackp, evt);
271 
272  if (trkQualityOk && hitsCheckOk ) result.push_back(trackp);
273  }
274  }
275 
276  return result;
277 }
bool detailedHitsCheck(const reco::Track *track, const edm::Event &evt) const
checking hit requirements beyond simple number of valid hits
double p() const
momentum vector magnitude
Definition: TrackBase.h:129
double d0() const
dxy parameter in perigee convention (d0 = - dxy)
Definition: TrackBase.h:123
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:111
std::vector< const reco::Track * > Tracks
double phi() const
azimuthal angle of momentum vector
Definition: TrackBase.h:139
unsigned short numberOfLostHits() const
number of cases where track crossed a layer without getting a hit.
Definition: TrackBase.h:234
T eta() const
double eta() const
pseudorapidity of momentum vector
Definition: TrackBase.h:141
double pt() const
track transverse momentum
Definition: TrackBase.h:131
tuple result
Definition: query.py:137
unsigned short numberOfValidHits() const
number of valid hits found
Definition: TrackBase.h:232
bool isOkTrkQuality(const reco::Track *track) const
double dz() const
dz parameter (= dsz/cos(lambda)). This is the track z0 w.r.t (0,0,0) only if the refPoint is close to...
Definition: TrackBase.h:127
const double ptMin_
if true, cut min/maxMultiplicity on input instead of on final result
tuple tracks
Definition: testEve_cfg.py:39
int charge() const
track electric charge
Definition: TrackBase.h:113
Definition: DDAxes.h:10
AlignmentTrackSelector::Tracks AlignmentTrackSelector::checkPrescaledHits ( const Tracks tracks,
const edm::Event evt 
) const
private

Definition at line 593 of file AlignmentTrackSelector.cc.

References SiStripRecHit1D::cluster(), SiPixelRecHit::cluster(), SiStripRecHit2D::cluster(), clusterValueMapTag_, cond::rpcobgas::detid, TrackingRecHit::geographicalId(), edm::Event::getByLabel(), AlignmentClusterFlag::isTaken(), TrackingRecHit::isValid(), kBPIX, kFPIX, minPrescaledHits_, reco::Track::recHitsBegin(), reco::Track::recHitsEnd(), query::result, and DetId::subdetId().

Referenced by select().

594 {
595  Tracks result;
596 
597  //take Cluster-Flag Assomap
599  evt.getByLabel( clusterValueMapTag_, fMap);
600  const AliClusterValueMap &flagMap=*fMap;
601 
602  //for each track loop on hits and count the number of taken hits
603  for (Tracks::const_iterator ittrk=tracks.begin(); ittrk != tracks.end(); ++ittrk) {
604  const reco::Track* trackp=*ittrk;
605  int ntakenhits=0;
606  // float pt=trackp->pt();
607 
608  for (trackingRecHit_iterator ith = trackp->recHitsBegin(), edh = trackp->recHitsEnd(); ith != edh; ++ith) {
609  const TrackingRecHit *hit = ith->get(); // ith is an iterator on edm::Ref to rechit
610  if(! hit->isValid())continue;
611  DetId detid = hit->geographicalId();
612  int subDet = detid.subdetId();
614 
615  bool isPixelHit=(subDet == kFPIX || subDet == kBPIX);
616 
617  if (!isPixelHit){
618  const std::type_info &type = typeid(*hit);
619 
620  if (type == typeid(SiStripRecHit2D)) {
621  const SiStripRecHit2D* striphit=dynamic_cast<const SiStripRecHit2D*>(hit);
622  if(striphit!=0){
623  SiStripRecHit2D::ClusterRef stripclust(striphit->cluster());
624  flag = flagMap[stripclust];
625  }
626  }
627  else if(type == typeid(SiStripRecHit1D)){
628  const SiStripRecHit1D* striphit = dynamic_cast<const SiStripRecHit1D*>(hit);
629  if(striphit!=0){
630  SiStripRecHit1D::ClusterRef stripclust(striphit->cluster());
631  flag = flagMap[stripclust];
632  }
633  }
634  else{
635  edm::LogError("AlignmentTrackSelector")<<"ERROR in <AlignmentTrackSelector::checkPrescaledHits>: Dynamic cast of Strip RecHit failed!"
636  <<" Skipping this hit.";
637  continue;
638  }
639 
640 
641 
642  }//end if hit in Strips
643  else{ // test explicitely BPIX/FPIX
644  const SiPixelRecHit* pixelhit= dynamic_cast<const SiPixelRecHit*>(hit);
645  if(pixelhit!=0){
646  SiPixelRecHit::ClusterRef pixclust(pixelhit->cluster());
647  flag = flagMap[pixclust];
648  }
649  else{
650  edm::LogError("AlignmentTrackSelector")<<"ERROR in <AlignmentTrackSelector::checkPrescaledHits>: Dynamic cast of Pixel RecHit failed! ";
651  }
652  }//end else hit is in Pixel
653 
654  if(flag.isTaken())ntakenhits++;
655 
656  }//end loop on hits
657  if(ntakenhits >= minPrescaledHits_)result.push_back(trackp);
658  }//end loop on tracks
659 
660  return result;
661 }//end checkPrescaledHits
ClusterRef cluster() const
type
Definition: HCALResponse.h:22
long int flag
Definition: mlp_lapack.h:47
std::vector< const reco::Track * > Tracks
const int kFPIX
tuple result
Definition: query.py:137
ClusterRef cluster() const
trackingRecHit_iterator recHitsBegin() const
Iterator to first hit on the track.
Definition: Track.h:63
int subdetId() const
get the contents of the subdetector field (not cast into any detector&#39;s numbering enum) ...
Definition: DetId.h:39
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:356
Definition: DetId.h:20
tuple tracks
Definition: testEve_cfg.py:39
const edm::InputTag clusterValueMapTag_
ClusterRef cluster() const
Definition: SiPixelRecHit.h:41
bool isValid() const
DetId geographicalId() const
Our base class.
Definition: SiPixelRecHit.h:22
const int kBPIX
trackingRecHit_iterator recHitsEnd() const
Iterator to last hit on the track.
Definition: Track.h:65
bool AlignmentTrackSelector::detailedHitsCheck ( const reco::Track track,
const edm::Event evt 
) const
private

checking hit requirements beyond simple number of valid hits

Definition at line 281 of file AlignmentTrackSelector.cc.

References abs, applyIsolation_, chargeCheck_, isHit2D(), isIsolated(), isOkCharge(), TIDDetId::isZMinusSide(), TECDetId::isZMinusSide(), TIDDetId::isZPlusSide(), TECDetId::isZPlusSide(), kBPIX, kFPIX, maxHitDiffEndcaps_, minHitsinBPIX_, minHitsinENDCAP_, minHitsinENDCAPminus_, minHitsinENDCAPplus_, minHitsinFPIX_, minHitsinFPIXminus_, minHitsinFPIXplus_, minHitsinPIX_, minHitsinTEC_, minHitsinTECminus_, minHitsinTECplus_, minHitsinTIB_, minHitsinTID_, minHitsinTIDminus_, minHitsinTIDplus_, minHitsinTOB_, nHitMin2D_, reco::Track::recHitsBegin(), reco::Track::recHitsEnd(), seedOnlyFromAbove_, PXFDetId::side(), DetId::subdetId(), SiStripDetId::TEC, SiStripDetId::TIB, SiStripDetId::TID, SiStripDetId::TOB, and DetId::Tracker.

Referenced by basicCuts().

282 {
283  // checking hit requirements beyond simple number of valid hits
284 
292  // any detailed hit cut is active, so have to check
293 
294  int nhitinTIB = 0, nhitinTOB = 0, nhitinTID = 0;
295  int nhitinTEC = 0, nhitinBPIX = 0, nhitinFPIX = 0, nhitinPIXEL=0;
296  int nhitinENDCAP = 0, nhitinENDCAPplus = 0, nhitinENDCAPminus = 0;
297  int nhitinTIDplus = 0, nhitinTIDminus = 0;
298  int nhitinFPIXplus = 0, nhitinFPIXminus = 0;
299  int nhitinTECplus = 0, nhitinTECminus = 0;
300  unsigned int nHit2D = 0;
301  unsigned int thishit = 0;
302 
303  for (trackingRecHit_iterator iHit = trackp->recHitsBegin(); iHit != trackp->recHitsEnd(); ++iHit) {
304  thishit++;
305  const DetId detId((*iHit)->geographicalId());
306  const int subdetId = detId.subdetId();
307 
308  // *** thishit == 1 means last hit in CTF ***
309  // (FIXME: assumption might change or not be valid for all tracking algorthms)
310  // ==> for cosmics
311  // seedOnlyFrom = 1 is TIB-TOB-TEC tracks only
312  // seedOnlyFrom = 2 is TOB-TEC tracks only
313  if (seedOnlyFromAbove_ == 1 && thishit == 1
314  && (subdetId == int(SiStripDetId::TOB) || subdetId == int(SiStripDetId::TEC))){
315  return false;
316  }
317  if (seedOnlyFromAbove_ == 2 && thishit == 1 && subdetId == int(SiStripDetId::TIB)) {
318  return false;
319  }
320 
321  if (!(*iHit)->isValid()) continue; // only real hits count as in trackp->numberOfValidHits()
322  if (detId.det() != DetId::Tracker) {
323  edm::LogError("DetectorMismatch") << "@SUB=AlignmentTrackSelector::detailedHitsCheck"
324  << "DetId.det() != DetId::Tracker (=" << DetId::Tracker
325  << "), but " << detId.det() << ".";
326  }
327  const TrackingRecHit* therechit = (*iHit).get();
328  if (chargeCheck_ && !(this->isOkCharge(therechit))) return false;
329  if (applyIsolation_ && (!this->isIsolated(therechit, evt))) return false;
330  if (SiStripDetId::TIB == subdetId) ++nhitinTIB;
331  else if (SiStripDetId::TOB == subdetId) ++nhitinTOB;
332  else if (SiStripDetId::TID == subdetId) {
333  ++nhitinTID;
334  ++nhitinENDCAP;
335  TIDDetId tidId(detId);
336  if (tidId.isZMinusSide()) {
337  ++nhitinTIDminus;
338  ++nhitinENDCAPminus;
339  }
340  else if (tidId.isZPlusSide()) {
341  ++nhitinTIDplus;
342  ++nhitinENDCAPplus;
343  }
344  }
345  else if (SiStripDetId::TEC == subdetId) {
346  ++nhitinTEC;
347  ++nhitinENDCAP;
348  TECDetId tecId(detId);
349  if (tecId.isZMinusSide()) {
350  ++nhitinTECminus;
351  ++nhitinENDCAPminus;
352  }
353  else if (tecId.isZPlusSide()) {
354  ++nhitinTECplus;
355  ++nhitinENDCAPplus;
356  }
357  }
358  else if ( kBPIX == subdetId) {++nhitinBPIX;++nhitinPIXEL;}
359  else if ( kFPIX == subdetId) {
360  ++nhitinFPIX;
361  ++nhitinPIXEL;
362  PXFDetId fpixId(detId);
363  if (fpixId.side()==1) ++nhitinFPIXminus;
364  else if (fpixId.side()==2) ++nhitinFPIXplus;
365  }
366  // Do not call isHit2D(..) if already enough 2D hits for performance reason:
367  if (nHit2D < nHitMin2D_ && this->isHit2D(**iHit)) ++nHit2D;
368  } // end loop on hits
369 
370  return (nhitinTIB >= minHitsinTIB_ && nhitinTOB >= minHitsinTOB_
371  && nhitinTID >= minHitsinTID_ && nhitinTEC >= minHitsinTEC_
372  && nhitinENDCAP >= minHitsinENDCAP_ && nhitinENDCAPplus >= minHitsinENDCAPplus_ && nhitinENDCAPminus >= minHitsinENDCAPminus_
373  && std::abs(nhitinENDCAPplus-nhitinENDCAPminus) <= maxHitDiffEndcaps_
374  && nhitinTIDplus >= minHitsinTIDplus_ && nhitinTIDminus >= minHitsinTIDminus_
375  && nhitinFPIXplus >= minHitsinFPIXplus_ && nhitinFPIXminus >= minHitsinFPIXminus_
376  && nhitinTECplus >= minHitsinTECplus_ && nhitinTECminus >= minHitsinTECminus_
377  && nhitinBPIX >= minHitsinBPIX_
378  && nhitinFPIX >= minHitsinFPIX_ && nhitinPIXEL>=minHitsinPIX_
379  && nHit2D >= nHitMin2D_);
380  } else { // no cuts set, so we are just fine and can avoid loop on hits
381  return true;
382  }
383 
384 }
bool isIsolated(const TrackingRecHit *therechit, const edm::Event &evt) const
#define abs(x)
Definition: mlp_lapack.h:159
const unsigned int nHitMin2D_
const int kFPIX
bool isHit2D(const TrackingRecHit &hit) const
int subdetId() const
get the contents of the subdetector field (not cast into any detector&#39;s numbering enum) ...
Definition: DetId.h:39
Definition: DetId.h:20
bool isOkCharge(const TrackingRecHit *therechit) const
if valid, check for minimum charge (currently only in strip), if invalid give true ...
const int kBPIX
bool AlignmentTrackSelector::isHit2D ( const TrackingRecHit hit) const
private

Definition at line 388 of file AlignmentTrackSelector.cc.

References countStereoHitAs2D_, TrackingRecHit::dimension(), TrackingRecHit::geographicalId(), TrackingRecHit::isValid(), kBPIX, kFPIX, ProjectedSiStripRecHit2D::originalHit(), SiStripDetId::stereo(), and DetId::Tracker.

Referenced by detailedHitsCheck().

389 {
390  // we count SiStrip stereo modules as 2D if selected via countStereoHitAs2D_
391  // (since they provide theta information)
392  if (!hit.isValid() ||
393  (hit.dimension() < 2 && !countStereoHitAs2D_ && !dynamic_cast<const SiStripRecHit1D*>(&hit))){
394  return false; // real RecHit1D - but SiStripRecHit1D depends on countStereoHitAs2D_
395  } else {
396  const DetId detId(hit.geographicalId());
397  if (detId.det() == DetId::Tracker) {
398  if (detId.subdetId() == kBPIX || detId.subdetId() == kFPIX) {
399  return true; // pixel is always 2D
400  } else { // should be SiStrip now
401  const SiStripDetId stripId(detId);
402  if (stripId.stereo()) return countStereoHitAs2D_; // stereo modules
403  else if (dynamic_cast<const SiStripRecHit1D*>(&hit)
404  || dynamic_cast<const SiStripRecHit2D*>(&hit)) return false; // rphi modules hit
405  //the following two are not used any more since ages...
406  else if (dynamic_cast<const SiStripMatchedRecHit2D*>(&hit)) return true; // matched is 2D
407  else if (dynamic_cast<const ProjectedSiStripRecHit2D*>(&hit)) {
408  const ProjectedSiStripRecHit2D* pH = static_cast<const ProjectedSiStripRecHit2D*>(&hit);
409  return (countStereoHitAs2D_ && this->isHit2D(pH->originalHit())); // depends on original...
410  } else {
411  edm::LogError("UnkownType") << "@SUB=AlignmentTrackSelector::isHit2D"
412  << "Tracker hit not in pixel, neither SiStripRecHit[12]D nor "
413  << "SiStripMatchedRecHit2D nor ProjectedSiStripRecHit2D.";
414  return false;
415  }
416  }
417  } else { // not tracker??
418  edm::LogWarning("DetectorMismatch") << "@SUB=AlignmentTrackSelector::isHit2D"
419  << "Hit not in tracker with 'official' dimension >=2.";
420  return true; // dimension() >= 2 so accept that...
421  }
422  }
423  // never reached...
424 }
virtual int dimension() const =0
const int kFPIX
bool isHit2D(const TrackingRecHit &hit) const
Detector identifier class for the strip tracker.
Definition: SiStripDetId.h:17
Definition: DetId.h:20
bool isValid() const
DetId geographicalId() const
const SiStripRecHit2D & originalHit() const
const int kBPIX
bool AlignmentTrackSelector::isIsolated ( const TrackingRecHit therechit,
const edm::Event evt 
) const
private

Definition at line 526 of file AlignmentTrackSelector.cc.

References TrackingRecHit::geographicalId(), edm::Event::getByLabel(), BaseTrackerRecHit::localPosition(), TrackingRecHit::localPosition(), mag(), matchedrecHitsTag_, minHitIsolation_, DetId::rawId(), and rphirecHitsTag_.

Referenced by detailedHitsCheck().

527 {
528  // FIXME:
529  // adapt to changes after introduction of SiStripRecHit1D...
530  //
531  // edm::ESHandle<TrackerGeometry> tracker;
534  // es.get<TrackerDigiGeometryRecord>().get(tracker);
535  evt.getByLabel( rphirecHitsTag_, rphirecHits );
536  evt.getByLabel( matchedrecHitsTag_, matchedrecHits );
537 
540  const SiStripRecHit2DCollection::DataContainer& stripcollSt = rphirecHits->data();
541  const SiStripMatchedRecHit2DCollection::DataContainer& stripcollStm = matchedrecHits->data();
542 
543  DetId idet = therechit->geographicalId();
544 
545  // FIXME: instead of looping the full hit collection, we should explore the features of
546  // SiStripRecHit2DCollection::rangeRphi = rphirecHits.get(idet) and loop
547  // only from rangeRphi.first until rangeRphi.second
548  for( istripSt=stripcollSt.begin(); istripSt!=stripcollSt.end(); ++istripSt ) {
549  const SiStripRecHit2D *aHit = &*(istripSt);
550  DetId mydet1 = aHit->geographicalId();
551  if (idet.rawId() != mydet1.rawId()) continue;
552  float theDistance = ( therechit->localPosition() - aHit->localPosition() ).mag();
553  // std::cout << "theDistance1 = " << theDistance << "\n";
554  if (theDistance > 0.001 && theDistance < minHitIsolation_) return false;
555  }
556 
557  // FIXME: see above
558  for( istripStm=stripcollStm.begin(); istripStm!=stripcollStm.end(); ++istripStm ) {
559  const SiStripMatchedRecHit2D *aHit = &*(istripStm);
560  DetId mydet2 = aHit->geographicalId();
561  if (idet.rawId() != mydet2.rawId()) continue;
562  float theDistance = (therechit->localPosition() - aHit->localPosition()).mag();
563  // std::cout << "theDistance1 = " << theDistance << "\n";
564  if (theDistance > 0.001 && theDistance < minHitIsolation_) return false;
565  }
566  return true;
567 }
boost::transform_iterator< IterHelp, const_IdIter > const_iterator
std::vector< data_type > DataContainer
T mag() const
The vector magnitude. Equivalent to sqrt(vec.mag2())
uint32_t rawId() const
get the raw id
Definition: DetId.h:45
const edm::InputTag matchedrecHitsTag_
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:356
Definition: DetId.h:20
virtual LocalPoint localPosition() const
const edm::InputTag rphirecHitsTag_
DetId geographicalId() const
virtual LocalPoint localPosition() const =0
bool AlignmentTrackSelector::isOkCharge ( const TrackingRecHit therechit) const
private

if valid, check for minimum charge (currently only in strip), if invalid give true

Definition at line 428 of file AlignmentTrackSelector.cc.

References TrackingRecHit::geographicalId(), errorMatrix2Lands_multiChannel::id, isOkChargeStripHit(), TrackingRecHit::isValid(), kBPIX, kFPIX, SiStripMatchedRecHit2D::monoHit(), ProjectedSiStripRecHit2D::originalHit(), SiStripMatchedRecHit2D::stereoHit(), and DetId::Tracker.

Referenced by detailedHitsCheck().

429 {
430  if (!hit || !hit->isValid()) return true;
431 
432  // check det and subdet
433  const DetId id(hit->geographicalId());
434  if (id.det() != DetId::Tracker) {
435  edm::LogWarning("DetectorMismatch") << "@SUB=isOkCharge" << "Hit not in tracker!";
436  return true;
437  }
438  if (id.subdetId() == kFPIX || id.subdetId() == kBPIX) {
439  return true; // might add some requirement...
440  }
441 
442  // We are in SiStrip now, so test normal hit:
443  const std::type_info &type = typeid(*hit);
444 
445 
446  if (type == typeid(SiStripRecHit2D)) {
447  const SiStripRecHit2D *stripHit2D = dynamic_cast<const SiStripRecHit2D*>(hit);
448  if (stripHit2D) {
449  return this->isOkChargeStripHit(*stripHit2D);
450  }
451  }
452  else if(type == typeid(SiStripRecHit1D)){
453  const SiStripRecHit1D *stripHit1D = dynamic_cast<const SiStripRecHit1D*>(hit);
454  if (stripHit1D) {
455  return this->isOkChargeStripHit(*stripHit1D);
456  }
457  }
458 
459  else if(type == typeid(SiStripMatchedRecHit2D)){ // or matched (should not occur anymore due to hit splitting since 20X)
460  const SiStripMatchedRecHit2D *matchedHit = dynamic_cast<const SiStripMatchedRecHit2D*>(hit);
461  if (matchedHit) {
462  return (this->isOkChargeStripHit(matchedHit->monoHit())
463  && this->isOkChargeStripHit(matchedHit->stereoHit()));
464  }
465  }
466  else if(type == typeid(ProjectedSiStripRecHit2D)){
467  // or projected (should not occur anymore due to hit splitting since 20X):
468  const ProjectedSiStripRecHit2D *projHit = dynamic_cast<const ProjectedSiStripRecHit2D*>(hit);
469  if (projHit) {
470  return this->isOkChargeStripHit(projHit->originalHit());
471  }
472  }
473  else{
474  edm::LogError("AlignmentTrackSelector")<< "@SUB=isOkCharge" <<"Unknown type of a valid tracker hit in Strips "
475  << " SubDet = "<<id.subdetId();
476  return false;
477  }
478 
479 
480 
481  // and now? SiTrackerMultiRecHit? Not here I guess!
482  // Should we throw instead?
483  edm::LogError("AlignmentTrackSelector")
484  << "@SUB=isOkCharge" << "Unknown valid tracker hit not in pixel, subdet " << id.subdetId()
485  << ", SiTrackerMultiRecHit " << dynamic_cast<const SiTrackerMultiRecHit*>(hit)
486  << ", BaseTrackerRecHit " << dynamic_cast<const BaseTrackerRecHit*>(hit);
487 
488  return true;
489 }
type
Definition: HCALResponse.h:22
bool isOkChargeStripHit(const SiStripRecHit1D &siStripRecHit1D) const
const int kFPIX
Definition: DetId.h:20
SiStripRecHit2D stereoHit() const
SiStripRecHit2D monoHit() const
const SiStripRecHit2D & originalHit() const
const int kBPIX
bool AlignmentTrackSelector::isOkChargeStripHit ( const SiStripRecHit1D siStripRecHit1D) const
private

Definition at line 510 of file AlignmentTrackSelector.cc.

References DeDxDiscriminatorTools::charge(), SiStripRecHit1D::cluster(), and minHitChargeStrip_.

Referenced by isOkCharge().

511 {
512  double charge = 0.;
513 
514  SiStripRecHit1D::ClusterRef cluster(siStripRecHit1D.cluster());
515  const std::vector<uint8_t> &amplitudes = cluster->amplitudes();
516 
517  for (size_t ia = 0; ia < amplitudes.size(); ++ia) {
518  charge += amplitudes[ia];
519  }
520 
521  return (charge >= minHitChargeStrip_);
522 }
ClusterRef cluster() const
double charge(const std::vector< uint8_t > &Ampls)
bool AlignmentTrackSelector::isOkChargeStripHit ( const SiStripRecHit2D siStripRecHit2D) const
private

Definition at line 494 of file AlignmentTrackSelector.cc.

References DeDxDiscriminatorTools::charge(), SiStripRecHit2D::cluster(), and minHitChargeStrip_.

495 {
496  double charge = 0.;
497 
498  SiStripRecHit2D::ClusterRef cluster(siStripRecHit2D.cluster());
499  const std::vector<uint8_t> &amplitudes = cluster->amplitudes();
500 
501  for (size_t ia = 0; ia < amplitudes.size(); ++ia) {
502  charge += amplitudes[ia];
503  }
504 
505  return (charge >= minHitChargeStrip_);
506 }
double charge(const std::vector< uint8_t > &Ampls)
ClusterRef cluster() const
bool AlignmentTrackSelector::isOkTrkQuality ( const reco::Track track) const
private

Definition at line 665 of file AlignmentTrackSelector.cc.

References reco::TrackBase::algo(), applyIterStepCheck_, applyTrkQualityCheck_, i, reco::TrackBase::quality(), trkQualities_, and trkSteps_.

Referenced by basicCuts().

666 {
667  bool qualityOk=false;
668  bool iterStepOk=false;
669 
670  //check iterative step
672  for (unsigned int i = 0; i < trkSteps_.size(); ++i) {
673  if (track->algo()==(trkSteps_[i])) {
674  iterStepOk=true;
675  }
676  }
677  }
678  else iterStepOk=true;
679 
680  //check track quality
682  for (unsigned int i = 0; i < trkQualities_.size(); ++i) {
683  if (track->quality(trkQualities_[i])) {
684  qualityOk=true;
685  }
686  }
687  }
688  else qualityOk=true;
689 
690  return qualityOk&&iterStepOk;
691 }//end check on track quality
int i
Definition: DBlmapReader.cc:9
std::vector< reco::TrackBase::TrackAlgorithm > trkSteps_
std::vector< reco::TrackBase::TrackQuality > trkQualities_
TrackAlgorithm algo() const
Definition: TrackBase.h:332
bool quality(const TrackQuality) const
Track quality.
Definition: TrackBase.h:377
AlignmentTrackSelector::Tracks AlignmentTrackSelector::select ( const Tracks tracks,
const edm::Event evt 
) const

select tracks

Definition at line 187 of file AlignmentTrackSelector.cc.

References applyBasicCuts_, applyMultiplicityFilter_, applyNHighestPt_, applyPrescaledHitsFilter_, basicCuts(), checkPrescaledHits(), maxMultiplicity_, minMultiplicity_, multiplicityOnInput_, query::result, theNHighestPtTracks(), and testEve_cfg::tracks.

Referenced by Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::find(), Vispa.Views.LineDecayView.LineDecayContainer::mousePressEvent(), Vispa.Gui.PortConnection.PointToPointConnection::mousePressEvent(), Vispa.Gui.VispaWidget.VispaWidget::mousePressEvent(), Vispa.Views.AbstractView.AbstractView::restoreSelection(), and TrackConfigSelector::select().

188 {
189 
191  (tracks.size() < static_cast<unsigned int>(minMultiplicity_)
192  || tracks.size() > static_cast<unsigned int>(maxMultiplicity_))) {
193  return Tracks(); // empty collection
194  }
195 
196  Tracks result = tracks;
197  // apply basic track cuts (if selected)
198  if (applyBasicCuts_) result= this->basicCuts(result, evt);
199 
200  // filter N tracks with highest Pt (if selected)
201  if (applyNHighestPt_) result = this->theNHighestPtTracks(result);
202 
203  // apply minimum multiplicity requirement (if selected)
205  if (result.size() < static_cast<unsigned int>(minMultiplicity_)
206  || result.size() > static_cast<unsigned int>(maxMultiplicity_) ) {
207 
208  result.clear();
209  }
210  }
211 
213  result = this->checkPrescaledHits(result, evt);
214  }
215 
216  return result;
217 }
Tracks basicCuts(const Tracks &tracks, const edm::Event &evt) const
apply basic cuts on pt,eta,phi,nhit
std::vector< const reco::Track * > Tracks
Tracks checkPrescaledHits(const Tracks &tracks, const edm::Event &evt) const
Tracks theNHighestPtTracks(const Tracks &tracks) const
filter the n highest pt tracks
tuple result
Definition: query.py:137
tuple tracks
Definition: testEve_cfg.py:39
AlignmentTrackSelector::Tracks AlignmentTrackSelector::theNHighestPtTracks ( const Tracks tracks) const
private

filter the n highest pt tracks

Definition at line 572 of file AlignmentTrackSelector.cc.

References n, nHighestPt_, ptComparator, query::result, python.multivaluedict::sort(), and testEve_cfg::tracks.

Referenced by select().

573 {
574  Tracks sortedTracks=tracks;
575  Tracks result;
576 
577  // sort in pt
578  std::sort(sortedTracks.begin(),sortedTracks.end(),ptComparator);
579 
580  // copy theTrackMult highest pt tracks to result vector
581  int n=0;
582  for (Tracks::const_iterator it=sortedTracks.begin();
583  it!=sortedTracks.end(); ++it) {
584  if (n<nHighestPt_) { result.push_back(*it); n++; }
585  }
586 
587  return result;
588 }
std::vector< const reco::Track * > Tracks
tuple result
Definition: query.py:137
tuple tracks
Definition: testEve_cfg.py:39
bool AlignmentTrackSelector::useThisFilter ( )

Member Data Documentation

const bool AlignmentTrackSelector::applyBasicCuts_
private

Definition at line 64 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), select(), and useThisFilter().

const bool AlignmentTrackSelector::applyIsolation_
private

Definition at line 65 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), and detailedHitsCheck().

bool AlignmentTrackSelector::applyIterStepCheck_
private

Definition at line 92 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), basicCuts(), and isOkTrkQuality().

const bool AlignmentTrackSelector::applyMultiplicityFilter_
private

Definition at line 64 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), select(), and useThisFilter().

const bool AlignmentTrackSelector::applyNHighestPt_
private

Definition at line 64 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), select(), and useThisFilter().

const bool AlignmentTrackSelector::applyPrescaledHitsFilter_
private

Definition at line 86 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), select(), and useThisFilter().

bool AlignmentTrackSelector::applyTrkQualityCheck_
private

Definition at line 91 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), basicCuts(), and isOkTrkQuality().

const bool AlignmentTrackSelector::chargeCheck_
private

Definition at line 65 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), and detailedHitsCheck().

const double AlignmentTrackSelector::chi2nMax_
private

Definition at line 69 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), and basicCuts().

const edm::InputTag AlignmentTrackSelector::clusterValueMapTag_
private

Definition at line 84 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), and checkPrescaledHits().

const bool AlignmentTrackSelector::countStereoHitAs2D_
private

Definition at line 74 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), and isHit2D().

const double AlignmentTrackSelector::d0Max_
private

Definition at line 69 of file AlignmentTrackSelector.h.

Referenced by basicCuts().

const double AlignmentTrackSelector::d0Min_
private

Definition at line 69 of file AlignmentTrackSelector.h.

Referenced by basicCuts().

const double AlignmentTrackSelector::dzMax_
private

Definition at line 69 of file AlignmentTrackSelector.h.

Referenced by basicCuts().

const double AlignmentTrackSelector::dzMin_
private

Definition at line 69 of file AlignmentTrackSelector.h.

Referenced by basicCuts().

const double AlignmentTrackSelector::etaMax_
private

Definition at line 68 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), and basicCuts().

const double AlignmentTrackSelector::etaMin_
private

Definition at line 68 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), and basicCuts().

const edm::InputTag AlignmentTrackSelector::matchedrecHitsTag_
private

Definition at line 73 of file AlignmentTrackSelector.h.

Referenced by isIsolated().

const double AlignmentTrackSelector::maxHitDiffEndcaps_
private

Definition at line 81 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), and detailedHitsCheck().

const int AlignmentTrackSelector::maxMultiplicity_
private

Definition at line 66 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), and select().

const double AlignmentTrackSelector::minHitChargeStrip_
private

Definition at line 71 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), and isOkChargeStripHit().

const double AlignmentTrackSelector::minHitIsolation_
private

Definition at line 71 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), and isIsolated().

const int AlignmentTrackSelector::minHitsinBPIX_
private

Definition at line 77 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), and detailedHitsCheck().

const int AlignmentTrackSelector::minHitsinENDCAP_
private

Definition at line 80 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), and detailedHitsCheck().

const int AlignmentTrackSelector::minHitsinENDCAPminus_
private

Definition at line 80 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), and detailedHitsCheck().

const int AlignmentTrackSelector::minHitsinENDCAPplus_
private

Definition at line 80 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), and detailedHitsCheck().

const int AlignmentTrackSelector::minHitsinFPIX_
private

Definition at line 77 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), and detailedHitsCheck().

const int AlignmentTrackSelector::minHitsinFPIXminus_
private

Definition at line 79 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), and detailedHitsCheck().

const int AlignmentTrackSelector::minHitsinFPIXplus_
private

Definition at line 79 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), and detailedHitsCheck().

const int AlignmentTrackSelector::minHitsinPIX_
private

Definition at line 77 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), and detailedHitsCheck().

const int AlignmentTrackSelector::minHitsinTEC_
private

Definition at line 76 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), and detailedHitsCheck().

const int AlignmentTrackSelector::minHitsinTECminus_
private

Definition at line 78 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), and detailedHitsCheck().

const int AlignmentTrackSelector::minHitsinTECplus_
private

Definition at line 78 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), and detailedHitsCheck().

const int AlignmentTrackSelector::minHitsinTIB_
private

Definition at line 76 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), and detailedHitsCheck().

const int AlignmentTrackSelector::minHitsinTID_
private

Definition at line 76 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), and detailedHitsCheck().

const int AlignmentTrackSelector::minHitsinTIDminus_
private

Definition at line 78 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), and detailedHitsCheck().

const int AlignmentTrackSelector::minHitsinTIDplus_
private

Definition at line 78 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), and detailedHitsCheck().

const int AlignmentTrackSelector::minHitsinTOB_
private

Definition at line 76 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), and detailedHitsCheck().

const int AlignmentTrackSelector::minMultiplicity_
private

Definition at line 66 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), and select().

const int AlignmentTrackSelector::minPrescaledHits_
private

Definition at line 85 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), and checkPrescaledHits().

const bool AlignmentTrackSelector::multiplicityOnInput_
private

Definition at line 67 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), and select().

const int AlignmentTrackSelector::nHighestPt_
private

Definition at line 66 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), and theNHighestPtTracks().

const double AlignmentTrackSelector::nHitMax_
private

Definition at line 69 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), and basicCuts().

const unsigned int AlignmentTrackSelector::nHitMin2D_
private

Definition at line 75 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), and detailedHitsCheck().

const double AlignmentTrackSelector::nHitMin_
private

Definition at line 69 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), and basicCuts().

const double AlignmentTrackSelector::nLostHitMax_
private

Definition at line 82 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), and basicCuts().

const double AlignmentTrackSelector::phiMax_
private

Definition at line 68 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), and basicCuts().

const double AlignmentTrackSelector::phiMin_
private

Definition at line 68 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), and basicCuts().

const double AlignmentTrackSelector::pMax_
private

Definition at line 68 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), and basicCuts().

const double AlignmentTrackSelector::pMin_
private

Definition at line 68 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), and basicCuts().

ComparePt AlignmentTrackSelector::ptComparator
private

Definition at line 62 of file AlignmentTrackSelector.h.

Referenced by theNHighestPtTracks().

const double AlignmentTrackSelector::ptMax_
private

Definition at line 68 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), and basicCuts().

const double AlignmentTrackSelector::ptMin_
private

if true, cut min/maxMultiplicity on input instead of on final result

Definition at line 68 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), and basicCuts().

const edm::InputTag AlignmentTrackSelector::rphirecHitsTag_
private

Definition at line 72 of file AlignmentTrackSelector.h.

Referenced by isIsolated().

const bool AlignmentTrackSelector::seedOnlyFromAbove_
private

Definition at line 65 of file AlignmentTrackSelector.h.

Referenced by detailedHitsCheck().

const int AlignmentTrackSelector::theCharge_
private

Definition at line 70 of file AlignmentTrackSelector.h.

Referenced by basicCuts().

std::vector<reco::TrackBase::TrackQuality> AlignmentTrackSelector::trkQualities_
private

Definition at line 88 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), and isOkTrkQuality().

std::vector<reco::TrackBase::TrackAlgorithm> AlignmentTrackSelector::trkSteps_
private

Definition at line 90 of file AlignmentTrackSelector.h.

Referenced by AlignmentTrackSelector(), and isOkTrkQuality().