#ifndef MUON_H
#define MUON_H

#define EULER_CONSTANT 0.57721

double get_range(double T, double rho);
double get_T(double T0, double x, double rho);
double get_dEdx(double T, double rho);
double get_expected_charge(double x, double T, double *pos, double *dir, double *pmt_pos, double *pmt_normal, double r);

#endif