hijing/src/hidata.f

0001 C
0002 C
0003 C***************************************************************
0004
0005         BLOCK DATA HIDATA
0006         REAL*8 XL(10),XU(10),ACC
0007         COMMON/BVEG1/XL,XU,ACC,NDIM,NCALL,ITMX,NPRN
0008         SAVE  /BVEG1/
0009         COMMON/SEEDVAX/NUM1
0010         SAVE  /SEEDVAX/
0011         COMMON/HIPARNT/HIPR1(100),IHPR2(50),HINT1(100),IHNT2(50)
0012         SAVE  /HIPARNT/
0013         COMMON/RANSEED/NSEED
0014         SAVE  /RANSEED/
0015         COMMON/HIMAIN1/ NATT,EATT,JATT,NT,NP,N0,N01,N10,N11,IERRSTAT
0016         SAVE  /HIMAIN1/
0017 C+++BAC
0018 C       COMMON/HIMAIN2/KATT(130000,4),PATT(130000,4)
0019 C       SAVE  /HIMAIN2/
0020 C
0021         COMMON/HIMAIN2/KATT(130000,4),PATT(130000,4), VATT(130000,4)
0022         SAVE  /HIMAIN2/
0023 C---BAC
0024
0025         COMMON/HISTRNG/NFP(300,15),PP(300,15),NFT(300,15),PT(300,15)
0026         SAVE  /HISTRNG/
0027         COMMON/HIJCRDN/YP(3,300),YT(3,300)
0028         SAVE  /HIJCRDN/
0029         COMMON/HIJJET1/NPJ(300),KFPJ(300,500),PJPX(300,500),
0030      &               PJPY(300,500),PJPZ(300,500),PJPE(300,500),
0031      &               PJPM(300,500),NTJ(300),KFTJ(300,500),
0032      &               PJTX(300,500),PJTY(300,500),PJTZ(300,500),
0033      &               PJTE(300,500),PJTM(300,500)
0034         SAVE  /HIJJET1/
0035         COMMON/HIJJET2/NSG,NJSG(900),IASG(900,3),K1SG(900,100),
0036      &          K2SG(900,100),PXSG(900,100),PYSG(900,100),
0037      &          PZSG(900,100),PESG(900,100),PMSG(900,100)
0038         SAVE  /HIJJET2/
0039         COMMON/HIJDAT/HIDAT0(10,10),HIDAT(10)
0040         SAVE  /HIJDAT/
0041         COMMON/HIPYINT/MINT4,MINT5,ATCO(200,20),ATXS(0:200)
0042         SAVE  /HIPYINT/
0043         DATA NUM1/30123984/,XL/10*0.D0/,XU/10*1.D0/
0044         DATA NCALL/1000/ITMX/100/ACC/0.01/NPRN/0/
0045 C...give all the switchs and parameters the default values
0046 c       DATA NSEED/74769375/
0047 c seed changed to zero, otherwise ran(seed) in gfortran will not generate random numbers
0048         DATA NSEED/74769375/
0049         DATA HIPR1/
0050      &  1.5,  0.35, 0.5,  0.9,  2.0,  0.1,  1.5,  2.0, -1.0, -2.25,
0051      &  2.0,  0.5,  1.0,  2.0,  0.2,  2.0,  2.5,  0.3,  0.1,  1.4,
0052      &  1.6,  1.0,  2.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.4,  57.0,
0053      &  28.5, 3.9,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
0054      &  3.141592654,
0055      &  0.0,  0.4,  0.1,  1.5,  0.1, 0.25, 0.0,  0.5,  0.0,  0.0,
0056      &  50*0.0/
0057
0058         DATA IHPR2/
0059      &  1,    3,    0,    1,    1,    1,    1,    10,    0,    0,
0060      &  1,    1,    1,    1,    0,    0,    1,     0,    0,    1,
0061      &  30*0/
0062
0063         DATA HINT1/100*0/
0064         DATA IHNT2/50*0/
0065
0066 C...initialize all the data common blocks
0067         DATA NATT/0/EATT/0.0/JATT/0/NT/0/NP/0/N0/0/N01/0/N10/0/N11/0/
0068 C+++BAC
0069 C
0070 C       DATA KATT/520000*0/PATT/520000*0.0/
0071 C
0072         DATA KATT/520000*0/PATT/520000*0.0/VATT/520000*0.0/
0073 C---BAC
0074
0075         DATA NFP/4500*0/PP/4500*0.0/NFT/4500*0/PT/4500*0.0/
0076
0077         DATA YP/900*0.0/YT/900*0.0/
0078
0079         DATA NPJ/300*0/KFPJ/150000*0/PJPX/150000*0.0/PJPY/150000*0.0/
0080      &  PJPZ/150000*0.0/PJPE/150000*0.0/PJPM/150000*0.0/
0081         DATA NTJ/300*0/KFTJ/150000*0/PJTX/150000*0.0/PJTY/150000*0.0/
0082      &  PJTZ/150000*0.0/PJTE/150000*0.0/PJTM/150000*0.0/
0083
0084         DATA NSG/0/NJSG/900*0/IASG/2700*0/K1SG/90000*0/K2SG/90000*0/
0085      &  PXSG/90000*0.0/PYSG/90000*0.0/PZSG/90000*0.0/PESG/90000*0.0/
0086      &  PMSG/90000*0.0/
0087         DATA MINT4/0/MINT5/0/ATCO/4000*0.0/ATXS/201*0.0/
0088         DATA (HIDAT0(1,I),I=1,10)/0.0,0.0,0.0,0.0,0.0,0.0,2.25,
0089      &          2.5,4.0,4.1/
0090         DATA (HIDAT0(2,I),I=1,10)/2.0,3.0,5.0,6.0,7.0,8.0,8.0,10.0,
0091      &          10.0,10.0/
0092         DATA (HIDAT0(3,I),I=1,10)/1.0,0.8,0.8,0.7,0.45,0.215,
0093      &          0.21,0.19,0.19,0.19/
0094         DATA (HIDAT0(4,I),I=1,10)/0.35,0.35,0.3,0.3,0.3,0.3,
0095      &          0.5,0.6,0.6,0.6/
0096         DATA (HIDAT0(5,I),I=1,10)/23.8,24.0,26.0,26.2,27.0,28.5,28.5,
0097      &          28.5,28.5,28.5/
0098         DATA ((HIDAT0(J,I),I=1,10),J=6,9)/40*0.0/
0099         DATA (HIDAT0(10,I),I=1,10)/5.0,20.0,53.0,62.0,100.0,200.0,
0100      &          546.0,900.0,1800.0,4000.0/
0101         DATA HIDAT/10*0.0/
0102         END