#ifndef PROTON_H
#define PROTON_H

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

/* 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 PROTON_PHOTONS_PER_MEV 400.0

double proton_get_range(double T, double rho);
double proton_get_dEdx_rad(double T, double rho);
double proton_get_dEdx(double T, double rho);

#endif