Ssth23Constraint.cxx
Go to the documentation of this file.
2 
3 #include "OscLib/IOscCalc.h"
4 
5 #include "Utilities/func/MathUtil.h"
6 
7 namespace ana
8 {
9  //----------------------------------------------------------------------
11  const SystShifts& /*syst*/) const
12  {
13 
14  const double ssth23 = util::sqr(sin(osc->GetTh23()));
15 
16  double chiA,chiB;
17  if(ssth23<fMeanA) chiA = util::sqr((ssth23-fMeanA)/fSigmaAMin);
18  else chiA = util::sqr((ssth23-fMeanA)/fSigmaAPlu);
19 
20  if(ssth23<fMeanB) chiB = util::sqr((ssth23-fMeanB)/fSigmaBMin);
21  else chiB = util::sqr((ssth23-fMeanB)/fSigmaBPlu);
22 
23  return std::min(chiA,chiB);
24  }
25 }
Cuts and Vars for the 2020 FD DiF Study.
Definition: vars.h:6
double ssth23
Simple record of shifts applied to systematic parameters.
Definition: SystShifts.h:20
virtual T GetTh23() const
Definition: IOscCalc.h:89
T sqr(T x)
More efficient square function than pow(x,2)
Definition: MathUtil.h:23
virtual double ChiSq(osc::IOscCalcAdjustable *osc, const SystShifts &shift=SystShifts::Nominal()) const override
Oscillation probability calculators.
Definition: Calcs.h:5
T sin(T number)
Definition: d0nt_math.hpp:132
T min(const caf::Proxy< T > &a, T b)