iceCube.uhe.interactions
Class GlashowResonanceHadronic

java.lang.Object
  extended by iceCube.uhe.interactions.GlashowResonanceLeptonic
      extended by iceCube.uhe.interactions.GlashowResonanceHadronic
All Implemented Interfaces:
java.io.Serializable, Function

public class GlashowResonanceHadronic
extends GlashowResonanceLeptonic
implements Function

The Glashow Resonance reaction with W into the hadronic decay

   \bar{nu_e} + e^{-1} -> hadrons
   
The super class GlashowResonanceLeptonic class is used in most of the calculation The inelasiticity parameter y is fixed to be always 0 because all the final states are hadrons that generates cascades at once.

See Also:
Serialized Form

Field Summary
protected  double ratioOfDecayWidth
          Ration Gamma_hadron/Gamma_muon = Gamma_hadron/Gamma_tau
 
Fields inherited from class iceCube.uhe.interactions.GlashowResonanceLeptonic
chargePerNucleon, G_F, gammaW, hbar_c, isPerNucleon, massE, massW
 
Constructor Summary
GlashowResonanceHadronic(ParticlePoint s)
           
 
Method Summary
 double getDSigmaDy(double y)
          return dSigma/dy [cm^2] where y = 1 - - E_{l^{-1}}/E_{\bar{\nu_e}}
 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()
          return 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 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
 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 iceCube.uhe.interactions.GlashowResonanceLeptonic
areaFactorByWeakCoupling, calculateCrossSectionAsPerElectron, calculateCrossSectionAsPerNucleon, getYmax, getYmin, isValidInelasticity, isValidParticle, wResonance
 
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
 

Field Detail

ratioOfDecayWidth

protected final double ratioOfDecayWidth
Ration Gamma_hadron/Gamma_muon = Gamma_hadron/Gamma_tau

See Also:
Constant Field Values
Constructor Detail

GlashowResonanceHadronic

public GlashowResonanceHadronic(ParticlePoint s)
Method Detail

getDSigmaDy

public double getDSigmaDy(double y)
return dSigma/dy [cm^2] where y = 1 - - E_{l^{-1}}/E_{\bar{\nu_e}}

Overrides:
getDSigmaDy in class GlashowResonanceLeptonic

getSigma

public double getSigma()
return total cross section [cm^2]

Overrides:
getSigma in class GlashowResonanceLeptonic

interactionName

public java.lang.String interactionName()
Show Name of the Interaction

Overrides:
interactionName in class GlashowResonanceLeptonic

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()

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