Provided by You-Hua Chu
A shortened version of the cookbook example

>> valid_times
> observation control = N11P1
> data format (1 for German, 2 for US, 3 for RDF) = 2 > time delay = 15
> maximum allowed MV rate = 170

>> rate
> observation control = N11P1
> data format (1 for German, 2 for US, 3 for RDF) = 2 > gain control (1=>high, 2=>low, 3=>low gain ph in high gain area) = 2

>> ao
> observation control = N11P1
> data format (1 for German, 2 for US, 3 for RDF) = 2 > program control (1 for simple treatment, 2 for full treatment) = 1 > first temp in log10 and scale factor = 6.2 0.8 > second temp in log10 and scale factor = 5.7 0.2

>> rate_fit
> band control (1=>R1, 2=>R2, etc) = 4 > enter value for I0 = 1
> enter 1 to fix value = 0
> enter value for I1 = 1
> enter 1 to fix value = 0
> enter value for I2 = 0.01
> enter 1 to fix value = 0
> enter value for I3 = 0.001
> enter 1 to fix value = 0
> enter value for I4 = 0.0001
> enter 1 to fix value = 0
> enter value for SSX scale ISO = 1
> enter 1 to fix value = 0
> enter value for SSX scale IS1 = 1
> enter 1 to fix value = 0

• SAVE SCREEN OUTPUT OF: TIME OFFSET: T0 = 77679596 I0 0.57323 +/- 0.06347 I1 0.05372 +/- 0.01418 I2 -0.01786 +/- 0.00515 I3 0.00237 +/- 0.00115 I4 0.00051 +/- 0.00004 IS0 0.38260 +/- 0.05251 IS1 -0.00613 +/- 0.00315

>> cp valid_times_all.dat valid_times.dat

>> lte
> observation control = N11P1
> data format (1 for German, 2 for US, 3 for RDF) = 2 > enter T0, I0, I1, I2, I3, I4 = 77679596 0.57323 0.05372 -0.01786 0.00237 0.00051 (numbers from rate_fit)

• SAVE SCREEN OUTPUT OF: NET COUNTS = 30128.0

>> cast_data
> observation control = N11P1
> data format (1 for German, 2 for US, 3 for RDF) = 2 > gain control (1=>high, 2=>low, 3=>low gain ph in high gain area) = 2 > enter nominal pointing RA & Dec in decimal deg = 74.139583 -66.48

>> cast_exp
> observation control = N11P1
> data format (1 for German, 2 for US, 3 for RDF) = 2 > detector control ( 1 => PSPC C, 2=> PSPC B) = 2 > detector control ( 1 => high gain, 2=> low gain) = 2 > enter input band (1=>R1, 2=>R2,3=>R3...) = 4 > enter the nomical pointing RA & Dec in decimal deg = 74.139583 -66.48

• SAVE SCREEN OUTPUT OF: AVERAGE MV = 80.3969

>> cast_part
> observation control = N11P1
> data format (1 for German, 2 for US, 3 for RDF) = 2 > enter detector control (1=>PSPC C, 2=>PSPC B) = 2 > enter gain control (1=>high gain, 2=>low gain) = 2 > enter the number of AP counts = 0
> enter the nomical pointing RA & Dec in decimal deg = 74.139583 -66.48

• SAVE SCREEN OUTPUT OF: TOTAL = 1803.66 (total PB events) >> <cntr>C to halt the program, will come back after # of AP determined

>> run rate_fit to get IS01, IS11, IS02, IS12

• SAVE SCREEN OUTPUT OF: T0 = 77679596 IS01,IS11 = 1.73701 -0.09298 IS02,IS12 = 0.84114 -0.01975 IS04,IS14 = 0.38260 -0.00613

>> tilt
> observation control = N11P1
> data format (1 for German, 2 for US, 3 for RDF) = 2 > enter program control (1 for simple treatment, 2 for full reatment) = 1 > enter T0, IS01, IS11, IS02, IS12, IS04, IS14 = > enter first temp in log10 and scale factor = 6.2 0.8 > enter second temp in log10 and scale factor = 5.7 0.2

• SAVE SCREEN OUTPUT OF: TOTAL COUNTS = 5450.29 R4 TOTAL = 1534.24 R4 THETA = 46.32075 R4 TILT = 3.79171E-02

>> fit_ap
> enter number of PB counts = 1803.66 (from the first run of cast_part) > enter number of SX counts = 5450.29 (total counts from tilt) > enter IL, IH = 11 30 (PH range to be fitted) > enter log T for the first temp = 6.1

• change T until the best reduced chi-square is obtained -----
  • SAVE SCREEN OUTPUT OF: AP = 2590.566 >> <cntr>C to halt the program.

>> cast_part
> observation control = N11P1
> detector control (1=> PSPC C, 2=> PSPC B) = 2 > gain control (1=> high gain, 2=> low gain) = 2 > enter number of AP counts = 2590.566 > enter nominal pointing RA & Dec in decimal deg = 74.139583 -66.48

>> cast_ssx
> enter the band to be cast = 4
> enter total SSX counts, tilt angle, and tilt ampliude = 1534.24 46.32075 0.037917 (R4 results from tilt)
> enter total other counts = 30128.0 (from lte)

>> detect
> enter range in bands for source detection: 4 4 > enter stat significance and count-rate threshold = 5 0

>> final_image
> enter detector control (1 => PSPC C, 2=> PSPC B) = 2 > enter gain control (1 => high gain, 2=> low gain) = 2 > enter band range = 4 4
> hit <CR> only to use the exposure map = mask_4.fits