aboutsummaryrefslogtreecommitdiff
path: root/src/random.h
diff options
context:
space:
mode:
authortlatorre <tlatorre@uchicago.edu>2018-08-27 10:12:17 -0500
committertlatorre <tlatorre@uchicago.edu>2018-08-27 10:12:17 -0500
commit779266ec72a5c76ee52043ab3ae17479ba6a9788 (patch)
tree9a03c335731aeff2cadb728bc3fe4309aafef874 /src/random.h
parent4e3eb54e364f578b3ea6608b0b3fc34125985ab8 (diff)
downloadsddm-779266ec72a5c76ee52043ab3ae17479ba6a9788.tar.gz
sddm-779266ec72a5c76ee52043ab3ae17479ba6a9788.tar.bz2
sddm-779266ec72a5c76ee52043ab3ae17479ba6a9788.zip
add code to expand the track of a particle using a KL expansion
To fit the path of muons and electrons I use the Karhunen-Loeve expansion of a random 2D walk in the polar angle in x and y. This allows you to decompose the path into a sum over sine functions whose coefficients become random variables. The nice thing about fitting the path in this way is that you can capture *most* of the variation in the path using a small number of variables by only summing over the first N terms in the expansion and it is easy to calculate the probability of the coefficients since they are all uncorrelated.
Diffstat (limited to 'src/random.h')
-rw-r--r--src/random.h1
1 files changed, 1 insertions, 0 deletions
diff --git a/src/random.h b/src/random.h
index 309cd71..733279c 100644
--- a/src/random.h
+++ b/src/random.h
@@ -4,5 +4,6 @@
#include "mt19937ar.h"
double randn(void);
+void rand_sphere(double *dir);
#endif