osc::OscCalcSterile Class Reference

Adapt the PMNS_Sterile calculator to standard interface. More...

#include "/cvmfs/nova.opensciencegrid.org/externals/osclib/v00.12/src/OscLib/OscCalcSterile.h"

Inheritance diagram for osc::OscCalcSterile:

Public Member Functions
	OscCalcSterile ()
	OscCalcSterile (const OscCalcSterile &calc)
virtual	~OscCalcSterile ()
void	SetNFlavors (int nflavors)
virtual IOscCalcAdjustable *	Copy () const override
virtual double	P (int flavBefore, int flavAfter, double E) override
	E in GeV; flavors as PDG codes (so, neg==>antinu) More...
virtual void	SetL (double L) override
virtual void	SetRho (double rho) override
void	SetAngle (int i, int j, double th)
void	SetDelta (int i, int j, double delta)
void	SetDm (int i, double dm)
void	SetState (std::vector< double > state)
int	GetNFlavors () const
double	GetL () const override
double	GetRho () const override
double	GetDm (int i) const
double	GetAngle (int i, int j) const
double	GetDelta (int i, int j) const
std::vector< double >	GetState () const
virtual TMD5 *	GetParamsHash () const override
	Use to check two calculators are in the same state. More...
virtual void	SetDmsq21 (const T &dmsq21)=0
virtual void	SetDmsq32 (const T &dmsq32)=0
virtual void	SetTh12 (const T &th12)=0
virtual void	SetTh13 (const T &th13)=0
virtual void	SetTh23 (const T &th23)=0
virtual void	SetdCP (const T &dCP)=0
virtual T	GetDmsq21 () const
virtual T	GetDmsq32 () const
virtual T	GetTh12 () const
virtual T	GetTh13 () const
virtual T	GetTh23 () const
virtual T	GetdCP () const
virtual void	InvalidateCache ()
	Invalidate any caching used internally by the calculator. More...
virtual Eigen::Array< T, Eigen::Dynamic, 1 >	P (int flavBefore, int flavAfter, const std::vector< double > &E)
virtual Eigen::Array< T, Eigen::Dynamic, 1 >	P (int flavBefore, int flavAfter, const Eigen::ArrayXd &E)

Protected Member Functions
virtual void	SetDmsq21 (const double &dmsq21) override
virtual void	SetDmsq32 (const double &dmsq32) override
virtual void	SetTh12 (const double &th12) override
virtual void	SetTh13 (const double &th13) override
virtual void	SetTh23 (const double &th23) override
virtual void	SetdCP (const double &dCP) override
TMD5 *	GetParamsHashDefault (const std::string &txt) const
	This is only a safe implementation if your calculator only depends on these eight parameters. More...

Protected Attributes
PMNS_Sterile *	fPMNS_Sterile
int	fNFlavors
double	fRho
bool	fDirty
double	fPrevE
int	fPrevAnti
int	fPrevFlavBefore
double	fL
T	fDmsq21
T	fDmsq32
T	fTh12
T	fTh13
T	fTh23
T	fdCP

Detailed Description

Adapt the PMNS_Sterile calculator to standard interface.

Adapt the PMNS_Sterile calculator (3+N with matter effects) to standard interface

Definition at line 23 of file OscCalcSterile.h.

Constructor & Destructor Documentation

osc::OscCalcSterile::OscCalcSterile

(

)

Definition at line 10 of file OscCalcSterile.cxx.

References fNFlavors, and fPMNS_Sterile.

Referenced by Copy().

    : fPMNS_Sterile(0), fNFlavors(3), fDirty(true), fPrevE(0), fPrevAnti(0), fPrevFlavBefore(0)
  {
    fPMNS_Sterile = new PMNS_Sterile(fNFlavors);
  }

osc::OscCalcSterile::OscCalcSterile

(

const OscCalcSterile &

calc

)

Definition at line 51 of file OscCalcSterile.cxx.

References GetState(), SetState(), and makeDatasetsPage::state.

    : OscCalcSterile()
  {
    std::vector<double> state = calc.GetState();
    SetState(state);
  }

osc::OscCalcSterile::~OscCalcSterile ( )

virtual

Definition at line 71 of file OscCalcSterile.cxx.

References fPMNS_Sterile.

  {
    delete fPMNS_Sterile;
  }

Member Function Documentation

IOscCalcAdjustable * osc::OscCalcSterile::Copy ( ) const

overridevirtual

Implements osc::_IOscCalcAdjustable< T >.

Reimplemented in osc::OscCalcSterileTrivial, and osc::OscCalcSterileBeam.

Definition at line 77 of file OscCalcSterile.cxx.

References OscCalcSterile().

Referenced by osc::OscCalcSterileTrivial::~OscCalcSterileTrivial().

  {
    return new OscCalcSterile(*this);
  }

double osc::OscCalcSterile::GetAngle ( int  i, int  j  ) const

inline

Definition at line 52 of file OscCalcSterile.h.

References fPMNS_Sterile, and osc::PMNS_Sterile::GetAngle().

Referenced by ana::LikelihoodCovMxExperiment::ChiSq(), Draw2DSurface(), osc::OscCalcSterileBeam::GetParamsHash(), GetState(), ana::FitTheta13Sterile::GetValue(), ana::FitTheta23Sterile::GetValue(), ana::FitTheta14Sterile::GetValue(), ana::FitTheta24Sterile::GetValue(), ana::FitTheta34Sterile::GetValue(), ana::FitTheta13InDegreesSterile::GetValue(), ana::FitTheta23InDegreesSterile::GetValue(), ana::FitTheta14InDegreesSterile::GetValue(), ana::FitTheta24InDegreesSterile::GetValue(), ana::FitTheta34InDegreesSterile::GetValue(), PrintOscParams(), SetCalcParam(), and SystsOscParamAna().

{ return fPMNS_Sterile->GetAngle(i, j); }

template<typename T>

virtual T osc::_IOscCalcAdjustable< T >::GetdCP ( ) const

inlinevirtualinherited

Reimplemented in osc::OscCalcCPT.

Definition at line 95 of file IOscCalc.h.

Referenced by osc::CopyParams(), osc::OscCalcCPT::GetdCP(), getHists_FNEX(), joint_fit_2017_make_fc_slice(), joint_fit_2017_make_fc_surf(), jointsensitivity(), ana::ResetOscCalcToTemplate(), ana::ResetSterileCalcToDefault(), ana::SaveTo(), ana::FitSinSq2Theta13::SetValue(), test_ana(), and test_nue2018_fitter().

{ return fdCP    ; }

double osc::OscCalcSterile::GetDelta ( int  i, int  j  ) const

inline

Definition at line 53 of file OscCalcSterile.h.

References fPMNS_Sterile, osc::PMNS_Sterile::GetDelta(), GetParamsHash(), and GetState().

Referenced by ana::LikelihoodCovMxExperiment::ChiSq(), Draw2DSurface(), osc::OscCalcSterileBeam::GetParamsHash(), GetState(), ana::FitDelta13InPiUnitsSterile::GetValue(), ana::FitDelta14InPiUnitsSterile::GetValue(), ana::FitDelta24InPiUnitsSterile::GetValue(), PrintOscParams(), SetCalcParam(), and SystsOscParamAna().

{ return fPMNS_Sterile->GetDelta(i, j); }

double osc::OscCalcSterile::GetDm ( int  i ) const

inline

Definition at line 51 of file OscCalcSterile.h.

References fPMNS_Sterile, and osc::PMNS_Sterile::GetDm().

Referenced by ana::LikelihoodCovMxExperiment::ChiSq(), Draw2DSurface(), osc::OscCalcSterileBeam::GetParamsHash(), GetState(), ana::FitDmSq32Sterile::GetValue(), ana::FitDmSq41Sterile::GetValue(), ana::FitDmSq43Sterile::GetValue(), PrintOscParams(), SetCalcParam(), ana::FitDmSq32Sterile::SetValue(), ana::FitDmSq43Sterile::SetValue(), SystsOscParamAna(), and ana::ToEigenSterile().

{ return fPMNS_Sterile->GetDm(i); }

template<typename T>

virtual T osc::_IOscCalcAdjustable< T >::GetDmsq21 ( ) const

inlinevirtualinherited

Reimplemented in osc::OscCalcCPT.

Definition at line 90 of file IOscCalc.h.

Referenced by osc::CopyParams(), ana::SolarConstraints::Dmsq21(), osc::OscCalcCPT::GetDmsq21(), ana::FitDmSq21::GetValue(), ana::ResetOscCalcToTemplate(), ana::ResetSterileCalcToDefault(), and ana::SaveTo().

{ return fDmsq21 ; }

template<typename T>

virtual T osc::_IOscCalcAdjustable< T >::GetDmsq32 ( ) const

inlinevirtualinherited

Reimplemented in osc::OscCalcCPT.

Definition at line 91 of file IOscCalc.h.

Referenced by CAF_makeCAFSensitivities_for_FNEX(), cc(), ana::Dmsq32Constraint::ChiSq(), ana::AtmConstraint::ChiSq(), osc::CopyParams(), demoFitContours(), demoFitSlices(), getBestFit(), osc::OscCalcCPT::GetDmsq32(), getHists_FNEX(), ana::FitDmSq32::GetValue(), ana::FitDmSq32Scaled::GetValue(), ana::FitDmSq32NormalHierarchy::GetValue(), ana::FitDmSq32ScaledNH::GetValue(), ana::FitDmSq32ScaledIH::GetValue(), ana::FitDmSq32InvertedHierarchy::GetValue(), joint_fit_2017_contours(), joint_fit_2017_slices(), joint_fit_2018_contours(), joint_fit_2018_slices(), joint_fit_2019_contours(), joint_fit_2019_slices(), joint_fit_future_contour_univ(), jointsensitivity(), MakeCAFSensitivities_for_FNEX(), plot_3NDvsFD(), plot_NDvsFD_REW(), plot_NDvsFD_weights(), plot_nueFD_Signal_REWvsNOM(), ana::ResetOscCalcToTemplate(), ana::ResetSterileCalcToDefault(), run_joint_fit_2020_contours(), run_joint_fit_2020_slices(), ana::SaveTo(), sensitivity2018(), sensitivity2020(), test_nue2018_fitter(), test_nueextrapsyst(), Tutorial2019FitContours(), and Tutorial2019FitSlices().

{ return fDmsq32 ; }

double osc::OscCalcSterile::GetL ( ) const

inlineoverridevirtual

Reimplemented from osc::_IOscCalcAdjustable< T >.

Definition at line 49 of file OscCalcSterile.h.

References osc::_IOscCalcAdjustable< T >::fL.

Referenced by ana::ToEigenSterile().

{ return fL; }

int osc::OscCalcSterile::GetNFlavors ( ) const

inline

Definition at line 48 of file OscCalcSterile.h.

References fPMNS_Sterile, and osc::PMNS_Sterile::GetNFlavors().

{ return fPMNS_Sterile->GetNFlavors(); }

TMD5 * osc::OscCalcSterile::GetParamsHash ( ) const

overridevirtual

Use to check two calculators are in the same state.

Returns

Null means not implemented for this calculator

Reimplemented from osc::_IOscCalc< T >.

Reimplemented in osc::OscCalcSterileBeam.

Definition at line 59 of file OscCalcSterile.cxx.

References GetState(), runNovaSAM::ret, and string.

Referenced by GetDelta().

  {
    TMD5* ret = new TMD5;
    std::string txt = "PMNSSterile";
    ret->Update((unsigned char*)txt.c_str(), txt.size());
    std::vector<double> buf = GetState();
    ret->Update((unsigned char*)&buf[0], sizeof(double)*buf.size());
    ret->Final();
    return ret;
  }

template<typename T >

TMD5 * osc::_IOscCalcAdjustable< T >::GetParamsHashDefault ( const std::string &  txt ) const

protectedinherited

This is only a safe implementation if your calculator only depends on these eight parameters.

Parameters

txt	A string to uniquely identify your calculator class

Definition at line 47 of file IOscCalc.cxx.

References runNovaSAM::ret, and T.

Referenced by osc::_OscCalcPMNS< T >::GetParamsHash(), osc::_OscCalcPMNSOpt< T >::GetParamsHash(), and osc::analytic::_OscCalc< T >::GetParamsHash().

  {
    TMD5* ret = new TMD5;
    ret->Update((unsigned char*)txt.c_str(), txt.size());
    const int kNumParams = 8;
    T buf[kNumParams] = {fRho, fL, fDmsq21, fDmsq32,
                         fTh12, fTh13, fTh23, fdCP};
    ret->Update((unsigned char*)buf, sizeof(T)*kNumParams);
    ret->Final();
    return ret;
  }

double osc::OscCalcSterile::GetRho ( ) const

inlineoverridevirtual

Reimplemented from osc::_IOscCalcAdjustable< T >.

Definition at line 50 of file OscCalcSterile.h.

References fRho.

{ return fRho; }

std::vector< double > osc::OscCalcSterile::GetState

(

)

const

Definition at line 17 of file OscCalcSterile.cxx.

References osc::_IOscCalcAdjustable< T >::fL, fNFlavors, fRho, GetAngle(), GetDelta(), GetDm(), MECModelEnuComparisons::i, calib::j, and makeDatasetsPage::state.

Referenced by GetDelta(), GetParamsHash(), OscCalcSterile(), osc::OscCalcSterileTrivial::OscCalcSterileTrivial(), and ana::SaveTo().

  {
    std::vector<double> state;
    state.push_back((double)fNFlavors);
    state.push_back(fL);
    state.push_back(fRho);
    for(int i = 2; i <= fNFlavors; ++i) state.push_back(GetDm(i));
    for(int j = 2; j <= fNFlavors; ++j) {
      for(int i = 1; i < j; ++i) {
        state.push_back(GetAngle(i, j));
        if(i+1 != j) state.push_back(GetDelta(i, j));
      }
    }
    return state;
  }

template<typename T>

virtual T osc::_IOscCalcAdjustable< T >::GetTh12 ( ) const

inlinevirtualinherited

Reimplemented in osc::OscCalcCPT.

Definition at line 92 of file IOscCalc.h.

Referenced by osc::CopyParams(), osc::OscCalcCPT::GetTh12(), ana::FitTanSqTheta12::GetValue(), ana::FitSinSq2Theta12::GetValue(), ana::ResetOscCalcToTemplate(), ana::ResetSterileCalcToDefault(), ana::SaveTo(), and ana::SolarConstraints::Th12().

{ return fTh12   ; }

template<typename T>

virtual T osc::_IOscCalcAdjustable< T >::GetTh13 ( ) const

inlinevirtualinherited

Reimplemented in osc::OscCalcCPT.

Definition at line 93 of file IOscCalc.h.

Referenced by osc::CopyParams(), demoFitContours(), osc::OscCalcCPT::GetTh13(), ana::FitTheta13::GetValue(), ana::FitSinSq2Theta13::GetValue(), joint_fit_2018_contours(), joint_fit_future_bestfit_univ(), monoprob(), ana::ResetOscCalcToTemplate(), ana::ResetSterileCalcToDefault(), ana::SaveTo(), sensitivity2020(), test_ana(), and Tutorial2019FitContours().

{ return fTh13   ; }

template<typename T>

virtual T osc::_IOscCalcAdjustable< T >::GetTh23 ( ) const

inlinevirtualinherited

Reimplemented in osc::OscCalcCPT.

Definition at line 94 of file IOscCalc.h.

Referenced by cc(), ana::AtmConstraint::ChiSq(), ana::Ssth23Constraint::ChiSq(), osc::CopyParams(), getBestFit(), osc::OscCalcCPT::GetTh23(), ana::FitSinSqTheta23::GetValue(), FitSinSqTheta23UO::GetValue(), ana::FitSinSqTheta23UpperOctant::GetValue(), FitSinSqTheta23LO::GetValue(), ana::FitSinSqTheta23LowerOctant::GetValue(), ana::FitSinSq2Theta23::GetValue(), FitSinSqTheta23MaxMix::GetValue(), FitSinSqTheta23LowerOctant::GetValue(), FitSinSqTheta23UpperOctant::GetValue(), jointsensitivity(), ana::ResetOscCalcToTemplate(), ana::ResetSterileCalcToDefault(), ana::SaveTo(), and test_nue2018_fitter().

{ return fTh23   ; }

template<typename T>

virtual void osc::_IOscCalcAdjustable< T >::InvalidateCache ( )

inlinevirtualinherited

Invalidate any caching used internally by the calculator.

Some calculators use a cache that can become stale in ways that the calculator may not know about (e.g., Stan var clearing). Default implementation does nothing.

Reimplemented in osc::_OscCalcDMP< T >.

Definition at line 102 of file IOscCalc.h.

{};

template<class T >

Eigen::Array< T, Eigen::Dynamic, 1 > osc::_IOscCalc< T >::P ( int  flavBefore, int  flavAfter, const std::vector< double > &  E  )

virtualinherited

Default implementation forwards to non-vector version using a simple loop. Override if your calculator has a more efficient implementation.

Reimplemented in osc::analytic::_OscCalc< T >, osc::OscCalcPMNSOptEigen, and osc::_OscCalcDMP< T >.

Definition at line 25 of file IOscCalc.cxx.

References MECModelEnuComparisons::i, osc::_IOscCalc< T >::P(), P, runNovaSAM::ret, and submit_hadd::u.

  {
    Eigen::Matrix<T,Eigen::Dynamic,1> ret(E.size());
    for(auto i = 0u; i < E.size(); i++) {
      ret(i) = this->P(flavBefore, flavAfter, E[i]);
    }
    return ret.array().isNaN().select(0, ret);
  }

template<class T >

Eigen::Array< T, Eigen::Dynamic, 1 > osc::_IOscCalc< T >::P ( int  flavBefore, int  flavAfter, const Eigen::ArrayXd &  E  )

virtualinherited

Default implementation forawrds to vector<double> version. Override if your calculator has a more efficient implementation.

Reimplemented in osc::analytic::_OscCalc< T >.

Definition at line 35 of file IOscCalc.cxx.

References P.

  {
    return P(flavBefore, flavAfter, std::vector<double>(&E[0], &E[0]+E.size()));
  }

double osc::OscCalcSterile::P ( int  flavBefore, int  flavAfter, double  E  )

overridevirtual

E in GeV; flavors as PDG codes (so, neg==>antinu)

Implements osc::_IOscCalc< T >.

Reimplemented in osc::OscCalcSterileTrivial.

Definition at line 155 of file OscCalcSterile.cxx.

References abs(), ana::assert(), E, e, fDirty, osc::_IOscCalcAdjustable< T >::fL, fPMNS_Sterile, fPrevAnti, fPrevE, fPrevFlavBefore, fRho, MECModelEnuComparisons::i, calib::j, osc::PMNS_Sterile::P(), osc::PMNS_Sterile::PropMatter(), and osc::PMNS_Sterile::ResetToFlavour().

Referenced by LOverEPlot(), ana::ToEigenSterile(), and osc::OscCalcSterileTrivial::~OscCalcSterileTrivial().

  {
    const int anti = (flavBefore > 0) ? +1 : -1;
    //anti must be +/- 1 but flavAfter can be zero
    assert(flavAfter/anti >= 0);
    if (anti != fPrevAnti) fDirty = true;
    fPrevAnti = anti;

    if (flavBefore != fPrevFlavBefore) fDirty = true;
    fPrevFlavBefore = flavBefore;

    if (abs(fPrevE - E) > 1e-8) fDirty = true;
    fPrevE = E;
    
    int i = -1, j = -1;
    if(abs(flavBefore) == 12) i = 0;
    if(abs(flavBefore) == 14) i = 1;
    if(abs(flavBefore) == 16) i = 2;
    if(abs(flavAfter) == 12) j = 0;
    if(abs(flavAfter) == 14) j = 1;
    if(abs(flavAfter) == 16) j = 2;
    if(abs(flavAfter) == 0)  j = 3;
    assert(i >= 0 && j >= 0);

    if (fDirty) {
      fPMNS_Sterile->ResetToFlavour(i);      
      // Assume Z/A=0.5
      const double Ne = fRho/2;
      fPMNS_Sterile->PropMatter(fL, E, Ne, anti);
    }
    
    // Return the active fraction
    if (j == 3) return fPMNS_Sterile->P(0) + fPMNS_Sterile->P(1) + fPMNS_Sterile->P(2);
    else        return fPMNS_Sterile->P(j);
  }

void osc::OscCalcSterile::SetAngle	(	int	i,
		int	j,
		double	th
	)

Definition at line 92 of file OscCalcSterile.cxx.

References fDirty, fPMNS_Sterile, and osc::PMNS_Sterile::SetAngle().

Referenced by BlessedPlotsAna(), BlessedPlotsAnaByPeriod(), FitSystEffectsAna(), LOverEPlot(), PlotNus17Prediction(), ana::PredictionSterile::PredictUnoscillated(), ResetAngles(), ana::ResetSterileCalcToDefault(), SetAngles(), SetCalcParam(), SetFakeSignalParams(), SetNus20Params(), SetRho(), SetState(), ana::FitTheta13Sterile::SetValue(), ana::FitSinSqTheta13Sterile::SetValue(), ana::FitSinSq2Theta13Sterile::SetValue(), ana::FitTheta23Sterile::SetValue(), ana::FitSinSqTheta23Sterile::SetValue(), ana::FitTheta14Sterile::SetValue(), ana::FitSinSqTheta14Sterile::SetValue(), ana::FitTheta24Sterile::SetValue(), ana::FitSinSqTheta24Sterile::SetValue(), ana::FitTheta34Sterile::SetValue(), ana::FitSinSqTheta34Sterile::SetValue(), ana::FitTheta13InDegreesSterile::SetValue(), ana::FitTheta23InDegreesSterile::SetValue(), ana::FitTheta14InDegreesSterile::SetValue(), ana::FitTheta24InDegreesSterile::SetValue(), ana::FitTheta34InDegreesSterile::SetValue(), sterile_demo(), SterileOscillator(), and SystsOscParamAna().

  {
    fDirty = true; 
    fPMNS_Sterile->SetAngle(i, j, th);
  }

void osc::OscCalcSterile::SetdCP ( const double &  dCP )

overrideprotectedvirtual

Definition at line 148 of file OscCalcSterile.cxx.

References ana::assert(), om::cerr, and allTimeWatchdog::endl.

  {
    std::cerr << "Must use SetDelta!" << std::endl;
    assert(false);
  }

template<typename T>

virtual void osc::_IOscCalcAdjustable< T >::SetdCP ( const T &  dCP )

pure virtualinherited

Implemented in osc::analytic::_OscCalc< T >, osc::OscCalcPMNSOptEigen, osc::_OscCalcDMP< T >, osc::OscCalc, osc::_OscCalcPMNSOpt< T >, osc::OscCalcGeneral, osc::_OscCalcPMNS< T >, and osc::OscCalcPMNS_NSI.

Referenced by BicountEllipse_dCP(), BiprobEllipse_dCP(), CalcRWithSystsNus17(), ComparePredictions(), osc::CopyParams(), default_calc(), genie_contours(), getHists_FNEX(), joint_fit_2017_make_fc_slice(), joint_fit_2017_make_fc_surf(), joint_fit_2018_slices(), joint_fit_2019_bestfit(), joint_fit_2019_slices(), joint_fit_future_bestfit_univ(), jointsensitivity(), ana::LoadFrom< osc::IOscCalcAdjustable >(), make_extrap_figure(), make_extrap_figure_hists(), monoprob(), plot_datamcpred(), plot_nue_xsec_pred(), Plotting_DataAndPrediction(), PredRatioToNom(), resetCalc(), mcmc_ana::ResetCalculator(), test::ResetCalculator(), mcmc::ResetCalculator(), ana::ResetOscCalcToDefault(), ResetOscCalcToMyDefault(), ana::ResetOscCalcToTemplate(), RestartCalculator(), run_joint_fit_2020_bestfit(), sensitivity2018(), sensitivity2020(), osc::OscCalcCPT::SetdCP(), SetOscillationParameters(), cmf::ShifterAndWeighter::SetOscillationVals(), ana::FitDeltaInPiUnits::SetValue(), systematics_extrap_comp_from_pred_interp(), systematics_summary_from_pred_interp(), and TheoryEllipse_dCP().

void osc::OscCalcSterile::SetDelta	(	int	i,
		int	j,
		double	delta
	)

Definition at line 99 of file OscCalcSterile.cxx.

References fDirty, fPMNS_Sterile, and osc::PMNS_Sterile::SetDelta().

Referenced by MakeSurface(), ana::PredictionSterile::PredictUnoscillated(), ResetAngles(), ana::ResetSterileCalcToDefault(), SetAngles(), SetCalcParam(), SetNus20Params(), SetRho(), SetState(), ana::FitDelta13InPiUnitsSterile::SetValue(), ana::FitDelta14InPiUnitsSterile::SetValue(), ana::FitDelta24InPiUnitsSterile::SetValue(), sterile_demo(), SterileOscillator(), and SystsOscParamAna().

  {
    fDirty = true; 
    fPMNS_Sterile->SetDelta(i, j, delta);
  }

void osc::OscCalcSterile::SetDm	(	int	i,
		double	dm
	)

Definition at line 106 of file OscCalcSterile.cxx.

References fDirty, fPMNS_Sterile, and osc::PMNS_Sterile::SetDm().

Referenced by BlessedPlotsAna(), BlessedPlotsAnaByPeriod(), FitSystEffectsAna(), LOverEPlot(), MakeSurface(), ana::PredictionSterile::PredictUnoscillated(), ResetAngles(), ana::ResetSterileCalcToDefault(), SetAngles(), SetCalcParam(), SetFakeSignalParams(), SetNus20Params(), cmf::ShifterAndWeighter::SetOscillationVals(), SetRho(), SetState(), ana::FitDmSq32Sterile::SetValue(), ana::FitDmSq41Sterile::SetValue(), ana::FitDmSq43Sterile::SetValue(), sterile_demo(), SterileOscillator(), and SystsOscParamAna().

  {
    fDirty = true; 
    fPMNS_Sterile->SetDm(i, dm);
  }

void osc::OscCalcSterile::SetDmsq21 ( const double &  dmsq21 )

overrideprotectedvirtual

Definition at line 113 of file OscCalcSterile.cxx.

References ana::assert(), om::cerr, and allTimeWatchdog::endl.

  {
    std::cerr << "Must use SetDm!" << std::endl;
    assert(false);
  }

template<typename T>

virtual void osc::_IOscCalcAdjustable< T >::SetDmsq21 ( const T &  dmsq21 )

pure virtualinherited

Implemented in osc::analytic::_OscCalc< T >, osc::OscCalcPMNSOptEigen, osc::_OscCalcDMP< T >, osc::OscCalc, osc::_OscCalcPMNSOpt< T >, osc::_OscCalcPMNS< T >, osc::OscCalcPMNS_NSI, and osc::OscCalcGeneral.

Referenced by osc::CopyParams(), default_calc(), ana::LoadFrom< osc::IOscCalcAdjustable >(), resetCalc(), mcmc_ana::ResetCalculator(), test::ResetCalculator(), mcmc::ResetCalculator(), ana::ResetOscCalcToDefault(), ResetOscCalcToMyDefault(), ana::ResetOscCalcToTemplate(), RestartCalculator(), osc::OscCalcCPT::SetDmsq21(), cmf::ShifterAndWeighter::SetOscillationVals(), and ana::FitDmSq21::SetValue().

void osc::OscCalcSterile::SetDmsq32 ( const double &  dmsq32 )

overrideprotectedvirtual

Definition at line 120 of file OscCalcSterile.cxx.

References ana::assert(), om::cerr, and allTimeWatchdog::endl.

  {
    std::cerr << "Must use SetDm!" << std::endl;
    assert(false);
  }

template<typename T>

virtual void osc::_IOscCalcAdjustable< T >::SetDmsq32 ( const T &  dmsq32 )

pure virtualinherited

Implemented in osc::analytic::_OscCalc< T >, osc::OscCalcPMNSOptEigen, osc::_OscCalcDMP< T >, osc::OscCalc, osc::_OscCalcPMNSOpt< T >, osc::_OscCalcPMNS< T >, osc::OscCalcPMNS_NSI, and osc::OscCalcGeneral.

Referenced by CAF_makeCAFSensitivities_for_FNEX(), CalcRWithSystsNus17(), ComparePredictions(), osc::CopyParams(), default_calc(), demoFitContours(), demoFitSlices(), demoPlotSystBands(), genie_contours(), getContProf(), getContProf_Sensitivity(), getHists_FNEX(), joint_fit_2017_contours(), joint_fit_2017_slices(), joint_fit_2018_contours(), joint_fit_2018_slices(), joint_fit_2019_bestfit(), joint_fit_2019_contours(), joint_fit_2019_slices(), joint_fit_future_contour_univ(), jointsensitivity(), ana::LoadFrom< osc::IOscCalcAdjustable >(), make_extrap_figure(), make_extrap_figure_hists(), make_prediction_extrap(), make_prediction_noextrap(), make_prediction_tables(), make_surfprof(), make_surfprof_sensitivity(), MakeCAFSensitivities_for_FNEX(), makePredTables(), monoprob(), plot_datamcpred(), plot_nue_xsec_pred(), plot_oscprob(), plot_predictions(), plotDataPred(), plotDataPred_SplitCanvas(), Plotting_DataAndPrediction(), PredRatioToNom(), resetCalc(), ResetCalc(), mcmc_ana::ResetCalculator(), test::ResetCalculator(), mcmc::ResetCalculator(), ana::ResetOscCalcToDefault(), ana::ResetOscCalcToDefaultIH(), ResetOscCalcToMyDefault(), ana::ResetOscCalcToTemplate(), RestartCalculator(), run_joint_fit_2020_bestfit(), run_joint_fit_2020_contours(), run_joint_fit_2020_slices(), saveFDMCHists(), sensitivity2018(), sensitivity2020(), osc::OscCalcCPT::SetDmsq32(), SetHierarchy(), SetOscillationParameters(), cmf::ShifterAndWeighter::SetOscillationVals(), ana::FitDmSq32::SetValue(), ana::FitDmSq32Scaled::SetValue(), ana::FitDmSq32NormalHierarchy::SetValue(), ana::FitDmSq32ScaledNH::SetValue(), ana::FitDmSq32ScaledIH::SetValue(), ana::FitDmSq32InvertedHierarchy::SetValue(), systematics_extrap_comp_from_pred_interp(), systematics_summary_from_pred_interp(), Tutorial2019FitContours(), and Tutorial2019FitSlices().

virtual void osc::OscCalcSterile::SetL ( double  L )

inlineoverridevirtual

Implements osc::_IOscCalcAdjustable< T >.

Definition at line 38 of file OscCalcSterile.h.

References fDirty, osc::_IOscCalcAdjustable< T >::fL, and CLHEP::L.

Referenced by LOverEPlot(), ana::PredictionSterile::PredictUnoscillated(), ana::ResetSterileCalcToDefault(), SetState(), sterile_demo(), SterileOscillator(), and ana::ToEigenSterile().

{fDirty = true; fL   = L;}

void osc::OscCalcSterile::SetNFlavors

(

int

nflavors

)

Definition at line 83 of file OscCalcSterile.cxx.

References fDirty, fNFlavors, and fPMNS_Sterile.

Referenced by BlessedPlotsAna(), BlessedPlotsAnaByPeriod(), ana::DefaultSterileCalc(), Draw2DSurface(), FitParamEffectsAna(), FitSystEffectsAna(), make_fc_nus_surfs_nersc_2018(), make_nus17_fc_surfs(), MakeExtrapSurface(), PlotNus17Prediction(), PlotNus17PredSystsData(), ana::PredictionSterile::PredictUnoscillated(), SetState(), SpectrumParamEffectsAna(), sterile_demo(), and SystsOscParamAna().

  {
    fDirty = true;
    delete fPMNS_Sterile;
    fNFlavors = nflavors;
    fPMNS_Sterile = new PMNS_Sterile(fNFlavors);
  }

virtual void osc::OscCalcSterile::SetRho ( double  rho )

inlineoverridevirtual

Implements osc::_IOscCalcAdjustable< T >.

Definition at line 39 of file OscCalcSterile.h.

References delta, fDirty, fRho, MECModelEnuComparisons::i, calib::j, SetAngle(), SetDelta(), SetDm(), SetState(), and makeDatasetsPage::state.

Referenced by ana::PredictionSterile::PredictUnoscillated(), ana::ResetSterileCalcToDefault(), SetState(), and sterile_demo().

{fDirty = true; fRho = rho;}

void osc::OscCalcSterile::SetState

(

std::vector< double >

state

)

Definition at line 34 of file OscCalcSterile.cxx.

References fDirty, fNFlavors, MECModelEnuComparisons::i, calib::j, SetAngle(), SetDelta(), SetDm(), SetL(), SetNFlavors(), and SetRho().

Referenced by OscCalcSterile(), osc::OscCalcSterileTrivial::OscCalcSterileTrivial(), and SetRho().

  {
    int iState(0);
    fDirty = true ;
    SetNFlavors((int)state[iState++]);
    SetL(state[iState++]);
    SetRho(state[iState++]);
    for (int i = 2; i <= fNFlavors; ++i) SetDm(i, state[iState++]);
    for(int j = 2; j <= fNFlavors; ++j) {
      for(int i = 1; i < j; ++i) {
        SetAngle(i, j, state[iState++]);
        if(i+1 != j) SetDelta(i, j, state[iState++]);
      }
    }
  }

void osc::OscCalcSterile::SetTh12 ( const double &  th12 )

overrideprotectedvirtual

Definition at line 127 of file OscCalcSterile.cxx.

References ana::assert(), om::cerr, and allTimeWatchdog::endl.

  {
    std::cerr << "Must use SetAngle!" << std::endl;
    assert(false);
  }

template<typename T>

virtual void osc::_IOscCalcAdjustable< T >::SetTh12 ( const T &  th12 )

pure virtualinherited

Implemented in osc::analytic::_OscCalc< T >, osc::OscCalcPMNSOptEigen, osc::_OscCalcDMP< T >, osc::OscCalc, osc::_OscCalcPMNSOpt< T >, osc::OscCalcGeneral, osc::_OscCalcPMNS< T >, and osc::OscCalcPMNS_NSI.

Referenced by osc::CopyParams(), default_calc(), ana::LoadFrom< osc::IOscCalcAdjustable >(), resetCalc(), mcmc_ana::ResetCalculator(), test::ResetCalculator(), mcmc::ResetCalculator(), ana::ResetOscCalcToDefault(), ResetOscCalcToMyDefault(), ana::ResetOscCalcToTemplate(), RestartCalculator(), cmf::ShifterAndWeighter::SetOscillationVals(), osc::OscCalcCPT::SetTh12(), ana::FitTanSqTheta12::SetValue(), and ana::FitSinSq2Theta12::SetValue().

void osc::OscCalcSterile::SetTh13 ( const double &  th13 )

overrideprotectedvirtual

Definition at line 134 of file OscCalcSterile.cxx.

References ana::assert(), om::cerr, and allTimeWatchdog::endl.

  {
    std::cerr << "Must use SetAngle!" << std::endl;
    assert(false);
  }

template<typename T>

virtual void osc::_IOscCalcAdjustable< T >::SetTh13 ( const T &  th13 )

pure virtualinherited

Implemented in osc::analytic::_OscCalc< T >, osc::OscCalcPMNSOptEigen, osc::_OscCalcDMP< T >, osc::OscCalc, osc::_OscCalcPMNSOpt< T >, osc::OscCalcGeneral, osc::_OscCalcPMNS< T >, and osc::OscCalcPMNS_NSI.

Referenced by ComparePredictions(), osc::CopyParams(), default_calc(), joint_fit_2018_slices(), joint_fit_2019_bestfit(), joint_fit_2019_slices(), ana::LoadFrom< osc::IOscCalcAdjustable >(), plot_nue_xsec_pred(), PredRatioToNom(), resetCalc(), mcmc_ana::ResetCalculator(), test::ResetCalculator(), mcmc::ResetCalculator(), ana::ResetOscCalcToDefault(), ResetOscCalcToMyDefault(), ana::ResetOscCalcToTemplate(), RestartCalculator(), run_joint_fit_2020_bestfit(), SetOscillationParameters(), cmf::ShifterAndWeighter::SetOscillationVals(), osc::OscCalcCPT::SetTh13(), ana::FitTheta13::SetValue(), ana::FitSinSq2Theta13::SetValue(), and systematics_summary_from_pred_interp().

void osc::OscCalcSterile::SetTh23 ( const double &  th23 )

overrideprotectedvirtual

Definition at line 141 of file OscCalcSterile.cxx.

References ana::assert(), om::cerr, and allTimeWatchdog::endl.

  {
    std::cerr << "Must use SetAngle!" << std::endl;
    assert(false);
  }

template<typename T>

virtual void osc::_IOscCalcAdjustable< T >::SetTh23 ( const T &  th23 )

pure virtualinherited

Implemented in osc::analytic::_OscCalc< T >, osc::OscCalcPMNSOptEigen, osc::_OscCalcDMP< T >, osc::OscCalc, osc::_OscCalcPMNSOpt< T >, osc::OscCalcGeneral, osc::_OscCalcPMNS< T >, and osc::OscCalcPMNS_NSI.

Referenced by CalcRWithSystsNus17(), ComparePredictions(), osc::CopyParams(), default_calc(), demoPlotSystBands(), genie_contours(), getContProf(), getContProf_Sensitivity(), joint_fit_2018_slices(), joint_fit_2019_bestfit(), joint_fit_2019_slices(), jointsensitivity(), ana::LoadFrom< osc::IOscCalcAdjustable >(), make_extrap_figure(), make_extrap_figure_hists(), make_prediction_extrap(), make_prediction_noextrap(), make_prediction_tables(), make_surfprof(), make_surfprof_sensitivity(), makePredTables(), monoprob(), plot_datamcpred(), plot_nue_xsec_pred(), plot_oscprob(), plot_predictions(), plotDataPred(), plotDataPred_SplitCanvas(), Plotting_DataAndPrediction(), PredRatioToNom(), resetCalc(), ResetCalc(), mcmc_ana::ResetCalculator(), test::ResetCalculator(), mcmc::ResetCalculator(), ana::ResetOscCalcToDefault(), ResetOscCalcToMyDefault(), ana::ResetOscCalcToTemplate(), RestartCalculator(), run_joint_fit_2020_bestfit(), saveFDMCHists(), sensitivity2018(), sensitivity2020(), SetOscillationParameters(), cmf::ShifterAndWeighter::SetOscillationVals(), osc::OscCalcCPT::SetTh23(), ana::FitSinSqTheta23::SetValue(), FitSinSqTheta23UO::SetValue(), ana::FitSinSqTheta23UpperOctant::SetValue(), FitSinSqTheta23LO::SetValue(), ana::FitSinSqTheta23LowerOctant::SetValue(), ana::FitSinSq2Theta23::SetValue(), FitSinSqTheta23MaxMix::SetValue(), FitSinSqTheta23LowerOctant::SetValue(), FitSinSqTheta23UpperOctant::SetValue(), systematics_extrap_comp_from_pred_interp(), and systematics_summary_from_pred_interp().

Member Data Documentation

template<typename T>

T osc::_IOscCalcAdjustable< T >::fdCP

protectedinherited

Definition at line 119 of file IOscCalc.h.

Referenced by osc::analytic::_OscCalc< T >::_P(), osc::OscCalcPMNS_CPT::P(), osc::_OscCalcPMNS< T >::P(), osc::_OscCalcPMNSOpt< T >::P(), osc::OscCalcPMNS_CPT::SetdCP(), osc::_OscCalcPMNS< T >::SetdCP(), osc::_OscCalcPMNSOpt< T >::SetdCP(), and osc::analytic::_OscCalc< T >::SetdCP().

bool osc::OscCalcSterile::fDirty

protected

Definition at line 69 of file OscCalcSterile.h.

Referenced by P(), SetAngle(), SetDelta(), SetDm(), SetL(), SetNFlavors(), SetRho(), and SetState().

template<typename T>

T osc::_IOscCalcAdjustable< T >::fDmsq21

protectedinherited

Definition at line 114 of file IOscCalc.h.

Referenced by osc::OscCalcPMNS_CPT::P(), osc::_OscCalcPMNS< T >::P(), osc::_OscCalcPMNSOpt< T >::P(), osc::OscCalcPMNS_CPT::SetDmsq21(), osc::_OscCalcPMNS< T >::SetDmsq21(), osc::_OscCalcPMNSOpt< T >::SetDmsq21(), osc::analytic::_OscCalc< T >::SetDmsq21(), and osc::analytic::SolveCubic().

template<typename T>

T osc::_IOscCalcAdjustable< T >::fDmsq32

protectedinherited

Definition at line 115 of file IOscCalc.h.

Referenced by osc::OscCalcPMNS_CPT::P(), osc::_OscCalcPMNS< T >::P(), osc::_OscCalcPMNSOpt< T >::P(), osc::OscCalcPMNS_CPT::SetDmsq32(), osc::_OscCalcPMNS< T >::SetDmsq32(), osc::_OscCalcPMNSOpt< T >::SetDmsq32(), osc::analytic::_OscCalc< T >::SetDmsq32(), and osc::analytic::SolveCubic().

template<typename T>

double osc::_IOscCalcAdjustable< T >::fL

protectedinherited

Definition at line 113 of file IOscCalc.h.

Referenced by osc::analytic::_OscCalc< T >::_P(), GetL(), osc::OscCalcSterileBeam::GetParamsHash(), osc::OscCalcPMNS_CPT::GetParamsHashDefaultBar(), GetState(), osc::OscCalcPMNS_CPT::P(), osc::_OscCalcPMNS< T >::P(), osc::_OscCalcPMNSOpt< T >::P(), P(), osc::OscCalcPMNS_CPT::SetL(), osc::_OscCalcPMNS< T >::SetL(), osc::_OscCalcPMNSOpt< T >::SetL(), SetL(), and osc::analytic::_OscCalc< T >::SetL().

int osc::OscCalcSterile::fNFlavors

protected

Definition at line 67 of file OscCalcSterile.h.

Referenced by osc::OscCalcSterileBeam::GetParamsHash(), GetState(), OscCalcSterile(), SetNFlavors(), and SetState().

PMNS_Sterile* osc::OscCalcSterile::fPMNS_Sterile

protected

Definition at line 58 of file OscCalcSterile.h.

Referenced by GetAngle(), GetDelta(), GetDm(), GetNFlavors(), OscCalcSterile(), P(), SetAngle(), SetDelta(), SetDm(), SetNFlavors(), and ~OscCalcSterile().

int osc::OscCalcSterile::fPrevAnti

protected

Definition at line 71 of file OscCalcSterile.h.

Referenced by P().

double osc::OscCalcSterile::fPrevE

protected

Definition at line 70 of file OscCalcSterile.h.

Referenced by P().

int osc::OscCalcSterile::fPrevFlavBefore

protected

Definition at line 72 of file OscCalcSterile.h.

Referenced by P().

double osc::OscCalcSterile::fRho

protected

Definition at line 68 of file OscCalcSterile.h.

Referenced by osc::OscCalcSterileBeam::GetParamsHash(), GetRho(), GetState(), P(), and SetRho().

template<typename T>

T osc::_IOscCalcAdjustable< T >::fTh12

protectedinherited

Definition at line 116 of file IOscCalc.h.

Referenced by osc::analytic::_OscCalc< T >::_P(), osc::OscCalcPMNS_CPT::P(), osc::_OscCalcPMNS< T >::P(), osc::_OscCalcPMNSOpt< T >::P(), osc::OscCalcPMNS_CPT::SetTh12(), osc::_OscCalcPMNS< T >::SetTh12(), osc::_OscCalcPMNSOpt< T >::SetTh12(), and osc::analytic::_OscCalc< T >::SetTh12().

template<typename T>

T osc::_IOscCalcAdjustable< T >::fTh13

protectedinherited

Definition at line 117 of file IOscCalc.h.

Referenced by osc::analytic::_OscCalc< T >::_P(), osc::OscCalcPMNS_CPT::P(), osc::_OscCalcPMNS< T >::P(), osc::_OscCalcPMNSOpt< T >::P(), osc::OscCalcPMNS_CPT::SetTh13(), osc::_OscCalcPMNS< T >::SetTh13(), osc::_OscCalcPMNSOpt< T >::SetTh13(), and osc::analytic::_OscCalc< T >::SetTh13().

template<typename T>

T osc::_IOscCalcAdjustable< T >::fTh23

protectedinherited

Definition at line 118 of file IOscCalc.h.

Referenced by osc::analytic::_OscCalc< T >::_P(), osc::OscCalcPMNS_CPT::P(), osc::_OscCalcPMNS< T >::P(), osc::_OscCalcPMNSOpt< T >::P(), osc::OscCalcPMNS_CPT::SetTh23(), osc::_OscCalcPMNS< T >::SetTh23(), osc::_OscCalcPMNSOpt< T >::SetTh23(), and osc::analytic::_OscCalc< T >::SetTh23().

---

The documentation for this class was generated from the following files:

• OscLib/OscCalcSterile.h
• OscLib/OscCalcSterile.cxx