power_law_distribution

romanisim.cr.power_law_distribution(x, slope=-4.33)

Return unnormalized power-law distribution parameterized by a log-log slope, used to describe the cosmic ray path lengths.

Parameters:

x1-d array of floats

An array of cosmic ray path lengths (units: micron).

slopefloat

The log-log slope of the distribution, default based on Miles et al. (2021).

Returns:

power_law1-d array of floats

Power-law distribution (pdf) evaluated on x grid of points.