gen
Class EarthModel

java.lang.Object
  extended by mmc.PhysicsModel
      extended by gen.EarthModel
All Implemented Interfaces:
FunctionInt, FunctionInt2, FunctionOfx

public class EarthModel
extends PhysicsModel

Earth density model implementation.


Field Summary
 Interpolate J
           
 boolean jt
           
 double[] mRa
           
 double[] mRo
           
 int num
           
 double px
           
 double py
           
 double pz
           
static double R0
           
 double tf
           
 double ti
           
 
Fields inherited from class mmc.PhysicsModel
Alpha, bigEnergy, C, Cmon, computerPrecision, De2, Ds2, ebig, elow, Gf, Gw, Gz, halfPrecision, imaxs, iprec, iprec2, iromb, K, Lmu, Log10, Lstau, Ltau, Ma1, Me, Mmon, Mmu, Mn, Mp, Mpi, Mrh, Mrs, Mstau, Mtau, Mw, Mz, Na, nlow, num1, num2, num3, Pi, Re, Ry, sqrt2, sqrt3, sqrtE, St2, Tt2, xres, Xw
 
Constructor Summary
EarthModel(int model, double z0, double b0, double D)
          Initialize class with atmospheric model, ground elevation z0 in [km], bedrock elevation b0 in [km], and detector depth D.
 
Method Summary
 double function(double x)
          Returns the density in [g/cm^3] as a function of the path length in [km] - interface to Integral.
 double functionInt(double x, double t)
          2d parametrization - interface to Interpolate
 void interpolate()
          Parametrizes the integral of this class.
 double intRho(double x)
          Returns the density in [g/cm^3] as a function of the path length in [km], gives density of the internal medium.
static void main(java.lang.String[] args)
          Prints the density profile.
 double r(double x, double y, double z)
          Calculates path length in [km] from the entry point.
 double rho(double R)
          Calculates density in [g/cm^3] as a function of distance to the center R in [km] according to the preliminary Earth model.
 void sett(double px, double py, double pz)
          Calculates and stores path length in [km] to the entry and exit points.
 double t()
          Returns the path length in [km] corresponding to the mass overburden X0 in [g/cm^2] set by a call to X(t1, t2, X0).
 double X(double t1, double t2)
          Calculates mass overburden in [g/cm^2] as a function of path length in [km] to the entry and exit points.
 double X(double t1, double t2, double X0)
          Calculates mass overburden in [g/cm^2] as a function of path length in [km] to the entry and exit points.
 
Methods inherited from class mmc.PhysicsModel
functionInt
 
Methods inherited from class java.lang.Object
clone, equals, finalize, getClass, hashCode, notify, notifyAll, toString, wait, wait, wait
 

Field Detail

R0

public static final double R0
See Also:
Constant Field Values

num

public int num

mRa

public double[] mRa

mRo

public double[] mRo

px

public double px

py

public double py

pz

public double pz

ti

public double ti

tf

public double tf

J

public Interpolate J

jt

public boolean jt
Constructor Detail

EarthModel

public EarthModel(int model,
                  double z0,
                  double b0,
                  double D)
Initialize class with atmospheric model, ground elevation z0 in [km], bedrock elevation b0 in [km], and detector depth D.

Method Detail

rho

public double rho(double R)
Calculates density in [g/cm^3] as a function of distance to the center R in [km] according to the preliminary Earth model.


sett

public void sett(double px,
                 double py,
                 double pz)
Calculates and stores path length in [km] to the entry and exit points.


r

public double r(double x,
                double y,
                double z)
Calculates path length in [km] from the entry point. sett must be called first.


X

public double X(double t1,
                double t2)
Calculates mass overburden in [g/cm^2] as a function of path length in [km] to the entry and exit points.


X

public double X(double t1,
                double t2,
                double X0)
Calculates mass overburden in [g/cm^2] as a function of path length in [km] to the entry and exit points. Also calculates path length corresponding to the mass overburden X0 in [g/cm^2] from the entry point.


t

public double t()
Returns the path length in [km] corresponding to the mass overburden X0 in [g/cm^2] set by a call to X(t1, t2, X0).


function

public double function(double x)
Returns the density in [g/cm^3] as a function of the path length in [km] - interface to Integral.

Specified by:
function in interface FunctionOfx
Overrides:
function in class PhysicsModel

intRho

public double intRho(double x)
Returns the density in [g/cm^3] as a function of the path length in [km], gives density of the internal medium.


main

public static void main(java.lang.String[] args)
Prints the density profile.


interpolate

public void interpolate()
Parametrizes the integral of this class.


functionInt

public double functionInt(double x,
                          double t)
2d parametrization - interface to Interpolate

Specified by:
functionInt in interface FunctionInt2
Overrides:
functionInt in class PhysicsModel