def mtrack(V, jet): #Inv. mass using all tracks in iso and core region, assuming each with pi0 mass mpi0=0.139 jetsum = LParticle() for i in range(len(V)): if (track_check(V[i])): r = delR(V[i], jet) if (r < 0.4): E = sqrt(mpi0**2 + V[i].x()**2 + V[i].y()**2 + V[i].z()**2) W = LParticle(V[i].x(), V[i].y(), V[i].z(), E) jetsum.add(W) ### print jetsum.calcMass() return jetsum.calcMass()
l = len(tau_d) if (l > 1): #when tau starts decaying for jj in range(len(tau_d)): daughter = tau_d[jj].getPDG() if ((abs(daughter)==15)): photorecoil=True # print "photon recoil" break # tau just emits photon if ((abs(daughter)>=11) and (abs(daughter)<=14)): flaglep=True #looking for lepton decay tau_lepdk = tau_lepdk + 1 break if ((photorecoil==False)): taucount=taucount+1 if ((flaglep==False)and(photorecoil==False)): leptau=LParticle('tauhad-dir') leptau.setPxPyPz(par.getMomentum()[0], par.getMomentum()[1], par.getMomentum()[2]) truejet=[] for aa in range(len(tau_d)): daughtau=tau_d[aa] if ((abs(daughtau.getPDG()))!=16): if (((abs(daughtau.getPDG()))!=113) and ((abs(daughtau.getPDG()))!=213) and ((abs(daughtau.getPDG()))!=323) and ((abs(daughtau.getPDG()))!=313)): truejetmem=LParticle(str(daughtau.getPDG())) truejetmem.setPxPyPzE(daughtau.getMomentum()[0], daughtau.getMomentum()[1], daughtau.getMomentum()[2], daughtau.getEnergy()) truejetmem.setCharge(daughtau.getCharge()) truejet.append(truejetmem) else: for reso in daughtau.getDaughters(): truejetmem=LParticle(str(reso.getPDG())) truejetmem.setPxPyPzE(reso.getMomentum()[0], reso.getMomentum()[1], reso.getMomentum()[2], reso.getEnergy()) truejetmem.setCharge(reso.getCharge())
tau_d = par.getDaughters() l = len(tau_d) if (l > 1): #when tau starts decaying for jj in range(len(tau_d)): daughter = tau_d[jj].getPDG() if ((abs(daughter)==15)): photorecoil=True # print "photon recoil" break # tau just emits photon if ((abs(daughter)>=11) and (abs(daughter)<=14)): flaglep=True #looking for lepton decay break # else: print "D: ", daughter if ((flaglep==False)and(photorecoil==False)): # print par.getMomentum()[0], par.getMomentum()[1], par.getMomentum()[2] leptau=LParticle('tauhad-dir') leptau.setPxPyPz(par.getMomentum()[0], par.getMomentum()[1], par.getMomentum()[2]) truejet=[] for aa in range(len(tau_d)): daughtau=tau_d[aa] if ((abs(daughtau.getPDG()))!=16): truejetmem=LParticle(str(daughtau.getPDG())) truejetmem.setPxPyPzE(daughtau.getMomentum()[0], daughtau.getMomentum()[1], daughtau.getMomentum()[2], daughtau.getEnergy()) truejetmem.setCharge(daughtau.getCharge()) truejet.append(truejetmem) # print leptau if ((PTCUT(truejet)) and (prong(truejet)==N)): FCENT_t = fcent(truejet, leptau) FTRACK_t = ftrack(truejet, leptau) RTRACK_t = Rtrack(truejet, leptau) NISO_t = Niso_track(truejet, leptau)