Age | Commit message (Collapse) | Author | |
---|---|---|---|
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 | 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 | move everything to src directory | tlatorre | |