#ifndef ELECTRON_H
#define ELECTRON_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. */
#define ELECTRON_PHOTONS_PER_MEV 400.0

void electron_get_position_distribution_parameters(double T0, double *a, double *b);
double electron_get_angular_distribution_alpha(double T0);
double electron_get_angular_distribution_beta(double T0);
double electron_get_angular_pdf_delta_ray(double cos_theta, double alpha, double beta, double mu);
double electron_get_angular_pdf(double cos_theta, double alpha, double beta, double mu);
double electron_get_range(double T, double rho);
double electron_get_dEdx_rad(double T, double rho);
double electron_get_dEdx(double T, double rho);

#endif