KLMMuonIDDNNExpertModule.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/core/Module.h>
#include <framework/database/DBObjPtr.h>
#include <framework/datastore/StoreObjPtr.h>
#include <framework/datastore/StoreArray.h>
 
#include <klm/dataobjects/bklm/BKLMElementNumbers.h>
#include <klm/dataobjects/eklm/EKLMElementNumbers.h>
 
 
#include <string>
 
namespace Belle2 {
  class DatabaseRepresentationOfWeightfile;
  class KLMHit2d;
  class KLMMuidLikelihood;
  class KLMMuonIDDNNInputVariable;
  class Track;
 
  namespace MVA {
    class Weightfile;
    class Expert;
    class SingleDataset;
  }
  class KLMMuonIDDNNExpertModule : public Module {
 
  public:

    KLMMuonIDDNNExpertModule();

    ~KLMMuonIDDNNExpertModule() override;

    void initialize() override;

    void beginRun() override;

    void event() override;

    void terminate() override;
 
  private:

    void initializeMVA(MVA::Weightfile& weightfile);

    float getNNmuProbability(const Track* track, const KLMMuidLikelihood* klmll);

    float getHitWidth(const KLMHit2d* klmhit);

    StoreArray<Track> m_tracks;

    StoreArray<KLMMuonIDDNNInputVariable> m_inputVariable;

    const std::string m_identifier = "KLMMuonIDDNNWeightFile";

    DBObjPtr<DatabaseRepresentationOfWeightfile> m_weightfile_representation{m_identifier};

    std::unique_ptr<MVA::Expert> m_expert;

    std::unique_ptr<MVA::SingleDataset> m_dataset;

    static constexpr int m_maxBKLMLayers = BKLMElementNumbers::getMaximalLayerNumber();

    static constexpr int m_maxEKLMLayers = EKLMElementNumbers::getMaximalLayerNumber();

    std::array < float, m_maxBKLMLayers > m_BarrelPhiStripWidth{};

    std::array < float, m_maxBKLMLayers > m_BarrelZStripWidth{};

    float m_EndcapScintWidth = 0.;

    std::array < float, m_maxBKLMLayers + m_maxEKLMLayers > m_hitpattern_width{};

    std::array < float, m_maxBKLMLayers + m_maxEKLMLayers > m_hitpattern_steplength{};

    std::array < float, m_maxBKLMLayers + m_maxEKLMLayers > m_hitpattern_chi2{};

    std::array < bool, m_maxBKLMLayers + m_maxEKLMLayers > m_hitpattern_hasext{};
 
  };

}