* * $Id: eicht2.F,v 1.1.1.1 1996/01/11 14:14:36 mclareni Exp $ * * $Log: eicht2.F,v $ * Revision 1.1.1.1 1996/01/11 14:14:36 mclareni * Cojets * * #include "cojets/pilot.h" FUNCTION EICHT2(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 2, 1986 DATA C1/ * 2.36660, 0.44530, 0.36600, 0.09467, 0.13412, 0.01661, * -3.16964, -1.79481, 0.03313, -0.28736, -0.09827, -0.07119, * 1.82343, 1.45736, -0.24652, 0.03739, 0.00609, 0.01814, * -1.03315, -0.98274, 0.21364, 0.11686, 0.05001, 0.01684, * 0.51327, 0.52594, -0.11727, -0.11388, -0.04988, -0.02021, * -0.28808, -0.31452, 0.05667, 0.09161, 0.04568, 0.01951/ DATA C2/ * 0.72374, -0.21888, -0.29951, -0.01909, -0.01477, 0.00025, * -0.53139, -0.24246, 0.32831, 0.11193, 0.02223, 0.00707, * 0.22891, 0.18904, -0.09859, -0.06900, -0.01747, -0.00508, * -0.10415, -0.10839, 0.02108, 0.02975, 0.00983, 0.00283, * 0.04394, 0.05116, -0.00141, -0.01055, -0.00423, -0.00127, * -0.01991, -0.02539, -0.00278, 0.00343, 0.00172, 0.00055/ DATA C3/ * 0.35776, -0.08622, -0.14799, -0.01840, -0.00782, -0.00045, * -0.29250, -0.13038, 0.16959, 0.06243, 0.01353, 0.00375, * 0.13182, 0.10408, -0.05486, -0.03872, -0.01038, -0.00285, * -0.06162, -0.06143, 0.01303, 0.01740, 0.00594, 0.00167, * 0.02643, 0.02957, -0.00149, -0.00645, -0.00263, -0.00077, * -0.01218, -0.01497, -0.00126, 0.00224, 0.00112, 0.00035/ DATA C4/ * 0.10077, -0.07100, 0.01973, -0.00571, 0.00293, -0.00099, * -0.05271, -0.01823, 0.01792, -0.00658, 0.00175, -0.00155, * 0.01220, 0.01763, -0.00869, -0.00088, -0.00116, -0.00021, * -0.00119, -0.00718, 0.00236, 0.00189, 0.00077, 0.00041, * -0.00091, 0.00204, -0.00031, -0.00105, -0.00040, -0.00024, * 0.00119, -0.00017, -0.00020, 0.00042, 0.00017, 0.00010/ DATA C5/ * 0.06478, -0.04537, 0.01643, -0.00349, 0.00271, -0.00067, * -0.02223, -0.02126, 0.01247, -0.00629, 0.00112, -0.00144, * -0.00134, 0.01362, -0.00613, -0.00079, -0.00090, -0.00020, * 0.00508, -0.00361, 0.00170, 0.00183, 0.00068, 0.00040, * -0.00358, 0.00006, -0.00026, -0.00105, -0.00038, -0.00023, * 0.00242, 0.00093, -0.00010, 0.00045, 0.00017, 0.00011/ DATA C6/ * 0.00998, -0.01945, 0.01055, -0.00687, 0.00186, -0.00156, * 0.00570, -0.01203, 0.00625, -0.00486, 0.00131, -0.00137, * -0.00449, 0.00799, -0.00417, 0.00205, -0.00044, 0.00033, * 0.00147, -0.00248, 0.00146, -0.00057, 0.00012, -0.00001, * -0.00009, 0.00015, -0.00032, 0.00012, -0.00006, -0.00004, * -0.00042, 0.00076, -0.00014, 0.00004, 0.00007, 0.00005/ DATA C7/ * 0.00898, -0.01459, 0.00751, -0.00441, 0.00131, -0.00107, * 0.00597, -0.00944, 0.00480, -0.00302, 0.00091, -0.00085, * -0.00305, 0.00444, -0.00210, 0.00085, -0.00024, 0.00014, * 0.00053, -0.00130, 0.00056, -0.00027, 0.00003, -0.00002, * 0.00020, 0.00014, -0.00011, 0.00010, 0.00000, 0.00000, * -0.00026, 0.00032, 0.00000, -0.00003, 0.00001, -0.00001/ DATA C8/ * 0.00426, -0.00753, 0.00383, -0.00268, 0.00076, -0.00073, * 0.00364, -0.00605, 0.00303, -0.00209, 0.00059, -0.00060, * -0.00092, 0.00171, -0.00082, 0.00050, -0.00012, 0.00010, * -0.00005, -0.00016, 0.00013, -0.00009, 0.00003, 0.00000, * 0.00013, -0.00021, -0.00001, -0.00002, -0.00002, -0.00001, * -0.00008, 0.00018, -0.00005, 0.00002, 0.00000, 0.00000/ DATA D1/ * 30.35544,-40.62373, 15.77920, -3.69929, 0.60204, -0.07031, * 26.99749,-41.67159, 17.69848, -4.80436, 0.78620, -0.10596, * -1.90878, 1.35687, 1.12665, -0.71813, 0.22324, -0.02481, * -0.24877, 0.97811, -0.81272, 0.20935, -0.02997, -0.00471, * 0.25061, -0.54273, 0.26717, -0.03103, -0.00180, 0.00287, * -0.11277, 0.20870, -0.06972, -0.00248, 0.00263, -0.00084/ DATA D2/ * 0.24099, 0.28839, 0.09369, 0.01900, 0.00253, 0.00024, * 0.01765, -0.00922, -0.03037, -0.02085, -0.00844, -0.00281, * -0.00645, -0.00526, 0.00172, 0.00311, 0.00183, 0.00087, * 0.00212, 0.00232, 0.00026, -0.00049, -0.00039, -0.00023, * -0.00069, -0.00082, -0.00020, 0.00007, 0.00009, 0.00006, * 0.00024, 0.00031, 0.00011, 0.00000, -0.00002, -0.00002/ DATA D3/ * 0.12633, 0.13341, 0.03732, 0.00707, 0.00126, 0.00034, * 0.00366, -0.01357, -0.01795, -0.01031, -0.00388, -0.00128, * -0.00210, -0.00036, 0.00205, 0.00192, 0.00098, 0.00044, * 0.00077, 0.00054, -0.00024, -0.00039, -0.00024, -0.00013, * -0.00026, -0.00023, 0.00002, 0.00009, 0.00006, 0.00004, * 0.00009, 0.00010, 0.00002, -0.00002, -0.00002, -0.00001/ DATA D4/ * 1.08105, -1.18897, 0.38685, -0.08617, 0.01115, -0.00118, * 0.99169, -1.39649, 0.49978, -0.11585, 0.01674, -0.00172, * 0.05099, -0.13384, 0.09173, -0.02885, 0.00589, -0.00065, * -0.03178, 0.05703, -0.02070, 0.00244, 0.00011, -0.00009, * 0.00897, -0.01392, 0.00205, 0.00065, -0.00023, 0.00002, * -0.00234, 0.00301, 0.00050, -0.00039, 0.00006, 0.00001/ DATA D5/ * 0.98680, -1.17126, 0.39398, -0.08459, 0.01124, -0.00125, * 1.00058, -1.38342, 0.50438, -0.11521, 0.01658, -0.00183, * 0.04928, -0.13678, 0.09021, -0.02935, 0.00580, -0.00066, * -0.03133, 0.05785, -0.02023, 0.00263, 0.00016, -0.00008, * 0.00884, -0.01416, 0.00190, 0.00058, -0.00025, 0.00001, * -0.00230, 0.00308, 0.00055, -0.00037, 0.00007, 0.00001/ DATA D6/ * 0.86983, -1.13053, 0.38364, -0.08111, 0.01048, -0.00130, * 0.96262, -1.32081, 0.48537, -0.10910, 0.01583, -0.00170, * 0.03057, -0.10878, 0.08022, -0.02676, 0.00559, -0.00056, * -0.02845, 0.05164, -0.01918, 0.00221, -0.00004, -0.00015, * 0.01311, -0.01751, 0.00331, 0.00051, -0.00012, 0.00005, * -0.00859, 0.00838, -0.00092, -0.00026, 0.00001, -0.00001/ DATA D7/ * 0.86718, -1.17370, 0.42652, -0.09252, 0.01244, -0.00146, * 0.85004, -1.19366, 0.46304, -0.10831, 0.01614, -0.00183, * -0.02241, -0.00563, 0.02815, -0.01425, 0.00352, -0.00043, * -0.00073, 0.00803, -0.00578, 0.00138, -0.00013, -0.00004, * 0.00346, -0.00538, 0.00196, -0.00021, 0.00001, 0.00001, * -0.00485, 0.00595, -0.00189, 0.00026, -0.00003, 0.00000/ DATA D8/ * 0.71455, -1.00730, 0.39320, -0.09246, 0.01366, -0.00154, * 0.68561, -0.98276, 0.39768, -0.09795, 0.01540, -0.00179, * -0.03053, 0.02758, 0.00215, -0.00488, 0.00164, -0.00025, * 0.00092, 0.00042, -0.00134, 0.00046, -0.00008, -0.00001, * 0.00423, -0.00566, 0.00214, -0.00043, 0.00006, 0.00000, * -0.00389, 0.00500, -0.00174, 0.00033, -0.00004, 0.00000/ DATA ICALL/0/ IF(ICALL.GT.0) GO TO 10 ICALL=1 ALAMB=.290 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)=7.4474 TM(7,2)=TM(7,1) TM(8,1)=10.8097 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 EICHT2=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) EICHT2=0. DO 101 MX=1,6 DO 101 MQ=1,6 101 EICHT2=EICHT2+C(MX,MQ,IP)*TCX(MX)*TCQ(MQ) EICHT2=EICHT2*(1.-X)**LA(IP) IF(EICHT2.LT.0.) EICHT2=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) EICHT2=0. DO 201 MX=1,6 DO 201 MQ=1,6 201 EICHT2=EICHT2+D(MX,MQ,IP)*TCX(MX)*TCQ(MQ) EICHT2=EICHT2*(1.-X)**LA(IP) IF(EICHT2.LT.0.) EICHT2=0. RETURN END