[u/mrichter/AliRoot.git] / PDF / spdf / mrseb.F

*
* $Id$
*
* $Log$
* Revision 1.1.1.2  1996/10/30 08:28:53  cernlib
* Version 7.04
*
* Revision 1.1.1.1  1996/04/12 15:29:25  plothow
* Version 7.01
*
*
#include "pdf/pilot.h"
C-----------------------------------------------------------------------
      SUBROUTINE MRSEB(X,SCALE,MODE,UPV,DNV,SEA,STR,CHM,BOT,GL)
C***************************************************************C
C                                                               C
C  MODE 1 CORRESPONDS TO                                        C
C  MARTIN, ROBERTS, STIRLING (EMC FIT)    WITH LAMBDA= 100 MEV  C
C                                                               C
C  MODE 2  CORRESPONDS TO                                       C
C  MARTIN, ROBERTS, STIRLING (BCDMS FIT)  WITH LAMBDA= 200 MEV  C
C                                                               C
C  (  SOFT GLUE :  X G(X,Q0) = A (1-X)**5 )                     C
C                                                               C
C                         -*-                                   C
C                                                               C
C    (NOTE THAT X TIMES THE PARTON DISTRIBUTION FUNCTION        C
C    IS RETURNED I.E. G(X) = GLU/X ETC, AND THAT "SEA"          C
C    IS THE LIGHT QUARK SEA I.E. UBAR(X)=DBAR(X)=  ...          C
C    = SEA/X FOR A PROTON.  IF IN DOUBT, CHECK THE              C
C    MOMENTUM SUM RULE! NOTE ALSO THAT SCALE=Q IN GEV)          C
C                                                               C
C                         -*-                                   C
C                                                               C
C     (THE RANGE OF APPLICABILITY IS CURRENTLY:                 C
C     10**-4 < X < 1  AND  5 < Q**2 < 1.31 * 10**6              C
C     HIGHER Q**2 VALUES CAN BE SUPPLIED ON REQUEST             C
C     - PROBLEMS, COMMENTS ETC TO SRG$T3@GEN                    C
C                                                               C
C                                                               C
C***************************************************************C
#include "pdf/impdp.inc"
      IF(MODE.EQ.1) CALL STRUCE(X,SCALE,UPV,DNV,SEA,STR,CHM,BOT,GL)
      IF(MODE.EQ.2) CALL STRUCB(X,SCALE,UPV,DNV,SEA,STR,CHM,BOT,GL)
      RETURN
      END
Commit	Line	Data
21886bb6	1	*
	2	* $Id$
	3	*
	4	* $Log$
	5	* Revision 1.1.1.2 1996/10/30 08:28:53 cernlib
	6	* Version 7.04
	7	*
	8	* Revision 1.1.1.1 1996/04/12 15:29:25 plothow
	9	* Version 7.01
	10	*
	11	*
	12	#include "pdf/pilot.h"
	13	C-----------------------------------------------------------------------
	14	SUBROUTINE MRSEB(X,SCALE,MODE,UPV,DNV,SEA,STR,CHM,BOT,GL)
	15	C***************************************************************C
	16	C C
	17	C MODE 1 CORRESPONDS TO C
	18	C MARTIN, ROBERTS, STIRLING (EMC FIT) WITH LAMBDA= 100 MEV C
	19	C C
	20	C MODE 2 CORRESPONDS TO C
	21	C MARTIN, ROBERTS, STIRLING (BCDMS FIT) WITH LAMBDA= 200 MEV C
	22	C C
	23	C ( SOFT GLUE : X G(X,Q0) = A (1-X)**5 ) C
	24	C C
	25	C -*- C
	26	C C
	27	C (NOTE THAT X TIMES THE PARTON DISTRIBUTION FUNCTION C
	28	C IS RETURNED I.E. G(X) = GLU/X ETC, AND THAT "SEA" C
	29	C IS THE LIGHT QUARK SEA I.E. UBAR(X)=DBAR(X)= ... C
	30	C = SEA/X FOR A PROTON. IF IN DOUBT, CHECK THE C
	31	C MOMENTUM SUM RULE! NOTE ALSO THAT SCALE=Q IN GEV) C
	32	C C
	33	C -*- C
	34	C C
	35	C (THE RANGE OF APPLICABILITY IS CURRENTLY: C
	36	C 10-4 < X < 1 AND 5 < Q2 < 1.31 * 10**6 C
	37	C HIGHER Q**2 VALUES CAN BE SUPPLIED ON REQUEST C
	38	C - PROBLEMS, COMMENTS ETC TO SRG$T3@GEN C
	39	C C
	40	C C
	41	C***************************************************************C
	42	#include "pdf/impdp.inc"
	43	IF(MODE.EQ.1) CALL STRUCE(X,SCALE,UPV,DNV,SEA,STR,CHM,BOT,GL)
	44	IF(MODE.EQ.2) CALL STRUCB(X,SCALE,UPV,DNV,SEA,STR,CHM,BOT,GL)
	45	RETURN
	46	END