* EVGEN/AliGenEMCocktail.cxx (& .h) is the class that generates the cocktail follow...

[u/mrichter/AliRoot.git] / PYTHIA6 / pdf_alice.F

diff --git a/PYTHIA6/pdf_alice.F b/PYTHIA6/pdf_alice.F
index 6e6e7546cb41b34b459689193fe169895d4303bb..ca1ca70bff2ed2c8ce9cb0585348af559bca746a 100644 (file)
--- a/PYTHIA6/pdf_alice.F
+++ b/PYTHIA6/pdf_alice.F
@@ -19,6 +19,15 @@ C...  MINT(124) : side (1 or 2)
 C...
       IMPLICIT DOUBLE PRECISION(A-H, O-Z)
       IMPLICIT INTEGER(I-N)
+C...Interface to PDFLIB.
+      COMMON/LW50512/QCDL4,QCDL5
+      SAVE /LW50512/
+      DOUBLE PRECISION QCDL4,QCDL5
+      COMMON/LW50513/XMIN,XMAX,Q2MIN,Q2MAX
+      SAVE /LW50513/
+      DOUBLE PRECISION XMIN,XMAX,Q2MIN,Q2MAX
+C...
+      COMMON/PYDAT1/MSTU(200),PARU(200),MSTJ(200),PARJ(200)
       COMMON/PYPARS/MSTP(200),PARP(200),MSTI(200),PARI(200)  
       DOUBLE PRECISION VALUE(20)
       CHARACTER*20 PARM(20)
@@ -33,10 +42,21 @@ C...
          VALUE(6)=MSTP(191)/1000
          PARM(7)='NASET'
          VALUE(7)=MOD(MSTP(191),1000)
-         CALL PDFSET(PARM,VALUE)
+         CALL PDFSET(PARM,VALUE,
+     >        MSTU(11),MSTP(51),MSTP(53),MSTP(55),
+     >        QCDL4,QCDL5,
+     >        XMIN,XMAX,Q2MIN,Q2MAX)
+         IF (MSTP(194) .EQ. 0) THEN 
+            CALL SETLHAPARM("EKS98")
+         ELSE
+            CALL SETLHAPARM("EPS08")
+         ENDIF
       ELSE 
          write(6,*) "-> pdfset"
-         CALL PDFSET(PARM,VALUE)
+         CALL PDFSET(PARM,VALUE,
+     >        MSTU(11),MSTP(51),MSTP(53),MSTP(55),
+     >        QCDL4,QCDL5,
+     >        XMIN,XMAX,Q2MIN,Q2MAX)
       ENDIF
       write(6,*) "done"
       END