Two_Time_Correlation_Function.py 文件源码

def plot_aged_g2( g2_aged, tau=None,timeperframe=1, ylim=None, xlim=None):
    ''''A plot of g2 calculated from two-time'''
    fig = plt.figure(figsize=(8,10)) 
    age_center = list( sorted( g2_aged.keys() ) )
    gs = gridspec.GridSpec(len(age_center),1 ) 
    for n,i in enumerate( age_center):        
        ax = plt.subplot(gs[n]) 
        if tau is None:
            gx= np.arange(len(g2_aged[i])) * timeperframe
        else:
            gx=tau[i]
        marker = markers[n]
        c = colors[n]
        ax.plot( gx,g2_aged[i],  '-%s'%marker, c=c, label=r"$age= %.1f s$"%(i*timeperframe))
        ax.set_xscale('log')
        ax.legend(fontsize='large', loc='best' ) 
        ax.set_xlabel(r"$\tau $ $(s)$", fontsize=18) 
        ax.set_ylabel("g2")    
        if ylim is not None:
            ax.set_ylim( ylim )
        if xlim is not None:
            ax.set_ylim( xlim )   

#####################################
#get fout-time