CDECK  ID>, QQLMAT.

*CMZ :-        -28/01/92  12.34.44  by  Mike Seymour

*-- Author :    Luca Stanco

C-----------------------------------------------------------------------

      INTEGER FUNCTION QQLMAT(IDL,NDIR)

C-----------------------------------------------------------------------

C. QQLMAT - Given a particle flavor (KF), converts it to QQ particle number

C.          (KF = IDPDG code)

C.

C. Inputs    : IDL    (input  particle code)

C              NDIR = 1   LUND --> QQ

C              NDIR = 2   QQ   --> LUND

C

C. Outputs   : QQLMAT (output particle code)

C.

C-----------------------------------------------------------------------

      IMPLICIT NONE

C-- Calling variable

      INTEGER IDL,NDIR

C-- External declaration

C-- Local variables

      INTEGER AKF(321),I

      DATA (AKF(I), I=1,151) /

     +    0,    0,    0,    0,    0,    0,    0,   21,   -6,   -5,

     +   -4,   -3,   -1,   -2,    6,    5,    4,    3,    1,    2,

     +    0,

     +   22,   23,   24,  -24,   90,    0,   11,  -11,   12,  -12,

     +   13,  -13,   14,  -14,   15,  -15,   16,  -16,20313,-20313,

     +  211, -211,  321, -321,  311, -311,  421, -421,  411, -411,

     +  431, -431, -521,  521, -511,  511, -531,  531, -541,  541,

     +  621, -621,  611, -611,  631, -631,  641, -641,  651, -651,

     +  111,  221,  331,  441,20551,  661,  310,  130,10313,-10313,

     +  213, -213,  323, -323,  313, -313,  423, -423,  413, -413,

     +  433, -433, -523,  523, -513,  513, -533,  533, -543,  543,

     +  623, -623,  613, -613,  633, -633,  643, -643,  653, -653,

     +  113,  223,  333,  443,  553,  136,  20553, 30553, 40553, 551,

     +  10553, 555, 10551,70553,10555, 0, 20213, 20113, -20213, 10441,

     +  10443, 445, 8*0,

     +  3122, -3122, 4122, -4122, 4232, -4232, 4132, -4132, 3212, -3212/

      DATA (AKF(I), I=152,321) /

     +  4212, -4212, 4322, -4322, 4312, -4312, 2212, -2212, 3222, -3222,

     +  4222, -4222, 2112, -2112, 3112, -3112, 4112, -4112, 3322, -3322,

     +  3312, -3312, 4332, -4332, 6*0,

     +  3214, -3214, 4214, -4214, 4324, -4324, 4314, -4314, 2214, -2214,

     +  3224, -3224, 4224, -4224, 2114, -2114, 3114, -3114, 4114, -4114,

     +  3324, -3324, 3314, -3314, 4334, -4334, 4*0,

     +  0, 0,  2224, -2224, 1114, -1114, 3334, -3334, 0, 0,

     +  10323, -10323, 20323, -20323, 6*0,

     +  30443, 0, 0, 0, 70443, 50553, 60553, 80553, 20443, 0,

     +  10411, 20413, 10413, 415,

     + -10411,-20413,-10413,-415,

     +  10421, 20423, 10423, 425,

     + -10421,-20423,-10423,-425,

     +  10431, 20433, 10433, 435,

     + -10431,-20433,-10433,-435, 0,0,0,0,0,0,

     +  10111, 10211,-10211, 115, 215, -215,10221,10331,20223,20333,

     +  225, 335, 10223, 10333, 10113, 10213,-10213, 33*0 /

      IF(NDIR.EQ.1) THEN

        DO 10 I=1,321

        IF (IDL.EQ.AKF(I)) THEN

          QQLMAT=I-21

          RETURN

        ENDIF

  10    CONTINUE

        WRITE(6,20) IDL

  20    FORMAT(1X,'Lund code particle ',I6,' not recognized')

      ELSEIF(NDIR.EQ.2) THEN

        QQLMAT = AKF(IDL+21)

      ENDIF

      RETURN

      END