osc::OscCalcDumb Class Reference

Simple oscillation probability calculator that has no solar term or mass hierarchy or delta so it's some kind of average of all of those. More...

#include "/cvmfs/nova.opensciencegrid.org/externals/osclib/v00.07/src/OscLib/OscCalcDumb.h"

Inheritance diagram for osc::OscCalcDumb:

Public Member Functions
virtual IOscCalc *	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 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)
virtual TMD5 *	GetParamsHash () const
	Use to check two calculators are in the same state. More...

Protected Member Functions
template<class T >
T	sqr (T x) const

Detailed Description

Simple oscillation probability calculator that has no solar term or mass hierarchy or delta so it's some kind of average of all of those.

Definition at line 16 of file OscCalcDumb.h.

Member Function Documentation

virtual IOscCalc* osc::OscCalcDumb::Copy ( ) const

inlineoverridevirtual

Implements osc::_IOscCalc< T >.

Definition at line 19 of file OscCalcDumb.h.

References osc::_IOscCalc< T >::P().

Referenced by futureSig_getNueNumuSyst(), genie_preds_make(), GetCalculator(), make_prediction(), make_prediction_rhc(), plot_joint_fit_2020_contours(), plot_spectra_2dplots(), syst_table_fit(), syst_table_fit_new(), and test_numupred2017().

{return new OscCalcDumb;}

template<typename T>

virtual TMD5* osc::_IOscCalc< T >::GetParamsHash ( ) const

inlinevirtualinherited

Use to check two calculators are in the same state.

Returns

Null means not implemented for this calculator

Reimplemented in osc::analytic::_OscCalc< T >, osc::OscCalcCPT, osc::OscCalcPMNSOptEigen, osc::_OscCalcDMP< T >, osc::OscCalc, osc::_NoOscillations< T >, osc::OscCalcSterile, osc::OscCalcPMNS_CPT, osc::OscCalcGeneral, osc::_OscCalcPMNSOpt< T >, osc::OscCalcSterileBeam, and osc::_OscCalcPMNS< T >.

Definition at line 39 of file IOscCalc.h.

Referenced by ana::OscillatableSpectrum::_Oscillated(), ana::PredictionInterp::_PredictComponentSyst(), and osc::OscCalcCPT::GetParamsHash().

{return 0;}

virtual double osc::OscCalcDumb::P ( int  flavBefore, int  flavAfter, double  E  )

inlineoverridevirtual

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

Implements osc::_IOscCalc< T >.

Definition at line 21 of file OscCalcDumb.h.

References abs(), CLHEP::L, std::sin(), sqr(), and ssth23.

Referenced by DrawCCOscillations(), ana::CovMxSurface::DrawSurfacePoint(), NuMu2020_TrimCAFs(), prob_biprob_pedagogical(), and caf::CAFMaker::SimpleOscProb().

    {
      const double L      = 810.;    // km
      const double ldm    = 2.40e-3; // large delta m^2, m23
      const double ss2t13 = 0.1;     // sin^2(2theta_13)
      const double ss2t23 = 1.;      // maximal sin^2(2theta_23)
      const double ssth23 = 0.5;     // sin^2(theta_23) corresponding to above

      // Signal
      if(abs(flavAfter) == 12 && abs(flavBefore) == 14){
        // Nue appearance
        return ss2t13*ssth23*sqr(sin(1.267*ldm*L/E));
      }
      if(abs(flavAfter) == 14 && abs(flavBefore) == 14){
        // CC mu disappearance
        return 1-ss2t23*sqr(sin(1.267*ldm*L/E));
      }

      // Background
      if(abs(flavAfter) == 12 && abs(flavBefore) == 12){
        // Beam nue
        return 1-ss2t13*sqr(sin(1.267*ldm*L/E));
      }
      if(abs(flavAfter) == 14 && abs(flavBefore) == 12){
        // CC mu appearance
        return ss2t13*ssth23*sqr(sin(1.267*ldm*L/E));
      }
      if(abs(flavAfter) == 16 && abs(flavBefore) == 14){
        //numu to nutau CC appearance
        return (1-ss2t13)*ss2t23*sqr(sin(1.267*ldm*L/E));
      }
      if(abs(flavAfter) == 16 && abs(flavBefore) == 12){
        //nue to nutau CC appearance
        return ss2t13*(1-ssth23)*sqr(sin(1.267*ldm*L/E));
      }
      if(abs(flavAfter) == 16 && abs(flavBefore) == 16){
        //nutau to nutau CC disappearance
        return 1-(ss2t23-ss2t23*ss2t13+ss2t13-ss2t13*ssth23)*sqr(sin(1.267*ldm*L/E));
      }

      // Don't know what this is
      return 0;
    }

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 27 of file IOscCalc.cxx.

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

Referenced by osc::_IOscCalc< T >::P().

  {
    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 37 of file IOscCalc.cxx.

References P.

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

template<class T >

T osc::OscCalcDumb::sqr ( T  x ) const

inlineprotected

Definition at line 66 of file OscCalcDumb.h.

References submit_syst::x.

Referenced by P().

{return x*x;}

---

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

• OscLib/OscCalcDumb.h