| 1 | SUBROUTINE UTQVAL(Q,NEND)
|
|---|
| 2 |
|
|---|
| 3 | C-----------------------------------------------------------------------
|
|---|
| 4 | C UT(ILITY ROUTINE) VAL(ENCE QUARK STRUCTURE FUNCTION)
|
|---|
| 5 | C
|
|---|
| 6 | C VALENCE QUARK STRUCTURE FUNCTION
|
|---|
| 7 | C RETURNS INTEGRAL (XVA(1)->XVA(I)) FU(Z) DZ
|
|---|
| 8 | C THIS INTEGRAL IS ONLY CALCULATED FOR SMALL VALUES OF XVA UP TO 25
|
|---|
| 9 | C TIMES THE VALUE OF XCUT. FOR LARGER VALUES THE TABULATED VALUES OF
|
|---|
| 10 | C DATASET 'VENUSDAT' ARE TAKEN AND CORRECTED BY THE CONSTANT SHIFT
|
|---|
| 11 | C DELTA0 (FOR HADRONS) OR DELTA1 (FOR PIONS).
|
|---|
| 12 | C
|
|---|
| 13 | C THIS SUBROUTINE IS CALLED FROM VENLNK
|
|---|
| 14 | C
|
|---|
| 15 | C DESIGN : D. HECK IK3 FZK KARLSRUHE
|
|---|
| 16 | C-----------------------------------------------------------------------
|
|---|
| 17 |
|
|---|
| 18 | *KEEP,RUNPAR.
|
|---|
| 19 | COMMON /RUNPAR/ FIXHEI,THICK0,HILOECM,HILOELB,
|
|---|
| 20 | * STEPFC,NRRUN,NSHOW,PATAPE,MONIIN,
|
|---|
| 21 | * MONIOU,MDEBUG,NUCNUC,
|
|---|
| 22 | * CETAPE,
|
|---|
| 23 | * SHOWNO,ISHW,NOPART,NRECS,NBLKS,MAXPRT,NDEBDL,
|
|---|
| 24 | * N1STTR,MDBASE,
|
|---|
| 25 | * DEBDEL,DEBUG,FDECAY,FEGS,FIRSTI,FIXINC,FIXTAR,
|
|---|
| 26 | * FIX1I,FMUADD,FNKG,FPRINT,FDBASE
|
|---|
| 27 | * ,GHEISH,GHESIG
|
|---|
| 28 | COMMON /RUNPAC/ DSN,HOST,USER
|
|---|
| 29 | DOUBLE PRECISION FIXHEI,THICK0,HILOECM,HILOELB
|
|---|
| 30 | REAL STEPFC
|
|---|
| 31 | INTEGER NRRUN,NSHOW,PATAPE,MONIIN,MONIOU,MDEBUG,NUCNUC,
|
|---|
| 32 | * SHOWNO,ISHW,NOPART,NRECS,NBLKS,MAXPRT,NDEBDL,
|
|---|
| 33 | * N1STTR,MDBASE
|
|---|
| 34 | INTEGER CETAPE
|
|---|
| 35 | CHARACTER*79 DSN
|
|---|
| 36 | CHARACTER*20 HOST,USER
|
|---|
| 37 |
|
|---|
| 38 | LOGICAL DEBDEL,DEBUG,FDECAY,FEGS,FIRSTI,FIXINC,FIXTAR,
|
|---|
| 39 | * FIX1I,FMUADD,FNKG,FPRINT,FDBASE
|
|---|
| 40 | * ,GHEISH,GHESIG
|
|---|
| 41 | *KEND.
|
|---|
| 42 |
|
|---|
| 43 | PARAMETER (NSTRU=2049)
|
|---|
| 44 | COMMON /CIPIO/ IPIO
|
|---|
| 45 | COMMON /FILES/ IFCH,IFDT,IFHI,IFMT,IFOP
|
|---|
| 46 | COMMON /PARO1/ AMPRIF,AMSIAC,BMAXIM,BMINIM,CORE,CUTMSQ,CUTMSS
|
|---|
| 47 | * ,DELMSS,DELREM,FCTRMX,GAUMX,OVERLP,PAREA,PDIQUA
|
|---|
| 48 | * ,PHARD,PSPINL,PSPINH,PISPN,PTF,PTH,PTMX,PTQ,PUD
|
|---|
| 49 | * ,PVALEN,QSEPC,QSETC,QMUST,QVAPC,QVATC,RADIAC
|
|---|
| 50 | * ,RADIAS,RSTRAS,SIGJ,SIGPPI,TAUMAX,TAUMIN
|
|---|
| 51 | * ,TAUMX,TAUNLL,TENSN,THEMAS,WPROJ,WTARG,WTMINI
|
|---|
| 52 | * ,WTSTEP,XCUT
|
|---|
| 53 | * ,IAQU,IFRADE,IOJINT,IOPBRK,IOPENT,IOPENU
|
|---|
| 54 | * ,IOPTF,IOPTQ,IRESCL,IWCENT,KENTRO,KO1KO2
|
|---|
| 55 | * ,LABSYS,MAXRES,NCLEAN,NCOLMX,NDECAW,NEQMN,NEQMX
|
|---|
| 56 | * ,NSTTAU,NTRYMX,NUMTAU
|
|---|
| 57 | COMMON /STRU2/ DELTA0,DELTA1,QSEH(NSTRU),QSEPI(NSTRU)
|
|---|
| 58 | * ,QVAH(NSTRU),QVAPI(NSTRU),XSE(NSTRU),XVA(NSTRU)
|
|---|
| 59 |
|
|---|
| 60 | DIMENSION Y0(9),Y1(9),Q(NEND)
|
|---|
| 61 | C-----------------------------------------------------------------------
|
|---|
| 62 |
|
|---|
| 63 | IF ( DEBUG ) WRITE(MDEBUG,*) 'UTQVAL: IPIO,NEND=',IPIO,NEND
|
|---|
| 64 |
|
|---|
| 65 | XCUT2 = XCUT**2
|
|---|
| 66 | Q(1) = 0.
|
|---|
| 67 | Z = XVA(1)
|
|---|
| 68 | DENOMI = 1. / SQRT(Z**2 + XCUT2)
|
|---|
| 69 |
|
|---|
| 70 | IF ( IPIO .EQ. 0 ) THEN
|
|---|
| 71 | C CALCULATE THE FIRST NEND VALUES OF STRUCTURE FUNCTION FOR HADRONS
|
|---|
| 72 | Y0(1) = 0.
|
|---|
| 73 | DO 3 I=2,NEND
|
|---|
| 74 | FACT = (XVA(I) - Z) * 0.125
|
|---|
| 75 | DO 2 J=2,8
|
|---|
| 76 | Z = Z + FACT
|
|---|
| 77 | DENOMI = 1. / SQRT(Z**2 + XCUT2)
|
|---|
| 78 | Y0(J) = (1.-Z)**3.46 * Z**.419 * (2.74793064*Z + 0.62452969)
|
|---|
| 79 | * * DENOMI
|
|---|
| 80 | 2 CONTINUE
|
|---|
| 81 | Z = XVA(I)
|
|---|
| 82 | DENOMI = 1. / SQRT(Z**2 + XCUT2)
|
|---|
| 83 | Y0(9) = (1.-Z)**3.46 * Z**.419 * (2.74793064*Z + 0.62452969)
|
|---|
| 84 | * * DENOMI
|
|---|
| 85 | C INTEGRATION AFTER BODE'S RULE (ABRAMOWITZ + STEGUN, HANDBOOK OF
|
|---|
| 86 | C MATHEMATICAL FUNCTIONS, DOVER PUBLICATIONS (1970), FORMULA 25.4.18)
|
|---|
| 87 | Q(I) = 2.8218694E-4 * FACT * ( 989. * (Y0(1) + Y0(9))
|
|---|
| 88 | * + 5888. * (Y0(2) + Y0(8)) - 928. * (Y0(3) + Y0(7))
|
|---|
| 89 | * + 10496. * (Y0(4) + Y0(6)) - 4540. * Y0(5) )
|
|---|
| 90 | * + Q(I-1)
|
|---|
| 91 | Y0(1) = Y0(9)
|
|---|
| 92 | 3 CONTINUE
|
|---|
| 93 |
|
|---|
| 94 | ELSE
|
|---|
| 95 | C CALCULATE THE FIRST NEND VALUES OF STRUCTURE FUNCTION FOR PIONS
|
|---|
| 96 | Y1(1) = 0.
|
|---|
| 97 | DO 5 I=2,NEND
|
|---|
| 98 | FACT = (XVA(I) - Z) * 0.125
|
|---|
| 99 | DO 4 J=2,8
|
|---|
| 100 | Z = Z + FACT
|
|---|
| 101 | DENOMI = 1. / SQRT(Z**2 + XCUT2)
|
|---|
| 102 | Y1(J) = (1.-Z)**0.7 * Z**.4 * DENOMI
|
|---|
| 103 | 4 CONTINUE
|
|---|
| 104 | Z = XVA(I)
|
|---|
| 105 | DENOMI = 1. / SQRT(Z**2 + XCUT2)
|
|---|
| 106 | Y1(9) = (1.-Z)**0.7 * Z**.4 * DENOMI
|
|---|
| 107 | C INTEGRATION AFTER BODE'S RULE (ABRAMOWITZ + STEGUN, HANDBOOK OF
|
|---|
| 108 | C MATHEMATICAL FUNCTIONS, DOVER PUBLICATIONS (1970), FORMULA 25.4.18)
|
|---|
| 109 | Q(I) = 2.8218694E-4 * FACT * ( 989. * (Y1(1) + Y1(9))
|
|---|
| 110 | * + 5888. * (Y1(2) + Y1(8)) - 928. * (Y1(3) + Y1(7))
|
|---|
| 111 | * + 10496. * (Y1(4) + Y1(6)) - 4540. * Y1(5) )
|
|---|
| 112 | * * 0.1730725 + Q(I-1)
|
|---|
| 113 | Y1(1) = Y1(9)
|
|---|
| 114 | 5 CONTINUE
|
|---|
| 115 | ENDIF
|
|---|
| 116 |
|
|---|
| 117 | RETURN
|
|---|
| 118 | END
|
|---|
