| Age | Commit message (Collapse) | Author | |
|---|---|---|---|
| 2020-09-07 | update plot-muons script so it works outside the utils/ directory | tlatorre | |
| 2020-09-07 | update python scripts to not call plt.show() when run with --save | tlatorre | |
| 2020-09-07 | update Makefile | tlatorre | |
| 2020-09-06 | update chi2 to correct for energy bias | tlatorre | |
| This commit updates the chi2 script to correct for the energy bias of the reconstruction relative to Monte Carlo. | |||
| 2020-09-06 | fix best fit line by changing axhline -> hlines | tlatorre | |
| 2020-09-06 | update plot-muons to fit energy scale and resolution difference | tlatorre | |
| 2020-08-31 | add a command line argument to print the nll values | tlatorre | |
| 2020-08-31 | don't call plt.show() when running with --save | tlatorre | |
| 2020-08-31 | add GENIE weights to chi2 analysis | tlatorre | |
| 2020-08-31 | add estimate_errors to chi2 analysis | tlatorre | |
| This commit updates the estimate_errors() function so that it works without a list of constraints and uses arrays of low and high bounds passed in instead of hardcoded constraints. I can now call this function from the chi2 analysis to get the stepsizes before running the MCMC. | |||
| 2020-08-31 | add a script to convert genie reweight file from ROOT -> HDF5 | tlatorre | |
| 2020-08-30 | don't sample more instrumentals than there are in the data | tlatorre | |
| 2020-08-30 | update dc scripts | tlatorre | |
| - delete instrumental tags in both scripts since we already tag them in get_events() - apply the same cuts to the MC as to the data - add a comment about the fitted fraction - quit if we don't have at least 1 type of each instrumental | |||
| 2020-08-30 | increase energy resolution uncertainty when running with --coverage | tlatorre | |
| 2020-08-30 | treat muons like Monte Carlo | tlatorre | |
| 2020-08-30 | update energy resolution -> max(EPSILON,resolution) | tlatorre | |
| This commit updates the chi2 script to make sure that the energy resolution is nonzero otherwise we get nans. | |||
| 2020-08-30 | update bins when adding energy resolution | tlatorre | |
| This commit updates the first bins that we use before applying the energy resolution in the chi2 script. The reason we do this is because when adding the energy resolution, events outside of the original bins may migrate into the first and last bin. | |||
| 2020-08-30 | small updates to chi2 script | tlatorre | |
| 2020-08-25 | fix a stupid bug in the last commit | tlatorre | |
| 2020-08-25 | speed up get_multinomial_prob() | tlatorre | |
| 2020-08-25 | fix stupid bug in last commit | tlatorre | |
| 2020-08-25 | speed up nll function in chi2 | tlatorre | |
| 2020-08-25 | update chi2 script to plot p-value coverage | tlatorre | |
| 2020-08-20 | add external muons to the chi2 fit | tlatorre | |
| 2020-08-19 | update chi2 analysis to include bias and resolution terms | tlatorre | |
| 2020-08-17 | add a new script to plot the michel energy distribution | tlatorre | |
| 2020-08-17 | update stoppig muon criteria | tlatorre | |
| 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-08-17 | simplify likelihood calculation in chi2 | tlatorre | |
| 2020-08-17 | simplify nllr() function | tlatorre | |
| 2020-08-17 | update function to tag atmospherics | tlatorre | |
| Previously the function to tag atmospherics was looking at the *first* event to come after a prompt event and checking to see if it was a neutron. However, this has a huge issue in that for large energy events, there is often secondary events caused by afterpulsing. I've now updated the algorithm to look for any events in which there is *any* follower event that passes the neutron criteria. | |||
| 2020-08-17 | fix bug in get_events() | tlatorre | |
| ev.r -> ev_single_particle.r | |||
| 2020-07-28 | only plot Michels from stopping muons | tlatorre | |
| 2020-07-28 | catch Exception when zdab-cat fails and continue | tlatorre | |
| 2020-07-27 | update plot-muons | tlatorre | |
| 2020-07-27 | update plot-muons | tlatorre | |
| - 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-27 | identify muons in the MC by looking for 'cosmic' in the filename | tlatorre | |
| 2020-07-27 | update how the negative log likelihood ratio is calculated | tlatorre | |
| I found a really simple form for the log likelihood ratio of a Poisson and multinomial likelihood. | |||
| 2020-07-27 | update how we calculte the p-value in chi2 | tlatorre | |
| 2020-07-27 | add 20 MeV cut to dc, dc-closure-test, and plot-dc | tlatorre | |
| 2020-07-27 | loop over runs in dc, dc-closure-test, and plot-dc to prevent using too much ↵ | tlatorre | |
| memory | |||
| 2020-07-27 | only merge fit info for events with at least 10 fits | tlatorre | |
| This commit updates get_events() to only merge fit info for events with at least 10 events. The reason for this is that when analyzing recent data where not all the fits have completed we don't want to plot the data for events which haven't completely finished being fit. | |||
| 2020-07-27 | update how the p-value is calculated in get_multinomial_prob() | tlatorre | |
| 2020-07-07 | add zdab-reprocess-orphans and plot-muons to Makefile | tlatorre | |
| 2020-07-07 | add a script to reprocess just junk events | tlatorre | |
| 2020-07-06 | add function to print particle probabilities in plot-muons | tlatorre | |
| 2020-07-06 | update plot-muons to plot energy distributions for each particle ID for ↵ | tlatorre | |
| stopping muons | |||
| 2020-07-06 | small updates to plot-muons and plot_energy.py | tlatorre | |
| - 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-06 | update plot-muons | tlatorre | |
| - 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-06 | add sddm/stats.py | tlatorre | |
| This commit adds the new file sddm/stats.py to and adds a function to correctly sample a Monte Carlo histogram when computing p-values. In particular, I now take into account the uncertainty on the total number of expected events by drawing from a gamma distribution, which is the posterior of the Poisson likelihood function with a prior of 1/lambda. | |||
| 2020-07-06 | only look at Michel electrons where the muon had < 2500 nhit | tlatorre | |
 |
index : sddm | |
| Code for the self-destructing dark matter search in SNO | Anthony LaTorre |
| aboutsummaryrefslogtreecommitdiff |