mk_hk_rsp

To find the response function for GS, use the mk_hk_rsp program. Calling the program with -help as argument will list the arguments to the program, which at the time this documentation was written were: 
$Id: mk_hk_rsp.c,v 1.19 1999/04/15 18:40:59 tl Exp $
__running on Wed Jul 19 09:23:58 2000
$Id: mk_hk_rsp.c,v 1.19 1999/04/15 18:40:59 tl Exp $

-nm #       ; number of blocks (default: 1000)
-ns #       ; number of events to initially skip
-p #        ; number of events to print (default: 25)
-input s    ; input drive
-tge  lo hi ; ge  time gate limits
-tbgo lo hi ; BGO time gate limits
-egate lo hi; hi res energy gate limits
-line_e #   ; response energy [kev]
-gecal # #  ; off and gain for low res ge
-BGOcal # # ; off and gain for BGO detectors
-nfiles #   ; number of files to sort
-hveto      ; honeycomb suppress gating gamma ray
-addtac2    ; add tac2 (undo on-line sub)

-v          ; version
-help       ; help
A typical script file to run the program will look like this 
mk_hk_rsp \
-input /dev/rmt/1mbn \
-nfiles 1 \
-tge 3816 4200 \
-tbgo 1800 2092 \
-gate 2751 2759 \
-line_e 898.0 \
-gecal -7.809905 4.001766 \
-bgocal 18.32 2.6248 \
-nm 100000 \
-ns 1000000
Which is appropriate for 88Y. The data format should be ususal GS data format WITH EVERYTHING WRITTEN TO TAPE (we need the lo-res signals!). The time and energy gates should in the the 'tape' units. The resulting H and K response functions 
h_rsp.mascii
k_rsp.mascii
matrices will be in ascii format, i.e., the k_rsp.mascii file will contain lines like: 
1 0 2399
1 1 3650
1 2 26
1 3 17
1 4 10
1 5 1
1 6 3
1 7 1
1 8 6
1 9 1
2 0 909
2 1 2932
2 2 2150
2 3 56
2 4 33
2 5 15
.
.
.
110 75 269
110 76 248
110 77 213
110 78 150
110 79 133
110 80 106
110 81 77
110 82 75
110 83 51
110 84 26
110 85 30
110 86 19
110 87 13
110 88 9
110 89 7
110 90 5
110 91 3
110 92 5
110 93 1
110 94 1
110 98 1
110 116 1
Notice that zeros are suppressed. The accompanying description file, k_rsp.mascii.dm, will contain this information 
128
128
sun
ascii
destribing how the data in k_rsp.mascii should be read, in this case: x-dim=128, y-dim=128, date in 'sun' format, and values are in ascii.

(You may change 128 in this file to 256 by hand afterwards if you need to)

You can display the matrices with the 2dplot utility.

[email protected]