PredictiveThrower.h

Go to the documentation of this file.

#pragma once
 
#include "Fitters/FitterBase.h"
 
class ParameterHandlerGeneric;
 
 
class PredictiveThrower : public FitterBase {
 public:
  PredictiveThrower(manager * const fitMan);
  virtual ~PredictiveThrower();
 
  void ProduceToys();
 
  void RunPredictiveAnalysis();
 
  void RunMCMC() override {
    MACH3LOG_ERROR("{} is not supported in {}", __func__, GetName());
    throw MaCh3Exception(__FILE__ , __LINE__ );
  };
 
 private:
  void SetParamters();
 
  void SetupToyGeneration();
 
  bool LoadToys();
 
  void SetupSampleInformation();
 
  std::unique_ptr<TH1D> MakePredictive(const std::vector<std::unique_ptr<TH1D>>& Toys,
                                       const std::string& Sample_Name,
                                       const std::string& suffix,
                                       const bool DebugHistograms);
 
  std::vector<std::unique_ptr<TH2D>> ProduceSpectra(const std::vector<std::vector<std::unique_ptr<TH1D>>>& Toys,
                                                    const std::string suffix);
 
  void PosteriorPredictivepValue(const std::vector<std::unique_ptr<TH1D>>& PostPred_mc,
                                 //const std::vector<std::unique_ptr<TH1D>>& PostPred_w2,
                                 const std::vector<TDirectory*>& SampleDir);
 
 
  double GetLLH(const std::unique_ptr<TH1D>& DatHist,
                                   const std::unique_ptr<TH1D>& MCHist,
                                   const std::unique_ptr<TH1D>& W2Hist,
                                   SampleHandlerBase* SampleHandler);
 
  void MakeChi2Plots(const std::vector<std::vector<double>>& Chi2_x,
                     const std::string& Chi2_x_title,
                     const std::vector<std::vector<double>>& Chi2_y,
                     const std::string& Chi2_y_title,
                     const std::vector<TDirectory*>& SampleDir,
                     const std::string Title);
 
  bool FullLLH;
  int NModelParams;
  bool Is_PriorPredictive;
 
  int TotalNumberOfSamples;
  std::vector<std::string> SampleNames;
  std::vector<int> SampleObjectMap;
 
  int Ntoys;
  std::vector<std::string> ParameterGroupsNotVaried;
  std::unordered_set<int> ParameterOnlyToVary;
 
  ParameterHandlerGeneric* ModelSystematic;
 
  std::vector<std::unique_ptr<TH1D>> Data_Hist;
  std::vector<std::unique_ptr<TH1D>> MC_Nom_Hist;
  std::vector<std::unique_ptr<TH1D>> W2_Nom_Hist;
 
  std::vector<std::vector<std::unique_ptr<TH1D>>> MC_Hist_Toy;
  std::vector<std::vector<std::unique_ptr<TH1D>>> W2_Hist_Toy;
 
  std::vector<double> ReweightWeight;
  std::vector<double> PenaltyTerm;
};