Commit ab1d5ba1 authored by Plaszczynski Stephane's avatar Plaszczynski Stephane
Browse files

No commit message

No commit message
parent b6d275d7
.r parnames
if (n_elements(write) eq 0) then write=0
fileout=dir+'.ps'
if (n_elements(t) eq 0) then begin &$
file=dir+'/samples_cut.fits' &$
print,"reading : ",file &$
t=mrdfits(file,1,hdr)&$
......@@ -13,7 +11,6 @@ n=sz[0] &$
;name=['chi2',name] &$
name=strtrim(hdr[5:n_elements(hdr)-4],2) &$
derived,t,name &$
endif
;noms/val
;if (write eq 0) then
if (write ne 0) then openplot,fileout,xs=30,ys=30,charsize=2 else cgdisplay,900,700,tit=dir,/free
......@@ -47,9 +44,9 @@ if (n_elements(name2) ne 0) then begin &$
ibf=strindex(name2,name[i],/silent) &$
cgplot,[bf[ibf],bf[ibf]],!y.crange,/over,line=2 &$
endif &$
if (i le 6 ) then dec=5 else dec=3 &$
if (i le 6 ) then dec=3 else dec=2 &$
ss=cgNumber_Formatter(m,dec=dec)+'+-'+cgNumber_Formatter(s,dec=dec) &$
cgtext,m-2*s,0.6,ss,/data,charsize=1,col='blue' &$
cgtext,m-2*s,0.6,ss,/data,charsize=.8,col='blue' &$
if (ibf ge 0) then begin&$
print,'||',name[i],'||',cgNumber_Formatter(bf[ibf],dec=3),'||',ss,'||' &$
endif else begin &$
......
delvar,t
t1=mrdfits(dir+"/samples.fits",1,hdr)
name=strtrim(hdr[5:n_elements(hdr)-4],2)
t2=mrdfits(dir+"/samples.fits",2)
t3=mrdfits(dir+"/samples.fits",3)
t4=mrdfits(dir+"/samples.fits",4)
in=''
print,"chosse burn_in"
read,in
if (strcmp(in,'s') eq 1) then stop
burn_in=long(in)
print,burn_in
t=[t1[burn_in:*,*],t2[burn_in:*,*],t3[burn_in:*,*],t4[burn_in:*,*]]
help,t
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment