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

Age	Commit message (Collapse)	Author
2018-10-21	set default number of points to 500	tlatorre

2018-10-19	epsrel -> npoints	tlatorre

2018-10-18	update fit to fit for electrons and protons	tlatorre

2018-10-12	skip PMTs which weren't hit for the fast likelihood calculation	tlatorre

2018-10-06	prune output of mcvx and mctk banks	tlatorre

2018-10-05	epsrel = 1e-2 to speed up fit	tlatorre

2018-10-02	update optics to use LETA constants for the D2O phase	tlatorre
	This commit updates the PMT response bank and absorption lengths for H2O and D2O based on the values in mcprod which were used for the LETA analysis.
2018-10-01	flag PMTs with the KPF_DIS bit set	tlatorre
	Doh! I was previously using the tube number variable!
2018-09-25	increase maxeval to 20 for the "quick" minimization phase	tlatorre
	When testing out the fitter on 500 MeV muons, there was at least one event which started out at a position very far from the true position. This event had a secondary electron like ring which is what I suspect caused the fit to start out in a position far from the true position. This fix correctly starts the minimization close to the true position. In the future I should look at updating get_direction() so that it finds the largest ring direction instead of just doing a weighted average of all the vectors from the position to the PMTs.
2018-09-20	add particle id code to output file	tlatorre

2018-09-20	add git SHA1 hash to output file	tlatorre

2018-09-20	add time elapsed to the output file	tlatorre

2018-09-20	add a command line option to only fit the first event	tlatorre
	This commit adds the ability to run the fit program with the --skip-second-event flag to only fit the first event after a MAST bank. This way we avoid fitting secondaries like Michel electrons when fitting MC events.
2018-09-20	add absorption lengths for D2O and H2O weighted by the Cerenkov spectrum and ↵	tlatorre
	the quantum efficiency
2018-09-19	change output file format to YAML	tlatorre

2018-09-18	stop fitting when the likelihood difference is less than 1e-2	tlatorre

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

2018-09-18	add a ctrl-c handler	tlatorre

2018-09-13	fix width formatting for t0	tlatorre

2018-09-12	small updates to speed things up	tlatorre

2018-09-12	update the starting parameters	tlatorre

2018-09-11	only print the likelihood value once for each iteration during the "quick" ↵	tlatorre
	minimization phase
2018-09-11	fix the energy and position when doing the "quick" minimizations	tlatorre

2018-09-11	switch order of expressions to avoid a valgrind warning	tlatorre

2018-09-11	add PMT response	tlatorre
	This commit adds code to read in the PMT response from the PMTR bank from SNOMAN. This file was used for the grey disk model in SNOMAN and was created using a full 3D simulation of the PMT and concentrator. Since the PMT response in SNOMAN included the quantum efficiency of the PMT, we have to divide that out to get just the PMT response independent of the quantum efficiency. I also updated the likelihood calculation to use the pmt response. Currently the energy is being fit too high which I think will improve when we update the solid angle calculation to use the radius of the concentrator instead of the PMT.
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-04	update kinetic energy step size to 2% of initial kinetic energy guess	tlatorre

2018-09-04	update fit to guess energy, direction, and t0	tlatorre
	This commit updates the initial guess for the energy using a simple heuristic of ~6 hits/MeV. I also updated the initial phase where we do a bunch of "quick" minimizations to loop over a series of starting positions and automatically calculate the approximate direction and t0 for the event.
2018-08-31	start by doing a series of "quick" minimizations	tlatorre
	This commit updates the fit_event function to first do a series of "quick" minimizations to try and find a good set of starting parameters for the "real" fit. I also updated the bounds on theta and phi since having hard bounds on these angle coordinates could prevent the fitter from finding the minimum.
2018-08-31	add epsrel argument to likelihood function	tlatorre

2018-08-31	update the lower bound for the energy in the fit	tlatorre

2018-08-31	update the criterion for the fit convergence	tlatorre

2018-08-31	update likelihood check to 1e-5 since that's what we pass to nlopt	tlatorre

2018-08-31	fit in a do while loop until the fit converges to the same likelihood value	tlatorre

2018-08-31	add option to save fit results to a text file	tlatorre

2018-08-31	update printf arguments to keep output aligned	tlatorre

2018-08-31	print out how long the likelihood function takes	tlatorre

2018-08-28	add path to the likelihood fit	tlatorre
	This commit updates the likelihood fit to use the KL path expansion. Currently, I'm just using one coefficient for the path in both x and y.
2018-08-14	fix ev pointer bug	tlatorre

2018-08-14	update pmt hit array in event struct to be MAX_PMTS long	tlatorre

2018-08-14	add lower and upper bounds for the fit parameters	tlatorre

2018-08-14	add function to fit event and clear event after each fit	tlatorre

2018-08-14	set stopping critera	tlatorre

2018-08-14	update the printf format string for the fit parameters	tlatorre

2018-08-14	update fit to use nlopt for minimization	tlatorre
	The GSL library only has the Nelder Mead Simplex algorithm for doing multidimensional minimization without gradient information. The nlopt library has lots of different minimization algorithms so it's easier to switch between them to see which one works best.
2018-08-14	fix how the RMS scattering angle is calculated	tlatorre
	The RMS scattering angle calculation comes from Equation 33.15 in the PDG article on the passage of particles through matter. It's not entirely obvious if this equation is correct for a long track. It seems like it should be integrated along the track to add up the contributions at different energies, but it's not obvious how to do that with the log term. In any case, the way I was previously calculating it (by using the momentum and velocity at each point along the track) was definitely wrong. I will try this out and perhaps try to integrate it later.
2018-08-14	move everything to src directory	tlatorre