MuidBuilder.h

/**************************************************************************
 * basf2 (Belle II Analysis Software Framework)                           *
 * Author: The Belle II Collaboration                                     *
 *                                                                        *
 * See git log for contributors and copyright holders.                    *
 * This file is licensed under LGPL-3.0, see LICENSE.md.                  *
 **************************************************************************/
 
#pragma once
 
/* KLM headers. */
#include <klm/dbobjects/KLMLikelihoodParameters.h>
#include <klm/muid/MuidElementNumbers.h>
 
/* Basf2 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);

    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;
 
  };

}