| Age | Commit message (Collapse) | Author | |
|---|---|---|---|
| 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 | use the index of refraction from snoman when computing the angular PDF for ↵ | tlatorre | |
| photons | |||
| 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 | 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 | |
 |
index : sddm | |
| Code for the self-destructing dark matter search in SNO | Anthony LaTorre |
| aboutsummaryrefslogtreecommitdiff |