* * $Id: eicht1.F,v 1.1.1.1 1996/01/11 14:14:35 mclareni Exp $ * * $Log: eicht1.F,v $ * Revision 1.1.1.1 1996/01/11 14:14:35 mclareni * Cojets * * #include "cojets/pilot.h" FUNCTION EICHT1(X,Q2,IP) C ************************ C-- EICHTEN ET AL. EVOLVED PARTON DENSITIES TIMES X -- SET 2 C-- IP=1,8 -- GLUON, UV, DV, US, S, C, B, T C-- REV.MOD.PHYS., 56 (1984) 579, ERRATUM 58 (1986) 1065. #if defined(CERNLIB_SINGLE) IMPLICIT REAL (A-H,O-Z) #endif #if defined(CERNLIB_DOUBLE) IMPLICIT DOUBLE PRECISION (A-H,O-Z) #endif #include "cojets/itapes.inc" DIMENSION TCX(6),TCQ(6) DIMENSION LA(8),C(6,6,8),D(6,6,8),TM(8,2),AM(8) DIMENSION C1(6,6),C2(6,6),C3(6,6),C4(6,6),C5(6,6),C6(6,6),C7(6,6) DIMENSION C8(6,6),D1(6,6),D2(6,6),D3(6,6),D4(6,6),D5(6,6),D6(6,6) DIMENSION D7(6,6),D8(6,6) DATA LA/6,3,4,5*7/,AM/6*0.,5.5,30./ C-- DATA EHLQ SET 1, 1986 DATA C1/ * 0.94819, -0.95779, 0.10085, -0.10510, 0.03456, -0.03054, * -0.96265, 0.53790, 0.33684, -0.09525, 0.01488, -0.02051, * 0.43004, -0.08306, -0.33719, 0.04902, -0.00916, 0.01041, * -0.19249, -0.01790, 0.21830, 0.00749, 0.00414, -0.00186, * 0.08183, 0.01926, -0.10718, -0.01944, -0.00277, -0.00052, * -0.03884, -0.01234, 0.05410, 0.01879, 0.00335, 0.00104/ DATA C2/ * 0.76772, -0.20874, -0.33026, -0.02517, -0.01570, -0.00010, * -0.53259, -0.26612, 0.32007, 0.11918, 0.02434, 0.00762, * 0.21618, 0.18812, -0.08375, -0.06515, -0.01743, -0.00504, * -0.09211, -0.09952, 0.01373, 0.02506, 0.00877, 0.00255, * 0.03670, 0.04409, 0.00096, -0.00796, -0.00342, -0.00105, * -0.01549, -0.02026, -0.00306, 0.00222, 0.00124, 0.00041/ DATA C3/ * 0.38130, -0.08090, -0.16336, -0.02185, -0.00843, -0.00062, * -0.29475, -0.14348, 0.16650, 0.06638, 0.01473, 0.00408, * 0.12518, 0.10422, -0.04722, -0.03683, -0.01038, -0.00286, * -0.05478, -0.05678, 0.00890, 0.01484, 0.00534, 0.00152, * 0.02220, 0.02567, -0.00003, -0.00497, -0.00216, -0.00065, * -0.00953, -0.01204, -0.00151, 0.00151, 0.00083, 0.00027/ DATA C4/ * 0.06870, -0.06861, 0.02973, -0.00540, 0.00378, -0.00097, * -0.01802, 0.00014, 0.00649, -0.00854, 0.00122, -0.00175, * -0.00465, 0.00148, -0.00593, 0.00060, -0.00103, -0.00008, * 0.00644, 0.00257, 0.00283, 0.00115, 0.00071, 0.00033, * -0.00393, -0.00254, -0.00116, -0.00077, -0.00036, -0.00019, * 0.00234, 0.00193, 0.00053, 0.00037, 0.00016, 0.00009/ DATA C5/ * 0.04968, -0.04173, 0.02102, -0.00327, 0.00324, -0.00067, * -0.00615, -0.01294, 0.00674, -0.00689, 0.00090, -0.00151, * -0.00858, 0.00505, -0.00490, -0.00016, -0.00094, -0.00015, * 0.00784, 0.00151, 0.00222, 0.00140, 0.00070, 0.00035, * -0.00441, -0.00222, -0.00089, -0.00085, -0.00036, -0.00020, * 0.00252, 0.00184, 0.00041, 0.00039, 0.00016, 0.00009/ DATA C6/ * 0.00927, -0.01817, 0.00959, -0.00639, 0.00169, -0.00154, * 0.00571, -0.01188, 0.00609, -0.00465, 0.00124, -0.00131, * -0.00396, 0.00710, -0.00359, 0.00184, -0.00039, 0.00034, * 0.00112, -0.00196, 0.00112, -0.00048, 0.00010, -0.00004, * 0.00004, -0.00003, -0.00018, 0.00009, -0.00005, -0.00002, * -0.00042, 0.00073, -0.00016, 0.00005, 0.00005, 0.00005/ DATA C7/ * 0.00901, -0.01401, 0.00715, -0.00413, 0.00126, -0.00104, * 0.00628, -0.00932, 0.00478, -0.00289, 0.00091, -0.00082, * -0.00293, 0.00409, -0.00189, 0.00076, -0.00023, 0.00014, * 0.00039, -0.00120, 0.00044, -0.00025, 0.00002, -0.00002, * 0.00026, 0.00014, -0.00008, 0.00010, 0.00001, 0.00001, * -0.00026, 0.00032, 0.00001, -0.00001, 0.00001, -0.00001/ DATA C8/ * 0.00441, -0.00748, 0.00377, -0.00258, 0.00073, -0.00071, * 0.00384, -0.00605, 0.00303, -0.00203, 0.00058, -0.00059, * -0.00088, 0.00166, -0.00075, 0.00047, -0.00010, 0.00010, * -0.00008, -0.00015, 0.00012, -0.00009, 0.00003, 0.00000, * 0.00013, -0.00022, -0.00002, -0.00002, -0.00002, -0.00002, * -0.00007, 0.00019, -0.00004, 0.00002, 0.00000, 0.00000/ DATA D1/ * 29.47734, -39.02468, 14.63570, -3.33516, 0.50538, -0.05915, * 25.58960, -39.54527, 16.61420, -4.29861, 0.69036, -0.08243, * -1.66291, 1.17624, 1.11844, -0.70986, 0.19481, -0.02404, * -0.21679, 0.81705, -0.71688, 0.18507, -0.01924, -0.00325, * 0.20880, -0.43547, 0.22391, -0.02446, -0.00362, 0.00191, * -0.09097, 0.16009, -0.05681, -0.00250, 0.00258, -0.00047/ DATA D2/ * 0.23946, 0.29055, 0.09778, 0.02149, 0.00344, 0.00050, * 0.01751, -0.00609, -0.02687, -0.01916, -0.00797, -0.00275, * -0.00576, -0.00504, 0.00108, 0.00249, 0.00153, 0.00075, * 0.00174, 0.00196, 0.00030, -0.00034, -0.00029, -0.00018, * -0.00053, -0.00064, -0.00017, 0.00004, 0.00006, 0.00004, * 0.00017, 0.00022, 0.00008, 0.00001, -0.00001, -0.00001/ DATA D3/ * 0.12613, 0.13542, 0.03958, 0.00824, 0.00166, 0.00045, * 0.00389, -0.01159, -0.01625, -0.00961, -0.00371, -0.00126, * -0.00191, -0.00056, 0.00159, 0.00159, 0.00084, 0.00039, * 0.00064, 0.00049, -0.00015, -0.00029, -0.00018, -0.00010, * -0.00020, -0.00019, 0.00000, 0.00006, 0.00004, 0.00003, * 0.00007, 0.00008, 0.00002, -0.00001, -0.00001, -0.00001/ DATA D4/ * 1.01386, -1.10585, 0.33739, -0.07444, 0.00885, -0.00087, * 0.92334, -1.28541, 0.44755, -0.09786, 0.01419, -0.00112, * 0.04888, -0.12708, 0.08606, -0.02608, 0.00478, -0.00060, * -0.02691, 0.04887, -0.01771, 0.00162, 0.00025, -0.00006, * 0.00704, -0.01113, 0.00159, 0.00070, -0.00020, 0.00000, * -0.00171, 0.00229, 0.00038, -0.00035, 0.00004, 0.00001/ DATA D5/ * 0.92351, -1.08483, 0.34642, -0.07210, 0.00914, -0.00091, * 0.93146, -1.27376, 0.45122, -0.09775, 0.01380, -0.00131, * 0.04739, -0.12960, 0.08482, -0.02642, 0.00476, -0.00057, * -0.02653, 0.04953, -0.01735, 0.00175, 0.00028, -0.00006, * 0.00694, -0.01132, 0.00148, 0.00065, -0.00021, 0.00000, * -0.00168, 0.00234, 0.00042, -0.00034, 0.00005, 0.00001/ DATA D6/ * 0.80983, -1.04168, 0.33980, -0.06824, 0.00876, -0.00090, * 0.89606, -1.21708, 0.43386, -0.09287, 0.01304, -0.00129, * 0.03058, -0.10402, 0.07604, -0.02415, 0.00460, -0.00050, * -0.02451, 0.04432, -0.01651, 0.00143, 0.00012, -0.00010, * 0.01122, -0.01457, 0.00268, 0.00058, -0.00012, 0.00003, * -0.00773, 0.00733, -0.00076, -0.00024, 0.00001, 0.00000/ DATA D7/ * 0.80288, -1.07532, 0.37920, -0.07843, 0.01007, -0.00109, * 0.79033, -1.09887, 0.41532, -0.09301, 0.01317, -0.00141, * -0.01704, -0.01130, 0.02882, -0.01341, 0.00304, -0.00036, * -0.00072, 0.00723, -0.00516, 0.00108, -0.00005, -0.00004, * 0.00305, -0.00461, 0.00166, -0.00013, -0.00001, 0.00001, * -0.00436, 0.00523, -0.00161, 0.00020, -0.00002, 0.00000/ DATA D8/ * 0.66233, -0.92481, 0.35193, -0.07930, 0.01110, -0.00118, * 0.63797, -0.90619, 0.35816, -0.08479, 0.01265, -0.00139, * -0.02581, 0.02125, 0.00419, -0.00498, 0.00149, -0.00021, * 0.00071, 0.00053, -0.00127, 0.00039, -0.00005, -0.00001, * 0.00385, -0.00506, 0.00186, -0.00035, 0.00004, 0.00000, * -0.00353, 0.00446, -0.00150, 0.00027, -0.00003, 0.00000/ DATA ICALL/0/ IF(ICALL.GT.0) GO TO 10 ICALL=1 ALAMB=.200 ALAMB2=ALAMB*ALAMB TX=LOG(1.E8/ALAMB2) TMM=LOG(5./ALAMB2) DO 1 M1=1,6 DO 1 M2=1,2 1 TM(M1,M2)=TMM TM(7,1)=8.1905 TM(7,2)=TM(7,1) TM(8,1)=11.5528 TM(8,2)=TM(8,1) DO 2 M1=1,6 DO 2 M2=1,6 C(M1,M2,1)=C1(M1,M2) C(M1,M2,2)=C2(M1,M2) C(M1,M2,3)=C3(M1,M2) C(M1,M2,4)=C4(M1,M2) C(M1,M2,5)=C5(M1,M2) C(M1,M2,6)=C6(M1,M2) C(M1,M2,7)=C7(M1,M2) C(M1,M2,8)=C8(M1,M2) D(M1,M2,1)=D1(M1,M2) D(M1,M2,2)=D2(M1,M2) D(M1,M2,3)=D3(M1,M2) D(M1,M2,4)=D4(M1,M2) D(M1,M2,5)=D5(M1,M2) D(M1,M2,6)=D6(M1,M2) D(M1,M2,7)=D7(M1,M2) D(M1,M2,8)=D8(M1,M2) 2 CONTINUE 10 CONTINUE IF(IP.LT.7) GO TO 20 EICHT1=0. B2=1.-4.*AM(IP)**2/(Q2*(1.-X)) IF(B2.LT.0.) RETURN 20 CONTINUE T=LOG(Q2/ALAMB2) IF(X.LE..1) GO TO 200 X1=(2.*X-1.1)/0.9 T1=(2.*T-(TX+TM(IP,1)))/(TX-TM(IP,1)) CALL CJEB(X1,TCX) CALL CJEB(T1,TCQ) EICHT1=0. DO 101 MX=1,6 DO 101 MQ=1,6 101 EICHT1=EICHT1+C(MX,MQ,IP)*TCX(MX)*TCQ(MQ) EICHT1=EICHT1*(1.-X)**LA(IP) IF(EICHT1.LT.0.) EICHT1=0. RETURN 200 Y=(2.*LOG(X)+11.51293)/6.90776 T1=(2.*T-(TX+TM(IP,2)))/(TX-TM(IP,2)) CALL CJEB(Y,TCX) CALL CJEB(T1,TCQ) EICHT1=0. DO 201 MX=1,6 DO 201 MQ=1,6 201 EICHT1=EICHT1+D(MX,MQ,IP)*TCX(MX)*TCQ(MQ) EICHT1=EICHT1*(1.-X)**LA(IP) IF(EICHT1.LT.0.) EICHT1=0. RETURN END