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/StoreArray.h>
 
#include <klm/dataobjects/KLMHit2d.h>
#include <klm/dataobjects/KLMMuidLikelihood.h>
#include <klm/dataobjects/bklm/BKLMElementNumbers.h>
#include <klm/dataobjects/eklm/EKLMElementNumbers.h>
 
#include <mdst/dataobjects/Track.h>
 
#include <reconstruction/dataobjects/KLMMuonIDDNNInputVariable.h>
 
#include <mva/dataobjects/DatabaseRepresentationOfWeightfile.h>
#include <mva/interface/Expert.h>
#include <mva/interface/Weightfile.h>
 
#include <string>
 
namespace Belle2 {
  class KLMMuonIDDNNExpertModule : public Module {
 
  public:

    KLMMuonIDDNNExpertModule();

    ~KLMMuonIDDNNExpertModule();

    void initialize() override;

    void beginRun() override;

    void event() override;

    void endRun() 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;

    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;
 
  };

}