# Time-stamp: <2001-02-06 18:54:45 dph> 
# MIT Directory: ~dph/h3/CXC/TG/AGfCHRS/Sl/
# File: chips_scripts.sl
# Author: D. Huenemoerder
# Original version: 2000.12.18
#====================================================================

%;;; plot the field

f_evt="~/d1/CXO_Data/II_Peg/Evt/acisf01451_000N003_evt2.fits";
f_get="/home/dph/h3/CXC/TG/Doc/AGfCHRS/Sl/get_evt_data.sl";

e=readfile(f_evt);
r=readfile(f_evt+"[REGION]");

evalfile(f_get);

% get event arrays for different tg_part:

(xm,ym)=get_evt_part(e,2,"sky");  
(xh,yh)=get_evt_part(e,1,"sky");
(xz,yz)=get_evt_part(e,0,"sky");
(xo,yo)=get_evt_part(e,99,"sky");

# get region arrays:

# zero-order region, % HEG box,   % MEG box
(xzo,yzo)=reg_cir_toxy(r,0); 
(xbh,ybh)=reg_box_toxy(r,1);    
(xbm,ybm)=reg_box_toxy(r,2);  

redraw off
pack on

curve x xo y yo
curve x xm y ym
curve x xh y yh
curve x xz y yz

curve x xzo y yzo
curve x xbm y ybm
curve x xbh y ybh

# 1:other, 2:MEG  3:HEG  4:ZO   5:ZOreg   6: MEGreg  7:HEGreg
c 1 symbol size 0.1
c 2 symbol size 0.1
c 3 symbol size 0.1
c 4 symbol size 0.1
c 5 symbol size 0.1
c 6 symbol size 0.1
c 7 symbol size 0.1

c 2 symbol red     
c 3 symbol green   
c 4 symbol blue    

c 5 curve simpleline
c 6 curve simpleline
c 7 curve simpleline

c 5 curve yellow  
c 6 curve magenta 
c 7 curve cyan    


xlabel "sky x [pixel]"
ylabel "sky y [pixel]"
title "Sky field view"

#limits auto auto auto auto
limits 1000 7000 1000 7000
redraw

print postfile Sky_xy.ps
$ pstopnm -portrait -xsize 512 -ysize 512 Sky_xy.ps
$ pnmtopng Sky_xy001.ppm > Sky_xy.png
$ rm Sky_xy.ps Sky_xy001.ppm

# zoom on z.o.

limits 3600 4600 3600 4600
redraw

print postfile Sky_xy_zoom.ps
$ pstopnm -portrait -xsize 512 -ysize 512 Sky_xy_zoom.ps
$ pnmtopng Sky_xy_zoom001.ppm > Sky_xy_zoom.png
$ rm Sky_xy_zoom.ps Sky_xy_zoom001.ppm


##############
# plot tg_lam, tg_d

# MEG

redraw off

c 1 del
c 2 del
c 3 del
c 4 del
d 1 del

clear
redraw

(xm,ym) = get_evt_part_order(e, 2, 1, "tglam");

curve x xm y ym
c 1 symbol size .1
limits auto auto -0.002 0.002

xlabel "Wavelength [Angstrom]"
ylabel "tg_d [degrees]"
title "MEG part photons, +-1st order"

redraw
print postfile MEG_lamd.ps
$ pstopnm -portrait -xsize 512 -ysize 512 MEG_lamd.ps
$ pnmtopng MEG_lamd001.ppm > MEG_lamd.png
$ rm MEG_lamd.ps MEG_lamd001.ppm

# HEG

c 1 del
d 1 del

clear
redraw

(xh,yh) = get_evt_part_order(e, 1, 1, "tglam");
curve x xh y yh
c 1 symbol size .1
limits auto auto -0.002 0.002

xlabel "Wavelength [Angstrom]"
ylabel "tg_d [degrees]"
title "HEG part photons, +-1st orders"

redraw

print postfile HEG_lamd.ps
$ pstopnm -portrait -xsize 512 -ysize 512 HEG_lamd.ps
$ pnmtopng HEG_lamd001.ppm > HEG_lamd.png
$ rm HEG_lamd.ps HEG_lamd001.ppm


##### order-ratio plots

d 1 del
clear
redraw

### MEG
(xm1,ym1) = get_evt_part_order(e, 2, 1, "tgmlamccd");
(xm2,ym2) = get_evt_part_order(e, 2, 2, "tgmlamccd");
(xm3,ym3) = get_evt_part_order(e, 2, 3, "tgmlamccd");
(xmo,ymo) = get_evt_part(e, 2, "tgmlamccd");

yy1=abs(xm1/ym1)
yy2=abs(xm2/ym2)
yy3=abs(xm3/ym3)
yyo=abs(xmo/ymo)

curve x xmo y yyo
curve x xm1 y yy1
curve x xm2 y yy2
curve x xm3 y yy3

c 1 symbol size .1
c 2 symbol size .1
c 3 symbol size .1
c 4 symbol size .1

c 2 symbol red
c 3 symbol green
c 4 symbol blue

limits -30 30 0 5

xlabel "Wavelength*order [Angstrom]"
ylabel "abs(order)~abs(tg_mlam)/ccd_lam [1]"
title "MEG part photons, order-sorting"

redraw
print postfile MEG_osort.ps
$ pstopnm -portrait -xsize 512 -ysize 512 MEG_osort.ps
$ pnmtopng MEG_osort001.ppm > MEG_osort.png
$ rm MEG_osort.ps MEG_osort001.ppm


############ HEG order-sorting
d 1 del
clear
redraw

(xh1,yh1) = get_evt_part_order(e, 1, 1, "tgmlamccd");
(xh2,yh2) = get_evt_part_order(e, 1, 2, "tgmlamccd");
(xh3,yh3) = get_evt_part_order(e, 1, 3, "tgmlamccd");
(xho,yho) = get_evt_part(e, 1, "tgmlamccd");

yy1=abs(xh1/yh1)
yy2=abs(xh2/yh2)
yy3=abs(xh3/yh3)
yyo=abs(xho/yho)

curve x xho y yyo
curve x xh1 y yy1
curve x xh2 y yy2
curve x xh3 y yy3

c 1 symbol size .1
c 2 symbol size .1
c 3 symbol size .1
c 4 symbol size .1

c 2 symbol red
c 3 symbol green
c 4 symbol blue

limits -15 15 0 5

xlabel "Wavelength*order [Angstrom]"
ylabel "abs(order)~abs(tg_mlam)/ccd_lam [1]"
title "HEG part photons, order-sorting"

redraw

print postfile HEG_osort.ps
$ pstopnm -portrait -xsize 512 -ysize 512 HEG_osort.ps
$ pnmtopng HEG_osort001.ppm > HEG_osort.png
$ rm HEG_osort.ps HEG_osort001.ppm