#include "MvdPid/PndMvdPidCalc.h"
#include "MvdData/PndMvdPidCand.h"

void PndMvdPidCalc::CalcMean(PndMvdPidCand* cand) {
  cand->fmeanMomentum=0;
  Double_t dE=0;
  Double_t dx=0;
  for (Int_t k=0;k<cand->fdE.size(); k++) {
    dE+=cand->fdE[k];
    dx+=cand->fdx[k];
    cand->fmeanMomentum+=cand->fmomentum[k];
  }
  cand->fmeanEnergyloss=dE/dx;
  cand->fmeanMomentum/=cand->fdE.size();
}