Warning, /coresoftware/generators/hijing/doc/readme.txt is written in an unsupported language. File is not indexed.
0001
0002
0003 Version 2.03
0004
0005 ESK parameterization of nuclear shadowing of parton distributions
0006 is replaced by a user implemented subroutine SHADOW. Users can replace
0007 this subroutine with their own parametrization
0008
0009 Versio 2.02
0010
0011 in order to fit the total cross section and dn/dy in pp or ppbar
0012 collisions the pt cut-off now varies with energy roughly as
0013 p0=1.3+0.09*log(sqrt(s)/10)+0.065*log(sqrt(s)/10)**2
0014
0015
0016 Version 2.01
0017
0018 Major modification started on Feb. 22, 1999:
0019 Include options for different parametrization of parton distributions.
0020 IHPR2(7):(D=4)
0021 =1 or 2 Duke-Owns set 1 or 2
0022 =3 MRS D-' parametrization
0023 =4 GRV HO
0024 =5 CTEQ 2pM
0025
0026
0027 IHPR2(6):(D=1)
0028 =1 the original HIJING parametrization of shadowing
0029 =2 Eskola-Kolhinen-Salgado parametrization
0030
0031 *****************
0032 * *
0033 * Version 2.00 *
0034 * *
0035 *****************
0036
0037 Modification Oct. 8, 1998. In hijing, log(ran(nseed)) occasionally
0038 causes overfloat. It is modified to log(max(ran(nseed),1.0e-20)).
0039
0040 replaces all the generic random number generator ran(nseed) by the
0041 portable random number generator in jetset RLU(0).
0042
0043 add the following common block to store the formation time and position
0044 of each hadron
0045 COMMON/HIMAIN3/VATT(130000,4)
0046 For directly produced or massless particles (gamma, dilepton, etc),
0047 we set VATT(I,4)=-10, and VATT(I,1-3)=0
0048 HIPR1(50)=particle formation time in rest-frame with default value
0049 1.0 fm
0050
0051
0052
0053 Version 1.36
0054
0055
0056 (5/30/98) Some statements in HIJSFT ae modified to avoid questionable
0057 branch into another IF block (new compiler on alpha machines does not
0058 like this).
0059
0060 Nothing important has been changed here. A few 'garbage' has been
0061 cleaned up here, like common block HIJJET3 for the sea quark strings
0062 which were originally created to implement the DPM scheme which
0063 later was abadoned in the final version. The lines which operate
0064 on these data are also deleted in the program.
0065
0066
0067
0068 *************************
0069 * *
0070 * Version 1.35 *
0071 * *
0072 *************************
0073
0074 Last modification on Feb. 24, 1997.
0075 There are some changes in the program: subroutine HARDJET is now
0076 consolidated with HIJHRD. HARDJET is used to re-initiate PYTHIA
0077 for the triggered hard processes. Now that is done altogether
0078 with other normal hard processes in modified JETINI. In the new
0079 version one calls JETINI every time one calls HIJHRD. In the new
0080 version the effect of the isospin of the nucleon on hard processes,
0081 especially direct photons is correctly considered.
0082 For A+A collisions, one has to initilize pythia
0083 separately for each type of collisions, pp, pn,np and nn,
0084 or hp and hn for hA collisions. In JETINI we use the following
0085 catalogue for different types of collisions:
0086 h+h: h+h (I_TYPE=1)
0087 h+A: h+p (I_TYPE=1), h+n (I_TYPE=2)
0088 A+h: p+h (I_TYPE=1), n+h (I_TYPE=2)
0089 A+A: p+p (I_TYPE=1), p+n (I_TYPE=2), n+p (I_TYPE=3), n+n (I_TYPE=4)
0090
0091
0092 ****************************
0093 * *
0094 * Version 1.34 *
0095 * *
0096 ****************************
0097
0098 Last modification on January 5, 1998. Two misstakes are corrected in
0099 function G. A Misstake in the subroutine Parton is also corrected.
0100 (These are pointed out by Yasushi Nara).
0101
0102
0103 Last modifcation on April 10, 1996. To conduct final
0104 state radiation, PYTHIA reorganize the two scattered
0105 partons and their final momenta will be a little
0106 different. The summed total momenta of the partons
0107 from the final state radiation are stored in HINT1(26-29)
0108 and HINT1(36-39) which are little different from
0109 HINT1(21-24) and HINT1(41-44).
0110
0111
0112
0113
0114
0115 *************************
0116 * *
0117 * Version 1.33 *
0118 * *
0119 *************************
0120
0121 Last modfication on September 11, 1995. When HIJING and
0122 PYTHIA are initialized, the shadowing is evaluated at
0123 b=0 which is the maximum. This will cause overestimate
0124 of shadowing for peripheral interactions. To correct this
0125 problem, shadowing is set to zero when initializing. Then
0126 use these maximum cross section without shadowing as a
0127 normalization of the Monte Carlo. This however increase
0128 the computing time. IHNT2(16) is used to indicate whether
0129 the sturcture function is called for (IHNT2(16)=1) initialization
0130 or for (IHNT2(16)=0)normal collisions simulation
0131
0132 Last modification on Aagust 28, 1994. Two bugs associate
0133 with the impact parameter dependence of the shadowing is
0134 corrected.
0135
0136 *********************************
0137 * HIJING version 1.31 *
0138 *********************************
0139
0140 No major change has been made in this version. A bug
0141 associated with triggering direct photon is corrected. This
0142 problem was reported by Jim Carroll and Mike Beddo. Another
0143 bug in tracing decay trees reported by Matt Bloomer was also
0144 corrected.
0145
0146 Another feature of triggered jet production (including
0147 direct photon) in HIJING needs to be clarified here. The energy
0148 and momentum stored in HINT1(21)-HINT1(25); HINT1(31)-HINT1(35);
0149 HINT1(41)-HINT1(45); and HINT1(51)-HINT1(55) are the momenta
0150 of the hard scattered partons in the lowest order calculation.
0151 To include the higher order correction, final state radiations
0152 are implemented. In this implementation, momentum and energy
0153 are reshuffled so that the momentum and energy you find in the
0154 finally produced parton or particle might not match those in
0155 HINT1(21)-HINT1(55). But the momentum is still balanced for
0156 the back-to-back jets (or jet and photon). So information in
0157 HINT1(21)-HINT1(55) still give the right direction of the two
0158 high pt jets.
0159
0160 *********************************
0161 * HIJING Version 1.30 *
0162 *********************************
0163
0164 The option to trig on heavy quark production (charm IHPR2(18)=0
0165 or beauty IHPR2(18)=1) is added. To do this, set IHPR2(3)=3. For inclusive
0166 production, one should reset HIPR1(10)=0.0. One can also trig larger pt
0167 QQbar production by giving HIPR1(10) a nonvanishing value. The mass
0168 of the heavy quark in the calculation of the cross section (HINT1(59)-
0169 -HINT1(65)) is given by HIPR1(7) (the default is the charm mass D=1.5).
0170 We also include a separate K-factor for heavy quark and direct photon
0171 production by HIPR1(23)(D=2.0).
0172
0173 *********************************
0174 * HIJING version 1.20 *
0175 *********************************
0176
0177 Beside some modifications to correct some minor errors, there three
0178 major changes since the first version.
0179
0180 (1) An option is added to keep the information of all
0181 particles including the decay chains. The switch of this option
0182 is IHPR(21)=1 (default=0). The array KATT(130000,2) in
0183 COMMON/HIMAIN2/ is expanded to KATT(130000,4). K(I,3) is the line
0184 number of the parent particle of the current line which is produced
0185 via a decay. KATT(I,4) is the status number of the particle: 11=particle
0186 which has decayed; 1=finally or directly produced particle.
0187
0188 (2)Heavy flavor production processes through q+\bar{q}-->c+\bar{c}
0189 and g+g-->c+\bar{c} have been switched on. Charm production is the default,
0190 B-quark option is IHPR2(18)=1 (default=0). Note that when B-quark production
0191 is switched on, charm quark is automatically off.
0192
0193 (3) A subroutine PROFILE(XB) is added to assist the simulation
0194 of minimum biased events of proton-nucleus and nucleus-nucleus collisions.
0195 Refer to the examples in hijing.tex for the simulation.
0196
0197
0198