]>
Commit | Line | Data |
---|---|---|
1 | * $Id$ | |
2 | ||
3 | C********************************************************************* | |
4 | ||
5 | FUNCTION ULALPS_HIJING(Q2) | |
6 | ||
7 | C...Purpose: to give the value of alpha_strong. | |
8 | #include "ludat1_hijing.inc" | |
9 | #include "ludat2_hijing.inc" | |
10 | ||
11 | C...Constant alpha_strong trivial. | |
12 | IF(MSTU(111).LE.0) THEN | |
13 | ULALPS_HIJING=PARU(111) | |
14 | MSTU(118)=MSTU(112) | |
15 | PARU(117)=0. | |
16 | PARU(118)=PARU(111) | |
17 | RETURN | |
18 | ENDIF | |
19 | ||
20 | C...Find effective Q2, number of flavours and Lambda. | |
21 | Q2EFF=Q2 | |
22 | IF(MSTU(115).GE.2) Q2EFF=MAX(Q2,PARU(114)) | |
23 | NF=MSTU(112) | |
24 | ALAM2=PARU(112)**2 | |
25 | 100 IF(NF.GT.MAX(2,MSTU(113))) THEN | |
26 | Q2THR=PARU(113)*PMAS(NF,1)**2 | |
27 | IF(Q2EFF.LT.Q2THR) THEN | |
28 | NF=NF-1 | |
29 | ALAM2=ALAM2*(Q2THR/ALAM2)**(2./(33.-2.*NF)) | |
30 | GOTO 100 | |
31 | ENDIF | |
32 | ENDIF | |
33 | 110 IF(NF.LT.MIN(8,MSTU(114))) THEN | |
34 | Q2THR=PARU(113)*PMAS(NF+1,1)**2 | |
35 | IF(Q2EFF.GT.Q2THR) THEN | |
36 | NF=NF+1 | |
37 | ALAM2=ALAM2*(ALAM2/Q2THR)**(2./(33.-2.*NF)) | |
38 | GOTO 110 | |
39 | ENDIF | |
40 | ENDIF | |
41 | IF(MSTU(115).EQ.1) Q2EFF=Q2EFF+ALAM2 | |
42 | PARU(117)=SQRT(ALAM2) | |
43 | ||
44 | C...Evaluate first or second order alpha_strong. | |
45 | B0=(33.-2.*NF)/6. | |
46 | ALGQ=LOG(Q2EFF/ALAM2) | |
47 | IF(MSTU(111).EQ.1) THEN | |
48 | ULALPS_HIJING=PARU(2)/(B0*ALGQ) | |
49 | ELSE | |
50 | B1=(153.-19.*NF)/6. | |
51 | ULALPS_HIJING=PARU(2)/(B0*ALGQ)*(1.-B1*LOG(ALGQ)/(B0**2*ALGQ)) | |
52 | ENDIF | |
53 | MSTU(118)=NF | |
54 | PARU(118)=ULALPS_HIJING | |
55 | ||
56 | RETURN | |
57 | END |