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)):
FCENT_t = fcent(truejet, leptau)
