cres 0.9.1

Cell resampling for collider events
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

// minimal example for cell resampling
// run with `cargo run --release --example minimal -- IN.hepmc OUT.hepmc`
// set the environment variable `RUST_LOG=info` for command-line output
use std::error::Error;

use cres::prelude::*;

use env_logger;

fn main() -> Result<(), Box<dyn Error>> {
    // initialise logging from the RUST_LOG environment variable
    env_logger::init();

    // access command line arguments, ignoring the program name
    let mut args = std::env::args().skip(1);
    let infile = args.next().unwrap().into();
    let outfile = args.next().unwrap().into();

    // How to access (read & write) events
    let event_io = IOBuilder::default();
    let event_io = event_io.build_from_files(infile, outfile)?;

    // How to convert into internal event format
    let converter = Converter::new();

    // Cluster outgoing particles into IRC safe objects
    // To actually perform clustering use `DefaultClustering` instead
    let clustering = NO_CLUSTERING;

    // Resample with default settings
    let resampler = ResamplerBuilder::default().build();

    let mut cres = CresBuilder {
        event_io,
        converter,
        clustering,
        resampler,
        unweighter: NO_UNWEIGHTING, // disable unweighting
    }
    .build();

    // Run the resampler
    cres.run()?;
    Ok(())
}