OscProcessor.h

Go to the documentation of this file.

#pragma once
 
// MaCh3 includes
#include "Fitters/MCMCProcessor.h"
#include "Samples/HistogramUtils.h"
 
class OscProcessor : public MCMCProcessor {
  public:
    OscProcessor(const std::string &InputFile);
    virtual ~OscProcessor();
 
    void PerformJarlskogAnalysis();
 
    void MakePiePlot();
 
  protected:
    void LoadAdditionalInfo() override;
 
    void MakeJarlskogPlot(const std::unique_ptr<TH1D>& jarl,
                          const std::unique_ptr<TH1D>& jarl_flatsindcp,
                          const std::unique_ptr<TH1D>& jarl_NH,
                          const std::unique_ptr<TH1D>& jarl_NH_flatsindcp,
                          const std::unique_ptr<TH1D>& jarl_IH,
                          const std::unique_ptr<TH1D>& jarl_IH_flatsindcp);
 
    double CalcJarlskog(const double s2th13, const double s2th23, const double s2th12, const double dcp) const;
    double SamplePriorForParam(const int paramIndex,
                               const std::unique_ptr<TRandom3>& randGen,
                               const std::vector<double>& FlatBounds) const;
 
    bool PlotJarlskog;
 
    bool OscEnabled;
 
    std::string Sin2Theta13Name;
    std::string Sin2Theta12Name;
    std::string Sin2Theta23Name;
    std::string DeltaCPName;
    std::string DeltaM2_23Name;
 
    int Sin2Theta13Index;
    int Sin2Theta12Index;
    int Sin2Theta23Index;
    int DeltaCPIndex;
    int DeltaM2_23Index;
};