The program select

The program select, introduced in section 2.6, selectively extracts a Ginga cube, time series or spectral file from an existing data cube. Any number of housekeeping (HK) parameters can be selected upon inside or outside a specified window. Thus, if you applied some ``loose'' selection criteria during the initial data extraction using sortac, you can "tighten" those loose limits with select (but see warning in section 2.8).

1. Extract a cube from another cube.

   qcl>select
    
   ** SELECT Version 2.4 Thu Sep  4 16:06:06 1997 **
    
   SEL_DATIN  Input data filename =mkn335_91a.sub_corr
   SEL_OPT  Extract  a TIME / PHAD / CUBE  file? = PHAD =cube
   SEL_DATOUT  Output data filename = mkn335_91a.cube =
    
   SEL_NPHA  How many PHA ranges? = 48 =
   SEL_NUM_SEL  Number of HK parameters to select on  (-1 = default list)
   = 1 =
    
   SEL_HK_NAME  Name of HK parameter =TRANSMISSION
   SEL_WINDOW  Normal or inverse  window = Norm =
   SEL_PAR_LO  Lower limit on TRANSMISSION = 0.0 =0.7
   SEL_PAR_HI  Upper limit on TRANSMISSION = 0.0 =1.0
    
   SEL_NUM_PER  Number of time windows = 1 =
   SEL_ST_PER001  Start time of window - secs from begining = 0.0 =
   SEL_END_PER001  End time of window - secs from begining = 50304.0 =
    
        DATIN  mkn335_91a.sub_corr
        Shape :      393 x    48   x     2  =    37728
                  TIME    PHA     IDENT   
    
        Time binning is 128.0000 secs
        ... searching    393 time bins from range of    393
    
         ID_MASK( 2)
        0101010101010101
        1010101010101010
    
    
        OUTPUT  mkn335_91a.cube
        Shape :      393 x    48   x     2  =    37728
                  TIME    PHA     IDENT   
    
            88 time bins are occupied 
   qcl>

   The program select prompts for the number of spectral channels in the output cube (the default is the maximum number of channels available), the number and names of HK parameters to select on, and the number, start period and stop period of time windows. If the number of spectral channels is less than the maximum number, then select will prompt for the lower and upper spectral channel contributing to each output channel. The start and stop times are expressed in seconds.

2. Extract a spectral file from an entire cube, or selected time period. This is the same as option 3 in section  5.1 except that HK parameters can be selected upon inside or outside a specified window.

   qcl>select
   
   ** SELECT Version 2.4 Mon Aug 11 14:16:18 1997 **
    
   SEL_DATIN  Input data filename =mkn335_91a.sub_corr
   SEL_OPT  Extract  a TIME / PHAD / CUBE  file? = PHAD =
   SEL_DATOUT  Output data filename = mkn335_91a.phad =
   SEL_NUM_SEL  Number of HK parameters to select on  (-1 = default list)
   = 1 =
    
   SEL_HK_NAME  Name of HK parameter =sud
   SEL_WINDOW  Normal or inverse  window = Norm =
   SEL_PAR_LO  Lower limit on SUD = 0.0 =
   SEL_PAR_HI  Upper limit on SUD = 0.0 =8.0
   
   SEL_NUM_PER  Number of time windows = 1 =
   SEL_ST_PER001  Start time of window - secs from begining = 0.0 =
   SEL_END_PER001  End time of window - secs from begining = 50304.0 =
    
        DATIN  mkn335_91a.sub_corr
        Shape :      393 x    48   x     2  =    37728
                  TIME    PHA     IDENT   
    
        Time binning is 128.0000 secs
        ... searching    393 time bins from range of    393
    
         ID_MASK( 2)
        0101010101010101
        1010101010101010
    
    
        OUTPUT  mkn335_91a.phad
        Shape :        1 x    48   x     2  =       96
                  TIME    PHA     IDENT   
    
            88 time bins contributed to the output spectrum
   qcl>

3. Extract a time series, covering a specified range of spectral channels, from an existing cube. This is the same as option 4 in section 5.1 except that HK parameters can be selected upon inside or outside a specified window.

   qcl>select
    
   ** SELECT Version 2.4 Thu Sep  4 16:15:25 1997 **
    
   SEL_DATIN  Input data filename =mkn335_91a.sub_corr
   SEL_OPT  Extract  a TIME / PHAD / CUBE  file? = PHAD =time
   SEL_DATOUT  Output data filename = mkn335_91a.time =
   SEL_PH_LO001  Low PHA channel of range = 1 =4
   SEL_PH_HI001  High PHA channel of range = 48 =30
   SEL_NUM_SEL  Number of HK parameters to select on  (-1 = default list)
   = 1 =
    
   SEL_HK_NAME  Name of HK parameter =sud
   SEL_WINDOW  Normal or inverse  window = Norm =
   SEL_PAR_LO  Lower limit on SUD = 0.0 =
   SEL_PAR_HI  Upper limit on SUD = 0.0 =8.0
    
   SEL_NUM_PER  Number of time windows = 1 =
   SEL_ST_PER001  Start time of window - secs from begining = 0.0 =
   SEL_END_PER001  End time of window - secs from begining = 50304.0 =
    
        DATIN  mkn335_91a.sub_corr
        Shape :      393 x    48   x     2  =    37728
                  TIME    PHA     IDENT   
    
        Time binning is 128.0000 secs
        ... searching    393 time bins from range of    393
    
         ID_MASK( 2)
        0101010101010101
        1010101010101010
    
    
        OUTPUT  mkn335_91a.time
        Shape :      393 x     1   x     2  =      786
                  TIME    PHA     IDENT   
    
            88 time bins are occupied 
   qcl>