*
* $Id$
*
* $Log$
* Revision 1.1  2000/06/19 20:00:37  eugenio
* Initial revision
*
* Revision 1.2  1996/04/02 15:33:26  zfiles
* Updates from Peter.
* Mass/decay width/branching ratios of all particles are checked and
* updated to match PDG.
*
* Revision 1.1.1.1  1994/10/08  02:21:34  zfiles
* first version of qqlib in CVS
*
*
#include "sys/CLEO_machine.h"
#include "pilot.h"
*CMZ :  1.03/23 23/09/91  15.15.19  by  Peter C Kim
*CMZ :  1.03/22 19/09/91  12.10.04  by  Peter C Kim
*CMZ :  1.03/03 04/04/91  10.55.53  by  Peter C Kim
*CMZ :  1.03/00 02/04/91  13.56.07  by  Peter C Kim
*CMZ :  1.02/00 21/02/91  11.50.51  by  Peter C Kim
*-- Author :    Peter C Kim   12/02/91
      INTEGER FUNCTION LQPMAT(IDL,NDIR)
C.......................................................................
C.
C. LQPMAT - Given a particle flavor (KF), converts it to QQ particle number
C.
C. Inputs    : IDL    (input  particle code)
C              NDIR = 1   LUND --> QQ
C              NDIR = 2   QQ   --> LUND
C
C. Outputs   : LQPMAT (output particle code)
C.
C. Calls     : IUCOMP
C. Called    : DECADL
C.
C. Author    : Peter C Kim  12/02/91  11.12.54
C.......................................................................
#if defined(CLEO_TYPECHEK)
      IMPLICIT NONE
#endif

C-- Calling variable

      INTEGER IDL,NDIR

C-- External declaration

      INTEGER IUCOMP
      EXTERNAL IUCOMP

C-- Local variables

      INTEGER AKF(521),NLIST,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,  551,  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,
     +  621, -621,  613, -613,  633, -633,  643, -643,  653, -653,
     +  113,  223,  333,  443,  553,  663,  30553, 50553, 60553, 10551,
     +  20553, 555, 20551,40553,10555,40113,20213,20113,-20213,10441,
     +  20443, 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, 40213, -40213, 4*0,
     +  40443,50443,60443,70443,10443,70553,80553,10553,30443,90553,
     +  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, 10311,-10311,10321,
     +  -10321, 315,-315, 325,-325, 25*0 /

      DATA (AKF(I), I=322,521) / 200*0.0 /


      IF(NDIR.EQ.1) THEN

        NLIST = 521
        LQPMAT = IUCOMP(IDL, AKF, NLIST) - 21

      ELSEIF(NDIR.EQ.2) THEN

        LQPMAT = AKF(IDL+21)

      ENDIF

      RETURN
      END