iceCube.uhe.muonModel
public class PropagatingAtmMuonFlux extends java.lang.Object implements Function
Constructor and Description |
---|
PropagatingAtmMuonFlux()
Constructor
|
PropagatingAtmMuonFlux(double muonThresholdEnergy,
boolean cutoffExists)
Constructor
|
PropagatingAtmMuonFlux(double alpha,
double muonThresholdEnergy,
boolean cutoffExists)
Constructor
|
Modifier and Type | Method and Description |
---|---|
protected double |
getAverageMuonEnergyLossAfterPropagation()
Calculate and return the averaged relative energy loss of a muon
after its propagation in the Earth.
|
double |
getDFMuDLogCREDLogE(double logCosmicRayEnergy,
double halfWidthOfLogE,
double logE,
double cos_th)
Calculate dF^2/dLogEcrDLogE [/cm^2 sec sr] for nu-mu
Because the cosmic ray follows the rapidly falling spectrum, the flux is
averaged over a given log(cosmic ray energy [GeV]) in a half width given
by halfWidthOfLogE.
|
double |
getDFMuDLogE(double logE,
double cos_th)
Calculate dF/dLogE [/cm^2 sec sr] for nu-mu
|
double |
getFunction(int functionIndex,
double[] parameters,
double x)
Implementation of interface numRecipes.Function
used for numerical integration of the d^F/dLogEcrDlogEmu over
logEcr, the cosmic ray energy.
|
protected void |
listFluxes(double logEnergy,
double cosTheta)
integrate getDFDLogCREDLogE(logEcr,logEmu,cosTheta) over logEcr
and compare with getDFDLogE(logEmu,cosTheta) for the consistency check
|
void |
readMatrix(java.io.DataInputStream in)
Read the calculated propagation matrix
|
void |
setConfidenceLevelOfFluctuation(double prob)
Set Confidence level of the Energy rato parameter
R = (E_muon/E0)/ Bar(E_muon/E0).
|
void |
setCutOffFeature(boolean cutoffExists)
set the GZK cutoff feature or not
|
void |
setFluxCalculationMode(boolean includeFluctuationEffects,
boolean fluctuateEventByEvent)
set the mode on the flux caulation
|
void |
whetherPropagationMatrixWithGlashowResonance(boolean flag)
or Glashow Resonance
|
public PropagatingAtmMuonFlux()
public PropagatingAtmMuonFlux(double muonThresholdEnergy, boolean cutoffExists)
public PropagatingAtmMuonFlux(double alpha, double muonThresholdEnergy, boolean cutoffExists)
public void readMatrix(java.io.DataInputStream in) throws java.io.IOException
java.io.IOException
public void setFluxCalculationMode(boolean includeFluctuationEffects, boolean fluctuateEventByEvent)
boolean includeFluctuationEffects : true default. calculate the flux taking into account flucuation of muon energies due to the EAS cascading. CascadeFluctuationFactory class does this part of calculation. false calculate the bare flux given by the Elbert formula. an energy of muon bundles is associated with the energy of primary cosmic ray by one-on-one relation. boolean fluctuateEventByEvent : Valid when includeFluctuationEffects = true true : Flux is calculated such that an event by event flucuation is included. Useful to evaluate flucuation of number of background events false: default. the flucuation effects is included in average way. i.e., The systematic factor of the flux shift due the cascade flucuation is calculated and multiplied in the flux calculation. Useful to the MC-data fitting, calculation of average number of events
public void setConfidenceLevelOfFluctuation(double prob)
public void setCutOffFeature(boolean cutoffExists)
public double getDFMuDLogE(double logE, double cos_th)
public double getDFMuDLogCREDLogE(double logCosmicRayEnergy, double halfWidthOfLogE, double logE, double cos_th)
double logCosmicRayEnergy : log10(Cosmic Ray Energy [GeV]) producing muon bundles with logE double halfWidthOfLogE : a half width of integration for averaging the cosmic ray flux double logE : log10(Muon Bundle Energy [GeV]) double cos_th : cos(zenith angle)
protected double getAverageMuonEnergyLossAfterPropagation()
AtmMuonBundle.getEffectiveEnergyOfCRs(double muonBundleEnergy, double cosTheta_ice3)
public void whetherPropagationMatrixWithGlashowResonance(boolean flag)
protected void listFluxes(double logEnergy, double cosTheta)
public double getFunction(int functionIndex, double[] parameters, double x)
double x = logCosmicRayEnergy parameter[0] : log10(Muon Energy [GeV]) parameter[1] : cosine(Zenith angle)
getFunction
in interface Function