CDCDedxLikelihood.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 CDCDEDXLIKELIHOOD_H
#define CDCDEDXLIKELIHOOD_H
 
#include <framework/datastore/RelationsObject.h>
#include <framework/gearbox/Const.h>
 
#include <cmath>
 
namespace Belle2 {
  class CDCDedxLikelihood : public RelationsObject {
  public:
    CDCDedxLikelihood(): RelationsObject()
    {
      //for all particles
      for (unsigned int i = 0; i < Const::ChargedStable::c_SetSize; i++)
        m_cdcLogl[i] = 0.0;
    }
 
    explicit CDCDedxLikelihood(const double* cdcLogl): RelationsObject()
    {
      //for all particles
      for (unsigned int i = 0; i < Const::ChargedStable::c_SetSize; i++)
        m_cdcLogl[i] = cdcLogl[i];
    }
 
    double getLogL(const Const::ChargedStable& type) const { return m_cdcLogl[type.getIndex()]; }
 
    double getCDCLikelihood(const Const::ChargedStable& type) const { return exp((double)m_cdcLogl[type.getIndex()]); }
 
  private:
    double m_cdcLogl[Const::ChargedStable::c_SetSize]; 
    ClassDef(CDCDedxLikelihood, 1); 
  };
}
#endif