ECLClusterPSD.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
 
/* Basf2 headers. */
#include <framework/core/Module.h>
#include <framework/database/DBObjPtr.h>
#include <framework/datastore/StoreArray.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