From be9255a0d769e0d06b1d72c66b284f42ba859b00 Mon Sep 17 00:00:00 2001 From: JCGoran <goran.jelic-cizmek@unige.ch> Date: Thu, 28 Jun 2018 20:38:38 +0200 Subject: [PATCH] Added non-source files --- PkL_CLASS.dat | 159 +++++++++++++++++++++++++++++++++++ separations.dat | 184 +++++++++++++++++++++++++++++++++++++++++ settings.cfg | 215 ++++++++++++++++++++++++++++++++++++++++++++++++ 3 files changed, 558 insertions(+) create mode 100644 PkL_CLASS.dat create mode 100644 separations.dat create mode 100644 settings.cfg diff --git a/PkL_CLASS.dat b/PkL_CLASS.dat new file mode 100644 index 0000000..d77268d --- /dev/null +++ b/PkL_CLASS.dat @@ -0,0 +1,159 @@ + 1.046002468616e-07 6.235364087170e-01 + 1.316839088540e-07 7.781293869189e-01 + 1.657802191807e-07 9.710481518263e-01 + 2.087049306994e-07 1.211799689987e+00 + 2.627439408242e-07 1.512240402026e+00 + 3.307750238985e-07 1.887186149825e+00 + 4.164210831727e-07 2.355091179122e+00 + 5.242430836130e-07 2.938987403334e+00 + 6.599829399177e-07 3.667640613187e+00 + 8.308692944130e-07 4.576937720508e+00 + 1.046002468616e-06 5.711684797774e+00 + 1.316839088540e-06 7.127797422125e+00 + 1.657802191807e-06 8.894993767927e+00 + 2.087049306994e-06 1.110039754482e+01 + 2.627439408242e-06 1.385238500634e+01 + 3.307750238985e-06 1.728689331607e+01 + 4.164210831727e-06 2.157269994273e+01 + 5.242430836129e-06 2.692138005360e+01 + 6.599829399177e-06 3.359570573321e+01 + 8.308692944130e-06 4.192510972798e+01 + 1.046002468616e-05 5.231953251804e+01 + 1.316839088540e-05 6.529098007757e+01 + 1.657802191807e-05 8.147833032282e+01 + 2.087049306994e-05 1.016787939908e+02 + 2.627439408242e-05 1.268871368768e+02 + 3.307750238985e-05 1.583445906013e+02 + 4.164210831727e-05 1.975997213788e+02 + 5.242430836130e-05 2.465842771463e+02 + 6.599829399177e-05 3.077074737843e+02 + 8.308692944130e-05 3.839729381223e+02 + 1.046002468616e-04 4.791232322596e+02 + 1.316839088540e-04 5.978175666464e+02 + 1.657802191807e-04 7.458483460835e+02 + 2.087049306994e-04 9.304015200118e+02 + 2.627439408242e-04 1.160363097424e+03 + 3.307750238985e-04 1.446668624527e+03 + 4.164210831727e-04 1.802682995344e+03 + 5.242430836130e-04 2.244585193304e+03 + 6.599829399177e-04 2.791718483337e+03 + 8.308692944130e-04 3.466843869212e+03 + 1.046002468616e-03 4.296160591409e+03 + 1.316839088540e-03 5.308791813948e+03 + 1.657802191807e-03 6.535188754858e+03 + 2.087049306994e-03 8.004803666436e+03 + 2.627439408242e-03 9.741203518106e+03 + 3.307750238985e-03 1.175527796161e+04 + 4.164210831727e-03 1.403598949228e+04 + 5.242430836130e-03 1.653770957442e+04 + 6.599829399177e-03 1.916470680453e+04 + 8.308692944130e-03 2.175693061375e+04 + 1.046002468616e-02 2.407037967204e+04 + 1.316839088540e-02 2.576916386540e+04 + 1.657802191804e-02 2.644066661190e+04 + 2.087049290115e-02 2.568269081177e+04 + 2.627426607523e-02 2.330574353176e+04 + 3.305983861456e-02 1.967028021559e+04 + 4.105589834407e-02 1.602355992606e+04 + 4.779681971948e-02 1.392626514416e+04 + 5.268562702534e-02 1.288880123643e+04 + 5.671721594501e-02 1.223759727526e+04 + 6.031687829181e-02 1.174365662217e+04 + 6.367418442744e-02 1.131141199246e+04 + 6.688946261170e-02 1.089205686404e+04 + 7.002286650976e-02 1.046031012785e+04 + 7.311390692949e-02 1.000446612186e+04 + 7.619042858984e-02 9.520258933464e+03 + 7.927320502998e-02 9.012490350506e+03 + 8.237849179503e-02 8.490012576078e+03 + 8.551953925527e-02 7.965379106689e+03 + 8.870753630877e-02 7.452866103636e+03 + 9.195222427013e-02 6.966656922411e+03 + 9.526231024366e-02 6.519762304196e+03 + 9.864575353763e-02 6.121786086448e+03 + 1.021099688283e-01 5.778453921247e+03 + 1.056619730371e-01 5.491472744261e+03 + 1.093084931093e-01 5.256628633898e+03 + 1.130560459802e-01 5.065780809561e+03 + 1.169109983473e-01 4.906439520419e+03 + 1.208796115317e-01 4.763158483537e+03 + 1.249680751995e-01 4.618879035664e+03 + 1.291825327147e-01 4.458647792313e+03 + 1.335291002364e-01 4.271387181131e+03 + 1.380138812293e-01 4.053559592995e+03 + 1.426429777658e-01 3.809515635752e+03 + 1.474224998042e-01 3.551298500134e+03 + 1.523585735031e-01 3.295395783048e+03 + 1.574573495511e-01 3.059396256620e+03 + 1.627250124421e-01 2.856661250781e+03 + 1.681677915969e-01 2.693462239672e+03 + 1.737919752073e-01 2.567549830606e+03 + 1.796039276714e-01 2.468144036854e+03 + 1.856101114764e-01 2.379834974201e+03 + 1.918171143817e-01 2.285726308966e+03 + 1.982316827589e-01 2.172737501207e+03 + 2.048607619578e-01 2.039108418104e+03 + 2.117115446000e-01 1.894369971192e+03 + 2.187915277619e-01 1.753620388901e+03 + 2.261085801025e-01 1.630316189727e+03 + 2.336710201441e-01 1.531112366691e+03 + 2.414877071275e-01 1.453842042334e+03 + 2.495681461696e-01 1.387905170938e+03 + 2.579226098714e-01 1.320134555909e+03 + 2.665622790882e-01 1.242906350006e+03 + 2.754994063270e-01 1.157703473246e+03 + 2.847475062343e-01 1.073446906898e+03 + 2.943215789532e-01 1.000161847404e+03 + 3.042383738781e-01 9.404478634107e+02 + 3.145167036585e-01 8.890953092738e+02 + 3.251778214267e-01 8.382522377623e+02 + 3.362458784590e-01 7.834337869853e+02 + 3.477484852899e-01 7.277216952219e+02 + 3.597174073793e-01 6.767921568778e+02 + 3.721894378311e-01 6.330801940196e+02 + 3.852075059997e-01 5.932430307158e+02 + 3.988221046420e-01 5.532692664077e+02 + 4.130931536463e-01 5.130955397824e+02 + 4.280924719868e-01 4.756274674146e+02 + 4.439071126178e-01 4.417256922022e+02 + 4.606439469169e-01 4.095987886090e+02 + 4.784361004430e-01 3.778085203090e+02 + 4.974522035406e-01 3.475017783359e+02 + 5.179100495760e-01 3.193258926116e+02 + 5.400973908098e-01 2.919260169042e+02 + 5.644047482072e-01 2.653798353486e+02 + 5.913793669410e-01 2.399835927855e+02 + 6.218183726548e-01 2.150674574747e+02 + 6.569391134310e-01 1.906853493852e+02 + 6.987123569887e-01 1.663926167344e+02 + 7.505660696221e-01 1.418939207055e+02 + 8.189917989690e-01 1.166537176031e+02 + 9.173402034529e-01 9.018541068760e+01 + 1.072860309668e+00 6.291194875862e+01 + 1.320414800762e+00 3.871991604207e+01 + 1.659546108755e+00 2.247636909241e+01 + 2.089192449348e+00 1.287803244810e+01 + 2.630137270126e+00 7.317815913718e+00 + 3.311146645767e+00 4.127008815959e+00 + 4.168486654533e+00 2.311755560598e+00 + 5.247813778117e+00 1.287002807355e+00 + 6.606606121635e+00 7.124349059768e-01 + 8.317224332241e+00 3.923399364081e-01 + 1.047076506745e+01 2.150334916342e-01 + 1.318191222434e+01 1.173322685326e-01 + 1.659504427526e+01 6.375816692002e-02 + 2.089192294798e+01 3.451098286821e-02 + 2.630137270045e+01 1.861093164120e-02 + 3.311146645767e+01 9.999853410280e-03 + 4.168486654533e+01 5.353185944966e-03 + 5.247813778117e+01 2.854539593210e-03 + 6.606606121635e+01 1.515512297386e-03 + 8.317224332241e+01 8.004430803193e-04 + 1.047076506745e+02 4.199626525805e-04 + 1.318191222434e+02 2.184026684970e-04 + 1.659504427526e+02 1.122067043641e-04 + 2.089192294798e+02 5.668850302366e-05 + 2.630137270045e+02 2.800241111863e-05 + 3.311146645767e+02 1.344544980356e-05 + 4.168486654533e+02 6.260444081632e-06 + 5.247813778117e+02 2.845605974333e-06 diff --git a/separations.dat b/separations.dat new file mode 100644 index 0000000..efc588d --- /dev/null +++ b/separations.dat @@ -0,0 +1,184 @@ +0.1 +0.2 +0.4 +0.8 +1. +1.5 +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. +32. +34. +36. +38. +40. +42. +44. +46. +48. +50. +52. +54. +56. +58. +60. +62. +64. +66. +68. +70. +72. +74. +76. +78. +80. +81. +82. +83. +84. +85. +86. +87. +88. +89. +90. +91. +92. +93. +94. +95. +95.5 +96. +96.5 +97. +97.5 +98. +98.5 +99. +99.5 +100. +100.5 +101. +101.5 +102. +102.5 +103. +103.5 +104. +104.5 +105. +106. +107. +108. +109. +110. +112. +114. +116. +118. +120. +124. +128. +132. +136. +140. +144. +148. +152. +156. +160. +164. +168. +172. +176. +180. +185. +190. +195. +200. +205. +210. +215. +220. +225. +230. +235. +240. +250. +260. +270. +280. +290. +300. +320. +340. +360. +380. +400. +420. +440. +460. +480. +500. +520. +540. +560. +580. +600. +620. +640. +660. +680. +700. +750. +800 +850 +900 +950 +1000 +1050 +1100 +1150 +1200 +1250 +1300 +1350 +1400 +1450 +1500 +1550 +1600 +1650 +1700 +1750 +1800 +1850 +1900 +1950 +2000 diff --git a/settings.cfg b/settings.cfg new file mode 100644 index 0000000..0a75858 --- /dev/null +++ b/settings.cfg @@ -0,0 +1,215 @@ +#*~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*# +#* COFF-E parameter file *# +#*~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*# + + +############### +#(1): Input # +############### + +### (1.a) +# file containing the separations for which the +# correlation function is computed (not needed if "output_type" is 0, 4, 5, 6) +# NOTE: separations must be in units Mpc/h + +input_separations = "separations.dat"; + +### (1.b) +# file containing the power spectrum +# NOTE: file containing spectrum must have two columns; +# allowed separators are ' ', ',', '\t' (tabs), or ':'; +# k must be in h/Mpc, and P(k) must be in (Mpc/h)^3 + +input_power_spectrum = "PkL_CLASS.dat"; + +### (1.c) +# matter density parameter and radiation density parameter +# NOTE: dark energy density parameter is computed from +# the closure relation sum(Omega) + Omega(DE) = 1 + +omega_m = 0.3128; +omega_gamma = 0.0000916; + + +### (1.d) +# equation of state for dark energy, parametrised as w(z) = w0 + wa z/(1+z) + +w0 = -1.0; +wa = 0.0; + +### (1.e) +# galaxy bias, magnification bias and evolution bias +# either a constant value or read a file containing the bias as a function of redshift (in ascending redshift) +# "1" and "2" label the galaxies populations (with only one population only "1" is read) + +matter_bias1 = 1.; +read_matter_bias1 = 0; +input_matter_bias1 = ""; + +matter_bias2 = 1.; +read_matter_bias2 = 0; +input_matter_bias2 = ""; + +magnification_bias1 = 0.; +read_magnification_bias1 = 0; +input_magnification_bias1 = ""; + +magnification_bias2 = 0.; +read_magnification_bias2 = 0; +input_magnification_bias2 = ""; + +evolution_bias1 = 0.; +read_evolution_bias1 = 0; +input_evolution_bias1 = ""; + +evolution_bias2 = 0.; +read_evolution_bias2 = 0; +input_evolution_bias2 = ""; + +### (1.f) +# parameter for the covariance +# respectively: the pixel size (in Mpc/h), the mean number density at z_mean (in (h/Mpc)^3) and the sky coverage of the catalog + +covariance_z_mean = [0.5, 1.0, 1.5]; +covariance_deltaz = [0.1, 0.1, 0.2]; +covariance_density = [0.1, 0.2, 0.01]; +covariance_fsky = [0.3, 0.5, 0.2]; +covariance_pixelsize = 10.0; + +covariance_zmin = [2.0, 2.2, 2.3]; +covariance_zmax = [2.5, 2.8, 2.5]; + +############### +#(2): Output # +############### + +### (2.a) +# the relative or absolute path of output; if directory doesn't exist, COFFE will try to create it +# a prefix can be appended to the standard outputs; if it's "$TIME", then it appends the timestamp +# in the format "MM-DD-HH-mm_". Note that COFFE may overwrite previous results, it is up to the user +# to make sure the path doesn't already exist! + +output_path = "results/"; + +output_prefix = "$TIME"; + +### (2.b) +# which projection effect to take into account (see 1708.00492 for details), possible values are: +# rsd = redshift space distortion +# den = density +# d1, d2 = nonintegrated doppler terms +# g1, g2, g3 = nonintegrated gravitational potential terms +# g4, g5 = integrated gravitational potential terms +# len = integrated lensing term +# NOTE: the correlation calculated for types A and B is automatically A*A + A*B + B*A + B*B + +correlation_contributions = ["den", "rsd"]; + +### (2.c) +# for output - specify if you want: +# (see the User guide for details) +# (0) angular correlation function +# (1) correlation function as function of mean redshift, angles, and separation +# (2) multipoles as function of l, mean redshift, and separation +# (3) average multipoles as function of l and separation +# (4) covariance of multipoles +# (5) covariance of redshift averaged multipoles +# (6) 2D correlation function as a function of r_parallel and r_perpendicular + +output_type = 2; + +### (2.d) +# if output_type is not 3, the mean redshift for which to calculate the output + +z_mean = 1.0; + +### (2.e) +# the thickness of the redshift bin for which we calculate +# the output; cannot be negative or larger than z_mean for obvious reasons + +deltaz = 0.5; + +z_min = 0.7; +z_max = 1.3; + +### (2.f) +# needed if output_type = 1 + +mu = [0.7]; + +### (2.g) +# needed if output_type = 2, 3, 4 or 5 +# multipoles you want to compute: 0=monopole, 2=quadrupole, 4=hexadecapole, 6=tetrahexadecapole, etc. + +multipoles = [0, 2, 4]; + +### (2.h) +# optional: which background quantities to output; possible values are: +# z = redshift +# a = scale factor +# H = Hubble rate +# conformal_H = conformal Hubble rate +# conformal_H_prime = derivative of conformal Hubble rate w.r.t. conformal time +# D1 = growth function +# f = derivative d(ln D1)/d(ln a) +# comoving_distance +# NOTE: they will be written into the file in the same order as given here + +output_background = ["z", "a", "H", "conformal_H", "conformal_H_prime", "D1", "f", "comoving_distance"]; + +########################### +#(3): Precision settings # +########################### + +### (3.a) +# the sampling rate for the background; about 1 second to sample 10000 points + +background_sampling = 10000; + +### (3.b) +# for how many points to compute the integral of P(k) k^2 j_l(kr) (NOTE: runtime is <1 s for 10000 points) + +bessel_sampling = 10000; + +### (3.c) +# the sampling for the angular correlation function (between 0 and pi/2) + +theta_sampling = 3000; + +### (3.d) +# double integrated terms are computed using monte carlo methods from GSL; +# the available methods are: +# 0 - standard random sampling +# 1 - MISER algorithm of Press and Farrar; based on recursive stratified sampling +# 2 - VEGAS algorithm of Lepage; based on importance sampling +# NOTE: if CUBA is used, only the integration_len parameter is needed +# reference: about 60000 for correlation function, +# 300000 for multipoles, more for redshift-averaged multipoles + +integration_method = 2; +integration_sampling = 750000; + +### (3.e) +# optional: the range of integration for the integral +# over the power spectrum +# NOTE: the range should be -smaller- than the one +# read from the file, otherwise the range from the file +# will be used instead +# NOTE: must be in h/Mpc + +k_min = 0.00001; +k_max = 300.; + +### (3.f) +# the interpolation type for all interpolation variables +# possible choices (taken from GSL manual 2.4): +# 1 - linear +# 2 - polynomial +# 3 - cubic spline with natural boundary conditions +# 4 - cubic spline with periodic boundary conditions +# 5 - non-rounded Akima spline with natural boundary conditions +# 6 - non-rounded Akima spline with periodic boundary conditions +# 7 - monotone cubic spline, described in: +# M. Steffen, A simple method for monotonic interpolation in one dimension, Astron. Astrophys. 239, 443-450, 1990. + +interpolation = 5; -- GitLab