aboutsummaryrefslogtreecommitdiff
path: root/utils/plot-muons
AgeCommit message (Collapse)Author
2020-09-07update python scripts to not call plt.show() when run with --savetlatorre
2020-09-06fix best fit line by changing axhline -> hlinestlatorre
2020-09-06update plot-muons to fit energy scale and resolution differencetlatorre
2020-08-17update stoppig muon criteriatlatorre
This commit updates the criteria for selecting stopping muons from: - calibrated nhit < 4000 - udotr < -0.5 to - reconstructed kinetic energy < 10 GeV The previous criteria were intended to remove through going atmospheric events but produced a strong bias in the comparison due to the nhit cut and an energy bias in the data relative to the Monte Carlo. The new cut does a good job of cutting through going muons but doesn't produce the same bias.
2020-07-28only plot Michels from stopping muonstlatorre
2020-07-27update plot-muonstlatorre
2020-07-27update plot-muonstlatorre
- added a cos(theta) cut - plot the energy and angular distribution of stopping muons - fix bug in calculating Michel normalization constant - plot legend for energy resolution plot
2020-07-27identify muons in the MC by looking for 'cosmic' in the filenametlatorre
2020-07-06add function to print particle probabilities in plot-muonstlatorre
2020-07-06update plot-muons to plot energy distributions for each particle ID for ↵tlatorre
stopping muons
2020-07-06small updates to plot-muons and plot_energy.pytlatorre
- use pd.Series.where() instead of DataFrame.loc() to speed things up in tag_michels - don't set y limits when plotting bias and resolution for stopping muons
2020-07-06update plot-muonstlatorre
- add get_multinomial_prob() function to stats.py - add plot_hist2_data_mc() function to do the normal particle id plot but also print p values - other small bug fixes
2020-07-06only look at Michel electrons where the muon had < 2500 nhittlatorre
2020-07-06update plot-muonstlatorre
- only look at muons with nhit < 4000 and udotr < -0.5 - switch from energy1 -> ke
2020-07-06add first draft of plot-muonstlatorre
This commit adds a first draft of a script to plot the michel energy distribution and particle id histograms for data and Monte Carlo and to plot the energy bias and resolution for stopping muons.