 |
Public Member Functions | |
| Taumodel () | |
| double | GetTauWeight (Primaries *primary1, Settings *settings1, IceModel *antarctica1, Interaction *interaction1, double pnu, int nu_nubar, double &ptauf, int &crust_entered) |
| GetTauWeight is the function that will calculate the probability that a tau neutrino will interact along its path through the earth,and the tau will survive the rest of the journey and decay in the ice. This probability is calculated for final energies from 10^15.5 to the energy of the neutrino. More... | |
| ClassDef (Taumodel, 1) | |
Public Attributes | |
| double | ptauf |
| Final energy of the tau. | |
| double | weight_tau_prob |
| Weight for tau neutrino to interact, create a tau, tau survives and decays in the ice. | |
| int | inu |
| double | weight_nu_prob |
Definition at line 29 of file Taumodel.hh.
| Taumodel::Taumodel | ( | ) |
For Total Tau Survival probability equation n.b. not in SI units. //from Tau neutrino propagaiton and tau energy loss 2005 Dutta, Huang, & Reno. Equation 16 & used in Equation 30.
| m^2/kg |
| m^2/kg | }parameterization using a logarithmic dependence on energy for B,
| eV | the tau elecromagnetic energy loss parameter.
| eV |Mass of Tau
| m |Tau Decay length (86.93 microMeters) | |
| g |nucleon/ proton mass in grams,also equal to 0.938 GeV.
| none |constant that sets the total probability to unity
these last two constants from Connolly Calc 2011, used in d_dzPsurvNu().
Definition at line 50 of file Taumodel.cc.
| double Taumodel::GetTauWeight | ( | Primaries * | primary1, |
| Settings * | settings1, | ||
| IceModel * | antarctica1, | ||
| Interaction * | interaction1, | ||
| double | pnu, | ||
| int | nu_nubar, | ||
| double & | ptauf, | ||
| int & | crust_entered | ||
| ) |
GetTauWeight is the function that will calculate the probability that a tau neutrino will interact along its path through the earth,and the tau will survive the rest of the journey and decay in the ice. This probability is calculated for final energies from 10^15.5 to the energy of the neutrino.
GetTauWeight is the function that will calculate the probability that a tau neutrino will interact along its path through the earth,and the tau will survive the rest of the journey and decay in the ice. This probability is calculated for final energies from 10^15.5 to the energy of the neutrino.
vector from earth_in to "interaction point". Sets the path direction
normalized chord3
Bring in useful variables from other classes
Find the chord, its length and its unit vector.
posnu-earth_in;
normalized chord3
total tau survival probability.
how big is the step size
set up stuff to be used later.
where the particle enters the ice.
Get the density vectors this function needs
set energy vector
set tau surv vector
Set the starting position for tau. First possible interaction point for neutrino
step number to get the correct initial energy at that point;
vector pointing to the step we are currently on.
average density neutrino has seen to that step
density of the current step
Energy vector is filled backwards (final to initial), pull out correct energy
chance tau survives from creation to ice
inelasticity
if neutrino interacts in ice, the tau will already be in ice.
prob neutrino survives to this point
how the initial tau energy is related to final
probability that tau survives to ice
chance to decay in ice
chance to decay in ice if already in ice
total prob for neutrino to make tau and for tau to decay in ice at this energy
dP/dE ->dP/d(log10(E))
dP/dlog10(E) * dlog10(E)
We have the weight. Now use a PDF to find the final energy of the tau.///
Definition at line 82 of file Taumodel.cc.
1.8.11
