iceCube.uhe.interactions
Class PairCreationFit

java.lang.Object
  iceCube.uhe.interactions.PairCreationFit

All Implemented Interfaces:

java.io.Serializable, Function

public class PairCreationFit
extends java.lang.Object
implements Function

    The Pair Creation interactions
    taking place in undergound rock and ice are calculated in this subclass.

    The class variables  Particle, and Point are necessary for the methods
    described here because the cross section depends on the particle propaty
    and the medium like Z and A(atomic number).

    This class, different from PairCreation.class, reads the numerical table
    generated by MakePairCreationTable.java and calculate the relevant cross sections
    with numericaly interpolating the table.
    

See Also:

Serialized Form

Constructor Summary

PairCreationFit(Particle p, ParticlePoint s, int Flavor)
Constructor: Register the Particle and ParticlePoint classes.

Method Summary

double	getDSigmaDy(double y) Differential cross section dsigma/dy [cm^2] y = 1 - Erecoiling/Eincoming -- inelasticity parameter
double	getDSigmaDz(double z) Differential cross section dsigma/dz [cm^2] z = Erecoiling/Eincoming -- inelasticity parameter
double	getEnergyCut()
double	getFunction(int functionIndex, double[] parameters, double x) Method for interface .
double	getIncidentParticleEnergy()
double	getSigma() Total cross section [cm^2]
double	getYDSigmaDy(double lowerY, double upperY) Integral y*dSigma/dy over a given range to obtain the inelasticity distribution.
double	getYmax()
double	getYmin() Getting the range of allowed y for a given interaction
double	getZDSigmaDZ(double lowerZ, double upperZ) Integral z*dSigma/dz over a given range to obtain the inelasticity distribution.
double	integralDSigmaDy(double lowerY, double upperY) Integral dSigma/dy over a given range to obtain a partial cross section
double	integralDSigmaDz(double lowerZ, double upperZ) Integral dSigma/dz over a given range to obtain a partial cross section
java.lang.String	interactionName() Show Name of the Interaction
boolean	isValidInelasticity(double y) Checking the range of the given inelasticity y that is determined in an individual interaction channel.
boolean	isValidParticle(Particle p) Checking the particle kind involved with a given interaction.
void	setEnergyCut(double cutEnergy) Energy Cut Parameter in integration to obtain the total cross section.
void	setIncidentParticleEnergy(double energy) Set The incident particle energy [GeV].
void	setIncidentParticleEnergy(int iLogE)
void	showIntegralErrorMessage(double lowerY, double upperY) Error message utility

Methods inherited from class java.lang.Object

clone, equals, finalize, getClass, hashCode, notify, notifyAll, toString, wait, wait, wait

Methods inherited from interface numRecipes.Function

getFunction

Constructor Detail

PairCreationFit

public PairCreationFit(Particle p,
                       ParticlePoint s,
                       int Flavor)
                throws java.io.IOException

Constructor: Register the Particle and ParticlePoint classes. It also reads the pre-calculated y/E * dsigma/dy from the data file.

Throws:

java.io.IOException

Method Detail

getDSigmaDy

public double getDSigmaDy(double y)

Differential cross section dsigma/dy [cm^2] y = 1 - Erecoiling/Eincoming -- inelasticity parameter

isValidInelasticity

public boolean isValidInelasticity(double y)

Checking the range of the given inelasticity y that is determined in an individual interaction channel.

getYmin

public double getYmin()

Getting the range of allowed y for a given interaction

getYmax

public double getYmax()

isValidParticle

public boolean isValidParticle(Particle p)

Checking the particle kind involved with a given interaction. Only muons and taus are allowed to be involved with the pair creation in the medium.

interactionName

public java.lang.String interactionName()

Show Name of the Interaction

setIncidentParticleEnergy

public void setIncidentParticleEnergy(double energy)

        Set The incident particle energy [GeV].
        The default value has been given by the constructor
        Interactions( ) with p.getEnergy( ) where p is 
        the Particle class. You might want to set, however,
        a different value such as 
        logE = logEnergyMinimum + deltaLogE*ilogE
        where ilogE is i'th index of the logEnergyMatrix
        in the Particle class p. This method provides you
        with a way to put the incident particle energy.
        

setIncidentParticleEnergy

public void setIncidentParticleEnergy(int iLogE)

getIncidentParticleEnergy

public double getIncidentParticleEnergy()

getDSigmaDz

public double getDSigmaDz(double z)

        Differential cross section dsigma/dz [cm^2]
        z = Erecoiling/Eincoming -- inelasticity parameter
        

setEnergyCut

public void setEnergyCut(double cutEnergy)

Energy Cut Parameter in integration to obtain the total cross section.

getEnergyCut

public double getEnergyCut()

getSigma

public double getSigma()

Total cross section [cm^2]

integralDSigmaDy

public double integralDSigmaDy(double lowerY,
                               double upperY)

Integral dSigma/dy over a given range to obtain a partial cross section

integralDSigmaDz

public double integralDSigmaDz(double lowerZ,
                               double upperZ)

Integral dSigma/dz over a given range to obtain a partial cross section

getYDSigmaDy

public double getYDSigmaDy(double lowerY,
                           double upperY)

Integral y*dSigma/dy over a given range to obtain the inelasticity distribution. The energy transfer probability would also require this value.

getZDSigmaDZ

public double getZDSigmaDZ(double lowerZ,
                           double upperZ)

Integral z*dSigma/dz over a given range to obtain the inelasticity distribution. The energy transfer probability would also require this value. z = 1-y

getFunction

public double getFunction(int functionIndex,
                          double[] parameters,
                          double x)

        Method for interface . 
        Interface the differntial cross sections given here
        to the utility methods such as the Romberg
        Integration code that is desinged for a genereal
        function in form of Func(x).

        functionIndex     1     dsigma/dy
        functionIndex     2     dsigma/dz  z = y-1
        functionIndex     3     y x dsigma/dy 
        functionIndex     4     z x dsigma/dz 
        

Specified by:

getFunction in interface Function

showIntegralErrorMessage

public void showIntegralErrorMessage(double lowerY,
                                     double upperY)

Error message utility