Hi, I'm trying to plot the correlation function of the experimental D_\ell^{TT} from planck, but the function needs a 4 dimensional array and I only have a one dimensional array, so I did this:And the curve I got looks wierd comparing with the expected one, has anyone attempted something similar
Thanks for reading.
Code:
import camb.correlationsell = unbin_cl['ell']xvals = np.linspace(-0.999,0.999,len(ell))theta=np.arccos(xvals)corrs = camb.correlations.cl2corr(np.tile(np.array(D_ell).reshape(-1, 1), 4),xvals)fig, ax = plt.subplots(2,2, figsize = (15,15))ax[0,0].plot(theta,corrs[:,0], color='k')ax[0,0].set_title(r'$TT\, [\mu K^2]$')ax[0,1].plot(theta, corrs[:,2]);ax[0,1].set_title(r'BB')ax[1,0].plot(theta,corrs[:,1], color='k')ax[1,0].set_title(r'$EE\, [\mu K^2]$')ax[1,1].plot(theta,corrs[:,3], color='k')ax[1,1].set_title(r'$TE\, [\mu K^2]$');for ax in ax.reshape(-1): ax.set_xlabel(r'$\theta [rad]$');Thanks for reading.
Statistics: Posted by Alv Men — February 28 2025