- extend debug output for GTU simulation
[u/mrichter/AliRoot.git] / TRD / AliTRDgtuParam.cxx
1 /**************************************************************************
2  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3  *                                                                        *
4  * Author: The ALICE Off-line Project.                                    *
5  * Contributors are mentioned in the code where appropriate.              *
6  *                                                                        *
7  * Permission to use, copy, modify and distribute this software and its   *
8  * documentation strictly for non-commercial purposes is hereby granted   *
9  * without fee, provided that the above copyright notice appears in all   *
10  * copies and that both the copyright notice and this permission notice   *
11  * appear in the supporting documentation. The authors make no claims     *
12  * about the suitability of this software for any purpose. It is          *
13  * provided "as is" without express or implied warranty.                  *
14  **************************************************************************/
15
16 /* $Id: AliTRDgtuParam.cxx 28397 2008-09-02 09:33:00Z cblume $ */
17
18 ////////////////////////////////////////////////////////////////////////////
19 //                                                                        //
20 //  Parameters for GTU simulation                                         //
21 //                                                                        //
22 //  Author: J. Klein (Jochen.Klein@cern.ch)                               //
23 //                                                                        //
24 ////////////////////////////////////////////////////////////////////////////
25
26 #include "TROOT.h"
27 #include "TMath.h"
28 #include "TMatrix.h"
29 #include "TDecompLU.h"
30 #include "TGraphAsymmErrors.h"
31 #include "TCanvas.h"
32
33 #include "AliLog.h"
34 #include "AliTRDgtuParam.h"
35 #include "AliTRDgeometry.h"
36 #include "AliTRDpadPlane.h"
37
38 ClassImp(AliTRDgtuParam)
39
40 AliTRDgtuParam *AliTRDgtuParam::fgInstance = 0;
41 Bool_t AliTRDgtuParam::fgUseGTUconst = kTRUE;
42
43 // ----- matching windows -----
44       Int_t     AliTRDgtuParam::fgDeltaY     = 19;
45       Int_t     AliTRDgtuParam::fgDeltaAlpha = 21;
46 // ----- reference layers -----
47       Int_t     AliTRDgtuParam::fgRefLayers[] = { 3, 2, 1 };
48
49 // ----- Bin widths (granularity) -----
50 const Float_t   AliTRDgtuParam::fgkBinWidthY  = 160e-4;
51 const Float_t   AliTRDgtuParam::fgkBinWidthdY = 140e-4;
52
53 // ----- Bit widths (used for internal representation) -----
54 const Int_t     AliTRDgtuParam::fgkBitWidthY      = 13;
55 const Int_t     AliTRDgtuParam::fgkBitWidthdY     = 7;
56 const Int_t     AliTRDgtuParam::fgkBitWidthYProj  = 10;
57 const Int_t     AliTRDgtuParam::fgkBitExcessY     = 4;
58 const Int_t     AliTRDgtuParam::fgkBitExcessAlpha = 10;
59 const Int_t     AliTRDgtuParam::fgkBitExcessYProj = 2;
60
61 // ----- z-channel tables -----
62 const Bool_t    AliTRDgtuParam::fgZChannelMap[5][16][6][16] = {
63
64 {  /* --- Stack 0 --- */
65
66 /*  . x x . . . . . . . . . . . . .  */
67 /*  x . . . . . . . . . . . . . . .  */
68 /*  X . . . . . . . . . . . . . . .  */
69 /*  x x . . . . . . . . . . . . . .  */
70 /*  x x . . . . . . . . . . . . . .  */
71 /*  x . . . . . . . . . . . . . . .  */
72
73 {{0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
74  {1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
75  {1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
76  {1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
77  {1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
78  {1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}},
79
80 /*  . . x x . . . . . . . . . . . .  */
81 /*  x x . . . . . . . . . . . . . .  */
82 /*  . X . . . . . . . . . . . . . .  */
83 /*  . x x . . . . . . . . . . . . .  */
84 /*  . x x . . . . . . . . . . . . .  */
85 /*  x x . . . . . . . . . . . . . .  */
86
87 {{0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
88  {1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
89  {0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
90  {0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
91  {0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
92  {1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}},
93
94 /*  . . . x x . . . . . . . . . . .  */
95 /*  . x x . . . . . . . . . . . . .  */
96 /*  . . X . . . . . . . . . . . . .  */
97 /*  . . x x . . . . . . . . . . . .  */
98 /*  . . x x . . . . . . . . . . . .  */
99 /*  . x x . . . . . . . . . . . . .  */
100
101 {{0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
102  {0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
103  {0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
104  {0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
105  {0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
106  {0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}},
107
108 /*  . . . . x x . . . . . . . . . .  */
109 /*  . . x x . . . . . . . . . . . .  */
110 /*  . . . X . . . . . . . . . . . .  */
111 /*  . . . x x . . . . . . . . . . .  */
112 /*  . . . x x . . . . . . . . . . .  */
113 /*  . . x x . . . . . . . . . . . .  */
114
115 {{0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
116  {0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
117  {0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
118  {0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
119  {0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
120  {0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}},
121
122 /*  . . . . . x x . . . . . . . . .  */
123 /*  . . . x x . . . . . . . . . . .  */
124 /*  . . . . X . . . . . . . . . . .  */
125 /*  . . . . x x . . . . . . . . . .  */
126 /*  . . . . x x . . . . . . . . . .  */
127 /*  . . . x x . . . . . . . . . . .  */
128
129 {{0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0},
130  {0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
131  {0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
132  {0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
133  {0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
134  {0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}},
135
136 /*  . . . . . . x x . . . . . . . .  */
137 /*  . . . . x x . . . . . . . . . .  */
138 /*  . . . . . X . . . . . . . . . .  */
139 /*  . . . . . x x . . . . . . . . .  */
140 /*  . . . . . x x . . . . . . . . .  */
141 /*  . . . . x x . . . . . . . . . .  */
142
143 {{0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0},
144  {0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
145  {0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
146  {0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0},
147  {0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0},
148  {0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}},
149
150 /*  . . . . . . . x x . . . . . . .  */
151 /*  . . . . . x x . . . . . . . . .  */
152 /*  . . . . . . X . . . . . . . . .  */
153 /*  . . . . . . x x . . . . . . . .  */
154 /*  . . . . . . x x . . . . . . . .  */
155 /*  . . . . . x x . . . . . . . . .  */
156
157 {{0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0},
158  {0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0},
159  {0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0},
160  {0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0},
161  {0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0},
162  {0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0}},
163
164 /*  . . . . . . . . x x . . . . . .  */
165 /*  . . . . . . x x x . . . . . . .  */
166 /*  . . . . . . . X . . . . . . . .  */
167 /*  . . . . . . . x x . . . . . . .  */
168 /*  . . . . . . . x x . . . . . . .  */
169 /*  . . . . . . x x . . . . . . . .  */
170
171 {{0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0},
172  {0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0},
173  {0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0},
174  {0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0},
175  {0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0},
176  {0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0}},
177
178 /*  . . . . . . . . x x x . . . . .  */
179 /*  . . . . . . . x x x . . . . . .  */
180 /*  . . . . . . . . X . . . . . . .  */
181 /*  . . . . . . . x x x . . . . . .  */
182 /*  . . . . . . . x x x . . . . . .  */
183 /*  . . . . . . . x x . . . . . . .  */
184
185 {{0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0},
186  {0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0},
187  {0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0},
188  {0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0},
189  {0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0},
190  {0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0}},
191
192 /*  . . . . . . . . . x x x . . . .  */
193 /*  . . . . . . . . x x x . . . . .  */
194 /*  . . . . . . . . . X . . . . . .  */
195 /*  . . . . . . . . x x x . . . . .  */
196 /*  . . . . . . . . x x x . . . . .  */
197 /*  . . . . . . . . x x . . . . . .  */
198
199 {{0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0},
200  {0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0},
201  {0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0},
202  {0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0},
203  {0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0},
204  {0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0}},
205
206 /*  . . . . . . . . . . x x x . . .  */
207 /*  . . . . . . . . . x x x . . . .  */
208 /*  . . . . . . . . . . X . . . . .  */
209 /*  . . . . . . . . . x x x . . . .  */
210 /*  . . . . . . . . . x x x . . . .  */
211 /*  . . . . . . . . . x x . . . . .  */
212
213 {{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0},
214  {0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0},
215  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0},
216  {0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0},
217  {0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0},
218  {0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0}},
219
220 /*  . . . . . . . . . . . x x x . .  */
221 /*  . . . . . . . . . . x x x . . .  */
222 /*  . . . . . . . . . . . X . . . .  */
223 /*  . . . . . . . . . . x x x . . .  */
224 /*  . . . . . . . . . . x x x . . .  */
225 /*  . . . . . . . . . . x x . . . .  */
226
227 {{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0},
228  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0},
229  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0},
230  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0},
231  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0},
232  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0}},
233
234 /*  . . . . . . . . . . . . x x x .  */
235 /*  . . . . . . . . . . . x x x . .  */
236 /*  . . . . . . . . . . . . X . . .  */
237 /*  . . . . . . . . . . . x x x . .  */
238 /*  . . . . . . . . . . . x x x . .  */
239 /*  . . . . . . . . . . . x x . . .  */
240
241 {{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0},
242  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0},
243  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0},
244  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0},
245  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0},
246  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0}},
247
248 /*  . . . . . . . . . . . . . x x x  */
249 /*  . . . . . . . . . . . . x x x .  */
250 /*  . . . . . . . . . . . . . X . .  */
251 /*  . . . . . . . . . . . . x x x .  */
252 /*  . . . . . . . . . . . . x x x .  */
253 /*  . . . . . . . . . . . . x x . .  */
254
255 {{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1},
256  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0},
257  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0},
258  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0},
259  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0},
260  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0}},
261
262 /*  . . . . . . . . . . . . . . x x  */
263 /*  . . . . . . . . . . . . . x x x  */
264 /*  . . . . . . . . . . . . . . X .  */
265 /*  . . . . . . . . . . . . . x x x  */
266 /*  . . . . . . . . . . . . . x x x  */
267 /*  . . . . . . . . . . . . . x x .  */
268
269 {{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1},
270  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1},
271  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0},
272  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1},
273  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1},
274  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0}},
275
276 /*  . . . . . . . . . . . . . . . x  */
277 /*  . . . . . . . . . . . . . . x x  */
278 /*  . . . . . . . . . . . . . . . X  */
279 /*  . . . . . . . . . . . . . . x x  */
280 /*  . . . . . . . . . . . . . . x x  */
281 /*  . . . . . . . . . . . . . . x x  */
282
283 {{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1},
284  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1},
285  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1},
286  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1},
287  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1},
288  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1}}},
289
290 {  /* --- Stack 1 --- */
291
292 /*  x x x . . . . . . . . . . . . .  */
293 /*  x x . . . . . . . . . . . . . .  */
294 /*  X . . . . . . . . . . . . . . .  */
295 /*  x x . . . . . . . . . . . . . .  */
296 /*  x x . . . . . . . . . . . . . .  */
297 /*  x . . . . . . . . . . . . . . .  */
298
299 {{1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
300  {1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
301  {1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
302  {1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
303  {1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
304  {1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}},
305
306 /*  . x x x . . . . . . . . . . . .  */
307 /*  x x x . . . . . . . . . . . . .  */
308 /*  . X . . . . . . . . . . . . . .  */
309 /*  x x x . . . . . . . . . . . . .  */
310 /*  x x x . . . . . . . . . . . . .  */
311 /*  x x . . . . . . . . . . . . . .  */
312
313 {{0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
314  {1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
315  {0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
316  {1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
317  {1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
318  {1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}},
319
320 /*  . . x x x . . . . . . . . . . .  */
321 /*  . x x x . . . . . . . . . . . .  */
322 /*  . . X . . . . . . . . . . . . .  */
323 /*  . x x x . . . . . . . . . . . .  */
324 /*  . x x x . . . . . . . . . . . .  */
325 /*  . x x . . . . . . . . . . . . .  */
326
327 {{0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
328  {0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
329  {0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
330  {0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
331  {0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
332  {0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}},
333
334 /*  . . . x x x . . . . . . . . . .  */
335 /*  . . x x x . . . . . . . . . . .  */
336 /*  . . . X . . . . . . . . . . . .  */
337 /*  . . x x x . . . . . . . . . . .  */
338 /*  . . x x x . . . . . . . . . . .  */
339 /*  . . x x . . . . . . . . . . . .  */
340
341 {{0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
342  {0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
343  {0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
344  {0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
345  {0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
346  {0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}},
347
348 /*  . . . . x x x . . . . . . . . .  */
349 /*  . . . x x x . . . . . . . . . .  */
350 /*  . . . . X . . . . . . . . . . .  */
351 /*  . . . x x x . . . . . . . . . .  */
352 /*  . . . x x x . . . . . . . . . .  */
353 /*  . . . x x . . . . . . . . . . .  */
354
355 {{0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0},
356  {0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
357  {0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
358  {0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
359  {0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
360  {0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}},
361
362 /*  . . . . . x x x . . . . . . . .  */
363 /*  . . . . x x x . . . . . . . . .  */
364 /*  . . . . . X . . . . . . . . . .  */
365 /*  . . . . x x . . . . . . . . . .  */
366 /*  . . . . x x . . . . . . . . . .  */
367 /*  . . . . x x . . . . . . . . . .  */
368
369 {{0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0},
370  {0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0},
371  {0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
372  {0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
373  {0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
374  {0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}},
375
376 /*  . . . . . . x x x . . . . . . .  */
377 /*  . . . . . . x x . . . . . . . .  */
378 /*  . . . . . . X . . . . . . . . .  */
379 /*  . . . . . x x . . . . . . . . .  */
380 /*  . . . . . x x . . . . . . . . .  */
381 /*  . . . . . x x . . . . . . . . .  */
382
383 {{0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0},
384  {0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0},
385  {0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0},
386  {0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0},
387  {0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0},
388  {0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0}},
389
390 /*  . . . . . . . x x . . . . . . .  */
391 /*  . . . . . . . x x . . . . . . .  */
392 /*  . . . . . . . X . . . . . . . .  */
393 /*  . . . . . . x x . . . . . . . .  */
394 /*  . . . . . . x x . . . . . . . .  */
395 /*  . . . . . . x x . . . . . . . .  */
396
397 {{0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0},
398  {0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0},
399  {0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0},
400  {0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0},
401  {0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0},
402  {0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0}},
403
404 /*  . . . . . . . . x x . . . . . .  */
405 /*  . . . . . . . . x x . . . . . .  */
406 /*  . . . . . . . . X . . . . . . .  */
407 /*  . . . . . . . x x . . . . . . .  */
408 /*  . . . . . . . x x . . . . . . .  */
409 /*  . . . . . . . x x . . . . . . .  */
410
411 {{0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0},
412  {0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0},
413  {0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0},
414  {0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0},
415  {0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0},
416  {0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0}},
417
418 /*  . . . . . . . . . x x . . . . .  */
419 /*  . . . . . . . . . x x . . . . .  */
420 /*  . . . . . . . . . X . . . . . .  */
421 /*  . . . . . . . . x x . . . . . .  */
422 /*  . . . . . . . . x x . . . . . .  */
423 /*  . . . . . . . . x x . . . . . .  */
424
425 {{0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0},
426  {0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0},
427  {0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0},
428  {0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0},
429  {0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0},
430  {0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0}},
431
432 /*  . . . . . . . . . . x x . . . .  */
433 /*  . . . . . . . . . . x x . . . .  */
434 /*  . . . . . . . . . . X . . . . .  */
435 /*  . . . . . . . . . x x . . . . .  */
436 /*  . . . . . . . . . x x . . . . .  */
437 /*  . . . . . . . . . x x . . . . .  */
438
439 {{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0},
440  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0},
441  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0},
442  {0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0},
443  {0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0},
444  {0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0}},
445
446 /*  . . . . . . . . . . . x x . . .  */
447 /*  . . . . . . . . . . . x x . . .  */
448 /*  . . . . . . . . . . . X . . . .  */
449 /*  . . . . . . . . . . x x . . . .  */
450 /*  . . . . . . . . . . x x . . . .  */
451 /*  . . . . . . . . . . x x . . . .  */
452
453 {{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0},
454  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0},
455  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0},
456  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0},
457  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0},
458  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0}},
459
460 /*  . . . . . . . . . . . . x x . .  */
461 /*  . . . . . . . . . . . . x x . .  */
462 /*  . . . . . . . . . . . . X . . .  */
463 /*  . . . . . . . . . . . x x . . .  */
464 /*  . . . . . . . . . . . x x . . .  */
465 /*  . . . . . . . . . . . x x . . .  */
466
467 {{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0},
468  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0},
469  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0},
470  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0},
471  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0},
472  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0}},
473
474 /*  . . . . . . . . . . . . . x x .  */
475 /*  . . . . . . . . . . . . . x x .  */
476 /*  . . . . . . . . . . . . . X . .  */
477 /*  . . . . . . . . . . . . x x . .  */
478 /*  . . . . . . . . . . . . x x . .  */
479 /*  . . . . . . . . . . . . x x . .  */
480
481 {{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0},
482  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0},
483  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0},
484  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0},
485  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0},
486  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0}},
487
488 /*  . . . . . . . . . . . . . . x x  */
489 /*  . . . . . . . . . . . . . . x x  */
490 /*  . . . . . . . . . . . . . . X .  */
491 /*  . . . . . . . . . . . . . x x .  */
492 /*  . . . . . . . . . . . . . x x .  */
493 /*  . . . . . . . . . . . . . x x .  */
494
495 {{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1},
496  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1},
497  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0},
498  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0},
499  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0},
500  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0}},
501
502 /*  . . . . . . . . . . . . . . . x  */
503 /*  . . . . . . . . . . . . . . . x  */
504 /*  . . . . . . . . . . . . . . . X  */
505 /*  . . . . . . . . . . . . . . x x  */
506 /*  . . . . . . . . . . . . . . x x  */
507 /*  . . . . . . . . . . . . . . x x  */
508
509 {{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1},
510  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1},
511  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1},
512  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1},
513  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1},
514  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1}}},
515
516 {  /* --- Stack 2 --- */
517
518 /*  x x . . . . . . . . . .          */
519 /*  x x . . . . . . . . . .          */
520 /*  X . . . . . . . . . . .          */
521 /*  x . . . . . . . . . . .          */
522 /*  x . . . . . . . . . . .          */
523 /*  x . . . . . . . . . . .          */
524
525 {{1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
526  {1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
527  {1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
528  {1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
529  {1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
530  {1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}},
531
532 /*  . x x . . . . . . . . .          */
533 /*  . x x . . . . . . . . .          */
534 /*  . X . . . . . . . . . .          */
535 /*  x x . . . . . . . . . .          */
536 /*  x x . . . . . . . . . .          */
537 /*  x x . . . . . . . . . .          */
538
539 {{0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
540  {0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
541  {0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
542  {1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
543  {1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
544  {1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}},
545
546 /*  . . x x . . . . . . . .          */
547 /*  . . x x . . . . . . . .          */
548 /*  . . X . . . . . . . . .          */
549 /*  . x x . . . . . . . . .          */
550 /*  . x x . . . . . . . . .          */
551 /*  . x x . . . . . . . . .          */
552
553 {{0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
554  {0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
555  {0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
556  {0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
557  {0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
558  {0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}},
559
560 /*  . . . x x . . . . . . .          */
561 /*  . . . x x . . . . . . .          */
562 /*  . . . X . . . . . . . .          */
563 /*  . . x x x . . . . . . .          */
564 /*  . . x x x . . . . . . .          */
565 /*  . . x x x . . . . . . .          */
566
567 {{0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
568  {0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
569  {0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
570  {0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
571  {0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
572  {0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}},
573
574 /*  . . . x x x . . . . . .          */
575 /*  . . . x x x . . . . . .          */
576 /*  . . . . X . . . . . . .          */
577 /*  . . . x x x . . . . . .          */
578 /*  . . . x x x . . . . . .          */
579 /*  . . . x x x . . . . . .          */
580
581 {{0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
582  {0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
583  {0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
584  {0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
585  {0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
586  {0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}},
587
588 /*  . . . . x x x . . . . .          */
589 /*  . . . . x x x . . . . .          */
590 /*  . . . . . X . . . . . .          */
591 /*  . . . . x x x . . . . .          */
592 /*  . . . . x x x . . . . .          */
593 /*  . . . . x x x . . . . .          */
594
595 {{0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0},
596  {0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0},
597  {0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
598  {0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0},
599  {0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0},
600  {0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0}},
601
602 /*  . . . . . x x x . . . .          */
603 /*  . . . . . x x x . . . .          */
604 /*  . . . . . . X . . . . .          */
605 /*  . . . . . x x x . . . .          */
606 /*  . . . . . x x x . . . .          */
607 /*  . . . . . x x x . . . .          */
608
609 {{0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0},
610  {0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0},
611  {0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0},
612  {0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0},
613  {0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0},
614  {0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0}},
615
616 /*  . . . . . . x x x . . .          */
617 /*  . . . . . . x x x . . .          */
618 /*  . . . . . . . X . . . .          */
619 /*  . . . . . . x x x . . .          */
620 /*  . . . . . . x x x . . .          */
621 /*  . . . . . . x x x . . .          */
622
623 {{0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0},
624  {0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0},
625  {0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0},
626  {0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0},
627  {0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0},
628  {0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0}},
629
630 /*  . . . . . . . x x . . .          */
631 /*  . . . . . . . x x . . .          */
632 /*  . . . . . . . . X . . .          */
633 /*  . . . . . . . x x x . .          */
634 /*  . . . . . . . x x x . .          */
635 /*  . . . . . . . x x x . .          */
636
637 {{0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0},
638  {0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0},
639  {0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0},
640  {0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0},
641  {0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0},
642  {0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0}},
643
644 /*  . . . . . . . . x x . .          */
645 /*  . . . . . . . . x x . .          */
646 /*  . . . . . . . . . X . .          */
647 /*  . . . . . . . . . x x .          */
648 /*  . . . . . . . . . x x .          */
649 /*  . . . . . . . . . x x .          */
650
651 {{0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0},
652  {0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0},
653  {0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0},
654  {0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0},
655  {0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0},
656  {0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0}},
657
658 /*  . . . . . . . . . x x .          */
659 /*  . . . . . . . . . x x .          */
660 /*  . . . . . . . . . . X .          */
661 /*  . . . . . . . . . . x x          */
662 /*  . . . . . . . . . . x x          */
663 /*  . . . . . . . . . . x x          */
664
665 {{0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0},
666  {0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0},
667  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0},
668  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0},
669  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0},
670  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0}},
671
672 /*  . . . . . . . . . . x x          */
673 /*  . . . . . . . . . . x x          */
674 /*  . . . . . . . . . . . X          */
675 /*  . . . . . . . . . . . x          */
676 /*  . . . . . . . . . . . x          */
677 /*  . . . . . . . . . . . x          */
678
679 {{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0},
680  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0},
681  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0},
682  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0},
683  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0},
684  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0}},
685
686 /*  . . . . . . . . . . . .          */
687 /*  . . . . . . . . . . . .          */
688 /*  . . . . . . . . . . . .          */
689 /*  . . . . . . . . . . . .          */
690 /*  . . . . . . . . . . . .          */
691 /*  . . . . . . . . . . . .          */
692
693 {{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
694  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
695  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
696  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
697  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
698  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}},
699
700 /*  . . . . . . . . . . . .          */
701 /*  . . . . . . . . . . . .          */
702 /*  . . . . . . . . . . . .          */
703 /*  . . . . . . . . . . . .          */
704 /*  . . . . . . . . . . . .          */
705 /*  . . . . . . . . . . . .          */
706
707 {{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
708  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
709  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
710  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
711  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
712  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}},
713
714 /*  . . . . . . . . . . . .          */
715 /*  . . . . . . . . . . . .          */
716 /*  . . . . . . . . . . . .          */
717 /*  . . . . . . . . . . . .          */
718 /*  . . . . . . . . . . . .          */
719 /*  . . . . . . . . . . . .          */
720
721 {{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
722  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
723  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
724  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
725  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
726  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}},
727
728 /*  . . . . . . . . . . . .          */
729 /*  . . . . . . . . . . . .          */
730 /*  . . . . . . . . . . . .          */
731 /*  . . . . . . . . . . . .          */
732 /*  . . . . . . . . . . . .          */
733 /*  . . . . . . . . . . . .          */
734
735 {{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
736  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
737  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
738  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
739  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
740  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}}},
741
742 {  /* --- Stack 3 --- */
743
744 /*  x . . . . . . . . . . . . . . .  */
745 /*  x . . . . . . . . . . . . . . .  */
746 /*  X . . . . . . . . . . . . . . .  */
747 /*  x x . . . . . . . . . . . . . .  */
748 /*  x x . . . . . . . . . . . . . .  */
749 /*  x x . . . . . . . . . . . . . .  */
750
751 {{1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
752  {1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
753  {1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
754  {1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
755  {1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
756  {1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}},
757
758 /*  x x . . . . . . . . . . . . . .  */
759 /*  x x . . . . . . . . . . . . . .  */
760 /*  . X . . . . . . . . . . . . . .  */
761 /*  . x x . . . . . . . . . . . . .  */
762 /*  . x x . . . . . . . . . . . . .  */
763 /*  . x x . . . . . . . . . . . . .  */
764
765 {{1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
766  {1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
767  {0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
768  {0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
769  {0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
770  {0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}},
771
772 /*  . x x . . . . . . . . . . . . .  */
773 /*  . x x . . . . . . . . . . . . .  */
774 /*  . . X . . . . . . . . . . . . .  */
775 /*  . . x x . . . . . . . . . . . .  */
776 /*  . . x x . . . . . . . . . . . .  */
777 /*  . . x x . . . . . . . . . . . .  */
778
779 {{0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
780  {0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
781  {0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
782  {0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
783  {0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
784  {0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}},
785
786 /*  . . x x . . . . . . . . . . . .  */
787 /*  . . x x . . . . . . . . . . . .  */
788 /*  . . . X . . . . . . . . . . . .  */
789 /*  . . . x x . . . . . . . . . . .  */
790 /*  . . . x x . . . . . . . . . . .  */
791 /*  . . . x x . . . . . . . . . . .  */
792
793 {{0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
794  {0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
795  {0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
796  {0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
797  {0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
798  {0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}},
799
800 /*  . . . x x . . . . . . . . . . .  */
801 /*  . . . x x . . . . . . . . . . .  */
802 /*  . . . . X . . . . . . . . . . .  */
803 /*  . . . . x x . . . . . . . . . .  */
804 /*  . . . . x x . . . . . . . . . .  */
805 /*  . . . . x x . . . . . . . . . .  */
806
807 {{0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
808  {0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
809  {0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
810  {0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
811  {0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
812  {0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}},
813
814 /*  . . . . x x . . . . . . . . . .  */
815 /*  . . . . x x . . . . . . . . . .  */
816 /*  . . . . . X . . . . . . . . . .  */
817 /*  . . . . . x x . . . . . . . . .  */
818 /*  . . . . . x x . . . . . . . . .  */
819 /*  . . . . . x x . . . . . . . . .  */
820
821 {{0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
822  {0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
823  {0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
824  {0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0},
825  {0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0},
826  {0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0}},
827
828 /*  . . . . . x x . . . . . . . . .  */
829 /*  . . . . . x x . . . . . . . . .  */
830 /*  . . . . . . X . . . . . . . . .  */
831 /*  . . . . . . x x . . . . . . . .  */
832 /*  . . . . . . x x . . . . . . . .  */
833 /*  . . . . . . x x . . . . . . . .  */
834
835 {{0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0},
836  {0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0},
837  {0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0},
838  {0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0},
839  {0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0},
840  {0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0}},
841
842 /*  . . . . . . x x . . . . . . . .  */
843 /*  . . . . . . x x . . . . . . . .  */
844 /*  . . . . . . . X . . . . . . . .  */
845 /*  . . . . . . . x x . . . . . . .  */
846 /*  . . . . . . . x x . . . . . . .  */
847 /*  . . . . . . . x x . . . . . . .  */
848
849 {{0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0},
850  {0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0},
851  {0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0},
852  {0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0},
853  {0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0},
854  {0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0}},
855
856 /*  . . . . . . . x x . . . . . . .  */
857 /*  . . . . . . . x x . . . . . . .  */
858 /*  . . . . . . . . X . . . . . . .  */
859 /*  . . . . . . . . x x . . . . . .  */
860 /*  . . . . . . . . x x . . . . . .  */
861 /*  . . . . . . . . x x . . . . . .  */
862
863 {{0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0},
864  {0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0},
865  {0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0},
866  {0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0},
867  {0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0},
868  {0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0}},
869
870 /*  . . . . . . . x x x . . . . . .  */
871 /*  . . . . . . . . x x . . . . . .  */
872 /*  . . . . . . . . . X . . . . . .  */
873 /*  . . . . . . . . . x x . . . . .  */
874 /*  . . . . . . . . . x x . . . . .  */
875 /*  . . . . . . . . . x x . . . . .  */
876
877 {{0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0},
878  {0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0},
879  {0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0},
880  {0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0},
881  {0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0},
882  {0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0}},
883
884 /*  . . . . . . . . x x x . . . . .  */
885 /*  . . . . . . . . . x x x . . . .  */
886 /*  . . . . . . . . . . X . . . . .  */
887 /*  . . . . . . . . . . x x . . . .  */
888 /*  . . . . . . . . . . x x . . . .  */
889 /*  . . . . . . . . . . x x . . . .  */
890
891 {{0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0},
892  {0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0},
893  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0},
894  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0},
895  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0},
896  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0}},
897
898 /*  . . . . . . . . . x x x . . . .  */
899 /*  . . . . . . . . . . x x x . . .  */
900 /*  . . . . . . . . . . . X . . . .  */
901 /*  . . . . . . . . . . x x x . . .  */
902 /*  . . . . . . . . . . x x x . . .  */
903 /*  . . . . . . . . . . . x x . . .  */
904
905 {{0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0},
906  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0},
907  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0},
908  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0},
909  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0},
910  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0}},
911
912 /*  . . . . . . . . . . x x x . . .  */
913 /*  . . . . . . . . . . . x x x . .  */
914 /*  . . . . . . . . . . . . X . . .  */
915 /*  . . . . . . . . . . . x x x . .  */
916 /*  . . . . . . . . . . . x x x . .  */
917 /*  . . . . . . . . . . . . x x . .  */
918
919 {{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0},
920  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0},
921  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0},
922  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0},
923  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0},
924  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0}},
925
926 /*  . . . . . . . . . . . x x x . .  */
927 /*  . . . . . . . . . . . . x x x .  */
928 /*  . . . . . . . . . . . . . X . .  */
929 /*  . . . . . . . . . . . . x x x .  */
930 /*  . . . . . . . . . . . . x x x .  */
931 /*  . . . . . . . . . . . . . x x .  */
932
933 {{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0},
934  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0},
935  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0},
936  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0},
937  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0},
938  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0}},
939
940 /*  . . . . . . . . . . . . x x x .  */
941 /*  . . . . . . . . . . . . . x x x  */
942 /*  . . . . . . . . . . . . . . X .  */
943 /*  . . . . . . . . . . . . . x x x  */
944 /*  . . . . . . . . . . . . . x x x  */
945 /*  . . . . . . . . . . . . . . x x  */
946
947 {{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0},
948  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1},
949  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0},
950  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1},
951  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1},
952  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1}},
953
954 /*  . . . . . . . . . . . . . x x x  */
955 /*  . . . . . . . . . . . . . . x x  */
956 /*  . . . . . . . . . . . . . . . X  */
957 /*  . . . . . . . . . . . . . . x x  */
958 /*  . . . . . . . . . . . . . . x x  */
959 /*  . . . . . . . . . . . . . . . x  */
960
961 {{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1},
962  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1},
963  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1},
964  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1},
965  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1},
966  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1}}},
967
968 {  /* --- Stack 4 --- */
969
970 /*  x . . . . . . . . . . . . . . .  */
971 /*  x x . . . . . . . . . . . . . .  */
972 /*  X . . . . . . . . . . . . . . .  */
973 /*  x x . . . . . . . . . . . . . .  */
974 /*  x x . . . . . . . . . . . . . .  */
975 /*  x x . . . . . . . . . . . . . .  */
976
977 {{1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
978  {1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
979  {1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
980  {1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
981  {1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
982  {1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}},
983
984 /*  x x . . . . . . . . . . . . . .  */
985 /*  x x x . . . . . . . . . . . . .  */
986 /*  . X . . . . . . . . . . . . . .  */
987 /*  x x x . . . . . . . . . . . . .  */
988 /*  x x x . . . . . . . . . . . . .  */
989 /*  . x x . . . . . . . . . . . . .  */
990
991 {{1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
992  {1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
993  {0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
994  {1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
995  {1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
996  {0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}},
997
998 /*  x x x . . . . . . . . . . . . .  */
999 /*  . x x x . . . . . . . . . . . .  */
1000 /*  . . X . . . . . . . . . . . . .  */
1001 /*  . x x x . . . . . . . . . . . .  */
1002 /*  . x x x . . . . . . . . . . . .  */
1003 /*  . . x x . . . . . . . . . . . .  */
1004
1005 {{1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
1006  {0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
1007  {0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
1008  {0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
1009  {0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
1010  {0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}},
1011
1012 /*  . x x x . . . . . . . . . . . .  */
1013 /*  . . x x x . . . . . . . . . . .  */
1014 /*  . . . X . . . . . . . . . . . .  */
1015 /*  . . x x x . . . . . . . . . . .  */
1016 /*  . . x x x . . . . . . . . . . .  */
1017 /*  . . . x x . . . . . . . . . . .  */
1018
1019 {{0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
1020  {0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
1021  {0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
1022  {0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
1023  {0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
1024  {0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}},
1025
1026 /*  . . x x x . . . . . . . . . . .  */
1027 /*  . . . x x x . . . . . . . . . .  */
1028 /*  . . . . X . . . . . . . . . . .  */
1029 /*  . . . x x x . . . . . . . . . .  */
1030 /*  . . . x x x . . . . . . . . . .  */
1031 /*  . . . . x x . . . . . . . . . .  */
1032
1033 {{0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
1034  {0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
1035  {0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
1036  {0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
1037  {0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
1038  {0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}},
1039
1040 /*  . . . x x x . . . . . . . . . .  */
1041 /*  . . . . x x x . . . . . . . . .  */
1042 /*  . . . . . X . . . . . . . . . .  */
1043 /*  . . . . x x x . . . . . . . . .  */
1044 /*  . . . . x x x . . . . . . . . .  */
1045 /*  . . . . . x x . . . . . . . . .  */
1046
1047 {{0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
1048  {0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0},
1049  {0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
1050  {0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0},
1051  {0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0},
1052  {0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0}},
1053
1054 /*  . . . . x x x . . . . . . . . .  */
1055 /*  . . . . . x x x . . . . . . . .  */
1056 /*  . . . . . . X . . . . . . . . .  */
1057 /*  . . . . . x x x . . . . . . . .  */
1058 /*  . . . . . x x x . . . . . . . .  */
1059 /*  . . . . . . x x . . . . . . . .  */
1060
1061 {{0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0},
1062  {0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0},
1063  {0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0},
1064  {0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0},
1065  {0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0},
1066  {0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0}},
1067
1068 /*  . . . . . x x x . . . . . . . .  */
1069 /*  . . . . . . x x x . . . . . . .  */
1070 /*  . . . . . . . X . . . . . . . .  */
1071 /*  . . . . . . x x x . . . . . . .  */
1072 /*  . . . . . . x x x . . . . . . .  */
1073 /*  . . . . . . . x x . . . . . . .  */
1074
1075 {{0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0},
1076  {0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0},
1077  {0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0},
1078  {0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0},
1079  {0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0},
1080  {0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0}},
1081
1082 /*  . . . . . . x x . . . . . . . .  */
1083 /*  . . . . . . . x x x . . . . . .  */
1084 /*  . . . . . . . . X . . . . . . .  */
1085 /*  . . . . . . . x x . . . . . . .  */
1086 /*  . . . . . . . x x . . . . . . .  */
1087 /*  . . . . . . . . x x . . . . . .  */
1088
1089 {{0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0},
1090  {0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0},
1091  {0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0},
1092  {0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0},
1093  {0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0},
1094  {0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0}},
1095
1096 /*  . . . . . . . x x . . . . . . .  */
1097 /*  . . . . . . . . . x x . . . . .  */
1098 /*  . . . . . . . . . X . . . . . .  */
1099 /*  . . . . . . . . x x . . . . . .  */
1100 /*  . . . . . . . . x x . . . . . .  */
1101 /*  . . . . . . . . . x x . . . . .  */
1102
1103 {{0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0},
1104  {0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0},
1105  {0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0},
1106  {0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0},
1107  {0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0},
1108  {0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0}},
1109
1110 /*  . . . . . . . . x x . . . . . .  */
1111 /*  . . . . . . . . . . x x . . . .  */
1112 /*  . . . . . . . . . . X . . . . .  */
1113 /*  . . . . . . . . . x x . . . . .  */
1114 /*  . . . . . . . . . x x . . . . .  */
1115 /*  . . . . . . . . . . x x . . . .  */
1116
1117 {{0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0},
1118  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0},
1119  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0},
1120  {0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0},
1121  {0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0},
1122  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0}},
1123
1124 /*  . . . . . . . . . x x . . . . .  */
1125 /*  . . . . . . . . . . . x x . . .  */
1126 /*  . . . . . . . . . . . X . . . .  */
1127 /*  . . . . . . . . . . x x . . . .  */
1128 /*  . . . . . . . . . . x x . . . .  */
1129 /*  . . . . . . . . . . . x x . . .  */
1130
1131 {{0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0},
1132  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0},
1133  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0},
1134  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0},
1135  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0},
1136  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0}},
1137
1138 /*  . . . . . . . . . . x x . . . .  */
1139 /*  . . . . . . . . . . . . x x . .  */
1140 /*  . . . . . . . . . . . . X . . .  */
1141 /*  . . . . . . . . . . . x x . . .  */
1142 /*  . . . . . . . . . . . x x . . .  */
1143 /*  . . . . . . . . . . . . x x . .  */
1144
1145 {{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0},
1146  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0},
1147  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0},
1148  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0},
1149  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0},
1150  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0}},
1151
1152 /*  . . . . . . . . . . . x x . . .  */
1153 /*  . . . . . . . . . . . . . x x .  */
1154 /*  . . . . . . . . . . . . . X . .  */
1155 /*  . . . . . . . . . . . . x x . .  */
1156 /*  . . . . . . . . . . . . x x . .  */
1157 /*  . . . . . . . . . . . . . x x .  */
1158
1159 {{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0},
1160  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0},
1161  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0},
1162  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0},
1163  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0},
1164  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0}},
1165
1166 /*  . . . . . . . . . . . . x x . .  */
1167 /*  . . . . . . . . . . . . . . x x  */
1168 /*  . . . . . . . . . . . . . . X .  */
1169 /*  . . . . . . . . . . . . . x x .  */
1170 /*  . . . . . . . . . . . . . x x .  */
1171 /*  . . . . . . . . . . . . . . x x  */
1172
1173 {{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0},
1174  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1},
1175  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0},
1176  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0},
1177  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0},
1178  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1}},
1179
1180 /*  . . . . . . . . . . . . . x x .  */
1181 /*  . . . . . . . . . . . . . . . x  */
1182 /*  . . . . . . . . . . . . . . . X  */
1183 /*  . . . . . . . . . . . . . . x x  */
1184 /*  . . . . . . . . . . . . . . x x  */
1185 /*  . . . . . . . . . . . . . . . x  */
1186
1187 {{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0},
1188  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1},
1189  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1},
1190  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1},
1191  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1},
1192  {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1}}}
1193
1194 };
1195
1196 AliTRDgtuParam::AliTRDgtuParam() :
1197   fVertexSize(20.0),
1198   fCurrTrackletMask(0),
1199   fMagField(0.5),
1200   fGeo(0x0)
1201 {
1202   // default ctor
1203   fGeo = new AliTRDgeometry();
1204   for (Int_t iLayer = 0; iLayer < 6; iLayer++) {
1205     fAki[iLayer] = 0.;
1206     fBki[iLayer] = 0.;
1207     fCki[iLayer] = 0.;
1208   }
1209
1210   GenerateZChannelMap();
1211 }
1212
1213 AliTRDgtuParam::~AliTRDgtuParam()
1214 {
1215   // dtor
1216
1217   delete fGeo;
1218 }
1219
1220 AliTRDgtuParam* AliTRDgtuParam::Instance()
1221 {
1222   // get (or create) the single instance
1223
1224   if (fgInstance == 0)
1225     fgInstance = new AliTRDgtuParam();
1226
1227   return fgInstance;
1228 }
1229
1230 void AliTRDgtuParam::Terminate()
1231 {
1232   // destruct the instance
1233
1234   if (fgInstance != 0) {
1235     delete fgInstance;
1236     fgInstance = 0x0;
1237   }
1238 }
1239
1240 Bool_t AliTRDgtuParam::IsInZChannel(Int_t stack, Int_t layer, Int_t zchannel, Int_t zpos) const
1241 {
1242   return (fZSubChannel[stack][zchannel][layer][zpos] != 0);
1243 }
1244
1245 Int_t AliTRDgtuParam::GetZSubchannel(Int_t stack, Int_t layer, Int_t zchannel, Int_t zpos) const
1246 {
1247   return fZSubChannel[stack][zchannel][layer][zpos];
1248 }
1249
1250 Int_t AliTRDgtuParam::GetRefLayer(Int_t refLayerIdx)
1251 {
1252   // returns the reference layer indexed by refLayerIdx
1253
1254   if (refLayerIdx >= 0 && refLayerIdx < fgkNRefLayers)
1255     return fgRefLayers[refLayerIdx];
1256   else
1257     return -1;
1258 }
1259
1260 Int_t AliTRDgtuParam::GenerateZChannelMap()
1261 {
1262   // generate the z-channel map
1263   // assuming that the tracks come from the vertex
1264   // +/- fVertexSize in z-direction
1265
1266   if (fgUseGTUconst) {
1267     for (Int_t iStack = 0; iStack < fGeo->Nstack(); iStack++) {
1268       for (Int_t iChannel = 0; iChannel < fGeo->GetRowMax(fgkFixLayer, iStack, 0); iChannel++) {
1269         for (Int_t iLayer = 0; iLayer < fGeo->Nlayer(); iLayer++) {
1270           for (Int_t iRow = 0; iRow < fGeo->GetRowMax(iLayer, iStack, 0); iRow++) {
1271             if (fgZChannelMap[iStack][iChannel][iLayer][iRow] != 0) {
1272               fZChannelMap[iStack][iChannel][iLayer][iRow] = 1;
1273               fZSubChannel[iStack][iChannel % fgkNZChannels][iLayer][iRow] = iChannel / fgkNZChannels + 1;
1274             }
1275           }
1276         }
1277       }
1278     }
1279
1280     return kTRUE;
1281   }
1282   else {
1283     Int_t iSec = 0; // sector is irrelevant
1284     Bool_t collision = kFALSE;
1285
1286     for (Int_t iStack = 0; iStack < fGeo->Nstack(); iStack++) {
1287
1288       Float_t x[6] = { 0 };
1289       Float_t z[6][16] = {{ 0 }};
1290       Float_t dZ[6][16] = {{ 0 }};
1291
1292       for (Int_t iLayer = 0; iLayer < fGeo->Nlayer(); iLayer++) {
1293         AliTRDpadPlane *pp = fGeo->GetPadPlane(iLayer, iStack);
1294         x[iLayer]  = fGeo->GetTime0(iLayer) - fGeo->CdrHght(); // ???
1295         for (Int_t iRow = 0; iRow < fGeo->GetRowMax(iLayer, iStack, iSec); iRow++) {
1296           z[iLayer][iRow]  = pp->GetRowPos(iRow); // this is the right (pos. z-direction) border of the pad
1297           dZ[iLayer][iRow] = pp->GetRowSize(iRow); // length of the pad in z-direction
1298           for (Int_t i = 0; i < fgkNZChannels; i++)
1299             fZSubChannel[iStack][i][iLayer][iRow] = 0;
1300         }
1301       }
1302
1303       for (Int_t fixRow = 0; fixRow < fGeo->GetRowMax(fgkFixLayer, iStack, iSec); fixRow++) {
1304
1305         Double_t fixZmin = z[fgkFixLayer][fixRow] - dZ[fgkFixLayer][fixRow];
1306         Double_t fixZmax = z[fgkFixLayer][fixRow];
1307         Double_t fixX    = x[fgkFixLayer] + 1.5; // ??? 1.5 from where?
1308
1309         for (Int_t iLayer = 0; iLayer < fGeo->Nlayer(); iLayer++) {
1310           Double_t leftZ, rightZ;
1311
1312           if (iLayer <= fgkFixLayer) {
1313             leftZ  = (fixZmin + fVertexSize) * (x[iLayer] + 1.5) / fixX - fVertexSize;
1314             rightZ = (fixZmax - fVertexSize) * (x[iLayer] + 1.5) / fixX + fVertexSize;
1315           }
1316           else {
1317             leftZ  = (fixZmin - fVertexSize) * (x[iLayer] + 1.5) / fixX + fVertexSize;
1318             rightZ = (fixZmax + fVertexSize) * (x[iLayer] + 1.5) / fixX - fVertexSize;
1319           }
1320
1321           Double_t epsilon = 0.001;
1322           for (Int_t iRow = 0; iRow < fGeo->GetRowMax(iLayer, iStack, iSec); iRow++) {
1323             if ( (z[iLayer][iRow] )                    > (leftZ  + epsilon) &&
1324                  (z[iLayer][iRow] - dZ[iLayer][iRow] ) < (rightZ - epsilon) ) {
1325               fZChannelMap[iStack][fixRow][iLayer][iRow] = 1;
1326               if (fZSubChannel[iStack][fixRow % fgkNZChannels][iLayer][iRow] != 0) {
1327                 AliError("Collision in Z-Channel assignment occured! No reliable tracking!!!");
1328                 collision = kTRUE;
1329               }
1330               else
1331                 fZSubChannel[iStack][fixRow % fgkNZChannels][iLayer][iRow] = fixRow / fgkNZChannels + 1;
1332             }
1333
1334           }
1335         }
1336       }
1337     }
1338
1339     return ~collision;
1340   }
1341 }
1342
1343 Bool_t AliTRDgtuParam::DisplayZChannelMap(Int_t zchannel, Int_t subchannel) const
1344 {
1345   // display the z-channel map
1346
1347   if (zchannel >= fgkNZChannels) {
1348     AliError("Invalid Z channel!");
1349     return kFALSE;
1350   }
1351
1352   Int_t zchmin = zchannel >= 0 ? zchannel : 0;
1353   Int_t zchmax = zchannel >= 0 ? zchannel + 1 : fgkNZChannels;
1354   Int_t i = 0;
1355   Int_t j = 0;
1356   TCanvas *c = new TCanvas("zchmap", "Z-Chhannel Mapping");
1357   c->cd();
1358   TGraph **graphz = new TGraph*[fgkNZChannels];
1359   for (Int_t zch = zchmin; zch < zchmax; zch++)
1360     graphz[zch] = new TGraph;
1361   TGraphAsymmErrors *graph = new TGraphAsymmErrors();
1362   graph->SetTitle("Z-Channel Map");
1363   graph->SetPoint(i, 0, 0); // vertex
1364   graph->SetPointError(i++, 20, 20, 0, 0);
1365   //  graph->SetRange //????
1366   for (Int_t iLayer = 0; iLayer < fGeo->Nlayer(); iLayer++) {
1367     for (Int_t iStack = 0; iStack < fGeo->Nstack(); iStack++) {
1368       AliTRDpadPlane *pp = fGeo->GetPadPlane(iLayer, iStack);
1369       for (Int_t iRow = 0; iRow < fGeo->GetRowMax(iLayer, iStack, 0); iRow++) {
1370         graph->SetPoint(i, pp->GetRowPos(iRow), fGeo->GetTime0(iLayer) - fGeo->CdrHght());
1371         graph->SetPointError(i++, pp->GetRowSize(iRow), 0, 0, 0);
1372         for (Int_t zch = zchmin; zch < zchmax; zch++)
1373           if (fZSubChannel[iStack][zch][iLayer][iRow] != 0)
1374             if (subchannel == 0 || fZSubChannel[iStack][zch][iLayer][iRow] == subchannel)
1375               graphz[zch]->SetPoint(j++, pp->GetRowPos(iRow)  - pp->GetRowSize(iRow)/2, fGeo->GetTime0(iLayer) - fGeo->CdrHght());
1376       }
1377     }
1378   }
1379   graph->SetMarkerStyle(kDot);
1380   graph->Draw("AP");
1381   gROOT->Add(graph);
1382   for (Int_t zch = zchmin; zch < zchmax; zch++) {
1383     graphz[zch]->SetMarkerStyle(kCircle);
1384     graphz[zch]->SetMarkerColor(zch+2);
1385     graphz[zch]->SetMarkerSize(0.3 + zch*0.2);
1386     graphz[zch]->Draw("P");
1387     gROOT->Add(graphz[zch]);
1388   }
1389   delete [] graphz;
1390   return kTRUE;
1391 }
1392
1393 Int_t AliTRDgtuParam::GetCiAlpha(Int_t layer) const
1394 {
1395   // get the constant for the calculation of alpha
1396
1397   Int_t ci = TMath::Nint(GetChamberThickness() / fGeo->GetTime0(layer) * GetBinWidthY() / GetBinWidthdY() * (1 << (GetBitExcessAlpha() + GetBitExcessY() + 1)) );
1398   return ci;
1399 }
1400
1401 Int_t AliTRDgtuParam::GetCiYProj(Int_t layer) const
1402 {
1403   // get the constant for the calculation of y_proj
1404
1405   Float_t xmid = (fGeo->GetTime0(0) + fGeo->GetTime0(fGeo->Nlayer()-1)) / 2.;
1406   Int_t ci = TMath::Nint(- (fGeo->GetTime0(layer) - xmid) / GetChamberThickness() * GetBinWidthdY() / GetBinWidthY() * (1 << GetBitExcessYProj()) );
1407   return ci;
1408 }
1409
1410 Int_t AliTRDgtuParam::GetYt(Int_t stack, Int_t layer, Int_t zrow) const
1411 {
1412     return (Int_t) (- ( (layer % 2 ? 1 : -1) *
1413                         (GetGeo()->GetPadPlane(layer, stack)->GetRowPos(zrow) - GetGeo()->GetPadPlane(layer, stack)->GetRowSize(zrow) / 2) *
1414                         TMath::Tan(- 2.0 / 180.0 * TMath::Pi()) ) / 0.016 );
1415 }
1416
1417 Bool_t AliTRDgtuParam::GenerateRecoCoefficients(Int_t trackletMask)
1418 {
1419   // calculate the coefficients for the straight line fit
1420   // depending on the mask of contributing tracklets
1421
1422   fCurrTrackletMask = trackletMask;
1423
1424   TMatrix a(GetNLayers(), 3);
1425   TMatrix b(3, GetNLayers());
1426   TMatrix c(3, 3);
1427
1428   for (Int_t layer = 0; layer < GetNLayers(); layer++) {
1429       if ( (trackletMask & (1 << layer)) == 0) {
1430           a(layer, 0) = 0;
1431           a(layer, 1) = 0;
1432           a(layer, 2) = 0;
1433       }
1434       else {
1435           a(layer, 0) = 1;
1436           a(layer, 1) = fGeo->GetTime0(layer);
1437           a(layer, 2) = (layer % 2 ? 1 : -1) * fGeo->GetTime0(layer);
1438       }
1439   }
1440
1441   b.Transpose(a);
1442   c = b * a;
1443   c.InvertFast();
1444   b = c * b;
1445
1446   for (Int_t layer = 0; layer < GetNLayers(); layer++) {
1447       fAki[layer] = b.GetMatrixArray()[layer];
1448       fBki[layer] = b.GetMatrixArray()[GetNLayers() + layer];
1449       fCki[layer] = b.GetMatrixArray()[2 * GetNLayers() + layer];
1450     }
1451   return kTRUE;
1452 }
1453
1454 Float_t AliTRDgtuParam::GetAki(Int_t k, Int_t i)
1455 {
1456   // get A_ki for the calculation of the tracking parameters
1457   if (fCurrTrackletMask != k)
1458     GenerateRecoCoefficients(k);
1459
1460   return -fAki[i];
1461 }
1462
1463 Float_t AliTRDgtuParam::GetBki(Int_t k, Int_t i)
1464 {
1465   // get B_ki for the calculation of the tracking parameters
1466
1467   if (fCurrTrackletMask != k)
1468     GenerateRecoCoefficients(k);
1469
1470   return fBki[i];
1471 }
1472
1473 Float_t AliTRDgtuParam::GetCki(Int_t k, Int_t i)
1474 {
1475   // get B_ki for the calculation of the tracking parameters
1476
1477   if (fCurrTrackletMask != k)
1478     GenerateRecoCoefficients(k);
1479
1480   return fCki[i];
1481 }
1482
1483 /*
1484 Float_t AliTRDgtuParam::GetD(Int_t k) const
1485 {
1486   // get the determinant for the calculation of the tracking parameters
1487
1488   TMatrix t(3, 3);
1489   for (Int_t i = 0; i < GetNLayers(); i++) {
1490     if ( !((k >> i) & 0x1) )
1491       continue;
1492     Float_t xi = fGeo->GetTime0(i);
1493     t(0,0) += 1;
1494     t(1,0) += xi;
1495     t(2,0) += TMath::Power(-1, i) * xi;
1496     t(0,1) += xi;
1497     t(1,1) += TMath::Power(xi, 2);
1498     t(2,1) += TMath::Power(-1, i) * TMath::Power(xi, 2);
1499     t(0,2) += TMath::Power(-1, i) * xi;
1500     t(1,2) += TMath::Power(-1, i) * TMath::Power(xi, 2);
1501     t(2,2) += TMath::Power(xi, 2);
1502   }
1503   return t.Determinant();
1504 }
1505
1506 Bool_t AliTRDgtuParam::GetFitParams(TVectorD& rhs, Int_t k)
1507 {
1508   // calculate the fitting parameters
1509   // will be changed!
1510
1511   TMatrix t(3,3);
1512   for (Int_t i = 0; i < GetNLayers(); i++) {
1513     if ( !((k >> i) & 0x1) )
1514       continue;
1515     Float_t xi = fGeo->GetTime0(i);
1516     t(0,0) += 1;
1517     t(1,0) += xi;
1518     t(2,0) += TMath::Power(-1, i) * xi;
1519     t(0,1) += xi;
1520     t(1,1) += TMath::Power(xi, 2);
1521     t(2,1) += TMath::Power(-1, i) * TMath::Power(xi, 2);
1522     t(0,2) -= TMath::Power(-1, i) * xi;
1523     t(1,2) -= TMath::Power(-1, i) * TMath::Power(xi, 2);
1524     t(2,2) -= TMath::Power(xi, 2);
1525   }
1526   TDecompLU lr(t);
1527   lr.Solve(rhs);
1528   return lr.Decompose();
1529 }
1530 */
1531
1532 Bool_t AliTRDgtuParam::GetIntersectionPoints(Int_t k, Float_t &x1, Float_t &x2)
1533 {
1534   // get the x-coord. of the assumed circle/straight line intersection points
1535
1536   Int_t l1 = -1;
1537   Int_t l2 = -1;
1538   Int_t nHits = 0;
1539   for (Int_t layer = 0; layer < GetNLayers(); layer++) {
1540     if ( (k >> layer) & 0x1 ) {
1541       if (l1 < 0)
1542         l1 = layer;
1543       l2 = layer;
1544       nHits++;
1545     }
1546   }
1547
1548   if ( (l1 >= 0) && (l2 >= 0) ) {
1549     x1 = fGeo->GetTime0(l1) + 10./6 * (nHits -1);
1550     x2 = fGeo->GetTime0(l2) - 10./6 * (nHits -1);
1551     return kTRUE;
1552   }
1553   else
1554     return kFALSE;
1555 }
1556
1557 Int_t AliTRDgtuParam::GetPt(Int_t layerMask, Int_t a, Float_t /* b */, Float_t x1, Float_t x2, Float_t magField)
1558 {
1559   // returns 0.3 * B * 1/a (1/128 GeV/c)
1560   // a : offset, b : slope (not used)
1561
1562   if (fgUseGTUconst) {
1563     //----- calculation as in the GTU ----
1564     const Int_t maskIdLut[64] = {
1565       -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1,  0,
1566       -1, -1, -1, -1, -1, -1, -1,  1, -1, -1, -1,  2, -1,  3,  4,  5,
1567       -1, -1, -1, -1, -1, -1, -1,  6, -1, -1, -1,  7, -1,  8,  9, 10,
1568       -1, -1, -1, 11, -1, 12, 13, 14, -1, 15, 16, 17, 18, 19, 20, 21
1569     };
1570
1571     const Int_t c1Lut[32] = {
1572       -2371, -2474, -2474, -2474, -2563, -2448, -2578, -2578,
1573       -2578, -2670, -2557, -2578, -2578, -2670, -2557, -2578,
1574       -2670, -2557, -2763, -2557, -2644, -2523,    -1,    -1,
1575       -1,    -1,    -1,    -1,    -1,    -1,    -1,    -1
1576     };
1577
1578     Int_t layerMaskId = maskIdLut[layerMask];
1579     Int_t c1 = c1Lut[layerMaskId];
1580     Int_t c1Ext = c1 << 8;
1581     Int_t ptRawStage4 = c1Ext / (a >> 2);
1582     Int_t ptRawComb4 = ptRawStage4;
1583     Int_t ptExtComb4 = (ptRawComb4 > 0) ? ptRawComb4 + 33 : ptRawComb4 - 30;
1584
1585     return ((Int_t) ptExtComb4/2);
1586   }
1587   else {
1588     //----- simple calculation -----
1589     Float_t c1 = x1 * x2 / 2. / 10000.; // conversion cm to m
1590     Float_t r = 0;
1591     if ( (a >> 1) != 0)
1592       r = (0.3 * magField / 2. / (fgkBinWidthY/100.)) * (((Int_t) c1) << 8) / (a >> 1); //??? why shift of a?
1593
1594     Int_t pt = (Int_t) (2 * r);
1595     if (pt >= 0)
1596       pt += 32;
1597     else
1598       pt -= 29;
1599     pt /= 2;
1600     return pt;
1601   }
1602 }