Skip to main content

compare_with_convolution

Function compare_with_convolution 

Source
pub fn compare_with_convolution(
    direct: &DirectDose,
    convolution: &Exposure,
    conv_rel_err: f64,
    tolerance_sigmas: f64,
    floor_fraction: f64,
) -> Result<Vec<PointComparison>, ExposureError>
Expand description

The direct and convolution doses side by side at every point of direct.grid; a point is flagged when the two differ by more than tolerance_sigmas combined standard errors (see the module docs). conv_rel_err is the relative uncertainty assigned to the convolution.

Resolution floor. A transport run cannot resolve a dose far below the grid’s peak: with few or no deposits in a cell the replica variance is zero or itself noise. A point where both doses are below floor_fraction times the largest dose of either method on the grid is reported but never flagged (0 judges every point). Choose it with the histories per shot in mind: it should be well above 1 / (histories per replica) of the peak.