]>
Commit | Line | Data |
---|---|---|
fe4da5cc | 1 | * |
2 | * $Id$ | |
3 | * | |
4 | * $Log$ | |
5 | * Revision 1.1.1.2 1996/10/30 08:29:19 cernlib | |
6 | * Version 7.04 | |
7 | * | |
8 | * Revision 1.1.1.1 1996/04/12 15:29:29 plothow | |
9 | * Version 7.01 | |
10 | * | |
11 | * | |
12 | #include "pdf/pilot.h" | |
13 | C---------------------------------------------------------------------- | |
14 | FUNCTION SEADIS(XX) | |
15 | DIMENSION | |
16 | + SQ2(2:3) | |
17 | COMMON/W5051I1/ | |
18 | + GPSI(2:3),GAA(2:3),GPSIA(2:3),GAPSI(2:3),GMINUS(2:3), | |
19 | 1 AL(2:3),BET(2:3),GPLUS(2:3), | |
20 | 2 VNQ02(2:3),SNQ02(2:3),GNQ02(2:3) | |
21 | COMMON/W5051I2/ Q2L | |
22 | *- sea quarks' momenta | |
23 | PARAMETER (Q02=0.5,ALAM2=0.045*0.045) | |
24 | ||
25 | DATA ZERO/0./,HALF/0.5/,ONE/1./,TWO/2./ | |
26 | DATA PTF/0.25/,PSF/0.75/ | |
27 | DATA IRAZ/0/, Q2OLD/-1./ | |
28 | SAVE SQ2, BSQ2 | |
29 | ||
30 | IF(Q2L.NE.Q2OLD) THEN | |
31 | Q2OLD = Q2L | |
32 | S = LOG( LOG( MAX(Q02,Q2L)/ALAM2)/ LOG( Q02/ALAM2)) | |
33 | DO I=2,3 | |
34 | D1 = SNQ02(I)*EXP(-GPSI(I)*S) | |
35 | D2PLUS=(ONE-AL(I))*(SNQ02(I)+VNQ02(I))-BET(I)*GNQ02(I) | |
36 | D2MINU=AL(I)*(SNQ02(I)+VNQ02(I))+BET(I)*GNQ02(I) | |
37 | D2 =D2PLUS*EXP(-GPLUS (I)*S) + | |
38 | 1 D2MINU*EXP(-GMINUS(I)*S) - | |
39 | 2 VNQ02(I)*EXP(-GPSI(I)*S) | |
40 | SQ2(I)=PSF*D2 + PTF*D1 | |
41 | ENDDO | |
42 | BSQ2=SQ2(2)/SQ2(3) - TWO | |
43 | ENDIF | |
44 | SEADIS=SQ2(2) * (ONE-XX)**BSQ2 * (BSQ2+ONE) | |
45 | RETURN | |
46 | END |