SampleHandlerNuDockBase.h

Go to the documentation of this file.

 
#include "Parameters/ParameterHandlerGeneric.h"
#include "Samples/SampleHandlerInterface.h"
#include "NuDockFactory.h"
_MaCh3_Safe_Include_Start_ //{
#include "nudock.hpp"
_MaCh3_Safe_Include_End_ //}
 
class SampleHandlerNuDockBase : public SampleHandlerInterface {
public:
  SampleHandlerNuDockBase(std::string configFile, ParameterHandlerGeneric* xsec_cov);
 
  virtual ~SampleHandlerNuDockBase();
 
  void Reweight() override;
 
  double GetLikelihood() const override;
 
  std::string GetSampleTitle([[maybe_unused]] const int Sample) const override { return "NuDockSample"; };
 
  std::string GetName() const override { return "NuDockSample"; };
 
  double GetSampleLikelihood([[maybe_unused]] const int isample) const override { return GetLikelihood(); };
 
  void PrintRates([[maybe_unused]] const bool DataOnly = false) override { MACH3LOG_INFO("No rates to print for NuDock sample handler"); };
 
  int GetNOscChannels([[maybe_unused]] const int iSample) const override { return 0; };
 
  std::string GetKinVarName([[maybe_unused]] const int iSample, [[maybe_unused]] const int Dimension) const override { return ""; };
 
  const TH1* GetDataHist([[maybe_unused]] const int Sample) override { return nullptr; };
  const TH1* GetMCHist([[maybe_unused]] const int Sample) override { return nullptr; };
  const TH1* GetW2Hist([[maybe_unused]] const int Sample) override { return nullptr; };
 
  int GetNDim([[maybe_unused]] const int Sample) const override { return 0; };
  std::string GetFlavourName([[maybe_unused]] const int iSample, [[maybe_unused]] const int iChannel) const override { return ""; };
 
  std::vector<double> ReturnKinematicParameterBinning([[maybe_unused]] const int Sample, [[maybe_unused]] const std::string &KinematicParameter) const override { return {}; };
  std::unique_ptr<TH1> Get1DVarHistByModeAndChannel([[maybe_unused]] const int iSample, [[maybe_unused]] const std::string& ProjectionVar_Str,
                                            [[maybe_unused]] int kModeToFill = -1, [[maybe_unused]] int kChannelToFill = -1,
                                            [[maybe_unused]] int WeightStyle = 0) override { return nullptr; };
  std::unique_ptr<TH2> Get2DVarHistByModeAndChannel([[maybe_unused]] const int iSample, [[maybe_unused]] const std::string& ProjectionVar_StrX,
                                            [[maybe_unused]] const std::string& ProjectionVar_StrY, [[maybe_unused]] int kModeToFill = -1,
                                            [[maybe_unused]] int kChannelToFill = -1, [[maybe_unused]] int WeightStyle = 0) override { return nullptr;};
  std::unique_ptr<TH1> Get1DVarHist([[maybe_unused]] const int iSample, [[maybe_unused]] const std::string &ProjectionVar,
                           [[maybe_unused]] const std::vector<KinematicCut> &EventSelectionVec = {}, [[maybe_unused]] int WeightStyle = 0,
                           [[maybe_unused]] const std::vector<KinematicCut> &SubEventSelectionVec = {}) override {return nullptr;};
  std::unique_ptr<TH2> Get2DVarHist([[maybe_unused]] const int iSample, [[maybe_unused]] const std::string& ProjectionVarX,
                                    [[maybe_unused]] const std::string& ProjectionVarY, [[maybe_unused]] const std::vector< KinematicCut >& EventSelectionVec = {},
                                    [[maybe_unused]] int WeightStyle = 0, [[maybe_unused]] const std::vector< KinematicCut >& SubEventSelectionVec = {}) override { return nullptr; };
 
protected:
  void CleanMemoryBeforeFit() override {};
 
  void Init();
 
  std::unique_ptr<NuDock> nudock_ptr;
 
  std::unique_ptr<Manager> SampleManager;
 
  bool verbose;
 
  std::vector<int> nudockParamInds;
 
  ParameterHandlerGeneric* ParHandler;
};