Back to home page

sPhenix code displayed by LXR

 
 

    


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