Meeting ID: 958 9367 1004
For MC-hits, I use Poisson hit probability, within 80 m from the track.
For real hits, use both hit and no-hit information.
I use two tests to see whether there is a Bremsstrahlung shower.
Poisson-test shows 5 Bremsstrahlung shower candidates in 21 events.
If I scan over a range of shower distances, the likelihood contains a lot of local minima.
I'm investigating whether the scipy-optimizer finds the global minimum correctly here.
If I fix the distance, but vary the energy, the likelihood looks much better.
For the real hits, the time-residuals are not very well approximated by the PDF.
Usually, I see less hits than I would expect around a time-residual of 0.
Ronald: Can you look at whether the extra spread in the time residuals derives from the amount of misreconstruction?
Ronald: Are you constructing the PDF from the reco-tracks or MC-truth tracks?
Ronald: Using the reco-tracks may bias the PDF.
Maarten: The reconstruction looks for the optimal track, so you don't expect a perfect match here. Additionally, your hypothesis may be more perfect than reality here. The events you show seem not very well reconstructed and maybe also do not contain as perfect a shower as ...
If your procedure can be used to pick out events with a huge stochastic energy loss, this would be extremely cool and we can use your tool to
Would love to see identification of events with 2 or more showers and plot the energy ratio of the first and the last shower. There should be a bias there.
Ronald: You can actually fit the energy loss.
A small addition to what I showed yesterday:
Deviation around x-position of vertex gives very sharp likelihood minimum.
Maarten: What do you expect? Shower fit angular resolution is several dozen of degrees. So an improvement here is probably not expected.
Jordan: You can maybe estimate the resolution from the width of the minimum and compare to the angular resolution you expect.
Ronald: The x-, y- and z-axis are fixed to the detector geometry here, so these plots are sometimes difficult to interpret.
Maarten: At these energies, the only information that depends on the shower direction is the hit- no-hit info.
Can we do a signal+background and background only hypothesis and take the difference?
Then the number of no-hit PMTs included should not matter a lot anymore.
Ronald: One of the issues here is that we may be including a large amount of hits also from the hadronic shower in the nue-CC events. These should be removed to make a clear comparison.
Thijs: In energy-deviation right-side plot, we see the minimum located at -10 GeV. If you would have a hadronic shower that was not taken into account, we would expect the minimum to be at positive energy here.
Ronald: That depends upon how the deviation is defined here.
Jordan: Might also be a good idea to include energy-elongation.