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

Age	Commit message (Collapse)	Author
2019-08-05	update sno_charge to prevent nan when charge is negative	tlatorre
	This commit updates the charge spectrum initialization to convolve the charge distribution over a wider range to make sure that we have a nonzero probability of observing charges all the way down to qlo. The reason for this is that previously for some negative charges the probability of observing that charge was zero, which was causing the likelihood to return nan.
2019-07-29	flag PMT charges below qlo	tlatorre
	This commit updates get_event() to flag PMT charges below qlo, which is the minimum charge value that we compute the charge PDFs for. This is to prevent the likelihood from returning nan.
2019-06-14	add a function to compute the most likely number of PE given an observed charge	tlatorre
	Also, call this function when computing the psi parameter in nll_best().
2019-03-26	small update to the charge likelihood calculation	tlatorre
	This commit updates the charge likelihood calculation to calculate: P(hit,q\|n) = P(q\|hit,n)*P(hit\|n) This has almost no effect on the fit results, but is technically correct.
2019-03-16	add GPLv3 license	tlatorre

2019-03-04	log(norm(...)) -> log_norm(...)	tlatorre

2018-11-28	update sno_charge.c	tlatorre
	This commit adds lots of comments to sno_charge.c and makes a couple of other changes: - use interp1d() instead of the GSL interpolation routines - increase MAX_PE to 100 I increased MAX_PE because I determined that it had a rather large impact on the likelihood function for 500 MeV electrons. This unfortunately slows down the initialization by a lot. I think I could speed this up by convolving the single PE charge distribution with a gaussian before convolving the charge distributions to compute the charge distributions for multiple PE.
2018-11-27	a bunch of small changes to speed things up	tlatorre

2018-09-18	add free_charge() to free memory used to interpolate the charge distributions	tlatorre

2018-09-17	fix bug in single PE charge	tlatorre
	This commit makes sure that when we conolve the single PE charge distribution with a gaussian we integrate starting at zero since the PDF is zero for q < 0.
2018-09-10	add a fast likelihood function	tlatorre
	This commit adds a fast function to calculate the expected number of PE at a PMT without numerically integrating over the track. This calculation is much faster than integrating over the track (~30 ms compared to several seconds) and so we use it during the "quick" minimization phase of the fit to quickly find the best position.
2018-09-09	fix bug in charge PDF calculation	tlatorre
	This commit fixes a bug in the charge PDF calculation for n > MAX_PE. The standard deviation should scale like sqrt(n)*qstd where qstd is the standard deviation of the single PE charge distribution.
2018-08-14	move everything to src directory	tlatorre