ECLClusterPSD.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2018 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors:  Savino Longo (longos@uvic.ca)                           *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
 
#pragma once
 
// FRAMEWORK
#include <framework/core/Module.h>
#include <framework/datastore/StoreArray.h>
#include <framework/database/DBObjPtr.h>
 
namespace Belle2 {
  class DatabaseRepresentationOfWeightfile;
  class ECLShower;
  class ECLCalDigit;
 
  namespace MVA {
    class Expert;
    class SingleDataset;
  }
 
  class ECLClusterPSDModule: public Module {
 
  public:
 
    ECLClusterPSDModule();
 
    ~ECLClusterPSDModule();
 
    virtual void initialize() override;
 
    virtual void beginRun() override;
 
    virtual void event() override;
 
    virtual void endRun() override;
 
    virtual void terminate() override;
 
    virtual const char* eclCalDigitArrayName() const
    { return "ECLCalDigits" ; }
 
    virtual const char* eclShowerArrayName() const
    { return "ECLShowers" ; }
 
  private:
    double m_CrystalHadronEnergyThreshold; 
    double m_CrystalHadronIntensityThreshold; 
    StoreArray<ECLCalDigit> m_eclCalDigits;  
    StoreArray<ECLShower> m_eclShowers;  
    const unsigned int maxdigits = 20; 
    const unsigned int m_numMVAvariables = 160; 
    std::string m_MVAidentifier; 
    std::unique_ptr<DBObjPtr<DatabaseRepresentationOfWeightfile>>
                                                               m_weightfile_representation; 
    std::unique_ptr<MVA::Expert> m_expert; 
    std::unique_ptr<MVA::SingleDataset>  m_dataset; 
    void initializeMVAweightFile(const std::string& identifier,
                                 std::unique_ptr<DBObjPtr<DatabaseRepresentationOfWeightfile>>& weightFileRepresentation);
 
    void initializeMVA(const std::string& identifier,
                       std::unique_ptr<DBObjPtr<DatabaseRepresentationOfWeightfile>>& weightFileRepresentation, std::unique_ptr<MVA::Expert>& expert);
 
    double evaluateMVA(const ECLShower* cluster); 
  };
 
  class ECLClusterPSDPureCsIModule : public ECLClusterPSDModule {
  public:
    virtual const char* eclCalDigitArrayName() const override
    { return "ECLCalDigitsPureCsI" ; }
 
    virtual const char* eclShowerArrayName() const override
    { return "ECLShowersPureCsI" ; }
 
  };
 
} // end Belle2 namespace