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()) truejet.append(truejetmem) if((prong(truejet)==1)): mctau1p = mctau1p + 1 #@#print "1p", mctau1p if((prong(truejet)==3)): mctau3p = mctau3p + 1 # print leptau if ((PTCUT(truejet)) and (prong(truejet)==N)):