Age | Commit message (Collapse) | Author | |
---|---|---|---|
2018-09-26 | increase number of points in cos(theta) interpolation to 1000 | tlatorre | |
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-08-27 | fix how multiple Coulomb scattering is treated | tlatorre | |
Previously I had been assuming that a particle undergoing many small angle Coulomb scatters had a track direction whose polar angle was a Gaussian. However, this was just due to a misunderstanding of the PDG section "Multiple scattering through small angles" in the "Passage of particles through matter" article. In fact, what is described by a Gaussian is the polar angle projected onto a plane. Therefore the distribution of the polar angle is actually: (1/(sqrt(2*pi)*theta0**2))*theta*exp(-theta**2/(2*theta0)) This commit updates the code in scattering.c to correctly calculate the probability that a photon is emitted at a particular angle. I also updated test-likelihood.c to simulate a track correctly. | |||
2018-08-14 | use the index of refraction from snoman when computing the angular PDF for ↵ | tlatorre | |
photons | |||
2018-08-14 | move everything to src directory | tlatorre | |