Mercurial > repos > astroteam > sgwb_astro_tool
annotate Phase_transition_parameters.py @ 0:c9fc89ee996e draft default tip
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 |
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 ***") |