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