Welcome to ulysses’s documentation!

Indices and tables

• Index
• Module Index
• Search Page

ULSBase

This is the base class.

class ulysses.ULSBase(*args, **kwargs)

__init__(*args, **kwargs)

Base class for all EtaB calculators.

Set global constants here.

Keyword Arguments:  

• vev (float) – Higgs VEV in GeV
• mhiggs (float) – Higgs mass in GeV
• mz (float) – Z-boson mass in GeV
• gstar (float) – Relativistic degrees of freedom at high temperature
• mplanck (float) – Planck mass in GeV
• mstar (float) – Neutrino cosmological mass in GeV

__call__(x)

Operator that returns EtaB for a given parameter point.

Arguments:

• x (dict) – parameter dictionary

NOTE — this operator is intended to be used with derived classes where EtaB is implemented

D1(k1, z)

Decay term for Boltzmann equation

D2(k, z)

Decay term for Boltzmann equation with two decaying steriles.

D3(k, z)

Decay term for Boltzmann equation with three decaying steriles.

DM

Diagonal heavy mass matrix.

DS(k1, z)

Function that replaces decay term when scattering is also present

Gamma1

Decay rate of N1.

Gamma2

Decay rate of N2.

Gamma3

Decay rate of N3.

N1Eq(z)

Equilibrium number density with one decaying sterile.

N2Eq(z)

Equilibrium number density with two decaying steriles.

N3Eq(z)

Equilibrium number density with three decaying steriles.

R

Orthogonal matrix R = R1.R2.R3

SqrtDM

Square root of diagonal heavy mass matrix. TODO: is this np.sqrt(self.DM) ???

SqrtDm

Square root of diagonal light mass matrix. Everything is in GeV.

U

PMNS matrix using the PDG parametrisation convention.

W1(k1, z)

Washout parameter with one decaying sterile.

W2(k, z)

Washout parameter with two decaying steriles.

W3(k, z)

Washout parameter with three decaying steriles.

c1a(a)

Probability coefficient for 1 a

c2a(a)

Probability coefficient for 2 a

c3a(a)

Probability coefficient for 3 a

epsilon(i, j, k, m)

CP asymmetry parameter. i,j,k,m denote indices in the heavy neutrino mass matrix.

epsilon1ab(a, b)

Off-diagonal CP asymmetry parameter for decays of N1. a and b denote lepton flavour.

epsilon2ab(a, b)

Off-diagonal CP asymmetry parameter for decays of N2. a and b denote lepton flavour.

epsilon3ab(a, b)

Off-diagonal CP asymmetry parameter for decays of N3. a and b denote lepton flavour.

epsilonaaRES(a)

CP asymmetry for resonant Leptogenesis.

epsiloniaaRES(a, i, j)

CP asymmetry for resonant leptogenesis in terms of i and j.

evolData

Getter:Return a 4-D array of the evolution data.

The first column is the evolution variable The second column corresponds to Ntautau, the third to Nmumu and the last columnd to Nee

f1(x)

Loop function in epsilon.

f2(x)

Loop function in epsilon.

fMLoop(x)

The loop function.

fMLoopHelper

Helper function.

fMR

This function returns the diagonl heavy neutrino mass matrix taking radiative corrections into account. (see h_loop) I.e. equivalent to self.SqrtDM in loop case

h

YUKAWA matrix.

h_loop

Yukawa matrix (LOOP + Tree).

h_tree

Yukawa matrix, tree-level.

hterm(a, b)

Projection probability projecting onto certain directions of flavour space indicated by the indices a and b.

a … [0,1,2]

0 = e 1 = mu 2 = tau

b … [0,1]

0 = term1 1 = term2

isPerturbative

Check for perturbative nature of Yukawas

k1

Decay parameter 1

k2

Decay parameter 2

k3

Decay parameter 3

m_loop

One-loop-level mass matrix.

m_tree

Tree-level mass matrix.

meff1

Effective mass 1 used for decay ans washout.

meff2

Effective mass 2 used for decay ans washout.

meff3

Effective mass 3 used for decay and washout.

printParams()

Print current parameters.

resonance(z)

calculate decay rate Gamma in terms of the total epsilon (epstot)

scat(z)

Function that multiplies washouts to incorporate scatterings.

setParams(pdict)

This set the model parameters. pdict is expected to be a dictionary

ulysses.tools

ulysses.tools.loadPlugin(model, **kwargs)

Plugin loader. This attemts to load and return and instance of the class CLASS in the file FILENAME. Both are given as a single string using : as separator. The kwargs are passed on to the base class.

Arguments:

• model (str) – FILENAME:CLASS

ulysses.tools.selectModel(model, **kwargs)

This function loads and returns an instance of one of the built-in models. The kwargs are passed on to the base class.

Arguments:

• model (str) – The model. Must be in [“1DME”, “2DME”, “3DME”, “1BE”, “2BE”, “2resonant”, “3DMEsct”, “3DMErhtau”]

Built-in models

class ulysses.EtaB_1DME(*args, **kwargs)

Density matrix equation (DME) with one decaying sterile. See arxiv:1112.4528.

RHS(y0, z, epstt, epsmm, epsee, epstm, epste, epsme, c1t, c1m, c1e, k)

EtaB()

class ulysses.EtaB_2DME(*args, **kwargs)

Density matrix equation (DME) with two decaying steriles. See arxiv:1112.4528.

RHS(y0, zzz, ETA, C, K, W)

EtaB()

class ulysses.EtaB_3DME(*args, **kwargs)

Density matrix equation (DME) with three decaying steriles. See arxiv:1112.4528.

RHS(y0, zzz, ETA, C, K, W)

EtaB()

class ulysses.EtaB_2Resonant(*args, **kwargs)

Resonant equations with two steriles and three lepton flavours. See arxiv:0705.2183.

RHS(y0, zzz, ETA, C, K)

EtaB()

class ulysses.EtaB_3DME_Scattering(*args, **kwargs)

Density matrix equation (DME) with three decaying steriles including scatterings. See hep-ph/0401240.

RHS(y0, zzz, ETA, C, K, W)

EtaB()

class ulysses.EtaB_3DS_Scattering_RHtaur(*args, **kwargs)

Density matrix equation (DME) with three decaying steriles including scatterings and effect of right-handed tau participation in kinetic equations. See hep-ph/0401240 and arxiv:0807.1408.

RHS(y0, zzz, ETA, C, K, W)

EtaB()