def test_CrossCorrelation(): ns_astro = NeutrinoSample() ns_astro.inputData(astropath) gs = GALAXY_LIBRARY.get_sample('analy') w_cross = ns_astro.getCrossCorrelation(gs.overdensityalm) if MAKE_TEST_PLOTS: figs = FigureDict() color = [ 'r', 'orange', 'limegreen', 'skyblue', 'mediumslateblue', 'purple', 'grey' ] o_dict = figs.setup_figure('Wcross', xlabel="l", ylabel=r'$w$', figsize=(6, 8)) axes = o_dict['axes'] for i in range(Defaults.NEbin): axes.plot(Defaults.ell, gs.analyCL[0:3 * Defaults.NSIDE] * 10**(i * 2), color='k', lw=2) w_cross[i] *= 10**(i * 2) figs.plot_cl('Wcross', Defaults.ell, np.abs(w_cross), xlabel="l", ylabel=r'$C_{l}$', colors=color, ymin=1e-7, ymax=1e10, lw=3) figs.save_all(testfigpath, 'pdf')
def atmBG_coszenith(energyBin=0): eg = EventGenerator() N_coszenith = eg.atm_gen.coszenith()[energyBin] recovered_values = eg.atmBG_coszenith(int(np.sum(N_coszenith[:, 1])), energyBin) index = np.where(np.abs(recovered_values) > 1) if len(index) > 1: print(index, recovered_values[index]) if MAKE_TEST_PLOTS: figs = FigureDict() figkey = 'N_coszenith' + str(energyBin) o_dict = figs.setup_figure(figkey, xlabel=r'$\cos\,\theta$', ylabel='Number of counts', figsize=(8, 6)) fig = o_dict['fig'] axes = o_dict['axes'] axes.plot(N_coszenith[:, 0], N_coszenith[:, 1], lw=2, label='data') axes.hist(recovered_values, N_coszenith[:, 0], label='mock') fig.legend() figs.save_all(testfigpath, 'pdf')