VXDDedxLikelihood.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2012 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Jake Bennett, Christian Pulvermacher
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
 
#ifndef VXDDEDXLIKELIHOOD_H
#define VXDDEDXLIKELIHOOD_H
 
#include <framework/datastore/RelationsObject.h>
#include <framework/gearbox/Const.h>
 
#include <cmath>
 
namespace Belle2 {
  class VXDDedxLikelihood : public RelationsObject {
  public:
    VXDDedxLikelihood(): RelationsObject()
    {
      //for all particles
      for (unsigned int i = 0; i < Const::ChargedStable::c_SetSize; i++)
        m_vxdLogl[i] = 0.0;
    }
 
    explicit VXDDedxLikelihood(const double* vxdLogl): RelationsObject()
    {
      //for all particles
      for (unsigned int i = 0; i < Const::ChargedStable::c_SetSize; i++)
        m_vxdLogl[i] = vxdLogl[i];
    }
 
    double getLogL(const Const::ChargedStable& type) const { return m_vxdLogl[type.getIndex()]; }
 
    double getVXDLikelihood(const Const::ChargedStable& type) const { return exp((double)m_vxdLogl[type.getIndex()]); }
 
  private:
    double m_vxdLogl[Const::ChargedStable::c_SetSize]; 
    ClassDef(VXDDedxLikelihood, 1); 
  };
}
#endif