DataDriven.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2016 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Thomas Keck                                              *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
 
#pragma once
#ifndef INCLUDE_GUARD_BELLE2_MVA_SPLOT_HEADER
#define INCLUDE_GUARD_BELLE2_MVA_SPLOT_HEADER
 
#include <mva/interface/Dataset.h>
#include <mva/utility/Binning.h>
#include <mva/interface/Expert.h>
 
 
namespace Belle2 {
  namespace MVA {
 
    class ReweightingDataset : public Dataset {
 
    public:
      ReweightingDataset(const GeneralOptions& general_options, Dataset& dataset, const std::vector<float>& weights);
 
      virtual unsigned int getNumberOfFeatures() const override { return m_dataset.getNumberOfFeatures(); }
 
      virtual unsigned int getNumberOfSpectators() const override { return m_dataset.getNumberOfSpectators(); }
 
      virtual unsigned int getNumberOfEvents() const override { return m_dataset.getNumberOfEvents(); };
 
      virtual std::vector<float> getFeature(unsigned int iFeature) override { return m_dataset.getFeature(iFeature); }
 
      virtual std::vector<float> getSpectator(unsigned int iSpectator) override { return m_dataset.getSpectator(iSpectator); }
 
      virtual std::vector<float> getWeights() override
      {
        auto weights = m_dataset.getWeights();
        for (unsigned int iEvent = 0; iEvent < weights.size(); ++iEvent) {
          weights[iEvent] *= m_weights[iEvent];
        }
        return weights;
      }
 
      virtual std::vector<float> getTargets() override { return m_dataset.getTargets(); }
 
      virtual std::vector<bool> getSignals() override { return m_dataset.getSignals(); }
 
      virtual void loadEvent(unsigned int event) override;
 
    private:
      Dataset& m_dataset; 
      std::vector<float> m_weights; 
    };
 
    class SidebandDataset : public Dataset {
 
    public:
      SidebandDataset(const GeneralOptions& general_options, Dataset& dataset, Dataset& mc_dataset, const std::string& sideband_variable);
 
      virtual unsigned int getNumberOfFeatures() const override { return m_dataset.getNumberOfFeatures(); }
 
      virtual unsigned int getNumberOfSpectators() const override { return m_dataset.getNumberOfSpectators(); }
 
      virtual unsigned int getNumberOfEvents() const override { return m_dataset.getNumberOfEvents(); };
 
      virtual std::vector<float> getFeature(unsigned int iFeature) override { return m_dataset.getFeature(iFeature); }
 
      virtual std::vector<float> getSpectator(unsigned int iSpectator) override { return m_dataset.getSpectator(iSpectator); }
 
      virtual void loadEvent(unsigned int event) override;
 
    private:
      Dataset& m_dataset; 
      int m_spectator_index; 
      double m_signal_weight; 
      double m_background_weight; 
      double m_negative_signal_weight; 
    };
 
    class SPlotDataset : public Dataset {
 
    public:
      SPlotDataset(const GeneralOptions& general_options, Dataset& dataset, const std::vector<float>& weights, float signalFraction);
 
      virtual unsigned int getNumberOfFeatures() const override { return m_dataset.getNumberOfFeatures(); }
 
      virtual unsigned int getNumberOfSpectators() const override { return m_dataset.getNumberOfSpectators(); }
 
      virtual unsigned int getNumberOfEvents() const override { return 2 * m_dataset.getNumberOfEvents(); };
 
      virtual void loadEvent(unsigned int event) override;
 
      virtual float getSignalFraction() override;
 
    private:
      Dataset& m_dataset; 
      std::vector<float> m_weights; 
      float m_signalFraction; 
    };
 
    std::vector<float> getSPlotWeights(Dataset& dataset, Binning& binning);
 
    std::vector<float> getBoostWeights(Dataset& dataset, Binning& binning);
 
    std::vector<float> getAPlotWeights(Dataset& dataset, Binning& binning, const std::vector<float>& boost_prediction);
 
  }
}
#endif