gp("set xrange[" + str(b_low) + ":" + str(b_high) + "]")
# gp("sb=sa")
gp("fit b(x) '" + src + str(v_int) + "/e_files/" + str(v_int) + "_step_" +
str(eq_step) + "_e" + str(echannel) + "' u 1:2:3 via B,xb,sb")
# gp("fit b(x) '"+src+str(v_int)+"/e_files/"+str(v_int)+"_step_"+str(eq_step)+"_e_all' u 1:2:3 via B,xb,sb")
gp("set xrange[" + str(pr_low) + ":" + str(pr_high) + "]")
gp("plot '" + src + str(v_int) + "/e_files/" + str(v_int) + "_step_" +
str(eq_step) + "_e" + str(echannel) +
"' u 1:2:3 with errorbars,a(x),b(x),a(x)+b(x)")
# gp("plot '"+src+str(v_int)+"/e_files/"+str(v_int)+"_step_"+str(eq_step)+"_e_all' u 1:2:3 with errorbars , a(x),b(x),a(x)+b(x)")
# gp("set term postscript color 18")
# gp("set output '"+src+str(v_int)+"/diff_files/"+"plot_"+str(eq_step)+"_"+str(i)+".ps'")
# gp("plot '"+src+str(v_int)+"/e_files/"+str(v_int)+"_step_"+str(eq_step)+"_e"+str(echannel)+"' u 1:2:3 with errorbars , a(x),b(x)")
# gp("set term X11")
a[l - 1, 0] = float(str(gp.getvar("xa")))
b[l - 1, 0] = float(str(gp.getvar("xb")))
c[l - 1, 0] = float(str(gp.getvar("xc")))
d[l - 1, 0] = float(str(gp.getvar("xd")))
a[l - 1, 1] = float(str(gp.getvar("sa")))
b[l - 1, 1] = float(str(gp.getvar("sb")))
c[l - 1, 1] = float(str(gp.getvar("sc")))
d[l - 1, 1] = float(str(gp.getvar("sd")))
a[l - 1, 2] = float(str(gp.getvar("A")))
b[l - 1, 2] = float(str(gp.getvar("B")))
c[l - 1, 2] = float(str(gp.getvar("C")))
d[l - 1, 2] = float(str(gp.getvar("D")))
for j in range(3):
cfile.write("%f\t" % mousexpos[k])
cfile.write("%f\t" % mouseypos[k])
if (mousexpos[k]>0. and mouseypos[k]>0.):
cfile.write("%f\t" % 1.)
else:
cfile.write("%f\t" % 0.)
cfile.write("\n")
cfile.close()
gp("set xrange[%i:%i]" % (xrange[0],xrange[1]))
gp("set yrange[%i:%i]" % (yrange[0],yrange[1]))
gp("set title '%s'" %(labelnames[i]))
gp("splot 'tmpetmat.dat' w pm3d t '', 'tmppointssrimmat.dat' w p lt 5 pt 2 ps 3 t '', 'tmppointsmat.dat' w p lt 3 pt 1 ps 3 t''")
gp("pause mouse '%s' "%(labelnames[i]))
tmpx=float(gp.getvar("MOUSE_X"))
mousexpos[i]=tmpx
tmpy=float(gp.getvar("MOUSE_Y"))
mouseypos[i]=tmpy
print mousexpos[i],mouseypos[i]
gp("set terminal x11 close")
f=open("/home/asterix/berger/projects/ionpositions/test_ion_pos.dat")
w=open("/home/asterix/berger/projects/ionpositions/tmp_ion_pos_end.dat","w")
print "loading test_ion_pos.dat"
for i in range(step-1):
s=f.readline()
w.write(s)
k=s.split()
