alpha_g_physics 0.1.5

Data analysis library for the ALPHA-g experiment
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
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101

use super::*;
use alpha_g_detector::alpha16::aw_map::TPC_ANODE_WIRES;

#[test]
fn try_wire_gain_wire_gain_map_error() {
    for run_number in 0..=9276 {
        for i in 0..TPC_ANODE_WIRES {
            let wire = TpcWirePosition::try_from(i).unwrap();
            assert!(try_wire_gain(run_number, wire).is_err());
        }
    }
}

#[test]
fn try_wire_gain_correctness_9277() {
    for i in 0..TPC_ANODE_WIRES {
        let wire = TpcWirePosition::try_from(i).unwrap();
        let gain = try_wire_gain(9277, wire);
        assert!(gain.is_ok());
    }

    assert_eq!(
        try_wire_gain(9277, TpcWirePosition::try_from(0).unwrap()).unwrap(),
        1.1678439324476453
    );
    assert_eq!(
        try_wire_gain(9277, TpcWirePosition::try_from(32).unwrap()).unwrap(),
        1.0491624717198518
    );
    assert_eq!(
        try_wire_gain(9277, TpcWirePosition::try_from(64).unwrap()).unwrap(),
        1.0873347124407264
    );
    assert_eq!(
        try_wire_gain(9277, TpcWirePosition::try_from(96).unwrap()).unwrap(),
        1.290917387251699
    );
    assert_eq!(
        try_wire_gain(9277, TpcWirePosition::try_from(128).unwrap()).unwrap(),
        1.3295441389947555
    );
    assert_eq!(
        try_wire_gain(9277, TpcWirePosition::try_from(160).unwrap()).unwrap(),
        1.2381314918699662
    );
    assert_eq!(
        try_wire_gain(9277, TpcWirePosition::try_from(192).unwrap()).unwrap(),
        1.0467826211111837
    );
    assert_eq!(
        try_wire_gain(9277, TpcWirePosition::try_from(224).unwrap()).unwrap(),
        1.1702237830563136
    );
}

#[test]
fn try_wire_gain_ok_11084() {
    for i in 0..TPC_ANODE_WIRES {
        let wire = TpcWirePosition::try_from(i).unwrap();
        let gain = try_wire_gain(11084, wire);
        if i == 111 {
            assert!(gain.is_err());
        } else {
            assert!(gain.is_ok());
        }
    }
}

#[test]
fn try_wire_gain_correctness_11084() {
    assert_eq!(
        try_wire_gain(11084, TpcWirePosition::try_from(224).unwrap()).unwrap(),
        1.0729483132909592
    );
}

#[test]
fn try_wire_gain_ok_11356() {
    for i in 0..TPC_ANODE_WIRES {
        let wire = TpcWirePosition::try_from(i).unwrap();
        let gain = try_wire_gain(11356, wire);
        assert!(gain.is_ok());
    }
}

#[test]
fn try_wire_gain_correctness_11356() {
    assert_eq!(
        try_wire_gain(11356, TpcWirePosition::try_from(111).unwrap()).unwrap(),
        1.0839053881577092
    );
}

#[test]
fn try_wire_gain_correctness_sim() {
    for i in 0..TPC_ANODE_WIRES {
        let wire = TpcWirePosition::try_from(i).unwrap();
        let gain = try_wire_gain(u32::MAX, wire).unwrap();
        assert_eq!(gain, 1.0);
    }
}