avery 0.7.1

A scattering event format for converting between other formats
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
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116

use stripper_xml::{Id, Momentum};

use crate::{
    event::{Particle, Reweight, Scales, Status, WeightInfo},
    Event,
};

impl From<stripper_xml::SubEvent> for Event {
    fn from(ev: stripper_xml::SubEvent) -> Self {
        let weight = WeightInfo {
            weight: Some(ev.weight),
            ..Default::default()
        };
        Self {
            weights: vec![weight],
            scales: Scales {
                mu_r: Some(ev.mu_r),
                mu_f: Some(ev.mu_f),
                mu_ps: None,
            },
            particles: ev.particles.into_iter().map(|p| p.into()).collect(),
            reweights: ev.reweight.into_iter().map(|rw| rw.into()).collect(),
            ..Default::default()
        }
    }
}

impl From<Event> for stripper_xml::SubEvent {
    fn from(ev: Event) -> Self {
        let particles = ev
            .particles
            .into_iter()
            .filter_map(|p| match p.status {
                Some(Status::Incoming | Status::Outgoing) => Some(p.into()),
                _ => None,
            })
            .collect();
        Self {
            weight: ev
                .weights
                .first()
                .map(|w| w.weight.unwrap_or_default())
                .unwrap_or_default(),
            mu_r: ev.scales.mu_r.unwrap_or_default(),
            mu_f: ev.scales.mu_r.unwrap_or_default(),
            particles,
            reweight: ev.reweights.into_iter().map(|rw| rw.into()).collect(),
        }
    }
}

impl From<Particle> for stripper_xml::Particle {
    fn from(particle: Particle) -> Self {
        Self {
            id: Id {
                status: particle.status.unwrap().into(),
                pdg_id: particle.id.unwrap(),
            },
            momentum: Momentum(particle.p.unwrap_or_default()),
        }
    }
}

impl From<stripper_xml::Particle> for Particle {
    fn from(particle: stripper_xml::Particle) -> Self {
        Self {
            id: Some(particle.id.pdg_id),
            p: Some(particle.momentum.0),
            status: Some(particle.id.status.into()),
            ..Default::default()
        }
    }
}

impl From<Reweight> for stripper_xml::Reweight {
    fn from(rw: Reweight) -> Self {
        Self {
            channel: rw.channel,
            reweights: stripper_xml::Reweights {
                x1: rw.x1,
                x2: rw.x2,
                log_coeff: rw.log_coeff,
            },
        }
    }
}

impl From<stripper_xml::Reweight> for Reweight {
    fn from(rw: stripper_xml::Reweight) -> Self {
        Self {
            channel: rw.channel,
            x1: rw.reweights.x1,
            x2: rw.reweights.x2,
            log_coeff: rw.reweights.log_coeff,
        }
    }
}

impl From<Status> for stripper_xml::Status {
    fn from(status: Status) -> Self {
        match status {
            Status::Incoming => stripper_xml::Status::Incoming,
            Status::Outgoing => stripper_xml::Status::Outgoing,
            _ => panic!("Can only convert incoming and outgoing particles to STRIPPER-XML format")
        }
    }
}

impl From<stripper_xml::Status> for Status {
    fn from(status: stripper_xml::Status) -> Self {
        match status {
            stripper_xml::Status::Incoming => Status::Incoming,
            stripper_xml::Status::Outgoing => Status::Outgoing,
        }
    }
}