*********************************************************************** * * SUBROUTINE MULTBOOK * * * MODIFIED BY: R. Holzmann DATE: 10/07/97 * * * * book histograms for multiplcity array of July 1997 test exp. * * * *********************************************************************** IMPLICIT NONE #include "multtups.inc" CHARACTER*63 NTUPLE_FILE, NT_FILE_NAME INTEGER*4 ICYCLE, ISTAT, I, IDTUP CHARACTER*40 NAME C WRITE(6,'(/'' MULTBOOK is called'')') CALL HCDIR('//PAWC',' ') CALL HMDIR('MULT','S') c NTUPLE_FILE='MULT.RZT' c WRITE(6,'(//A,$)') ' NTUPLE_FILE_NAME (=''MULT.RZT''): ' c READ(5,'(A)') NT_FILE_NAME c IF (NT_FILE_NAME.NE.' ') NTUPLE_FILE = NT_FILE_NAME c CALL HROPEN(81,'MULTTUP',NTUPLE_FILE,'N',4096,ISTAT) CALL HCDIR('//PAWC/MULT',' ') DO I=1,NDET ! book energy spectra NAME = 'Multiplicity Array E'//CHAR(48+I/10)// 1 CHAR(48+I-(I/10)*10)//' (MeV)$' CALL HBOOK1(5000+I,NAME,1000,0.0,MAXENERGY,0.0) END DO CALL HBOOK1(5000+NDET+1,'Multiplicity Array Hit Pattern$', 1 NDET+2,-0.5,FLOAT(NDET)+1.5,0.0) CALL HBOOK1(5000+NDET+2,'Multiplictiy Array Hit Multiplcity$', 1 NDET+6,-0.5,FLOAT(NDET)+5.5,0.0) IDTUP = 5000+NDET+3 CALL HCDIR('//MULTTUP',' ') CALL HBNT(IDTUP,'Multiplicity Array N-tuple',' ') ! book CWN n-tuple CALL HBNAME(IDTUP,'MULT',NHIT, 1 'NHIT[0,45],ID(NHIT)[0,45],E(NHIT)') CALL HCDIR('//PAWC',' ') RETURN ************************************************************ ENTRY MULTBOOKLAST ************************************************************ C WRITE(6,'(/''MULTBOOKLAST is called''/)') CALL HCDIR('//PAWC/MULT',' ') CALL HLDIR(' ',' ') CALL HCDIR('//MULTTUP',' ' ) DO I=1,NDET+3 CALL HROUT(5000+I,ICYCLE,' ') ! save histograms and n-tuple to file END DO c CALL HREND('MULTTUP') c CLOSE(81) c CALL HCDIR('//PAWC',' ') RETURN END