#ifndef TRACK_HH
#define TRACK_HH

#include <iostream>
#include <iomanip>
#include <math.h>

// generate uniformy distributed random number between min and max
double unirandom(float min, float max);

// state of a particle
struct state {
  float x;
  float y;
  float theta;

// extrapolate a state to a new x position (=detector plane)
state extrapolate(state st, float x);

// scatter particle with probablity 50% by +/- 2 deg
void scatter(state *st);