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

Age	Commit message (Collapse)	Author
2018-08-14	add refractive index for heavy and light water from snoman	tlatorre

2018-08-14	update muon.c to use the new table for heavy water	tlatorre

2018-08-14	add a table for muon dE/dx and range for D2O	tlatorre
	There was no entry for heavy water at http://pdg.lbl.gov/2018/AtomicNuclearProperties/index.html, so I emailed Don Groom who maintains the website. Amazingly he agreed to rerun the code to add an entry for D2O. Apparently all of the old Fortran code was not set up to deal with isotopes, but he updated everything and reran the code for heavy water! The new results are at http://pdgprod.lbl.gov/~deg/AtomicNuclearProperties/ but they should make it to the main page soon.
2018-08-14	move everything to src directory	tlatorre

2018-08-14	initial commit of likelihood fit for muons	tlatorre
	This commit contains code to fit for the energy, position, and direction of muons in the SNO detector. Currently, we read events from SNOMAN zebra files and fill an event struct containing the PMT hits and fit it with the Nelder Mead simplex algorithm from GSL. I've also added code to read in ZEBRA title bank files to read in the DQXX files for a specific run. Any problems with channels in the DQCH and DQCR banks are flagged in the event struct by masking in a bit in the flags variable and these PMT hits are not included in the likelihood calculation. The likelihood for an event is calculated by integrating along the particle track for each PMT and computing the expected number of PE. The charge likelihood is then calculated by looping over all possible number of PE and computing: P(q\|n)*P(n\|mu) where q is the calibrated QHS charge, n is the number of PE, and mu is the expected number of photoelectrons. The latter is calculated assuming the distribution of PE at a given PMT follows a Poisson distribution (which I think should be correct given the track, but is probably not perfect for tracks which scatter a lot). The time part of the likelihood is calculated by integrating over the track for each PMT and calculating the average time at which the PMT is hit. We then assume the PDF for the photons to arrive is approximately a delta function and compute the first order statistic for a given time to compute the probability that the first photon arrived at a given time. So far I've only tested this with single tracks but the method was designed to be easy to use when you are fitting for multiple particles.
2018-07-05	analyze_genie_mc.py -> analyze-genie-mc	tlatorre

2018-07-05	refractive_index.c -> optics.c	tlatorre

2018-07-05	add a function to compute the energy along a muon track	tlatorre

2018-07-05	add a simple test for the various muon functions	tlatorre

2018-07-05	add some comments and divide by the density	tlatorre

2018-07-05	add functions to compute the range and dE/dx for muons in water	tlatorre

2018-07-05	add a text file with values for dE/dx and ranges for muons in water	tlatorre

2018-07-04	only require a 1% agreement, because I don't know how Schiebener converted ↵	tlatorre
	temperature and pressure to density in his tables
2018-07-04	add a function to compute the refractive index of water as a function of ↵	tlatorre
	wavelength
2018-07-04	update	tlatorre

2018-07-04	initial commit of a file to calculate the quantum efficiency	tlatorre

2018-07-04	add a file with the measured quantum efficiency of the sno PMTs from snoman	tlatorre

2018-07-04	update test	tlatorre

2018-07-04	update solid angle test and fix a bug in the solid angle calculation	tlatorre

2018-07-04	update test	tlatorre

2018-07-04	update	tlatorre

2018-07-04	fix a bug in get_solid_angle()	tlatorre

2018-07-04	add a test of the get_solid_angle function	tlatorre

2018-07-04	update TODO	tlatorre

2018-07-04	initial commit of a function to calculate the solid angle subtended by a ↵	tlatorre
	circular disk
2018-07-04	update index script to fit the inverse of the index of refraction to a ↵	tlatorre
	polynomial
2018-07-04	initial commit of a script to calculate the index of refraction of water	tlatorre

2018-05-30	add an energy resolution parameter and add a signal to the visible energy plot	tlatorre

2018-05-30	plot kinetic energy instead of calresp0 variable. Also update the expected ↵	tlatorre
	number of events to 230 per year since that's what's in Richie's thesis
2018-05-30	forgot to check for nc events	tlatorre

2018-05-23	update TODO	tlatorre

2018-05-23	update script to analyze GENIE Monte Carlo	tlatorre
	- plot energy for just the neutron follower cut - update pions to produce 2 muon like rings - add a check to see if the pion is above the cerenkov threshold commit 924d4277be5ff67938247cd1babc83eaf85e7852 Author: tlatorre <tlatorre@uchicago.edu> Date: Wed May 23 12:15:32 2018 -0400 update
2018-05-23	fix small bug in convert-genie-to-gst script	tlatorre

2018-05-22	add TODO text file	tlatorre

2018-05-22	initial commit of script to plot GENIE atmospheric Monte Carlo information	tlatorre

2018-05-15	add a python script to calculate the event rate for SNO and dark matter ↵	tlatorre
	experiments
2018-05-01	add some phase space factors for the integral over the earth	tlatorre

2018-05-01	update get_event_rate to return an event rate assuming some simplified model ↵	tlatorre
	for the dark matter in the galaxy
2018-05-01	update makefile to compile with -O4	tlatorre

2018-05-01	make epsabs and epsrel global variables	tlatorre

2018-05-01	add the integration over the gamma production point	tlatorre

2018-05-01	start working on the code to numerically integrate to find the event rate	tlatorre

2018-05-01	add some vector functions and add a function to rotate from earth centered ↵	tlatorre
	coordinates into sno coordinates
2018-05-01	initial commit	tlatorre