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)
