Project bayesian-priors

Invoke with: import icecube.bayesian_priors

C++ ServiceFactorys

I3ConstantZenithWeightServiceFactory

"I3ConstantZenithWeightServiceFactory" (C++ ServiceFactory)

<undocumented>

Param CosZenithRange:

Default = [-1.0, 1.0], Minimum and maximum value of cos(zenith) (before any flipping).

Param FlipTrack:

Default = False, Boolean option: if you set this to true, then the weight will be computed as if the track pointed in the opposite direction.

Param PenaltySlope:

Default = 0.0, If your goal is to force a reconstruction to end up within within a specific zenith range, but you don’t like setting bounds on the parameters, then you can stimulate the reconstruction to ‘roll’ back to the allowed range via a sloped penalty value, i.e. a penalty which gets worse farther away from the allowed zenith range, then you should set a negitive value here, e.g. -1000.

Param PenaltyValue:

Default = -200.0, LogLikelihood value outside the cos(zenith) range. Should be a large negative number.

I3PolynomialZenithWeightServiceFactory

"I3PolynomialZenithWeightServiceFactory" (C++ ServiceFactory)

<undocumented>

Param Coefficients:

Default = [-0.007355395666666667, 0.1899537877777778, -1.5003446333333332, 1.6376264333333335, 23.163495444444443, -61.62309222222222, 60.70669777777778, -21.576655], Coefficients for a polynomial of cos(zenith). The default uses the same coefficients as Ty’s P7 zenith weight function for AMANDA data.

Param CosZenithRange:

Default = [-1.0, 1.0], Minimum and maximum value of cos(zenith) (before any flipping).

Param DefaultWeight:

Default = 1.111111111111111e-06, Default weight value for cos(zen) < MinCosZenith

Param FlipTrack:

Default = False, Boolean option: if you set this to true, then the weight will be computed as if the track pointed in the opposite direction.

Param MinCosZenith:

Default = 0.1, Minimum cos(zenith) value; for lower values the default weight is used.

Param PenaltySlope:

Default = 0.0, If your goal is to force a reconstruction to end up within within a specific zenith range, but you don’t like setting bounds on the parameters, then you can stimulate the reconstruction to ‘roll’ back to the allowed range via a sloped penalty value, i.e. a penalty which gets worse farther away from the allowed zenith range, then you should set a negitive value here, e.g. -1000.

Param PenaltyValue:

Default = -200.0, LogLikelihood value outside the cos(zenith) range. Should be a large negative number.

I3PowExpZenithWeightServiceFactory

"I3PowExpZenithWeightServiceFactory" (C++ ServiceFactory)

This service provides a gulliver I3EventLogLikelihoodBase which includes a prior probability to the liklihood as a function of zenith angle to provide a “Bayesian” reconstruction.

Param Amplitude:

Default = 2.49655e-07, Amplitude factor (a0 in w=a0*pow(z,a1)*exp(-a2/z) with z=cos(zenith))

Param CosZenithRange:

Default = [-1.0, 1.0], Minimum and maximum value of cos(zenith) (before any flipping).

Param DefaultWeight:

Default = 1.3838965267367376e-87, Default weight for to low cos(zenith) values

Param ExponentFactor:

Default = 0.778393, Numerator in the exponent of the exponent factor (a2 in w=a0*pow(z,a1)*exp(-a2/z) with z=cos(zenith))

Param FlipTrack:

Default = False, Boolean option: if you set this to true, then the weight will be computed as if the track pointed in the opposite direction.

Param PenaltySlope:

Default = 0.0, If your goal is to force a reconstruction to end up within within a specific zenith range, but you don’t like setting bounds on the parameters, then you can stimulate the reconstruction to ‘roll’ back to the allowed range via a sloped penalty value, i.e. a penalty which gets worse farther away from the allowed zenith range, then you should set a negitive value here, e.g. -1000.

Param PenaltyValue:

Default = -200.0, LogLikelihood value outside the cos(zenith) range. Should be a large negative number.

Param Power:

Default = 1.67721, Exponent in the power factor (a1 in w=a0*pow(z,a1)*exp(-a2/z) with z=cos(zenith))

I3TableZenithWeightServiceFactory

"I3TableZenithWeightServiceFactory" (C++ ServiceFactory)

<undocumented>

Param CosZenithRange:

Default = [-1.0, 1.0], Minimum and maximum value of cos(zenith) (before any flipping).

Param FlipTrack:

Default = False, Boolean option: if you set this to true, then the weight will be computed as if the track pointed in the opposite direction.

Param MinCosZenith:

Default = 0.1, Minimum cos(zenith) value; for lower values the default weight is used.

Param PenaltySlope:

Default = 0.0, If your goal is to force a reconstruction to end up within within a specific zenith range, but you don’t like setting bounds on the parameters, then you can stimulate the reconstruction to ‘roll’ back to the allowed range via a sloped penalty value, i.e. a penalty which gets worse farther away from the allowed zenith range, then you should set a negitive value here, e.g. -1000.

Param PenaltyValue:

Default = -200.0, LogLikelihood value outside the cos(zenith) range. Should be a large negative number.

Param TableValues:

Default = [1.0, 200.0, 5000.0, 50000.0, 100000.0, 200000.0, 400000.0, 500000.0, 600000.0, 800000.0, 900000.0], Table values, for equidistant values of cos(zenith). The table should contain at least 2 values, the first the first gives the weight at the minimum cos(zenith) value, the first gives the weight at cos(zenith)==1. The default uses values originally optimized for AMANDA data.