PndMvdAdvancedPidAlgo.cxx File Reference

#include "TMath.h"
#include "TRandom3.h"
#include "PndMvdAdvancedPidAlgo.h"
#include "PndMvdPidCand.h"

Go to the source code of this file.

Functions
	ClassImp (PndMvdAdvancedPidAlgo) void PndMvdAdvancedPidAlgo

Function Documentation

ClassImp

(

PndMvdAdvancedPidAlgo

)

Definition at line 8 of file PndMvdAdvancedPidAlgo.cxx.

References dE, and dx.

                                                              {
  //computes pid likelihoods for pi, k, p, mu, e and stores them in the likelihood
  //map with the map key used as pdg id. This is basically done by integrating the
  //respective energy loss distribution within the chosen selector limits

  double weight[kNPidType];
  double dE=0;
  double dx=0;
  fmomentum=0;

  for (int k=0;k<cand->GetMvdHits(); k++) {
    dE+=cand->GetMvdHitdE(k);
    dx+=cand->GetMvdHitdx(k);
    fmomentum+=cand->GetMvdHitMomentum(k);
  }

  if (dx>0) {
    fmomentum/=cand->GetMvdHits();
    fenergyloss=dE/dx;
  } else
    fmomentum=-1; // whateva

  CalcLikelihood(weight);

  cand->SetLikelihood(2212, weight[kProton]);
  cand->SetLikelihood(211, weight[kPion]);
  cand->SetLikelihood(321, weight[kKaon]);
  cand->SetLikelihood(13, weight[kMuon]);
  cand->SetLikelihood(11, weight[kElectron]);
}