#ifndef SCATTERING_H
#define SCATTERING_H

/* Maximum value for the scattering RMS `theta0`.
 *
 * This is to prevent the interpolation from going out of bounds. */
#define MAX_THETA0 1.0

void init_interpolation(void);
double get_probability(double beta, double cos_theta, double theta0);
double get_probability2(double beta);
double get_weighted_quantum_efficiency(void);
void free_interpolation(void);

#endif