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.

1 files changed, 4 insertions, 2 deletions

diff --git a/src/Makefile b/src/Makefile
index ae8807b..1ea1685 100644
--- a/src/Makefile
+++ b/src/Makefile
@@ -1,7 +1,7 @@
 CFLAGS=-O0 -Wall -g -DSWAP_BYTES
 LDLIBS=-lm -lgsl -lgslcblas -lnlopt_cxx -lstdc++
 
-all: test test-vector test-likelihood fit test-charge
+all: test test-vector test-likelihood fit test-charge test-path
 
 Makefile.dep:
 	-$(CC) -MM *.c > Makefile.dep
@@ -12,12 +12,14 @@ calculate_limits: calculate_limits.c
 
 solid_angle.o: solid_angle.c
 
-test: test.o solid_angle.o optics.o muon.o vector.o quantum_efficiency.o pdg.o scattering.o misc.o mt19937ar.o sno_charge.o
+test: test.o solid_angle.o optics.o muon.o vector.o quantum_efficiency.o pdg.o scattering.o misc.o mt19937ar.o sno_charge.o path.o random.o
 
 test-vector: test-vector.o vector.o mt19937ar.o
 
 test-likelihood: test-likelihood.o muon.o random.o optics.o quantum_efficiency.o mt19937ar.o pdg.o vector.o solid_angle.o scattering.o
 
+test-path: test-path.o mt19937ar.o vector.o path.o random.o
+
 test-charge: test-charge.o sno_charge.o misc.o
 
 test-zebra: test-zebra.o zebra.o pack2b.o