ReactorExperiment.h

Go to the documentation of this file.

#pragma once

#include "CAFAna/Experiment/IExperiment.h"

#include <memory>

namespace ana
{
  /// Very simple model allowing inclusion of reactor constraints
  class ReactorExperiment : public IExperiment
  {
  public:
    ReactorExperiment(double bestFit, double sigma)
      : fBestFit(bestFit), fSigma(sigma)
    {
    }

    virtual double ChiSq(osc::IOscCalcAdjustable* osc,
                         const SystShifts& shift = SystShifts::Nominal()) const override;

    void SaveTo(TDirectory* dir, const std::string& name) const override;
    static std::unique_ptr<ReactorExperiment> LoadFrom(TDirectory* dir, const std::string& name);
  protected:
    virtual double SSTh13(osc::IOscCalcAdjustable* osc) const;

    double fBestFit, fSigma;
  };

  /// For use in CPT analysis
  class ReactorExperimentNu: public ReactorExperiment
  {
  public:
    ReactorExperimentNu(double bestFit, double sigma)
      : ReactorExperiment(bestFit, sigma)
    {
    }

  protected:
    virtual double SSTh13(osc::IOscCalcAdjustable* osc) const override;
  };

  /// For use in CPT analysis
  class ReactorExperimentNubar: public ReactorExperiment
  {
  public:
    ReactorExperimentNubar(double bestFit, double sigma)
      : ReactorExperiment(bestFit, sigma)
    {
    }

  protected:
    virtual double SSTh13(osc::IOscCalcAdjustable* osc) const override;
  };


  /// A \ref ReactorExperiment initialized with the Nu2014 Daya Bay constraints
  const ReactorExperiment* DayaBayConstraint2014();

  /// Weighted average of all experiments as of first nue paper writing
  const ReactorExperiment* WorldReactorConstraint2015();

  /// Updated value for SecondAna based on the latest PDG
  const ReactorExperiment* WorldReactorConstraint2016();

  /// Reactor constraint from PDG 2017
  const ReactorExperiment* WorldReactorConstraint2017();

  /// Reactor constraint from NuFit 4.0 2018 CPT nu
  const ReactorExperiment* ReactorConstraintNuFit2018nu();

  /// Reactor constraint from NuFit 4.0 2018 CPT bar
  const ReactorExperiment* ReactorConstraintNuFit2018bar();

  /// Reactor constraint from PDG 2019
  const ReactorExperiment* WorldReactorConstraint2019();

}