Struct IkedaCarpenter 

pub struct IkedaCarpenter { /* private fields */ }

Analytical Ikeda–Carpenter resolution model.

Synthesizes a TabulatedResolution at construction and applies it through the same broadening path as a Monte-Carlo file. Cloning is cheap (the synthesized table is the only large field). Re-synthesize (construct anew) when fitting changes the parameters.

Implementations

impl IkedaCarpenter

pub fn new( params: IkedaCarpenterParams, flight_path_m: f64, grid: &SynthesisGrid, ) -> Result<Self, ResolutionParseError>

Build the model, synthesizing its kernel table over grid.

Errors

Returns ResolutionParseError::InvalidFormat for a non-positive flight path, a degenerate grid (n_energies < 2, n_tau < 8, e_min ≤ 0, e_max ≤ e_min), a non-positive β(E), a parameter/grid combination whose τ-grid cannot resolve the prompt core and requested folds within the MAX_TAU_SAMPLES cap at some reference energy (see tau_geometry — remedy: larger β, R = 0, or a wider/disabled fold), or if the synthesized kernels fail TabulatedResolution::from_kernels validation.

pub fn tabulated(&self) -> &TabulatedResolution

The synthesized tabulated kernel set (the broadening engine).

pub fn params(&self) -> &IkedaCarpenterParams

The IC parameters.

pub fn flight_path_m(&self) -> f64

Flight-path length (m).

pub fn with_flight_path( &self, flight_path_m: f64, ) -> Result<Self, ResolutionParseError>

The same pulse read against a different flight path.

The synthesized kernels are emission-time distributions built from α(E), β(E) and R(E) — no flight path appears in the synthesis. The flight path enters only the TOF↔energy map and the nominal arrival time, so rebinding it is exact and does not resynthesize the table.

Errors

Returns ResolutionParseError::InvalidFormat if flight_path_m is not positive and finite.

pub fn ref_energies(&self) -> &[f64]

Reference energies (eV, ascending) the table was synthesized on.

pub fn kernel_at( &self, energy_ev: f64, ) -> Result<(Vec<f64>, Vec<f64>), ResolutionParseError>

Evaluate the (burst/channel-folded) IC kernel at one energy.

Returns ascending TOF-offsets (µs, 0 = pulse start) and peak-normalized weights (max = 1), matching the TabulatedResolution storage convention.

Errors

Returns ResolutionParseError::InvalidFormat when the requested energy or an energy-dependent rate/fraction is non-physical, or when the τ-grid cannot resolve the prompt core and requested folds within MAX_TAU_SAMPLES at this energy. Construction validates every reference energy, but a probe outside [e_min, e_max] can still leave the physical or resolvable region.

pub fn delays_us( &self, energy_ev: f64, ) -> Result<(f64, f64), ResolutionParseError>

The first and last delay, in µs after the nominal arrival, of a neutron of energy_ev, outside which Self::detector_bin_probabilities gives it less than NEGLIGIBLE_ARRIVAL_PROBABILITY chance of arriving: without a fold, 0 and the delay where 1 − ic_cdf falls to that chance; with a fold, the fold’s reach before 0 and a bound on the end of the sampled pulse those probabilities integrate.

Errors

ResolutionParseError::InvalidFormat when energy_ev is not positive and finite, or a law is singular or out of range there.

pub fn rise_us(&self, energy_ev: f64) -> Result<f64, ResolutionParseError>

The time in µs from the first sample of the pulse at energy_ev to its peak.

Errors

As Self::source_pulse_at.

pub fn source_pulse_at( &self, energy_ev: f64, ) -> Result<(Vec<f64>, Vec<f64>), ResolutionParseError>

Evaluate the physical source pulse at one true neutron energy.

Returns sampled moderator-delay coordinates in µs and peak-normalized densities. Unlike Self::kernel_at, this method does not move the pulse mode to zero. With no symmetric proton/channel fold, the delay is causal and starts at zero. A symmetric fold may extend the sampled support below zero relative to its stated time origin.

Errors

Returns ResolutionParseError::InvalidFormat if energy_ev is not a positive finite true energy, if an energy law is unphysical at that energy, or if the requested pulse cannot be resolved by the configured sampling limits.

pub fn detector_bin_probabilities( &self, true_energy_ev: f64, detector_time_edges_us: &[f64], timing_offset_us: f64, ) -> Result<Vec<f64>, ResolutionParseError>

Probability that a neutron of known true energy is recorded in each supplied detector-time bin.

detector_time_edges_us are the actual measured bin edges. The nominal arrival time is timing_offset_us + TOF_FACTOR * flight_path_m / sqrt(true_energy_ev). timing_offset_us represents the shared clock/detector offset; it does not absorb or remove the moderator pulse’s physical mode.

The returned vector has one entry per adjacent edge pair and is not renormalized to the supplied window: bins that do not cover the full pulse correctly sum to less than one.

With no burst or channel fold, probabilities come directly from the analytical IC CDF. With either optional fold, the continuous pulse is represented on the configured synthesis grid and integrated as a piecewise-linear density. The finite numerical support is not silently renormalized; omitted physical tail probability remains omitted.

Errors

Returns ResolutionParseError::InvalidFormat unless the true energy is physical, the timing offset is finite, and at least two finite bin edges are supplied in strictly increasing order.

Trait Implementations

impl Clone for IkedaCarpenter

fn clone(&self) -> IkedaCarpenter

Returns a duplicate of the value.

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

Performs copy-assignment from source.

impl Debug for IkedaCarpenter

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

Formats the value using the given formatter.