sddm - Code for the self-destructing dark matter search in SNO

Age	Commit message (Collapse)	Author
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-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	simplify likelihood calculation in chi2	tlatorre

2020-07-27	update how we calculte the p-value in chi2	tlatorre

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-06-29	add energy scale to chi2 fit	tlatorre

2020-06-23	first draft of a script to do the chi2 analysis	tlatorre
	This commit adds a first draft of a script called chi2. This script calculates a chi2 for the null hypothesis test if the events in the energy range 20 MeV - 10 GeV match what we expect from atmospheric neutrino events.