laddu 0.19.2

Amplitude analysis made short and sweet
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

use std::collections::HashMap;

use laddu::generation::{
    CompositeGenerator, EventGenerator, GeneratedParticle, GeneratedReaction, InitialGenerator,
    MandelstamTDistribution, Reconstruction, StableGenerator,
};

#[test]
fn umbrella_crate_exposes_generation_api() {
    let beam = GeneratedParticle::initial(
        "beam",
        InitialGenerator::beam_with_fixed_energy(0.0, 8.0),
        Reconstruction::Stored,
    );
    let target = GeneratedParticle::initial(
        "target",
        InitialGenerator::target(0.938272),
        Reconstruction::Missing,
    );
    let kshort1 = GeneratedParticle::stable(
        "kshort1",
        StableGenerator::new(0.497611),
        Reconstruction::Stored,
    );
    let kshort2 = GeneratedParticle::stable(
        "kshort2",
        StableGenerator::new(0.497611),
        Reconstruction::Stored,
    );
    let kk = GeneratedParticle::composite(
        "kk",
        CompositeGenerator::new(1.1, 1.6),
        (&kshort1, &kshort2),
        Reconstruction::Composite,
    );
    let recoil = GeneratedParticle::stable(
        "recoil",
        StableGenerator::new(0.938272),
        Reconstruction::Stored,
    );
    let reaction = GeneratedReaction::two_to_two(
        beam,
        target,
        kk,
        recoil,
        MandelstamTDistribution::Exponential { slope: 0.1 },
    )
    .unwrap();
    let reconstructed = reaction.reconstructed_reaction().unwrap();
    let generator = EventGenerator::new(reaction, HashMap::new(), Some(12345));
    let dataset = generator.generate_dataset(4).unwrap();

    assert_eq!(dataset.n_events(), 4);
    assert_eq!(
        dataset.metadata().p4_names(),
        &["beam", "target", "kk", "kshort1", "kshort2", "recoil"]
    );
    assert_eq!(
        reconstructed.decay("kk").unwrap().daughters(),
        ["kshort1", "kshort2"]
    );
}