def plotSNfit(result, vc, pidx, digi='D1'): ''' result : fitting results vc, variable carrier pidx power index digi = 'D1' or 'D2' ''' vc.Tz = result.params['Tz'].value if digi == 'D1': f = vc.f1 T = result.params['T1'].value G = result.params['G1'].value Tn = result.params['Tn1'].value data = np.array(vc.cPD1[pidx]) * 1.0 if digi == 'D2': f = vc.f2 T = result.params['T2'].value G = result.params['G2'].value Tn = result.params['Tn2'].value data = np.array(vc.cPD2[pidx]) * 1.0 if digi == 'D12': f = vc.f1 T = result.params['T12'].value G = result.params['G12'].value Tn = result.params['Tn12'].value data = np.array(vc.cPD3[pidx]) * 1.0 if digi == 'D12c': f = vc.f1 T = result.params['T12c'].value G = result.params['G12c'].value Tn = result.params['Tn12c'].value data = np.array(vc.cPD4[pidx]) * 1.0 factor = f*h*vc.B*G # factor to photon # SNf = fitfunc(G, Tn, T, f, vc) Amp = G*vc.B*kB*Tn title2 = (digi + ', RF-Drive ' + str(vc.d2.lin[pidx]) + ' G ' + str(np.round(G/1e7, 2)) + 'e7 = ' + str(np.round(np.log10(G)*10.0, 2)) + 'dB T ' + str(np.round(Tn, 2)) + 'K') dataname = digi + '_' + str(vc.d2.lin[pidx]) + '.dat' gp.c('') gp.figure() # gp.c('clear') gp.s([vc.I*1e6, (data)/factor, (SNf)/factor, np.ones_like(data)*Amp/factor], filename=vc.resultfolder+dataname) gp.c('set title "' + title2 + '"') gp.c('set xrange[-19:19]') gp.c('set key center top') gp.c('plot "' + dataname + '" u 1:2 w l t "Data" ') gp.c('replot "' + dataname + '" u 1:3 w l t "Fit" ') gp.c('replot "' + dataname + '" u 1:4 w l t "Amplifier Noise" ') gp.c('save "' + dataname[:-3] + 'gn"') gp.pdf(dataname[:-3]+'pdf') print dataname[:-3]+'pdf'
def plot3dMSR(uc_3): hyperPoints = cdf.getMSR_ND(uc_s, N=190) gp.s(np.transpose(hyperPoints), filename="plot3dMSR.dat") gp.c('set ticslevel 0') gp.c('unset key') gp.figure('0') c_plot = 'splot "plot3dMSR.dat" with lines ' print c_plot gp.c(c_plot)