PndMvdIdealPidAlgo.h

Go to the documentation of this file.

#ifndef PndMvdIdealPidAlgo_h
#define PndMvdIdealPidAlgo_h

#include <math.h>
#include <map>
#include "Rtypes.h"

class PndMvdPidCand;

class PndMvdIdealPidAlgo {
public:

  //Write likelihoods to PndMvdPidCand
  static void CalcLikelihood(PndMvdPidCand* cand);

private:

  //Lower boundary of energy loss distribution, as given by Bethe-Bloch formula
  static double LowerBoundary(double p, double m);

  //These functions represent the lower boundary of the energy loss distribution
  //at a given momentum. The offset numbers in both momentum and energy loss
  //axis are adjusted that at least 99% of all events have are gathered.
  static double LowerProtonBoundary(double momentum);

  static double LowerKaonBoundary(double momentum);

  static double LowerPionBoundary(double momentum);

  //Naive approximation of the Landau distribution, integrated from -inf to x
  static double LandauIntegral(double x);

  //Constants
  static float fpiMass;
  static float fkMass;
  static float fpMass;
  static float feMass;
  static float fc;
  static float feb;

  //Landau distribution parameters
  static double fkShift;
  static double fkScale;

  static double fpiShift;
  static double fpiScale;

  //Integral remainders
  static double fPRemainder;
  static double fKRemainder;
  static double fPiRemainder;
};

#endif