NuRadioMC.utilities.fluxes module

NuRadioMC.utilities.fluxes.get_limit_from_aeff(energy, aeff, livetime, signalEff=1.0, energyBinsPerDecade=1.0, upperLimOnEvents=2.44)[source]

Limit from effective volume

Parameters:
energy: array of floats

neutrino energy

veff: array of floats

effective volumes

livetime: float

time used

signalEff: float

efficiency of signal reconstruction

energyBinsPerDecade: float

1 for decade bins, 2 for half-decade bins, etc.

upperLimOnEvents: float

2.3 for Neyman UL w/ 0 background, 2.44 for F-C UL w/ 0 background, etc

NuRadioMC.utilities.fluxes.get_limit_flux(energy, veff_sr, livetime, signalEff=1.0, energyBinsPerDecade=1.0, upperLimOnEvents=2.44, nuCrsScn='ctw', inttype='total')[source]

Limit from effective volume

Parameters:
energy: array of floats

neutrino energy

veff_sr: array of floats

effective volume x solid angle

livetime: float

time used

signalEff: float

efficiency of signal reconstruction

energyBinsPerDecade: float

1 for decade bins, 2 for half-decade bins, etc.

upperLimOnEvents: float

2.3 for Neyman UL w/ 0 background, 2.44 for F-C UL w/ 0 background, etc

nuCrsScn: str

type of neutrino cross-section

NuRadioMC.utilities.fluxes.get_limit_e1_flux(energy, veff_sr, livetime, signalEff=1.0, energyBinsPerDecade=1.0, upperLimOnEvents=2.44, nuCrsScn='ctw', inttype='total')[source]

Limit from effective volume on E^1 flux plot

Parameters:
energy: array of floats

neutrino energy

veff_sr: array of floats

effective volume x solid angle

livetime: float

time used

signalEff: float

efficiency of signal reconstruction

energyBinsPerDecade: float

1 for decade bins, 2 for half-decade bins, etc.

upperLimOnEvents: float

2.3 for Neyman UL w/ 0 background, 2.44 for F-C UL w/ 0 background, etc

nuCrsScn: str

type of neutrino cross-section

NuRadioMC.utilities.fluxes.get_limit_e2_flux(energy, veff_sr, livetime, signalEff=1.0, energyBinsPerDecade=1.0, upperLimOnEvents=2.44, nuCrsScn='ctw', inttype='total')[source]

Limit from effective volume on E^2 flux plot

Parameters:
energy: array of floats

neutrino energy

veff_sr: array of floats

effective volumes x solid angle

livetime: float

time used

signalEff: float

efficiency of signal reconstruction

energyBinsPerDecade: float

1 for decade bins, 2 for half-decade bins, etc.

upperLimOnEvents: float

2.3 for Neyman UL w/ 0 background, 2.44 for F-C UL w/ 0 background, etc

nuCrsScn: str

type of neutrino cross-section

NuRadioMC.utilities.fluxes.get_number_of_events_for_flux(energies, flux, Veff, livetime, nuCrsScn='ctw')[source]

calculates the number of expected neutrinos for a certain flux assumption

Parameters:
energies: array of floats

energies (the bin centers), the binning in log10(E) must be equidistant!

flux: array of floats

the flux at energy logE

Veff: array of floats

the effective volume per energy logE

livetime: float

the livetime of the detector (including signal efficiency)

Returns:
array of floats: number of events per energy bin
NuRadioMC.utilities.fluxes.get_exposure(energy, Veff, field_of_view=6.283185307179586)[source]

calculate exposure from effective volume

Parameters:
energy: float

neutrino energy (needed to calculate interaction length)

Veff: float

effective volume

field_of_view: float

the field of view of the detector

Returns:
float: exposure
NuRadioMC.utilities.fluxes.get_integrated_exposure(exp_func, E_low, E_high)[source]

calculates the integral int(E^-2 * exposure(E) dE)

integration is performed in log space

NuRadioMC.utilities.fluxes.get_fluence_limit(int_exp)[source]

calculates the fluence limit for a integrated exposure