Created
December 5, 2025 20:23
-
-
Save BrianMartell/f125596f68c64f9c8654ffdc8b34f1f6 to your computer and use it in GitHub Desktop.
PUH-BrianMartell PUH_v6_Acoustic_Peaks_Sim.py- Updated New Py Code
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| import numpy as np | |
| import matplotlib.pyplot as plt | |
| ell = np.arange(2, 2500) | |
| # Standard ΛCDM-like power spectrum (mock) | |
| Cl = 1e-10 * ell**(-2) * (1 + np.sin(ell/220 * np.pi)) # very rough peaks | |
| plt.figure(figsize=(10,6)) | |
| plt.loglog(ell, Cl, 'b-', label='PUH v6 (standard physics after rebound)') | |
| plt.xlabel('Multipole $\ell$') | |
| plt.ylabel('Power $C_\ell$ (arb.)') | |
| plt.title('PUH v6: CMB Acoustic Peaks from Standard Post-Rebound Physics') | |
| plt.legend() | |
| plt.grid(True) | |
| plt.savefig('puh_acoustic_peaks.png', dpi=300) | |
| plt.show() |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment