MACRO figure7
*
* This is an example how to read the extract_log - output file
*    extract_log_500-833. 
*
* It generates figure7 of the report
*
*
* Ntuple variables:   
* 
*    For the acsii file extarct_log_500-833.out, 
*       the rows are just r1,r2, ..., a1...b1,,,x1... . 
*
*    Also, see variable names with extract6.cpp and the report.ps.
*
*     r1   =   run_year      year
*     r2   =   run_day       day
*     r3   =   run_no        run number 
*     r4   =   run_file      file number in current run
*
*     a1   =   A_total       total number of Am-A events
*     a2   =   A_1time       File start time (1st event),   GPS-secs
*     a3   =   A_detim       File duration,                 GPS secs
*     a4   =   A_spase       number of Am-A Spase trigger
*     a5   =   A_resyn       number of resynched events  in Am-A 
*
*     b1   =   B_total       total number of Am-B events
*     b2   =   B_1time       File start tike (1st event),   GPS-secs
*     b3   =   B_detim       File duration,                 GPS secs
*     b4   =   B_spase       number of Am-B Spase trigger
*     b5   =   B_resyn       number of resynched events  in Am-B 
*
*     x1   =   rawcoin       raw coincidences as found by extraction program
*     x2   =   coin          coincidences
*     x3   =   updwncnt      up-down filter events
*     x4   =   evtswdp       events with duplicate channels
*     x5   =   overlap       overlap time of Am-A and Am-B file
*     x6   =   gaspcnt       number of Am-B GASP coincidences
*     c    =   check         quality flag  (not uesd currently)


  if [1].ne.1 then
    nt/cre 1 'a' 21 ' ' 2000 r1 r2 r3 r4 a1 a2 a3 a4 a5 b1 b2 b3 b4 b5 x1 x2 x3 x4 x5 x6 c
    nt/rea 1 extract_log_500-833.out  ! ! ! -/*/
  endif 
  

* if you want to plot...  
*
*   for/file 66 figure7.eps
*   meta   66 -113

  opt *
  title 'Amanda-B trigger rate vs. time'
  set stat 000010
  opt tic   
  opt stat 

  2dhisto 2000 '' 1151 7500000 19000000 250 0 250  
  ntuple/proj 2000 1.b1/b3%b2 b3>0
  histo/plo 2000
  atit 'time in GPS seconds                  .'  'B-rate (Hz)'

* if you want to plot...  
*
*   close 66
  histo/delete 2000
  set *
  opt *

RETURN