c******************************************************************* subroutine showbook c c c c open 6 SHOWER n-tuples with layout: c c IDEVT, NTRK, SHOWTRK SHOWDET, SHOWDE, SHOWX, SHOWY, c ^ ^ ^ ^ ^ ^ ^ c | | | | | | | c ev. ID nb. tracks track nb. det. nb. energy x pos y pos c c ... SHOWTHETA, SHOWPHI, SHOWBETA c ^ ^ ^ c | | | c incidence azimuth beta (p/E) c c last modified on: 06/07/99 by J. Otwinowski c c******************************************************************* implicit none #include "geant321/gcflag.inc" #include "showtups.inc" integer i, ntupid, icycle, err call hbset('BSIZE',1024,err) do i=1,nsect ntupid = 4000+i call hbnt(ntupid,'SHOWER'//char(48+i),' ') call hbname(ntupid,'IEVT',idevt,'IDEVT:U') call hbname(ntupid,'SHOW'//char(48+i),ntrk, + 'NHIT[0,600]:U,'// ! nb. of tracks in this sector + 'SHOWTRK(NHIT)[0,4000]:U,'// ! track number + 'SHOWDET(NHIT)[0,3]:U,'// ! plane number + 'SHOWDE(NHIT),'// ! deposited energy (in MeV) + 'SHOWX(NHIT),'// ! x position (in cm) + 'SHOWY(NHIT),'// ! y position (in cm) + 'SHOWTHETA(NHIT),'// ! angle of incidence (0-180) + 'SHOWPHI(NHIT),'// ! azimuthal angle (0-360) + 'SHOWBETA(NHIT)') ! beta of particle enddo return entry showbooklast ! save ntuples to file do i=1,nsect c call hprint(4000+i) call hrout(4000+i,icycle,' ') enddo return end