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

use super::*;
use uom::si::angle::radian;
use uom::si::frequency::megahertz;
use uom::si::length::millimeter;

#[test]
fn anode_wires_radius() {
    let r = Length::new::<millimeter>(182.0);
    assert_eq!(ANODE_WIRES_RADIUS, r);
}

#[test]
fn trigger_clock_frequency() {
    let f = Frequency::new::<megahertz>(62.5);
    assert_eq!(TRG_CLOCK_FREQ, f);
}

#[test]
fn spacepoint_x_coordinate() {
    let p = SpacePoint {
        r: Length::new::<millimeter>(5.0),
        phi: Angle::new::<radian>(4.0f64.atan2(3.0)),
        z: Length::new::<millimeter>(0.0),
    };

    let diff = (p.x() - Length::new::<millimeter>(3.0)).abs();
    assert!(diff < Length::new::<millimeter>(1e-6));
}

#[test]
fn spacepoint_y_coordinate() {
    let p = SpacePoint {
        r: Length::new::<millimeter>(5.0),
        phi: Angle::new::<radian>(4.0f64.atan2(3.0)),
        z: Length::new::<millimeter>(0.0),
    };

    let diff = (p.y() - Length::new::<millimeter>(4.0)).abs();
    assert!(diff < Length::new::<millimeter>(1e-6));
}

#[test]
fn spacepoint_distance() {
    let p1 = SpacePoint {
        r: Length::new::<millimeter>(10.0),
        phi: Angle::new::<radian>(1.5),
        z: Length::new::<millimeter>(0.5),
    };
    let p2 = SpacePoint {
        r: Length::new::<millimeter>(5.0),
        phi: Angle::new::<radian>(0.5),
        z: Length::new::<millimeter>(-1.0),
    };

    assert_eq!(p1.distance(p2), p2.distance(p1));

    let diff = p1.distance(p2) - Length::new::<millimeter>(8.55685511232);
    assert!(diff.abs() < Length::new::<millimeter>(1e-6));

    let p3 = SpacePoint {
        r: Length::new::<millimeter>(15.0),
        phi: Angle::new::<radian>(1.5),
        z: Length::new::<millimeter>(0.5),
    };

    let diff = p1.distance(p3) - Length::new::<millimeter>(5.0);
    assert!(diff.abs() < Length::new::<millimeter>(1e-6));
}