ATLAS Offline Software
Loading...
Searching...
No Matches
TrigTrackSelector.cxx
Go to the documentation of this file.
1
9
10
11#include "TrigInDetAnalysisUtils/TrigTrackSelector.h"
12#include "TruthUtils/HepMCHelpers.h"
13
14#include "xAODTruth/TruthVertexContainer.h"
15
16#include <stdexcept>
17
21const double TrigTrackSelector::s_default_radius = 47;
22
23
24TrigTrackSelector::TrigTrackSelector( TrackFilter* selector, double radius, int selectPdgId, int selectParentPdgId ) :
25 TrackSelector(selector), m_id(0), m_xBeam(0), m_yBeam(0), m_zBeam(0),
26 m_correctTrkTracks(false),
27 m_radius(radius), m_selectPdgId(selectPdgId), m_selectParentPdgId(selectParentPdgId) {
28}
29
30
31
38
39const xAOD::TruthParticle* TrigTrackSelector::fromAncestor(const int pdg_id, const xAOD::TruthParticle *p) const {
40 if ( p==nullptr ) return nullptr;
41 if ( MC::isElectron(p) || MC::isMuon(p) ) return nullptr; //don't want light leptons from eg tau decays - they are found directly
42 if ( p->absPdgId()==pdg_id ) {
43 return p; // recursive stopping conditions
44 }
45 auto vertex = p->prodVtx();
46 if ( vertex == nullptr ) {
47 return nullptr; // has no production vertex !!!
48 }
49 if ( vertex->nIncomingParticles() < 1 ) {
50 return nullptr; // recursive stopping conditions
51 }
52 for( unsigned ip = 0; ip < vertex->nIncomingParticles(); ip++ ) {
53 auto* in = vertex->incomingParticle(ip);
54 auto parent = fromAncestor( pdg_id, in);
55 if ( parent!=nullptr ) {
56 if (parent->absPdgId()==pdg_id) return parent;
57 }
58 }
59
60 return nullptr;
61}
62
63
72
73const xAOD::TruthParticle* TrigTrackSelector::fromAncestor( const std::vector<int>& ids, const xAOD::TruthParticle *p ) const {
74 if ( p==nullptr ) return nullptr;
75 if (MC::isElectron(p) || MC::isMuon(p)) return nullptr; //don't want light leptons from eg tau decays - they are found directly
76 for ( size_t i=ids.size() ; i-- ; ) {
77 if ( p->absPdgId()==ids[i] ) return p; // recursive stopping conditions
78 }
79
80 auto vertex = p->prodVtx();
81 if ( vertex == nullptr ) return nullptr; // has no production vertex !!!
82
83 if ( vertex->nIncomingParticles()<1 ) return nullptr; // recursive stopping conditions
84
85 for( unsigned ip = 0; ip < vertex->nIncomingParticles(); ip++ ) {
86 auto* in = vertex->incomingParticle(ip);
87 auto parent = fromAncestor( ids, in);
88 if ( parent!=nullptr ) {
89 for ( size_t i=ids.size() ; i-- ; ) {
90 if ( parent->absPdgId()==ids[i] ) return parent;
91 }
92 }
93 }
94
95 return nullptr;
96}
97
98
99
106
107// const xAOD::TruthParticle* TrigTrackSelector::fromAncestor(const int pdg_id, const xAOD::TruthParticle *p) const {
108// return fromAncestor( std::vector<int>(1,pdg_id), p );
109// }
110
111
112
113
114
115
116bool TrigTrackSelector::selectTrack( const TrigInDetTrack* track, const TrigInDetTrackTruthMap* truthMap ) {
117 // do the track extraction stuff here....
118 if ( track ) {
119
120 double eta = track->param()->eta();
121 double phi = track->param()->phi0();
122 double z0 = track->param()->z0();
123 double pT = track->param()->pT();
124 double d0 = track->param()->a0();
125
126 double deta = track->param()->eeta();
127 double dphi = track->param()->ephi0();
128 double dz0 = track->param()->ez0();
129 double dpT = track->param()->epT();
130 double dd0 = track->param()->ea0();
131
132 double theta = 2*std::atan( exp( (-1)*eta ) );
133 correctToBeamline( z0, dz0, d0, dd0, theta, phi );
134
135 int algoid = track->algorithmId();
136
137 int nBlayerHits = (track->HitPattern() & 0x1);
138 int nPixelHits = 2 * track->NPixelSpacePoints(); // NB: for comparison with offline
139 int nSctHits = 2 * track->NSCT_SpacePoints(); // a spacepoint is 2 "hits"
140 int nStrawHits = track->NStrawHits();
141 int nTrHits = track->NTRHits();
142
143 int nSiHits = nPixelHits + nSctHits;
144
145 bool expectBL = false; //not filled in
146
147 unsigned long id = (unsigned long)track;
148
149 unsigned hitPattern = track->HitPattern();
150 unsigned multiPattern = 0;
151
152 double chi2 = track->chi2();
153 double dof = 0;
154
155 bool truth = false;
156 int match_uniqueID = HepMC::INVALID_PARTICLE_ID;
157
158 if ( truthMap ) {
159 const TrigInDetTrackTruth* trackTruth = truthMap->truth(track);
160 if (trackTruth!=0 && trackTruth->nrMatches() > 0) {
161 match_uniqueID = HepMC::uniqueID(trackTruth->bestSiMatch());
162 truth = true;
163 }
164 }
165
166
167 TIDA::Track* t = new TIDA::Track( eta, phi, z0, d0, pT, chi2, dof,
168 deta, dphi, dz0, dd0, dpT,
169 nBlayerHits, nPixelHits, nSctHits, nSiHits,
170 nStrawHits, nTrHits,
171 hitPattern, multiPattern,
172 algoid, truth, -1, match_uniqueID,
173 expectBL, id) ;
174
175 // std::cout << "SUTT ID track " << *t << "\t0x" << std::hex << track->HitPattern() << std::dec << std::endl;
176
177 if ( !addTrack( t ) ){
178 delete t;
179 return false;
180 }
181 return true;
182 }
183 return false;
184}
185
186
187// extract all the tracks from a TrigInDetTrack collection and associated TruthMap and convert them
188void TrigTrackSelector::selectTracks( const TrigInDetTrackCollection* trigtracks, const TrigInDetTrackTruthMap* truthMap ) {
189 // do the track extraction stuff here....
190 TrigInDetTrackCollection::const_iterator trackitr = trigtracks->begin();
191 TrigInDetTrackCollection::const_iterator trackend = trigtracks->end();
192 while ( trackitr!=trackend ) {
193 selectTrack( *trackitr, truthMap );
194 ++trackitr;
195 }
196}
197
198
199// add a TrackParticle
200bool TrigTrackSelector::selectTrack( const Rec::TrackParticle* track ) {
201
202 // do the track extraction stuff here....
203
204 static const int hpmap[20] = { 0, 1, 2, 7, 8, 9, 3, 4, 5, 6, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19 };
205
206 if ( track ) {
207
208#ifdef TRKPARAMETERS_MEASUREDPERIGEE_H
209 const Trk::MeasuredPerigee* measPer = track->measuredPerigee();
210#else
211 const Trk::Perigee* measPer = track->measuredPerigee();
212#endif
213 // CLHEP::HepVector perigeeParams = measPer->parameters();
214
215 double pT = measPer->pT();
216 double eta = measPer->eta();
217 double phi = measPer->parameters()[Trk::phi0];
218 double z0 = measPer->parameters()[Trk::z0] + m_zBeam;
219 double d0 = measPer->parameters()[Trk::d0];
220
221 double theta = measPer->parameters()[Trk::theta];
222 double p = 1/measPer->parameters()[Trk::qOverP];
223
224 // AAARCH!!!!! the TrackParticle pT is NOT SIGNED!!!! ( I ask you! )
225 if ( measPer->parameters()[Trk::qOverP]<0 && pT>0 ) pT *= -1;
226
227#ifdef TRKPARAMETERS_MEASUREDPERIGEE_H
228 const Trk::ErrorMatrix err = measPer->localErrorMatrix();
229 double dtheta = err.error(Trk::theta);
230 double dqovp = err.error(Trk::qOverP);
231 double covthetaOvP = err.covValue(Trk::qOverP,Trk::theta);
232#else
233 double dtheta = std::sqrt((*measPer->covariance())(Trk::theta,Trk::theta));
234 double dqovp = std::sqrt((*measPer->covariance())(Trk::qOverP,Trk::qOverP));
235 double covthetaOvP = (*measPer->covariance())(Trk::qOverP,Trk::theta);
236#endif
237
238
239 double deta = 0.5*dtheta/(std::cos(0.5*theta)*std::cos(0.5*theta)*std::tan(0.5*theta));
240
241#ifdef TRKPARAMETERS_MEASUREDPERIGEE_H
242 double dphi = err.error(Trk::phi0);
243 double dz0 = err.error(Trk::z0);
244 double dd0 = err.error(Trk::d0);
245#else
246 double dphi = std::sqrt((*measPer->covariance())(Trk::phi0,Trk::phi0));
247 double dz0 = std::sqrt((*measPer->covariance())(Trk::z0,Trk::z0));
248 double dd0 = std::sqrt((*measPer->covariance())(Trk::d0,Trk::d0));
249#endif
250
251 double dpT = 0;
252
253
254 double sintheta = std::sin(theta);
255 double costheta = std::cos(theta);
256 double dpt2 = (p*p*sintheta)*(p*p*sintheta)*dqovp*dqovp + (p*costheta)*(p*costheta)*dtheta*dtheta - 2*(p*p*sintheta)*(p*costheta)*covthetaOvP;
257
258 if ( dpt2>0 ) dpT = std::sqrt( dpt2 );
259
260 // Check number of hits
261 // NB: a spacepoint is two offline "hits", so a pixel spacepoint is really
262 // 2 "hits" and an offline SCT "hit" is really a 1D cluster, so two intersecting
263 // stereo clusters making a spacepoint are two "hits"
264 const Trk::TrackSummary *summary = track->trackSummary();
265 int nBlayerHits = 2*summary->get(Trk::numberOfBLayerHits);
266 int nPixelHits = 2*summary->get(Trk::numberOfPixelHits);
267 int nSctHits = summary->get(Trk::numberOfSCTHits);
268 int nStrawHits = summary->get(Trk::numberOfTRTHits);
269 int nTrHits = summary->get(Trk::numberOfTRTHighThresholdHits);
270
271 int nSiHits = nPixelHits + nSctHits;
272 bool expectBL = false; // Not stored for Rec::TrackParticle
273
274 const Trk::FitQuality *quality = track->fitQuality();
275 double chi2 = quality->chiSquared();
276 double dof = quality->numberDoF();
277
278 unsigned bitmap = 0;
279
280
281 unsigned long id = (unsigned long)track;
282
283 for ( int ih=0 ; ih<20 ; ih++ ) {
284 if ( summary->isHit(Trk::DetectorType(ih)) ) bitmap |= ( 1<<hpmap[ih] );
285 }
286
289
290 // std::cout << "fetching author info :" << track->info().trackFitter() << ":"
291 // << track->info().dumpInfo() << ": bm 0x" << std::hex << bitmap << std::dec << std::endl;
292
293 int fitter = track->info().trackFitter();
294 std::string dumpinfo = track->info().dumpInfo();
295
296 int trackAuthor = -1;
297 if ( fitter>0 && fitter<Trk::TrackInfo::NumberOfTrackFitters ) {
298 if ( dumpinfo.find("TRTStandalone")!=std::string::npos) trackAuthor = 2;
299 else if ( dumpinfo.find("TRTSeededTrackFinder")!=std::string::npos) trackAuthor = 1;
300 else trackAuthor = 0;
301 }
302
303#if 0
304 std::cout << "\t\t\tSUTT TP track"
305 << "\teta=" << eta // << " +- " << (*trackitr)->params()->deta()
306 << "\tphi=" << phi // << " +- " << (*trackitr)->params()->dphi()
307 << "\tz0=" << z0
308 << "\tpT=" << pT // << "\t( " << 1/qoverp << ")"
309 << "\td0=" << d0
310 << "\tNsi=" << nSiHits
311 << "\tNtrt=" << nTrHits
312 << "\tNstr=" << nStrawHits
313 << "\tauthor=" << trackAuthor
314 << std::endl;
315#endif
316
317 // Create and save Track
318
319 TIDA::Track* t = new TIDA::Track(eta, phi, z0, d0, pT, chi2, dof,
320 deta, dphi, dz0, dd0, dpT,
321 nBlayerHits, nPixelHits, nSctHits, nSiHits,
322 nStrawHits, nTrHits, bitmap, 0,
323 trackAuthor, false, -1, -1,
324 expectBL, id) ;
325
326 // std::cout << "SUTT TP track " << *t << "\t0x" << std::hex << bitmap << std::dec << std::endl;
327
328 if ( !addTrack( t ) ){
329 delete t;
330 return false;
331 }
332 return true;
333
334 }
335 return false;
336}
337
338
339// extract all the tracks from a TrackParticle collection and add them
340void TrigTrackSelector::selectTracks( const Rec::TrackParticleContainer* trigtracks ) {
341
342 // std::cout << "\t\t\tSUTT \tTrackParticleContainer->size() = " << trigtracks->size() << std::endl;
343
344 Rec::TrackParticleContainer::const_iterator trackitr = trigtracks->begin();
345 Rec::TrackParticleContainer::const_iterator trackend = trigtracks->end();
346
347 while ( trackitr!=trackend ) {
348
349 selectTrack( *trackitr );
350
351 ++trackitr;
352
353 } // loop over tracks
354
355}
356
357
358
359// extract all the tracks from a TruthParticle collection and add them
360void TrigTrackSelector::selectTracks( const TruthParticleContainer* truthtracks ) {
361 // std::cout << "\t\t\tSUTT \tTrackParticleContainer->size() = " << trigtracks->size() << std::endl;
362
363 TruthParticleContainer::const_iterator trackitr = truthtracks->begin();
364 TruthParticleContainer::const_iterator trackend = truthtracks->end();
365
366 while ( trackitr!=trackend ) {
367
368 selectTrack( *trackitr );
369
370 ++trackitr;
371
372 } // loop over tracks
373
374}
375
376
380
381void TrigTrackSelector::truthBeamline( const xAOD::TruthParticleContainer* truthtracks, double& x0, double& y0 ) {
382
383 x0 = 0;
384 y0 = 0;
385
387
388 int Nx = 300;
389 int Ny = 300;
390
392 std::vector<double> xpos(Nx,0);
393 std::vector<double> ypos(Ny,0);
394
396 std::vector<int> xn(Nx,0);
397 std::vector<int> yn(Ny,0);
398
399 int xoffset = Nx/2;;
400 int yoffset = Ny/2;
401
402 double deltax = 3.0/Nx;
403 double deltay = 3.0/Ny;
404
406
407 xAOD::TruthParticleContainer::const_iterator trackitr = truthtracks->begin();
408 xAOD::TruthParticleContainer::const_iterator trackend = truthtracks->end();
409
410 for ( ; trackitr!=trackend ; ++trackitr ) {
411
412 const xAOD::TruthParticle* track = (*trackitr);
413
414 if ( !MC::isStable(track) || !track->hasProdVtx() ) continue;
415
417
418 double xp[3] = { track->prodVtx()->x(), track->prodVtx()->y(), track->prodVtx()->z() };
419
421
422 int ix = xp[0]/deltax + xoffset;
423 int iy = xp[1]/deltay + yoffset;
424
425 if ( ix<0 || ix>=Nx || iy<0 || iy>=Nx ) continue;
426
427 xpos[ix] += xp[0];
428 ypos[iy] += xp[1];
429
430 xn[ix]++;
431 yn[iy]++;
432
433 } // loop over tracks
434
436
437 int imx = 0;
438 int imy = 0;
439
440 for ( size_t i=0 ; i<xpos.size() ; i++ ) {
441 if ( xn[i]>xn[imx] ) imx = i;
442 if ( yn[i]>yn[imy] ) imy = i;
443 }
444
448 if ( xn[imx]>1 ) x0 = xpos[imx]/xn[imx];
449 if ( yn[imy]>1 ) y0 = ypos[imy]/yn[imy];
450
451}
452
453
454// extract all the tracks from a xAOD::TruthParticle collection and add them
455void TrigTrackSelector::selectTracks( const xAOD::TruthParticleContainer* truthtracks) {
456
457 xAOD::TruthParticleContainer::const_iterator trackitr = truthtracks->begin();
458 xAOD::TruthParticleContainer::const_iterator trackend = truthtracks->end();
459
461 double x0 = 0;
462 double y0 = 0;
463
464 truthBeamline( truthtracks, x0, y0 );
465
466
467 for ( ; trackitr!=trackend; ++trackitr) {
468
469
470 // Only select charged final state particles
471 double q = (*trackitr)->charge();
472
474 static const particleType ptype;
475 if ( q==-999 ) q = ptype.charge( (*trackitr)->pdgId() );
476
477 if (q == 0 || !MC::isStable(*trackitr) ) continue;
478
479 // If looking for tau parents, don't select mu or e children
480
481 // select based on the pdg of final state particle
482 bool gotPdgId = true;
483 if (m_selectPdgId!=0) gotPdgId = (*trackitr)->absPdgId()==m_selectPdgId;
484
485 // select based on the pdg of the parent or ancestor
486 bool gotParentPdgId = true;
487 if ( gotPdgId && m_selectParentPdgId!=0 ) gotParentPdgId = fromAncestor(m_selectParentPdgId, (*trackitr))!=nullptr;
490
491 if ( gotParentPdgId && gotPdgId ) selectTrack( *trackitr, x0, y0);
492
493 } // loop over tracks
494
495}
496
497
498
499// add a TruthParticle from a GenParticle - easy, bet it doesn't work
500bool TrigTrackSelector::selectTrack( HepMC::ConstGenParticlePtr track ) {
501
503 if ( !MC::isStable(track) ) return false;
504
506//AV Using memory to get some value is not a good idea. This is not a repruducible/portable way, but I leave it as is.
507 m_id = (unsigned long)(track.get());
508 bool sel;
509 sel = selectTrack( TruthParticle(track) );
510 m_id = 0;
511
512 return sel;
513
514}
515
516
517// add a TruthParticle
518bool TrigTrackSelector::selectTrack( const TruthParticle& track ) {
519
520 return selectTrack( &track );
521
522}
523
524
525// add a TruthParticle
526bool TrigTrackSelector::selectTrack( const TruthParticle* track ) {
527 TIDA::Track* t = makeTrack( track, m_id );
528 if ( t == 0 ) return false;
529 if ( !addTrack(t) ) {
530 delete t;
531 return false;
532 }
533 return true;
534}
535
536
537
538
539// add an xAOD::TruthParticle
540bool TrigTrackSelector::selectTrack( const xAOD::TruthParticle* track, double x0, double y0) {
541 if ( track ) {
542
543 if (!MC::isStable(track) ) return false;
544
546 const xAOD::TruthParticle* measPer = track;
547
548 double pT = measPer->pt();
549 double eta = measPer->eta();
550 double phi = measPer->phi();
551
552 // AAARCH!!!!! the TrackParticle pT is NOT SIGNED!!!! ( I ask you! )
553 if ( measPer->charge()<0 && pT>0 ) pT *= -1;
554 double q = track->charge();
555
556 static const particleType ptype;
557
559 if ( q==-999 ) q = ptype.charge( track->pdgId() );
560
562 if ( q==0 ) return 0;
563
567 // double xbeam = getBeamX(); // track->vx();
568 // double ybeam = getBeamY(); // track->vy();
569 // double zbeam = getBeamZ(); // track->vz();
570
571 if ( !track->hasProdVtx() ) return false;
572
574
575 double xp[3] = { measPer->prodVtx()->x(), measPer->prodVtx()->y(), measPer->prodVtx()->z() };
576 double xb[3] = { xp[0]-x0, xp[1]-y0, measPer->prodVtx()->z() };
577 double xd[3] = { 0, 0, 0 };
578
579 if ( track->hasDecayVtx() ) {
580 xd[0] = track->decayVtx()->x();
581 xd[1] = track->decayVtx()->y();
582 xd[2] = track->decayVtx()->z();
583 }
584
585 double rp = std::sqrt( xp[0]*xp[0] + xp[1]*xp[1] );
586 double rd = std::sqrt( xd[0]*xd[0] + xd[1]*xd[1] );
587
590 double theta = 2*std::atan( std::exp( -eta ) );
591 double z0 = xb[2] - (xb[0]*std::cos(phi) + xb[1]*std::sin(phi))/std::tan(theta);
592 double d0 = xb[1]*std::cos(phi) - xb[0]*std::sin(phi);
593
594 bool final_state = false;
595
607 const double inner_radius = m_radius;
608 const double outer_radius = m_radius;
609
610 if ( ( track->hasProdVtx() && rp<=inner_radius ) &&
611 ( !track->hasDecayVtx() || rd>outer_radius ) ) final_state = true;
612
613 if ( !final_state ) return false;
614
615 double deta = 0;
616 double dphi = 0;
617 double dz0 = 0;
618 double dd0 = 0;
619
620 double dpT = 0;
621
622 int nBlayerHits = 0;
623 int nPixelHits = 0;
624 int nSctHits = 0;
625 int nStrawHits = 0;
626 int nTrtHits = 0;
627
628 double chi2 = 0;
629 double dof = 0;
630
631 bool expectBL = false;
632
633 nSctHits += 0;
634 nPixelHits += 0;
635
637 int nSiHits = 0;
638
639 unsigned long id = (unsigned long)track;
640
641 unsigned bitmap = 0;
642
643 int trackAuthor = track->pdgId();
644 int uniqueID = HepMC::uniqueID(track);
645
646#if 0
647 std::cout << "\t\t\tSUTT TP track"
648 << "\teta=" << eta
649 << "\tphi=" << phi
650 << "\tz0=" << z0
651 << "\tpT=" << pT
652 << "\td0=" << d0
653 << "\tauthor=" << trackAuthor
654 << "\tVTX x " << xp[0]<< "\ty " << xp[1] << "\tz " << xp[2]
655 << std::endl;
656#endif
657
658 // Create and save Track
659
660 TIDA::Track* t = new TIDA::Track( eta, phi, z0, d0, pT, chi2, dof,
661 deta, dphi, dz0, dd0, dpT,
662 nBlayerHits, nPixelHits, nSctHits, nSiHits,
663 nStrawHits, nTrtHits, bitmap, 0,
664 trackAuthor, false, uniqueID, -1,
665 expectBL, id) ;
666
668 // std::cout << "SUTT TP track " << *t << "\t0x" << std::hex << bitmap << std::dec << std::endl;
669
670 // addTrack applies additional cuts using the Filter
671 if ( !addTrack( t ) ){
672 delete t;
673 return false;
674 }
675 }
676 return false;
677
678}
679
680
681
682
683// make a TIDA::Track from a GenParticle
684TIDA::Track* TrigTrackSelector::makeTrack(HepMC::ConstGenParticlePtr track ) {
685//AV Using memory to get some value is not a good idea. This is not a repruducible/portable way, but I leave it as is.
686 unsigned long id = (unsigned long)(track.get());
687 TruthParticle t = TruthParticle(track);
688 return makeTrack( &t, id );
689}
690
691// make a TIDA::Track from a TruthParticle
692TIDA::Track* TrigTrackSelector::makeTrack( const TruthParticle* track, unsigned long tid ) {
693
694 if ( track==0 ) return 0;
695 if ( !MC::isStable(track)) return 0;
696
697
698 double phi = track->phi();
699 double eta = track->eta();
700
704 double xp[3] = { 0, 0, 0 };
705
706 if ( track->genParticle()->production_vertex() ) {
707 xp[0] = track->genParticle()->production_vertex()->position().x();
708 xp[1] = track->genParticle()->production_vertex()->position().y();
709 xp[2] = track->genParticle()->production_vertex()->position().z();
710 }
711
712 // CHANGED BY JK - z0 with respect to (0,0)
713 // double z0 = xp[2];
714 double theta = 2*std::atan( exp( (-1)*eta ) );
715 double z0 = xp[2] - (xp[0]*std::cos(phi) + xp[1]*std::sin(phi))/std::tan(theta);
716
717 double xd[3] = { 0, 0, 0 };
718
719 if ( track->genParticle()->end_vertex() ) {
720 xd[0] = track->genParticle()->end_vertex()->position().x();
721 xd[1] = track->genParticle()->end_vertex()->position().y();
722 xd[2] = track->genParticle()->end_vertex()->position().z();
723 }
724
725 double rp = std::sqrt( xp[0]*xp[0] + xp[1]*xp[1] );
726 double rd = std::sqrt( xd[0]*xd[0] + xd[1]*xd[1] );
727
728
729 bool final_state = false;
730
742 const double inner_radius = m_radius;
743 const double outer_radius = m_radius;
744 if ( ( track->genParticle()->production_vertex() && rp<=inner_radius ) &&
745 ( track->genParticle()->end_vertex()==0 || rd>outer_radius ) ) final_state = true;
746
747
748 if ( !final_state ) return 0;
749
760
761 double q = track->charge();
762
763 static const particleType ptype;
764
766 if ( q==-999 ) q = ptype.charge( track->pdgId() );
767
769 if ( q==0 ) return 0;
770
771 double pT = q*track->pt();
772
773
774 double d0 = 0;
775
780
782
783
784 // CHANGED BY JK - d0 with respect to (0,0)
785 // d0 = q*rp*std::sin(phi);
786 d0 = xp[1]*std::cos(phi) - xp[0]*std::sin(phi);
787
788
790
791 double dz0 = 0;
792 double dd0 = 0;
793
794 correctToBeamline( z0, dz0, d0, dd0, theta, phi );
795
796
797
799
801 int author = track->pdgId();
802 int uniqueID = HepMC::uniqueID(track);
803
804
805 unsigned long id = (unsigned long)track;
806 if ( tid!=0 ) id = tid;
807
811
812
814
815 // std::cout << "\t\t\tSUTT Truth track"
816 // << "\teta=" << eta // << " +- " << (*trackitr)->params()->deta()
817 // << "\tphi=" << phi // << " +- " << (*trackitr)->params()->dphi()
818 // << "\tz0=" << z0
819 // << "\tpT=" << pT // << "\t( " << 1/qoverp << ")"
820 // << "\td0=" << d0
821 // << "\tauthor=" << author
822 // << std::endl;
823
824
825
826 TIDA::Track* t = new TIDA::Track(eta, phi, z0, d0, pT, 0, 0,
827 0, 0, 0, 0, 0,
828 0, 0, 0, 0,
829 0, 0, 0, 0,
830 author, false, uniqueID, -1,
831 false,
832 id ) ;
833
834 return t;
835
836}
837
838
839
840// add a Trk::Track
841bool TrigTrackSelector::selectTrack( const Trk::Track* track ) {
842
843 // do the track extraction stuff here....
844
845 static const int hpmap[20] = { 0, 1, 2, 7, 8, 9, 3, 4, 5, 6, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19 };
846 // std::cout << __FILE__<< " " <<__LINE__<<std::endl;
847
848 if ( track ) {
849
850 // const Trk::Perigee* startPerigee = track->perigeeParameters();
851
852#ifdef TRKPARAMETERS_MEASUREDPERIGEE_H
853 const Trk::MeasuredPerigee* startPerigee = dynamic_cast<const Trk::MeasuredPerigee*>(track->perigeeParameters());
854 // const Trk::MeasuredPerigee* measPer = startPerigee; // just out of laziness
855#else
856 const Trk::Perigee* startPerigee = track->perigeeParameters();
857 const Trk::Perigee* measPer = startPerigee; // just out of laziness
858#endif
859
860
861 // CLHEP::HepVector perigeeParams = measPer->parameters();
862 // double pT = measPer->pT();
863 // double eta = measPer->eta();
864 // double phi = perigeeParams[Trk::phi0];
865 // double z0 = perigeeParams[Trk::z0];
866 // double d0 = perigeeParams[Trk::d0];
867 // // AAARCH!!!!! the TrackParticle pT is NOT SIGNED!!!! ( I ask you! )
868 // if ( perigeeParams[Trk::qOverP]<0 ) pT *= -1;
869 // std::cout <<pT1<<" pt1vspT "<<pT<<std::endl;
870
871 if (startPerigee){
872
873 double theta = startPerigee->parameters()[Trk::theta];
874 double p = 1/startPerigee->parameters()[Trk::qOverP];
875 double qOverPt = startPerigee->parameters()[Trk::qOverP]/std::sin(theta);
876 double charge = startPerigee->charge();
877 double eta = startPerigee->eta();
878 double phi = startPerigee->parameters()[Trk::phi0];
879 double z0 = startPerigee->parameters()[Trk::z0];
880 double d0 = startPerigee->parameters()[Trk::d0];
881 // double pT = (1./qOverPt)*(charge);
882 double pT = (1./qOverPt); // always use signed PT
883
884 if ( charge<0 && pT>0 ) pT *= -1;
885 if ( charge<0 && p>0 ) p *= -1;
886
887
888
889#ifdef TRKPARAMETERS_MEASUREDPERIGEE_H
890 const Trk::ErrorMatrix err = startPerigee->localErrorMatrix();
891 double dtheta = err.error(Trk::theta);
892 double dqovp = err.error(Trk::qOverP);
893 double covthetaOvP = err.covValue(Trk::qOverP,Trk::theta);
894
895 double dphi = err.error(Trk::phi0);
896 double dz0 = err.error(Trk::z0);
897 double dd0 = err.error(Trk::d0);
898#else
899 double dtheta = std::sqrt((*measPer->covariance())(Trk::theta,Trk::theta));
900 double dqovp = std::sqrt((*measPer->covariance())(Trk::qOverP,Trk::qOverP));
901 double covthetaOvP = (*measPer->covariance())(Trk::qOverP,Trk::theta);
902
903 double dphi = std::sqrt((*measPer->covariance())(Trk::phi0,Trk::phi0));
904 double dz0 = std::sqrt((*measPer->covariance())(Trk::z0,Trk::z0));
905 double dd0 = std::sqrt((*measPer->covariance())(Trk::d0,Trk::d0));
906#endif
907
908 double deta = 0.5*dtheta/(std::cos(0.5*theta)*std::cos(0.5*theta)*std::tan(0.5*theta));
909
910
911 if ( m_correctTrkTracks ) correctToBeamline( z0, dz0, d0, dd0, theta, phi );
912
913 double dpT = 0;
914
915
916 double sintheta = std::sin(theta);
917 double costheta = std::cos(theta);
918 double dpT2 = (p*p*sintheta)*(p*p*sintheta)*dqovp*dqovp + (p*costheta)*(p*costheta)*dtheta*dtheta - 2*(p*p*sintheta)*(p*costheta)*covthetaOvP;
919
920 if ( dpT2>0 ) dpT = std::sqrt( dpT2 );
921
922 // Check number of hits
923 // NB: a spacepoint is two offline "hits", so a pixel spacepoint is really
924 // 2 "hits" and an offline SCT "hit" is really a 1D cluster, so two intersetcting
925 // stereo clusters making a spacepoint are two "hits"
926 // const Trk::TrackSummary *summary = dynamic_cast<const Trk::TrackSummary*>(track->trackSummary());
927 //ToolHandle< Trk::ITrackSummaryTool > m_trackSumTool;
928 //m_trackSumTool = ToolHandle<Trk::ITrackSummaryTool>("Trk::TrackSummaryTool/InDetTrackSummaryTool");
929 //const Trk::TrackSummary* summary = NULL;
930 //summary = m_trackSumTool->createSummary(*track);
931
932 const Trk::TrackSummary * summary = track->trackSummary();
933 int nBlayerHits = 0;
934 int nPixelHits = 0;
935 int nSctHits = 0;
936 int nStrawHits = 0;
937 int nTrHits = 0;
938 int nSiHits = 0;
939 bool expectBL = false; // Not stored for Trk::Track
940 unsigned bitmap = 0;
941
942 if(summary==0){
943 std::cerr << "Could not create TrackSummary - Track will likely fail hits requirements" << std::endl;
944 }
945 else{
946 nBlayerHits = 2*summary->get(Trk::numberOfBLayerHits);
947 nPixelHits = 2*summary->get(Trk::numberOfPixelHits);
948 nSctHits = summary->get(Trk::numberOfSCTHits);
949 nStrawHits = summary->get(Trk::numberOfTRTHits);
950 nTrHits = summary->get(Trk::numberOfTRTHighThresholdHits);
951 nSiHits = nPixelHits + nSctHits;
952
953 for ( int ih=0 ; ih<20 ; ih++ ) {
954 if ( summary->isHit(Trk::DetectorType(ih)) ) bitmap |= ( 1<<hpmap[ih] );
955 }
956 }
957
958 unsigned long id = (unsigned long)track;
959 double chi2 = 0;
960 double dof = 0;
961 //const Trk::FitQuality *quality = dynamic_cast<const Trk::FitQuality*>(track->fitQuality());
962 const Trk::FitQuality *quality = (track->fitQuality());
963 if(quality==0) std::cerr << "Could not create FitQuality - Track will likely fail hits requirements" << std::endl;
964 else{
965 chi2 = quality->chiSquared();
966 dof = quality->numberDoF();
967 }
968
969 int trackAuthor = -1;
970
973
974 // std::cout << "fetching author info :" << track->info().trackFitter() << ":"
975 // << track->info().dumpInfo() << ": bm 0x" << std::hex << bitmap << std::dec << std::endl;
976
977 int fitter = track->info().trackFitter();
978 // std::string dumpinfo = track->info().dumpInfo();
979
980 if ( fitter>0 && fitter<Trk::TrackInfo::NumberOfTrackFitters ) {
981 if ((track->info().dumpInfo()).find("TRTStandalone") != std::string::npos) trackAuthor = 2;
982 else if ((track->info().dumpInfo()).find("TRTSeededTrackFinder") != std::string::npos) trackAuthor = 1;
983 else trackAuthor = 0;
984 }
985
986 #if 0
987 std::cout << "\t\t\tSUTT TP track"
988 << "\teta=" << eta // << " +- " << (*trackitr)->params()->deta()
989 << "\tphi=" << phi // << " +- " << (*trackitr)->params()->dphi()
990 << "\tz0=" << z0
991 << "\tpT=" << pT // << "\t( " << 1/qoverp << ")"
992 << "\td0=" << d0
993 << "\tNsi=" << nSiHits
994 << "\tNtrt=" << nTrHits
995 << "\tNstr=" << nStrawHits
996 << "\tauthor=" << trackAuthor
997 << std::endl;
998 #endif
999 // Create and save Track
1000 TIDA::Track* t = new TIDA::Track(eta, phi, z0, d0, pT, chi2, dof,
1001 deta, dphi, dz0, dd0, dpT,
1002 nBlayerHits, nPixelHits, nSctHits, nSiHits,
1003 nStrawHits, nTrHits, bitmap, 0,
1004 trackAuthor, false, -1, -1,
1005 expectBL, id) ;
1006
1007 if ( !addTrack( t ) ){
1008 delete t;
1009 return false;
1010 }
1011 return true;
1012
1013 //std::cout << "SUTT TP track " << *t << "\t0x" << std::hex << bitmap << std::dec << std::endl;
1014 }
1015 }
1016
1017 return false;
1018}
1019
1020// extract all the tracks from a TrackCollection and add them
1021void TrigTrackSelector::selectTracks( const TrackCollection* trigtracks ) {
1022
1023 // std::cout << "\t\t\tSUTT \tTrackContainer->size() = " << trigtracks->size() << std::endl;
1024
1025 TrackCollection::const_iterator trackitr = trigtracks->begin();
1026 TrackCollection::const_iterator trackend = trigtracks->end();
1027
1028 while ( trackitr!=trackend ) {
1029 selectTrack( *trackitr );
1030 ++trackitr;
1031 } // loop over tracks
1032
1033}
1034
1035
1036
1037bool TrigTrackSelector::selectTrack( const xAOD::TrackParticle* track, void* ) {
1038 // do the track extraction stuff here....
1039
1040 // static int hpmap[20] = { 0, 1, 2, 7, 8, 9, 3, 4, 5, 6, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19 };
1041
1042 if ( track ) {
1043
1045 const xAOD::TrackParticle* measPer = track;
1046
1047 // CLHEP::HepVector perigeeParams = measPer->parameters();
1048
1049 double pT = measPer->pt();
1050 double eta = measPer->eta();
1051 double phi = measPer->phi0();
1052 double z0 = measPer->z0() + measPer->vz();
1053 double d0 = measPer->d0();
1054
1055 double theta = measPer->theta();
1056
1061 // if ( measPer->qOverP()==0 ) return false;
1062 if ( measPer->qOverP()==0 ) throw std::runtime_error( "probable corrupted track - this should never happen" );
1063 double p = 1/measPer->qOverP();
1064
1065 // AAARCH!!!!! the TrackParticle pT is NOT SIGNED!!!! ( I ask you! )
1066 if ( measPer->qOverP()<0 && pT>0 ) pT *= -1;
1067
1068 double dtheta = std::sqrt(measPer->definingParametersCovMatrix()(Trk::theta,Trk::theta));
1069 double dqovp = std::sqrt(measPer->definingParametersCovMatrix()(Trk::qOverP,Trk::qOverP));
1070 double covthetaOvP = measPer->definingParametersCovMatrix()(Trk::qOverP,Trk::theta);
1071
1072 double deta = 0.5*dtheta/(std::cos(0.5*theta)*std::cos(0.5*theta)*std::tan(0.5*theta)); // ???? CHECK THIS <<--
1073
1074 double dphi = std::sqrt(measPer->definingParametersCovMatrix()(Trk::phi0,Trk::phi0));
1075 double dz0 = std::sqrt(measPer->definingParametersCovMatrix()(Trk::z0,Trk::z0));
1076 double dd0 = std::sqrt(measPer->definingParametersCovMatrix()(Trk::d0,Trk::d0));
1077
1078 double dpT = 0;
1079
1080
1082 // if ( m_xBeam!=0 || m_yBeam!=0 ) correctToBeamline( z0, dz0, d0, dd0, theta, phi );
1083
1084
1085 double sintheta = std::sin(theta);
1086 double costheta = std::cos(theta);
1087 double dpt2 = (p*p*sintheta)*(p*p*sintheta)*dqovp*dqovp + (p*costheta)*(p*costheta)*dtheta*dtheta - 2*(p*p*sintheta)*(p*costheta)*covthetaOvP;
1088
1089 if ( dpt2>0 ) dpT = std::sqrt( dpt2 );
1090
1091 // Check number of hits
1092 // NB: a spacepoint is two offline "hits", so a pixel spacepoint is really
1093 // 2 "hits" and an offline SCT "hit" is really a 1D cluster, so two intersetcting
1094 // stereo clusters making a spacepoint are two "hits"
1095
1096 uint8_t sum_nBlayerHits = 0;
1097 track->summaryValue( sum_nBlayerHits, xAOD::numberOfInnermostPixelLayerHits);
1098 int nBlayerHits = 2*sum_nBlayerHits;
1099
1100 uint8_t sum_nPixelHits = 0;
1101 track->summaryValue( sum_nPixelHits, xAOD::numberOfPixelHits);
1102 int nPixelHits = 2*sum_nPixelHits;
1103
1104 uint8_t sum_nSctHits = 0;
1105 track->summaryValue( sum_nSctHits, xAOD::numberOfSCTHits);
1106 int nSctHits = sum_nSctHits;
1107
1108 uint8_t sum_nStrawHits = 0;
1109 track->summaryValue( sum_nStrawHits, xAOD::numberOfTRTHits);
1110 int nStrawHits = sum_nStrawHits;
1111
1112 uint8_t sum_nTrtHits = 0;
1113 track->summaryValue( sum_nTrtHits, xAOD::numberOfTRTHighThresholdHits);
1114 int nTrtHits = sum_nTrtHits;
1115
1116
1117 uint8_t sum_expectBL = 0;
1118 track->summaryValue( sum_expectBL, xAOD::expectInnermostPixelLayerHit);
1119 bool expectBL = ( sum_expectBL ? true : false );
1120
1122
1123 uint8_t sum_sctholes = 0;
1124 track->summaryValue( sum_sctholes, xAOD::numberOfSCTHoles);
1125
1126 uint8_t sum_pixholes = 0;
1127 track->summaryValue( sum_pixholes, xAOD::numberOfPixelHoles);
1128
1132
1133 nSctHits += 1000*sum_sctholes;
1134 nPixelHits += 1000*sum_pixholes;
1135
1137 int nSiHits = nPixelHits + nSctHits;
1138
1140
1141 double chi2 = track->chiSquared();
1142 double dof = track->numberDoF();
1143
1144 unsigned long id = (unsigned long)track;
1145
1146 unsigned bitmap = track->hitPattern();
1147
1148
1149
1150 double xbeam = track->vx();
1151 double ybeam = track->vy();
1152 double zbeam = track->vz();
1153
1154 if ( xbeam!=getBeamX() || ybeam!=getBeamY() || zbeam!=getBeamZ() ) setBeamline( xbeam, ybeam, zbeam );
1155
1156 int trackAuthor = 0;
1157
1158 int fitter = track->trackFitter();
1159 std::bitset<xAOD::NumberOfTrackRecoInfo> patternrec = track->patternRecoInfo();
1160
1161 //int icount = 0;<- never used if section below is commented
1162 for ( unsigned ipr=patternrec.size() ; ipr-- ; ) {
1163 if ( patternrec[ipr] ) {
1164 //icount++; <- never used if section below is commented
1165 trackAuthor |= (ipr >> 16);
1166 // static bool first = true;
1167 // if ( first && icount>1 ) {
1168 // std::cerr << "more than one pattern rec strategy " << ipr << "\t(suppressing further output)" << std::endl;
1169 // first = false;
1170 // }
1171 }
1172 }
1173
1174 trackAuthor |= fitter;
1175
1176 // if ( fitter>0 && fitter<Trk::TrackInfo::NumberOfTrackFitters ) {
1177 // if ( dumpinfo.find("TRTStandalone")!=std::string::npos) trackAuthor = 2;
1178 // else if ( dumpinfo.find("TRTSeededTrackFinder")!=std::string::npos) trackAuthor = 1;
1179 // else trackAuthor = 0;
1180 // }
1181
1182#if 0
1183 std::cout << "\t\t\tSUTT TP track"
1184 << "\teta=" << eta // << " +- " << (*trackitr)->params()->deta()
1185 << "\tphi=" << phi // << " +- " << (*trackitr)->params()->dphi()
1186 << "\tz0=" << z0
1187 << "\tpT=" << pT // << "\t( " << 1/qoverp << ")"
1188 << "\td0=" << d0
1189 << "\tNsi=" << nSiHits
1190 // << "\tNtrt=" << nTrtHits
1191 // << "\tNstr=" << nStrawHits
1192 << "\tfitter=" << fitter
1193 << "\tauthor=" << trackAuthor
1194 << "\tVTX x " << track->vx() << "\ty " << track->vy() << "\tz " << track->vz()
1195 << std::endl;
1196#endif
1197
1198 // Create and save Track
1199
1200 TIDA::Track* t = new TIDA::Track( eta, phi, z0, d0, pT, chi2, dof,
1201 deta, dphi, dz0, dd0, dpT,
1202 nBlayerHits, nPixelHits, nSctHits, nSiHits,
1203 nStrawHits, nTrtHits, bitmap, 0,
1204 trackAuthor, false, -1, -1,
1205 expectBL, id) ;
1206
1207 // std::cout << "SUTT TP track " << *t << "\t0x" << std::hex << bitmap << std::dec << std::endl;
1208
1209 if ( !addTrack( t ) ){
1210 delete t;
1211 return false;
1212 }
1213 return true;
1214
1215 }
1216 return false;
1217
1218
1219}
1220
1221
1222void TrigTrackSelector::selectTracks( const xAOD::TrackParticleContainer* tracks, void* ) {
1223 // std::cout << "\t\t\tSUTT \tTrackContainer->size() = " << trigtracks->size() << std::endl;
1224 xAOD::TrackParticleContainer::const_iterator trackitr = tracks->begin();
1225 xAOD::TrackParticleContainer::const_iterator trackend = tracks->end();
1226 while ( trackitr!=trackend ) {
1227 selectTrack( *trackitr );
1228 ++trackitr;
1229 } // loop over tracks
1230}
1231
1232
1233
1234void TrigTrackSelector::selectTracks( xAOD::TrackParticleContainer::const_iterator trackitr,
1235 xAOD::TrackParticleContainer::const_iterator trackend, void* ) {
1238 // std::cout << "\t\t\tSUTT \tTrackContainer->size() = " << trigtracks->size() << std::endl;
1239 while ( trackitr!=trackend ) {
1240 selectTrack( *trackitr );
1241 ++trackitr;
1242 } // loop over tracks
1243}
1244
1245
1246void TrigTrackSelector::selectTracks( const TrackParticleLinks_t& tracks ) {
1247 for( const auto& track : tracks ) selectTrack( *track );
1248}
1249
1250
1251
1252
1253
1254void TrigTrackSelector::correctToBeamline( double& z0, double& dz0,
1255 double& d0, double& dd0,
1256 double theta, double phi ) {
1257
1259
1260 // if ( m_first ) {
1261 // if ( m_xBeam==0 && m_yBeam==0 ) {
1262 // std::cerr << "TrigTrackSelector::correctToBeamline() WARNING -- Beamline set to (0,0) -- WARNING" << std::endl;
1263 // }
1264 // else {
1265 // std::cout << "TrigTrackSelector::correctToBeamline() Beamline set to " << m_xBeam << " " << m_yBeam << std::endl;
1266 // }
1267 // m_first = false;
1268 // }
1269
1270
1271 // double theta = 2*std::atan( exp( (-1)*eta ) );
1272 double z0t = z0 + ((std::cos(phi)*m_xBeam + std::sin(phi)*m_yBeam)/std::tan(theta));
1273 double a0t = d0 + std::sin(phi)*m_xBeam - std::cos(phi)*m_yBeam;
1274
1276 double dz0t = dz0 + ((std::cos(phi)*m_xBeam + std::sin(phi)*m_yBeam)/std::tan(theta));
1277 double da0t = dd0 + std::sin(phi)*m_xBeam - std::cos(phi)*m_yBeam;
1278
1279 z0 = z0t;
1280 d0 = a0t;
1281
1282 dz0 = dz0t;
1283 dd0 = da0t;
1284}
1285
eta
Scalar eta() const
pseudorapidity method
Definition AmgMatrixBasePlugin.h:83
phi
Scalar phi() const
phi method
Definition AmgMatrixBasePlugin.h:67
theta
Scalar theta() const
theta method
Definition AmgMatrixBasePlugin.h:75
charge
double charge(const T &p)
Definition AtlasPID.h:997
HepMCHelpers.h
ATLAS-specific HepMC functions.
rp
ReadCards * rp
Definition IReadCards.cxx:26
TrackCollection
DataVector< Trk::Track > TrackCollection
This typedef represents a collection of Trk::Track objects.
Definition TrackCollection.h:19
DataVector< TrigInDetTrack >::const_iterator
DataModel_detail::const_iterator< DataVector > const_iterator
Definition DataVector.h:838
DataVector::end
const_iterator end() const noexcept
Return a const_iterator pointing past the end of the collection.
DataVector::begin
const_iterator begin() const noexcept
Return a const_iterator pointing at the beginning of the collection.
TrackSelector::addTrack
virtual bool addTrack(TIDA::Track *t, bool(*f)(const TIDA::Track *)=0)
Definition Trigger/TrigAnalysis/TrigInDetAnalysis/TrigInDetAnalysis/TrackSelector.h:35
TrigInDetTrackTruthMap::truth
const TrigInDetTrackTruth * truth(const TrigInDetTrack *p_trig_trk) const
Definition TrigInDetTrackTruthMap.cxx:82
TrigInDetTrackTruth::bestSiMatch
const HepMcParticleLink * bestSiMatch() const
returns best match according to the number of hits
Definition TrigInDetTrackTruth.cxx:186
TrigInDetTrackTruth::nrMatches
unsigned int nrMatches() const
returns number of matching particles
Definition TrigInDetTrackTruth.cxx:240
TrigInDetTrack
represents a LVL2 ID track
Definition TrigInDetTrack.h:34
TrigTrackSelector::setBeamline
void setBeamline(double x, double y, double z=0)
Definition TrigTrackSelector.h:94
TrigTrackSelector::fromAncestor
const xAOD::TruthParticle * fromAncestor(const int pdg_id, const xAOD::TruthParticle *p) const
recursive functions to identify whether a particle comes from some particle of a specific PDG ID,...
Definition TrigTrackSelector.cxx:39
TrigTrackSelector::selectTrack
bool selectTrack(const TrigInDetTrack *track, const TrigInDetTrackTruthMap *truthMap=0)
neater code to make use of vector function also for a single ancestor pdgid, instead of the full code...
Definition TrigTrackSelector.cxx:116
TrigTrackSelector::getBeamY
double getBeamY() const
Definition TrigTrackSelector.h:97
TrigTrackSelector::truthBeamline
void truthBeamline(const xAOD::TruthParticleContainer *truthtracks, double &x0, double &y0)
extract all the tracks from a xAOD::TruthParticle collection and histogram the x and y production coo...
Definition TrigTrackSelector.cxx:381
TrigTrackSelector::getBeamX
double getBeamX() const
Definition TrigTrackSelector.h:96
TrigTrackSelector::TrigTrackSelector
TrigTrackSelector(TrackFilter *selector)
use a radius of 47 mm corresponding to the Run 1 pixel inner radius For the IBL it should be 32 mm,...
Definition TrigTrackSelector.h:70
TrigTrackSelector::selectTracks
void selectTracks(const TrigInDetTrackCollection *trigtracks, const TrigInDetTrackTruthMap *truthMap=0)
Definition TrigTrackSelector.cxx:188
TrigTrackSelector::makeTrack
TIDA::Track * makeTrack(HepMC::ConstGenParticlePtr track)
Definition TrigTrackSelector.cxx:684
TrigTrackSelector::getBeamZ
double getBeamZ() const
Definition TrigTrackSelector.h:98
TrigTrackSelector::correctToBeamline
void correctToBeamline(double &z0, double &dz0, double &d0, double &dd0, double theta, double phi)
Definition TrigTrackSelector.cxx:1254
TrigTrackSelector::s_default_radius
static const double s_default_radius
NB: This was 47 for Run 2, but with the addition of the IBL it should be 32 It was kept at 47 for all...
Definition TrigTrackSelector.h:205
TrigTrackSelector::TrackParticleLinks_t
std::vector< ElementLink< xAOD::TrackParticleContainer > > TrackParticleLinks_t
Definition TrigTrackSelector.h:62
Trk::FitQuality
Class to represent and store fit qualities from track reconstruction in terms of and number of degre...
Definition FitQuality.h:97
Trk::FitQuality::numberDoF
int numberDoF() const
returns the number of degrees of freedom of the overall track or vertex fit as integer
Definition FitQuality.h:60
Trk::FitQuality::chiSquared
double chiSquared() const
returns the of the overall track fit
Definition FitQuality.h:56
Trk::ParametersBase::eta
double eta() const
Access method for pseudorapidity - from momentum.
Trk::ParametersBase::charge
double charge() const
Returns the charge.
Trk::ParametersBase::pT
double pT() const
Access method for transverse momentum.
Trk::TrackSummary
A summary of the information contained by a track.
Definition Tracking/TrkEvent/TrkTrackSummary/TrkTrackSummary/TrackSummary.h:287
particleType::charge
double charge(int id) const
Definition particleType.h:121
xAOD::TrackParticle_v1::z0
float z0() const
Returns the parameter.
xAOD::TrackParticle_v1::theta
float theta() const
Returns the parameter, which has range 0 to .
xAOD::TrackParticle_v1::definingParametersCovMatrix
const ParametersCovMatrix_t definingParametersCovMatrix() const
Returns the 5x5 symmetric matrix containing the defining parameters covariance matrix.
Definition TrackParticle_v1.cxx:260
xAOD::TrackParticle_v1::vz
float vz() const
The z origin for the parameters.
xAOD::TrackParticle_v1::d0
float d0() const
Returns the parameter.
xAOD::TrackParticle_v1::qOverP
float qOverP() const
Returns the parameter.
xAOD::TrackParticle_v1::pt
virtual double pt() const override final
The transverse momentum ( ) of the particle.
Definition TrackParticle_v1.cxx:75
xAOD::TrackParticle_v1::eta
virtual double eta() const override final
The pseudorapidity ( ) of the particle.
Definition TrackParticle_v1.cxx:79
xAOD::TrackParticle_v1::phi0
float phi0() const
Returns the parameter, which has range to .
Definition TrackParticle_v1.cxx:150
xAOD::TruthParticle_v1::pt
virtual double pt() const override final
The transverse momentum ( ) of the particle.
Definition TruthParticle_v1.cxx:161
xAOD::TruthParticle_v1::prodVtx
const TruthVertex_v1 * prodVtx() const
The production vertex of this particle.
Definition TruthParticle_v1.cxx:75
xAOD::TruthParticle_v1::charge
double charge() const
Physical charge.
xAOD::TruthParticle_v1::eta
virtual double eta() const override final
The pseudorapidity ( ) of the particle.
Definition TruthParticle_v1.cxx:169
xAOD::TruthParticle_v1::phi
virtual double phi() const override final
The azimuthal angle ( ) of the particle.
Definition TruthParticle_v1.cxx:176
xAOD::TruthVertex_v1::z
float z() const
Vertex longitudinal distance along the beam line form the origin.
xAOD::TruthVertex_v1::y
float y() const
Vertex y displacement.
xAOD::TruthVertex_v1::x
float x() const
Vertex x displacement.
zbeam
Definition zbeam.h:24
chi2
double chi2(TH1 *h0, TH1 *h1)
Definition comparitor.cxx:525
HepMC::INVALID_PARTICLE_ID
constexpr int INVALID_PARTICLE_ID
Definition MagicNumbers.h:58
HepMC::uniqueID
int uniqueID(const T &p)
Definition MagicNumbers.h:89
HepMC::ConstGenParticlePtr
HepMC3::ConstGenParticlePtr ConstGenParticlePtr
Definition GenParticle.h:20
MC::isElectron
bool isElectron(const T &p)
Definition HepMCHelpers.h:203
MC::isStable
bool isStable(const T &p)
Identify if the particle is stable, i.e. has not decayed.
Definition HepMCHelpers.h:51
MC::isMuon
bool isMuon(const T &p)
Definition HepMCHelpers.h:206
Trk::Perigee
ParametersT< TrackParametersDim, Charged, PerigeeSurface > Perigee
Definition Tracking/TrkEvent/TrkParameters/TrkParameters/TrackParameters.h:33
Trk::MeasuredPerigee
Perigee MeasuredPerigee
Definition MeasuredPerigeeCnv_p1.h:23
Trk::phi0
@ phi0
Definition ParamDefs.h:65
Trk::theta
@ theta
Definition ParamDefs.h:66
Trk::qOverP
@ qOverP
perigee
Definition ParamDefs.h:67
Trk::d0
@ d0
Definition ParamDefs.h:63
Trk::z0
@ z0
Definition ParamDefs.h:64
Trk::DetectorType
DetectorType
enumerates the various detector types currently accessible from the isHit() method.
Definition Tracking/TrkEvent/TrkTrackSummary/TrkTrackSummary/TrackSummary.h:231
Trk::numberOfPixelHits
@ numberOfPixelHits
number of pixel layers on track with absence of hits
Definition Tracking/TrkEvent/TrkTrackSummary/TrkTrackSummary/TrackSummary.h:57
Trk::numberOfBLayerHits
@ numberOfBLayerHits
these are the hits in the 0th pixel layer?
Definition Tracking/TrkEvent/TrkTrackSummary/TrkTrackSummary/TrackSummary.h:51
Trk::numberOfTRTHighThresholdHits
@ numberOfTRTHighThresholdHits
total number of TRT hits which pass the high threshold
Definition Tracking/TrkEvent/TrkTrackSummary/TrkTrackSummary/TrackSummary.h:87
xAOD::TrackParticle
TrackParticle_v1 TrackParticle
Reference the current persistent version:
Definition Event/xAOD/xAODTracking/xAODTracking/TrackParticle.h:13
xAOD::TruthParticle
TruthParticle_v1 TruthParticle
Typedef to implementation.
Definition Event/xAOD/xAODTruth/xAODTruth/TruthParticle.h:15
xAOD::TrackParticleContainer
TrackParticleContainer_v1 TrackParticleContainer
Definition of the current "TrackParticle container version".
Definition Event/xAOD/xAODTracking/xAODTracking/TrackParticleContainer.h:14
xAOD::expectInnermostPixelLayerHit
@ expectInnermostPixelLayerHit
Do we expect a 0th-layer barrel hit for this track?
Definition TrackingPrimitives.h:237
xAOD::numberOfPixelHoles
@ numberOfPixelHoles
number of pixel layers on track with absence of hits [unit8_t].
Definition TrackingPrimitives.h:262
xAOD::numberOfTRTHits
@ numberOfTRTHits
number of TRT hits [unit8_t].
Definition TrackingPrimitives.h:276
xAOD::numberOfSCTHits
@ numberOfSCTHits
number of hits in SCT [unit8_t].
Definition TrackingPrimitives.h:269
xAOD::numberOfInnermostPixelLayerHits
@ numberOfInnermostPixelLayerHits
these are the hits in the 0th pixel barrel layer
Definition TrackingPrimitives.h:238
xAOD::numberOfPixelHits
@ numberOfPixelHits
these are the pixel hits, including the b-layer [unit8_t].
Definition TrackingPrimitives.h:260
xAOD::numberOfTRTHighThresholdHits
@ numberOfTRTHighThresholdHits
number of TRT hits which pass the high threshold (only xenon counted) [unit8_t].
Definition TrackingPrimitives.h:279
xAOD::numberOfSCTHoles
@ numberOfSCTHoles
number of SCT holes [unit8_t].
Definition TrackingPrimitives.h:271
xAOD::TruthParticleContainer
TruthParticleContainer_v1 TruthParticleContainer
Declare the latest version of the truth particle container.
Definition Event/xAOD/xAODTruth/xAODTruth/TruthParticleContainer.h:17
Generated on for ATLAS Offline Software by doxygen 1.16.1