Add UT to VeloSciFi matching
Description
Need lhcb-datapkg/ParamFiles!80 to be merged first to run CI/CD test in GitLab
Add the last missing piece for the Hybridseeding project.
Add UT hits for Velo-SciFi matched long tracks, the idea is inspired by HLT1 downstream: after Velo-SciFi matching, a second-order polynomial trajectory model is used to extrapolate the state from Velo to UT, open a wide search window to find the first hit, for each first hit create a long track candidate and use the position of it to correct the trajectory, finally find the rest of hits with a very narrow search window to confirm the candidate.
Since we use the first hit information to correct trajectory, the algorithm requires at least to find 2 UT hits in order to build a long track, and both the number of UT hits and the chi2 are used to evaluate a global ghost probability through a single layer neural network. This probability is then used to reject clone/ghost tracks.
The MR also modified the VeloSciFi matching (no UT) to use ghost probability to replace the chi2 comparison during the clone killing stage, it improves both the efficiency and ghost rate of the algorithm. Proposed at https://indico.cern.ch/event/1370606/#5-studies-on-hlt1-nn-ghost-kil
Performance
Validator output: Bs2PhiPhi (1k events)
forward validation:
TrackChecker output : 1630/ 44944 3.63% ghosts
for P>3GeV,Pt>0.5GeV : 906/ 30332 2.99% ghosts
01_long : 40748/ 78155 52.14% ( 53.86%), 708 ( 1.71%) clones, pur 98.75%, hit eff 98.20%
02_long_P>5GeV : 34326/ 51059 67.23% ( 69.06%), 610 ( 1.75%) clones, pur 98.88%, hit eff 98.51%
03_long_strange : 1294/ 3858 33.54% ( 34.34%), 19 ( 1.45%) clones, pur 98.34%, hit eff 97.90%
04_long_strange_P>5GeV : 966/ 1834 52.67% ( 54.22%), 15 ( 1.53%) clones, pur 98.54%, hit eff 98.33%
05_long_fromB : 3379/ 4585 73.70% ( 74.83%), 58 ( 1.69%) clones, pur 98.87%, hit eff 98.63%
06_long_fromB_P>5GeV : 3115/ 3752 83.02% ( 83.45%), 55 ( 1.74%) clones, pur 98.94%, hit eff 98.81%
07_long_electrons : 799/ 5658 14.12% ( 14.70%), 27 ( 3.27%) clones, pur 97.55%, hit eff 97.79%
08_long_electrons_P>5GeV : 694/ 2953 23.50% ( 24.71%), 22 ( 3.07%) clones, pur 97.43%, hit eff 97.87%
09_long_fromB_electrons : 81/ 240 33.75% ( 37.55%), 2 ( 2.41%) clones, pur 98.34%, hit eff 98.21%
10_long_fromB_electrons_P>5GeV : 76/ 168 45.24% ( 48.69%), 1 ( 1.30%) clones, pur 98.21%, hit eff 98.40%
long_P>5GeV_AND_Pt>1GeV : 11210/ 13336 84.06% ( 86.34%), 171 ( 1.50%) clones, pur 98.83%, hit eff 98.66%
long_fromB_P>5GeV_AND_Pt>1GeV : 2168/ 2500 86.72% ( 86.58%), 35 ( 1.59%) clones, pur 98.88%, hit eff 98.93%
11_noVelo_UT : 0/ 8699 0.00% ( 0.00%), 0 ( 0.00%) clones, pur -nan%, hit eff -nan%
12_noVelo_UT_P>5GeV : 0/ 3556 0.00% ( 0.00%), 0 ( 0.00%) clones, pur -nan%, hit eff -nan%
13_long_PT>2GeV : 2809/ 3268 85.95% ( 86.75%), 38 ( 1.33%) clones, pur 98.85%, hit eff 98.90%
14_long_from_B_PT>2GeV : 1017/ 1172 86.77% ( 86.65%), 8 ( 0.78%) clones, pur 98.86%, hit eff 99.17%
15_long_strange_P>5GeV : 966/ 1834 52.67% ( 54.22%), 15 ( 1.53%) clones, pur 98.54%, hit eff 98.33%
16_long_strange_P>5GeV_PT>500MeV : 590/ 738 79.95% ( 80.53%), 7 ( 1.17%) clones, pur 98.60%, hit eff 98.46%
17_long_fromSignal : 2162/ 2621 82.49% ( 81.75%), 37 ( 1.68%) clones, pur 98.85%, hit eff 98.73%
trackmatching validation:
TrackChecker output : 2127/ 53303 3.99% ghosts
for P>3GeV,Pt>0.5GeV : 498/ 28548 1.74% ghosts
01_long : 48489/ 78155 62.04% ( 65.30%), 0 ( 0.00%) clones, pur 98.55%, hit eff 98.33%
02_long_P>5GeV : 38974/ 51059 76.33% ( 80.64%), 0 ( 0.00%) clones, pur 98.65%, hit eff 98.61%
03_long_strange : 1942/ 3858 50.34% ( 52.76%), 0 ( 0.00%) clones, pur 97.97%, hit eff 98.17%
04_long_strange_P>5GeV : 1309/ 1834 71.37% ( 74.85%), 0 ( 0.00%) clones, pur 98.23%, hit eff 98.56%
05_long_fromB : 3513/ 4585 76.62% ( 77.43%), 0 ( 0.00%) clones, pur 98.88%, hit eff 98.71%
06_long_fromB_P>5GeV : 3201/ 3752 85.31% ( 85.63%), 0 ( 0.00%) clones, pur 98.92%, hit eff 98.89%
07_long_electrons : 1526/ 5658 26.97% ( 28.57%), 0 ( 0.00%) clones, pur 96.49%, hit eff 98.43%
08_long_electrons_P>5GeV : 1188/ 2953 40.23% ( 42.55%), 0 ( 0.00%) clones, pur 96.55%, hit eff 98.54%
09_long_fromB_electrons : 107/ 240 44.58% ( 48.17%), 0 ( 0.00%) clones, pur 98.20%, hit eff 99.02%
10_long_fromB_electrons_P>5GeV : 101/ 168 60.12% ( 63.69%), 0 ( 0.00%) clones, pur 98.28%, hit eff 99.04%
long_P>5GeV_AND_Pt>1GeV : 10845/ 13336 81.32% ( 85.17%), 0 ( 0.00%) clones, pur 98.86%, hit eff 98.72%
long_fromB_P>5GeV_AND_Pt>1GeV : 2181/ 2500 87.24% ( 87.08%), 0 ( 0.00%) clones, pur 98.92%, hit eff 98.96%
11_noVelo_UT : 0/ 8699 0.00% ( 0.00%), 0 ( 0.00%) clones, pur -nan%, hit eff -nan%
12_noVelo_UT_P>5GeV : 0/ 3556 0.00% ( 0.00%), 0 ( 0.00%) clones, pur -nan%, hit eff -nan%
13_long_PT>2GeV : 2767/ 3268 84.67% ( 86.45%), 0 ( 0.00%) clones, pur 98.87%, hit eff 98.97%
14_long_from_B_PT>2GeV : 1054/ 1172 89.93% ( 89.31%), 0 ( 0.00%) clones, pur 98.87%, hit eff 99.08%
15_long_strange_P>5GeV : 1309/ 1834 71.37% ( 74.85%), 0 ( 0.00%) clones, pur 98.23%, hit eff 98.56%
16_long_strange_P>5GeV_PT>500MeV : 590/ 738 79.95% ( 81.60%), 0 ( 0.00%) clones, pur 98.66%, hit eff 98.75%
17_long_fromSignal : 2182/ 2621 83.25% ( 83.09%), 0 ( 0.00%) clones, pur 98.92%, hit eff 98.82%Validator plots: Bs2PhiPhi (10k events)
| Efficiency | Ghost rate |
|---|---|
 |
 |
 |
 |
 |
 |
 |
 |
Comparison between matching (no UT) and matching (with UT)
| Efficiency | Ghost rate |
|---|---|
 |
 |
 |
 |
 |
 |
 |
 |
Note: As we expect to see, adding UT hits to matching should have small impact in efficiency and great improvement in ghost rate
Throughput
The throughput is measured with RTX 4070Ti (faster than 3090) and then scaled down to the expected A5000 throughput.
FYI: @adeoyang @bjashal @cagapopo @dovombru @mveghel @lohenry