| About ROSAT |
ROSAT Home Page | ROSAT Images |
|---|
Sample run through US data set RP900148 on a unix system. The symbol >> indicates commands to the operating system. The > symbol indicates user input to the programs. The ! symbol indicates annotation for this listing.
>> ls ! list out data sets in work
! directory
uma_ATTITUDE.FITS uma_EVENTRATES.FITS uma_ORBIT.FITS uma_ATTITUDE1.FITS uma_EVENTS.FITS
>> valid_times.o ! invoke the program VALID_TIMES
***** VALID_TIMES *****
Enter observation control
> uma ! observation control is "uma"
OBS =uma
Enter 1 for German data
Enter 2 for US data
Enter 3 for RDF data
> 2 ! US data set
OPEN
uma_EVENTS.FITS
1 44435430.000000 44435432.000000 ! These are the accepted times
2 44435460.000000 44436663.000000 ! (GTI) in the observation data
3 44436937.000000 44436967.000000 ! set
4 44441078.000000 44442716.000000
5 44491486.000000 44492387.000000
6 44492614.000000 44492618.000000
7 44508959.000000 44511695.000000
8 44521431.000000 44523193.000000
9 44526918.000000 44528939.000000
10 44577681.000000 44578653.000000
11 44607652.000000 44609310.000000
12 44609312.000000 44609412.000000
13 44663932.000000 44664901.000000
14 44669788.000000 44672641.000000
OPEN
uma_EVENTRATES.FITS
Enter a time delay in seconds for the counter to be ready
> 15 ! input time delay of 15 seconds
Time Delay = 15.000000000000
Enter the maximum allowed Master Veto rate
> 170 ! input master veto threshold of
! 170 counts/second
Maximum Master Veto = 170
1 44435463 44436663 ! list of selected time
2 44441111 44442715 ! intervals
3 44491518 44492122
4 44492358 44492387
5 44508989 44511695
6 44521466 44523192
7 44527018 44528939
8 44577712 44578417
9 44578485 44578653
10 44607685 44609310
11 44609312 44609411
12 44663965 44664602
13 44664670 44664901
14 44669821 44672641
>> rate.o ! invoke the program RATE
Enter observation control
> uma ! observation control is "uma"
OBS =uma
Enter 1 for German data
Enter 2 for US data
Enter 3 for RDF data
> 2 ! US data set
Enter the gain control
1 => high gain
2 => low gain
3 => low gain pulse height, high gain area
> 2 ! data collected in low gain mode
uma_EVENTS.FITS
BEGIN = 44435463 END = 44436663
BEGIN = 44441111 END = 44442715
BEGIN = 44491518 END = 44492122
BEGIN = 44492358 END = 44492387
BEGIN = 44508989 END = 44511695
BEGIN = 44521466 END = 44523192
BEGIN = 44527018 END = 44528939
BEGIN = 44577712 END = 44578417
BEGIN = 44578485 END = 44578653
BEGIN = 44607685 END = 44609310
BEGIN = 44609312 END = 44609411
BEGIN = 44663965 END = 44664602
BEGIN = 44664670 END = 44664901
BEGIN = 44669821 END = 44672641
TOTAL NUMBER OF INTERVALS IS 14
10000
20000
30000
40000
50000
60000
70000
80000
90000
100000
110000
120000
130000
140000
150000
160000
170000
180000
190000
200000
210000
220000
230000
240000
250000
260000
229453
uma_EVENTRATES.FITS
>> ao.o ! invoke the program AO
Enter observation control
> uma ! observation control is "uma"
OBS =uma
Enter 1 for German data
Enter 2 for US data
Enter 3 for RDF data
> 2 ! US data set
Enter program control
1 => simple treatment
2 => full treatment
> 1 ! use the simple method
uma_ORBIT.FITS
uma_ATTITUDE1.FITS
uma_EVENTRATES.FITS
BEGIN = 44435463 END = 44436663
BEGIN = 44441111 END = 44442715
BEGIN = 44491518 END = 44492122
BEGIN = 44492358 END = 44492387
BEGIN = 44508989 END = 44511695
BEGIN = 44521466 END = 44523192
BEGIN = 44527018 END = 44528939
BEGIN = 44577712 END = 44578417
BEGIN = 44578485 END = 44578653
BEGIN = 44607685 END = 44609310
BEGIN = 44609312 END = 44609411
BEGIN = 44663965 END = 44664602
BEGIN = 44664670 END = 44664901
BEGIN = 44669821 END = 44672641
BEGIN = 0 END = 0
TOTAL NUMBER OF INTERVALS IS 14
Enter first temperature in log10 and scale factor
> 6.2 0.8 ! first component has log10 T=6.2
! and scale factor of 0.8
25 26 6.15000 6.20000
Enter second temperature in log10 and scale factor
> 5.7 0.2 ! second component has log10 T=5.7
! and scale factor of 0.8
15 16 5.65000 5.70000
558.17847 0.130 0.274 0.314 0.534 0.137 0.0070 0.0000 1.396
558.17914 0.130 0.275 0.314 0.535 0.138 0.0070 0.0000 1.399
558.17987 0.136 0.288 0.329 0.561 0.144 0.0073 0.0000 1.466
558.18054 0.143 0.303 0.346 0.590 0.151 0.0077 0.0000 1.541
.
.
! The program lists out the nominal count rates for the seven bands plus
! the total count rate.
! Time (d) R1 R2 R3 R4 R5 R6 R7 TOTAL
558.19165 0.761 1.614 1.844 3.137 0.806 0.0409 0.0002 8.203
558.19238 0.936 1.985 2.268 3.858 0.991 0.0503 0.0002 10.089
558.19305 1.140 2.418 2.760 4.695 1.206 0.0612 0.0003 12.280
558.24377 0.139 0.294 0.336 0.572 0.147 0.0075 0.0000 1.496
558.24445 0.138 0.291 0.333 0.566 0.145 0.0074 0.0000 1.480
558.24512 0.141 0.297 0.340 0.578 0.149 0.0075 0.0000 1.512
.
.
560.92224 2.329 4.954 5.662 9.630 2.473 0.1254 0.0005 25.174
560.92291 3.620 7.694 8.767 14.902 3.826 0.1941 0.0008 39.003
560.92358 6.207 13.177 14.926 25.343 6.505 0.3303 0.0014 66.489
>> rate_fit.o ! invoke the program RATE_FIT
Enter band control
1 => R1
2 => R2
3 => R3
4 => R4
5 => R5
6 => R6
7 => R7
> 4 ! fit the band R4 light curve
BEGIN = 44435463 END = 44436663
BEGIN = 44441111 END = 44442715
BEGIN = 44491518 END = 44492122
BEGIN = 44492358 END = 44492387
BEGIN = 44508989 END = 44511695
BEGIN = 44521466 END = 44523192
BEGIN = 44527018 END = 44528939
BEGIN = 44577712 END = 44578417
BEGIN = 44578485 END = 44578653
BEGIN = 44607685 END = 44609310
BEGIN = 44609312 END = 44609411
BEGIN = 44663965 END = 44664602
BEGIN = 44664670 END = 44664901
BEGIN = 44669821 END = 44672641
TOTAL NUMBER OF INTERVALS IS 14
YMIN = 0.267000 YMAX= 14.4270
529 ! number of data points
Enter the value for I0
> 1 ! guess for I0 is 1.0
Enter 1 to fix the value
> 0 ! let the parameter vary
Enter the value for I1
> 1 ! guess for I1 is 1.0, fixed
Enter 1 to fix the value
> 0 ! let the parameter vary
Enter the value for I2
> 1 ! guess for I2 is 1.0, fixed
Enter 1 to fix the value
> 0 ! let the parameter vary
Enter the value for I3
> 1 ! guess for I3 is 1.0, fixed
Enter 1 to fix the value
> 0 ! let the parameter vary
Enter the value for I4
> 1 ! guess for I4 is 1.0, fixed
Enter 1 to fix the value
> 0 ! let the parameter vary
Enter the value for SSX scale IS0
> 1 ! guess for IS0 is 1.0
Enter 1 to fix the value
> 0 ! let the parameter vary
Enter the value for SSX scale IS1
> 1 ! guess for IS1 is 1.0, fixed
Enter 1 to fix the value
> 0 ! let the parameter vary
Initial Model Value Settings:
*I0= 1.00000000 *I1= 1.00000000 *I2= 1.00000000
*I3= 1.00000000
*I4= 1.00000000 *IS0= 1.00000000 *IS0= 1.00000000
ITER = 2 ALAMDA = 1.00000E-04 M.L.E. = 1867.70 Intermediate step..
*I0= 1.06111407 *I1= -0.80486083 *I2= 0.18138754
*I3= -0.13042748
*I4= 0.10288483 *IS0= 0.29994959 *IS0= -0.05306673
.
.
.
ITER = 18 ALAMDA = 100.0000 M.L.E. = 1793.08 Intermediate step..
*I0= 1.19334090 *I1= -1.58024240 *I2= 1.13640761
*I3= -0.47476155
*I4= 0.12959373 *IS0= 0.29634619 *IS0= -0.05163757
Final Model Values: (*=Floating)
*I0= 1.19334090 *I1= -1.58024240 *I2= 1.13640761
*I3= -0.47476155
*I4= 0.12959373 *IS0= 0.29634619 *IS0= -0.05163757
Np = 7 N = 529 M.L.E. = 1793.
Final parameter values and the formal standard errors
I0 1.19334 +/- 0.03145 I1 -1.58024 +/- 0.27440 I2 1.13641 +/- 0.56276 I3 -0.47476 +/- 0.37396 I4 0.12959 +/- 0.07789 IS0 0.29635 +/- 0.00738 IS1 -0.05164 +/- 0.00585
TIME OFFSET: T0 = 44435376
LONG TERM: I0,I1,I2,I3,I4= 1.193 -1.580 1.136 -0.475 0.130
SSX: IS0,IS1 = 0.2963 -0.0516
NO = 529 NP = 7 CHISQ =1793.1 REDCHISQ = 3.4350
! This is a messy case where all parameters are used and still the fit ! is terrible. Examination of the plot produced (plot_4_full.qdp) shows ! that a time interval towards the end of the observation data set is ! strongly affected by short-term enhancements. By removing that interval ! and the times of greatest excesses from scattered solar X-rays, a ! considerably better fit is achieved.
>> rate_fit.o ! invoke the program RATE_FIT
Enter band control
1 => R1
2 => R2
3 => R3
4 => R4
5 => R5
6 => R6
7 => R7
> 4 ! fit the band R4 light curve
BEGIN = 44435463 END = 44436663
BEGIN = 44441111 END = 44442328
BEGIN = 44491518 END = 44492122
BEGIN = 44492358 END = 44492387
BEGIN = 44508989 END = 44511584
BEGIN = 44521466 END = 44522980
BEGIN = 44527018 END = 44528712
BEGIN = 44577712 END = 44578417
BEGIN = 44578485 END = 44578653
BEGIN = 44607685 END = 44609232
BEGIN = 44609312 END = 44609114
BEGIN = 44671308 END = 44672388
TOTAL NUMBER OF INTERVALS IS 12
YMIN = 0.267000 YMAX= 3.88600
415 ! number of data points
Enter the value for I0
> 1 ! guess for I0 is 1.0
Enter 1 to fix the value
> 0 ! let the parameter vary
Enter the value for I1
> 1 ! guess for I1 is 1.0, fixed
Enter 1 to fix the value
> 0 ! let the parameter vary
Enter the value for I2
> 1 ! guess for I2 is 1.0, fixed
Enter 1 to fix the value
> 0 ! let the parameter vary
Enter the value for I3
> 1 ! guess for I3 is 1.0, fixed
Enter 1 to fix the value
> 0 ! let the parameter vary
Enter the value for I4
> 1 ! guess for I4 is 1.0, fixed
Enter 1 to fix the value
> 0 ! let the parameter vary
Enter the value for SSX scale IS0
> 1 ! guess for IS0 is 1.0
Enter 1 to fix the value
> 0 ! let the parameter vary
Enter the value for SSX scale IS1
> 0 ! guess for IS1 is 0.0, fixed
Initial Model Value Settings:
*I0= 1.00000000 *I1= 1.00000000 *I2= 1.00000000
*I3= 1.00000000
*I4= 1.00000000 *IS0= 1.00000000 IS0= 0.00000000
ITER = 2 ALAMDA = 1.00000E-04 M.L.E. = 717.224 Intermediate step..
*I0= 1.17301834 *I1= -1.08642650 *I2= 0.37974840
*I3= -0.06773710
*I4= 0.04088503 *IS0= 0.21684825 IS0= 0.00000000
.
.
.
ITER = 21 ALAMDA = 1000.000 M.L.E. = 668.732 Intermediate step..
*I0= 1.30878007 *I1= -2.48289609 *I2= 2.84655619
*I3= -1.51039016
*I4= 0.30858064 *IS0= 0.21624272 IS0= 0.00000000
Final Model Values: (*=Floating)
*I0= 1.30878007 *I1= -2.48289609 *I2= 2.84655619
*I3= -1.51039016
*I4= 0.30858064 *IS0= 0.21624272 IS0= 0.00000000
Np = 6 N = 415 M.L.E. = 668.7
Final parameter values and the formal standard errors
I0 1.30878 +/- 0.03057 I1 -2.48290 +/- 0.27820 I2 2.84656 +/- 0.57429 I3 -1.51039 +/- 0.38234 I4 0.30858 +/- 0.07970 IS0 0.21624 +/- 0.00646
TIME OFFSET: T0 = 44435376
LONG TERM: I0,I1,I2,I3,I4= 1.309 -2.483 2.847 -1.510 0.309
SSX: IS0,IS1 = 0.2162 0.0000
NO = 415 NP = 6 CHISQ = 668.7 REDCHISQ = 1.6350
! This was a relatively simple case where a reasonable (as opposed to ! statistically significant fit) fit was achieved using most of the ! parameters and exclusion of not a great amount of data (~20%) .
>> cp valid_times_all.dat valid_times.dat ! copy the valid_times_all.dat
! accepted times file to
! valid_times.dat:given the
! quality of the fit, there is
! no apparent reason to exclude
! any further data
>> lte.o ! invoke the program LTE
Enter observation control
> uma ! observation control is "uma"
OBS =uma
Enter 1 for German data
Enter 2 for US data
Enter 3 for RDF data
> 2 ! US data set
ENTER T0, I0, I1, I2, I3, I4
> 44435376 1.30878007 -2.48289609 2.84655619 -1.51039016 0.30858064
! fitted values from RATE_FIT
! when the higher-order elements
! of the quadratic are used,
! more digit should be used to
! avoid round-off error
uma_EVENTRATES.FITS
BEGIN = 44435463 END = 44436663
BEGIN = 44441111 END = 44442328
BEGIN = 44491518 END = 44492122
BEGIN = 44492358 END = 44492387
BEGIN = 44508989 END = 44511584
BEGIN = 44521466 END = 44522980
BEGIN = 44527018 END = 44528712
BEGIN = 44577712 END = 44578417
BEGIN = 44578485 END = 44578653
BEGIN = 44607685 END = 44609232
BEGIN = 44609312 END = 44609114
BEGIN = 44671308 END = 44672388
TOTAL NUMBER OF INTERVALS IS 12
TOTAL TIME = 12061.0 ! observation live time
TOTAL COUNTS = 8246.98 ! number of counts in polynomial,
! in this case only I0 component
NET COUNTS = 2785.43 ! integral (number of counts) of
! polynomial above minimum value
! of polynomial
>> cast_data.o ! invoke program CAST_DATA
Enter observation control
> uma ! observation control is "uma"
OBS =uma
Enter 1 for German data
Enter 2 for US data
Enter 3 for RDF data
> 2 ! US data set
Enter the gain control
1 => high gain
2 => low gain
3 => low gain pulse height, high gain area
> 2 ! data collected in low gain mode
Enter the nominal pointing direction (RA, Dec)
in decimal degrees
> 172.96 63.84 ! this information will be used
! in the output fits files
DSCS= 44435430.000000 DSCSO= 44435432.000000
DSCS= 44435460.000000 DSCSO= 44436663.000000
DSCS= 44436937.000000 DSCSO= 44436967.000000
DSCS= 44441078.000000 DSCSO= 44442716.000000
DSCS= 44491486.000000 DSCSO= 44492387.000000
DSCS= 44492614.000000 DSCSO= 44492618.000000
DSCS= 44508959.000000 DSCSO= 44511695.000000
DSCS= 44521431.000000 DSCSO= 44523193.000000
DSCS= 44526918.000000 DSCSO= 44528939.000000
DSCS= 44577681.000000 DSCSO= 44578653.000000
DSCS= 44607652.000000 DSCSO= 44609310.000000
DSCS= 44609312.000000 DSCSO= 44609412.000000
DSCS= 44663932.000000 DSCSO= 44664901.000000
DSCS= 44669788.000000 DSCSO= 44672641.000000
BEGIN = 44435463 END = 44436663
BEGIN = 44441111 END = 44442328
BEGIN = 44491518 END = 44492122
BEGIN = 44492358 END = 44492387
BEGIN = 44508989 END = 44511584
BEGIN = 44521466 END = 44522980
BEGIN = 44527018 END = 44528712
BEGIN = 44577712 END = 44578417
BEGIN = 44578485 END = 44578653
BEGIN = 44607685 END = 44609232
BEGIN = 44609312 END = 44609114
BEGIN = 44671308 END = 44672388
TOTAL NUMBER OF INTERVALS IS 12
IE= 10000 ! event counter
IE= 20000
IE= 30000
IE= 40000
IE= 50000
IE= 60000
IE= 70000
IE= 80000
IE= 90000
IE= 100000
IE= 110000
IE= 120000
IE= 130000
IE= 140000
IE= 150000
IE= 160000
IE= 170000
IE= 180000
IE= 190000
IE= 200000
IE= 210000
IE= 220000
IE= 230000
IE= 240000
IE= 250000
IE= 260000
IEE= 145377 ! number of events in accepted
! time
Number of counts in each of the seven bands
8-201 (or 11-201): 136928
8-19 (or 11-19): 39066 ! band R1L
20-41: 51788 ! band R2
42-51: 6980 ! band R3
52-69: 12387 ! band R4
70-90: 10423 ! band R5
91-131: 10859 ! band R6
132-201: 5425 ! band R7
Number of rejected counts
PH < 8 (or 11): 6437 ! out of accepted pulse-height range
PH > 201: 1502 ! out of accepted pulse-height range
After Pulse: 490 ! events within 0.35 ms of
! preceeding event
Bad Region: 11556 ! out of accepted detector region
Not Valid: 106379 ! out of accepted time intervals
Outside of image: 0 ! outside of 512x512 array
File to store diagnostic information:
uma_diag.out ! file for diagnostic output
File to store PHDs for afterpulse analysis:
spec.dat ! file for PHD output
Writing output file
b_1l.fits
Writing output file
b_2.fits
Writing output file
b_3.fits
Writing output file
b_4.fits
Writing output file
b_5.fits
Writing output file
b_6.fits
Writing output file
b_7.fits
>> cast_exp.o ! invoke program CAST_EXP
Enter observation control
> uma ! observation control is "uma"
OBS =uma
Enter 1 for German data
Enter 2 for US data
Enter 3 for RDF data
> 2 ! US data set
Enter detector control
1 => PSPC C before 26 January 1991
2 => PSPC B after 26 January 1991
> 2 ! PSPC B
IDET = 2
Enter detector control
1 => high gain before 11 October 1991
2 => low gain after 11 October 1991
> 2 ! low gain
IG = 2
Enter input band
1 => CH 8-19 R1 2 => CH 20-41 R2 3 => CH 42-51 R3 4 => CH 52-69 R4 5 => CH 70-90 R5
6 => CH 91-131 R6
7 => CH 132-201 R7
8 => CH 11-19 R1L
9 => CH 8-41 SASS soft band
10 => CH 52-90
11 => CH 91-201
12 => CH 42-131 SASS exposure band
13 => CH 42-201 SASS hard band
> 4 ! select band R4
IB = 4
Enter the nominal pointing direction (RA, Dec)
in decimal degrees
172.96 63.84
Selection: 4
/caldb/data/rosat/pspc/cpf/detmaps/det_52_69_b.fits
! detector map
Reading input file
/caldb/data/rosat/pspc/cpf/detmaps/det_52_69_b.fits
uma_ATTITUDE.FITS
uma_EVENTRATES.FITS
BEGIN = 44441111 END = 44442328
BEGIN = 44491518 END = 44492122
BEGIN = 44492358 END = 44492387
BEGIN = 44508989 END = 44511584
BEGIN = 44521466 END = 44522980
BEGIN = 44527018 END = 44528712
BEGIN = 44577712 END = 44578417
BEGIN = 44578485 END = 44578653
BEGIN = 44607685 END = 44609232
BEGIN = 44609312 END = 44609114
BEGIN = 44671308 END = 44672388
TOTAL NUMBER OF INTERVALS IS 12
1000
2000
3000
4000
5000
6000
7000
8000
9000
10000
11000
12000
LIST LENGTH = 442 ! number of independent entries
! in the offset list
GOOD STEPS = 12361
BAD STEPS = 2
Total Time = 12362.0 Total Exposure = 12059.0
! total number of seconds and
! total livetime
AVERAGE MV = 80.9570 ! average master veto rate
! during accepted time intervals
File for storing diagnostic information:
uma_diag.out ! file to write diagnostic output
PROCESSING 10% COMPLETE
PROCESSING 20% COMPLETE
PROCESSING 30% COMPLETE
PROCESSING 40% COMPLETE
PROCESSING 50% COMPLETE
PROCESSING 60% COMPLETE
PROCESSING 70% COMPLETE
PROCESSING 80% COMPLETE
PROCESSING 90% COMPLETE
PROCESSING 100% COMPLETE
Writing output file
expose_4.fits
>> tilt.o ! invoke program TILT
Enter observation control
> uma ! observation control is "uma"
OBS =uma
Enter 1 for German data
Enter 2 for US data
Enter 3 for RDF data
> 2 ! US data set
Enter program control
1 => simple treatment
2 => full treatment
> 1 ! simple treatment
uma_ORBIT.FITS
uma_ATTITUDE1.FITS
uma_EVENTRATES.FITS
Enter T0, IS01, IS11, IS02, IS12, IS04, IS14
> 44435376 0.1 0 0.1 0 0.21624 0 ! parameters from RATE_FIT
! use nonzero values for IS01
! and IS02 to avoid divide by 0
BEGIN = 44435463 END = 44436663
BEGIN = 44441111 END = 44442328
BEGIN = 44491518 END = 44492122
BEGIN = 44492358 END = 44492387
BEGIN = 44508989 END = 44511584
BEGIN = 44521466 END = 44522980
BEGIN = 44527018 END = 44528712
BEGIN = 44577712 END = 44578417
BEGIN = 44578485 END = 44578653
BEGIN = 44607685 END = 44609232
BEGIN = 44609312 END = 44609114
BEGIN = 44671308 END = 44672388
TOTAL NUMBER OF INTERVALS IS 12
Enter first temperature in log10 and scale factor
> 6.2 0.8 ! same spectrum as for AO
25 26 6.15000 6.20000
Enter second temperature in log10 and scale factor
> 5.7 0.2 ! same spectrum as for AO
15 16 5.65000 5.70000
558.179138 0.013 0.028 0.051 0.119 0.031 0.002 0.000 0.245
558.179138 0.014 0.029 0.052 0.121 0.031 0.002 0.000 0.247
558.179871 0.014 0.029 0.053 0.124 0.032 0.002 0.000 0.254
558.179871 0.014 0.030 0.054 0.126 0.032 0.002 0.000 0.258
58.61 0.50 0.00 1.05 0.01 4.43 0.05
58.61 1.14 0.01 0.06 0.00 0.00 0.00
58.61 0.82 0.01 1.73 0.02 7.28 0.09
.
.
.
560.920837 0.121 0.258 0.468 1.089 0.280 0.014 0.000 2.230
560.920837 0.134 0.284 0.516 1.200 0.308 0.016 0.000 2.457
65.20 8.02 0.81 17.05 1.74 71.87 7.39
65.20 18.46 1.90 0.94 0.10 0.00 0.00
TOTAL TIME = 12030.4 ! total live time
TOTAL COUNTS = 7704.04 ! total counts
R1 TOTAL = 421.323
R1 THETA = 92.642518628301
R1 TILT = 5.81998E-02
R2 TOTAL = 893.790
R2 THETA = 92.625297174948
R2 TILT = 5.94798E-02
R4 TOTAL = 3759.01 ! counts in R4 band R4 THETA = 92.554456238150 ! rotation angle for gradient R4 TILT = 6.04035E-02 ! gradient
R5 TOTAL = 965.336
R5 THETA = 92.552689049694
R5 TILT = 6.03616E-02
R6 TOTAL = 48.9829
R6 THETA = 92.560059038743
R6 TILT = 6.01161E-02
R7 TOTAL = 0.205248
R7 THETA = 92.651093294314
R7 TILT = 5.72581E-02
>> cast_part.o ! invoke program CAST_PART
Enter observation control
> uma ! observation control is "uma"
OBS =uma
Enter 1 for German data
Enter 2 for US data
Enter 3 for RDF data
2 ! US data set
Enter detector control
1 => PSPC C before 26 January 1991
2 => PSPC B after 26 January 1991
2 ! PSPC B
IDET = 2
Enter gain control
1 => high gain before 11 October 1991
2 => low gain after 11 October 1991
2 ! low gain
IG = 2
Enter the number of AP counts
0 ! zero AP events (we haven't
! fit them yet)
Enter the nominal pointing direction (RA, Dec)
in decimal degrees
172.96 63.84
Selection:
/caldb/data/rosat/pspc/cpf/detp_i_l_b.fits ! detector map
Reading input file
/caldb/data/rosat/pspc/cpf/detmaps/detp_i_l_b.fits
uma_ATTITUDE.FITS
uma_EVENTRATES.FITS
BEGIN = 44435463 END = 44436663
BEGIN = 44441111 END = 44442328
BEGIN = 44491518 END = 44492122
BEGIN = 44492358 END = 44492387
BEGIN = 44508989 END = 44511584
BEGIN = 44521466 END = 44522980
BEGIN = 44527018 END = 44528712
BEGIN = 44577712 END = 44578417
BEGIN = 44578485 END = 44578653
BEGIN = 44607685 END = 44609232
BEGIN = 44609312 END = 44609114
BEGIN = 44671308 END = 44672388
TOTAL NUMBER OF INTERVALS IS 12
1000
2000
3000
4000
5000
6000
7000
8000
9000
10000
11000
12000
File for storing diagnostic information:
uma_diag.out ! diagnostic file name
LIST LENGTH = 444
GOOD STEPS = 12361
BAD STEPS = 2
TOTAL TIME = 12362.0
TOTAL EXP = 12059.0
AVERAGE MV = 80.9570
TOTAL = 1438.70 ! total PB events
TOTAL INT = 954.244 ! internally produced PB events
TOTAL EX1 = 305.833 ! externally produced PB events
! continuum component
TOTAL EX2 = 178.626 ! externally produced PB events
! Al K alpha
>> <cntr> C ! halt program this time
! will come back to it after
! the number of AP counts has
! been determined
>> fit_ap.o ! invoke program FIT_AP
Enter number of PB counts
> 1438.70 ! total counts from first run
! of CAST_PART
Enter number of SX counts
> 7704.04 ! total counts from TILT
Enter IL, IH
> 11 30 ! pulse-height range to be fit
Enter log T for first temperature
> 6.1 ! guess for temperature
23 24 6.05000 6.10000
CHISQ = 37.9, DOF = 16, REDCHISQ = 2.366
PB = 1438.700 AP = 1050.479
Enter log T for first temperature
> 6.09 ! guess for temperature
23 24 6.05000 6.10000
CHISQ = 35.0, DOF = 16, REDCHISQ = 2.185 ! better
PB = 1438.700 AP = 596.087
Enter log T for first temperature
> 6.08 ! guess for temperature
23 24 6.05000 6.10000
CHISQ = 36.3, DOF = 16, REDCHISQ = 2.267 ! overshot minimum
PB = 1438.700 AP = 192.006
Enter log T for first temperature
> 6.09 ! guess for temperature
23 24 6.05000 6.10000
CHISQ = 35.0, DOF = 16, REDCHISQ = 2.185
PB = 1438.700 AP = 596.087 ! best guess for number of AP
! events
>> <cntr> C ! halt program when best
! fit is achieved, plot is found
! in plot_ap.qdp
>> cast_part.o ! reinvocation of program CAST_PART
Enter observation control
> uma ! observation control is "uma"
OBS =uma
Enter 1 for German data
Enter 2 for US data
Enter 3 for RDF data
2 ! US data set
Enter detector control
1 => PSPC C before 26 January 1991
2 => PSPC B after 26 January 1991
2 ! PSPC B
IDET = 2
Enter gain control
1 => high gain before 11 October 1991
2 => low gain after 11 October 1991
2 ! low gain
IG = 2
Enter the number of AP counts
596.087 ! number of AP counts from FIT_AP
Enter the nominal pointing direction (RA, Dec)
in decimal degrees
172.96 63.84
Selection:
/caldb/data/rosat/pspc/cpf/detmaps/detp_ap_l.fits
Reading input file
/caldb/data/rosat/pspc/cpf/detmaps/detp_ap_l.fits
uma_ATTITUDE.FITS
uma_EVENTRATES.FITS
BEGIN = 44435463 END = 44436663
BEGIN = 44441111 END = 44442328
BEGIN = 44491518 END = 44492122
BEGIN = 44492358 END = 44492387
BEGIN = 44508989 END = 44511584
BEGIN = 44521466 END = 44522980
BEGIN = 44527018 END = 44528712
BEGIN = 44577712 END = 44578417
BEGIN = 44578485 END = 44578653
BEGIN = 44607685 END = 44609232
BEGIN = 44609312 END = 44609114
BEGIN = 44671308 END = 44672388
TOTAL NUMBER OF INTERVALS IS 12
1000
2000
3000
4000
5000
6000
7000
8000
9000
10000
11000
12000
File for storing diagnostic information:
uma_diag.out ! diagnostic file name
LIST LENGTH = 444
GOOD STEPS = 12361
BAD STEPS = 2
TOTAL TIME = 12362.0
TOTAL EXP = 12059.0
AVERAGE MV = 80.9570
TOTAL = 1438.70 ! total PB events
TOTAL INT = 954.244 ! internally produced PB events
TOTAL EX1 = 305.833 ! externally produced PB events
! continuum component
TOTAL EX2 = 178.626 ! externally produced PB events
! Al K alpha
PROCESSING 10% COMPLETE ! cast afterpulse component
PROCESSING 20% COMPLETE
PROCESSING 30% COMPLETE
PROCESSING 40% COMPLETE
PROCESSING 50% COMPLETE
PROCESSING 60% COMPLETE
PROCESSING 70% COMPLETE
PROCESSING 80% COMPLETE
PROCESSING 90% COMPLETE
PROCESSING 100% COMPLETE
Writing output file
partap.fits
Selection:
/caldb/data/rosat/pspc/cpf/detmaps/detp_i_l_b.fits ! detector map for
! internally-produced component
! cast internally-produced component
Reading input file
/caldb/data/rosat/pspc/cpf/detmaps/detp_i_l_b.fits
PROCESSING 10% COMPLETE
PROCESSING 20% COMPLETE
PROCESSING 30% COMPLETE
PROCESSING 40% COMPLETE
PROCESSING 50% COMPLETE
PROCESSING 60% COMPLETE
PROCESSING 70% COMPLETE
PROCESSING 80% COMPLETE
PROCESSING 90% COMPLETE
PROCESSING 100% COMPLETE
Writing output file
partin.fits
Selection:
/caldb/data/rosat/pspc/cpf/detmaps/detp_e_l.fits ! detector map for
! externally-produced component
! cast externally-produced
! continuum component
Reading input file
/caldb/data/rosat/pspc/cpf/detmaps/detp_e_l.fits
PROCESSING 10% COMPLETE
PROCESSING 20% COMPLETE
PROCESSING 30% COMPLETE
PROCESSING 40% COMPLETE
PROCESSING 50% COMPLETE
PROCESSING 60% COMPLETE
PROCESSING 70% COMPLETE
PROCESSING 80% COMPLETE
PROCESSING 90% COMPLETE
PROCESSING 100% COMPLETE
Writing output file
partex1.fits
Reading input file
/home/lhea1/snowden/code/ssx/detp_e_l.fits
PROCESSING 10% COMPLETE ! cast externally-produced PROCESSING 20% COMPLETE ! Al K alpha component
PROCESSING 30% COMPLETE
PROCESSING 40% COMPLETE
PROCESSING 50% COMPLETE
PROCESSING 60% COMPLETE
PROCESSING 70% COMPLETE
PROCESSING 80% COMPLETE
PROCESSING 90% COMPLETE
PROCESSING 100% COMPLETE
Writing output file
partex2.fits
>> cast_ssx.o ! invoke program CAST_SSX
Enter the band to be cast
> 4 ! process band R4
expose_4.fits
Reading input file
expose_4.fits
expose_4.fits
Reading input file
expose_4.fits
Enter the total SSX counts, tilt angle, and tilt amplitude
> 3759.01 92.55 0.0604 ! results from program TILT
Enter the total other counts (stuff other than particle and SSX, e.g. LTE)
2785.43 ! results from program LTE
lte_4.fits
Writing output file
lte_4.fits
ssx_4.fits
Writing output file
ssx_4.fits
>> detect.o ! invoke program DETECT
Enter range in bands for the source detection Usually 1(8) - 2 or 4 - 7
> 4 4 ! since we only have the exposure
! map for band 4, we can only use
! band 4
4 4
Weighted Energy = 600.000 ! energy for the PSF calculation
Enter statistical significance and count-rate threshold
> 5 0 ! 5 sigma, o count rate
Statistical significance 5.0000000000000
SIG: 0.99999942657000
b_4.fits ! read in the count image
Reading input file
b_4.fits
0.0000 0.0000 0.0000 1.0000 0.0000 0.0000 0.0000
! band weighting for exposure maps
expose_4.fits ! read in the exposure image
Reading input file
expose_4.fits
EXPAVE = 6928.15 ! average exposure mask_4.fits ! output mask name diag_4.dat ! diagnostic data file
! the following lines occur when a source is detected
ADD! 0.99999999999963 1.1177158007008D-02 124.00000000000 250.0000000
124.00000000000 250.00000000000 1.9927663027031D-02 9.868221083
131.00000000000 6632940.3413086 302.00000000000 30878716.367126
64.871478440333
18 44 1009
ADD! 0.99999999325521 4.2425965930909D-03 286.00000000000 424.0000000
286.00000000000 424.00000000000 6.1922272483099D-03 2.373890413
48.000000000000 1961168.3345947 97.000000000000 10337882.4219971
18.401575940824
9 23 253
ADD! 0.99999999465804 1.5889453804094D-03 366.00000000000 336.0000000
366.00000000000 336.00000000000 1.7703464406142D-03 3.402956364
20.000000000000 553563.96777344 21.000000000000 3023841.3037109
3.8443959704419
4 10 49
ADD! 1.0000000000000 6.9409516810909D-03 392.00000000000 250.0000000
392.00000000000 250.00000000000 8.2636893419773D-03 2.016832993
66.000000000000 1190025.3740234 91.000000000000 6722916.5910645
16.107935829528
7 18 149
ADD! 1.0000000000000 8.5340192436083D-03 424.00000000000 279.0000000
424.00000000000 279.00000000000 1.0530590141792D-02 2.849973023
103.000000000000 2474600.1552734 154.00000000000 13671006.583008
27.875666769540
9 23 253
Writing output file ! write the output file
mask_4.fits
>> final_image.o ! invoke program FINAL_IMAGE
Enter 1 => PCPC C
Enter 2 => PCPC B
> 2 ! PSPC B
Enter 1 => High gain
Enter 2 => Low gain
> 2 ! low gain
Enter the band range
> 4 4 ! cast band R4
PB scale factors: 8.10000E-02 8.40000E-02 0.
b_4.fits ! input count image
Reading input file
b_4.fits
Total counts = 12387.0 ! total number of events
1.00000 ! exposure weightings for the
! bands (only one in this case)
expose_4.fits ! input exposure image
partin.fits
Reading input file
partin.fits
partex1.fits
Reading input file
partex1.fits
! The file partex2 is not included since the Al K alpha contamination ! does not extend down to band R4
! If band R1 or R1L were being processed, the program would query here ! if an afterpulse file should be subtracted
Total modelled PB counts = 102.9804 ! PB events in band R4
Writing output file
part_4.fits ! output PB image
lte_4.fits
Reading input file
lte_4.fits ! input LTE image
ssx_4.fits
Reading input file
ssx_4.fits ! input SSX image
Total modelled background counts = 6650.94 ! total number of background
! events in band R4
Writing output file
back_4.fits ! output background image
Enter 1 to use the exposure map for the mask
> 0 ! use a mask file
Enter mask filename
Hit <CR> only to use the exposure map
> mask_4.fits ! mask to use for band R4,
! the exposure map will include
! all data
mask_4.fits
Reading input file
mask_4.fits
Writing output file
rate_bin_4.fits ! binned rate image
! x10^-6 counts s^-1 arcmin^-2
sig_bin_4.fits ! binned sigma image
radial_4.qdp ! QDP plot of radial profile
73.3528 1.49346 ! average count rate and uncertainty
! for the masked field is
! 73.4 +/- 1.5 x 10^-6
! counts s^-1 arcmin^-2
rate_4.fits ! unbinned rate image
Writing output file
rate_4.fits
sig_4.fits ! unbinned sigma image
Writing output file
sig_4.fits
>> adapt.o ! invoke program ADAPT Enter number of counts for smoothing
> 50 ! use 50 counts as the smoothing
! limit
Enter scale factor for output
> 1 ! allow for scaling of output
! image, units are
! 10e-6 counts/s/arcmin^2
Enter fractional limit for template
> 0.02 ! use only regions of the image
! where the mask value is greater
! than 0.02 of the average value
Enter the binning factor (1, 2, or 4)
> 4 ! bin by 4 in each dimension
Enter the band range to be smoothed
> 4 4 ! process band R4
Reading input file
b_4.fits ! input count file
Reading input file
back_4.fits ! input background file
Reading input file
expose_4.fits ! input exposure file
Enter mask filename
Hit <CR> only to use the exposure map
> mask_4.fits
Reading input file
mask_4.fits
START THE BINNING
START THE SMOOTHING
25
50
75
100
125
Writing output file
adapt_4.fits
*** Program completed successfully ***
>> ls
AO_COL.OUT expose_4.fits sig_4.fits AO_SSX_SC.OUT lte_4.fits sig_bin_4.fits RATE1.DAT mask_4.fits spec.dat RATE2.DAT part_4.fits ssx_4.fits adapt_4.fits partap.fits uma_ATTITUDE.FITS b_1l.fits partex1.fits uma_ATTITUDE1.FITS b_2.fits partex2.fits uma_EVENTRATES.FITS b_3.fits partin.fits uma_EVENTS.FITS b_4.fits plot_4.qdp uma_ORBIT.FITS b_5.fits plot_4_all.qdp uma_diag.out b_6.fits plot_ap.qdp valid_times.dat b_7.fits radial_4.qdp valid_times_all.dat back_4.fits rate_4.fits valid_times_ful.dat diag_4.dat rate_bin_4.fits