**********************************************************
*
* printcal.kum
* this kumac has been modified 8/7/97 by D. Vite
* please do not modify without consulting me. Thanks.
*
***********************************************************
 
* 
macro printcal file=calib\calib.cal printfile=dump\printcal.ps
* 111macro printcal file=calib.cal printfile=printcal.ps
* 
  
  
 exec davmodif
 
 
* 
for/file 80 [printfile]
meta 80 -111
*
* set stat 1110
*
vect/crea nva(1) r
vect/crea nchannel(1) r
*
nt/crea 1 calib 5 ! ! channel pedestal sigma gain status
vect/read nva,nchannel,name [file] 'f2.0,1x,f3.0,1x,f10.0' 'o'
k=nva(1)
*
vect/crea cns([k])
vect/read cns [file] 'f12.4' ' '
vect/pri nva
vect/pri nchannel
vect/pri name
*
kname=name(1)
title ladder:temp.sens.no..[kname]
* 
nt/read 1 [file] 'f5.0,1x,f12.4,1x,f12.4,1x,f12.4,1x,f2.0' ' '
* this command reopens the file so the first 17 lines are junk
* it is fortunate that they can read them at all
*
*zone 1 2 
* 
igset mtyp 29 
igset mscf 0.5
*set htyp -1

* 
* 
zone 2 2
nt/plot 1.pedestal channel<641 ! 18
nt/plot 1.pedestal channel>640 ! 18
zone 1 2 2 s
* 
* option nsta
 set stat 110000
nt/plot 1.pedestal%channel  ! ! 18
* exe linesvas 
 set stat 111111
* option stat
*
*zone 2 2
*
wait
*
zone 2 2
ntuple/plot 1.sigma sigma<20.and.channel<641 ! 18
ntuple/plot 1.sigma sigma<20.and.channel>640 ! 18
zone 1 2 2 s
* 
* option nsta
 set stat 110000
nt/plo 1.sigma%channel sigma<20  ! 18
* exe linesvas 
 set stat 111111
* option stat
* 
*zone 2 2
* 
wait
* 
* 
zone 2 2
nt/plot 1.gain channel<641 ! 18
nt/plot 1.gain channel>640 ! 18
zone 1 2 2 s
* 
* option nsta
 set stat 110000
nt/plot 1.gain%channel    ! ! 18
* exe linesvas 
 set stat 111111
* option stat
* 
* 
wait
* 
* 
zone 2 2
nt/plot 1.status channel<641 ! 18
nt/plot 1.status channel>640 ! 18
zone 1 2 2 s
* 
* option nsta
 set stat 110000
nt/plot 1.status%channel ! ! 18
* exe linesvas 
 set stat 111111
* option stat
* 
 
 set htyp 0
 
 meta 0
close 80
* 
*