TOPPixelLikelihood.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
 
#include <framework/datastore/RelationsObject.h>
#include <framework/gearbox/Const.h>
#include <framework/logging/Logger.h>
 
#include <map>
#include <vector>
#include <tuple>
#include <utility>
 
namespace Belle2 {
  class TOPPixelLikelihood : public RelationsObject {
 
  public:
 
    typedef std::array<float, 512> PixelArray_t;
 
    TOPPixelLikelihood() { }
 
    bool addHypothesisLikelihoods(const PixelArray_t& plkhs, const int hypothesis)
    {
      auto result = m_logls.insert(std::make_pair(hypothesis, plkhs));
      if (not result.second) {
        B2WARNING("Likelihoods already exist for this track");
      }
      return result.second;
    }
 
    bool addHypothesisSignalPhotons(const PixelArray_t& sfots, const int hypothesis)
    {
      auto result = m_sfots.insert(std::make_pair(hypothesis, sfots));
      if (not result.second) {
        B2WARNING("Signal photons already exist for this track");
      }
      return result.second;
    }
 
    const PixelArray_t& getHypothesisLikelihoods(const int hypothesis) const
    {
      return m_logls.at(hypothesis);
    }
 
    const PixelArray_t& getHypothesisSignalPhotons(const int hypothesis) const
    {
      return m_sfots.at(hypothesis);
    }
 
    const PixelArray_t& getPixelLogL_e() const { return m_logls.at(11); }
 
    const PixelArray_t& getPixelLogL_mu() const { return m_logls.at(13); }
 
    const PixelArray_t& getPixelLogL_pi() const { return m_logls.at(211); }
 
    const PixelArray_t& getPixelLogL_K() const { return m_logls.at(321); }
 
    const PixelArray_t& getPixelLogL_p() const { return m_logls.at(2212); }
 
    const PixelArray_t& getPixelSigPhot_e() const { return m_sfots.at(11); }
 
    const PixelArray_t& getPixelSigPhot_mu() const { return m_sfots.at(13); }
 
    const PixelArray_t& getPixelSigPhot_pi() const { return m_sfots.at(211); }
 
    const PixelArray_t& getPixelSigPhot_K() const { return m_sfots.at(321); }
 
    const PixelArray_t& getPixelSigPhot_p() const { return m_sfots.at(2212); }
 
    void setModuleID(int moduleID)
    {
      m_moduleID = moduleID;
    }
 
    int getModuleID() const { return m_moduleID; }
 
  private:
    std::map<int, PixelArray_t> m_logls; 
    std::map<int, PixelArray_t> m_sfots; 
    int m_moduleID = 0; 
    ClassDef(TOPPixelLikelihood, 3); 
  };
} // end namespace Belle2