macro dsm_sorter

*kumac to sort dsm data from akio's trg2ntp program  which converts
*the dsm info in a raw daq file to an ntuple.  The format is not
*to our liking so we will process further here

*do the following file
*hi/file 1 data/runst_physics_4076012_raw_0020001.daq.rzd
*hi/file 1 data/runst_physics_4079058_raw_0010001.daq.rzd


  rzd=data/runst_physics_4079058_raw_001000
  do j=1,4
    r=[j]
    run=[rzd][j].daq.rzd
    message j=[j] r=400[r]  run=[run]
    chain c[j] [run]
  enddo
  chain all c1 c2 c3 c4 
  zebra/rz/cd //all

*some test histos to see we got something
zone 2 2
* data from any input and any bXing
n/pl 111.eemclayer0 eemclayer0.gt.1

* data from  input=75 and any bXing
*n/pl 111.eemclayer0(75) eemclayer0.gt.1

* data from  input=75 and trig bXing
*n/pl 111.eemclayer0(75) eemclayer0.gt.1.and.ipre.eq.0 

* data from  input=75 and bXing=-1
*n/pl 111.eemclayer0(75) eemclayer0.gt.1.and.ipre.eq.1 

hi/crea/1d 1 'dummy' 100 0 100 0 0 0 
hi/create/2d 10 'high tower' 65 0 65 100 0 100
hi/create/2d 12 'trig patch' 65 0 65 100 0 100
hi/create/2d 20 'high tower masked' 66 0 66 100 0 100
hi/create/2d 22 'trig patch masked' 66 0 66 100 0 100

*dummy line to call the fortran program 
nt/proj 1 111.BBC(1) dsmsort.f    

zone 1 2
opt logz
hi/pl 10(,0:57) col
hi/pl 12(,0:57) col
return