iceCube.uhe.interactions
Class PhotoNuclear

java.lang.Object
  iceCube.uhe.interactions.PhotoNuclear

All Implemented Interfaces:

java.io.Serializable, Function

public class PhotoNuclear
extends java.lang.Object
implements Function

See Also:

Serialized Form

Constructor Summary

PhotoNuclear(Particle p, ParticlePoint s)
Constructor

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 the name of 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

PhotoNuclear

public PhotoNuclear(Particle p,
                    ParticlePoint s)
             throws java.io.IOException

Constructor

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 muon and tauon but no neutrinos and pions are allowed to be involved with the photo-nuclear interaction in the medium.

interactionName

public java.lang.String interactionName()

Show the name of 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