Ionization

java.lang.Object
- iceCube.uhe.interactions.Ionization

All Implemented Interfaces:

java.io.Serializable, Function
```
public class Ionization
extends java.lang.Object
```
The Ionization loss formula for JULIeT application. This code introduces some tweaks to handle the coinuous processes like the ionization as the Interactions class has been originaly designed for stochastic interaction process which dominates anyway in ultra-high energy region. 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).

See Also:
Serialized Form

Field Summary

Fields
Modifier and Type	Field and Description
`protected static double[]`	`a_screen`
`protected static double[]`	`c_screen` The coefficinets in the screening effects
`protected static double[]`	`IonizedPotential` Ion potential table
`protected static double`	`Kion` The dimension constatnt in the dE/dX.
`protected static double`	`log10_screen`
`protected static double[]`	`m_screen`
`protected static double[]`	`x0_screen`
`protected static double[]`	`x1_screen`

Constructor Summary

Constructors
Constructor and Description

Ionization(Particle p, ParticlePoint s)
Constructor.

Constructors
Constructor and Description
`Ionization(Particle p, ParticlePoint s)` Constructor.

Method Summary

Methods
Modifier and Type	Method and Description
`double`	`getDEDX()` dE/dX - A modified Bethe-Bloch formula [GeV/g/cm^2]
`double`	`getDNDyDX(double y)` the ionization y distribution
`double`	`getDSigmaDy(double y)` Differential Cross section.
`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()`
`protected double`	`getScreenFactor()`
`double`	`getSigma()` Total cross section [cm^2]
`double`	`getYcut()`
`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`	`getYupper()`
`double`	`getZDSigmaDZ(double lowerZ, double upperZ)` Integral z*dSigma/dz over a given range to obtain the inelasticity distribution.
`double`	`inelasticCorrection(double y)` inelasticity correction
`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

Field Detail

IonizedPotential

protected static final double[] IonizedPotential

Ion potential table

c_screen
```
protected static final double[] c_screen
```
The coefficinets in the screening effects

a_screen

protected static final double[] a_screen

m_screen

protected static final double[] m_screen

log10_screen

protected static final double log10_screen

x0_screen

protected static final double[] x0_screen

x1_screen

protected static final double[] x1_screen

Kion
```
protected static final double Kion
```
The dimension constatnt in the dE/dX.

See Also:
Constant Field Values

Constructor Detail
- Ionization
```
public Ionization(Particle p,
          ParticlePoint s)
```
  Constructor. Setting the coefficients depending on the propagation medium.

Method Detail

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)

getYcut
```
public double getYcut()
```

getYmin
```
public double getYmin()
```
Getting the range of allowed y for a given interaction

getYmax
```
public double getYmax()
```

getYupper
```
public double getYupper()
```

getDEDX
```
public double getDEDX()
```
dE/dX - A modified Bethe-Bloch formula [GeV/g/cm^2]

getDNDyDX
```
public double getDNDyDX(double y)
```
the ionization y distribution

inelasticCorrection

public double inelasticCorrection(double y)

inelasticity correction

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 
        functionIndex     5    dN/dYdX(y=exp x) * coorection(y=exp x) y

Specified by:: getFunction in interface Function

getScreenFactor
```
protected double getScreenFactor()
```

getDSigmaDy
```
public double getDSigmaDy(double y)
```
Differential Cross section. This is not a true cross section but defined to give dE/dX, because the ionization is continuous process.

isValidInelasticity
```
public boolean isValidInelasticity(double y)
```
Checking the range of the given inelasticity y that is determined in an individual interaction channel.

isValidParticle
```
public boolean isValidParticle(Particle p)
```
Checking the particle kind involved with a given interaction. Only charged particles but not neutrinos are allowed to be involved with the ionization in the medium.

interactionName
```
public java.lang.String interactionName()
```
Show Name of the Interaction. Returns the Sting where the interation name is written.

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

showIntegralErrorMessage

public void showIntegralErrorMessage(double lowerY,
                            double upperY)

Error message utility

Class Ionization

Field Summary

Constructor Summary

Method Summary

Methods inherited from class java.lang.Object

Field Detail

IonizedPotential

c_screen

a_screen

m_screen

log10_screen

x0_screen

x1_screen

Kion

Constructor Detail

Ionization

Method Detail

setIncidentParticleEnergy

setIncidentParticleEnergy

getYcut

getYmin

getYmax

getYupper

getDEDX

getDNDyDX

inelasticCorrection

getFunction

getScreenFactor

getDSigmaDy

isValidInelasticity

isValidParticle

interactionName

getIncidentParticleEnergy

getDSigmaDz

setEnergyCut

getEnergyCut

getSigma

integralDSigmaDy

integralDSigmaDz

getYDSigmaDy

getZDSigmaDZ

showIntegralErrorMessage