*********************************************************************** * * subroutine kinebook * * * * * book ntuple for events kinematics * * * * last modified: 02/03/2000 by R. Holzmann * *********************************************************************** implicit none #include "kinetups.inc" #include "user.inc" integer*4 icycle, ih, ierr c write(6,'(/'' KINEBOOK is called'')') call hcdir('//PAWC',' ') call hmdir('KINE','S') ! make a directory KINE call hcdir('//PAWC/KINE',' ') ! change to it call hcdir('//KINETUP',' ') ! change to rzt file call hbset('BSIZE',1024,ierr) call hbnt(idkinetup,'HGEANT kinematics N-tuple',' ') ! book CWN n-tuple call hbname(idkinetup,'EVENT',eventid,'IDEVT:U,IDRUN,'// + 'EBEAM,BPAR,PHI') call hbname(idkinetup,'TRACKS',ntrk,'NTRK[0,4000],'// + 'TRK(NTRK)[0,4000]:U,'// + 'PARTRK(NTRK)[0,4000]:U,'// + 'TVERT(3,NTRK),MECHMED(NTRK):U,'// + 'PID(NTRK)[0,127]:U,PMOM(3,NTRK),'// + 'GENE(NTRK),WGHT(NTRK)') call hbook1(idkinehisto,'KINE multiplicity',5000,-0.5,4999.5,0.0) call hbook1(idkinehisto+1,'KINE particle ID',100,-0.5,99.5,0.0) call hbook1(idkinehisto+2,'KINE mechanism ID',100,-0.5,99.5,0.0) call hbook2(idkinehisto+3,'KINE xy vertex', + 50,-25.,25.,50,-25.,25.,0.0) call hbook1(idkinehisto+4,'KINE z vertex',1000,-100.,100.,0.0) call hcdir('//PAWC',' ') return ********************************************************************** entry kinebooklast ********************************************************************** c write(6,'(/''KINEBOOKLAST is called''/)') call hcdir('//PAWC/KINE',' ') call hldir(' ',' ') call hcdir('//KINETUP',' ' ) call hrout(idkinetup,icycle,' ') ! save n-tuple to file do ih=0,nhisto-1 call hrout(idkinehisto+ih,icycle,' ') ! save histograms enddo call hrend('KINETUP') call hcdir('//PAWC',' ') return end