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

No commit message

No commit message
parent 5e0fd08e
......@@ -18,10 +18,10 @@ print,+"# free params="+strn(n_elements(names))
print,names
; prof2
ichi=sz[1]-12
;ichi=sz[1]-12
;Mnu
;ichi=sz[1]-11
ichi=sz[1]-11
;
imu=ichi-1
;imu=ichi-2 ; Alens class
......
!x.margin=[6,5]
xfit=0.2
yup=18
yup=7
xr=[-.5,.5]
yr=[0,yup]
xr=[-.1,.3]
yr=[-.2,yup]
dirs=["../hlpTT_bflike_vhl_LCDM_nuscan","../hlpTT_bflike_vhl2_LCDM_nuscan"]
;dirs=["hlpTT_bflike_vhl_LCDMClass_nuscan","hlpTT_bflike_vhl_LCDMClass_nuscan0"]
;dirs=["hlpTT_bflike_vhl_LCDMClass_nuscan","hlpTT_bflike_vhl_JLA_LCDMClass_nuscan"]
dirs=["../hlpTT_bflike_vhl2012_LCDM_nuscan","../hlpTT_bflike_vhl2_LCDM_nuscan"]
;dirs=["../hlpTT_bflike_vhl2012_JLA_BAO_LCDM_nuscan","../hlpTT_bflike_vhl2012_lensing_JLA_BAO_LCDM_nuscan"]
names=dirs
cols=['black','red','blue','green','orange','magenta','violet','purple','navy']
......@@ -36,14 +38,15 @@ sz=size(t,/dim) &$
ichi=1 &$
itau=0 &$
chi2=t[*,ichi] &$
dchi2=chi2 &$
dchi2=chi2-min(chi2,imin) &$
tau=t[*,itau] &$
iord=sort(tau) &$
tau=tau[iord] &$
dchi2=dchi2[iord] &$
diff=dchi2-shift(dchi2,-1) &$
w=where(dchi2 lt 5 and abs(diff) lt 1) &$ ; points used
;w=where(dchi2 gt 0) &$
;w=where(dchi2 lt 5 and abs(diff) lt 1) &$ ; points used
w=where(dchi2 gt 0) &$
x=tau &$
chi2min1=min(chi2) &$
chi2-=chi2min1 &$
......@@ -52,8 +55,8 @@ par=svdfit(x[cut],chi2[cut],3,/double) &$
xmle=-par[1]/(2*par[2]) &$
chi2min=poly(xmle,par) &$
sig=1/sqrt(par[2]) &$
print,dirs[i]+' xmle=',+strn(xmle),' sigma='+strn(sig)+' nsig='+strn(xmle/sig) &$
cgplot,x,chi2-chi2min,psym=16,col=cols[i],/over &$
print," *"+dirs[i]+' xmle=',+strn(xmle),' sigma='+strn(sig)+' #sig='+strn(xmle/sig)+" chi2min="+strn(chi2min1)," raw: "+strn(t[imin,ichi])+" @M="+strn(tau[imin]) &$
cgplot,x,chi2-chi2min,psym=16,col=cols[i],/over,symsize=0.8 &$
xx=range(xr[0],xr[1],N=100) &$
cgplot,xx,poly(xx,par)-chi2min,linestyle=2,col=cols[i],/over &$
xx=range(0,xfit,N=100) &$
......@@ -63,8 +66,8 @@ xf=[-.5,-.4,-.3,-.2,-.1,0.,.1,.2,.3,.4,.5] &$
yf95=[1.49,1.58,1.67,1.77,1.86,1.96,2.06,2.16,2.26,2.36,2.46] &$
mu2=interpol(yf95,xf,xmle/sig) &$
x0=mu2*sig &$
print, 'FC=',x0 &$
y0=poly(x0,par)-chi2min &$
print, 'FC='+strn(y0)+"/ 95\%CL="+strn(x0) &$
cgplot,[0,x0],[y0,y0],col=cols[i],thick=2,/over &$
cgplot,[x0,x0],[0,y0],col=cols[i],thick=2,/over &$
endfor
......
......@@ -48,22 +48,21 @@ for i=1,N-1 do begin &$
if (total(t[*,i]-t[0,i]) eq 0) then continue &$
m=mean(t[*,i],/double) &$
s=stddev(t[*,i]) &$
print,name[i] &$
ff="(f20."+strn(decs[i])+")"&$
print,m,format=ff &$
print,s,format=ff &$
hist_plot,t[*,i],/noplot,x=x,y=y,nbin=101 &$
;hist_plot,t[*,i],/noplot,x=x,y=y,nbin=50 &$
yy=smooth(y,15,/edge_trunc)&$
ff="(e20."+strn(decs[i])+")"&$
;ss=strn(m,format=ff)+"+-"+strn(s,format=ff) &$
ss=cgNumber_Formatter(m,dec=decs[i])+'+-'+cgNumber_Formatter(s,dec=decs[i]) &$
;hist_plot,t[*,i],/noplot,x=x,y=y,nbin=101 &$
hist_plot,t[*,i],/noplot,x=x,y=y,nbin=50 &$
yy=smooth(y,5,/edge_trunc)&$
yy/=max(yy) &$
cgplot,x,yy,xtit=name[i],yr=[0,1],/xs,xticks=3,yticks=3,ytickv=[0.,0.5,1] &$
cgplot,x,yy,xtit=name[i],yr=[0,1],yticks=3,ytickv=[0.,0.5,1] &$
ibf=-1 &$
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 &$
ss=cgNumber_Formatter(m,dec=decs[i])+'+-'+cgNumber_Formatter(s,dec=decs[i]) &$
cgtext,m-2*s,0.6,ss,/data,charsize=.8,col='blue' &$
;cgtext,m-2*s,0.6,ss,/data,charsize=.5,col='blue' &$
al_legend,ss,/top,/right,box=0,textcol=blue,chars=0.6 &$
if (ibf ge 0) then begin&$
print,'||',name[i],'||',cgNumber_Formatter(bf[ibf],dec=3),'||',ss,'||' &$
endif else begin &$
......
......@@ -10,8 +10,8 @@ for j=0,i do cor[i,j]=0 &$
endfor
zone
if keyword_set(write) then openplot,filename
cgDisplay, 800, 800, color = 'white',/free
if keyword_set(write) then openplot,filename,xs=20,ys=20
;cgDisplay, 800, 800, color = 'white',/free
LoadCT, 1
cgloadct, 1, ncolors = 128, bottom = 1, clip = 150
cgloadct, 3, /reverse, ncolors = 128, bottom = 128, clip = 100
......
......@@ -120,8 +120,8 @@ EOBATCH
#####################################################################################################################################
qsub -P P_planck_prod -t 1-5 -pe multicores 8 -q mc_long camelrun
#qsub -P P_planck_prod -t 6-10 -pe openmpi_8 8 -q pa_long camelrun
#qsub -P P_planck_prod -t 1-10 -pe multicores 8 -q mc_long camelrun
qsub -P P_planck_prod -t 1-10 -pe openmpi_8 8 -q pa_long camelrun
#sps=1
#qsub -t 1-10 -pe multicores 8 -l sps=1 -q mc_long camelrun
......
......@@ -79,7 +79,7 @@ algo=ada
length=$NSAMPLES
proposal_cov=$cov
bunchSize=1000
t0=200
t0=2000
ts=10000
scale=0.001
do_move=false
......
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