Struct MomentumTransferSampler 

pub struct MomentumTransferSampler { /* private fields */ }

Samples the momentum transfer q of an inelastic collision with incident kinetic energy E and loss W (module docs). Built once per (E, W).

Implementations

impl MomentumTransferSampler

pub fn new( penn: &SinglePolePenn, energy_ev: f64, loss_ev: f64, ) -> Result<Self, InelasticTableError>

Build the sampler with the default tolerance DEFAULT_Q_TOLERANCE.

pub fn with_tolerance( penn: &SinglePolePenn, energy_ev: f64, loss_ev: f64, tol: f64, ) -> Result<Self, InelasticTableError>

Build the sampler with density-model tolerance tol in [1e-9, 1e-2]. Fails if 0 < W <= E does not hold or the slice has no support (W below the ELF table, or no q allowed).

pub fn sample(&self, u: f64) -> f64

The momentum transfer (as a wavenumber, m⁻¹) at cumulative probability u in [0, 1] (clamped).

pub fn q_range_per_m(&self) -> (f64, f64)

The support (q_min, q_max) of the sampler, m⁻¹: the kinematic limits q± cut to the tabulated ELF range.

pub fn kinematic_limits_per_m(&self) -> (f64, f64)

The kinematic limits q- = k - k' and q+ = k + k', m⁻¹.

pub fn cos_theta(&self, q_per_m: f64) -> f64

The cosine of the primary’s deflection angle for momentum transfer q (m⁻¹), (k² + k'² - q²)/(2 k k'), clamped to [-1, 1].

pub fn diimfp_per_m_ev(&self) -> f64

The DIIMFP p(E, W) (m⁻¹ eV⁻¹) of the sampler’s own density model: the integral of the density over ln q, divided by π T'. It agrees with SinglePolePenn::diimfp_per_m_ev to the model tolerance.

Trait Implementations

impl Clone for MomentumTransferSampler

fn clone(&self) -> MomentumTransferSampler

Returns a duplicate of the value.

fn clone_from(&mut self, source: &Self)

Performs copy-assignment from source.

impl Debug for MomentumTransferSampler

fn fmt(&self, f: &mut Formatter<'_>) -> Result

Formats the value using the given formatter.

impl PartialEq for MomentumTransferSampler

fn eq(&self, other: &MomentumTransferSampler) -> bool

Equality operator ==.

fn ne(&self, other: &Rhs) -> bool

Inequality operator !=.

impl StructuralPartialEq for MomentumTransferSampler