From 8e6895f1e55290fb5ce6b7a6ff77b3819806c51f Mon Sep 17 00:00:00 2001 From: tlatorre Date: Tue, 14 Aug 2018 11:29:22 -0500 Subject: fix how the RMS scattering angle is calculated 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. --- src/fit.c | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) (limited to 'src/fit.c') diff --git a/src/fit.c b/src/fit.c index b6e441a..1f90f12 100644 --- a/src/fit.c +++ b/src/fit.c @@ -134,7 +134,7 @@ int main(int argc, char **argv) /* New event, so we need to fit the old event. */ iter = 0; - gsl_vector_set(x,0,1000.0); + gsl_vector_set(x,0,500.0); gsl_vector_set(x,1,100.0); gsl_vector_set(x,2,100.0); gsl_vector_set(x,3,100.0); -- cgit