5 * Revision 1.1.1.2 1996/10/30 08:31:56 cernlib
8 * Revision 1.2 1996/10/28 15:13:59 pdflib
9 * PDFLIB new version 7.04
11 * Revision 1.1.1.1 1996/04/12 15:30:06 plothow
15 #include "pdf/pilot.h"
17 C-----------------------------------------------------------------------
19 SUBROUTINE STRUCTM(X,SCALE,UPV,DNV,USEA,DSEA,STR,CHM,BOT,TOP,GL)
21 C *********************************************************************
23 C * Main steering routine for all sets of structure functions *
26 C * Input: X = x value of parton *
27 C * SCALE = QCD scale in GeV *
29 C * Output: UPV = up valence quark *
30 C * DNV = down valence quark *
31 C * USEA = sea (up_bar) *
32 C * DSEA = sea (down_bar) *
33 C * STR = strange quark *
34 C * CHM = charm quark *
35 C * BOT = bottom quark *
40 C * The variables NPTYPE, NGROUP and NSET should be, *
41 C * the variables NFL, LO, TMAS and *
42 C * QCDL4, QCDL5, XMIN, XMAX, Q2MIN, Q2MAX *
43 C * could be provided by the user via a call to the *
44 C * subroutine PDFSET at the initialization phase, where *
46 C * NPTYPE = Particle type *
47 C * (number or character string, *
48 C * 1,2,3 or 'NU','PI','PH') *
49 C * of desired structure functions set *
50 C * (Default: NPTYPE = 1 or 'NU' *
51 C * NGROUP = author group *
52 C * (number or character string, *
53 C * 1 to 7 or i.e. 'DFLM','MRS','MT','GRV',etc) *
54 C * of desired structure functions set *
55 C * (Default: NGROUP = 5 or 'GRV') *
56 C * NSET = number of desired structure functions set *
57 C * (Default: NSET = 3) *
58 C * NFL = desired number of flavours for alpha(s) *
59 C * (Default: NFL = 5) *
60 C * LO = order of alpha(s) calculation *
61 C * (Default: LO = 2) *
62 C * TMAS = top quark mass in GeV/c**2 (optional) *
63 C * (Default: TMAS = 100.0D0) *
64 C * QCDL4 = QCD scale in GeV for four flavours *
65 C * QCDL5 = QCD scale in GeV for five flavours *
66 C * corresponding to QCDL4 *
67 C * XMIN = minimal allowed x value *
68 C * XMAX = maximal allowed x value *
69 C * Q2MIN = minimal allowed Q**2 value *
70 C * Q2MAX = maximal allowed Q**2 value *
72 C * for each set of structure fuction. *
74 C * The internal COMMON blocks *
76 C * COMMON/W50511/ NPTYPE,NGROUP,NSET,MODE,NFL,LO,TMAS *
77 C * COMMON/W50512/ QCDL4,QCDL5 *
78 C * COMMON/W50513/ XMIN,XMAX,Q2MIN,Q2MAX *
80 C * are filled then by the subroutine PDFSET. *
83 C * Note: STRUCTM returns X * parton distribution function ! *
86 C * Author: H. Plothow-Besch *
87 C * CERN-PPE, CH - 1211 Geneva 23, Switzerland *
89 C * Please return any problems, questions, suggestions *
92 C *********************************************************************
95 #include "pdf/impdp.inc"
97 #include "pdf/w5051p1.inc"
98 #include "pdf/w5051p2.inc"
99 #include "pdf/w5051p7.inc"
100 #include "pdf/w50510.inc"
101 #include "pdf/w50511.inc"
102 #include "pdf/w50512.inc"
103 #include "pdf/w50513.inc"
104 #include "pdf/w50514.inc"
105 #include "pdf/w50514w.inc"
106 #include "pdf/w50515.inc"
107 #include "pdf/w50516.inc"
108 #include "pdf/w50517.inc"
109 #include "pdf/w50519.inc"
110 #include "pdf/w505120.inc"
111 #include "pdf/w505121.inc"
112 CHARACTER*20 PARM(NCHDIM)
113 #include "pdf/expdp.inc"
115 DATA ZEROD/0.D0/,ONED/1.D0/,TWOD/2.D0/
116 SAVE /W50514/, /W50514W/, /W50516/
118 #include "pdf/w50511c.inc"
120 C User wants new version (4.0 or bigger) of PDFLIB format
121 IF (IFLSET.NE.1) THEN
123 WRITE(N6,*) ' Warning : NO initialisation via PDFSET made !!'
125 + ' ALL Parameters set to default (Nucleon PDFs) !!'
138 CALL PDFSET(PARM,VAL)
141 C User wants old version (3.0 or less) of PDFLIB format
143 IF (MODE .GE.0 .AND. MODE .LE.MODEMX) THEN
149 NPTYPE = NPTYCR(MODE)
150 NGROUP = NGROCR(MODE)
152 C Check on validity of parameter values
153 IF (NPTYPE.LT.0 .OR. NGROUP.LT.0 .OR. NSET.LT.0) THEN
155 WRITE(N6,*) ' PDFLIB : MODE value INCORRECT, MODE = ',
157 WRITE(N6,*) ' Warning : MODE value set to OLD default !!'
194 IF (DX.LE.ZEROD .OR. DX.GE.ONED) THEN
195 WRITE(N6,*) ' PDFLIB : Nptype = ',NPTYPE,', Ngroup = ',
196 + NGROUP,', Nset = ',NSET
198 + ' Error : X value outside physical range , X = ',DX
201 WTXMIN = WTXMIN + PDFWGT
203 WXMIN = WXMIN + PDFWGT
205 WRITE(N6,*) ' PDFLIB : Nptype = ',NPTYPE,', Ngroup = ',
206 + NGROUP,', Nset = ',NSET
208 + ' Error : X value smaller allowed range , X = ',DX
211 WTXMAX = WTXMAX + PDFWGT
213 WXMAX = WXMAX + PDFWGT
215 WRITE(N6,*) ' PDFLIB : Nptype = ',NPTYPE,', Ngroup = ',
216 + NGROUP,', Nset = ',NSET
218 + ' Error : X value bigger allowed range , X = ',DX
222 WTQ2MIN = WTQ2MIN + PDFWGT
223 IF (DQ2.LT.Q2MIN) THEN
224 WQ2MIN = WQ2MIN + PDFWGT
226 WRITE(N6,*) ' PDFLIB : Nptype = ',NPTYPE,', Ngroup = ',
227 + NGROUP,', Nset = ',NSET
229 + ' Error : SCALE value smaller allowed range , Q**2 = ',DQ2
232 WTQ2MAX = WTQ2MAX + PDFWGT
233 IF (DQ2.GT.Q2MAX) THEN
234 WQ2MAX = WQ2MAX + PDFWGT
236 WRITE(N6,*) ' PDFLIB : Nptype = ',NPTYPE,', Ngroup = ',
237 + NGROUP,', Nset = ',NSET
239 + ' Error : SCALE value bigger allowed range , Q**2 = ',DQ2
243 IF (FIRST) WRITE(N6,*) ' '
244 IF (FIRST) WRITE(N6,*) ' '
245 C Call default structure functions set
246 IF (NPTYPE.EQ.0 .OR. NGROUP.EQ.0 .OR. NSET.EQ.0) THEN
247 C IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
248 C + ', Nset = ',NSET,
249 C + ', for GRV Set HO Structure Functions'
250 C CALL GRVHO(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DTOP,DGL)
253 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
255 + ', for MRS Set (G) (L255-MSb) Structure Functions'
256 CALL SFMRSG(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
260 C Start with NUCLEON structure functions
264 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
266 + ', for PRIVATE Structure Functions'
267 CALL STRPRIV(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
269 ELSEIF(NSET.EQ.2) THEN
270 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
272 + ', for BEBC Structure Functions'
273 CALL STRBEBC(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DGL)
279 ELSEIF(NSET.EQ.3) THEN
280 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
282 + ', for OR Structure Functions'
283 CALL STRUCOR(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DGL)
289 ELSEIF(NSET.EQ.4) THEN
290 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
292 + ', for BEP Structure Functions'
293 CALL STRBEP(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DGL)
299 ELSEIF(NSET.EQ.5) THEN
300 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
302 + ', for GHR Structure Functions'
303 CALL STRGHR(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DGL)
308 ELSEIF(NSET.EQ.6) THEN
309 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
311 + ', for DO Set 1 Structure Functions'
312 CALL STRDO1(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DGL)
317 ELSEIF(NSET.EQ.7) THEN
318 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
320 + ', for DO Set 2 Structure Functions'
321 CALL STRDO2(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DGL)
326 ELSEIF(NSET.EQ.8) THEN
327 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
329 + ', for EHLQ Set 1 Structure Functions'
330 CALL SFEHLQ1(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DTOP,DGL)
333 ELSEIF(NSET.EQ.9) THEN
334 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
336 + ', for EHLQ Set 2 Structure Functions'
337 CALL SFEHLQ2(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DTOP,DGL)
340 ELSEIF(NSET.EQ.10) THEN
341 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
343 + ', for New DO Set 1.1 Structure Functions'
344 CALL NEWDO1(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DGL)
349 ELSEIF(NSET.GT.11) THEN
350 WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
352 + ', Structure Functions not yet exsistent'
355 ELSEIF(NGROUP.EQ.2) THEN
357 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
359 + ', for DFLMsoft (Valence) Structure Functions'
360 CALL DFLM4(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DTOP,DGL)
363 ELSEIF(NSET.EQ.2) THEN
364 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
366 + ', for DFLMhard (Valence) Structure Functions'
367 CALL DFLM5(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DTOP,DGL)
370 ELSEIF(NSET.EQ.3) THEN
371 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
373 + ', for DFLMsoft (Gluon) Structure Functions'
374 CALL DFLM6(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DTOP,DGL)
377 ELSEIF(NSET.EQ.4) THEN
378 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
380 + ', for DFLMhard (Gluon) Structure Functions'
381 CALL DFLM7(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DTOP,DGL)
384 ELSEIF(NSET.EQ.5) THEN
385 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
387 + ', for DFLMaverage (LO) Structure Functions'
388 CALL DFLM8(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DTOP,DGL)
391 ELSEIF(NSET.EQ.6) THEN
392 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
394 + ', for DFLMaverage (DIS) Structure Functions'
395 CALL DFLM9(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DTOP,DGL)
398 ELSEIF(NSET.EQ.7) THEN
399 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
401 + ', for DFLM160 (DIS) Structure Functions'
402 CALL DFLM1(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DTOP,DGL)
405 ELSEIF(NSET.EQ.8) THEN
406 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
408 + ', for DFLM260 (DIS) Structure Functions'
409 CALL DFLM2(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DTOP,DGL)
412 ELSEIF(NSET.EQ.9) THEN
413 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
415 + ', for DFLM360 (DIS) Structure Functions'
416 CALL DFLM3(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DTOP,DGL)
419 ELSEIF(NSET.GT.10) THEN
420 WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
422 + ', Structure Functions not yet exsistent'
425 ELSEIF(NGROUP.EQ.3) THEN
427 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
429 + ', for MRS Set 1 (MSb) Structure Functions'
430 CALL STRMRS1(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DGL)
434 ELSEIF(NSET.EQ.2) THEN
435 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
437 + ', for MRS Set 2 (MSb) Structure Functions'
438 CALL STRMRS2(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DGL)
442 ELSEIF(NSET.EQ.3) THEN
443 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
445 + ', for MRS Set 3 (MSb) Structure Functions'
446 CALL STRMRS3(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DGL)
450 ELSEIF(NSET.EQ.4) THEN
451 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
453 + ', for MRS Set E Structure Functions'
454 CALL STRUCE(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DGL)
458 ELSEIF(NSET.EQ.5) THEN
459 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
461 + ', for MRS Set B Structure Functions'
462 CALL STRUCB(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DGL)
466 ELSEIF(NSET.EQ.6) THEN
467 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
469 + ', for MRS Set EP Structure Functions'
470 CALL STRCEP(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DGL)
474 ELSEIF(NSET.EQ.7) THEN
475 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
477 + ', for MRS Set BP Structure Functions'
478 CALL STRCBP(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DGL)
482 ELSEIF(NSET.EQ.8) THEN
483 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
485 + ', for HMRS1 Set E Structure Functions'
486 IF (FIRST) WRITE(N6,*)
487 +' Warning !! : in principle, set has been retracted by the auth
489 CALL STRH1E(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DGL)
493 C ELSEIF(NSET.EQ.9) THEN
494 C IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
495 C + ', Nset = ',NSET,
496 C + ', for HMRS1 Set B Structure Functions'
497 C IF (FIRST) WRITE(N6,*)
498 C +' Warning !! : in principle, set has been retracted by the auth
501 C CALL STRH1B(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DGL)
505 ELSEIF(NSET.EQ.9) THEN
506 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
508 + ', for KMRS Set B0 (190) Structure Functions (old)'
509 CALL STRKB0O(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DGL)
513 ELSEIF(NSET.EQ.10) THEN
514 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
516 + ', for HMRS2 Set E Structure Functions'
517 CALL STRH2E(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DGL)
521 ELSEIF(NSET.EQ.11) THEN
522 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
524 + ', for HMRS2 Set B Structure Functions'
525 CALL STRH2B(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DGL)
529 ELSEIF(NSET.EQ.12) THEN
530 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
532 + ', for HMRS Set E+ Structure Functions'
533 CALL STRH3EP(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DGL)
537 ELSEIF(NSET.EQ.13) THEN
538 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
540 + ', for HMRS Set E- Structure Functions'
541 CALL STRH3EM(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DGL)
545 ELSEIF(NSET.EQ.14) THEN
546 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
548 + ', for HMRS Set E Structure Functions'
549 CALL STRH3E(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DGL)
553 ELSEIF(NSET.EQ.15) THEN
554 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
556 + ', for HMRS Set B (190) Structure Functions'
557 CALL STRH3B(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DGL)
561 ELSEIF(NSET.EQ.16) THEN
562 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
564 + ', for HMRS Set B (100) Structure Functions'
565 CALL STRH3B1(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DGL)
569 ELSEIF(NSET.EQ.17) THEN
570 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
572 + ', for HMRS Set B (300) Structure Functions'
573 CALL STRH3B3(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DGL)
577 ELSEIF(NSET.EQ.18) THEN
578 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
580 + ', for KMRS Set B- Structure Functions'
581 CALL STRKBM(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DGL)
585 ELSEIF(NSET.EQ.19) THEN
586 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
588 + ', for KMRS Set B--R2-SH Structure Functions'
589 CALL SFKBMR2(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DGL)
593 ELSEIF(NSET.EQ.20) THEN
594 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
596 + ', for KMRS Set B--R5-SH Structure Functions'
597 CALL SFKBMR5(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DGL)
601 ELSEIF(NSET.EQ.21) THEN
602 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
604 + ', for KMRS Set B0 (190) updated Structure Functions'
605 CALL STRKB0(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
608 ELSEIF(NSET.EQ.22) THEN
609 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
611 + ', for MRS Set B0 (L135) Structure Functions'
612 CALL STRKB02(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DGL)
616 ELSEIF(NSET.EQ.23) THEN
617 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
619 + ', for MRS Set B0 (L160) Structure Functions'
620 CALL STRKB03(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DGL)
624 ELSEIF(NSET.EQ.24) THEN
625 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
627 + ', for MRS Set B0 (L200) Structure Functions'
628 CALL STRKB04(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DGL)
632 ELSEIF(NSET.EQ.25) THEN
633 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
635 + ', for MRS Set B0 (L235) Structure Functions'
636 CALL STRKB05(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DGL)
640 ELSEIF(NSET.EQ.26) THEN
641 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
643 + ', for MRS Set S0 (L215) Structure Functions'
644 IF (FIRST) WRITE(N6,*)
645 +' Warning !! : in principle, set has been retracted by the auth
647 CALL SFMRSS0(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
650 ELSEIF(NSET.EQ.27) THEN
651 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
653 + ', for MRS Set D0 (L215) Structure Functions'
654 IF (FIRST) WRITE(N6,*)
655 +' Warning !! : in principle, set has been retracted by the auth
657 CALL SFMRSD0(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
660 ELSEIF(NSET.EQ.28) THEN
661 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
663 + ', for MRS Set D- (L215) Structure Functions'
664 IF (FIRST) WRITE(N6,*)
665 +' Warning !! : in principle, set has been retracted by the auth
667 CALL SFMRSDM(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
670 ELSEIF(NSET.EQ.29) THEN
671 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
673 + ', for MRS Set S0p (L230-MSb) Structure Functions'
674 CALL SFMSBS0(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
677 ELSEIF(NSET.EQ.30) THEN
678 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
680 + ', for MRS Set D0p (L230-MSb) Structure Functions'
681 CALL SFMSBD0(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
684 ELSEIF(NSET.EQ.31) THEN
685 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
687 + ', for MRS Set D-p (L230-MSb) Structure Functions'
688 CALL SFMSBDM(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
691 ELSEIF(NSET.EQ.32) THEN
692 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
694 + ', for MRS Set S0p (L230-DIS) Structure Functions'
695 CALL SFDISS0(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
698 ELSEIF(NSET.EQ.33) THEN
699 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
701 + ', for MRS Set D0p (L230-DIS) Structure Functions'
702 CALL SFDISD0(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
705 ELSEIF(NSET.EQ.34) THEN
706 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
708 + ', for MRS Set D-p (L230-DIS) Structure Functions'
709 CALL SFDISDM(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
712 ELSEIF(NSET.EQ.35) THEN
713 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
715 + ', for MRS Set (H) (L230-MSb) Structure Functions'
716 CALL SFMSBH(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
719 ELSEIF(NSET.EQ.36) THEN
720 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
722 + ', for MRS Set (H) (L230-DIS) Structure Functions'
723 CALL SFDISH(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
726 ELSEIF(NSET.EQ.37) THEN
727 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
729 + ', for MRS Set (A) (L230-MSb) Structure Functions'
730 CALL SFMSBA(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
733 ELSEIF(NSET.EQ.38) THEN
734 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
736 + ', for MRS Fit (A) (L230-MSb) Structure Functions'
737 CALL SFMSFA(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
740 ELSEIF(NSET.EQ.39) THEN
741 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
743 + ', for MRS Set (Ap) (L231-MSb) Structure Functions'
744 CALL SFMSAP(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
747 ELSEIF(NSET.EQ.40) THEN
748 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
750 + ', for MRS Fit (Ap) (L231-MSb) Structure Functions'
751 CALL SFMFAP(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
754 ELSEIF(NSET.EQ.41) THEN
755 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
757 + ', for MRS Set (G) (L255-MSb) Structure Functions'
758 CALL SFMRSG(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
761 ELSEIF(NSET.EQ.42) THEN
762 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
764 + ', for MRS Fit (G) (L255-MSb) Structure Functions'
765 CALL SFMSFG(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
768 ELSEIF(NSET.EQ.43) THEN
769 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
771 + ', for MRS Set (A) (L230-MSb) and low Q**2 Structure
774 + CALL SFMSBA(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
776 + CALL SFMSQA(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
779 ELSEIF(NSET.EQ.44) THEN
780 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
782 + ', for MRS Set (A) (L230-DIS) and low Q**2 Structure
785 + CALL SFDISA(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
787 + CALL SFDSQA(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
790 ELSEIF(NSET.EQ.45) THEN
791 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
793 + ', for MRS Set (Ap) (L150-MSb) Structure Functions'
794 CALL SFMS105(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
797 ELSEIF(NSET.EQ.46) THEN
798 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
800 + ', for MRS Set (Ap) (L201-MSb) Structure Functions'
801 CALL SFMS110(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
804 ELSEIF(NSET.EQ.47) THEN
805 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
807 + ', for MRS Set (Ap) (L266-MSb) Structure Functions'
808 CALL SFMS115(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
811 ELSEIF(NSET.EQ.48) THEN
812 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
814 + ', for MRS Set (Ap) (L344-MSb) Structure Functions'
815 CALL SFMS120(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
818 ELSEIF(NSET.EQ.49) THEN
819 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
821 + ', for MRS Set (Ap) (L435-MSb) Structure Functions'
822 CALL SFMS125(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
825 ELSEIF(NSET.EQ.50) THEN
826 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
828 + ', for MRS Set (Ap) (L542-MSb) Structure Functions'
829 CALL SFMS130(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
832 ELSEIF(NSET.EQ.51) THEN
833 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
835 + ', for MRS Set (J) (L344-MSb) Structure Functions'
836 CALL SFMRSJ(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
839 ELSEIF(NSET.EQ.52) THEN
840 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
842 + ', for MRS Set (Jp) (L507-MSb) Structure Functions'
843 CALL SFMSJP(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
846 ELSEIF(NSET.EQ.53) THEN
847 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
849 + ', for MRS Set (R1) (L241-MSb) Structure Functions'
850 CALL SFMSR1(DX,DQ2,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
853 ELSEIF(NSET.EQ.54) THEN
854 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
856 + ', for MRS Set (R2) (L344-MSb) Structure Functions'
857 CALL SFMSR2(DX,DQ2,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
860 ELSEIF(NSET.EQ.55) THEN
861 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
863 + ', for MRS Set (R3) (L241-MSb) Structure Functions'
864 CALL SFMSR3(DX,DQ2,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
867 ELSEIF(NSET.EQ.56) THEN
868 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
870 + ', for MRS Set (R4) (L344-MSb) Structure Functions'
871 CALL SFMSR4(DX,DQ2,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
874 ELSEIF(NSET.EQ.57) THEN
875 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
877 + ', for MRS Fit (R1) (L241-MSb) Structure Functions'
878 CALL SFMFR1(DX,DQ2,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
881 ELSEIF(NSET.EQ.58) THEN
882 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
884 + ', for MRS Fit (R2) (L344-MSb) Structure Functions'
885 CALL SFMFR2(DX,DQ2,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
888 ELSEIF(NSET.GE.59) THEN
889 WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
891 + ', Structure Functions not yet exsistent'
894 ELSEIF(NGROUP.EQ.4) THEN
896 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
898 + ', for MT Set S1 (DIS) Structure Functions'
899 IF (FIRST) WRITE(N6,*)
900 +' Warning !! : in principle, set has been retracted by the auth
902 CALL SFMTU1(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
904 ELSEIF(NSET.EQ.2) THEN
905 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
907 + ', for MT Set B1 (DIS) Structure Functions'
908 IF (FIRST) WRITE(N6,*)
909 +' Warning !! : in principle, set has been retracted by the auth
911 CALL SFMTU2(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
913 ELSEIF(NSET.EQ.3) THEN
914 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
916 + ', for MT Set B2 (DIS) Structure Functions'
917 IF (FIRST) WRITE(N6,*)
918 +' Warning !! : in principle, set has been retracted by the auth
920 CALL SFMTU3(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
922 ELSEIF(NSET.EQ.4) THEN
923 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
925 + ', for MT Set E1 (DIS) Structure Functions'
926 IF (FIRST) WRITE(N6,*)
927 +' Warning !! : in principle, set has been retracted by the auth
929 CALL SFMTU4(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
931 ELSEIF(NSET.EQ.5) THEN
932 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
934 + ', for MT Set 6 (1/2s) (DIS) Structure Functions'
935 IF (FIRST) WRITE(N6,*)
936 +' Warning !! : in principle, set has been retracted by the auth
938 CALL SFMTU5(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
940 ELSEIF(NSET.EQ.6) THEN
941 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
943 + ', for MT Set S1 (MSb) Structure Functions'
944 IF (FIRST) WRITE(N6,*)
945 +' Warning !! : in principle, set has been retracted by the auth
947 CALL SFMTU6(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
949 ELSEIF(NSET.EQ.7) THEN
950 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
952 + ', for MT Set B1 (MSb) Structure Functions'
953 IF (FIRST) WRITE(N6,*)
954 +' Warning !! : in principle, set has been retracted by the auth
956 CALL SFMTU7(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
958 ELSEIF(NSET.EQ.8) THEN
959 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
961 + ', for MT Set B2 (MSb) Structure Functions'
962 IF (FIRST) WRITE(N6,*)
963 +' Warning !! : in principle, set has been retracted by the auth
965 CALL SFMTU8(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
967 ELSEIF(NSET.EQ.9) THEN
968 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
970 + ', for MT Set E1 (MSb) Structure Functions'
971 IF (FIRST) WRITE(N6,*)
972 +' Warning !! : in principle, set has been retracted by the auth
974 CALL SFMTU9(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
976 ELSEIF(NSET.EQ.10) THEN
977 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
979 + ', for MT Set 6 (1/2s) (MSb) Structure Functions'
980 IF (FIRST) WRITE(N6,*)
981 +' Warning !! : in principle, set has been retracted by the auth
983 CALL SFMT10(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
985 ELSEIF(NSET.EQ.11) THEN
986 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
988 + ', for MT Set LO Structure Functions'
989 IF (FIRST) WRITE(N6,*)
990 +' Warning !! : in principle, set has been retracted by the auth
992 CALL SFMT11(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
994 ELSEIF(NSET.EQ.12) THEN
995 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
997 + ', for CTEQ Set 1L (LO) Structure Functions'
998 IF (FIRST) WRITE(N6,*)
999 +' Warning !! : in principle, set has been retracted by the auth
1001 CALL SFCTQ11(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
1003 ELSEIF(NSET.EQ.13) THEN
1004 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1006 + ', for CTEQ Set 1M (MSb) Structure Functions'
1007 IF (FIRST) WRITE(N6,*)
1008 +' Warning !! : in principle, set has been retracted by the auth
1010 CALL SFCTQ12(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
1012 ELSEIF(NSET.EQ.14) THEN
1013 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1015 + ', for CTEQ Set 1MS (MSb) Structure Functions'
1016 IF (FIRST) WRITE(N6,*)
1017 +' Warning !! : in principle, set has been retracted by the auth
1019 CALL SFCTQ13(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
1021 ELSEIF(NSET.EQ.15) THEN
1022 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1024 + ', for CTEQ Set 1ML (MSb) Structure Functions'
1025 IF (FIRST) WRITE(N6,*)
1026 +' Warning !! : in principle, set has been retracted by the auth
1028 CALL SFCTQ14(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
1030 ELSEIF(NSET.EQ.16) THEN
1031 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1033 + ', for CTEQ Set 1D (DIS) Structure Functions'
1034 IF (FIRST) WRITE(N6,*)
1035 +' Warning !! : in principle, set has been retracted by the auth
1037 CALL SFCTQ15(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
1039 ELSEIF(NSET.EQ.17) THEN
1040 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1042 + ', for CTEQ Set 2L (LO) Structure Functions'
1043 IF (FIRST) WRITE(N6,*)
1044 +' Warning !! : in principle, set has been retracted by the auth
1046 CALL SFCTQ21(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
1048 ELSEIF(NSET.EQ.18) THEN
1049 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1051 + ', for CTEQ Set 2M (MSb) Structure Functions'
1052 IF (FIRST) WRITE(N6,*)
1053 +' Warning !! : in principle, set has been retracted by the auth
1055 CALL SFCTQ22(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
1057 ELSEIF(NSET.EQ.19) THEN
1058 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1060 + ', for CTEQ Set 2MS (MSb) Structure Functions'
1061 IF (FIRST) WRITE(N6,*)
1062 +' Warning !! : in principle, set has been retracted by the auth
1064 CALL SFCTQ23(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
1066 ELSEIF(NSET.EQ.20) THEN
1067 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1069 + ', for CTEQ Set 2MF (MSb) Structure Functions'
1070 IF (FIRST) WRITE(N6,*)
1071 +' Warning !! : in principle, set has been retracted by the auth
1073 CALL SFCTQ24(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
1075 ELSEIF(NSET.EQ.21) THEN
1076 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1078 + ', for CTEQ Set 2ML (MSb) Structure Functions'
1079 IF (FIRST) WRITE(N6,*)
1080 +' Warning !! : in principle, set has been retracted by the auth
1082 CALL SFCTQ25(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
1084 ELSEIF(NSET.EQ.22) THEN
1085 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1087 + ', for CTEQ Set 2D (DIS) Structure Functions'
1088 IF (FIRST) WRITE(N6,*)
1089 +' Warning !! : in principle, set has been retracted by the auth
1091 CALL SFCTQ26(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
1093 ELSEIF(NSET.EQ.23) THEN
1094 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1096 + ', for CTEQ Set 2pL (LO) Structure Functions'
1097 CALL SRCTQ21(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
1099 ELSEIF(NSET.EQ.24) THEN
1100 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1102 + ', for CTEQ Set 2pM (MSb) Structure Functions'
1103 CALL SRCTQ22(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
1105 ELSEIF(NSET.EQ.25) THEN
1106 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1108 + ', for CTEQ Set 2pMS (MSb) Structure Functions'
1109 CALL SRCTQ23(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
1111 ELSEIF(NSET.EQ.26) THEN
1112 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1114 + ', for CTEQ Set 2pMF (MSb) Structure Functions'
1115 CALL SRCTQ24(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
1117 ELSEIF(NSET.EQ.27) THEN
1118 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1120 + ', for CTEQ Set 2pML (MSb) Structure Functions'
1121 CALL SRCTQ25(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
1123 ELSEIF(NSET.EQ.28) THEN
1124 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1126 + ', for CTEQ Set 2pD (DIS) Structure Functions'
1127 CALL SRCTQ26(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
1129 ELSEIF(NSET.EQ.29) THEN
1130 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1132 + ', for CTEQ Set 3L (LO) Structure Functions'
1133 CALL SFCTQ31(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
1135 ELSEIF(NSET.EQ.30) THEN
1136 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1138 + ', for CTEQ Set 3M (MSb) Structure Functions'
1139 CALL SFCTQ32(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
1141 ELSEIF(NSET.EQ.31) THEN
1142 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1144 + ', for CTEQ Set 3D (DIS) Structure Functions'
1145 CALL SFCTQ33(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
1147 ELSEIF(NSET.EQ.32) THEN
1148 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1150 + ', for CTEQ Set 4L (LO) Structure Functions'
1151 CALL SFCTQ41(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
1153 ELSEIF(NSET.EQ.33) THEN
1154 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1156 + ', for CTEQ Set 4D (DIS) Structure Functions'
1157 CALL SFCTQ42(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
1159 ELSEIF(NSET.EQ.34) THEN
1160 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1162 + ', for CTEQ Set 4M (MSb) Structure Functions'
1163 CALL SFCTQ43(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
1165 ELSEIF(NSET.EQ.35) THEN
1166 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1168 + ', for CTEQ Set 4A1 (MSb) Structure Functions'
1169 CALL SFCTQ44(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
1171 ELSEIF(NSET.EQ.36) THEN
1172 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1174 + ', for CTEQ Set 4A2 (MSb) Structure Functions'
1175 CALL SFCTQ45(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
1177 ELSEIF(NSET.EQ.37) THEN
1178 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1180 + ', for CTEQ Set 4A3 (MSb) Structure Functions'
1181 CALL SFCTQ43(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
1183 ELSEIF(NSET.EQ.38) THEN
1184 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1186 + ', for CTEQ Set 4A4 (MSb) Structure Functions'
1187 CALL SFCTQ46(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
1189 ELSEIF(NSET.EQ.39) THEN
1190 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1192 + ', for CTEQ Set 4A5 (MSb) Structure Functions'
1193 CALL SFCTQ47(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
1195 ELSEIF(NSET.EQ.40) THEN
1196 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1198 + ', for CTEQ Set 4HJ (MSb) Structure Functions'
1199 CALL SFCTQ48(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
1201 ELSEIF(NSET.EQ.41) THEN
1202 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1204 + ', for CTEQ Set 4LQ (MSb) Structure Functions'
1205 CALL SFCTQ49(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
1207 ELSEIF(NSET.GE.42) THEN
1208 WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1210 + ', Structure Functions not yet exsistent'
1213 ELSEIF(NGROUP.EQ.5) THEN
1215 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1217 + ', for old GRV Set HO Structure Functions'
1218 IF (FIRST) WRITE(N6,*)
1219 +' Warning !! : in principle, set has been retracted by the auth
1221 CALL SFGRVH(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DTOP,DGL)
1224 ELSEIF(NSET.EQ.2) THEN
1225 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1227 + ', for old GRV Set LO Structure Functions'
1228 IF (FIRST) WRITE(N6,*)
1229 +' Warning !! : in principle, set has been retracted by the auth
1231 CALL SFGRVL(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DTOP,DGL)
1234 ELSEIF(NSET.EQ.3) THEN
1235 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1237 + ', for GRV Set HO Structure Functions'
1238 CALL GRVHO(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DTOP,DGL)
1241 ELSEIF(NSET.EQ.4) THEN
1242 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1244 + ', for GRV Set LO Structure Functions'
1245 CALL GRVLO(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DTOP,DGL)
1248 ELSEIF(NSET.EQ.5) THEN
1249 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1251 + ', for GRV-94 Set LO Structure Functions'
1252 CALL GRVLO(DX,DQ,ZUPV,ZDNV,ZUSEA,ZSTR,DCHM,DBOT,DTOP,ZGL)
1254 + DUPV,DDNV,DUSEA,DDSEA,DSTR,ZCHM,ZBOT,ZTOP,DGL)
1256 ELSEIF(NSET.EQ.6) THEN
1257 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1259 + ', for GRV-94 Set HO (MSb) Structure Functions'
1260 CALL GRVHO(DX,DQ,ZUPV,ZDNV,ZUSEA,ZSTR,DCHM,DBOT,DTOP,ZGL)
1262 + DUPV,DDNV,DUSEA,DDSEA,DSTR,ZCHM,ZBOT,ZTOP,DGL)
1264 ELSEIF(NSET.EQ.7) THEN
1265 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1267 + ', for GRV-94 Set HO (DIS) Structure Functions'
1268 CALL GRVHO(DX,DQ,ZUPV,ZDNV,ZUSEA,ZSTR,DCHM,DBOT,DTOP,ZGL)
1270 + DUPV,DDNV,DUSEA,DDSEA,DSTR,ZCHM,ZBOT,ZTOP,DGL)
1272 ELSEIF(NSET.EQ.8) THEN
1273 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1275 + ', for GRSV Set LO polarized Structure Functions (standard)'
1276 CALL GRSVL1(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DTOP,DGL)
1279 ELSEIF(NSET.EQ.9) THEN
1280 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1282 + ', for GRSV Set LO polarized Structure Functions (valence)'
1283 CALL GRSVL2(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DTOP,DGL)
1286 ELSEIF(NSET.EQ.10) THEN
1287 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1289 + ', for GRSV Set NLO polarized Structure Functions (standard)'
1290 CALL GRSVH1(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DTOP,DGL)
1293 ELSEIF(NSET.EQ.11) THEN
1294 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1296 + ', for GRSV Set NLO polarized Structure Functions (valence)'
1297 CALL GRSVH2(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DTOP,DGL)
1300 ELSEIF(NSET.GE.12) THEN
1301 WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1303 + ', Structure Functions not yet exsistent'
1306 ELSEIF(NGROUP.EQ.6) THEN
1308 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1310 + ', for ABFOW Structure Functions'
1311 CALL SFABFOW(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DGL)
1315 ELSEIF(NSET.GE.2) THEN
1316 WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1318 + ', Structure Functions not yet exsistent'
1321 ELSEIF(NGROUP.EQ.7) THEN
1323 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1325 + ', for BM Set A Structure Functions'
1326 CALL SFBMA(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
1329 ELSEIF(NSET.EQ.2) THEN
1330 IF (FIRST) WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1332 + ', for BM Set B Structure Functions'
1333 CALL SFBMB(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
1336 ELSEIF(NSET.GE.3) THEN
1337 WRITE(N6,*) ' Nucleon PDFs : Ngroup = ',NGROUP,
1339 + ', Structure Functions not yet exsistent'
1343 WRITE(N6,*)' NUCLEON Structure function call for unknown GROUP '
1344 + ,NGROUP,', and unknown NSET ',NSET
1348 C Now start PION structure functions
1349 IF(NPTYPE.EQ.2) THEN
1350 IF(NGROUP.EQ.1) THEN
1352 IF (FIRST) WRITE(N6,*) ' Pion PDFs : Ngroup = ',NGROUP,
1354 + ', for OW-P Set 1 Structure Functions'
1355 CALL STROWP1(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DGL)
1360 ELSEIF(NSET.EQ.2) THEN
1361 IF (FIRST) WRITE(N6,*) ' Pion PDFs : Ngroup = ',NGROUP,
1363 + ', for OW-P Set 2 Structure Functions'
1364 CALL STROWP2(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DGL)
1369 ELSEIF(NSET.GE.3) THEN
1370 WRITE(N6,*) ' Pion PDFs : Ngroup = ',NGROUP,
1372 + ', Structure Functions not yet exsistent'
1375 ELSEIF(NGROUP.EQ.2) THEN
1376 WRITE(N6,*) ' Pion PDFs : Ngroup = ',NGROUP,
1378 + ', Structure Functions not yet exsistent'
1380 ELSEIF(NGROUP.EQ.3) THEN
1382 IF (FIRST) WRITE(N6,*) ' Pion PDFs : Ngroup = ',NGROUP,
1384 + ', for SMRS-P Set 1 Structure Functions'
1385 CALL SMRSP31(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DGL)
1389 ELSEIF(NSET.EQ.2) THEN
1390 IF (FIRST) WRITE(N6,*) ' Pion PDFs : Ngroup = ',NGROUP,
1392 + ', for SMRS-P Set 2 Structure Functions'
1393 CALL SMRSP32(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DGL)
1397 ELSEIF(NSET.EQ.3) THEN
1398 IF (FIRST) WRITE(N6,*) ' Pion PDFs : Ngroup = ',NGROUP,
1400 + ', for SMRS-P Set 3 Structure Functions'
1401 CALL SMRSP33(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DGL)
1405 ELSEIF(NSET.GE.4) THEN
1406 WRITE(N6,*) ' Pion PDFs : Ngroup = ',NGROUP,
1408 + ', Structure Functions not yet exsistent'
1411 ELSEIF(NGROUP.EQ.4) THEN
1412 WRITE(N6,*) ' Pion PDFs : Ngroup = ',NGROUP,
1414 + ', Structure Functions not yet exsistent'
1416 ELSEIF(NGROUP.EQ.5) THEN
1418 IF (FIRST) WRITE(N6,*) ' Pion PDFs : Ngroup = ',NGROUP,
1420 + ', for GRV-P Set HO Structure Functions'
1421 CALL GRVPIHO(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DTOP,DGL)
1424 ELSEIF(NSET.EQ.2) THEN
1425 IF (FIRST) WRITE(N6,*) ' Pion PDFs : Ngroup = ',NGROUP,
1427 + ', for GRV-P Set LO Structure Functions'
1428 CALL GRVPILO(DX,DQ,DUPV,DDNV,DUSEA,DSTR,DCHM,DBOT,DTOP,DGL)
1431 ELSEIF(NSET.GE.3) THEN
1432 WRITE(N6,*) ' Pion PDFs : Ngroup = ',NGROUP,
1434 + ', Structure Functions not yet exsistent'
1437 ELSEIF(NGROUP.EQ.6) THEN
1439 IF (FIRST) WRITE(N6,*) ' Pion PDFs : Ngroup = ',NGROUP,
1441 + ', for ABFKW-P Set 1 Structure Functions'
1442 CALL ABFKW1(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DGL)
1446 ELSEIF(NSET.EQ.2) THEN
1447 IF (FIRST) WRITE(N6,*) ' Pion PDFs : Ngroup = ',NGROUP,
1449 + ', for ABFKW-P Set 2 Structure Functions'
1450 CALL ABFKW2(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DGL)
1454 ELSEIF(NSET.EQ.3) THEN
1455 IF (FIRST) WRITE(N6,*) ' Pion PDFs : Ngroup = ',NGROUP,
1457 + ', for ABFKW-P Set 3 Structure Functions'
1458 CALL ABFKW3(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DGL)
1462 ELSEIF(NSET.GE.4) THEN
1463 WRITE(N6,*) ' Pion PDFs : Ngroup = ',NGROUP,
1465 + ', Structure Functions not yet exsistent'
1469 WRITE(N6,*) ' PION Structure function call for unknown GROUP ',
1470 + NGROUP,', and unknown NSET ',NSET
1474 C Now start PHOTON structure functions
1475 IF(NPTYPE.EQ.3) THEN
1476 IF(NGROUP.EQ.1) THEN
1478 IF (FIRST) WRITE(N6,*) ' Photon PDFs : Ngroup = ',NGROUP,
1480 + ', for DO-G Set LO Structure Functions'
1481 CALL DOPHO1(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
1484 ELSEIF(NSET.EQ.2) THEN
1485 IF (FIRST) WRITE(N6,*) ' Photon PDFs : Ngroup = ',NGROUP,
1487 + ', for DO-G Set NLL Structure Functions'
1488 CALL DOPHO2(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
1491 ELSEIF(NSET.GE.3) THEN
1492 WRITE(N6,*) ' Photon PDFs : Ngroup = ',NGROUP,
1494 + ', Structure Functions not yet exsistent'
1497 ELSEIF(NGROUP.EQ.2) THEN
1499 IF (FIRST) WRITE(N6,*) ' Photon PDFs : Ngroup = ',NGROUP,
1501 + ', for DG-G Set 1 Structure Functions'
1502 CALL DGPHO1(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
1505 ELSEIF(NSET.EQ.2) THEN
1506 IF (FIRST) WRITE(N6,*) ' Photon PDFs : Ngroup = ',NGROUP,
1508 + ', for DG-G Set 2 Structure Functions'
1509 CALL DGPHO2(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
1512 ELSEIF(NSET.EQ.3) THEN
1513 IF (FIRST) WRITE(N6,*) ' Photon PDFs : Ngroup = ',NGROUP,
1515 + ', for DG-G Set 3 Structure Functions'
1516 CALL DGPHO3(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
1520 ELSEIF(NSET.EQ.4) THEN
1521 IF (FIRST) WRITE(N6,*) ' Photon PDFs : Ngroup = ',NGROUP,
1523 + ', for DG-G Set 4 Structure Functions'
1524 CALL DGPHO4(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
1528 ELSEIF(NSET.GE.5) THEN
1529 WRITE(N6,*) ' Photon PDFs : Ngroup = ',NGROUP,
1531 + ', Structure Functions not yet exsistent'
1534 ELSEIF(NGROUP.EQ.3) THEN
1536 IF (FIRST) WRITE(N6,*) ' Photon PDFs : Ngroup = ',NGROUP,
1538 + ', for LAC-G Set 1 Structure Functions'
1539 CALL SFLACG1(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
1541 ELSEIF(NSET.EQ.2) THEN
1542 IF (FIRST) WRITE(N6,*) ' Photon PDFs : Ngroup = ',NGROUP,
1544 + ', for LAC-G Set 2 Structure Functions'
1545 CALL SFLACG2(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
1547 ELSEIF(NSET.EQ.3) THEN
1548 IF (FIRST) WRITE(N6,*) ' Photon PDFs : Ngroup = ',NGROUP,
1550 + ', for LAC-G Set 3 Structure Functions'
1551 CALL SFLACG3(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
1553 ELSEIF(NSET.GE.4) THEN
1554 WRITE(N6,*) ' Photon PDFs : Ngroup = ',NGROUP,
1556 + ', Structure Functions not yet exsistent'
1559 ELSEIF(NGROUP.EQ.4) THEN
1561 IF (FIRST) WRITE(N6,*) ' Photon PDFs : Ngroup = ',NGROUP,
1563 + ', for GS-G Set HO Structure Functions'
1564 CALL SFGSHO(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
1567 ELSEIF(NSET.EQ.2) THEN
1568 IF (FIRST) WRITE(N6,*) ' Photon PDFs : Ngroup = ',NGROUP,
1570 + ', for GS-G LO Set 1 Structure Functions'
1571 CALL SFGSLO1(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
1574 ELSEIF(NSET.EQ.3) THEN
1575 IF (FIRST) WRITE(N6,*) ' Photon PDFs : Ngroup = ',NGROUP,
1577 + ', for GS-G LO Set 2 Structure Functions'
1578 CALL SFGSLO2(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
1581 ELSEIF(NSET.EQ.4) THEN
1582 IF (FIRST) WRITE(N6,*) ' Photon PDFs : Ngroup = ',NGROUP,
1584 + ', for GS-96-G HO Structure Functions'
1585 CALL GS96HO(DX,DQ,DUPV,DDNV,DSTR,DCHM,DBOT,DGL)
1590 ELSEIF(NSET.EQ.5) THEN
1591 IF (FIRST) WRITE(N6,*) ' Photon PDFs : Ngroup = ',NGROUP,
1593 + ', for GS-96-G LO Structure Functions'
1594 CALL GS96LO(DX,DQ,DUPV,DDNV,DSTR,DCHM,DBOT,DGL)
1599 ELSEIF(NSET.GE.6) THEN
1600 WRITE(N6,*) ' Photon PDFs : Ngroup = ',NGROUP,
1602 + ', Structure Functions not yet exsistent'
1605 ELSEIF(NGROUP.EQ.5) THEN
1607 IF (FIRST) WRITE(N6,*) ' Photon PDFs : Ngroup = ',NGROUP,
1609 + ', for GRV-G Set L-HO Structure Functions'
1610 CALL GRVGAH0(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
1613 ELSEIF(NSET.EQ.2) THEN
1614 IF (FIRST) WRITE(N6,*) ' Photon PDFs : Ngroup = ',NGROUP,
1616 + ', for GRV-G Set HO Structure Functions'
1617 CALL GRVGAHO(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
1620 ELSEIF(NSET.EQ.3) THEN
1621 IF (FIRST) WRITE(N6,*) ' Photon PDFs : Ngroup = ',NGROUP,
1623 + ', for GRV-G Set LO Structure Functions'
1624 CALL GRVGALO(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
1627 ELSEIF(NSET.EQ.4) THEN
1629 WRITE(N6,*) ' Photon PDFs : Ngroup = ',NGROUP,
1631 + ', for GRS-G Set LO Photon Structure Functions'
1632 WRITE(N6,*) ' !! You MUST call STRUCTP '
1633 WRITE(N6,*) ' !! to get the INPUT parameters right !!!'
1636 ELSEIF(NSET.GE.5) THEN
1637 WRITE(N6,*) ' Photon PDFs : Ngroup = ',NGROUP,
1639 + ', Structure Functions not yet exsistent'
1642 ELSEIF(NGROUP.EQ.6) THEN
1644 IF (FIRST) WRITE(N6,*) ' Photon PDFs : Ngroup = ',NGROUP,
1646 + ', for ACFGP-G Set HO Structure Functions'
1647 CALL ACFGP1(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DGL)
1651 ELSEIF(NSET.EQ.2) THEN
1652 IF (FIRST) WRITE(N6,*) ' Photon PDFs : Ngroup = ',NGROUP,
1654 + ', for ACFGP-G Set HO-mc Structure Functions'
1655 CALL ACFGP2(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DGL)
1659 ELSEIF(NSET.EQ.3) THEN
1660 IF (FIRST) WRITE(N6,*) ' Photon PDFs : Ngroup = ',NGROUP,
1662 + ', for AFG-G Set HO Structure Functions'
1663 CALL SFAFG1(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DGL)
1667 ELSEIF(NSET.GE.4) THEN
1668 WRITE(N6,*) ' Photon PDFs : Ngroup = ',NGROUP,
1670 + ', Structure Functions not yet exsistent'
1673 ELSEIF(NGROUP.EQ.7) THEN
1674 WRITE(N6,*) ' Photon PDFs : Ngroup = ',NGROUP,
1676 + ', Structure Functions not yet exsistent'
1678 ELSEIF(NGROUP.EQ.8) THEN
1680 IF (FIRST) WRITE(N6,*) ' Photon PDFs : Ngroup = ',NGROUP,
1682 + ', for WHIT-G Set 1 Structure Functions'
1683 CALL SFWHI1(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DGL)
1687 ELSEIF(NSET.EQ.2) THEN
1688 IF (FIRST) WRITE(N6,*) ' Photon PDFs : Ngroup = ',NGROUP,
1690 + ', for WHIT-G Set 2 Structure Functions'
1691 CALL SFWHI2(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DGL)
1695 ELSEIF(NSET.EQ.3) THEN
1696 IF (FIRST) WRITE(N6,*) ' Photon PDFs : Ngroup = ',NGROUP,
1698 + ', for WHIT-G Set 3 Structure Functions'
1699 CALL SFWHI3(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DGL)
1703 ELSEIF(NSET.EQ.4) THEN
1704 IF (FIRST) WRITE(N6,*) ' Photon PDFs : Ngroup = ',NGROUP,
1706 + ', for WHIT-G Set 4 Structure Functions'
1707 CALL SFWHI4(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DGL)
1711 ELSEIF(NSET.EQ.5) THEN
1712 IF (FIRST) WRITE(N6,*) ' Photon PDFs : Ngroup = ',NGROUP,
1714 + ', for WHIT-G Set 5 Structure Functions'
1715 CALL SFWHI5(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DGL)
1719 ELSEIF(NSET.EQ.6) THEN
1720 IF (FIRST) WRITE(N6,*) ' Photon PDFs : Ngroup = ',NGROUP,
1722 + ', for WHIT-G Set 6 Structure Functions'
1723 CALL SFWHI6(DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DGL)
1727 ELSEIF(NSET.GE.7) THEN
1728 WRITE(N6,*) ' Photon PDFs : Ngroup = ',NGROUP,
1730 + ', Structure Functions not yet exsistent'
1733 ELSEIF(NGROUP.EQ.9) THEN
1736 WRITE(N6,*) ' Photon PDFs : Ngroup = ',NGROUP,
1738 + ', for SAS Photon Structure Functions'
1739 WRITE(N6,*) ' !! You MUST call STRUCTP '
1740 WRITE(N6,*) ' !! to get the INPUT parameters right !!!'
1743 ELSEIF(NSET.GE.9) THEN
1744 WRITE(N6,*) ' Photon PDFs : Ngroup = ',NGROUP,
1746 + ', Structure Functions not yet exsistent'
1750 WRITE(N6,*)' PHOTON Structure function call for unknown GROUP ',
1751 + NGROUP,', and unknown NSET ',NSET
1755 100 IF(FIRST) WRITE(N6,*)
1756 +' ----------------------------------------------------------------
1757 +------------------------------------------------'
1758 IF (FIRST) FIRST = .FALSE.
1760 UPV = MAX( ZEROD,DUPV)
1761 DNV = MAX( ZEROD,DDNV)
1762 USEA = MAX( ZEROD,DUSEA)
1763 DSEA = MAX( ZEROD,DDSEA)
1764 STR = MAX( ZEROD,DSTR)
1765 CHM = MAX( ZEROD,DCHM)
1766 BOT = MAX( ZEROD,DBOT)
1767 TOP = MAX( ZEROD,DTOP)
1768 GL = MAX( ZEROD,DGL)
1772 200 IF(FIRST) WRITE(N6,*)
1773 +' ----------------------------------------------------------------
1774 +------------------------------------------------'
1775 IF (FIRST) FIRST = .FALSE.