| Age | Commit message (Collapse) | Author | 
 | 
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.
 | 
 | 
 | 
 | 
 | 
 | 
- 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
 | 
 | 
 | 
 | 
 | 
 | 
stopping muons
 | 
 | 
- 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
 | 
 | 
- 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
 | 
 | 
 | 
 | 
- only look at muons with nhit < 4000 and udotr < -0.5
- switch from energy1 -> ke
 | 
 | 
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.
 |