MVAFilter.dcl.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 <tracking/trackFindingCDC/filters/base/FilterOnVarSet.dcl.h>
#include <tracking/trackFindingCDC/filters/base/Filter.fwd.h>
 
#include <tracking/trackFindingCDC/varsets/BaseVarSet.h>
 
#include <tracking/trackFindingCDC/numerics/Weight.h>
 
#include <framework/database/DBObjPtr.h>
#include <tracking/dbobjects/TrackingMVAFilterParameters.h>
 
#include <memory>
#include <string>
#include <cmath>
 
namespace Belle2 {
  class ModuleParamList;
 
  namespace TrackFindingCDC {
    class MVAExpert;

    template<class AFilter>
    class MVA : public OnVarSet<AFilter> {
 
    private:
      using Super = OnVarSet<AFilter>;
 
    public:
      using Object = typename AFilter::Object;
 
    private:
      using AVarSet = BaseVarSet<Object>;
 
    public:
      explicit MVA(std::unique_ptr<AVarSet> varSet,
                   const std::string& identifier = "",
                   double defaultCut = NAN,
                   const std::string& dbObjectName = "");

      virtual ~MVA();

      void exposeParameters(ModuleParamList* moduleParamList, const std::string& prefix) override;

      void initialize() override;

      void beginRun() override;
 
    public:
      Weight operator()(const Object& obj) override;

      virtual double predict(const Object& obj);

      std::vector<float> predict(const std::vector<Object*>& objs);

      virtual std::vector<float> operator()(const std::vector <Object*>& objs) override;
    private:
      std::string m_identifier = "";

      double m_cutValue;

      std::string m_DBPayloadName = "";

      std::unique_ptr<MVAExpert> m_mvaExpert;

      std::vector<Named<Float_t*>> m_namedVariables;
    };

    template<class AVarSet>
    class MVAFilter: public MVA<Filter<typename AVarSet::Object> > {
 
    private:
      using Super = MVA<Filter<typename AVarSet::Object> >;
 
    public:
      using Object = typename AVarSet::Object;

      explicit MVAFilter(const std::string& defaultTrainingName = "",
                         double defaultCut = NAN,
                         const std::string& defaultDBObjectName = "");

      ~MVAFilter();
    };
  }
}