Esempio n. 1
0
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);

}