;Plot Parameters vs. Time

;hsi_read_areas reads ascii files writen by the HESSI GUI containing the areas for different times
;This program plots those areas vs. time for multiple files.

checkvar, ps, 0
initialize, ps

;One plot per page [no. empty sectors, no. in X/Y/Z direction, (0 for left to right, 1 for top to bottom)]
!p.multi = [0,1,2]

path = 'C:\Documents and Settings\lhaga\My Documents\2003 October 28\Flux output files\'
;ascii file containing the area information
;file ='flux_output_april21_6-12_c50.txt'
file = path + 'flux_output_oct28_12-25.txt'

hsi_read_areas,file,tarr,utb,areas

; Time range to plot
;s = 0 & e = n_elements(tarr) - 1
divisor = 1.e19

utplot, tarr-utb, areas/divisor, utb,  title='Area vs. Time', $
	background=0, color=255, $
	xticklen = 1, xtickname = [' ',' ',' ',' ',' ',' ',' ',' ',' ',' ',' '], xtitle = ' ',  $
	ymargin = [0,3], yticklen = 1,$
	yrange=[0.5,1.2],ytitle = 'Area in 10!u19!n cm!u2!n'

outplot,tarr-utb, areas/divisor, utb, color=2
;position=[0.15,0.10,0.95,0.90]
xyouts,0.25,0.85,'12 - 25 keV',color=2,/normal


;Plot quicklook light curves

obj = hsi_obs_summary()
obj -> set, obs_time_interval=['28-oct-03 11:06', '28-oct-03 11:30']

rates_data = obj -> getdata()

corr_rates_data = obj -> getdata(/corrected)

times_rate = obj -> getaxis(/ut)

;To plot the observing summary rates summed over energy, showing the flare, night, and SAA flags, the command is:
;obj -> plot, /saa, /flare, /night, /dim1_sum

;The following example plots the 6-12 keV band of corrected count rates:

;obj =hsi_obs_summary(obs_time= ['23-jul-2002 00:16','23-jul-2002 01:15'])
data = obj -> getdata(/corrected)
time = obj -> getdata(/time)
;utplot, time - time[0], data.countrate[1,*], time[0]
;or
factor = [1.,1.,1.,1.,1.,1.,1.]
addition = [0.,0.,0.,25.,60.,10.,10.]

utplot, tarr-utb,areas/divisor, utb,/ylog,yrange=[1.e-1,1.e5], $
	xticklen = 1, yticklen= 1, ymargin = [3,0], $
	ytitle = 'Corrected count rate in counts s!u-1!n detector!u-1!n'
for i = 0,4 do begin
;	addition = data.countrate(i,0)
	outplot,atime(time), (data.countrate[i,*] - addition(i))*factor(i),color=2*i
	endfor
;outplot,atime(time), data.countrate[1,*],color=2
;outplot,atime(time),data.countrate[2,*],color=4
;outplot,atime(time),data.countrate[3,*],color=6
;outplot,atime(time),data.countrate[4,*],color=8

;The first utplot example uses time in seconds relative to a base time
;(base time provided as the third argument).
;The second utplot example uses fully qualified time in ASCII.

if ps gt 0 then device,/close
set_plot,'win'

end