def main():
fj.ClusterSequence.print_banner()
print()
sconfig_pythia = [
"Beams:eCM = 8000.", "HardQCD:all = on", "PhaseSpace:pTHatMin = 20."
]
pythia = create_and_init_pythia(sconfig_pythia)
jet_R0 = 0.4
jet_def = fj.JetDefinition(fj.antikt_algorithm, jet_R0)
jet_selector = fj.SelectorPtMin(20.0) & fj.SelectorPtMax(
40.0) & fj.SelectorAbsEtaMax(1)
# sd = rt.SoftDrop(0, 0.1, 1.0)
all_jets = []
for iEvent in tqdm(range(1000), 'event'):
if not pythia.next():
continue
parts = pyfj.vectorize(pythia, True, -1, 1, True)
jets = jet_selector(jet_def(parts))
all_jets.extend(jets)
print_jets(all_jets)
pythia = create_and_init_pythia(sconfig_pythia)
# print the banner first
fj.ClusterSequence.print_banner()
print()
# set up our jet definition and a jet selector
jet_R0 = 0.4
jet_def = fj.JetDefinition(fj.antikt_algorithm, jet_R0)
jet_selector = fj.SelectorPtMin(100.0) & fj.SelectorPtMax(
200.0) & fj.SelectorAbsEtaMax(1)
sd = rt.SoftDrop(0, 0.1, 1.0)
all_jets = []
for iEvent in tqdm(range(10), 'event'):
if not pythia.next(): continue
parts = pyfj.vectorize(pythia, True, -1, 1, False)
jets = jet_selector(jet_def(parts))
all_jets.extend(jets)
all_sd_jets = [sd.result(j) for j in all_jets]
pts = [j.pt() for j in all_jets]
sd_pts = [j.pt() for j in all_sd_jets]
sd_delta_pt = [delta for delta in deltas(all_jets, all_sd_jets)]
nangs0 = [pyfj.angularity(j, 0.) for j in all_jets]
jet_def_lund = fj.JetDefinition(fj.cambridge_algorithm, 1.0)
lund_gen = lund.LundGenerator(jet_def_lund)
lunds = [lund_gen.result(j) for j in all_jets]
print('listing lund plane points... Delta, kt - for {} selected jets'.format(
len(all_jets)))