A struct holding useful combinations of trajectory variables. More...

#include <EarthModelBinned.h>

Public Member Functions
	TrajConstants (double cosT=0, double phi=0, double DetLat=0, double DetLon=0, double rDet=6368)
	Constructor. More...

void	UpdateNuAngles (double cosTheta, double phi)

void	UpdateDetPos (double rDet, double DetLat, double DetLon)

void	Recalculate ()

Public Attributes
double	cosT
	cosT More...

double	cosA
	cos(phi) More...

double	sinA
	sin(phi) More...

double	sinDetLat
	sin(DetLat) More...

double	cosDetLon
	cos(DetLon) More...

double	sinDetLon
	sin(DetLon) More...

double	DetRadius
	rDet More...

double	sinSqT
	sin^2(T) = 1 - (cosT)^2 More...

double	sinT
	sin(T) = sqrt(sinSqT) More...

double	sinTsinA
	sin(T)*sin(phi) More...

double	sinTcosA
	sin(T)*cos(phi) More...

double	cosDetLat
	cos(DetLat) More...

double	rDetSinDetLat
	rDet*sin(DetLat) More...

double	rDetCosDetLat
	rDet*cos(DetLat) More...

double	sinTsinAsinDetLon
	sin(T)sin(phi)sin(DetLon) More...

double	sinTsinAcosDetLon
	sin(T)sin(phi)cos(DetLon) More...

double	cosTcosDetLat
	cosT*cos(DetLat) More...

double	rDetCosT
	rDet*cosT More...

double	rDetSinT
	rDet*sin(T) More...

double	alpha
	sin(T)cos(phi)sin(DetLat)-cosT*cos(DetLat) More...

double	beta
	sin(T)sin(DetLat)-cos(phi)cosT*cos(DetLat) More...

double	gamma
	sin(T)cos(phi)cos(DetLat)+cosT*sin(DetLat) More...

double	gammaSq
	[sin(T)cos(phi)cos(DetLat)+cosT*sin(DetLat)]^2 More...

double	rDetGammaSinDetLat
	rDetsin(DetLat)[sin(T)cos(phi)cos(DetLat)+cosT*sin(DetLat)] More...

double	maxSinSqLat
	1 - [sin(phi)*cos(DetLat)]^2 More...

double	xlatextreme
	rDetcos(phi)cos(DetLat)/[sin(T)sin(DetLat)-cos(phi)cosT*cos(DetLat)], unless the denominator is 0, in which case it is set to 0 More...

Detailed Description

This struct holds combinations of variables that are useful when doing the calculations needed for FindPath within EarthModelBinned.

Author: rpestes@apc.in2p3.fr

Definition at line 80 of file EarthModelBinned.h.

Constructor & Destructor Documentation

◆ TrajConstants()

OscProb::TrajConstants::TrajConstants	(	double	cosT = `0`,
		double	phi = `0`,
		double	DetLat = `0`,
		double	DetLon = `0`,
		double	rDet = `6368`
	)

inline

Constructor.

By default, it sets cosT, Az, latD, and lonD equal to 0, and it uses 6368km for rD.

Parameters

cosT	- Cosine of the zenith angle for the neutrino trajectory
phi	- The azimuthal angle for the neutrino trajectory (in rad)
DetLat	- The latitude of the detector (in rad)
DetLon	- The longitude of the detector (in rad)
rDet	- The distance from the center of the Earth to the detector (in km)

Definition at line 96 of file EarthModelBinned.h.

      {
        UpdateNuAngles(cosT, phi);
        UpdateDetPos(DetLat, DetLon, rDet);
      }

References cosT, UpdateDetPos(), and UpdateNuAngles().

Member Function Documentation

◆ Recalculate()

void TrajConstants::Recalculate ( )

Calculate constants that use combinations of detector position variables and neutrino direction angles

Recalculate constants that use combinations of detector position variables and neutrino direction angles.

Definition at line 66 of file EarthModelBinned.cxx.

{
  cosTcosDetLat      = cosT * cosDetLat;
  sinTsinAsinDetLon  = sinTsinA * sinDetLon;
  sinTsinAcosDetLon  = sinTsinA * cosDetLon;
  rDetCosT           = DetRadius * cosT;
  rDetSinT           = DetRadius * sinT;
  alpha              = sinTcosA * sinDetLat - cosTcosDetLat;
  beta               = sinT * sinDetLat - cosA * cosTcosDetLat;
  gamma              = sinTcosA * cosDetLat + cosT * sinDetLat;
  gammaSq            = gamma * gamma;
  rDetGammaSinDetLat = rDetSinDetLat * gamma;
  maxSinSqLat        = 1 - pow(sinA * cosDetLat, 2);
  if(beta == 0) { xlatextreme = 0; }
  else { xlatextreme = cosA * rDetCosDetLat / beta; }
}

References alpha, beta, cosA, cosDetLat, cosDetLon, cosT, cosTcosDetLat, DetRadius, gamma, gammaSq, maxSinSqLat, rDetCosDetLat, rDetCosT, rDetGammaSinDetLat, rDetSinDetLat, rDetSinT, sinA, sinDetLat, sinDetLon, sinT, sinTcosA, sinTsinA, sinTsinAcosDetLon, sinTsinAsinDetLon, and xlatextreme.

Referenced by OscProb::EarthModelBinned::FillPath().

◆ UpdateDetPos()

void TrajConstants::UpdateDetPos	(	double	rDet,
		double	DetLat,
		double	DetLon
	)

Update detector position for neutrino trajectory calculations

Update detector position for neutrino trajectory calculations.

Parameters

rDet	- The distance from the center of the Earth to the detector (in km)
DetLat	- The latitude of the detector (in rad)
DetLon	- The longitude of the detector (in rad)

Definition at line 49 of file EarthModelBinned.cxx.

{
  sinDetLat = sin(DetLat);
  cosDetLat = cos(DetLat);
  cosDetLon = cos(DetLon);
  sinDetLon = sin(DetLon);
  DetRadius = rDet;
 
  rDetSinDetLat = DetRadius * sinDetLat;
  rDetCosDetLat = DetRadius * cosDetLat;
}

References cosDetLat, cosDetLon, DetRadius, rDetCosDetLat, rDetSinDetLat, sinDetLat, and sinDetLon.

Referenced by OscProb::EarthModelBinned::SetDetPos(), and TrajConstants().

◆ UpdateNuAngles()

void TrajConstants::UpdateNuAngles	(	double	cosTheta,
		double	phi
	)

Update values of zenith angle and azimuthal angle for neutrino trajectory

Update values of zenith angle and azimuthal angle for neutrino trajectory.

Parameters

cosTheta	- Cosine of the zenith angle
phi	- The azimuthal angle in radians

Definition at line 27 of file EarthModelBinned.cxx.

{
  cosT   = cosTheta;
  sinSqT = 1 - cosT * cosT;
  sinT   = sqrt(sinSqT);
 
  cosA = cos(phi);
  sinA = sin(phi);
 
  sinTsinA = sinT * sinA;
  sinTcosA = sinT * cosA;
}