pisa.stages.absorption package

Submodules

pisa.stages.absorption.earth_absorption module

PISA pi stage for the calculation of earth layers and survival probabilities.

The stage calculates first the depth of a water column that is mass-equivalent to the path traversed by the neutrino through the earth. This is done using the same Layers module that is also used for oscillation. The survival probability is then calculated from the average cross-section with protons and neutrons.

pisa.stages.absorption.earth_absorption.FLAV_BAR_STR_MAPPING = {(0, -1): 'e_bar', (0, 1): 'e', (1, -1): 'mu_bar', (1, 1): 'mu', (2, -1): 'tau_bar', (2, 1): 'tau'}

Mapping from flav and nubar container content to the string for this neutrino in the ROOT file.

pisa.stages.absorption.earth_absorption.calculate_survivalprob(int_rho, xsection, out)[source]

Calculate survival probability given layer distances, layer densities and (pre-computed) cross-sections.

Parameters:

  • int_rho (scalar) – depth of mass equivalent water column in cm

  • xsection (scalar) – cross-section per nucleon in cm^2

  • out (scalar) – Result is stored here

class pisa.stages.absorption.earth_absorption.earth_absorption(earth_model, xsec_file, detector_depth=<Quantity(2.0, 'kilometer')>, prop_height=<Quantity(20.0, 'kilometer')>, **std_kwargs)[source]

Bases: Stage

earth absorption PISA Pi class

Parameters:

  • earth_model (str) – PREM file path

  • xsec_file (str) – path to ROOT file containing cross-sections

  • detector_depth (quantity (distance), optional) – detector depth

  • prop_height (quantity (distance), optional) – height of neutrino production in the atmosphere

apply_function()[source]

Implement in services (subclasses of Stage)

calculate_xsections(flav, nubar, energy)[source]

Calculates the cross-sections on isoscalar targets. The result is returned in cm^2. The xsection on one target is calculated by taking the xsection for O16 and dividing it by 16.

compute_function()[source]

Implement in services (subclasses of Stage)

setup_function()[source]

Implement in services (subclasses of Stage)

pisa.stages.absorption.earth_absorption.init_test(**param_kwargs)[source]

Instantiation example

Module contents