PIDNeuralNetwork.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 <analysis/VariableManager/Manager.h>
#include <analysis/dbobjects/PIDNeuralNetworkParameters.h>
#include <framework/database/DBObjPtr.h>
#include <framework/logging/Logger.h>
 
#include <framework/gearbox/Const.h>
 
#define FDEEP_FLOAT_TYPE float
#include <fdeep/fdeep.hpp>
 
 
 
namespace Belle2 {
  class PIDNeuralNetwork {
 
  public:
 
    PIDNeuralNetwork():
      m_pidNeuralNetworkParametersDB(nullptr),
      m_model(nullptr)
    {
      loadParametersFromDB();
      check();
    };
 
 
    PIDNeuralNetwork(const std::string& parameterName):
      m_pidNeuralNetworkParametersName(parameterName),
      m_pidNeuralNetworkParametersDB(nullptr),
      m_model(nullptr)
    {
      loadParametersFromDB();
      check();
    }
 
    PIDNeuralNetwork(PIDNeuralNetwork&& other):
      m_pidNeuralNetworkParametersName(std::move(other.m_pidNeuralNetworkParametersName)),
      m_pidNeuralNetworkParametersDB(std::move(other.m_pidNeuralNetworkParametersDB)),
      m_model(std::move(other.m_model)),
      m_inputBasf2Names(std::move(other.m_inputBasf2Names)),
      m_extraInfoNames(std::move(other.m_extraInfoNames))
    {
    }
 
    std::map<int, double> predict(std::vector<float> inputs) const;
 
    size_t getInputSize() const { return (*m_pidNeuralNetworkParametersDB)->getInputSize(); }
 
    const std::vector<std::string>& getInputBasf2Names() const { return m_inputBasf2Names; }
 
    const std::vector<int>& getOutputSpeciesPdg() const {return (*m_pidNeuralNetworkParametersDB)->getOutputSpeciesPdg();}
 
    const std::string& getPIDNeuralNetworkParametersName() const {return m_pidNeuralNetworkParametersName;}
 
    bool hasPdgCode(const int pdg, const bool throwException = false) const {return (*m_pidNeuralNetworkParametersDB)->hasPdgCode(pdg, throwException);}
 
    const std::string& getExtraInfoName(const int pdg) const {return m_extraInfoNames.at(pdg);}
 
 
  private:
 
    void loadParametersFromDB();
    void check();
 
    std::string m_pidNeuralNetworkParametersName = "PIDNeuralNetworkParameters"; 
    std::unique_ptr<DBObjPtr<PIDNeuralNetworkParameters>> m_pidNeuralNetworkParametersDB; 
    std::unique_ptr<const fdeep::model> m_model; 
    std::vector<std::string> m_inputBasf2Names; 
    std::map<int, std::string> m_extraInfoNames; 
  };
 
} // Belle2 namespace
 
 
inline void Belle2::PIDNeuralNetwork::check()
{
  for (const auto& name : getInputBasf2Names()) {
    if (!Variable::Manager::Instance().getVariable(name))
      B2FATAL("PID neural network needs input '" + name + "', but this input is not available!");
  }
 
}