HERWIG/src/qqlmat.f

   1
   2 CDECK  ID>, QQLMAT.
   3
   4 *CMZ :-        -28/01/92  12.34.44  by  Mike Seymour
   5
   6 *-- Author :    Luca Stanco
   7
   8 C-----------------------------------------------------------------------
   9
  10       INTEGER FUNCTION QQLMAT(IDL,NDIR)
  11
  12 C-----------------------------------------------------------------------
  13
  14 C. QQLMAT - Given a particle flavor (KF), converts it to QQ particle number
  15
  16 C.          (KF = IDPDG code)
  17
  18 C.
  19
  20 C. Inputs    : IDL    (input  particle code)
  21
  22 C              NDIR = 1   LUND --> QQ
  23
  24 C              NDIR = 2   QQ   --> LUND
  25
  26 C
  27
  28 C. Outputs   : QQLMAT (output particle code)
  29
  30 C.
  31
  32 C-----------------------------------------------------------------------
  33
  34       IMPLICIT NONE
  35
  36 C-- Calling variable
  37
  38       INTEGER IDL,NDIR
  39
  40 C-- External declaration
  41
  42 C-- Local variables
  43
  44       INTEGER AKF(321),I
  45
  46       DATA (AKF(I), I=1,151) /
  47
  48      +    0,    0,    0,    0,    0,    0,    0,   21,   -6,   -5,
  49
  50      +   -4,   -3,   -1,   -2,    6,    5,    4,    3,    1,    2,
  51
  52      +    0,
  53
  54      +   22,   23,   24,  -24,   90,    0,   11,  -11,   12,  -12,
  55
  56      +   13,  -13,   14,  -14,   15,  -15,   16,  -16,20313,-20313,
  57
  58      +  211, -211,  321, -321,  311, -311,  421, -421,  411, -411,
  59
  60      +  431, -431, -521,  521, -511,  511, -531,  531, -541,  541,
  61
  62      +  621, -621,  611, -611,  631, -631,  641, -641,  651, -651,
  63
  64      +  111,  221,  331,  441,20551,  661,  310,  130,10313,-10313,
  65
  66      +  213, -213,  323, -323,  313, -313,  423, -423,  413, -413,
  67
  68      +  433, -433, -523,  523, -513,  513, -533,  533, -543,  543,
  69
  70      +  623, -623,  613, -613,  633, -633,  643, -643,  653, -653,
  71
  72      +  113,  223,  333,  443,  553,  136,  20553, 30553, 40553, 551,
  73
  74      +  10553, 555, 10551,70553,10555, 0, 20213, 20113, -20213, 10441,
  75
  76      +  10443, 445, 8*0,
  77
  78      +  3122, -3122, 4122, -4122, 4232, -4232, 4132, -4132, 3212, -3212/
  79
  80       DATA (AKF(I), I=152,321) /
  81
  82      +  4212, -4212, 4322, -4322, 4312, -4312, 2212, -2212, 3222, -3222,
  83
  84      +  4222, -4222, 2112, -2112, 3112, -3112, 4112, -4112, 3322, -3322,
  85
  86      +  3312, -3312, 4332, -4332, 6*0,
  87
  88      +  3214, -3214, 4214, -4214, 4324, -4324, 4314, -4314, 2214, -2214,
  89
  90      +  3224, -3224, 4224, -4224, 2114, -2114, 3114, -3114, 4114, -4114,
  91
  92      +  3324, -3324, 3314, -3314, 4334, -4334, 4*0,
  93
  94      +  0, 0,  2224, -2224, 1114, -1114, 3334, -3334, 0, 0,
  95
  96      +  10323, -10323, 20323, -20323, 6*0,
  97
  98      +  30443, 0, 0, 0, 70443, 50553, 60553, 80553, 20443, 0,
  99
 100      +  10411, 20413, 10413, 415,
 101
 102      + -10411,-20413,-10413,-415,
 103
 104      +  10421, 20423, 10423, 425,
 105
 106      + -10421,-20423,-10423,-425,
 107
 108      +  10431, 20433, 10433, 435,
 109
 110      + -10431,-20433,-10433,-435, 0,0,0,0,0,0,
 111
 112      +  10111, 10211,-10211, 115, 215, -215,10221,10331,20223,20333,
 113
 114      +  225, 335, 10223, 10333, 10113, 10213,-10213, 33*0 /
 115
 116       IF(NDIR.EQ.1) THEN
 117
 118         DO 10 I=1,321
 119
 120         IF (IDL.EQ.AKF(I)) THEN
 121
 122           QQLMAT=I-21
 123
 124           RETURN
 125
 126         ENDIF
 127
 128   10    CONTINUE
 129
 130         WRITE(6,20) IDL
 131
 132   20    FORMAT(1X,'Lund code particle ',I6,' not recognized')
 133
 134       ELSEIF(NDIR.EQ.2) THEN
 135
 136         QQLMAT = AKF(IDL+21)
 137
 138       ENDIF
 139
 140       RETURN
 141
 142       END