Warning, /geant4/examples/extended/hadronic/Hadr02/hijing1_383/README is written in an unsupported language. File is not indexed.
0001 HIJING interface
0002 ----------------
0003
0004 Khaled Abdel-Waged
0005 --------------------
0006
0007 version 0.0, 1-11-2012
0008 ----------------------
0009
0010 >Installation requirements
0011 --------------------------
0012
0013 The Geant4 interface to HIJING uses the following software tools and packages:
0014
0015 1. HIJING (available at fttp://nta0.lbl.gov/pub/xnwang/hijing)
0016 2. The following Geant4 versions:
0017 − ver.9.5 or above
0018
0019 The interface and original HIJING code have been compiled and tested using:
0020 • gcc 4.1.2 with gfortran (FORTRAN95) and GNUmake)
0021
0022 Operating systems used for this test is:
0023 Red Hat Linux 4.1.2-64
0024
0025
0026 >Changes in the Fortran code
0027 ----------------------------
0028
0029 1) The main subroutines
0030 SUBROUTINE HIJSET(EFRM, FRAME, PROJ, TARG, IAP, IZP, IAT, IZT)
0031 is replaced with
0032 SUBROUTINE HIJSET(EFRM)
0033
0034 SUBROUTINE HIJING (FRAME, BMIN, BMAX)
0035 is replaced with
0036 SUBROUTINE HIJING (BMIN, BMAX)
0037
0038
0039 The input variables are directly inserted into the following common block:
0040
0041 COMMON/HIPARNT/HIPR1(100),IHPR2(50),HINT1(100),IHNT2(50)
0042 Projectile:
0043 IHNT2(1)=IAP // Nucleus mass number
0044 or IHNT2(1)=1 //Hadron
0045 IHNT2(2)=IZP //charge
0046 IHNT2(5)=0 //id code
0047 Target:
0048 IHNT2(3)=IAT
0049 IHNT2(4)=IZT
0050 IHNT2(6)=0 //id Target (Fixed)
0051 Rest Mass:
0052 HINT1(8)= //projectile
0053 HINT1(9)= //Target
0054
0055 2) In subroutines HIJSET and HIJING,
0056
0057 Since Geant4 hadronic cascade models always works in the LABoratory frame
0058
0059 the statement Frame=”LAB” is inserted.
0060
0061
0062 3)Random number generator
0063
0064 The two random generator functions become one!
0065
0066 This is done by replacing
0067 RAN(NSEED) with RLU(0) in hijing1.383.f
0068
0069
0070 >Consequences of conversion from f77 (g77) to gfortran
0071 ------------------------------------------------------
0072
0073 When running HIJING code in gfortran directly, the execution is blocked.
0074 This problem is solved by the following changes
0075
0076 1) You have to add in the GNUmake file, the line
0077
0078 FFLAGS=-fno-automatic
0079
0080 This treats each program unit as if the SAVE statement were specified
0081 for every local variable and array referenced in it.
0082
0083
0084 2) Problem related to the function ROMG(x) in hijing1.383.f:
0085
0086 FUNCTION ROMG(X)
0087 C This gives the eikonal function from a table
0088 C calculated in the first call
0089 DIMENSION FR(0:1000)
0090 DATA I0/0/
0091
0092 COMMON/EIKONAL/FR !New->Khaled
0093
0094 IF(I0.NE.0) GO TO 100
0095 DO 50 I=1,1001
0096 XR=(I-1)*0.01
0097 FR(I-1)=OMG0(XR)
0098 50 CONTINUE
0099 100 I0=1
0100 IF(X.GE.10.0) THEN
0101 ROMG=0.0
0102 RETURN
0103 ENDIF
0104 IX=INT(X*100)
0105 ROMG=(FR(IX)*((IX+1)*0.01-X)+FR(IX+1)*(X-IX*0.01))/0.01
0106 RETURN
0107 END
0108
0109 Our analysis shows that the array FR() is not saved in the subsequent calls of the function ROMG().
0110 Therefore, the statement “COMMON/EIKONAL/FR” is inserted, as shown above.
0111
0112
0113 >Interface design
0114 -----------------
0115
0116 The use of HIJING physics in Geant4 has resulted in the introduction of a new event model (G4HIJING_Model).
0117 The G4HIJING_Model class is derived from G4HadronicInteraction, and is defined within the Geant4 user physics
0118 list if access to HIJING physics is required. It controls initialisation of HIJING through common block variables.
