SignalPDF.cc

/**************************************************************************
 * 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.                  *
 **************************************************************************/
 
#include <top/reconstruction_cpp/SignalPDF.h>
#include <cmath>
 
using namespace std;
 
namespace Belle2 {
  namespace TOP {
 
    SignalPDF::SignalPDF(int pixelID, const TOPNominalTTS& tts):
      m_pixelID(pixelID), m_tts(&tts)
    {}
 
    double SignalPDF::getPDFValue(double time, double timeErr, double sigt) const
    {
      double wid0 = timeErr * timeErr + sigt * sigt;
      double f = 0;
      for (const auto& peak : m_peaks) {
        for (const auto& gaus : m_tts->getTTS()) {
          double sigma = peak.wid + gaus.sigma * gaus.sigma + wid0; // sigma squared!
          double x = pow(time - peak.t0 - gaus.position, 2) / sigma;
          if (x > 100) continue;
          sigma = sqrt(2 * M_PI * sigma);
          f += peak.nph * gaus.fraction / sigma * exp(-x / 2);
        }
      }
      return f;
    }
 
    double SignalPDF::getIntegral(double minTime, double maxTime) const
    {
      double s = 0;
      for (const auto& peak : m_peaks) {
        if (peak.t0 < minTime or peak.t0 > maxTime) continue;
        s += peak.nph;
      }
      return s;
    }
 
 
  } // namespace TOP
} // namespace Belle2