def __init__(self):
pythia8.UserHooks.__init__(self)
self.slowJet = pythia8.SlowJet(-1, 0.7, 10., 5.)
self.pTHat = 0.
R = float(sys.argv[6])
pTjetMin = int(sys.argv[7]) #Min pT for all jets
# eta_max=2.5
# R=1.5
# pTjetMin=300 #Min pT for all jets
#Exclude neutrinos (and other invisible) from study:
nSel = 2
# calorimeter granularity
etaedges = np.arange(-eta_max, eta_max + 0.025, 10 / 206.25)
phiedges = np.arange(-np.pi, np.pi + 0.024, np.pi / 130.)
# print etaedges, phiedges
my_map = plt.cm.gray
#Set up SlowJet jet finder, with anti-kT clustering and pion mass assumed for non-photons..
slowJet = pythia8.SlowJet(p_value, R, pTjetMin, eta_max, nSel, 1)
jet_list = []
jet_parton_list = []
jet_mass = []
##----------------------------------------------------------------------------------------------------------------------
# Begin event loop. Generate event. Skip if error. List first one.
for iEvent in range(0, Nevents):
if not pythia.next(): continue
# for i in range(0,pythia.event.size()):
# print('i {},id {}'.format(i,pythia.event[i].id()))
# Find number of all final charged particles and fill histogram.
nCharged = 0
for prt in pythia.event: