9 unstable releases (3 breaking)
0.4.1 | Sep 20, 2023 |
---|---|
0.4.0 | Sep 20, 2023 |
0.3.2 | Jul 28, 2023 |
0.2.1 | Jun 19, 2023 |
0.1.0 | Feb 12, 2021 |
#383 in Science
Used in cres
56KB
1.5K
SLoC
jetty
Implementations of common inclusive jet clustering algorithms.
In the current version, the following distance measures are implemented:
For state-of-the-art implementations of many more jet algorithms, have a look at the excellent fastjet library.
Examples
Cluster a number of partons into jets using the anti-kt algorithm with radius 0.4:
use jetty::{anti_kt_f, pseudojet_f, Cluster, ClusterHistory, ClusterStep};
let partons = vec![
pseudojet_f(0.2626773221934335, -0.08809521946454194, -0.1141608706693822, -0.2195584284654444),
pseudojet_f(2.21902459329915, -0.7529973704809976, -0.9658189214109036, -1.850475321845671)
];
// get all jets
let all_jets = partons.clone().cluster(anti_kt_f(0.4));
assert_eq!(all_jets.len(), 1);
// get all jets with at least 40 GeV
let jets_40gev = partons.clone().cluster_if(
anti_kt_f(0.4),
|jet| jet.pt2() > 40. * 40.
);
assert_eq!(jets_40gev.len(), 0);
// go through the cluster history step-by-step
let history = ClusterHistory::new(partons, anti_kt_f(0.4));
for step in history {
match step {
ClusterStep::Jet(j) => println!("Found a jet: {j:?}"),
ClusterStep::Combine([_j1, _j2]) => println!("Combined two pseudojets"),
}
}
License: GPL-3.0-or-later
Dependencies
~1.5MB
~29K SLoC