e_pi = hypot4(sim.fPxPi, sim.fPyPi, sim.fPzPi,
AliPID.ParticleMass(AliPID.kPion))
deu.SetXYZM(sim.fPxDeu, sim.fPyDeu, sim.fPzDeu,
AliPID.ParticleMass(AliPID.kDeuteron))
p.SetXYZM(sim.fPxP, sim.fPyP, sim.fPzP,
AliPID.ParticleMass(AliPID.kProton))
pi.SetXYZM(sim.fPxPi, sim.fPyPi, sim.fPzPi,
AliPID.ParticleMass(AliPID.kPion))
hyp = deu + p + pi
decay_lenght = TVector3(sim.fDecayVtxX, sim.fDecayVtxY, sim.fDecayVtxZ)
dl = decay_lenght.Mag()
if hyp.Gamma() == 0 or hyp.Beta() == 0:
continue
ct = dl / (hyp.Gamma() * hyp.Beta())
if hyp.Pt() < 1. or hyp.Pt() > 10.:
continue
hist_ctsim.Fill(ct)
hist_ptsim.Fill(hyp.Pt())
hist_psim.Fill(hyp.P())
hist_etasim.Fill(hyp.Eta())
hist_phisim.Fill(hyp.Phi())
# rec - sim diff
# if sim.fRecoIndex >= 0:
