@@ -24,12 +24,12 @@ cs = jet_reconstruct(particles::Vector{T}; p = -1, R = 1.0, recombine = +, strat
24
24
- Any type that supplies the methods ` pt2() ` , ` phi() ` , ` rapidity() ` , ` px() ` , ` py() ` , ` pz() ` , ` energy() ` can be used
25
25
- These methods have to be defined in the namespace of this package, i.e., ` JetReconstruction.pt2(::T) `
26
26
- The ` PseudoJet ` type from this package, or a 4-vector from ` LorentzVectorHEP ` are suitable (and have the appropriate definitions)
27
- - ` p ` - the transverse momentum power used in the $d_ {ij}$ metric for deciding on closest jets, as $k^{2p}_ \text{T}$
28
- - ` -1 ` gives anti-${k}_ \text{T}$ clustering (default)
27
+ - ` p ` - the transverse momentum power used in the $d_ {ij}$ metric for deciding on closest jets, as $k^{2p}_ \text{T}$. Different values of $p$ then give different reconstruction algorithms:
28
+ - ` -1 ` gives anti-$` {k}_\text{T} ` $ clustering (default)
29
29
- ` 0 ` gives Cambridge/Achen
30
30
- ` 1 ` gives inclusive $k_ \text{T}$
31
31
- ` R ` - the cone size parameter; no particles more geometrically distance than ` R ` will be merged (default 1.0)
32
- - ` recombine ` - the function used to merge two pseudojets (default is a simple 4-vector addition of $(E, \mathbf{p})$)
32
+ - ` recombine ` - the function used to merge two pseudojets (default is a simple 4-vector addition of $` (E, \mathbf{p}) ` $)
33
33
- ` strategy ` - the algorithm strategy to adopt, as described below (default ` JetRecoStrategy.Best ` )
34
34
35
35
The object returned is a ` ClusterSequence ` , which internally tracks all merge steps.
@@ -65,7 +65,7 @@ Generally one can use the `jet_reconstruct` interface, shown above, as the *Best
65
65
Another option, if one wishes to use a specific strategy, is to call that strategy's interface directly, e.g.,
66
66
67
67
``` julia
68
- # N2Plain
68
+ # For N2Plain strategy called directly
69
69
plain_jet_reconstruct (particles:: Vector{T} ; p = - 1 , R = 1.0 , recombine = + )
70
70
```
71
71
@@ -75,7 +75,7 @@ Note that there is no `strategy` option in these interfaces.
75
75
76
76
See the ` examples/jetreco.jl ` script for a full example of how to call jet reconstruction.
77
77
78
- ``` julia
78
+ ``` sh
79
79
julia --project=. examples/jetreco.jl --maxevents=100 --nsamples=1 --strategy=N2Plain test/data/events.hepmc3
80
80
...
81
81
julia --project=. examples/jetreco.jl --maxevents=100 --nsamples=1 --strategy=N2Tiled test/data/events.hepmc3
0 commit comments