macro pall 
*for plotting CLF energy monitors  (all data on different time scales)
opt *
set *
nt/create 1 energy 23 ! ! spare1 n1 n2 emon rms emon2 srms rat rrms vset q azi elv volts temp wind y m dy hr mn s ns   
nt/read 1 temp.dat

*macro/trace on wait
exec pall#plot all.vstats.clf.eps 2003.5
exec pall#plot most.vstats.clf.eps 2006.6
return

macro plot [yr] [mo]
opt *
set *
set gsiz 0.45
set asiz 0.35
set vsiz 0.3
set ndvx 510
set ndvy 505
opt grid
set *fon 60
set ygti 0.8
set mscf 0.4
set ywin 2.5
set ymgu 0.5
set yval 0.3
set tsiz 0.35
set xmgr 0.5
set xwin 2.5
rat=1.17e+6

day=(dy-1+(hr+(mn+s/60)/60)/24)
year=(y+(dy/31+m-1)/12)

cuts 1 (y>0)
cuts 2 (m>0)

cuts 3 $1.and.$2
nt/pl 1.[year] $3 ! ! ! N
nn=$hinfo(1000000,'events')
kuip/message 'CLF energy probe data for '//[yr][mo]//' '//[nn]//' entries'
if ([nn]<1) then
exitm
endif
min=$hinfo(1000000,'xmin')
max=$hinfo(1000000,'xmax')
min=[min]+0.3
max=[max]-0.25

*title 'Pierre Auger CLF  Beam Monitoring'
fort/file 44 [1]
xval=year
min=[2]
meta 44 -111

view/zone 2 2
opt grid

gra/att/col 1 0 0 0
set yhti 1.1
hist/delete 201
hist/create/2d 201  'Monitor Energy ' 10 [min] [max] 10 0 1.3

set mtyp 25
nt/pl 1.emon*1e+9%[[xval]] $3 ! ! ! ! 201
atitle [xval] 'Energy'
gra/hplot/key [min]+0.05 0.75 25 'Probe1(photodiode) pJ'
gra/att/col 1 0 0 1
set mtyp 23
nt/pl 1.emon2*1e+3%[[xval]] $3 ! ! ! s
gra/hplot/key [min]+0.05 0.68 23 'Probe2(pyroelectric) mJ'
gra/att/col 1 0 0 0

hist/delete 101
hist/create/2d 101 'Monitor Energy RMS' 10 [min] [max] 10 0.0 25
set mtyp 25
nt/pl 1.rms%[[xval]] $3 ! ! ! ! 101
gra/hplot/key [min]+0.05 23.5 25 Probe1
gra/att/col 1 0 0 1
set mtyp 23
nt/pl 1.srms%[[xval]] $3 ! ! ! s
gra/hplot/key [min]+0.05 21.0 23 Probe2
gra/att/col 1 0.0 .0 .0
atitle [xval] 'Percent'

gra/att/col 1 0 0 0
hist/delete 202
set yhti 1.1
hist/create/2d 202 'Ratio (Normalized) ' 10 [min] [max] 10 0.6 1.2
set yval 0.3
*view/zone 1 4 3 s 
set ywin 0.6
set mtyp 4
nt/pl 1.(emon2/emon/[rat])%[[xval]] $3 ! ! ! ! 202
gra/att/col 1 1 0 0
set mtyp 8
nt/pl 1.(emon2/emon/[rat])%[[xval]] $3 ! ! ! s
gra/att/col 1 0 0 0
atitle [xval] 'E2/E1/'//[rat] 

hist/delete 101
hist/create/2d 101 'RMS of Ratio ' 10 [min] [max] 10 0.0 25
set mtyp 4
nt/pl 1.rrms%[[xval]] $3 ! ! ! ! 101
gra/hplot/key 1.0 8.0 4
gra/att/col 1 1 0 0
set mtyp 8
nt/pl 1.rrms%[[xval]] $3 ! ! ! s
gra/att/col 1 1 0 0
gra/hplot/key 1.0 8.0 8 RMS-E2/E1
gra/att/col 1 0 0 0
atitle [xval] 'Percent'

close 44
return