Mercurial > repos > astroteam > sgwb_astro_tool

annotate Phase_transition_parameters.py @ 0:c9fc89ee996e draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
author astroteam
date Wed, 17 Apr 2024 10:29:23 +0000
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
1 #!/usr/bin/env python
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
2 # coding: utf-8
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
3
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
4 # flake8: noqa
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
5
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
6 import json
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
7 import os
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
8 import shutil
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
9
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
10 from oda_api.json import CustomJSONEncoder
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
11
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
12 get_ipython().run_cell_magic( # noqa: F821
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
13 "bash",
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
14 "",
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
15 "wget https://gitlab.renkulab.io/astronomy/mmoda/sgwb/-/raw/master/EPTA.csv\nwget https://gitlab.renkulab.io/astronomy/mmoda/sgwb/-/raw/master/NANOGrav23.csv\nwget https://gitlab.renkulab.io/astronomy/mmoda/sgwb/-/raw/master/PPTA.csv\nwget https://gitlab.renkulab.io/astronomy/mmoda/sgwb/-/raw/master/RoperPol22.csv\nwget https://gitlab.renkulab.io/astronomy/mmoda/sgwb/-/raw/master/MAGIC.csv\nwget https://gitlab.renkulab.io/astronomy/mmoda/sgwb/-/raw/master/Ellis.csv\nwget https://gitlab.renkulab.io/astronomy/mmoda/sgwb/-/raw/master/NANO_alpha_T.csv\nwget https://gitlab.renkulab.io/astronomy/mmoda/sgwb/-/raw/master/NANO_alpha_T1.csv\nwget https://gitlab.renkulab.io/astronomy/mmoda/sgwb/-/raw/master/NANO_bubble.csv\nwget https://gitlab.renkulab.io/astronomy/mmoda/sgwb/-/raw/master/NANO_sound.csv\n",
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
16 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
17
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
18 import matplotlib.pyplot as plt
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
19
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
20 # Loading necessary packages
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
21 import numpy as np
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
22
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
23 get_ipython().run_line_magic("matplotlib", "inline") # noqa: F821
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
24 import astropy.units as u
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
25 from astropy.constants import c
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
26 from numpy import exp, log, pi, sqrt
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
27 from oda_api.api import ProgressReporter
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
28 from oda_api.data_products import PictureProduct
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
29
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
30 # Parameters of the phase transition
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
31
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
32 # fraction of turbulent energy that goes to gw N.B. arXiv:1004.4187 claims that epsilon_turb=0.05, but checks below show that it is rather 0.01
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
33 epsilon_turb = 1.0 # http://odahub.io/ontology#Float
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
34
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
35 _galaxy_wd = os.getcwd()
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
36
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
37 with open("inputs.json", "r") as fd:
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
38 inp_dic = json.load(fd)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
39 if "_data_product" in inp_dic.keys():
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
40 inp_pdic = inp_dic["_data_product"]
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
41 else:
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
42 inp_pdic = inp_dic
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
43
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
44 for vn, vv in inp_pdic.items():
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
45 if vn != "_selector":
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
46 globals()[vn] = type(globals()[vn])(vv)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
47
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
48 Tmax_yrs = 10.0 # http://odahub.io/ontology#Float
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
49 Tmin_yrs = 2.0 # http://odahub.io/ontology#Float
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
50 # terminal velocity of bubbles
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
51 v_w = 0.999 # http://odahub.io/ontology#Float
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
52 h = 0.7 # http://odahub.io/ontology#Float
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
53 # Numbers of relativistic degrees of freedom
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
54 g_star = 20 # http://odahub.io/ontology#Integer
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
55
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
56 c_s = 3 ** (-0.5) # speed of sound
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
57 F0_GW = 1.64e-5 / h**2 * (100 / g_star) ** (1 / 3.0)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
58 pr = ProgressReporter()
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
59 pr.report_progress(stage="Progress", progress=1.0)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
60
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
61 # Eq.20 of arXiv:1512.06239, corrected according to Appendix A of arXiv:1004.4187
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
62
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
63 def ka_v(alpha_v, v):
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
64 zeta = ((2.0 / 3.0 * alpha_v + alpha_v**2) ** 0.5 + (1 / 3.0) ** 0.5) / (
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
65 1 + alpha_v
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
66 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
67 kappa_a = (
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
68 6.9
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
69 * v ** (6.0 / 5.0)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
70 * alpha_v
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
71 / (1.36 - 0.037 * alpha_v**0.5 + alpha_v)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
72 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
73 kappa_b = alpha_v ** (2.0 / 5.0) / (
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
74 0.017 + (0.997 + alpha_v) ** (2.0 / 5.0)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
75 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
76 kappa_c = alpha_v**0.5 / (0.135 + (0.98 + alpha_v) ** 0.5)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
77 kappa_d = alpha_v / (0.73 + 0.083 * alpha_v**0.6 + alpha_v)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
78 deltak = -0.9 * log(alpha_v**0.5 / (1 + alpha_v**0.5))
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
79 if v < c_s:
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
80 return (
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
81 c_s ** (11.0 / 5.0)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
82 * kappa_a
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
83 * kappa_b
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
84 / (
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
85 (c_s ** (11.0 / 5.0) - v ** (11.0 / 5.0)) * kappa_b
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
86 + v * c_s ** (6.0 / 5.0) * kappa_a
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
87 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
88 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
89 elif v > zeta:
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
90 return (
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
91 (zeta - 1) ** 3
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
92 * (zeta / v) ** (5.0 / 2)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
93 * kappa_c
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
94 * kappa_d
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
95 / (
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
96 ((zeta - 1) ** 3 - (v - 1) ** 3)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
97 * zeta ** (5.0 / 2.0)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
98 * kappa_c
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
99 + (v - 1) ** 3 * kappa_d
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
100 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
101 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
102 else:
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
103 return (
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
104 kappa_b
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
105 + (v - c_s) * deltak
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
106 + (v - c_s) ** 3
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
107 / (zeta - c_s) ** 3
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
108 * (kappa_c - kappa_b - (zeta - c_s) * deltak)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
109 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
110
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
111 T_star = 0.2
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
112
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
113 # Comoving Hubble rate at the phase transition Eq. 11 of arXiv:1512.06239
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
114 def H_star(T_star):
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
115 return (
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
116 16.5e-6
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
117 * (T_star / (100.0 * u.GeV))
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
118 * (g_star / 100) ** (1.0 / 6.0)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
119 * u.Hz
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
120 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
121
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
122 Hstar = H_star(T_star * u.GeV)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
123 logHstar = np.log10(Hstar / u.Hz)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
124
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
125 fmin = logHstar.value - 5
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
126 fmax = logHstar.value + 5
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
127 ff = np.logspace(fmin, fmax, 101)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
128
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
129 # HL model formula
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
130 def GW_sound(f, T_star, alpha, beta_H, v_w):
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
131 Hstar = H_star(T_star)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
132 kappa_v = ka_v(alpha, v_w)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
133 K = kappa_v * alpha / (1 + alpha)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
134 lambda_star = (8 * pi) ** (1 / 3) * max([v_w, c_s]) / (beta_H * Hstar) * c
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
135 Delta_w = sqrt((v_w - c_s) ** 2) / v_w
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
136 s2 = Delta_w * lambda_star * f / c
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
137 # print((c/lambda_star).cgs,(c/(Delta_w*lambda_star)).cgs)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
138 s1 = lambda_star * f / c
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
139 M = (
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
140 16
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
141 * (1 + Delta_w ** (-3)) ** (2 / 3.0)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
142 * (Delta_w * s1) ** 3
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
143 / ((1 + s1**3) ** (2.0 / 3.0) * (3 + s2**2) ** 2)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
144 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
145 factor = (K * lambda_star * Hstar) ** 2 / (
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
146 sqrt(K) + lambda_star * Hstar / c
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
147 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
148 mu = 4.78 - 6.27 * Delta_w + 3.34 * Delta_w**2
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
149 ff = f / u.Hz
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
150 dlnf = (ff[1] - ff[0]) / ff[0]
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
151 mu = sum(M) * dlnf
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
152 B = 1e-2 / mu
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
153 return (3 * B * factor / c**2 * F0_GW * M).cgs
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
154
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
155 # Eq 1 of the new Overleaf, from Alberto
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
156
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
157 def GW_turb_Andrii(f, T_star, alpha, beta_H, v_w, epsilon_turb):
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
158 Hstar = H_star(T_star)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
159
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
160 # Eq. 1 of 2307.10744
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
161 lambda_star = (
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
162 (8 * pi) ** (1 / 3) * max([v_w, c_s]) / (beta_H * Hstar)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
163 ) * c # characteristic light-crossing distance scale
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
164
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
165 kappa_v = ka_v(alpha, v_w)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
166 K = kappa_v * alpha / (1 + alpha)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
167
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
168 # Eq. 10 of 2307.10744
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
169 Omega_star = epsilon_turb * K
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
170
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
171 # Eq. 13 of 2307.10744 and text after it
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
172 u_star = sqrt(0.75 * Omega_star)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
173 dt_fin = 2 * lambda_star / u_star / c
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
174 s3 = dt_fin * f
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
175 s1 = f * lambda_star / c
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
176 # print('u_star=',u_star,'lambda_star=',lambda_star.cgs)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
177
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
178 # Eq. 15 of 2307.10744
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
179 T_GW = np.log(1 + Hstar * dt_fin / (2 * pi)) * (s3 < 1) + np.log(
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
180 1 + lambda_star * Hstar / c / (2 * pi * s1)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
181 ) * (s3 >= 1)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
182 T_GW = T_GW**2
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
183
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
184 # Eq. 17 of 2307.10744
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
185 alpha_pi = 2.15
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
186 s_pi = 2.2
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
187 P_pi = (1 + (s1 / s_pi) ** alpha_pi) ** (-11 / (3 * alpha_pi))
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
188
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
189 # Eq 18,19,20 of 2307.10744
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
190 A = 2e-3 * 1.4 * 0.6
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
191
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
192 # Eq. 14 of 2307.10744
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
193 Sturb = (
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
194 4
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
195 * pi**2
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
196 * s1**3
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
197 * T_GW
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
198 / (lambda_star * Hstar / c) ** 2
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
199 * P_pi
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
200 / 1.4
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
201 / 0.6
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
202 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
203
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
204 # Eq 9 of 2307.10744
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
205 res = (
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
206 3 * A * Omega_star**2 * (lambda_star * Hstar / c) ** 2 * F0_GW * Sturb
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
207 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
208
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
209 Omega_B = Omega_star / 2.0
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
210 Omega_gamma = 2 / g_star
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
211 B = 3e-6 * (Omega_B / Omega_gamma) ** 0.5
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
212 return res, lambda_star, B
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
213
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
214 H0 = 70 * (u.km / u.s) / u.Mpc
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
215
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
216 d = np.genfromtxt("NANOGrav23.csv")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
217 gammas_nano = d[:, 0]
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
218 As_nano = 10 ** d[:, 1]
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
219 ps_nano = 5 - gammas_nano
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
220
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
221 d = np.genfromtxt("EPTA.csv")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
222 gammas_epta = d[:, 0]
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
223 As_epta = 10 ** d[:, 1]
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
224 ps_epta = 5 - gammas_epta
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
225
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
226 d = np.genfromtxt("PPTA.csv")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
227 gammas_ppta = d[:, 0]
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
228 As_ppta = 10 ** d[:, 1]
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
229 ps_ppta = 5 - gammas_ppta
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
230
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
231 fref = (1 / u.yr).cgs.value
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
232 lgfmin = np.log10(fref / Tmax_yrs)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
233 lgfmax = np.log10(fref / Tmin_yrs)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
234 fff = np.logspace(lgfmin, lgfmax, 10) * u.Hz
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
235 min_nano = np.ones(len(fff))
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
236 max_nano = np.zeros(len(fff))
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
237 for i in range(len(As_nano)):
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
238 spec = (
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
239 2
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
240 * pi**2
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
241 / 3
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
242 / H0**2
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
243 * fff**2
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
244 * As_nano[i] ** 2
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
245 * (fff / fref) ** (3 - gammas_nano[i])
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
246 ).cgs.value
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
247 min_nano = np.minimum(spec, min_nano)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
248 max_nano = np.maximum(spec, max_nano)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
249 # plt.plot(ff,spec)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
250 plt.fill_between(
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
251 fff.value, min_nano, max_nano, color="red", alpha=0.5, label="NANOGrav"
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
252 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
253 min_epta = np.ones(len(fff))
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
254 max_epta = np.zeros(len(fff))
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
255 for i in range(len(As_epta)):
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
256 spec = (
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
257 2
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
258 * pi**2
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
259 / 3
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
260 / H0**2
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
261 * fff**2
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
262 * As_epta[i] ** 2
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
263 * (fff / fref) ** (3 - gammas_epta[i])
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
264 ).cgs.value
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
265 min_epta = np.minimum(spec, min_epta)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
266 max_epta = np.maximum(spec, max_epta)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
267 # plt.plot(ff,spec)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
268 plt.fill_between(
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
269 fff.value, min_epta, max_epta, color="blue", alpha=0.5, label="EPTA"
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
270 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
271 min_ppta = np.ones(len(fff))
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
272 max_ppta = np.zeros(len(fff))
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
273 for i in range(len(As_ppta)):
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
274 spec = (
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
275 2
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
276 * pi**2
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
277 / 3
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
278 / H0**2
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
279 * fff**2
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
280 * As_ppta[i] ** 2
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
281 * (fff / fref) ** (3 - gammas_ppta[i])
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
282 ).cgs.value
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
283 min_ppta = np.minimum(spec, min_ppta)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
284 max_ppta = np.maximum(spec, max_ppta)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
285 # plt.plot(ff,spec)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
286 plt.fill_between(
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
287 fff.value, min_ppta, max_ppta, color="green", alpha=0.5, label="PPTA"
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
288 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
289
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
290 # PBH constraints
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
291 def PBH(alpha):
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
292 return (5.2 * (1.0 - exp(-1.1 * (abs(alpha - 1.0)) ** (1 / 3))) + 1.0) * (
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
293 alpha > 1
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
294 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
295
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
296 Tgrid = np.logspace(-2.6, 1, 81)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
297 agrid = np.logspace(-1, 2, 61)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
298 bgrid = np.logspace(0, 2, 41)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
299 amin_b_nano_eps1 = 100.0 * np.ones(len(bgrid))
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
300 amax_b_nano_eps1 = 0.0 * np.ones(len(bgrid))
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
301 amin_b_epta_eps1 = 100.0 * np.ones(len(bgrid))
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
302 amax_b_epta_eps1 = 0.0 * np.ones(len(bgrid))
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
303 amin_b_ppta_eps1 = 100.0 * np.ones(len(bgrid))
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
304 amax_b_ppta_eps1 = 0.0 * np.ones(len(bgrid))
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
305 amin_T_nano_eps1 = 100.0 * np.ones(len(Tgrid))
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
306 amax_T_nano_eps1 = 0.0 * np.ones(len(Tgrid))
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
307 amin_T_epta_eps1 = 100.0 * np.ones(len(Tgrid))
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
308 amax_T_epta_eps1 = 0.0 * np.ones(len(Tgrid))
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
309 amin_T_ppta_eps1 = 100.0 * np.ones(len(Tgrid))
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
310 amax_T_ppta_eps1 = 0.0 * np.ones(len(Tgrid))
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
311
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
312 Tmin_b_nano_eps1 = 100.0 * np.ones(len(bgrid))
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
313 Tmax_b_nano_eps1 = 0.0 * np.ones(len(bgrid))
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
314 Tmin_b_epta_eps1 = 100.0 * np.ones(len(bgrid))
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
315 Tmax_b_epta_eps1 = 0.0 * np.ones(len(bgrid))
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
316 Tmin_b_ppta_eps1 = 100.0 * np.ones(len(bgrid))
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
317 Tmax_b_ppta_eps1 = 0.0 * np.ones(len(bgrid))
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
318 B_nano_eps1 = []
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
319 lam_nano_eps1 = []
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
320 B_ppta_eps1 = []
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
321 lam_ppta_eps1 = []
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
322 B_epta_eps1 = []
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
323 lam_epta_eps1 = []
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
324 for i, T in enumerate(Tgrid):
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
325 print(T, len(B_nano_eps1), len(B_epta_eps1), len(B_ppta_eps1))
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
326 pr.report_progress(stage="Progress", progress=1 + 94.0 * (i / len(Tgrid)))
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
327 for j, a in enumerate(agrid):
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
328 for k, b in enumerate(bgrid):
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
329 if b > PBH(a):
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
330 GW_t = GW_turb_Andrii(
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
331 ff * u.Hz, T * u.GeV, a, b, v_w, epsilon_turb
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
332 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
333 lam = GW_t[1].cgs.value
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
334 B = GW_t[2]
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
335 GW_t = GW_t[0]
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
336 GW_s = GW_sound(ff * u.Hz, T * u.GeV, a, b, v_w)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
337 GW = GW_s + GW_t
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
338 GW_interp = np.interp(fff.value, ff, GW)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
339 if sum((GW_interp < max_nano) * (GW_interp > min_nano)) == len(
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
340 fff
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
341 ):
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
342 if a < amin_b_nano_eps1[k]:
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
343 amin_b_nano_eps1[k] = a
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
344 if a > amax_b_nano_eps1[k]:
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
345 amax_b_nano_eps1[k] = a
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
346 if a < amin_T_nano_eps1[i]:
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
347 amin_T_nano_eps1[i] = a
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
348 if a > amax_T_nano_eps1[i]:
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
349 amax_T_nano_eps1[i] = a
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
350 if T < Tmin_b_nano_eps1[k]:
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
351 Tmin_b_nano_eps1[k] = T
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
352 if T > Tmax_b_nano_eps1[k]:
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
353 Tmax_b_nano_eps1[k] = T
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
354 B_nano_eps1.append(B)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
355 lam_nano_eps1.append(lam)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
356 if sum((GW_interp < max_epta) * (GW_interp > min_epta)) == len(
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
357 fff
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
358 ):
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
359 if a < amin_b_epta_eps1[k]:
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
360 amin_b_epta_eps1[k] = a
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
361 if a > amax_b_epta_eps1[k]:
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
362 amax_b_epta_eps1[k] = a
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
363 if a < amin_T_epta_eps1[i]:
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
364 amin_T_epta_eps1[i] = a
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
365 if a > amax_T_epta_eps1[i]:
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
366 amax_T_epta_eps1[i] = a
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
367 if T < Tmin_b_epta_eps1[k]:
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
368 Tmin_b_epta_eps1[k] = T
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
369 if T > Tmax_b_epta_eps1[k]:
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
370 Tmax_b_epta_eps1[k] = T
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
371 B_epta_eps1.append(B)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
372 lam_epta_eps1.append(lam)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
373 if sum((GW_interp < max_ppta) * (GW_interp > min_ppta)) == len(
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
374 fff
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
375 ):
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
376 if a < amin_b_ppta_eps1[k]:
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
377 amin_b_ppta_eps1[k] = a
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
378 if a > amax_b_ppta_eps1[k]:
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
379 amax_b_ppta_eps1[k] = a
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
380 if a < amin_T_ppta_eps1[i]:
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
381 amin_T_ppta_eps1[i] = a
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
382 if a > amax_T_ppta_eps1[i]:
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
383 amax_T_ppta_eps1[i] = a
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
384 if T < Tmin_b_ppta_eps1[k]:
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
385 Tmin_b_ppta_eps1[k] = T
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
386 if T > Tmax_b_ppta_eps1[k]:
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
387 Tmax_b_ppta_eps1[k] = T
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
388 B_ppta_eps1.append(B)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
389 lam_ppta_eps1.append(lam)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
390
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
391 lam_bins = np.logspace(-7, -5, 41)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
392 B_bins = np.logspace(-7, -5, 41)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
393
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
394 h = plt.hist2d(
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
395 np.array(lam_nano_eps1) / 3e24,
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
396 B_nano_eps1,
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
397 bins=[lam_bins, B_bins],
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
398 vmax=1,
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
399 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
400 plt.colorbar()
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
401 cnt = plt.contour(
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
402 sqrt(lam_bins[1:] * lam_bins[:-1]),
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
403 sqrt(B_bins[1:] * B_bins[:-1]),
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
404 np.transpose(h[0]),
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
405 levels=[1],
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
406 colors="white",
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
407 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
408 cont = cnt.get_paths()[0].vertices
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
409 B_nanos_eps1 = cont[:, 1]
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
410 lam_nanos_eps1 = cont[:, 0]
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
411
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
412 plt.xscale("log")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
413 plt.yscale("log")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
414
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
415 lam_bins = np.logspace(-7, -5, 41)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
416 B_bins = np.logspace(-7, -5, 41)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
417
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
418 h = plt.hist2d(
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
419 np.array(lam_epta_eps1) / 3e24,
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
420 B_epta_eps1,
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
421 bins=[lam_bins, B_bins],
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
422 vmax=1,
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
423 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
424 plt.colorbar()
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
425 cnt = plt.contour(
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
426 sqrt(lam_bins[1:] * lam_bins[:-1]),
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
427 sqrt(B_bins[1:] * B_bins[:-1]),
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
428 np.transpose(h[0]),
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
429 levels=[1],
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
430 colors="white",
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
431 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
432 cont = cnt.get_paths()[0].vertices
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
433 B_eptas_eps1 = cont[:, 1]
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
434 lam_eptas_eps1 = cont[:, 0]
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
435
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
436 plt.xscale("log")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
437 plt.yscale("log")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
438
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
439 lam_bins = np.logspace(-7, -5, 41)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
440 B_bins = np.logspace(-7, -5, 41)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
441
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
442 h = plt.hist2d(
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
443 np.array(lam_ppta_eps1) / 3e24,
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
444 B_ppta_eps1,
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
445 bins=[lam_bins, B_bins],
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
446 vmax=1,
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
447 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
448 plt.colorbar()
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
449 cnt = plt.contour(
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
450 sqrt(lam_bins[1:] * lam_bins[:-1]),
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
451 sqrt(B_bins[1:] * B_bins[:-1]),
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
452 np.transpose(h[0]),
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
453 levels=[1],
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
454 colors="white",
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
455 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
456 cont = cnt.get_paths()[0].vertices
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
457 B_pptas_eps1 = cont[:, 1]
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
458 lam_pptas_eps1 = cont[:, 0]
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
459
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
460 plt.xscale("log")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
461 plt.yscale("log")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
462
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
463 import numpy as np
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
464 from matplotlib import pyplot as plt
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
465
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
466 fig, ax = plt.subplots(figsize=(7, 7))
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
467 x = [1e-8, 1e3]
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
468 y1 = [3e-6, 3e-6]
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
469 y2 = [3e-5, 3e-5]
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
470 plt.fill_between(x, y1, y2, color="grey", alpha=0.5)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
471 ax.fill(lam_nanos_eps1, B_nanos_eps1, color="red", alpha=0.3, label="NANOGrav")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
472 ax.fill(lam_eptas_eps1, B_eptas_eps1, color="blue", alpha=0.3, label="EPTA")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
473 ax.fill(lam_pptas_eps1, B_pptas_eps1, color="green", alpha=0.3, label="PPTA")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
474 ax.set_xscale("log")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
475 ax.set_yscale("log")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
476
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
477 lam_med = np.median(lam_nanos_eps1)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
478 B_med = np.median(B_nanos_eps1)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
479 lam_evol = np.logspace(np.log10(lam_med), np.log10(lam_med) + 10, 100)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
480 B_evol = B_med * (lam_evol / lam_med) ** (-5 / 4.0)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
481 mask = B_evol > 10 ** (-8.5) * lam_evol
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
482 lam_evol = lam_evol[mask]
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
483 B_evol = B_evol[mask]
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
484 ax.annotate(
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
485 "",
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
486 xy=(lam_evol[-1], B_evol[-1]),
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
487 xytext=(lam_evol[0], B_evol[0]),
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
488 arrowprops={"arrowstyle": "->", "color": "black"},
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
489 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
490
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
491 lam_evol = np.logspace(np.log10(lam_med), np.log10(lam_med) + 10, 100)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
492 B_evol = B_med * (lam_evol / lam_med) ** (-5 / 2.0)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
493 mask = B_evol > 10 ** (-6.8) * lam_evol
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
494 lam_evol = lam_evol[mask]
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
495 B_evol = B_evol[mask]
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
496 ax.annotate(
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
497 "",
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
498 xy=(lam_evol[-1], B_evol[-1]),
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
499 xytext=(lam_evol[0], B_evol[0]),
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
500 arrowprops={"arrowstyle": "->", "linestyle": "dashed", "color": "black"},
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
501 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
502
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
503 x = np.logspace(-8, 3, 10)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
504 y = 10 ** (-8.5) * x
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
505 ax.plot(x, y, color="grey", linewidth=4, linestyle="dashed")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
506 y = 10 ** (-6.8) * x
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
507 ax.plot(x, y, color="grey", linewidth=4)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
508
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
509 d = np.genfromtxt("MAGIC.csv")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
510 ax.fill_between(d[:, 0], np.zeros(len(d)), d[:, 1], color="grey", alpha=0.5)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
511 ax.set_xlim(1e-8, 1e3)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
512 ax.set_ylim(5e-18, 2e-5)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
513 ax.set_xlabel(r"$\lambda_B$, Mpc")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
514 ax.set_ylabel(r"$B$, G")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
515 plt.text(1, 3e-17, "MAGIC '22", color="black")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
516 plt.text(
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
517 1e-6,
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
518 0.1e-14,
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
519 "Alfvenic larges processed eddies",
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
520 color="black",
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
521 rotation=42,
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
522 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
523
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
524 plt.text(
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
525 0.15e-7,
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
526 0.01e-13,
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
527 "Helicity fluctuaitons / reconnection controlled",
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
528 color="black",
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
529 rotation=41.5,
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
530 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
531 ax.legend(loc="lower left")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
532
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
533 y = [1.5e-12 * 1e1, 1.5e-12, 1.5e-12]
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
534 x = [0.3e-3, 0.03, 1e3]
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
535 plt.plot(x, y, color="olive", linewidth=4, linestyle="dotted")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
536 ax.annotate(
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
537 "",
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
538 xytext=(1, 1.0e-13),
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
539 xy=(1, 1.5e-12),
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
540 arrowprops={"arrowstyle": "->", "color": "olive", "linewidth": 4},
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
541 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
542 plt.text(0.5, 2e-12, "CTA", color="olive")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
543
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
544 x = [0.7e-3, 1e3]
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
545 y = [3e-11, 3e-11]
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
546 plt.plot(x, y, color="olive", linewidth=4, linestyle="dotted")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
547 ax.annotate(
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
548 "",
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
549 xytext=(1, 3e-10),
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
550 xy=(1, 3e-11),
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
551 arrowprops={"arrowstyle": "->", "color": "olive", "linewidth": 4},
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
552 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
553 plt.text(0.5, 1.1e-11, "CMB", color="olive")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
554
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
555 plt.savefig("B_lambdaB.png", bbox_inches="tight")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
556
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
557 fig = plt.figure()
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
558 mask = Tmax_b_nano_eps1 > 0
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
559 plt.fill_between(
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
560 bgrid[mask],
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
561 Tmin_b_nano_eps1[mask],
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
562 Tmax_b_nano_eps1[mask],
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
563 alpha=0.3,
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
564 color="red",
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
565 label="NANOGrav",
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
566 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
567 mask = Tmax_b_epta_eps1 > 0
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
568 plt.fill_between(
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
569 bgrid[mask],
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
570 Tmin_b_epta_eps1[mask],
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
571 Tmax_b_epta_eps1[mask],
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
572 alpha=0.3,
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
573 color="blue",
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
574 label="EPTA",
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
575 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
576 mask = Tmax_b_ppta_eps1 > 0
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
577 plt.fill_between(
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
578 bgrid[mask],
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
579 Tmin_b_ppta_eps1[mask],
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
580 Tmax_b_ppta_eps1[mask],
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
581 alpha=0.3,
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
582 color="green",
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
583 label="PPTA",
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
584 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
585 d = np.genfromtxt("Ellis.csv")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
586 lgT = d[:, 0]
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
587 lgb = d[:, 1]
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
588 plt.plot(10**lgb, 10**lgT, color="black", label="2308.08546")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
589 d = np.genfromtxt("NANO_bubble.csv")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
590 lgT = d[:, 0]
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
591 lgb = -d[:, 1]
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
592 plt.plot(
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
593 10**lgb, 10**lgT, color="black", linestyle="dashed", label="2306.162196"
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
594 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
595 d = np.genfromtxt("NANO_sound.csv")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
596 lgT = d[:, 0]
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
597 lgb = -d[:, 1]
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
598 plt.plot(10**lgb, 10**lgT, color="black", linestyle="dashed")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
599 plt.axhline(0.160, color="black", linestyle="dotted")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
600
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
601 plt.xscale("log")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
602 plt.yscale("log")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
603 plt.ylabel(r"$T$ [GeV]")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
604 plt.xlabel(r"$\beta/H$")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
605 plt.xlim(0.8, 80)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
606 plt.ylim(0.001, 10)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
607 plt.legend(loc="lower left")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
608 plt.savefig("T_Beta.png", bbox_inches="tight")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
609
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
610 fig = plt.figure()
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
611 mask = amax_b_nano_eps1 > 0
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
612 plt.fill_between(
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
613 bgrid[mask],
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
614 amin_b_nano_eps1[mask],
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
615 amax_b_nano_eps1[mask],
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
616 alpha=0.3,
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
617 color="red",
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
618 label="NANOGrav",
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
619 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
620 mask = amax_b_epta_eps1 > 0
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
621 plt.fill_between(
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
622 bgrid[mask],
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
623 amin_b_epta_eps1[mask],
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
624 amax_b_epta_eps1[mask],
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
625 alpha=0.3,
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
626 color="blue",
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
627 label="EPTA",
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
628 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
629 mask = amax_b_ppta_eps1 > 0
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
630 plt.fill_between(
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
631 bgrid[mask],
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
632 amin_b_ppta_eps1[mask],
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
633 amax_b_ppta_eps1[mask],
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
634 alpha=0.3,
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
635 color="green",
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
636 label="PPTA",
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
637 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
638 bpbh = PBH(agrid)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
639 x = [2, 3]
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
640 y = [1, 1]
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
641 y1 = [2, 2]
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
642 plt.fill_between(x, y, y1, color="white", linewidth=0)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
643 plt.fill(bpbh, agrid, color="white")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
644 plt.xscale("log")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
645 plt.yscale("log")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
646 plt.xlim(0.8, 80)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
647 plt.ylim(0.2, 80)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
648 plt.xlabel(r"$\beta/H$")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
649 plt.ylabel(r"$\alpha$")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
650 plt.legend(loc="upper left")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
651 plt.savefig("Alpha_Beta.png", bbox_inches="tight")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
652
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
653 fig = plt.figure()
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
654 mask = amax_T_nano_eps1 > 0
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
655 plt.fill_between(
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
656 Tgrid[mask],
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
657 amin_T_nano_eps1[mask],
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
658 amax_T_nano_eps1[mask],
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
659 alpha=0.3,
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
660 color="red",
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
661 label="NANOGrav",
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
662 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
663 mask = amax_T_epta_eps1 > 0
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
664 plt.fill_between(
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
665 Tgrid[mask],
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
666 amin_T_epta_eps1[mask],
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
667 amax_T_epta_eps1[mask],
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
668 alpha=0.3,
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
669 color="blue",
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
670 label="EPTA",
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
671 linewidth=0,
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
672 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
673 mask = amax_T_ppta_eps1 > 0
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
674 plt.fill_between(
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
675 Tgrid[mask],
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
676 amin_T_ppta_eps1[mask],
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
677 amax_T_ppta_eps1[mask],
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
678 alpha=0.3,
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
679 color="green",
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
680 label="PPTA",
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
681 linewidth=0,
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
682 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
683 d = np.genfromtxt("NANO_alpha_T.csv")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
684 lgT = d[:, 0]
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
685 lga = d[:, 1]
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
686 plt.plot(
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
687 10**lgT, 10**lga, color="black", label="2306.16219", linestyle="dashed"
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
688 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
689 d = np.genfromtxt("NANO_alpha_T1.csv")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
690 lgT = d[:, 0]
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
691 lga = d[:, 1]
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
692 plt.plot(10**lgT, 10**lga, color="black", linestyle="dashed")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
693 plt.axvline(0.16, color="black", linestyle="dotted")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
694 plt.xscale("log")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
695 plt.yscale("log")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
696 plt.xlim(0.001, 10)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
697 plt.ylim(0.2, 80)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
698 plt.xlabel(r"$T$ [GeV]")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
699 plt.ylabel(r"$\alpha$")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
700 plt.legend(loc="upper left")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
701 plt.savefig("Alpha_T.png", bbox_inches="tight")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
702
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
703 bin_image1 = PictureProduct.from_file("B_lambdaB.png")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
704 bin_image2 = PictureProduct.from_file("T_Beta.png")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
705 bin_image3 = PictureProduct.from_file("Alpha_Beta.png")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
706 bin_image4 = PictureProduct.from_file("Alpha_T.png")
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
707
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
708 B_lambdaB_png = bin_image1 # http://odahub.io/ontology#ODAPictureProduct
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
709 T_Beta_png = bin_image2 # http://odahub.io/ontology#ODAPictureProduct
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
710 Alpha_Beta_png = bin_image3 # http://odahub.io/ontology#ODAPictureProduct
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
711 Alpha_T_png = bin_image4 # http://odahub.io/ontology#ODAPictureProduct
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
712
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
713 # output gathering
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
714 _galaxy_meta_data = {}
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
715 _oda_outs = []
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
716 _oda_outs.append(
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
717 (
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
718 "out_Phase_transition_parameters_B_lambdaB_png",
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
719 "B_lambdaB_png_galaxy.output",
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
720 B_lambdaB_png,
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
721 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
722 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
723 _oda_outs.append(
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
724 (
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
725 "out_Phase_transition_parameters_T_Beta_png",
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
726 "T_Beta_png_galaxy.output",
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
727 T_Beta_png,
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
728 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
729 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
730 _oda_outs.append(
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
731 (
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
732 "out_Phase_transition_parameters_Alpha_Beta_png",
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
733 "Alpha_Beta_png_galaxy.output",
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
734 Alpha_Beta_png,
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
735 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
736 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
737 _oda_outs.append(
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
738 (
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
739 "out_Phase_transition_parameters_Alpha_T_png",
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
740 "Alpha_T_png_galaxy.output",
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
741 Alpha_T_png,
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
742 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
743 )
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
744
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
745 for _outn, _outfn, _outv in _oda_outs:
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
746 _galaxy_outfile_name = os.path.join(_galaxy_wd, _outfn)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
747 if isinstance(_outv, str) and os.path.isfile(_outv):
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
748 shutil.move(_outv, _galaxy_outfile_name)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
749 _galaxy_meta_data[_outn] = {"ext": "_sniff_"}
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
750 elif getattr(_outv, "write_fits_file", None):
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
751 _outv.write_fits_file(_galaxy_outfile_name)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
752 _galaxy_meta_data[_outn] = {"ext": "fits"}
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
753 elif getattr(_outv, "write_file", None):
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
754 _outv.write_file(_galaxy_outfile_name)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
755 _galaxy_meta_data[_outn] = {"ext": "_sniff_"}
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
756 else:
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
757 with open(_galaxy_outfile_name, "w") as fd:
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
758 json.dump(_outv, fd, cls=CustomJSONEncoder)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
759 _galaxy_meta_data[_outn] = {"ext": "json"}
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
760
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
761 with open(os.path.join(_galaxy_wd, "galaxy.json"), "w") as fd:
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
762 json.dump(_galaxy_meta_data, fd)
c9fc89ee996e planemo upload for repository https://github.com/esg-epfl-apc/tools-astro/tree/main/tools commit d5be3628a0bdad9747451669fa195f1d2acaf3f2
astroteam
parents:
diff changeset
763 print("*** Job finished successfully ***")