Документ взят из кэша поисковой машины. Адрес оригинального документа : http://hea-www.harvard.edu/ChaMPlane/quantile/pre-release/example.pro Дата изменения: Sat Jun 25 01:11:24 2005 Дата индексирования: Tue Oct 2 08:18:16 2012 Кодировка:

; To run this program
;
; IDL> @example.pro
;
; JaeSub Hong, 2003-2005, version 1.5
; Please report any problem or suggestion at jaesub@head.cfa.harvard.edu
;

range = [0.3,8.0] ; 0.3 - 8.0 keV
frac = [0.25, 0.33, 0.50, 0.67, 0.75]
nofixerror = 0 ; set 0 to fix error and 1 to not fix error

;---------------------------------------------------------
; example of simple source photons
readcol,'src.txt', src, format='f'

Ex = quantile(frac, src, range=range, err_Ex=err_Ex, $
Qx=Qx, err_Qx=err_Qx, $
nofixerror=nofixerror, $
QDx=QDx, QDy=QDy, $
err_QDx=err_QDx, err_QDy=err_QDy)

print,'----------------------'
print,'src ',n_elements(src)
print,'frac ',frac
print,'Ex ',Ex
print,'err_Ex',err_Ex
print,'Qx ',Qx
print,'err_Qx',err_Qx
print,'QDx ',QDx, err_QDx
print,'QDy ',QDy, err_QDy

;---------------------------------------------------------
; example of bkgnd subtraction
readcol,'bkg.txt', bkg, format='f'

ratio = 0.2
Ex = quantile(frac, src, bkg, ratio=ratio, $
range=range, err_Ex=err_Ex, $
Qx=Qx, err_Qx=err_Qx, $
nofixerror=nofixerror, $
QDx=QDx, QDy=QDy, $
err_QDx=err_QDx, err_QDy=err_QDy)

print,'----------------------'
print,'src ',n_elements(src)
print,'bkg ',n_elements(bkg)
print,'net ',n_elements(src)-ratio*n_elements(bkg)
print,'ratio ',ratio
print,'frac ',frac
print,'Ex ',Ex
print,'err_Ex',err_Ex
print,'Qx ',Qx
print,'err_Qx',err_Qx
print,'QDx ',QDx, err_QDx
print,'QDy ',QDy, err_QDy

;---------------------------------------------------------
; example of plotting QCCD

; to compile the plotting routines
qdplot

; rainbow color scheme
loadct,39

; set the plot window
xr =[-1.5,0.5]
yr =[0.5,2.5]
xtick =['0.6','1','2','3','4','5',' 6 keV','7'] ; top x-axis
qd_setwin, xr=xr, yr=yr, range=range, xtick=xtick, stu=0.5

; read grid data
file_grid = 'example_powerlaw_grid.rdb' ; use grid.pl to generate it
grid=qd_readgrid(file_grid, comment=comment)

; label grid
align=[1.0,0.5] ; label alignment
ori =[-20,0] ; orientation
par =[0.01,4] ; location
xo =[0.0,0.0] ; x offset
yo =[0.0,-0.05] ; y offset
col = [80,250] ; color setting
prefix1= ['','','','','','N!dH!n(x10!u22!n)=']
prefix2= '!4C!3='
qd_labelgrid, par, grid, comment=comment, $
prefix1=prefix1, prefix2=prefix2, $
align=align, ori=ori, col=col, xo=xo, yo=yo

; plot grid
qd_plotgrid, grid, comment=comment, col=col

; plot the data point
plots, QDx, QDy, psym=3
; plot errors
plots, QDx+err_QDx, QDy*[1,1]
plots, QDx*[1,1], QDy+err_QDy