***********************************************************************
*
* printcla.kum: printclu modified for "straight" tracks from Sr
* this kumac has been modified by D. Vite on 27/6/97
* last modification: 25/7/97
* please do not modify and contact me directly for any info. Thanks. DV
*
***********************************************************************

*
macro printcla dir=dump file=cluster.lis cut=event>0 printfil=printcl2.ps
* 111macro printclu file=cluster.lis cut=event>0 printfil=printclu.ps
*

exec davmodif


*
fort/file 80 [dir]\[printfil]
* 222fort/file 80 [printfil]
*

 meta 80 -111
*
opt zfl
hist/del 1
*
 nt/crea 1 cluster 8 ! ! event clu start len cog pos charge sovern

*
nt/read 1 [dir]\[file] 'f6.0,1x,f3.0,1x,f4.0,1x,f3.0,1x,f6.2,1x,f6.2,1x,f7.2,1x,f7.2'
* 333nt/read 1 [file] 'f6.0,1x,f3.0,1x,f4.0,1x,f3.0,1x,f6.2,1x,f6.2,1x,f7.2,1x,f7.2'
*


* set htyp -1

*****************
nt/cut 1 len<2.1
nt/cut 2 len<2.1
*****************
* to be modified each time depending on the actual position of the source
* or on the type of cut desired
*****************

*
*
his/crea/1d 10 'Cluster size (S)' 10 0. 10.
nt/pro 10 1.len [cut].and.start<641.and.$1
his/crea/1d 11 'Cluster size (K)' 10 0. 10.
nt/pro 11 1.len [cut].and.start>640.and.$2
*
his/crea/1d 20 'Cluster position' 1024 0. 1024.
nt/pro 20 1.start [cut].and.($1.or.$2)
*
his/crea/1d 30 'Cluster charge (S)' 400 0. 400.
nt/pro 30 1.charge [cut].and.start<641.and.$1
his/crea/1d 31 'Cluster charge (K)' 400 0. 400.
nt/pro 31 1.charge [cut].and.start>640.and.$2
*
his/crea/1d 40 'Cl. Signal/Noise (S)' 40 0. 80.
nt/pro 40 1.sovern [cut].and.start<641.and.$1
his/crea/1d 41 'Cl. Signal/Noise (K)' 40 0. 80.
nt/pro 41 1.sovern [cut].and.start>640.and.$2
*

 option stat
* set csiz 0.43

 zone 2 2
 histo/plot 10
 histo/plot 11
 histo/plot 40
 histo/plot 41

 wait

*
 zone 1 2
 option liny
 histo/plo 20
 option logy
 option nsta
 histo/plot 20
* exec linesvas
 option liny
 option stat
*

 wait

*
 zone 2 2
 option liny
 hist/plot 30
 hist/plot 31
 option logy
 option nsta
 option grid
 histo/plot 30
 histo/plot 31
 option liny
 option ngri
 option stat

   wait

 histo/create/1d 50 'eta function (S)' 100 -0.5 1.5
 histo/create/1d 51 'eta function (K)' 100 -0.5 1.5
 ntuple/pro 50 1.pos-start len=2.and.[cut].and.start<641.and.$1
 ntuple/pro 51 1.pos-start len=2.and.[cut].and.start>640.and.$2
 zone 2 2
 option liny
 histo/plot 50
 histo/plot 51
 option logy
 option nsta
 option grid
 histo/plot 50
 histo/plot 51
 option liny
 option ngri
 option stat
*
* set ltyp 3
* line 0 -10 0 1000
* line 0.5 -10 0.5 1000
* line 1 -10 1 1000
* set ltyp


*

 set htyp 0

*
meta 0
*
close 80
*