Skip to content
Draft
Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
71 commits
Select commit Hold shift + click to select a range
b1acf6e
OpenKF for our case is building and predicting
anthoak13 Jun 9, 2025
c507853
Implement update portion of UKF test
anthoak13 Jun 9, 2025
f0f9ee2
Remove debug statements
anthoak13 Jun 9, 2025
0909519
Add check for positive definite cov matrix
anthoak13 Jun 9, 2025
87587b2
Add physics test cases to be filled
Jun 10, 2025
6b8698b
Run clang format
Jun 10, 2025
9efdcfc
Force is working and propogator with no fields is working
anthoak13 Jul 17, 2025
abae9bc
Add additional test
anthoak13 Jul 17, 2025
a7036b6
Add test for fixing energy loss
anthoak13 Jul 17, 2025
37fc28b
Move some code from tests to AtTools
anthoak13 Jul 17, 2025
b0bed98
Add Propagator class for charged particles in mat
anthoak13 Jul 21, 2025
72e8087
Update to using XYZPoint over XYZVector
anthoak13 Jul 21, 2025
27ca184
Use internal scaling factor
anthoak13 Jul 21, 2025
ef1dec0
Fixed step propagator stopping in correct place
anthoak13 Jul 21, 2025
57ea952
Non-Adaptive tests working
anthoak13 Jul 21, 2025
fdee72a
Adaptive working
anthoak13 Jul 21, 2025
8b38a50
Add test for small initial step size
anthoak13 Jul 21, 2025
6cda995
Fix unit issue in stopping plane calculation.
anthoak13 Jul 22, 2025
511e04f
PropagateToPoint uses AtStepper interface
anthoak13 Jul 22, 2025
a848c2b
All RK4 to point tests use stepper interface
anthoak13 Jul 22, 2025
08c7ccb
Adaptive test working with stepper
anthoak13 Jul 22, 2025
365e211
Plane propagation uses AtStepper interface
anthoak13 Jul 22, 2025
1d563a7
Propogate to point using common surface formalism
anthoak13 Jul 22, 2025
dcccdfd
Migrate fixed eloss to new interfaces
anthoak13 Jul 22, 2025
df4ec30
Plans using interface
anthoak13 Jul 22, 2025
c847df1
Refactor and only project if at surface
anthoak13 Jul 22, 2025
1a8bdcc
Refactor state struct
anthoak13 Jul 22, 2025
ec81eb1
Refactor to passing around state struct
anthoak13 Jul 22, 2025
7332716
Remove double counted variables from state
anthoak13 Jul 22, 2025
03e6d77
Remove debug and more refactoring
anthoak13 Jul 22, 2025
4801616
Refactor step size defaults
anthoak13 Jul 22, 2025
9d06ade
Fix comments that implied AtEnergyLoss density was mg/cm^3
anthoak13 Jul 22, 2025
caffb8b
Possibly make tests portable.
anthoak13 Jul 23, 2025
2d24c7a
Try to setup env in github runner
anthoak13 Jul 23, 2025
0f009dc
Second attempt at setting CI env
anthoak13 Jul 23, 2025
8a29bd0
Try 3 - didn't run at all
anthoak13 Jul 23, 2025
aacb527
Probably really works
anthoak13 Jul 23, 2025
9b41edb
Add integration test for drawing figures
anthoak13 Jul 23, 2025
b3efecb
Remove old UKF tests in prep for real first test
anthoak13 Jul 24, 2025
a711861
Reset KF changes
anthoak13 Jul 24, 2025
6f74ffd
Original tests building again. Need fresh slate
anthoak13 Jul 24, 2025
b288edd
Base class with virtual tests pass
anthoak13 Jul 24, 2025
73f294b
UKF w/ prop is instantiating
anthoak13 Jul 24, 2025
8d67056
Prediction code runs
anthoak13 Jul 24, 2025
9573b64
Seemingly working forward step
anthoak13 Jul 24, 2025
5e8471c
Add full test track. Decomp failing
anthoak13 Jul 26, 2025
479bb84
Fix decomposition, but UKF seems suspicious.
anthoak13 Jul 26, 2025
ab3b642
Fix issue with initial COV being too high. Filter not follows measure…
anthoak13 Jul 26, 2025
3c516c2
Move most propagator logging to debug
anthoak13 Jul 28, 2025
2b070b9
Fix bug in straggling.
anthoak13 Jul 28, 2025
796d893
Cleanup UKF doc and headers
anthoak13 Jul 28, 2025
bd62e90
Save the information needed for smoothing
anthoak13 Jul 29, 2025
322e738
Refactor LLT decomp into function
anthoak13 Jul 29, 2025
841a7cd
Adjust weights to use method that should protect PD of COV better
anthoak13 Jul 29, 2025
d73785f
Add function to ensure PD and only save PD matrices
anthoak13 Jul 29, 2025
22c624e
Remove debug. Add diagnostic plots
anthoak13 Jul 29, 2025
56bc621
Parameter tuning
anthoak13 Jul 29, 2025
06919b8
Add ability to build CATIMA offline if source is already downloaded.
anthoak13 Jul 29, 2025
520786f
Increase stability of filter
anthoak13 Jul 29, 2025
d9247ca
Optimize Kalman gain calculation
anthoak13 Jul 29, 2025
42994ef
Finish optimizing of smoothing gain
anthoak13 Jul 29, 2025
075ab92
Zero missing parts of filter
anthoak13 Jul 30, 2025
534bfe1
Add test for many tracks
anthoak13 Jul 30, 2025
59f4cb4
Add diagnostic and add COV correction
anthoak13 Jul 30, 2025
596d81b
Add controllable diagnostics
anthoak13 Jul 30, 2025
f14be20
Rename noises
anthoak13 Jul 30, 2025
e94e3f0
More refactoring
anthoak13 Jul 30, 2025
ce22ed3
More refactoring
anthoak13 Jul 30, 2025
75ab272
Remove inverse calls
anthoak13 Jul 30, 2025
4caa601
Incorporate model noise
anthoak13 Jul 30, 2025
a078937
Add more output to many track test macro
anthoak13 Jul 30, 2025
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
143 changes: 143 additions & 0 deletions AtReconstruction/AtFitter/OpenKF/kalman_filter/HinH.txt
Original file line number Diff line number Diff line change
@@ -0,0 +1,143 @@
==================================================================
SRIM version ---> SRIM-2013.00
Calc. date ---> julio 16. 2025
==================================================================

Disk File Name = SRIM Outputs\Hydrogen in Hydrogen (gas).txt

Ion = Hydrogen [1] . Mass = 1.008 amu

Target Density = 4.1906E-05 g/cm3 = 2.5036E+19 atoms/cm3
Target is a GAS
======= Target Composition ========
Atom Atom Atomic Mass
Name Numb Percent Percent
---- ---- ------- -------
H 1 100.00 100.00
====================================
Bragg Correction = 0.00%
Stopping Units = MeV / (mg/cm2)
See bottom of Table for other Stopping units

Ion dE/dx dE/dx Projected Longitudinal Lateral
Energy Elec. Nuclear Range Straggling Straggling
-------------- ---------- ---------- ---------- ---------- ----------
999.999 eV 7.504E-01 1.692E-01 267.71 um 61.83 um 55.80 um
1.10 keV 7.871E-01 1.600E-01 291.11 um 64.85 um 59.45 um
1.20 keV 8.221E-01 1.519E-01 313.96 um 67.62 um 62.90 um
1.30 keV 8.556E-01 1.447E-01 336.31 um 70.18 um 66.16 um
1.40 keV 8.879E-01 1.383E-01 358.18 um 72.56 um 69.25 um
1.50 keV 9.191E-01 1.325E-01 379.59 um 74.77 um 72.19 um
1.60 keV 9.492E-01 1.272E-01 400.57 um 76.83 um 74.99 um
1.70 keV 9.785E-01 1.224E-01 421.14 um 78.76 um 77.65 um
1.80 keV 1.007E+00 1.180E-01 441.32 um 80.57 um 80.20 um
2.00 keV 1.061E+00 1.101E-01 480.60 um 83.97 um 84.97 um
2.25 keV 1.124E+00 1.019E-01 527.87 um 87.75 um 90.40 um
2.50 keV 1.183E+00 9.491E-02 573.34 um 91.09 um 95.34 um
2.75 keV 1.239E+00 8.894E-02 617.19 um 94.06 um 99.86 um
3.00 keV 1.293E+00 8.377E-02 659.59 um 96.72 um 104.03 um
3.25 keV 1.344E+00 7.923E-02 700.68 um 99.13 um 107.88 um
3.50 keV 1.393E+00 7.521E-02 740.56 um 101.32 um 111.46 um
3.75 keV 1.440E+00 7.162E-02 779.35 um 103.33 um 114.80 um
4.00 keV 1.486E+00 6.840E-02 817.12 um 105.18 um 117.93 um
4.50 keV 1.572E+00 6.285E-02 889.92 um 108.70 um 123.65 um
5.00 keV 1.654E+00 5.821E-02 959.49 um 111.75 um 128.76 um
5.50 keV 1.731E+00 5.428E-02 1.03 mm 114.44 um 133.36 um
6.00 keV 1.804E+00 5.090E-02 1.09 mm 116.83 um 137.55 um
6.50 keV 1.873E+00 4.795E-02 1.15 mm 118.97 um 141.38 um
7.00 keV 1.939E+00 4.536E-02 1.21 mm 120.92 um 144.91 um
8.00 keV 2.064E+00 4.101E-02 1.33 mm 124.82 um 151.22 um
9.00 keV 2.178E+00 3.748E-02 1.44 mm 128.13 um 156.72 um
10.00 keV 2.285E+00 3.457E-02 1.54 mm 131.01 um 161.59 um
11.00 keV 2.384E+00 3.211E-02 1.64 mm 133.54 um 165.95 um
12.00 keV 2.477E+00 3.001E-02 1.74 mm 135.79 um 169.89 um
13.00 keV 2.564E+00 2.819E-02 1.83 mm 137.82 um 173.48 um
14.00 keV 2.646E+00 2.660E-02 1.92 mm 139.67 um 176.77 um
15.00 keV 2.723E+00 2.519E-02 2.01 mm 141.36 um 179.82 um
16.00 keV 2.796E+00 2.393E-02 2.09 mm 142.92 um 182.65 um
17.00 keV 2.864E+00 2.281E-02 2.17 mm 144.37 um 185.29 um
18.00 keV 2.929E+00 2.179E-02 2.26 mm 145.72 um 187.77 um
20.00 keV 3.048E+00 2.003E-02 2.41 mm 149.01 um 192.31 um
22.50 keV 3.179E+00 1.822E-02 2.60 mm 153.11 um 197.34 um
25.00 keV 3.292E+00 1.673E-02 2.79 mm 156.75 um 201.82 um
27.50 keV 3.389E+00 1.549E-02 2.96 mm 160.04 um 205.87 um
30.00 keV 3.473E+00 1.443E-02 3.14 mm 163.06 um 209.56 um
32.50 keV 3.544E+00 1.351E-02 3.30 mm 165.87 um 212.98 um
35.00 keV 3.604E+00 1.272E-02 3.47 mm 168.50 um 216.15 um
37.50 keV 3.654E+00 1.201E-02 3.63 mm 170.98 um 219.13 um
40.00 keV 3.695E+00 1.139E-02 3.79 mm 173.35 um 221.95 um
45.00 keV 3.753E+00 1.033E-02 4.11 mm 180.64 um 227.18 um
50.00 keV 3.784E+00 9.466E-03 4.43 mm 187.42 um 231.99 um
55.00 keV 3.794E+00 8.743E-03 4.74 mm 193.85 um 236.49 um
60.00 keV 3.785E+00 8.129E-03 5.06 mm 200.03 um 240.76 um
65.00 keV 3.763E+00 7.601E-03 5.37 mm 206.04 um 244.84 um
70.00 keV 3.730E+00 7.142E-03 5.69 mm 211.93 um 248.78 um
80.00 keV 3.638E+00 6.382E-03 6.33 mm 232.69 um 256.39 um
90.00 keV 3.526E+00 5.776E-03 7.00 mm 252.74 um 263.79 um
100.00 keV 3.403E+00 5.282E-03 7.69 mm 272.49 um 271.10 um
110.00 keV 3.276E+00 4.870E-03 8.40 mm 292.22 um 278.42 um
120.00 keV 3.150E+00 4.522E-03 9.14 mm 312.12 um 285.84 um
130.00 keV 3.026E+00 4.222E-03 9.91 mm 332.30 um 293.41 um
140.00 keV 2.907E+00 3.963E-03 10.72 mm 352.84 um 301.18 um
150.00 keV 2.794E+00 3.735E-03 11.55 mm 373.81 um 309.18 um
160.00 keV 2.687E+00 3.533E-03 12.42 mm 395.24 um 317.44 um
170.00 keV 2.587E+00 3.353E-03 13.33 mm 417.16 um 326.00 um
180.00 keV 2.492E+00 3.192E-03 14.26 mm 439.58 um 334.87 um
200.00 keV 2.320E+00 2.915E-03 16.25 mm 524.78 um 353.64 um
225.00 keV 2.134E+00 2.632E-03 18.92 mm 651.26 um 379.16 um
250.00 keV 1.975E+00 2.402E-03 21.83 mm 773.93 um 407.12 um
275.00 keV 1.840E+00 2.211E-03 24.95 mm 895.52 um 437.62 um
300.00 keV 1.722E+00 2.050E-03 28.30 mm 1.02 mm 470.70 um
325.00 keV 1.620E+00 1.911E-03 31.87 mm 1.14 mm 506.38 um
350.00 keV 1.530E+00 1.791E-03 35.65 mm 1.26 mm 544.65 um
375.00 keV 1.451E+00 1.686E-03 39.65 mm 1.39 mm 585.49 um
400.00 keV 1.381E+00 1.594E-03 43.86 mm 1.52 mm 628.85 um
450.00 keV 1.260E+00 1.437E-03 52.90 mm 1.99 mm 723.02 um
500.00 keV 1.162E+00 1.310E-03 62.75 mm 2.44 mm 826.75 um
550.00 keV 1.079E+00 1.205E-03 73.40 mm 2.88 mm 939.68 um
600.00 keV 1.008E+00 1.116E-03 84.83 mm 3.31 mm 1.06 mm
650.00 keV 9.477E-01 1.039E-03 97.02 mm 3.75 mm 1.19 mm
700.00 keV 8.947E-01 9.735E-04 109.96 mm 4.18 mm 1.33 mm
800.00 keV 8.064E-01 8.650E-04 138.03 mm 5.79 mm 1.63 mm
900.00 keV 7.355E-01 7.791E-04 169.00 mm 7.29 mm 1.96 mm
1.00 MeV 6.772E-01 7.095E-04 202.78 mm 8.74 mm 2.32 mm
1.10 MeV 6.322E-01 6.518E-04 239.22 mm 10.17 mm 2.71 mm
1.20 MeV 5.879E-01 6.031E-04 278.33 mm 11.61 mm 3.13 mm
1.30 MeV 5.513E-01 5.616E-04 320.21 mm 13.07 mm 3.57 mm
1.40 MeV 5.197E-01 5.256E-04 364.75 mm 14.54 mm 4.05 mm
1.50 MeV 4.918E-01 4.942E-04 411.91 mm 16.04 mm 4.54 mm
1.60 MeV 4.670E-01 4.665E-04 461.66 mm 17.55 mm 5.07 mm
1.70 MeV 4.448E-01 4.418E-04 513.97 mm 19.09 mm 5.62 mm
1.80 MeV 4.247E-01 4.197E-04 568.83 mm 20.65 mm 6.19 mm
2.00 MeV 3.899E-01 3.819E-04 686.03 mm 26.57 mm 7.41 mm
2.25 MeV 3.542E-01 3.436E-04 846.45 mm 35.05 mm 9.09 mm
2.50 MeV 3.250E-01 3.125E-04 1.02 m 43.08 mm 10.91 mm
2.75 MeV 3.005E-01 2.868E-04 1.21 m 50.96 mm 12.88 mm
3.00 MeV 2.797E-01 2.652E-04 1.42 m 58.81 mm 14.99 mm
3.25 MeV 2.618E-01 2.467E-04 1.64 m 66.70 mm 17.25 mm
3.50 MeV 2.462E-01 2.307E-04 1.87 m 74.66 mm 19.66 mm
3.75 MeV 2.325E-01 2.168E-04 2.12 m 82.73 mm 22.20 mm
4.00 MeV 2.203E-01 2.045E-04 2.39 m 90.91 mm 24.88 mm
4.50 MeV 1.997E-01 1.838E-04 2.95 m 121.70 mm 30.66 mm
5.00 MeV 1.829E-01 1.671E-04 3.58 m 150.69 mm 36.98 mm
5.50 MeV 1.689E-01 1.533E-04 4.26 m 179.04 mm 43.83 mm
6.00 MeV 1.570E-01 1.416E-04 4.99 m 207.27 mm 51.20 mm
6.50 MeV 1.468E-01 1.317E-04 5.78 m 235.63 mm 59.08 mm
7.00 MeV 1.379E-01 1.231E-04 6.61 m 264.26 mm 67.48 mm
8.00 MeV 1.232E-01 1.090E-04 8.44 m 370.83 mm 85.77 mm
9.00 MeV 1.116E-01 9.794E-05 10.48 m 470.37 mm 106.02 mm
10.00 MeV 1.021E-01 8.897E-05 12.71 m 567.83 mm 128.19 mm
11.00 MeV 9.416E-02 8.155E-05 15.15 m 665.14 mm 152.25 mm
-----------------------------------------------------------
Multiply Stopping by for Stopping Units
------------------- ------------------
4.1905E-04 eV / Angstrom
4.1905E-03 keV / micron
4.1905E-03 MeV / mm
1.0000E+00 keV / (ug/cm2)
1.0000E+00 MeV / (mg/cm2)
1.0000E+03 keV / (mg/cm2)
1.6738E+00 eV / (1E15 atoms/cm2)
5.5947E-01 L.S.S. reduced units
==================================================================
(C) 1984.1989.1992.1998.2008 by J.P. Biersack and J.F. Ziegler
Loading
Loading