X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;ds=sidebyside;f=PYTHIA6%2Fpdf_alice.F;h=ca1ca70bff2ed2c8ce9cb0585348af559bca746a;hb=190d081e2aff2c760ac35a8bbd4b66c1d454dc1a;hp=6e6e7546cb41b34b459689193fe169895d4303bb;hpb=2eef49c001b3c8ffa5c6cc35871f984ebe16a075;p=u%2Fmrichter%2FAliRoot.git diff --git a/PYTHIA6/pdf_alice.F b/PYTHIA6/pdf_alice.F index 6e6e7546cb4..ca1ca70bff2 100644 --- 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