*------------------------------------------------------------------------------ c SUBROUTINE EMCBOOK c c open 6 EMC n-tuples with layout: c c IDEVT, NTRKS, EMCTRK EMCDET, EMCE, EMCX, EMCY, EMCTOF c ^ ^ ^ ^ ^ ^ ^ ^ c | | | | | | | | c ev. ID nb. tracks track nb. det. nb. energy posx posy tof c c c last modified on: 06/10/2016 by R.Holzmann c *------------------------------------------------------------------------------ IMPLICIT NONE #include "emctuple.inc" #include "geant321/gcflag.inc" INTEGER I, NTUPID, ICYCLE, ERR, ISTAT CALL HBSET('BSIZE',1024,ERR) DO I=1,NSECT CALL HBNT(7000+I,'EMC'//CHAR(48+I),' ') CALL HBNAME(7000+I,'IEVT',IDEVT,'IDEVT:U') ! event number in 1st pos CALL HBNAME(7000+I,'EMC'//CHAR(48+I),NTRKS, & 'NHIT[0,200]:U,'// ! nb. of tracks in this EMC sector & 'EMCTRK(NHIT)[0,4000]:U,'// ! track number & 'EMCDET(NHIT)[0,2000]:U,'// ! detector number & 'EMCE(NHIT),'// ! deposited energy in MeV & 'EMCX(NHIT),'// ! x along crystal in cm & 'EMCY(NHIT),'// ! y along crystal in cm & 'EMCZ(NHIT),'// ! z along crystal in cm & 'EMCTOF(NHIT)') ! tof in ns ENDDO RETURN ENTRY EMCBOOKLAST ! save ntuples to file DO I=1,NSECT c CALL HPRINT(7000+I) CALL HROUT(7000+I,ICYCLE,' ') END DO RETURN END