SolarConstraints.h

Go to the documentation of this file.

#pragma once

#include "CAFAna/Experiment/IExperiment.h"

#include <memory>

namespace ana
{
  /// \brief Constraints on the parameters \f$ \Delta m^2_{21} \f$ and
  /// \f$ \sin^22\theta_{12} \f$ from solar experiments
  class SolarConstraints: public IExperiment
  {
  public:
    SolarConstraints();

    SolarConstraints(const double dmsq,  const double errorDmsq,
                     const double ss2th, const double errorSs2th);

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

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

    double fCentralDmsq, fErrorDmsq;
    double fCentralAngle, fErrorAngle;
  };

  /// For use in CPT analysis
  class SolarConstraintsNu: public SolarConstraints
  {
  public:
    SolarConstraintsNu(const double dmsq,  const double errorDmsq,
                       const double ss2th, const double errorSs2th);

  protected:
    virtual double Dmsq21(osc::IOscCalcAdjustable* calc) const override;
    virtual double Th12(osc::IOscCalcAdjustable* calc) const override;
  };

  /// For use in CPT analysis
  class SolarConstraintsNubar: public SolarConstraints
  {
  public:
    SolarConstraintsNubar(const double dmsq,  const double errorDmsq,
                          const double ss2th, const double errorSs2th);

  protected:
    virtual double Dmsq21(osc::IOscCalcAdjustable* calc) const override;
    virtual double Th12(osc::IOscCalcAdjustable* calc) const override;
  };

  extern const SolarConstraints kSolarConstraintsPDG2017;
  extern const SolarConstraintsNu kSolarConstraintsNuFit2018nu;
  extern const SolarConstraintsNubar kSolarConstraintsNuFit2018bar;

}