[u/mrichter/AliRoot.git] / HERWIG / src / qqlmat.f


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
Commit	Line	Data
3820ca8e	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