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