gammapbh 1.1.3__py3-none-any.whl

gammapbh/__init__.py

@@ -0,0 +1,42 @@
+# src/gammapbh/__init__.py
+"""
+gammapbh — utilities to load BlackHawk tables, compose PBH gamma spectra,
+and plot/average spectra for various PBH mass distributions.
+
+Units
+-----
+Mass: grams (g)
+Energy: MeV
+Spectral density: dN/dE [MeV^-1 s^-1]
+"""
+
+from importlib.metadata import version as _v
+
+__all__ = [
+    "load_spectra_components",
+    "discover_mass_folders",
+    "mass_function",
+    "mass_function_exact",
+    "mass_function_lognormal",
+    # add other public functions/classes you want users to import
+]
+__version__ = _v("gammapbh")
+
+# Re-export key APIs from submodules (so docs & users can do `from gammapbh import ...`)
+from .cli import (  # or better: move non-CLI helpers to a library module and import from there
+    discover_mass_folders,
+    load_spectra_components,
+    mass_function,
+    mass_function_exact,
+    mass_function_lognormal,
+)
+
+__all__ = ["cli"]
+__version__ = "1.1.3"
+
+
+
+
+
+
+

gammapbh/__main__.py

@@ -0,0 +1,29 @@
+# src/gammapbh/__main__.py
+"""
+Module entry point so `python -m gammapbh [args]` works.
+
+This simply forwards argv to `gammapbh.cli.main`, so things like:
+  - python -m gammapbh --help
+  - python -m gammapbh monochromatic --help
+  - python -m gammapbh  (interactive CLI, if your main supports it)
+behave exactly the same as running the `gammapbh` console script.
+"""
+
+from __future__ import annotations
+import sys
+
+
+def _entry() -> int:
+    # Import inside the function so import errors show as runtime errors
+    # when the module is invoked (mirrors how console_scripts work).
+    from .cli import main  # your argparse/click entry point
+
+    # Forward only the user args (skip the interpreter/module bits)
+    # Expect `main(argv: list[str] | None)` to handle parsing and return
+    # an int exit code or None (treated as 0 here).
+    rv = main(sys.argv[1:])
+    return int(rv) if isinstance(rv, int) else 0
+
+
+if __name__ == "__main__":
+    raise SystemExit(_entry())

gammapbh/blackhawk_data/1.0e+14/1.0e+14.txt

@@ -0,0 +1,59 @@
+destination_folder = 1.0e+14							# name of the output folder in results/
+full_output = 1										# quantity of information displayed (0=less, 1=more)
+interpolation_method = 0							# interpolation in the numerical tables (0=linear, 1=logarithmic)
+
+metric = 0											# BH metric: 0=Kerr, 1=polymerized, 2=charged, 3=higher-dimensional
+
+BH_number = 1										# number of BH masses (should be the number of tabulated masses if spectrum_choice=5)
+Mmin = 1e+14											# lowest BH mass in g (larger than the Planck mass)
+Mmax = 1e+14										# highest BH mass in g (larger than the Planck mass)
+param_number = 1									# number of Kerr spins
+amin = 0.											# lowest Kerr spin
+amax = 0.											# highest Kerr spin
+Qmin = 0.											# lowest Reissner-Norström charge
+Qmax = 0.											# highest Reissner-Norström charge
+
+epsilon_LQG = 1.5									# dimensionless epsilon parameter for the polymerized metric
+a0_LQG = 0.											# minimal area for the polymerized metric in GeV^(-2)
+n = 0.												# number of extra spatial dimensions in higher-dimensional metric
+
+spectrum_choice = 0									# form of the BH distribution: 0=Dirac, 1=log-normal for the mass, 11: log-normal for the number, 2=power-law, 3=critical collapse, 4=peak theory, 5=uniform -1=user-defined
+spectrum_choice_param = 0							# form of the spin dsitribution for each mass: 0=Dirac, 1=uniform, 2=Gaussian
+
+amplitude_lognormal = 1.							# amplitude of the log-normal (mass density) distribution in g.cm^-3
+amplitude_lognormal2 = 1.							# amplitude of the log-normal (number density) distribution in cm^-3
+stand_dev_lognormal = 1.							# dimensionless variance of the log-normal distribution 
+crit_mass_lognormal = 1.							# characteristic mass of the log-normal distribution in g
+
+amplitude_powerlaw = 1.								# amplitude of the power-law distribution in g^(gamma-1).cm^-3
+eqstate_powerlaw = 0.3333							# equation of state of the Universe at the BH formation time P = w.rho
+
+amplitude_critical_collapse = 1.					# amplitude of the critical collapse distribution in g^(-2.85).cm^-3
+crit_mass_critical_collapse = 1.					# characteristic mass of the critical collapse distribution in g
+
+amplitude_uniform = 1.								# amplitude of the uniform mass distribution in cm^(-3)
+
+stand_dev_param_Gaussian = 1.						# standard deviation of the Gaussian spin distribution
+mean_param_Gaussian = 0.5							# mean of the Gaussian spin distribution
+
+table = spin_distribution_BH.txt					# table containing the User's BH distribution
+
+tmin_manual = 0										# 1: user-defined tmin, 0:automatically set tmin
+tmin = 1.e-30										# initial integration time of the evolution of BH in s
+limit = 5000										# iteration limit when computing the time evolution of a single BH
+BH_remnant = 0										# 0: total evaporation, 1: BH relic at mass M_relic
+M_remnant = 1e-4									# BH relic mass in g
+
+E_number = 500										# number of primary particles energies to be simulated
+Emin = 1e-6											# minimal energy in GeV of the primary particles
+Emax = 1e+5											# maximal energy in GeV of the primary particles
+
+grav = 1											# 0=no graviton, 1=emission of gravitons
+add_DM = 1											# 0=no DM added, 1=one DM particle
+m_DM = 1.											# DM mass in GeV
+spin_DM = 0.										# DM spin
+dof_DM = 1.											# number of DM degrees of freedom
+
+primary_only = 0									# 1=no secondary spectrum, 0=secondary spectrum computed
+
+hadronization_choice = 2							# 0=PYTHIA at the BBN epoch, 1=HERWIG at the BBN epoch, 2=PYTHIA (new) at the present epoch, 3=HAZMA at the present epoch, 4=HDMSpectra at the present epoch

gammapbh/blackhawk_data/1.0e+14/BH_spectrum.txt

@@ -0,0 +1,4 @@
+Initial BH comoving number density as a function of their mass and parameter.
+
+      mass/spin    0.00000e+00
+    1.00000e+14    1.00000e+00