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