pro histoplot,a,amean,asig,title=maintitle,chi=deg index=where(finite(a) eq 1) z=moment(a(index)) Na=N_elements(a(index)) print,z b=histogram(a(index),binsize=sqrt(z(1))/5.) print,Na mean=strreal(z(0)) sig=strreal(sqrt(z(1))) Nb=N_elements(b) xx=findgen(Nb)*(max(a(index))-min(a(index)))/(Nb-1)+min(a(index)) if (N_elements(deg) ne 0) then begin if (deg eq 2) then begin yy=exp(-xx) endif else begin yy=xx^(deg/2.-1)*exp(-xx*deg/2.)/gamma(deg/2.)/sqrt(2./deg)^deg endelse endif If (n_elements(maintitle) gt 0) then begin plot,xx,b,title=maintitle+', mean='+mean+', sigma='+sig,xtitle='Data',ytitle='Number of points' if (N_elements(deg) ne 0) then begin plot,xx,b,title=maintitle,xtitle='Data',ytitle='Number of points' print,'chi 2 dof=',deg yy=yy*total(b)/total(yy) oplot,xx,yy,linestyle=2 zcumul=fltarr(Nb) for i=1,Nb-1 do begin zcumul(i)=int_tabulated(findgen(i+1),yy(0:i)) endfor zcumul=zcumul/zcumul(Nb-1) print,'Kolmogorov test=',kolmo(a,zcumul,Nb,0) endif else begin oplot,xx,exp(-(xx-z(0))^2/2./z(1))*Na/sqrt(2*3.141592654)/5.,linestyle=2 endelse endif else begin if (N_elements(deg) ne 0) then begin print,'chi 2 dof=',deg plot,xx,b,title='Histogram'+', mean='+mean+', sigma='+sig+', chi 2 dof='+strtrim(deg,1),xtitle='Data',ytitle='Number of points' yy=yy*total(b)/total(yy) oplot,xx,yy,linestyle=2 Nx=N_elements(xx) zcumul=fltarr(Nb) for i=1,Nx-1 do begin zcumul(i)=int_tabulated(findgen(i+1),yy(0:i)) endfor zcumul=zcumul/zcumul(Nb-1) print,'Kolmogorov test=',kolmo(a,zcumul,Nb,0) endif else begin plot,xx,b,title='Histogram'+', mean='+mean+', sigma='+sig,xtitle='Data',ytitle='Number of points' oplot,xx,exp(-(xx-z(0))^2/2./z(1))*Na/sqrt(2*3.141592654)/5.,linestyle=2 endelse endelse amean=z(0) asig=sqrt(z(1)) end