Muon Energy Resolution

The code lives in the EnergyResolution subdir in the kleykamp_validation branch of dune_tms. It takes in the validation plots, described here

Then take the output file and call:

python fit_resolution.py /exp/dune/data/users/kleykamp/dune-tms/Validation/Tracking_Validation/2025-04-10_stereo.root

It will do a bunch of fits, and then output the linear fit function.

#{Chi}^2/dof: 5.688	ndf: 10	chi2/ndf: 0.5688	mean: -0.07 +/- 0.45	sigma: 0.03 +/- 1.11	0
Fit function:
#define length_to_energy(l) (1.75*l*1.0020 + -0.7820)*1e-3

Example output plots: /exp/dune/app/users/kleykamp/tms_validation/dune-tms/scripts/EnergyResolution/2025-04-10_stereo

This function then can be placed in the validation scripts as the correct function for this sample (changes per geometry and with reco changes). Add it to the function, and assign the right enum based on your naming scheme here

Tech Note

Here's the energy resolution tech note

Location of plots

The code is here. It reco/true in slices of true energy, along with diagnostic plots.

Outputs to <basedir>/energy_resolution