120 {
121
122 static int hpmap[20] = { 0, 1, 2, 7, 8, 9, 3, 4, 5, 6, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19 };
123
126
127 while(muonitr!=muonend) {
128
129 bool hasInDet = (*muonitr)->hasInDetTrackParticle();
130 const Rec::TrackParticle *trackitr=NULL;
131 if(hasInDet) trackitr=(*muonitr)->inDetTrackParticle();
132 muonitr++;
133 if(!hasInDet) continue;
134
136
137 CLHEP::HepVector perigeeParams = measPer->parameters();
138 double pT = measPer->
pT();
139 double eta = measPer->
eta();
143
144
145 double deta = 1;
146 double dphi = 1;
147 double dz0 = 1;
148 double dd0 = 1;
149 double dpT = 1;
150
151
152
153
154
161
163
164 const Trk::FitQuality *quality = trackitr->
fitQuality();
166 double dof = quality->quality->
numberDoF();
167
168 unsigned bitmap = 0;
169
170 for ( int ih=0 ; ih<20 ; ih++ ) {
172 }
173
174
176
177
178 TIDA::Track*
t =
new TIDA::Track(
eta,
phi, z0, d0, pT,
chi2, dof,
179 deta, dphi, dz0, dd0, dpT,
180 nBlayerHits, nPixelHits, nSctHits, nSiHits,
181 nStrawHits, nTrHits, bitmap, 0,
182 -1) ;
183
185 }
186
187 }
void addTrack(TIDA::Track *t)
void ipCorr(double d0, double z0, double &d0c, double &z0c, double phi0, double eta, double pt)
DataModel_detail::const_iterator< DataVector > const_iterator
const_iterator end() const noexcept
Return a const_iterator pointing past the end of the collection.
const_iterator begin() const noexcept
Return a const_iterator pointing at the beginning of the collection.
const Trk::Perigee * measuredPerigee() const
Accessor method for Perigee.
int numberDoF() const
returns the number of degrees of freedom of the overall track or vertex fit as integer
double chiSquared() const
returns the of the overall track fit
double eta() const
Access method for pseudorapidity - from momentum.
double pT() const
Access method for transverse momentum.
const FitQuality * fitQuality() const
accessor function for FitQuality.
const TrackSummary * trackSummary() const
accessor function for TrackSummary.
double chi2(TH1 *h0, TH1 *h1)
float nSiHits(const U &p)
float nPixelHits(const U &p)
DetectorType
enumerates the various detector types currently accessible from the isHit() method.
@ numberOfSCTHits
number of SCT holes
@ numberOfPixelHits
number of pixel layers on track with absence of hits
@ numberOfBLayerHits
these are the hits in the 0th pixel layer?
@ numberOfTRTHits
number of TRT outliers
@ numberOfTRTHighThresholdHits
total number of TRT hits which pass the high threshold