Figure out pid for large ecut / vcut

[apizzuto]

For now, we are just going to put a recommended max ecut / vcut, but the problem seems to be the line decay_products = [p for i,p in zip(range(max(len(particles)-3,0),len(particles)), particles[-3:]) if int(p.type) <= 1000000001] for large ecut, because it just outputs a handful of particles (mostly nutau and pions?)