ana::OscCalcSterileApprox Class Reference

#include <OscCalcSterileApprox.h>

Inheritance diagram for ana::OscCalcSterileApprox:

Public Member Functions
virtual double	P (int from, int to, double E) override
double	P_range (int from, int to, double Elo, double Ehi)
double	P_LoverE (int from, int to, double LElo, double LEhi)
virtual OscCalcSterileApprox *	Copy () const override
void	SetDmsq (double d)
double	GetDmsq () const
void	SetSinSq2ThetaMuMu (double t)
double	GetSinSq2ThetaMuMu () const
void	SetSinSq2ThetaMuE (double t)
double	GetSinSq2ThetaMuE () const
void	SetSinSq2ThetaEE (double t)
double	GetSinSq2ThetaEE () const
void	SetL (double L)
double	GetL () const
TMD5 *	GetParamsHash () const override

Protected Member Functions
double	PFromDelta (int from, int to, double Delta) const

Protected Attributes
double	fDmsq
double	fSinSq2ThetaMuMu
double	fSinSq2ThetaMuE
double	fSinSq2ThetaEE
bool	fSinSq2ThetaMuMuSet = false
bool	fSinSq2ThetaMuESet = false
bool	fSinSq2ThetaEESet = false
double	fL

Detailed Description

Definition at line 14 of file OscCalcSterileApprox.h.

Member Function Documentation

OscCalcSterileApprox * ana::OscCalcSterileApprox::Copy ( ) const

overridevirtual

Definition at line 252 of file OscCalcSterileApprox.cxx.

  {
    OscCalcSterileApprox* ret = new OscCalcSterileApprox;

    ret->fDmsq = fDmsq;
    ret->fSinSq2ThetaMuMuSet = fSinSq2ThetaMuMuSet;
    ret->fSinSq2ThetaMuESet = fSinSq2ThetaMuESet;
    ret->fSinSq2ThetaEESet = fSinSq2ThetaEESet;
    if(fSinSq2ThetaMuMuSet)
      ret->fSinSq2ThetaMuMu = fSinSq2ThetaMuMu;
    if(fSinSq2ThetaMuESet)
      ret->fSinSq2ThetaMuE = fSinSq2ThetaMuE;
    if(fSinSq2ThetaEESet)
      ret->fSinSq2ThetaEE = fSinSq2ThetaEE;
    ret->fL = fL;

    return ret;
  }

double ana::OscCalcSterileApprox::GetDmsq ( ) const

inline

Definition at line 28 of file OscCalcSterileApprox.h.

{return fDmsq;}

double ana::OscCalcSterileApprox::GetL ( ) const

inline

Definition at line 41 of file OscCalcSterileApprox.h.

{return fL;}

TMD5 * ana::OscCalcSterileApprox::GetParamsHash ( ) const

override

Definition at line 272 of file OscCalcSterileApprox.cxx.

  {
    TMD5* ret = new TMD5;
    const std::string txt = "SterileApprox";
    ret->Update((unsigned char*)txt.c_str(), txt.size());
    const int kNumParams = 4;
    double angles[2] = {0, 0};
    int angles_set = 0;
    if(fSinSq2ThetaMuMuSet && angles_set < 2)
      angles[angles_set++] = fSinSq2ThetaMuMu;
    if(fSinSq2ThetaMuESet && angles_set < 2)
      angles[angles_set++] = fSinSq2ThetaMuE;
    if(fSinSq2ThetaEESet && angles_set < 2)
      angles[angles_set++] = fSinSq2ThetaEE;
    double buf[kNumParams] = {fDmsq, angles[0], angles[1], fL};
    ret->Update((unsigned char*)buf, sizeof(double)*kNumParams);
    ret->Final();
    return ret;
  }

double ana::OscCalcSterileApprox::GetSinSq2ThetaEE

(

)

const

Definition at line 24 of file OscCalcSterileApprox.cxx.

  {
    if(fSinSq2ThetaEESet) return fSinSq2ThetaEE;
    if(!fSinSq2ThetaMuMuSet && !fSinSq2ThetaMuESet) return 0;

    // The three angles are coupled via
    //
    // ss2thmm = 4*Um4^2*(1-Um4^2)
    // ss2thee = 4*Ue4^2*(1-Ue4^2)
    // ss2thme = 4*Um4^2*Ue4^2
    //
    // Solve for the final (nue survival angle). Emperically, choosing the
    // negative sign in the expression works well.

    // A couple of limits
    if(fSinSq2ThetaMuMu <= 0) return fSinSq2ThetaMuE / 4;
    if(fSinSq2ThetaMuMu >= 1) return fSinSq2ThetaMuE / 2;

    // General case
    const double Ue4sq = fSinSq2ThetaMuE/(2*fSinSq2ThetaMuMu)*(1-sqrt(1-fSinSq2ThetaMuMu));
    assert(!isinf(Ue4sq) && !isnan(Ue4sq));
    const double sinsq2thetaee = 4.0*Ue4sq*(1-Ue4sq);

    assert(sinsq2thetaee >= 0 && sinsq2thetaee <= 1);

    // Cross-check that we actually solved the equations correctly
    if(Ue4sq > 0){
      const double Umu4sq = fSinSq2ThetaMuE/(4*Ue4sq);
      assert(fabs(4*Umu4sq*(1-Umu4sq) - fSinSq2ThetaMuMu) < 1e-6);
    }

    return sinsq2thetaee;
  }

double ana::OscCalcSterileApprox::GetSinSq2ThetaMuE

(

)

const

Definition at line 110 of file OscCalcSterileApprox.cxx.

  {
    if(fSinSq2ThetaMuESet) return fSinSq2ThetaMuE;
    if(!fSinSq2ThetaMuMuSet && !fSinSq2ThetaEESet) return 0;

    // The three angles are coupled via
    //
    // ss2thmm = 4*Um4^2*(1-Um4^2)
    // ss2thee = 4*Ue4^2*(1-Ue4^2)
    // ss2thme = 4*Um4^2*Ue4^2
    //
    // Choose 1-sqrt() case instead of 1+sqrt() since in the later
    // if both ss2thmm and ss2thee are 0 we stil get ss2thme = 1,
    // which seems wrong.

    const double Um4sq = 0.5 * (1.0 - std::sqrt(1.0 - fSinSq2ThetaMuMu));
    const double Ue4sq = 0.5 * (1.0 - std::sqrt(1.0 - fSinSq2ThetaEE));

    assert(!isinf(Um4sq) && !isnan(Um4sq));
    assert(!isinf(Ue4sq) && !isnan(Ue4sq));
    const double sinsq2thetame = 4.0*Um4sq*Ue4sq;

    assert(sinsq2thetame >= 0 && sinsq2thetame <= 1);

    return sinsq2thetame;
  }

double ana::OscCalcSterileApprox::GetSinSq2ThetaMuMu

(

)

const

Definition at line 67 of file OscCalcSterileApprox.cxx.

  {
    if(fSinSq2ThetaMuMuSet) return fSinSq2ThetaMuMu;
    if(!fSinSq2ThetaEESet && !fSinSq2ThetaMuESet) return 0;

    // The three angles are coupled via
    //
    // ss2thmm = 4*Um4^2*(1-Um4^2)
    // ss2thee = 4*Ue4^2*(1-Ue4^2)
    // ss2thme = 4*Um4^2*Ue4^2
    //
    // Solve for the final (numu survival angle). Emperically, choosing the
    // negative sign in the expression works well.

    // A couple of limits
    if(fSinSq2ThetaEE <= 0) return fSinSq2ThetaMuE / 4;
    if(fSinSq2ThetaEE >= 1) return fSinSq2ThetaMuE / 2;

    // General case
    const double Um4sq = fSinSq2ThetaMuE/(2*fSinSq2ThetaEE)*(1-sqrt(1-fSinSq2ThetaEE));
    assert(!isinf(Um4sq) && !isnan(Um4sq));
    const double sinsq2thetamm = 4.0*Um4sq*(1-Um4sq);

    assert(sinsq2thetamm >= 0 && sinsq2thetamm <= 1);

    // Cross-check that we actually solved the equations correctly
    if(Um4sq > 0){
      const double Ue4sq = fSinSq2ThetaMuE/(4*Um4sq);
      assert(fabs(4*Ue4sq*(1-Ue4sq) - fSinSq2ThetaEE) < 1e-6);
    }

    return sinsq2thetamm;
  }

double ana::OscCalcSterileApprox::P ( int  from, int  to, double  E  )

overridevirtual

Definition at line 178 of file OscCalcSterileApprox.cxx.

  {
    return P_range(from, to, E, E);
  }

double ana::OscCalcSterileApprox::P_LoverE	(	int	from,
		int	to,
		double	LElo,
		double	LEhi
	)

Definition at line 232 of file OscCalcSterileApprox.cxx.

  {
    if(fDmsq == 0) return (from == to || to == 0) ? 1 : 0;

    LElo = std::max(0., LElo);

    LElo *= 1.267*fDmsq;
    LEhi *= 1.267*fDmsq;

    if(LElo == LEhi) return PFromDelta(from, to, util::sqr(sin(LElo)));

    // Average value of sin^2 over the range
    const double Delta = (.25*(sin(2*LElo) - sin(2*LEhi)) + .5*(LEhi - LElo))/(LEhi-LElo);
    assert(!isnan(Delta));

    return PFromDelta(from, to, Delta);
  }

double ana::OscCalcSterileApprox::P_range	(	int	from,
		int	to,
		double	Elo,
		double	Ehi
	)

Definition at line 216 of file OscCalcSterileApprox.cxx.

  {
    if(Ehi <= 0) return 0;

    if(fDmsq == 0) return (from == to || to == 0) ? 1 : 0;

    Elo = std::max(0., Elo);

    assert(fL != 0);
    const double Delta = AvgSinSq(1.267*fDmsq*fL, Elo, Ehi);
    assert(!isnan(Delta));

    return PFromDelta(from, to, Delta);
  }

double ana::OscCalcSterileApprox::PFromDelta ( int  from, int  to, double  Delta  ) const

protected

Definition at line 184 of file OscCalcSterileApprox.cxx.

  {
    if(abs(from) == 14 && abs(to) == 14){
      return 1-GetSinSq2ThetaMuMu()*Delta;
    }
    else if(abs(from) == 12 && abs(to) == 12){
      return 1-GetSinSq2ThetaEE()*Delta;
    }
    else if(abs(from) == 14 && abs(to) == 12){
      return GetSinSq2ThetaMuE()*Delta;
    }
    else if(abs(from) == 12 && abs(to) == 14){
      // TODO - this seems reasonable, is it right?
      return GetSinSq2ThetaMuE()*Delta;
    }
    else if(abs(to) == 16){ // no tau appearance
      return 0;
    }

    //Option to return the active fraction
    else if (abs(from) == 14 && to == 0) {
      return (1-GetSinSq2ThetaMuMu()*Delta) + (GetSinSq2ThetaMuE()*Delta);
    }
    else if (abs(from) == 12 && to == 0) {
      return (1-GetSinSq2ThetaEE()*Delta) + (GetSinSq2ThetaMuE()*Delta);
    }

    std::cout << "OscCalcSterileApprox: P(" << from << ", " << to << ") not implemented" << std::endl;
    abort();
  }

void ana::OscCalcSterileApprox::SetDmsq ( double  d )

inline

Definition at line 27 of file OscCalcSterileApprox.h.

{fDmsq = d;}

void ana::OscCalcSterileApprox::SetL ( double  L )

inline

Definition at line 40 of file OscCalcSterileApprox.h.

{fL = L;}

void ana::OscCalcSterileApprox::SetSinSq2ThetaEE

(

double

t

)

Definition at line 16 of file OscCalcSterileApprox.cxx.

  {
    assert(!fSinSq2ThetaMuMuSet || !fSinSq2ThetaMuESet);
    fSinSq2ThetaEESet = true;
    fSinSq2ThetaEE = t;
  }

void ana::OscCalcSterileApprox::SetSinSq2ThetaMuE

(

double

t

)

Definition at line 102 of file OscCalcSterileApprox.cxx.

  {
    assert(!fSinSq2ThetaMuMuSet || !fSinSq2ThetaEESet);
    fSinSq2ThetaMuESet = true;
    fSinSq2ThetaMuE = t;
  }

void ana::OscCalcSterileApprox::SetSinSq2ThetaMuMu

(

double

t

)

Definition at line 59 of file OscCalcSterileApprox.cxx.

  {
    assert(!fSinSq2ThetaEESet || ! fSinSq2ThetaMuESet);
    fSinSq2ThetaMuMuSet = true;
    fSinSq2ThetaMuMu = t;
  }

Member Data Documentation

double ana::OscCalcSterileApprox::fDmsq

protected

Definition at line 47 of file OscCalcSterileApprox.h.

double ana::OscCalcSterileApprox::fL

protected

Definition at line 54 of file OscCalcSterileApprox.h.

double ana::OscCalcSterileApprox::fSinSq2ThetaEE

protected

Definition at line 50 of file OscCalcSterileApprox.h.

bool ana::OscCalcSterileApprox::fSinSq2ThetaEESet = false

protected

Definition at line 53 of file OscCalcSterileApprox.h.

double ana::OscCalcSterileApprox::fSinSq2ThetaMuE

protected

Definition at line 49 of file OscCalcSterileApprox.h.

bool ana::OscCalcSterileApprox::fSinSq2ThetaMuESet = false

protected

Definition at line 52 of file OscCalcSterileApprox.h.

double ana::OscCalcSterileApprox::fSinSq2ThetaMuMu

protected

Definition at line 48 of file OscCalcSterileApprox.h.

bool ana::OscCalcSterileApprox::fSinSq2ThetaMuMuSet = false

protected

Definition at line 51 of file OscCalcSterileApprox.h.

---

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

• OscCalcSterileApprox.h
• OscCalcSterileApprox.cxx