MuidBuilder.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2010 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Leo Piilonen                                             *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
 
#pragma once
 
/* KLM headers. */
#include <klm/dbobjects/KLMLikelihoodParameters.h>
#include <klm/muid/MuidElementNumbers.h>
 
/* Belle 2 headers. */
#include <framework/database/DBObjPtr.h>
 
namespace Belle2 {
  class KLMMuidLikelihood;
 
  class MuidBuilder {
 
  public:
 
    explicit MuidBuilder(int pdg);
 
    ~MuidBuilder();
 
    double getPDF(const KLMMuidLikelihood* muid) const;
 
  private:
 
    MuidBuilder();
 
    MuidBuilder(MuidBuilder&);
 
    MuidBuilder& operator=(const MuidBuilder&);
 
    void fillPDFs(MuidElementNumbers::Hypothesis hypothesis);
 
    void spline(int n, double dx, double Y[], double B[], double C[], double D[]);
 
    double getLongitudinalPDF(const KLMMuidLikelihood* muid) const;
 
    double getTransversePDF(const KLMMuidLikelihood* muid) const;
 
    double m_LayerPDF[MuidElementNumbers::getMaximalOutcome() + 1][MuidElementNumbers::getMaximalBarrelLayer() +
        1][MuidElementNumbers::getMaximalBarrelLayer() + MuidElementNumbers::getMaximalEndcapForwardLayer() + 2];
 
    double m_ReducedChiSquaredThreshold[MuidElementNumbers::getMaximalDetector() + 1][MuidElementNumbers::getMaximalHalfNdof() + 1];
 
    double m_ReducedChiSquaredScaleX[MuidElementNumbers::getMaximalDetector() + 1][MuidElementNumbers::getMaximalHalfNdof() + 1];
 
    double m_ReducedChiSquaredScaleY[MuidElementNumbers::getMaximalDetector() + 1][MuidElementNumbers::getMaximalHalfNdof() + 1];
 
    double m_ReducedChiSquaredPDF[MuidElementNumbers::getMaximalDetector() + 1][MuidElementNumbers::getMaximalHalfNdof() +
        1][MuidElementNumbers::getSizeReducedChiSquared()];
 
    double m_ReducedChiSquaredD1[MuidElementNumbers::getMaximalDetector() + 1][MuidElementNumbers::getMaximalHalfNdof() +
        1][MuidElementNumbers::getSizeReducedChiSquared()];
 
    double m_ReducedChiSquaredD2[MuidElementNumbers::getMaximalDetector() + 1][MuidElementNumbers::getMaximalHalfNdof() +
        1][MuidElementNumbers::getSizeReducedChiSquared()];
 
    double m_ReducedChiSquaredD3[MuidElementNumbers::getMaximalDetector() + 1][MuidElementNumbers::getMaximalHalfNdof() +
        1][MuidElementNumbers::getSizeReducedChiSquared()];
 
    double m_ReducedChiSquaredDx;
 
    DBObjPtr<KLMLikelihoodParameters> m_LikelihoodParameters;
 
  };
 
}