new (particlesR[nump]) TParticle(it->Encoding(), type, //pdg,stat
im1, im2, id1, id2, //m1,m2,d1,d2
mom[0], mom[1], mom[2], TMath::Sqrt(mom.Mag2() + m1 * m1), //px,py,pz,e
new (particlesR[nump]) TParticle(it->Encoding(), type, //pdg,stat
im1, im2, id1, id2, //m1,m2,d1,d2
mom[0], mom[1], mom[2], TMath::Sqrt(mom.Mag2() + m1 * m1), //px,py,pz,e