-
Notifications
You must be signed in to change notification settings - Fork 3
/
Copy pathexample_input.ini
56 lines (33 loc) · 2.06 KB
/
example_input.ini
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
//This file specifies the parameters that RelicFast takes, and gives some fiducial values.
//Make sure that if you add new stuff the code is consistent, and that the first instance of every variable name is right before the equal sign.
//Halo parameters
N_Mhalo = 1 ;//number of halo masses, logspaced. If 1 it chooses the minimum
Mhalo_min = 1e13;//minimum halo mass, in Msun
Mhalo_max = 1e15;//maximum halo mass, in Msun
//the code is not really designed to have more than one z_collapse, but it can support it.
z_collapse_bot = 0.7;//minimum redshift of collapse on the list
z_collapse_top = 1.4;//maximum redshift of collapse on the list
N_zcoll = 1;//number of redshifts of collapse, linearly spaced. If 1 it chooses the minimum
kbot= 1e-4;//minimum k_long in Mpc-1
ktop= 0.7;//maximum k_long in Mpc-1
N_klong = 4;// number of k_longs we take. If 1 it chooses the minimmum.
//if thermal relic:
tag_thermal_relic = 0; //whether to include a thermal relic
m_TR = .35 ; //mass of TR in eV
omega_TR = 0.0014; //physical TR energy density: Omega h^2.
//if regular sterile neutrino
tag_sterile_nu = 0; //whether to include sterile neutrino. Also used for a third neutrino.
m_SN = 0.02; //SN mass, in eV
//Fiducial Cosmological parameters (LCDM, within Planck 2015 errorbars)
omegab=0.022; //baryon physical density
omegac=0.12; //DM physical density
hubble=0.67; //h = H0/100. Reduced Hubble parameter.
// primordial curvature power spectrum
A_s=2.2e-9;//amplitude of perturbations
n_s=0.964; //scalar tilt
//Additional
mnu1=0.0; //mass of neutrino species 1, in eV (*if you choose a neutrino to be massive it should NOT be relativistic today, i.e., mnu_i>~0.01 eV)
mnu2=0.0; //mass of neutrino species 2, in eV
N_eff_input = 3.05; //what is the Neff that you want.
;//It will subtract 1 per each massive neutrino (including SN, but not TR).
;//If you modify it beyond normal LCDM values (or if you add a large N_eff in relics) modify explanatory_base.ini in the CLASS folder to have YHe=0.24.