ChargedPidMVAModule.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
 
// FRAMEWORK
#include <framework/core/Module.h>
#include <framework/dataobjects/EventMetaData.h>
#include <framework/database/DBObjPtr.h>
#include <analysis/DecayDescriptor/DecayDescriptor.h>
#include <framework/datastore/StoreArray.h>
 
// MVA
#include <mva/interface/Expert.h>
#include <mva/interface/Dataset.h>
 
// ANALYSIS
#include <analysis/VariableManager/Manager.h>
#include <analysis/dbobjects/ChargedPidMVAWeights.h>
#include <analysis/VariableManager/Utility.h>
 
 
namespace Belle2 {
  class ChargedPidMVAModule : public Module {
 
    typedef std::vector<std::unique_ptr<MVA::Expert> > ExpertsList; 
    typedef std::vector<std::unique_ptr<MVA::SingleDataset> > DatasetsList; 
    typedef std::vector<std::unique_ptr<Variable::Cut>> CutsList; 
    typedef std::vector< std::vector<const Variable::Manager::Var*> > VariablesLists; 
  public:
 
    ChargedPidMVAModule();
 
    virtual ~ChargedPidMVAModule();
 
    virtual void initialize() override;
 
    virtual void beginRun() override;
 
    virtual void event() override;
 
  private:
 
    void initializeMVA();
 
  private:
    StoreArray<Particle> m_particles;
    int m_sig_pdg;
 
    int m_bkg_pdg;
 
    std::vector<std::string> m_decayStrings;
 
    std::string m_payload_name;
 
    bool m_charge_independent;
 
    bool m_ecl_only;
 
    std::string m_score_varname;
 
    StoreObjPtr<EventMetaData> m_event_metadata;
 
    std::unique_ptr<DBObjPtr<ChargedPidMVAWeights>> m_weightfiles_representation;
 
    ExpertsList m_experts;
 
    DatasetsList m_datasets;
 
    CutsList m_cuts;
 
    VariablesLists m_variables;
 
    VariablesLists m_spectators;
 
    std::map<int, std::tuple<std::string, std::string, int>> m_stdChargedInfo = {
      { Const::electron.getPDGCode(), std::make_tuple("e", "pi", Const::pion.getPDGCode()) },
      { Const::muon.getPDGCode(), std::make_tuple("mu", "pi", Const::pion.getPDGCode()) },
      { Const::pion.getPDGCode(), std::make_tuple("pi", "K", Const::kaon.getPDGCode()) },
      { Const::kaon.getPDGCode(), std::make_tuple("K", "pi", Const::pion.getPDGCode()) },
      { Const::proton.getPDGCode(), std::make_tuple("p", "pi", Const::pion.getPDGCode()) },
      { Const::deuteron.getPDGCode(), std::make_tuple("d", "pi", Const::pion.getPDGCode()) }
    };
 
    void registerAliasesLegacy();
 
    void registerAliases();
 
  };
}