PndHypAdvancedPidAlgo.cxx File Reference

#include "TMath.h"
#include "TRandom3.h"
#include "PndHypAdvancedPidAlgo.h"
#include "PndHypPidCand.h"
#include <iostream>
#include "TROOT.h"

Go to the source code of this file.

Functions
	ClassImp (PndHypAdvancedPidAlgo) void PndHypAdvancedPidAlgo

Function Documentation

ClassImp

(

PndHypAdvancedPidAlgo

)

Definition at line 9 of file PndHypAdvancedPidAlgo.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->GetHypHits(); k++) {
    dE+=cand->GetHypHitdE(k);
    dx+=cand->GetHypHitdx(k);
    fmomentum+=cand->GetHypHitMomentum(k);
  }

  if (dx>0) {
    fmomentum/=cand->GetHypHits();
    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]);
}