Example of reading monitor rate data for an interval:


    IDL> o=hsi_monitor_rate()

    IDL> d=o->getdata(obs_time_interval=['2002-Aug-18 13:00:00','2002-Aug-18 14:30:00'])

    IDL> s=o->get()

    IDL> help, d, /str

     ** Structure HSI_MONITOR_RATE_CYCLE, 8 tags, length=428, data length=428:

          TIME FLOAT 0.00000 ;relative from start of interval

          PARTICLE_LO LONG Array[8] ;particle counter sampled at 8 Hz

          PARTICLE_HI LONG Array[8] ;p.c. higher-energy threshold

          PREAMP_RESET LONG Array[18] ;CSA reset events, 18 segments, 1 Hz

          SHAPER_VALID LONG Array[18] ;Slow-channel (spectroscopy) events, 1 Hz

          SHAPER_OVER_ULD LONG Array[18] ;Cosmic ray counts, 1 Hz

          DELAY_LINE_VALID

          LONG Array[18] ;Fast channel counts, 1 Hz

          LIVE_TIME FLOAT Array[18] ;Uncorrected livetime, scale of 0.0-1.0

;Generate an absolute time axis for 1 Hz rates:

    IDL> t = s.mon_ut_ref + d.time

;Plot rear G8 good event rate vs. time

    IDL> utplot,t,d.shaper_valid[16],yrange=[0,1500]

;Generate an absolute time axis for 8 Hz particle detector rates:

    IDL> t8 = dblarr(n_elements(d)*8)

    IDL> for i=0l,n_elements(d)-1l do for j=0l,7l do t8[i*8l+j] = t[i] + double(j)*(1.d / 8.d)

;Generate a continuous list of particle detector count rates:

    IDL> pd_lo = lonarr(n_elements(d)*8)

    IDL> for i=0l,n_elements(d)-1l do for j=0l,7l do pd_lo[i*8l+j] = d[i].particle_lo[j]

    IDL> pd_hi = lonarr(n_elements(d)*8)

    IDL> for i=0l,n_elements(d)-1l do for j=0l,7l do pd_hi[i*8l+j] = d[i].particle_hi[j]

;Look at particle detector rates (here we see a minor precipitation event at the center of the scale and the SAA on the right)

    IDL> utplot,t8,pd_lo,yrange=[.5,500],/ylog