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(const int Sample) const override { (void)Sample; return "NuDockSample"; };
 
  std::string GetName() const override { return "NuDockSample"; };
 
  double GetSampleLikelihood(const int isample) const override { (void)isample; return GetLikelihood(); };
 
  void PrintRates(const bool DataOnly = false) override { (void)DataOnly; MACH3LOG_INFO("No rates to print for NuDock sample handler"); };
 
  int GetNOscChannels(const int iSample) const override { (void)iSample; return 0; };
 
  // Functions for posterior predictive - left unimplemented in the base class since they are not needed for likelihood evaluation and may require experiment-specific handling
  std::string GetKinVarName(const int iSample, const int Dimension) const override { (void)iSample; (void)Dimension; return ""; };
 
  const TH1* GetDataHist(const int Sample) override { (void)Sample; return nullptr; };
  const TH1* GetMCHist(const int Sample) override { (void)Sample; return nullptr; };
  const TH1* GetW2Hist(const int Sample) override { (void)Sample; return nullptr; };
 
  int GetNDim(const int Sample) const override { (void)Sample; return 0; };
  std::string GetFlavourName(const int iSample, const int iChannel) const override { (void)iSample; (void)iChannel; return ""; };
 
  std::vector<double> ReturnKinematicParameterBinning(const int Sample, const std::string &KinematicParameter) const override { (void)Sample; (void)KinematicParameter; return {}; };
 
  std::unique_ptr<TH1> Get1DVarHistByModeAndChannel(const int iSample, const std::string& ProjectionVar_Str,
                                            int kModeToFill = -1, int kChannelToFill = -1, int WeightStyle = 0) override {
    (void)iSample; (void)ProjectionVar_Str; (void)kModeToFill; (void)kChannelToFill; (void)WeightStyle;
    return nullptr;
  };
 
  std::unique_ptr<TH2> Get2DVarHistByModeAndChannel(const int iSample, const std::string& ProjectionVar_StrX,
                                            const std::string& ProjectionVar_StrY, int kModeToFill = -1,
                                            int kChannelToFill = -1, int WeightStyle = 0) override {
    (void)iSample; (void)ProjectionVar_StrX; (void)ProjectionVar_StrY; (void)kModeToFill; (void)kChannelToFill; (void)WeightStyle;
    return nullptr;
  };
 
  std::unique_ptr<TH1> Get1DVarHist(const int iSample, const std::string &ProjectionVar,
                           const std::vector<KinematicCut> &EventSelectionVec = {}, int WeightStyle = 0,
                           const std::vector<KinematicCut> &SubEventSelectionVec = {}) override {
    (void)iSample; (void)ProjectionVar; (void)EventSelectionVec; (void)WeightStyle; (void)SubEventSelectionVec;
    return nullptr;
  };
 
  std::unique_ptr<TH2> Get2DVarHist(const int iSample, const std::string& ProjectionVarX, const std::string& ProjectionVarY,
                            const std::vector< KinematicCut >& EventSelectionVec = {},
                            int WeightStyle = 0, const std::vector< KinematicCut >& SubEventSelectionVec = {}) override {
    (void)iSample; (void)ProjectionVarX; (void)ProjectionVarY; (void)EventSelectionVec; (void)WeightStyle; (void)SubEventSelectionVec;
    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;
};