DUNE / dune-tms

DUNE ND Temporary Muon Spectrometer
0 stars 1 forks source link

Correct energy and time based on 3d track y information #143

Open jdkio opened 1 week ago

jdkio commented 1 week ago

Once we have track info, we can adjust the energy and time information of all the hits.

The PE in the fibers go as exp(-alpha dist), where alpha is the attenuation length. It's a bit more complex because ~half go the "long way" and are reflected. And there's also the secondary short range attentuation that isn't simulated, but may be in the future. But we can correct once we know y by multiplying the hit energy by 1/exp(-alpha reco dist) (or the more complex thing you'd use when accounting for reflection and short-range attenuation).

The time is wrong by dist/(cng), where c is the speed of light, n is the index of refraction and g is the geometric correction because the light is literally bouncing at an angle. There's also time slew that we may want to adjust for based on PE.

Note that we also want to store the uncorrected energies as well.

Bonus points for adjusting info of "nearby" hits. Like hits that are close to tracks are likely at the same or similar y position. This would require some study