[ 1 ]  [ 2 ]  [ 3 ]  [ 4 ]  [ 5 ]  [ 6 ]  [ 7 ]  [ 8 ]  [ 9 ]  [ 10 ]  [ 11 ]  [ 12 ]  [ 13 ]  [ 14 ]  [ 15 ]  [ back ] 

HISTIN

The main attention in the case of the use of HBOOK is that it makes use of REAL*4 variables, while the FLUKA environment is by default in double precision (REAL*8). For our example we intend to book a simple disk-resident ntupla opening the hbk.his file, having 8 variables: the muon charge, the residual energy, the 3 coordinates and the 3 direction cosines at boundary crossing. Notice that we have chosen to assign the logical unit 63 to the HBOOK file. Any number larger than 20 is good, provided that it is not already used by the user himself for other scoring purposes.

C------------------------------------------------------------
      SUBROUTINE HISTIN
C------------------------------------------------------------
C     
C     HBOOK HISTOGRAM BOOKING
C
      PARAMETER (NNTUPLE=8)
      COMMON/PAWC/HMEMOR(200000)
      REAL*4 HMEMOR
      COMMON/QUEST/IQUEST(100)
      CHARACTER*8 CHTAG(NNTUPLE)
      DATA CHTAG /'MUTYP' ,'ERES' ,'XMU'  ,'YMU', 'ZMU', 'THX',
     $     'THY', 'THZ'/ 
      WRITE(*,*) 'ENTER INTO HIST INIT.'
      IQUEST(10)=65000
      CALL HLIMIT(200000)
C
      NLUNI1 = 63
      CALL HROPEN(NLUNI1, 'Muons','hbk.his','N',1024,ISTAT)
      WRITE(*,*) 'HROPEN: ISTAT=',ISTAT
      CALL HBOOKN(100,'Muons',NNTUPLE,'//Muons',NNTUPLE*100,CHTAG)
C
      RETURN
      END

Giuseppe Battistoni; INFN, Milano
Francesco Cerutti; CERN, Geneva

Last updated: 10th of December, 2008