]>
Commit | Line | Data |
---|---|---|
0795afa3 | 1 | #include "isajet/pilot.h" |
2 | FUNCTION STRUCW(XW,IW,IH) | |
3 | C | |
4 | C LONGITUDINAL W STRUCTURE FUNCTIONS CALCULATED BY CONVOLUTION | |
5 | C OF EHLQ STRUCTURE FUNCTIONS AT Q**2=AMW**2 WITH F(W/Q) FROM | |
6 | C DAWSON, N.P. B249, 42 (1985). | |
7 | C IW = 1 2 3 4 | |
8 | C GM W+ W- Z0 | |
9 | C WARNING: DEFAULT VALUES ONLY FOR LAMBDA, SIN2W, ETC. | |
10 | C | |
11 | #if defined(CERNLIB_DOUBLE) | |
12 | DOUBLE PRECISION X,XLOG | |
13 | #endif | |
14 | C FOLLOWING CONSTANTS ARE (CV**2+CA**2)/(4*PI**2) | |
15 | DATA CVAW/2.701E-3/,CVAZUU/1.017E-3/,CVAZDD/1.298E-3/ | |
16 | C | |
17 | C STATEMENT FUNCTIONS CALCULATED BY INTEGRATING EACH TERM IN | |
18 | C EHLQ PARAMETERIZATION USING SMP 1.5.0. | |
19 | C | |
20 | C FROM STRUCW2.EX | |
21 | FUVAL(X) = 5.769575427 - 10.13681547*X + 3.042561145*XLOG - | |
22 | $ 0.2798411214*(1./X) - 0.3404284678*(XLOG/X) + 8.296794608*(X**2) | |
23 | $ - 6.017283047*(X**3) + 3.548706099*(X**4) - 1.560232679*(X**5) | |
24 | $ + 0.4342326806*(X**6) - 0.05513649922*(X**7) | |
25 | C | |
26 | FDVAL(X) = 2.533753356 - 4.57001915*X + 0.9589133982*XLOG - | |
27 | $ 0.307417692*(1./X) - 0.1793690733*(XLOG/X) + 4.566485508*(X**2) | |
28 | $ - 3.937617129*(X**3) + 2.773480477*(X**4) - 1.540248513*(X**5) | |
29 | $ + 0.6183097947*(X**6) - 0.1542126064*(X**7) + | |
30 | $ 0.0174859553*(X**8) | |
31 | C | |
32 | FUSEA(X) = -5.503869566 + 15.20985662*X - 5.166205929*XLOG - | |
33 | $ 1.728208206*(1./X) - 0.4215945253*(XLOG/X) - 18.82164974*(X**2) | |
34 | $ + 24.27679709*(X**3) - 27.22445715*(X**4) + 24.98752188*(X**5) | |
35 | $ - 18.12394891*(X**6) + 10.06144542*(X**7) - 4.106763252*(X**8) | |
36 | $ + 1.15847115*(X**9) - 0.2014529514*(X**10) + | |
37 | $ 0.01625760301*(X**11) | |
38 | C | |
39 | FDSEA(X) = -5.478593775 + 15.1097528*X - 5.581750835*XLOG - | |
40 | $ 1.987547927*(1./X) - 0.4944113864*(XLOG/X) - 17.23528157*(X**2) | |
41 | $ + 20.59071315*(X**3) - 21.48775889*(X**4) + 18.40750425*(X**5) | |
42 | $ - 12.46901771*(X**6) + 6.45303238*(X**7) - 2.445316084*(X**8) + | |
43 | $ 0.6363027326*(X**9) - 0.1011660335*(X**10) + | |
44 | $ 0.007376682094*(X**11) | |
45 | C | |
46 | GUVAL1(X) = 6.650062246 + 1.480836233*X - 3.536793901*XLOG - | |
47 | $ 5.08758928*(1./X) - 2.458893299*(XLOG/X) - | |
48 | $ 0.4159265541*(XLOG**2/X) - 0.03949356966*(XLOG**3/X) - | |
49 | $ 0.002175877338*(XLOG**4/X) - 6.382577207E-5*(XLOG**5/X) - | |
50 | $ 7.324244818E-7*(XLOG**6/X) + 0.575121435*(X*XLOG) + | |
51 | $ 0.09704190061*(X*XLOG**2) + 791./90953.*(X *XLOG**3) | |
52 | GUVAL2(X) = | |
53 | $ 3.962955366E-4*(X*XLOG**4) + 6.591820335E-6*(X *XLOG**5) - | |
54 | $ 0.07706665525*(X**2*XLOG) - 0.01260088982*(X**2*XLOG**2) - | |
55 | $ 0.001080476729*(X**2*XLOG**3) - 4.647425231E-5*(X**2*XLOG**4) - | |
56 | $ 7.324244829E-7*(X**2*XLOG**5) - 0.2027508558*(X**2) - | |
57 | $ 0.6721633358*(XLOG**2) - 0.07339139987*(XLOG**3) - | |
58 | $ 0.004744695259*(XLOG**4) - 1.651100349E-4*XLOG**5 - | |
59 | $ 2.197273446E-6*XLOG**6 | |
60 | C | |
61 | GDVAL1(X) = 2.403815112 + 1.359504335*X - 2.144999226*XLOG - | |
62 | $ 2.357471591*(1./X) - 1.13339787*(XLOG/X) - | |
63 | $ 0.2003574358*(XLOG**2/X) - 0.02221455952*(XLOG**3/X) - | |
64 | $ 0.001601991722*(XLOG**4/X) - 6.878677010E-5*(XLOG**5/X) - | |
65 | $ 1.323789434E-6*(XLOG**6/X) + 0.4984901703*(X*XLOG) + | |
66 | $ 0.09080940556*(X*XLOG**2) + 0.01045471879*(X *XLOG**3) | |
67 | GDVAL2(X) = | |
68 | $ 7.339489291E-4*(X*XLOG**4) + 2.382820980E-5*(X *XLOG**5) - | |
69 | $ 0.1348586553*(X**2*XLOG) - 0.02441501092*(X**2*XLOG**2) - | |
70 | $ 0.002734669503*(X**2*XLOG**3) - 1.807502877E-4*(X**2*XLOG**4) - | |
71 | $ 5.295157736E-6*(X**2*XLOG**5) + 0.01826324487*(X**3*XLOG) + | |
72 | $ 0.003291926818*(X**3*XLOG**2) + 3.636001057E-4*(X**3*XLOG**3) + | |
73 | $ 2.342115438E-5*(X**3*XLOG**4) + 6.618947159E-7*(X**3*XLOG**5) - | |
74 | $ 0.3679347826*(X**2) + 0.04985470259*(X**3) - | |
75 | $ 0.3928663839*(XLOG**2) - 0.04732954832*(XLOG**3) - | |
76 | $ 0.003974205548*( XLOG**4) - 2.116051876E-4*XLOG**5 - | |
77 | $ 5.295157736E-6*XLOG**6 | |
78 | C | |
79 | GUSEA1(X) = -0.8251281831 + 1.555766474*X - 0.476618796*XLOG - | |
80 | $ 0.157877015*(1./X) + 0.3273497735*(XLOG/X) + | |
81 | $ 0.1184829659*(XLOG**2/X) + 0.01147973292*(XLOG**3/X) + | |
82 | $ 0.001370332595*(XLOG**4/X) + 4.084139287E-5*(XLOG**5/X) + | |
83 | $ 2.284079310E-6*(XLOG**6/X) - 1.681676555*(X*XLOG) + | |
84 | $ 0.01771802464*(X*XLOG**2) - 0.04546554244*(X *XLOG**3) - | |
85 | $ 3.454606694E-4*(X*XLOG**4) - 1.438969965E-4*(X *XLOG**5) + | |
86 | $ 0.8875664376*(X**2*XLOG) + 0.04051742981*(X**2 *XLOG**2) | |
87 | GUSEA2(X) = | |
88 | $ 0.0254151271*(X**2*XLOG**3) + 4.583985126E-4*(X**2*XLOG**4) + | |
89 | $ 7.994277584E-5*(X**2*XLOG**5) - 0.4492144518*(X**3*XLOG) - | |
90 | $ 0.02965496152*(X**3*XLOG**2) - 0.01288679853*(X**3*XLOG**3) - | |
91 | $ 2.791634921E-4*(X**3*XLOG**4) - 3.997138792E-5*(X**3*XLOG**5) + | |
92 | $ 0.1638328221*(X**4*XLOG) + 0.01250393016*(X**4*XLOG**2) + | |
93 | $ 0.004685172364*(X**4*XLOG**3) + 1.100919901E-4*(X**4*XLOG**4) + | |
94 | $ 1.438969967E-5*(X**4*XLOG**5) - 0.0367848506*(X**5*XLOG) - | |
95 | $ 0.003035305139*(X**5*XLOG**2) - 0.001048416117*(X**5*XLOG**3) - | |
96 | $ 2.579726667E-5*(X**5*XLOG**4) - 3.197711036E-6*(X**5*XLOG**5) + | |
97 | $ 0.003783780648*(X**6*XLOG) + 3.278296157E-4*(X**6*XLOG**2) + | |
98 | $ 1.075296502E-4*(X**6*XLOG**3) + 2.725601849E-6*(X**6*XLOG**4) | |
99 | GUSEA3(X) = | |
100 | $ 3.262970444E-7*(X**6*XLOG**5) - 0.2259054436*(X**2) + | |
101 | $ 0.003364712414*(X**3) + 0.01992787001*(X**4) - 0.007399430903* | |
102 | $ (X**5) + 9.652150086E-4*X**6 + 0.554652844*(XLOG**2) + | |
103 | $ 0.011217842*(XLOG**3) + 0.007692743973*(XLOG**4) + | |
104 | $ 9.402708800E-5*XLOG**5 + 1.598855517E-5*XLOG**6 | |
105 | C | |
106 | GDSEA1(X) = -0.9201807217 + 2.243479849*X - 0.899698589*XLOG - | |
107 | $ 0.3970657521*(1./X) + 0.2818290666*(XLOG/X) + | |
108 | $ 0.120664241*(XLOG**2/X) + 0.01043451714 *(XLOG**3/X) + | |
109 | $ 0.001191246128*(XLOG**4/X) + 3.001102810E-5*(XLOG**5/X) + | |
110 | $ 2.039937816E-6*(XLOG**6/X) - 1.787863927*(X*XLOG) + | |
111 | $ 826./63247.*(X*XLOG**2) - 0.04263009293*(X*XLOG**3) + | |
112 | $ 3.087205491E-5*( X*XLOG**4) - 1.285160824E-4*(X*XLOG**5) + | |
113 | $ 0.9517638436*(X**2*XLOG) + 0.04047432532*(X**2*XLOG**2) + | |
114 | $ 0.02332032497*(X**2*XLOG**3) + 2.208416029E-4*(X**2*XLOG**4) | |
115 | GDSEA2(X) = | |
116 | $ 7.139782355E-5*(X**2*XLOG**5) - 0.4816526149*(X**3*XLOG) - | |
117 | $ 0.02893293689*(X**3*XLOG**2) - 0.01172595961*(X**3*XLOG**3) - | |
118 | $ 1.550444419E-4*(X**3*XLOG**4) - 3.569891178E-5*(X**3*XLOG**5) + | |
119 | $ 0.1755025883*(X**4*XLOG) + 0.01209064681*(X**4*XLOG**2) + | |
120 | $ 0.004244259809*(X**4*XLOG**3) + 6.438373778E-5*(X**4*XLOG**4) + | |
121 | $ 1.285160826E-5*(X**4*XLOG**5) - 0.03937044174*(X**5*XLOG) - | |
122 | $ 0.002920614679*(X**5*XLOG**2) - 9.471446970E-4*(X**5*XLOG**3) - | |
123 | $ 1.549746135E-5*(X**5*XLOG**4) - 2.855912944E-6*(X**5*XLOG**5) | |
124 | GDSEA3(X) = | |
125 | $ 0.004046756409*(X**6*XLOG) + 3.144555106E-4*(X**6*XLOG**2) + | |
126 | $ 9.696129739E-5*(X**6*XLOG**3) + 1.664636351E-6*(X**6*XLOG**4) + | |
127 | $ 2.914196881E-7*(X**6*XLOG**5) - 0.5703745807*(X**2) + | |
128 | $ 0.1676594704*(X**3) - 0.03765961644*(X**4) + | |
129 | $ 0.005180399826*(X**5) - 3.032419889E-4*X**6 + | |
130 | $ 0.5916013402*(XLOG**2) + 0.0113078292*(XLOG**3) + | |
131 | $ 0.007094724813*( XLOG**4) + 3.872242009E-5*XLOG**5 + | |
132 | $ 1.427956471E-5*XLOG**6 | |
133 | C | |
134 | EUVAL(X) = -3.398748694 + 0.6266420937*(1./X) | |
135 | EDVAL(X) = -1.300464877 + 0.2267175031*(1./X) | |
136 | EUSEA(X) = -0.4281951222 + 0.0600001177*(1./X) | |
137 | EDSEA(X) = -0.5744690066 + 0.08143317382*(1./X) | |
138 | C | |
139 | C ENTRY | |
140 | C | |
141 | IF(XW.LE.0..OR.XW.GE.1.) THEN | |
142 | STRUCW=0. | |
143 | RETURN | |
144 | ENDIF | |
145 | X=XW | |
146 | XLOG=LOG(X) | |
147 | C | |
148 | IIW=IW | |
149 | IF(IH.EQ.-1120) THEN | |
150 | IF(IW.EQ.2) IIW=3 | |
151 | IF(IW.EQ.3) IIW=2 | |
152 | ELSEIF(IH.EQ.1220) THEN | |
153 | IF(IW.EQ.2) IIW=3 | |
154 | IF(IW.EQ.3) IIW=2 | |
155 | ENDIF | |
156 | C | |
157 | IF(XW.GT..1) GO TO 1000 | |
158 | C | |
159 | C STRUCW = XW*F(XW) FOR IIW=W+ IN PROTON, XW<.1 | |
160 | IF(IIW.EQ.2) THEN | |
161 | G1=GUVAL1(X) | |
162 | G1=G1+GUVAL2(X) | |
163 | G2=GUSEA1(X) | |
164 | G2=G2+GUSEA2(X) | |
165 | G2=G2+GUSEA3(X) | |
166 | G3=GDSEA1(X) | |
167 | G3=G3+GDSEA2(X) | |
168 | G3=G3+GDSEA3(X) | |
169 | SUM=G1+G2+G3+EUVAL(X)+EUSEA(X)+EDSEA(X) | |
170 | STRUCW=X*CVAW*SUM | |
171 | C STRUCW = XW*F(XW) FOR IIW=W- IN PROTON, XW<.1 | |
172 | ELSEIF(IIW.EQ.3) THEN | |
173 | G1=GDVAL1(X) | |
174 | G1=G1+GDVAL2(X) | |
175 | G2=GDSEA1(X) | |
176 | G2=G2+GDSEA2(X) | |
177 | G2=G2+GDSEA3(X) | |
178 | G3=GUSEA1(X) | |
179 | G3=G3+GUSEA2(X) | |
180 | G3=G3+GUSEA3(X) | |
181 | SUM=G1+G2+G3+EDVAL(X)+EDSEA(X)+EUSEA(X) | |
182 | STRUCW=X*CVAW*SUM | |
183 | C STRUCW=XW*F(XW) FOR IIW=Z0 IN PROTON, XW<.1 | |
184 | ELSEIF(IIW.EQ.4) THEN | |
185 | G1=GUVAL1(X) | |
186 | G1=G1+GUVAL2(X) | |
187 | G2=GUSEA1(X) | |
188 | G2=G2+GUSEA2(X) | |
189 | G2=G2+GUSEA3(X) | |
190 | SUMU=G1+2.*G2+EUVAL(X)+2.*EUSEA(X) | |
191 | G1=GDVAL1(X) | |
192 | G1=G1+GDVAL2(X) | |
193 | G2=GDSEA1(X) | |
194 | G2=G2+GDSEA2(X) | |
195 | G2=G2+GDSEA3(X) | |
196 | SUMD=G1+2.*G2+EDVAL(X)+2.*EDSEA(X) | |
197 | STRUCW=X*(CVAZUU*SUMU+CVAZDD*SUMD) | |
198 | ENDIF | |
199 | IF(STRUCW.LT.0.) STRUCW=0. | |
200 | RETURN | |
201 | C | |
202 | 1000 CONTINUE | |
203 | C | |
204 | C STRUCW=XW*F(XW) FOR IIW=W+ IN PROTON, XW>.1 | |
205 | IF(IIW.EQ.2) THEN | |
206 | F1=FUVAL(X) | |
207 | F2=FUSEA(X) | |
208 | F3=FDSEA(X) | |
209 | SUM=F1+F2+F3 | |
210 | STRUCW=X*CVAW*SUM | |
211 | C STRUCW=XW*F(XW) FOR IIW=W- IN PROTON, XW>.1 | |
212 | ELSEIF(IIW.EQ.3) THEN | |
213 | F1=FDVAL(X) | |
214 | F2=FDSEA(X) | |
215 | F3=FUSEA(X) | |
216 | SUM=F1+F2+F3 | |
217 | STRUCW=X*CVAW*SUM | |
218 | C STRUCW=XW*F(XW) FOR IIW=Z0 IN PROTON, XW>.1 | |
219 | ELSEIF(IIW.EQ.4) THEN | |
220 | F1=FUVAL(X) | |
221 | F2=FUSEA(X) | |
222 | SUMU=F1+2.*F2 | |
223 | F1=FDVAL(X) | |
224 | F2=FDSEA(X) | |
225 | SUMD=F1+2.*F2 | |
226 | STRUCW=X*(CVAZUU*SUMU+CVAZDD*SUMD) | |
227 | ENDIF | |
228 | IF(STRUCW.LT.0.) STRUCW=0. | |
229 | RETURN | |
230 | END |