next up previous
Next: Definition of Input Files Up: User's Guide to SPECTRAL Previous: Fitting and displaying results

Example

Here is a transcript of a program run. The user's responses are indicated by square brackets. The default file used here is not the same as that shown earlier.

 Reading data...
 Number of class A acceptance cones:           25
 Setting up energy bins...



  EGRET SOURCE OBSERVATION
  Start time:  8- 5-90  6:23:52
  Stop time:   8-19-90  6:23:51
  Source celestial position:
    Right ascension:  200.21
    Declination:       20.28



 Reading photons...
 Reading response matrix...
 Reading background...

 Rescaling bins...
Commands: Fit,Counts,Rebin,Parameters,Write,Display,Top,Spawn,Help,Quit.
Enter command: [counts]




 Bin   High edge   Meas. Counts    Background    Source Estimate
  C    -------        2816           2227.7           588.3
  1        26.         161            137.8            23.2
  2        35.         336            265.6            70.4
  3        46.         672            479.0           193.0
  4        60.         886            653.7           232.3
  5        80.        1016            617.0           399.0
  6       105.         801            455.6           345.4
  7       138.         639            303.2           335.8
  8       182.         481            191.0           290.0
  9       240.         351            116.0           235.0
 10       317.         262             68.2           193.8
 11       418.         199             38.8           160.2
 12       551.         166             21.5           144.5
 13       726.          95             11.7            83.3
 14       957.          71              6.3            64.7
 15      1262.          77              3.4            73.6
 16      1664.          43              1.8            41.2
 17      2193.          27              1.0            26.0
 18      2891.          24              0.5            23.5
 19      3811.          14              0.3            13.7
 20      5024.          12              0.1            11.9
 21      6623.           9              0.1             8.9
 22      8730.           5              0.0             5.0
 23     11509.           5              0.0             5.0
 24     15172.           3              0.0             3.0
 25     20000.           4              0.0             4.0


Commands: Fit,Counts,Rebin,Parameters,Write,Display,Top,Spawn,Help,Quit.
Enter command: [param]

FIT PARAMETERS
  Form of spectrum: straight
  Guess for power-law slope: 2.00
  Number of lines: 0
  Bins used in analysis:   0 - 25 of  25
  Extra iterations: 1
  Automatic normalization


Options: <F>orm,<S>lope,<N>umber of lines,<L>ine energies,<D>ata bins,
  <B>reak energy,<I>terations,<A>utonorm,<E>norm,<H>elp,<Q>uit options.
Choose option-> [quit]
Commands: Fit,Counts,Rebin,Parameters,Write,Display,Top,Spawn,Help,Quit.
Enter command: [fit]
 LEAST-SQUARES FIT TO DATA

 Analyzing data...



ITERATION 0:

Estimate of alpha:   2.0705
Error on alpha:  0.0196

Flux at   227. MeV:  5.317E-09
Error on that flux:  1.185E-10
Flux at 1 MeV (I):  4.021E-04


Correlation coefficient: 0.00

Reduced chi-squared:  0.89
Degrees of freedom: 24
 Hit <cr> when ready. [cr]



ITERATION 1:

Estimate of alpha:   2.0619
Error on alpha:  0.0209

Flux at   234. MeV:  5.059E-09
Error on that flux:  1.214E-10
Flux at 1 MeV (I):  3.871E-04


Correlation coefficient: 0.00

Reduced chi-squared:  1.05
Degrees of freedom: 24
 Hit <cr> when ready. [cr]
Commands: Fit,Counts,Rebin,Parameters,Write,Display,Top,Spawn,Help,Quit.
Enter command: d
  * = measured counts;  + = model counts;  # = both
BIN   ENERGY  SOURCE   MODEL   +---------+---------+---------+---------+
  0      20.  588.33  553.55   |                                      +*
  1      26.   23.17   24.80   |                                 #     |
  2      35.   70.41   64.02   |                                   +*  |
  3      46.  193.05  153.38   |                                      #|
  4      60.  232.30  280.23   |                                      *+
  5      80.  399.01  359.52   |                                       #
  6     105.  345.37  366.28   |                                     *+|
  7     138.  335.80  338.13   |                                    #  |
  8     182.  290.02  296.06   |                                   #   |
  9     240.  234.98  250.15   |                                 #     |
 10     317.  193.81  204.52   |                               #       |
 11     418.  160.18  161.98   |                              #        |
 12     551.  144.48  125.09   |                            #          |
 13     726.   83.28   94.86   |                         *+            |
 14     957.   64.70   71.07   |                       *+              |
 15    1262.   73.63   52.94   |                     + *               |
 16    1664.   41.20   39.51   |                   +*                  |
 17    2193.   26.04   29.49   |                 #                     |
 18    2891.   23.48   21.99   |               #                       |
 19    3811.   13.72   16.41   |            *+                         |
 20    5024.   11.85   12.25   |           #                           |
 21    6623.    8.92    9.17   |         #                             |
 22    8730.    4.96    6.88   |      *+                               |
 23   11509.    4.98    5.00   |     #                                 |
 24   15172.    2.99    3.19   |  #                                    |
 25   20000.    3.99    1.86   +  *                                    |
                               +---------+---------+---------+---------+
 Hit <cr> when ready. [cr]
   * = (OBSERVED COUNTS - MODEL) normalized to unit variance
BIN   ENERGY     CHI   CHI^2   +---------+---------+---------+---------+
  0      20.    0.66    0.43   |                   |    *              |
  1      26.   -0.13    0.02   |                  *|                   |
  2      35.    0.35    0.12   |                   |  *                |
  3      46.    1.58    2.49   |                   |          *        |
  4      60.   -1.57    2.46   |        *          |                   |
  5      80.    1.26    1.60   |                   |        *          |
  6     105.   -0.73    0.53   |              *    |                   |
  7     138.   -0.09    0.01   |                  *|                   |
  8     182.   -0.27    0.08   |                 * |                   |
  9     240.   -0.80    0.63   |             *     |                   |
 10     317.   -0.65    0.42   |              *    |                   |
 11     418.   -0.13    0.02   |                  *|                   |
 12     551.    1.61    2.60   |                   |           *       |
 13     726.   -1.13    1.28   |           *       |                   |
 14     957.   -0.73    0.53   |              *    |                   |
 15    1262.    2.79    7.76   |                   |                   *
 16    1664.    0.27    0.07   |                   | *                 |
 17    2193.   -0.63    0.40   |              *    |                   |
 18    2891.    0.32    0.10   |                   | *                 |
 19    3811.   -0.67    0.44   |              *    |                   |
 20    5024.   -0.12    0.01   |                  *|                   |
 21    6623.   -0.08    0.01   |                  *|                   |
 22    8730.   -0.75    0.56   |              *    |                   |
 23   11509.   -0.01    0.00   |                   *                   |
 24   15172.   -0.11    0.01   |                  *|                   |
 25   20000.    1.60    2.55   |                   |          *        |
                               +---------+---------+---------+---------+
Commands: Fit,Counts,Rebin,Parameters,Write,Display,Top,Spawn,Help,Quit.
Enter command: [rebin]
<R>escale, <C>ombine, or <Q>uit? [combine]
First bin to combine: [24]
Last bin to combine: [25]
 Combining bins...
<R>escale, <C>ombine, or <Q>uit? [c]
First bin to combine: [22]
Last bin to combine: [23]
 Combining bins...
<R>escale, <C>ombine, or <Q>uit? [q]
Commands: Fit,Counts,Rebin,Parameters,Write,Display,Top,Spawn,Help,Quit.
Enter command: [p]

FIT PARAMETERS
  Form of spectrum: straight
  Guess for power-law slope: 2.06
  Number of lines: 0
  Bins used in analysis:   0 - 23 of  23
  Extra iterations: 1
  Automatic normalization


Options: <F>orm,<S>lope,<N>umber of lines,<L>ine energies,<D>ata bins,
  <B>reak energy,<I>terations,<A>utonorm,<E>norm,<H>elp,<Q>uit options.
Choose option-> [n]
Number of lines (present value = 0)
-->[1]
Choose option-> [l]
Energy of line 1 -->[1200]
Choose option-> [auto]
User-specified normalization.
Choose option-> [enorm]
Normalization energy (present value =  234.)
 -->[200]
Choose option-> [q]
Commands: Fit,Counts,Rebin,Parameters,Write,Display,Top,Spawn,Help,Quit.
Enter command: [f]
 LEAST-SQUARES FIT TO DATA

 Analyzing data...



ITERATION 0:

Estimate of alpha:   2.0884
Error on alpha:  0.0202

Flux at   200. MeV:  6.837E-09
Error on that flux:  1.472E-10
Flux at 1 MeV (I):  4.369E-04

Strength of line at  1200. MeV:  3.558E-08
Error on line strength:  1.440E-08

Correlation coefficient:-0.03

Reduced chi-squared:  0.75
Degrees of freedom: 21
 Hit <cr> when ready. [cr]



ITERATION 1:

Estimate of alpha:   2.0843
Error on alpha:  0.0205

Flux at   200. MeV:  6.870E-09
Error on that flux:  1.518E-10
Flux at 1 MeV (I):  4.294E-04

Strength of line at  1200. MeV:  3.614E-08
Error on line strength:  1.463E-08

Correlation coefficient:-0.02

Reduced chi-squared:  0.80
Degrees of freedom: 21
 Hit <cr> when ready. [cr]
Commands: Fit,Counts,Rebin,Parameters,Write,Display,Top,Spawn,Help,Quit.
Enter command: [write]
Commands: Fit,Counts,Rebin,Parameters,Write,Display,Top,Spawn,Help,Quit.
Enter command: [top]
Commands: Fit,Counts,Rebin,Parameters,Write,Display,Top,Spawn,Help,Quit.
Enter command: [quit]
FORTRAN STOP



CGRO SSC
1998-06-29