#ifndef MUON_H
#define MUON_H

#include <stddef.h> /* for size_t */

#define EULER_CONSTANT 0.57721

/* Number of photons in the range 200 nm - 800 nm generated per MeV of energy
 * lost to radiation for electrons.
 *
 * FIXME: This is just a rough estimate, should use an energy dependent
 * quantity from simulation.
 *
 * FIXME: Actually determine what this is. */
#define MUON_PHOTONS_PER_MEV 7368.0

void muon_get_position_distribution_parameters(double T0, double *a, double *b);
double muon_get_angular_distribution_alpha(double T0);
double muon_get_angular_distribution_beta(double T0);
void muon_get_delta_ray_distribution_parameters(double T0, double *a, double *b);
double muon_get_delta_ray_photons(double T0);
double muon_get_range(double T, double rho);
double muon_get_dEdx_rad(double T, double rho);
double muon_get_dEdx(double T, double rho);

#endif