def Plot(fname1='run/out_chip2_nmos2_drain1_5v_vddat1_5v_gnadtgrounded_rad0Krad.txt',
         fname2='run/beforeIrradiation/test_chip2_pmos5_drain1_5v_vddat1_5v_gnadtfloating.txt'):
    # Process
    info = Parse(fname1)
    #info = Parse('run/out_chip1_pmos8_drain1_5v_vddat1_5v_gnadtgrounded_rad0Kradproper.txt')
    info2 = Parse(fname2)
    #info = Parse('run/beforeIrradiation/test_chip3_pmos5_drain1_5v_vddat1_5v_gnadtfloating.txt')
    #setPlotDefaults(ROOT)
    Style()
    c1 = ROOT.TCanvas("c1","testbeam efficiency",50,50,600,600);
    mg = ROOT.TMultiGraph();
    
    gr1,gr1sq = GetTGraph(info,color=1)
    gr2,gr2sq = GetTGraph(info2,'gr2',color=2)
    gr1sq.Draw();
    #gr2.Draw('same')
    
    leg = ROOT.TLegend(0.45, 0.2, 0.85, 0.4);
    leg.SetBorderSize(0);
    leg.SetFillStyle(0);
    #leg.AddEntry(gr3, "60V Bias");        
    #leg.AddEntry(gr2, "90V Bias");    
    leg.AddEntry(gr1, "PMOS");
    leg.Draw();
    
    c1.Update();
    c1.WaitPrimitive();
    
    data=[]
    for i in info:
        d = DataPoint(i[0],i[1],i[2],i[3])
        data+=[d]
    dset = DataSet(data,'data')
    dset.set_transistor_dimensions(4.0, 0.18)
    dset.calc_characteristics()
    
    print 'Vthr: ',dset._vthr
    print dset.get_characteristics()
    
    vis = Visualizer('title', '2 1')    
    vis.add_plot(dset, ty='rawData', ind='0 0')
    vis.add_plot(dset, ty='sqData', ind='1 0')
    vis.get_plot()
    #raw_input('wait')   
    my_type='NMOS'
    if PMOS:
        my_type = 'PMOS'
    #ta.max_slope(dset)
    dataa=np.loadtxt(fname1,unpack=True, skiprows=2)
    #print dataa
    if NSAMPLE>0:
        ths=[]
        rand = ROOT.TRandom3()
        rand.SetSeed(4)
        for i in range(0,NSAMPLE):
            sampled_data = copy.deepcopy(dataa)
            for n in range(0,len(sampled_data[1])):
                sampled_data[1][n] += rand.Gaus(0.0,sampled_data[1][3]/3.0)
                #print sampled_data[1][n],' ',dataa[1][n]
            ths += [ta.global_plot_single(sampled_data, transistor= my_type, no_show=NSAMPLE>0)]

        mean_ths = np.mean(ths)
        print 'mean: ',mean_ths
        print ths
        rms_ths = []
        for m in ths:
            rms_ths += [m-mean_ths]
        rms_th = math.sqrt(np.sum(np.square(rms_ths))/float(len(rms_ths)-1))
        print 'std. dev.:',rms_th
    else:
        th = ta.global_plot_single(dataa, transistor= my_type, no_show=(NSAMPLE>0))
        print th    
dev_hv = KeithleySMU2400Series(configuration_file,1)
dev_hv.disable_output()
dev_hv.enable_output()

data, number = dev_hv.sweep(dev_hv)

for d in data:
    print 'Vg, Id, Id,errV,errI,   Vd, Ig, Id,errVd,errIg'
    print(d.x, d.y,d.x_err, d.y_err,d.w, d.z,d.w_err, d.z_err)

dset = DataSet(data,'data')
dset.set_transistor_dimensions(4.0, 0.18)
dset.calc_characteristics()

vis = Visualizer('title', '2 1')    
vis.add_plot(dset, ty='rawData', ind='0 0')
vis.add_plot(dset, ty='sqData', ind='1 0')
vis.get_plot()

output_file = 'run/out_hv_iv_chip%s.txt' %(chip_number)
f=time.strftime(output_file)
f=open(f,'w+')
print >> f, time.strftime("[%Y-%m-%d %H:%M:%S]\n"), "<V>,  <Id>,  sd_V,  sd_I, <Vd>,  <Ig>,  sd_Vd,  sd_Ig\n",
for d in data:
    print >> f,d.x, d.y, d.x_err, d.y_err,d.w, d.z,d.w_err, d.z_err
   
f.close()
   
dev_hv.disable_output()
print 'end'