This commit was generated by cvs2svn to compensate for changes in r2,

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

diff --git a/PDF/spdf/srctq24.F b/PDF/spdf/srctq24.F
new file mode 100644 (file)
index 0000000..e53fd98
--- /dev/null
+++ b/PDF/spdf/srctq24.F
@@ -0,0 +1,60 @@
+*
+* $Id$
+*
+* $Log$
+* Revision 1.1.1.2  1996/10/30 08:30:56  cernlib
+* Version 7.04
+*
+* Revision 1.1.1.1  1996/04/12 15:29:48  plothow
+* Version 7.01
+*
+*
+#include "pdf/pilot.h"
+C-----------------------------------------------------------------------
+      SUBROUTINE SRCTQ24(DX,DSCAL,
+     +                     DUPV,DDNV,DSEA,DSEAD,DSTR,DCHM,DBOT,DTOP,DGL)
+C
+C   ********************************************************************
+C   *                                                                  *
+C   *        Interface to CTEQset of structure functions               *
+C   *                                                                  *
+C   *        Author:    H. Plothow-Besch (CERN-PPE)                    *
+C   *                                                                  *
+C   ********************************************************************
+C
+C :::::::::::: Structure functions from Morfin and Wu-Ki Tung
+C :::::::::::: Lambda = 0.208 GeV (MSb) (flat small-x) (CTEQ2MF)
+C
+#include "pdf/expdp.inc"
+     +       DX,DSCAL,DUPV,DDNV,DSEA,DSEAD,DSTR,DCHM,DBOT,DTOP,DGL
+      REAL   X, Q
+      PARAMETER (ISET=3)
+#include "pdf/w50511.inc"
+C
+      X = DX
+      Q = DSCAL
+C
+C     generate the individual structure fcn calls
+C
+      UPV = Ctq2pPd(ISET, 1,X,Q,IFL)
+      DUPV = UPV
+      DNV = Ctq2pPd(ISET, 2,X,Q,IFL)
+      DDNV = DNV
+      SEAU = Ctq2pPd(ISET,-1,X,Q,IFL)
+      DSEA = SEAU
+      SEAD = Ctq2pPd(ISET,-2,X,Q,IFL)
+      DSEAD = SEAD
+      STR = Ctq2pPd(ISET,-3,X,Q,IFL)
+      DSTR = STR
+      CHM = Ctq2pPd(ISET,-4,X,Q,IFL)
+      DCHM = CHM
+      BOT = Ctq2pPd(ISET,-5,X,Q,IFL)
+      DBOT = BOT
+      TOP = 0.
+      IF (DSCAL.GT.TMAS) TOP = Ctq2pPd(ISET,-6,X,Q,IFL)
+      DTOP = TOP
+      GL = Ctq2pPd(ISET, 0,X,Q,IFL)
+      DGL = GL
+C
+      RETURN
+      END