Double_t ACAna::MT(TLorentzVector *particle[], Int_t Number) { TLorentzVector *sparticle = new TLorentzVector(0,0,0,0); for (Int_t i=0;i<Number;++i) { *sparticle += *(particle[i]); } Double_t Mt=sparticle->Mt(); delete sparticle; return Mt; //Double_t Minvsq = sE*sE - sPx*sPx - sPy*sPy- sPz*sPz; //return Minvsq >= 0. ? sqrt(Minvsq) : -sqrt(-Minvsq); }