laddu 0.20.0

Amplitude analysis tools for Rust
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82

//! End-to-end statistical closure test for generation and fitting.

#[allow(dead_code)]
#[path = "../examples/common/mod.rs"]
mod common;

use std::f64::consts::PI;

use common::{
    closure::{ClosureConfig, run_closure, wrapped_phase_residual},
    ksks::{F2_MAGNITUDE_TRUTH, F2_PHASE_TRUTH},
};
use laddu::prelude::Execution;

#[test]
fn generated_two_wave_sample_recovers_injected_coupling() {
    let result = run_closure(
        ClosureConfig {
            data_events: 4_000,
            normalization_events: 40_000,
            pilot_proposals: 20_000,
            ..ClosureConfig::default()
        },
        Execution::default(),
    )
    .expect("generation-to-fit closure should succeed");

    assert_eq!(result.data_report.produced, 4_000);
    assert_eq!(result.normalization_report.produced, 40_000);
    assert!(result.initial_nll.is_finite());
    assert!(
        result
            .initial_gradient
            .iter()
            .all(|value| value.is_finite())
    );
    assert!(result.fit.fx.is_finite());
    assert!(result.fit.fx < result.initial_nll);

    let magnitude = result
        .fitted("f2_magnitude")
        .expect("fit should return f2_magnitude");
    let phase = result
        .fitted("f2_phase")
        .expect("fit should return f2_phase");
    assert!((0.0..=2.0).contains(&magnitude));
    assert!((-PI..PI).contains(&phase));
    assert!((magnitude - F2_MAGNITUDE_TRUTH).abs() < 0.15);
    assert!(wrapped_phase_residual(phase, F2_PHASE_TRUTH).abs() < 0.35);

    assert_eq!(result.projection.data.histogram.bin_edges().len(), 51);
    assert_eq!(result.projection.projections.len(), 3);
    let data_yield = result
        .projection
        .data
        .histogram
        .counts()
        .iter()
        .sum::<f64>();
    let fit_yield = result.projection.projections[0]
        .histogram
        .counts()
        .iter()
        .sum::<f64>();
    assert!((fit_yield - data_yield).abs() < 1e-9 * data_yield);
    assert!(
        result.projection.projections[1]
            .histogram
            .counts()
            .iter()
            .sum::<f64>()
            > 0.0
    );
    assert!(
        result.projection.projections[2]
            .histogram
            .counts()
            .iter()
            .sum::<f64>()
            > 0.0
    );
}