ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
hidata.f
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file hidata.f
1 C
2 C
3 C***************************************************************
4 
5  BLOCK DATA hidata
6  REAL*8 xl(10),xu(10),acc
7  common/bveg1/xl,xu,acc,ndim,ncall,itmx,nprn
8  SAVE /bveg1/
9  common/seedvax/num1
10  SAVE /seedvax/
11  common/hiparnt/hipr1(100),ihpr2(50),hint1(100),ihnt2(50)
12  SAVE /hiparnt/
13  common/ranseed/nseed
14  SAVE /ranseed/
15  common/himain1/ natt,eatt,jatt,nt,np,n0,n01,n10,n11,ierrstat
16  SAVE /himain1/
17 C+++BAC
18 C COMMON/HIMAIN2/KATT(130000,4),PATT(130000,4)
19 C SAVE /HIMAIN2/
20 C
21  common/himain2/katt(130000,4),patt(130000,4), vatt(130000,4)
22  SAVE /himain2/
23 C---BAC
24 
25  common/histrng/nfp(300,15),pp(300,15),nft(300,15),pt(300,15)
26  SAVE /histrng/
27  common/hijcrdn/yp(3,300),yt(3,300)
28  SAVE /hijcrdn/
29  common/hijjet1/npj(300),kfpj(300,500),pjpx(300,500),
30  & pjpy(300,500),pjpz(300,500),pjpe(300,500),
31  & pjpm(300,500),ntj(300),kftj(300,500),
32  & pjtx(300,500),pjty(300,500),pjtz(300,500),
33  & pjte(300,500),pjtm(300,500)
34  SAVE /hijjet1/
35  common/hijjet2/nsg,njsg(900),iasg(900,3),k1sg(900,100),
36  & k2sg(900,100),pxsg(900,100),pysg(900,100),
37  & pzsg(900,100),pesg(900,100),pmsg(900,100)
38  SAVE /hijjet2/
39  common/hijdat/hidat0(10,10),hidat(10)
40  SAVE /hijdat/
41  common/hipyint/mint4,mint5,atco(200,20),atxs(0:200)
42  SAVE /hipyint/
43  DATA num1/30123984/,xl/10*0.d0/,xu/10*1.d0/
44  DATA ncall/1000/itmx/100/acc/0.01/nprn/0/
45 C...give all the switchs and parameters the default values
46 c DATA NSEED/74769375/
47 c seed changed to zero, otherwise ran(seed) in gfortran will not generate random numbers
48  DATA nseed/74769375/
49  DATA hipr1/
50  & 1.5, 0.35, 0.5, 0.9, 2.0, 0.1, 1.5, 2.0, -1.0, -2.25,
51  & 2.0, 0.5, 1.0, 2.0, 0.2, 2.0, 2.5, 0.3, 0.1, 1.4,
52  & 1.6, 1.0, 2.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.4, 57.0,
53  & 28.5, 3.9, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
54  & 3.141592654,
55  & 0.0, 0.4, 0.1, 1.5, 0.1, 0.25, 0.0, 0.5, 0.0, 0.0,
56  & 50*0.0/
57 
58  DATA ihpr2/
59  & 1, 3, 0, 1, 1, 1, 1, 10, 0, 0,
60  & 1, 1, 1, 1, 0, 0, 1, 0, 0, 1,
61  & 30*0/
62 
63  DATA hint1/100*0/
64  DATA ihnt2/50*0/
65 
66 C...initialize all the data common blocks
67  DATA natt/0/eatt/0.0/jatt/0/nt/0/np/0/n0/0/n01/0/n10/0/n11/0/
68 C+++BAC
69 C
70 C DATA KATT/520000*0/PATT/520000*0.0/
71 C
72  DATA katt/520000*0/patt/520000*0.0/vatt/520000*0.0/
73 C---BAC
74 
75  DATA nfp/4500*0/pp/4500*0.0/nft/4500*0/pt/4500*0.0/
76 
77  DATA yp/900*0.0/yt/900*0.0/
78 
79  DATA npj/300*0/kfpj/150000*0/pjpx/150000*0.0/pjpy/150000*0.0/
80  & pjpz/150000*0.0/pjpe/150000*0.0/pjpm/150000*0.0/
81  DATA ntj/300*0/kftj/150000*0/pjtx/150000*0.0/pjty/150000*0.0/
82  & pjtz/150000*0.0/pjte/150000*0.0/pjtm/150000*0.0/
83 
84  DATA nsg/0/njsg/900*0/iasg/2700*0/k1sg/90000*0/k2sg/90000*0/
85  & pxsg/90000*0.0/pysg/90000*0.0/pzsg/90000*0.0/pesg/90000*0.0/
86  & pmsg/90000*0.0/
87  DATA mint4/0/mint5/0/atco/4000*0.0/atxs/201*0.0/
88  DATA (hidat0(1,i),i=1,10)/0.0,0.0,0.0,0.0,0.0,0.0,2.25,
89  & 2.5,4.0,4.1/
90  DATA (hidat0(2,i),i=1,10)/2.0,3.0,5.0,6.0,7.0,8.0,8.0,10.0,
91  & 10.0,10.0/
92  DATA (hidat0(3,i),i=1,10)/1.0,0.8,0.8,0.7,0.45,0.215,
93  & 0.21,0.19,0.19,0.19/
94  DATA (hidat0(4,i),i=1,10)/0.35,0.35,0.3,0.3,0.3,0.3,
95  & 0.5,0.6,0.6,0.6/
96  DATA (hidat0(5,i),i=1,10)/23.8,24.0,26.0,26.2,27.0,28.5,28.5,
97  & 28.5,28.5,28.5/
98  DATA ((hidat0(j,i),i=1,10),j=6,9)/40*0.0/
99  DATA (hidat0(10,i),i=1,10)/5.0,20.0,53.0,62.0,100.0,200.0,
100  & 546.0,900.0,1800.0,4000.0/
101  DATA hidat/10*0.0/
102  END