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