#########
## Read an event with 354 particles
#########
## Example of basic event shape analysis
#########
Jet parameters: R_jet=1, pTcut=200
Trimming parameters: R_sub=0.2, fcut=0.05
#####
N_jet=1.99971
N_jet_trim (using built-in)=1.97156
N_jet_trim (using trimmer)=2
#####
H_T=1893.67
H_T_trim (using built-in)=1867.45
H_T_trim (using trimmer)=1867.45
#####
Missing H_T=76.1162
Missing H_T_trim (using built-in)=91.8116
Missing H_T_trim (using trimmer)=91.8116
#####
#--------------------------------------------------------------------------
#                         FastJet release 3.0.3
#                 M. Cacciari, G.P. Salam and G. Soyez                  
#     A software package for jet finding and analysis at colliders      
#                           http://fastjet.fr                           
#                                                                       
# Please cite EPJC72(2012)1896 [arXiv:1111.6097] if you use this package
# for scientific work and optionally PLB641(2006)57 [hep-ph/0512210].   
#								      	   
# FastJet is provided without warranty under the terms of the GNU GPLv2.
# It uses T. Chan's closest pair algorithm, S. Fortune's Voronoi code
# and 3rd party plugin jet algorithms. See COPYING file for details.
#--------------------------------------------------------------------------
#########
## Example of different trimming methods
#########
Anti-kT jet parameters: R_jet=1, pTcut=200
Trimming parameters: R_sub=0.2, fcut=0.05
Untrimmed Mass = 39.9912
Tree Trimmed Mass = 22.6675
Jet Shape Trimmed Mass = 25.4948
Event Shape Trimmed Mass = 25.4948
#####
#########
## Example of variable pT analysis
#########
Jet parameters: R_jet=1, pTcut=200 (unless otherwise stated)
Trimming parameters: R_sub=0.2, fcut=0.05
#####
N_jet from full pT_cut function=1.99971
N_jet_trim from full pT_cut function=1.97156
#####
Variable pT_cut
pT_cut=10, N_jet=4.57942
pT_cut=20, N_jet=2.98989
pT_cut=50, N_jet=2.96953
pT_cut=100, N_jet=1.99971
pT_cut=150, N_jet=1.99971
#####
pT of the hardest jet=983.637
pT of the hardest trimmed jet=983.637
#####
#########
## Example of variable R analysis
#########
Jet parameters: R_jet=1 (unless otherwise stated), pTcut=200
Trimming parameters: R_sub=0.2, fcut=0.05
#####
N_jet from full R function=1.99971
N_jet_trim from full R function=1.97156
#####
Variable Rjet
R=0.2, N_jet=1.99997
R=0.4, N_jet=2.00011
R=0.6, N_jet=1.99988
R=0.8, N_jet=2.00023
R=1, N_jet=1.99971
#########
## Hardest jet axis and weights
#########
Without global consistency:
Rap_hardest_jet=-0.868431
Phi_hardest_jet=2.90578
pt_weight_hardest_jet=983.637
N_jet_weight_hardest_jet=0.999981
#
With global consistency:
Rap_hardest_jet=-0.868431
Phi_hardest_jet=2.90578
pt_weight_hardest_jet=983.637
N_jet_weight_hardest_jet=0.999981
#########
## Examples showing more general shapes
#########
Jet parameters: R_jet=1, pTcut=200
Trimming parameters: R_sub=0.2, fcut=0.05
#####
Summed Jet Count as event shape, R_jet=1, pT_cut=200
Njet=1.99971
(Storage array works? Yes)
#
Summed Jet Count as event shape, R_jet=1, pT_cut=200, trimming with R_sub=0.2, fcut=0.05
Njet_trim=1.97156
(Storage array works? Yes)
#
Summed Jet Scalar Pt as event shape, R_jet=1, pT_cut=200
HT=1893.67
(Storage array works? Yes)
#
Summed Jet Scalar Pt as event shape, R_jet=1, pT_cut=200, trimming with R_sub=0.2, fcut=0.05
HT_trim=1867.45
(Storage array works? Yes)
#
Missing Transverse Momentum as event shape, R_jet=1, pT_cut=200
MissHT=76.1162
(Storage array works? Yes)
#
Missing Transverse Momentum as event shape, R_jet=1, pT_cut=200, trimming with R_sub=0.2, fcut=0.05
MissHT_trim=91.8116
(Storage array works? Yes)
#
Summed Jet Mass as event shape, R_jet=1, pT_cut=200
SigmaM=129.472
(Storage array works? Yes)
#
Summed Jet Mass as event shape, R_jet=1, pT_cut=200, trimming with R_sub=0.2, fcut=0.05
SigmaM_trim=125.566
(Storage array works? Yes)
#
Summed Jet Mass^2 as event shape, R_jet=1, pT_cut=200
SigmaM2=9852.27
(Storage array works? Yes)
#
Summed Jet Mass^2 as event shape, R_jet=1, pT_cut=200, trimming with R_sub=0.2, fcut=0.05
SigmaM2_trim=9114.19
(Storage array works? Yes)
#
Trimmed Subjet multiplicity with R_jet=1, pT_cut=200, R_sub=0.2, and fcut=0.05
SigmaNsub_trim=2.04888
(Storage array works? Yes)
#
Summed Jet Scalar Pt^-2 as event shape, R_jet=1, pT_cut=200
GenHT_-2=2.24003e-06
(Storage array works? Yes)
#
Summed Jet Scalar Pt^-2 as event shape, R_jet=1, pT_cut=200, trimming with R_sub=0.2, fcut=0.05
GenHT_-2_trim=2.2075e-06
(Storage array works? Yes)
#
Summed Jet Scalar Pt^-1 as event shape, R_jet=1, pT_cut=200
GenHT_-1=0.00211487
(Storage array works? Yes)
#
Summed Jet Scalar Pt^-1 as event shape, R_jet=1, pT_cut=200, trimming with R_sub=0.2, fcut=0.05
GenHT_-1_trim=0.00208462
(Storage array works? Yes)
#
Summed Jet Scalar Pt^0 as event shape, R_jet=1, pT_cut=200
GenHT_0=1.99971
(Storage array works? Yes)
#
Summed Jet Scalar Pt^0 as event shape, R_jet=1, pT_cut=200, trimming with R_sub=0.2, fcut=0.05
GenHT_0_trim=1.97156
(Storage array works? Yes)
#
Summed Jet Scalar Pt^1 as event shape, R_jet=1, pT_cut=200
GenHT_1=1893.67
(Storage array works? Yes)
#
Summed Jet Scalar Pt^1 as event shape, R_jet=1, pT_cut=200, trimming with R_sub=0.2, fcut=0.05
GenHT_1_trim=1867.45
(Storage array works? Yes)
#
Summed Jet Scalar Pt^2 as event shape, R_jet=1, pT_cut=200
GenHT_2=1.79595e+06
(Storage array works? Yes)
#
Summed Jet Scalar Pt^2 as event shape, R_jet=1, pT_cut=200, trimming with R_sub=0.2, fcut=0.05
GenHT_2_trim=1.77151e+06
(Storage array works? Yes)
#
#####