コード例 #1
0
ファイル: var_plot.py プロジェクト: mwhashim/CDEPNGpy
 def update(val):
     f_NL = sf_NL.val
     gamma = Decimal(sgamma.val).quantize(Decimal('0.1'),rounding=ROUND_DOWN)
     z = Decimal(int(sz.val)).quantize(Decimal('0.0'),rounding=ROUND_DOWN)
     
     APg_Delta_m = var.galaxy_angular_power("wcdm_%s-%s-%s" %(cc.w_x, cc.c2_s, z), 1, f_NL, **cosmo)
     APg_Delta_m11 = var.galaxy_angular_power("gwcdm_%s-%s-%s-%s-%s" %(cc.w_x, cc.c2_s, gamma, cc.Inrct, z), 1, f_NL0, **cosmo)
     APg_Delta_mt = var.galaxy_angular_power("gwcdm_%s-%s-%s-%s-%s" %(cc.w_x, cc.c2_s, gamma, cc.Inrct, z), 1, f_NL, **cosmo)
     
     l.set_ydata(APg_Delta_m0)
     l1.set_ydata(APg_Delta_m)
     l2.set_ydata(APg_Delta_m11)
     l3.set_ydata(APg_Delta_mt)
     
     draw()
コード例 #2
0
ファイル: plot.py プロジェクト: mwhashim/CDEPNGpy
def AP_g(model, BP, f_NL, **cosmo):
    AP_Delta_m = var.galaxy_angular_power(model, BP, f_NL, **cosmo)
    semilogy(cc.l, AP_Delta_m, next(linecycler), linewidth = 2, label = model + " - $f_{NL}: %s$" %f_NL)
    legend(loc = 'best')
    xlabel('$l$')
    ylabel('$l(l+1)/(2\pi) C_l$(a=1)')



#---------------------- END ----------------------
コード例 #3
0
ファイル: var_plot.py プロジェクト: mwhashim/CDEPNGpy
def IN_APg_var_plot():
    #---------- Cosmo parameters -----------------
    cosmo = {'omega_M_0': 0.315, 'omega_b_0': 0.045, 'omega_lambda_0': 1 - 0.315}
    
    ax = subplot(111)
    subplots_adjust(left=0.25, bottom=0.25)
    
    f_NL0 = cc.f_NL[0]
    f_NL1  = cc.f_NL[1]
    
    APg_Delta_m0 = var.galaxy_angular_power("wcdm_%s-%s-%s" %(cc.w_x, cc.c2_s, cc.z), 1, f_NL0, **cosmo)
    APg_Delta_m1 = var.galaxy_angular_power("wcdm_%s-%s-%s" %(cc.w_x, cc.c2_s, cc.z), 1, f_NL1, **cosmo)
    APg_Delta_m01 = var.galaxy_angular_power("gwcdm_%s-%s-%s-%s-%s" %(cc.w_x, cc.c2_s, cc.gamma, cc.Inrct, cc.z), 1, f_NL0, **cosmo)
    
    PPg_Delta_m01 = var.galaxy_angular_power("gwcdm_%s-%s-%s-%s-%s" %(cc.w_x, cc.c2_s, cc.gamma, cc.Inrct, cc.z), 1, f_NL0, **cosmo)[0]
    
    
#    f_NL = linspace(f_NL0, f_NL1, 1000)
#    for i in range(len(f_NL)):
#        PPg_Delta_m1 = var.galaxy_angular_power("wcdm_%s-%s-%s" %(cc.w_x, cc.c2_s, cc.z), 1, f_NL[i], **cosmo)[0]
#        if PPg_Delta_m01 >= PPg_Delta_m1:
#            f_NL01  = f_NL[i]
#        else:
#            exit
#    print "Effective f_NL: ", f_NL01

    l, = semilogy(cc.l, APg_Delta_m0, '-', linewidth = 2)#, label = "wcdm - $f_{NL}$: %s" % cc.f_NL[0] )
    l1, = semilogy(cc.l, APg_Delta_m1, '--', linewidth = 2, label = "wcdm - $f_{NL}$")
    l2, = semilogy(cc.l, APg_Delta_m01, '-.', linewidth = 2, label = "$\Gamma$-wcdm")
    l3, = semilogy(cc.l, APg_Delta_m01, ':', linewidth = 2, label = "$\Gamma$-wcdm + $f_{NL}$" )
    
    #axvline(x = 0.17 * 10**-2, color='r', ls = '--') # vertical line at k_eq
    #axvline(x = 0.14 * 10**-1, color='b', ls = '-.') # vertical line at k_eq
    
    #xlim((cc.k_min, cc.k_max))
    #xlim((10**-3, cc.k_max))
    #axvline(x = 10**-3, color='g', ls = '--') # vertical line at k_eq
    legend(loc = 'best')
    
    xlabel('l')
    ylabel('$l(l+1)/(2\pi) C_l$(a=1)')
    
    axcolor = 'lightgoldenrodyellow'
    
    
    axf_NL = axes([0.25, 0.05, 0.65, 0.03], axisbg=axcolor)
    axgamma = axes([0.25, 0.1, 0.65, 0.03], axisbg=axcolor)
    axz = axes([0.25, 0.15, 0.65, 0.03], axisbg=axcolor)
    
    sf_NL = Slider(axf_NL, '$f_{NL}$', f_NL0, 10.0, valfmt='%1.3f', valinit = f_NL1)
    sgamma = Slider(axgamma, '$\Gamma$', 0.0, 1.0, valfmt='%1.1f', valinit = cc.gamma)
    sz = Slider(axz, '$z$', 0.0, 3.0, valfmt='%1.0f', valinit = cc.z)
    
    def update(val):
        f_NL = sf_NL.val
        gamma = Decimal(sgamma.val).quantize(Decimal('0.1'),rounding=ROUND_DOWN)
        z = Decimal(int(sz.val)).quantize(Decimal('0.0'),rounding=ROUND_DOWN)
        
        APg_Delta_m = var.galaxy_angular_power("wcdm_%s-%s-%s" %(cc.w_x, cc.c2_s, z), 1, f_NL, **cosmo)
        APg_Delta_m11 = var.galaxy_angular_power("gwcdm_%s-%s-%s-%s-%s" %(cc.w_x, cc.c2_s, gamma, cc.Inrct, z), 1, f_NL0, **cosmo)
        APg_Delta_mt = var.galaxy_angular_power("gwcdm_%s-%s-%s-%s-%s" %(cc.w_x, cc.c2_s, gamma, cc.Inrct, z), 1, f_NL, **cosmo)
        
        l.set_ydata(APg_Delta_m0)
        l1.set_ydata(APg_Delta_m)
        l2.set_ydata(APg_Delta_m11)
        l3.set_ydata(APg_Delta_mt)
        
        draw()
    
    sf_NL.on_changed(update)
    sgamma.on_changed(update)
    sz.on_changed(update)
    
    resetax = axes([0.8, 0.0, 0.1, 0.04])
    button = Button(resetax, 'Reset', color=axcolor, hovercolor='0.975')
    
    def reset(event):
        sf_NL.reset()
        sgamma.reset()
        sz.reset()
    
    button.on_clicked(reset)
    show()
    
    return