Bispectrum Files

import numpy as np

cosmo = 'LC'
space='irsd3'
massbin = 'mbin0'
run = 1

similarly for the bispectrum

bs_path= f'psbs/LH_{cosmo}/bispectrum_hhh_pmwd_{space}_{cosmo}_{massbin_tag}_run{run}_grid144_z0p503.dat'

kF = 2*np.pi/1000. # fundamental mode in units of h/Mpc
k1, k2, k3, Pk1, Pk2, Pk3, B0BSN, BSN, N_tr, B2, B4 = np.loadtxt(bs_path, unpack=True)

# to set k1,k2,k3 in units of h/Mpc, multiply by kF
k1, k2, k3 = k1*kF, k2*kF, k3*kF

B0 = B0BSN - BSN  # we subtract the Poisson shot noise to the bispectrum 'B0(k1,k2,k3) [Mpc/h]^6'

fig = plt.figure(figsize=(12, 4))
plt.scatter(np.arange(k3.shape[0]),B0)

plt.yscale('log')
plt.xlabel('Triangle number', fontsize=18)
plt.ylabel('B0(k1,k2,k3) [Mpc/h]^6', fontsize=18)

Text(0, 0.5, 'B0(k1,k2,k3) [Mpc/h]^6')