Dmsq32Constraint.h

Go to the documentation of this file.

#pragma once

#include "CAFAna/Experiment/IExperiment.h"

namespace ana
{
  /// A simple Gaussian constraint on dmsq32, dependent on the hierarchy.
  class Dmsq32Constraint : public IExperiment
  {
  public:
    Dmsq32Constraint(double meanNH, double sigmaNH,
                     double meanIH, double sigmaIH)
      : fMeanNH(meanNH), fSigmaNH(sigmaNH),
        fMeanIH(fabs(meanIH)), fSigmaIH(sigmaIH)
    {
    }

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

  protected:
    double fMeanNH, fSigmaNH;
    double fMeanIH, fSigmaIH;
  };

  const Dmsq32Constraint kDmsq32ConstraintPDG2015(2.44e-3, 0.06e-3,  // NH
                                                  2.49e-3, 0.06e-3); // IH

  // http://pdg.lbl.gov/2017/tables/rpp2017-sum-leptons.pdf
  const Dmsq32Constraint kDmsq32ConstraintPDG2017(2.45e-3, 0.05e-3,  // NH
                                                  2.52e-3, 0.05e-3); // IH

  // http://pdg.lbl.gov/2019/tables/rpp2019-sum-leptons.pdf
  const Dmsq32Constraint kDmsq32ConstraintPDG2019(2.444e-3, 0.034e-3,  // NH
                                                  2.55e-3, 0.04e-3); // IH
}