Kalman filter - resolution


Previous: LDRD occupancies   |   Index   |   Next: material budget


Throw 5 Gev muons into the N arm:
cfmgmc.input = (5. 5.0 18.0 90.0 14.0 100 0.0 0.0 0.0 0.0)


From rcas cvs3 files fvtx_mc_eval_20perc4.root, fvtx_mc_eval_50100_4.root, fvtx_mc_eval_perf4.root.

x0reco   track projection to the z=0 plane
y0reco   track projection to the z=0 plane
z0reco   z=0 or to z-axis


5.0 GeV μ x0reco y0reco z0reco
mean sigma mean sigma mean sigma
pi = poriginal ~0 47 ~0 54 - -
pi->pi±20% ~0 131 ~0 71 - -
px=py=50, pz=100 0.08 228 ~0 354 - -
straight-line fit ~0 231 ~0 356 ~0 166


Note that with perfect information, the Kalman filter projects to %50μm, which is roughly what we expect from multiple scattering.

On the other extreme, when no momentum information is available, and the best we can do is assume a straight line by feeding in very high momenta, the Kalman filter returns the same values as a straight-line fit.

The realistic case is one where we have partial knowledge of the momentum, simulated by smearing the input momentum by 20%.





Macros: xres_vs_p.C.txt, yres_vs_p.C.txt

2.5 GeV μ x0reco y0reco z0reco
mean sigma mean sigma mean sigma
pi = poriginal ~0 95 ~0 102 - -
pi->pi±20% ~0 202 ~0 121 - -
px=py=50, pz=100 0.14cm 478 ~0 714 - -
straight-line fit 0.14 476 ~0 720 ~0 306


5.0 GeV μ x0reco y0reco z0reco
mean sigma mean sigma mean sigma
pi = poriginal ~0 47 ~0 54 - -
pi->pi±20% ~0 131 ~0 71 - -
px=py=50, pz=100 0.08 228 ~0 354 - -
straight-line fit ~0 231 ~0 356 ~0 166


10 GeV μ x0reco y0reco z0reco
mean sigma mean sigma mean sigma
pi = poriginal ~0 23 ~0 26 - -
pi->pi±20% ~0 67 ~0 32 - -
px=py=50, pz=100 0.03 96 ~0 160 - -
straight-line fit 0.03 104 ~0 160 ~0 86


20 GeV μ x0reco y0reco z0reco
mean sigma mean sigma mean sigma
pi = poriginal ~0 13 ~0 13 - -
pi->pi±20% ~0 33 ~0 18 - -
px=py=50, pz=100 0.015 52 ~0 69 - -
straight-line fit 0.018 56 ~0 81 ~0 36




Previous: LDRD occupancies   |   Index   |   Next: material budget