                     -----------------------------------

                               IMP-8 ANALIMP-F

                                 User's Guide

                     -----------------------------------

1. GENERAL INFORMATION

   The IMP-8 ANALIMPF program is to analyze pha data from PHA files. This
   program produces A vs B or D vs E plots of IMP-8 pha data, compute the
   distribution of particles about a centroidal curve, and prints
   histogram to reflect this distribution.

2. INPUT PARAMETER CARDS

   Card 1 (Event Identification Card)

   Name     Column   Format  Description

   ISPEC     1        I1     Histogram Option:
                             =0, compute histogram (default)
                             =1, no histograms

   KGAIN     3        I1     Gain Indicator
                             =0, High Gain (default)
                             =1, Low Gain

   ICMP      5        I2     Compression Factor
                             =1, first 128 channels on each axis
                                 of matrix with no compression
                             =2, first 256 channels on each axis
                                 compressed to a 128xl28 matrix
                             =4, first 512 channels on each axis
                                 compressed to a 128xl28 matrix
                             =8, all 1024 channels on each axis
                                 compressed to a 128xl28 matrix

   IEV       7        I1     Event to be analyzed
                             =1, A.B.-C and A.-B.-C if KGAIN=0;
                                 (A+B)1.A.B.-C and (A+B)1.A.-B.-C
                                 if KGAIN=1
                             =2, (A+B)2.A.B.-C and (A+B)2.A.-B.-C
                                 (KGAIN=0 only)
                             =3, D.E.-F.-G (KGAIN=0 only)
                             =4, (D1+E1)1.E.-F.-G (KGAIN=0 only)
                             =5, (D1+E1)2.E.-F.-G Low (10) (KGAIN=1 only)
                             =6, (D1+E1)2.E.-F.-G Low (50) (KGAIN=1 only)

   LDMASK   9-11      I3     LED T2T3T4 selection mask
                             =I2I3I4 where

                             Ij=0, then the LED Tj event type bit must be 0
                                   for an event to be accumulated
                             Ij=1, then the LED Tj event type bit must be 1
                                   for an event to be accumulated
                             Ij=2, then the LED Tj event type bit is not
                                   tested in the accumulation decision
                                   (LDMASK is a requisite only if IEV=1 or 2)

   STD(1)  13-22     F10.0   If IEV less than 3, then A standard endpoint
                             for gain factor computation. Otherwise, D
                             standard endpoint for gain factor computation.

   STD(2)  24-33     F10.0   If IEV less than 3, then B standard endpoint
                             for gain factor computation. Otherwise, E
                             standard endpoint for gain factor computation.

   If ISPEC = 1, cards 2 through 14 are omitted. Time cards (card 15)
   follow immediately after card 1.

   Card 2

   Name     Column   Format  Description

   COM       1-16     A16    A user comment - written on page 1 of
                             output listing

   NZ       19-20     I2     Particle atomic number

   XA       21-30     F10.0  Particle atomic mass (in proton units)

   IDIST     32       I1     Distance selector
                             =0, Perpendicular distance
                             =1, Vertical distance

   IPRINT    34       I1     Debug printout option
                             =0, no printout
                             =1, diagnostic printout requested

   NCHOF    36-40     I5     Number of channels offset
                             (Center of histograms offset from curve)

   MXNCL    41-45     I5     Maximum number of iterations allowed to
                             compute perpendicular distance of any
                             given point to event line

   Card 3 (Scan Parameters Card)

   Name     Column   Format  Description

   NLX       1-5      I5     X channel for lower horizontal scanning line

   NUX       6-10     I5     X channel for upper horizontal scanning line

   NLY      16-20     I5     Y channel for left vertical scanning line

   NUY      21-25     I5     Y channel for right vertical scanning line

   NIX1     31-35     I5     X coordinate of point (X1,0)

   NIY1     36-40     I5     Y coordinate of point (0,Y1)
                             Points (X1,0) and (0,Y1) are connected
                             to form lower scanning line

   NIX2     46-50     I5     X coordinate of point (X2,0)

   NIY2     51-55     I5     Y coordinate of point (0,Y2)
                             Points (Y2,0) and (0,Y2) are connected
                             to form upper scanning line

   Card 4 (Histogram Parameters Card)

   Name     Column   Format  Description

   NB        1-2      I2     Number of histograms (also number of
                             histogram energy limit cards to follow)
                             maximum = 10

   CHW       3-10     F8.0   Histogram bin width - in channels

   Cards 5 through N (where N = 4 + NB of card 4, maximum N = 14)

   Name     Column   Format  Description

                             One card for each histogram (maximum of 10)

   ELOW      1-10     F10.0  Lower energy limit for histogram

   EUP      11-20     F10.0  Upper energy limit for histogram

   Time Indicator Card

   Name     Column   Format  Description

   IOPT       1       I1     0= Do not plot, accumulate for the time
                                period on this card (default)
                             1= Plot all accumulated data.

   IORB1     3-7      I5     Start orbit.

   IORB2     8-12     I5     Stop orbit.

   IYR1     13-16     I4     Start year.

   IDAY1    18-20     I3     Start Julian day.

   IHR1     22-23     I2     Start hour.

   IMIN1    25-26     I2     Start minute.

   IYR2     28-31     I4     Stop year.

   IDAY2    33-35     I3     Stop Julian d,

   IHR2     37-38     I2     Stop hour.

   IMIN2    40-41     I2     Stop minute.

   END1     60-64     F5.0   Endpoint used in computing the gain shift
                             along the ordinates
                             see Note 2

   END2     66-70     F5.0   Endpoint used in computing the gain shift
                             along the abscissa
                             see Note 3

   GTABLE    72       I1     Gain factor table switch
                             see Note 4

   If IOPT=1 in the time indicator card, then omit the pick option card

   Pick option Card

   Name     Column   Format  Description

   CHOICE    1-4      A4     Selection flag
                             =PICK (default)

   QPHA       6       L1     =true  the program is to use only PHA files
                                    (KGAIN must equal to zero)
                                    currently we restrict the program
                                    using PHA files only
                             =false the program can use PHA, MATR or LOWG
                                    files

   QFINE      8       L1     =true  finegain factors are to be used if they
                                    are available
                             =false access the gain factors indicated by
                                    variable GTABLE

----------------------------------------------------------------------------

Note 1 :

(a) If an MED plot is requested, then this is the D detector endpoint.
(b) If an LED plot is requested, then this is the A detector endpoint.

Note 2 :

(a) If an MED plot is requested, then this is the E detector endpoint.
(b) If an LED plot is requested, then this is the B detector endpoint.

Note 3 :

GTABLE = 0: A, B events must use only this option gains are calculated
from card input.

(a) If STD values are zero, the program defaults to the internal A-B
    or D-E endpoints and uses them along with the END values to calculate
    gain factors.
(b) If END1 and END2 are zero, the program will bypass gain correction
    (effectively gains = 1)
(c) If the STD and END values do not give reasonable gain factors
    (STD(I)/ END(I) <= 3) the gain factors are set to 1.0.

GTABLE = 1: D-E events only.  Neither STD nor END values need be coded
on the time card.  The gain tables are accessed for each interval in the
time span/orbit(s) requested.

GTABLE = 2: D-E events only.  Neither STD nor END values need be coded
on the time card.  The gain tables are accessed for gain factors, and
the data will be normalized to IMP-6 data.

Description of Parameter Cards

All data specified on cards 1 through 14 is held constant for the entire run.
Invalid combinations for gain and event on card 1 will cause the job to be
terminated.

Cards 2 through 14 pertain to histograms. If no histograms are required
(ISPEC = 1), they must be omitted.

Figure 1 contains a diagram of the scanning box derived from the data on
card 3. All values specified on card 3 should be greated than or equal
to zero.

On card 4, an NB of less than one or greater than ten will cause job
termination.  CHW must be greater than zero and, unless specifically
requested otherwise, should be the same value as the compression factor
(ICMP on card 1).

Card 15 (Time card) may be repeated as many times as are necessary to
analyze all requested data.  If partial orbits are requested, the time cards
should be in ascending chronological order. When IOPT = 1, all other data
on that card is ignored. For all further discussion of this card, it is
assumed that IOPT = 0.

By specifying IORB1=0, IORB2 > 0 a time span can be broken up into time
increments of IORB2 (hours). For each increment, data is accumulated and
a plot (with histogram, if requested) is printed.  Increment processing
stops when the start time of an increment is larger than the stop time on
the input data card. If IORB2 is zero, the entire time period is processed.


3. PROCESSING ENVIRONMENT

   The IMP-8 Analimp-f program (ANALIMPF) is developed on an Sun worksation
   Sun 4.1.4 operating system. The source codes are written in Sun FORTRAN
   and C languages. A local I/O system FTIO package was developed. A local
   shell program was developed to execute the IMP-8 ANALIMPF program.
