Home » Blogs » jwebb's blog

Quick ID truth Track vs Track matching

Updated on Mon, 2019-07-01 08:57. Originally created by jwebb on 2019-06-27 14:27.
Method:
 
  • Throw 200 GeV RuRu (or ZrZr) events through y2018a geometry.
  • Reconstruct FZD file using 32bit chain
  • Reconstruct resulting event.root file (hits created @ 32bit) using 64bit chain
  • Loop over all tracks in both MuDst files for first event
  • Build maps
    • MAP:              ID truth --> StMuMcTrack
    • MULTIMAP:   ID truth --> StMuTrack @ 32
    • MULTIMAP:   ID truth --> StMuTrack @ 64
  • Loop over all id truths
  • Print the number of reconstructed tracks whose dominant contributor was the given id truth...
BFC:INFO  - idtruth = 00008 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00017 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00019 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00022 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00023 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00027 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00029 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00030 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00039 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00040 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00060 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00063 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00064 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00068 ... N reco @32 = 2 ... N reco @64 = 2
BFC:INFO  - idtruth = 00070 ... N reco @32 = 2 ... N reco @64 = 1
BFC:INFO  - idtruth = 00074 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00076 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00079 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00080 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00087 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00089 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00091 ... N reco @32 = 2 ... N reco @64 = 2
BFC:INFO  - idtruth = 00092 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00102 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00113 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00114 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00115 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00117 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00121 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00125 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00131 ... N reco @32 = 2 ... N reco @64 = 1
BFC:INFO  - idtruth = 00140 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00141 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00143 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00147 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00153 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00155 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00159 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00168 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00172 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00187 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00191 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00192 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00194 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00198 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00203 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00204 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00211 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00220 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00221 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00223 ... N reco @32 = 2 ... N reco @64 = 2
BFC:INFO  - idtruth = 00228 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00229 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00236 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00239 ... N reco @32 = 4 ... N reco @64 = 4
BFC:INFO  - idtruth = 00240 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00244 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00264 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00306 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00347 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00349 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00357 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00375 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00376 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00387 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00388 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00418 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00426 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00445 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00509 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00526 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00531 ... N reco @32 = 2 ... N reco @64 = 2
BFC:INFO  - idtruth = 00536 ... N reco @32 = 1 ... N reco @64 = 0
BFC:INFO  - idtruth = 00540 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00541 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00543 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00548 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00651 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00653 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00662 ... N reco @32 = 2 ... N reco @64 = 2
BFC:INFO  - idtruth = 00681 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00699 ... N reco @32 = 3 ... N reco @64 = 3
BFC:INFO  - idtruth = 00705 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00721 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00731 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00743 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00751 ... N reco @32 = 1 ... N reco @64 = 1
BFC:INFO  - idtruth = 00799 ... N reco @32 = 1 ... N reco @64 = 1
Allows us to compare the following outcomes: 1) Same number of tracks reconstructed in 32 and 64 2) Different number of tracks reconstructed in 32 and 64 3) Track lost at 32 or 64 Outcome #2 example below...
BFC:INFO  - 32 bit track(s)
Gl id:   92 fl:  101 vx: -1 p:   0.175    0.649    0.034 q:-1 eta: 0.051 phi: 1.307 pT:  0.672 DCA [78]:37.339 -10.030 -1.096 Total hits:11 fitted:11 poss:27 Points F: 52.104 135.265  2.334 L: 48.314 172.094  4.237 idT   70 qa 100 idParentVx 1
Gl id:  104 fl:  101 vx:  0 p:   0.193    0.196    0.013 q:-1 eta: 0.048 phi: 0.794 pT:  0.275 DCA [88]:-0.429   0.457 -0.524 Total hits:17 fitted:17 poss:43 Points F: 34.162  49.559 -1.718 L: 42.763 187.195  4.661 idT   70 qa 100 idParentVx 1

BFC:INFO  - 64 bit track(s)
Gl id:   96 fl:  101 vx:  0 p:   0.193    0.197    0.013 q:-1 eta: 0.047 phi: 0.795 pT:  0.276 DCA [78]:-0.399   0.424 -0.486 Total hits:26 fitted:26 poss:42 Points F: 34.163  49.558 -1.714 L: 42.792 187.195  4.693 idT   70 qa 100 idParentVx 1
»