Fixing coding conventions (const arguments in functions)
[u/mrichter/AliRoot.git] / MUON / AliMUONLocalTriggerBoard.cxx
CommitLineData
be972727 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$ */
17
3d1463c8 18//-----------------------------------------------------------------------------
71a2d3aa 19/// \class AliMUONLocalTriggerBoard
20/// A local trigger board has as input a bit pattern and returns
21/// the local trigger response after comparison w/ a LUT
22/// \todo Change member functions comments in capital letters to normal text
23///
24/// \author Rachid Guernane (LPCCFd)
3d1463c8 25//-----------------------------------------------------------------------------
be972727 26
b09247a2 27#include <cstdlib>
be972727 28#include "AliMUONLocalTriggerBoard.h"
29#include "AliMUONTriggerLut.h"
f4678817 30
be972727 31#include "AliLog.h"
32
33#include <TBits.h>
be972727 34#include <Riostream.h>
35
71a2d3aa 36/// \cond CLASSIMP
37ClassImp(AliMUONLocalTriggerBoard)
38/// \endcond
39
e44235aa 40const Int_t AliMUONLocalTriggerBoard::fgkCircuitId[234] =
41{
42 111, 121, 131, 141, 151, 161, 171,
43 211, 212, 221, 222, 231, 232, 241, 242, 251, 252, 261, 262, 271,
44 311, 312, 321, 322, 331, 332, 341, 342, 351, 352, 361, 362, 371,
45 411, 412, 413, 421, 422, 423, 424, 431, 432, 433, 434, 441, 442, 451, 452, 461, 462, 471,
46 521, 522, 523, 524, 531, 532, 533, 534, 541, 542, 551, 552, 561, 562, 571,
47 611, 612, 613, 621, 622, 623, 624, 631, 632, 633, 634, 641, 642, 651, 652, 661, 662, 671,
48 711, 712, 721, 722, 731, 732, 741, 742, 751, 752, 761, 762, 771,
49 811, 812, 821, 822, 831, 832, 841, 842, 851, 852, 861, 862, 871,
50 911, 921, 931, 941, 951, 961, 971,
51 -111, -121, -131, -141, -151, -161, -171,
52 -211, -212, -221, -222, -231, -232, -241, -242, -251, -252, -261, -262, -271,
53 -311, -312, -321, -322, -331, -332, -341, -342, -351, -352, -361, -362, -371,
54 -411, -412, -413, -421, -422, -423, -424, -431, -432, -433, -434, -441, -442, -451, -452, -461, -462, -471,
55 -521, -522, -523, -524, -531, -532, -533, -534, -541, -542, -551, -552, -561, -562, -571,
56 -611, -612, -613, -621, -622, -623, -624, -631, -632, -633, -634, -641, -642, -651, -652, -661, -662, -671,
57 -711, -712, -721, -722, -731, -732, -741, -742, -751, -752, -761, -762, -771,
58 -811, -812, -821, -822, -831, -832, -841, -842, -851, -852, -861, -862, -871,
59 -911, -921, -931, -941, -951, -961, -971
60};
61
be972727 62//___________________________________________
63AliMUONLocalTriggerBoard::AliMUONLocalTriggerBoard()
1657f946 64 : AliMUONTriggerBoard(),
ac809573 65 fMpLocalBoard(0x0),
1657f946 66 fStripX11(0),
8aac6cce 67 fStripY11(15),
1657f946 68 fDev(0),
8aac6cce 69 fTrigY(1),
1657f946 70 fOutput(0),
e93b11d0 71 fLUT(0x0),
ac809573 72 fCoinc44(0)
be972727 73{
ac809573 74/// default constructor
71a2d3aa 75///
ac809573 76
be972727 77 for (Int_t i=0; i<2; i++)
78 for (Int_t j=0; j<4; j++)
79 {
80 fXY[i][j] = fXYU[i][j] = fXYD[i][j] = 0;
81
82 fMask[i][j] = 0xFFFF;
83 }
84
be972727 85 for (Int_t i=0; i<5; i++) fMinDevStrip[i] = fMinDev[i] = fCoordY[i] = 0;
86
c720a83c 87 for (Int_t i=0; i<2; i++) fLutLpt[i] = fLutHpt[i] = 0;
be972727 88}
89
90//___________________________________________
ac809573 91AliMUONLocalTriggerBoard::AliMUONLocalTriggerBoard(AliMpLocalBoard* mpLocalBoard)
92 : AliMUONTriggerBoard(mpLocalBoard->GetName(), mpLocalBoard->GetSlot()),
93 fMpLocalBoard(mpLocalBoard),
1657f946 94 fStripX11(0),
8aac6cce 95 fStripY11(15),
1657f946 96 fDev(0),
8aac6cce 97 fTrigY(1),
1657f946 98 fOutput(0),
ac809573 99 fLUT(0x0),
e93b11d0 100 fCoinc44(0)
be972727 101{
71a2d3aa 102/// Standard constructor
103///
be972727 104
105 for (Int_t i=0; i<2; i++)
106 for (Int_t j=0; j<4; j++)
107 {
108 fXY[i][j] = fXYU[i][j] = fXYD[i][j] = 0;
109
110 fMask[i][j] = 0xFFFF;
111 }
112
be972727 113 for (Int_t i=0; i<5; i++) fMinDevStrip[i] = fMinDev[i] = fCoordY[i] = 0;
114
c720a83c 115 for (Int_t i=0; i<2; i++) fLutLpt[i] = fLutHpt[i] = 0;
be972727 116}
117
71a2d3aa 118//___________________________________________
00d46f24 119AliMUONLocalTriggerBoard::AliMUONLocalTriggerBoard(const AliMUONLocalTriggerBoard& right) :
120 AliMUONTriggerBoard(right),
121 fMpLocalBoard(right.fMpLocalBoard),
122 fStripX11(right.fStripX11),
123 fStripY11(right.fStripY11),
124 fDev(right.fDev),
125 fTrigY(right.fTrigY),
126 fOutput(right.fOutput),
127 fLUT(right.fLUT),
128 fCoinc44(right.fCoinc44)
129{
130 //
131 /// Copy constructor
132 //
133 for (Int_t i=0; i<2; i++) {
134 for (Int_t j=0; j<4; j++) {
135 fXY[i][j] = right.fXY[i][j];
136 fXYU[i][j] = right.fXYU[i][j];
137 fXYD[i][j] = right.fXYD[i][j];
138
139 fMask[i][j] = right.fMask[i][j];
140 }
141 }
142
143 for (Int_t i=0; i<2; i++) {
144 fLutLpt[i] = right.fLutLpt[i];
145 fLutHpt[i] = right.fLutHpt[i];
146 }
147
148 for (Int_t i=0; i<5; i++) {
149 fMinDevStrip[i] = right.fMinDevStrip[i];
150 fMinDev[i] = right.fMinDev[i];
151 fCoordY[i] = right.fCoordY[i];
152 }
153}
154
155//___________________________________________
156AliMUONLocalTriggerBoard& AliMUONLocalTriggerBoard::operator=(const AliMUONLocalTriggerBoard& right)
157{
158/// Assigment operator;
159/// equal operator
160
161 if (this == &right)
162 return *this;
163
164 // base class assignement
165 AliMUONTriggerBoard::operator=(right);
166
167 fMpLocalBoard = right.fMpLocalBoard;
168 fStripX11 = right.fStripX11;
169 fStripY11 = right.fStripY11;
170 fDev = right.fDev;
171 fTrigY = right.fTrigY;
172 fOutput = right.fOutput;
173 fLUT = right.fLUT;
174 fCoinc44 = right.fCoinc44;
175
176 for (Int_t i=0; i<2; i++) {
177 for (Int_t j=0; j<4; j++) {
178 fXY[i][j] = right.fXY[i][j];
179 fXYU[i][j] = right.fXYU[i][j];
180 fXYD[i][j] = right.fXYD[i][j];
181
182 fMask[i][j] = right.fMask[i][j];
183 }
184 }
185
186 for (Int_t i=0; i<2; i++) {
187 fLutLpt[i] = right.fLutLpt[i];
188 fLutHpt[i] = right.fLutHpt[i];
189 }
190
191 for (Int_t i=0; i<5; i++) {
192 fMinDevStrip[i] = right.fMinDevStrip[i];
193 fMinDev[i] = right.fMinDev[i];
194 fCoordY[i] = right.fCoordY[i];
195 }
196
197 return *this;
198}
199
200//___________________________________________
71a2d3aa 201AliMUONLocalTriggerBoard::~AliMUONLocalTriggerBoard()
884a73f1 202{
71a2d3aa 203/// Destructor
204}
884a73f1 205
ac809573 206
207//___________________________________________
208Int_t AliMUONLocalTriggerBoard::GetNumber() const
209{
210/// return board number for notified boards
211
212 if (fMpLocalBoard->IsNotified())
213 return fMpLocalBoard->GetId();
214 else
215 return 0;
216}
217
be972727 218//___________________________________________
219void AliMUONLocalTriggerBoard::Reset()
220{
71a2d3aa 221/// reset board
222///
be972727 223 for (Int_t i=0; i<2; i++)
224 for (Int_t j=0; j<4; j++)
225 fXY[i][j] = fXYU[i][j] = fXYD[i][j] = 0;
226
00d46f24 227 ResetResponse();
228}
be972727 229
00d46f24 230//___________________________________________
231void AliMUONLocalTriggerBoard::ResetResponse()
232{
233/// reset board response
234//
235 fResponse = 0;
be972727 236
00d46f24 237 for (Int_t i=0; i<5; i++) fMinDevStrip[i] = fMinDev[i] = fCoordY[i] = 0;
238
239 fOutput = 0;
be972727 240
00d46f24 241 fStripX11 = 0;
242 fStripY11 = 15;
243 fDev = 0;
244 fTrigY = 1;
be972727 245
00d46f24 246 for (Int_t i=0; i<2; i++) fLutLpt[i] = fLutHpt[i] = 0;
be972727 247}
248
00d46f24 249
be972727 250//___________________________________________
251void AliMUONLocalTriggerBoard::Setbit(Int_t strip, Int_t cathode, Int_t chamber)
252{
71a2d3aa 253/// 0 .. LBS : N-1 .. MSB
be972727 254 TBits w, m;
255
256 UShort_t xy = fXY[cathode][chamber], mask = fMask[cathode][chamber];
257
258 w.Set(16,&xy);
259 m.Set(16,&mask);
260
261 Int_t s = strip - int(strip / 16) * 16;
262
263 w.SetBitNumber(s);
264
265 w &= m;
266
267 UShort_t value;
268
269 w.Get(&value);
270
271 fXY[cathode][chamber] = value;
272}
273
274//___________________________________________
8af755f4 275void AliMUONLocalTriggerBoard::SetbitM(Int_t strip, Int_t cathode, Int_t chamber)
276{
71a2d3aa 277/// 0 .. LBS : N-1 .. MSB
8af755f4 278 TBits w, m;
279
280 UShort_t xy = fXY[cathode][chamber], mask = fMask[cathode][chamber];
281
282 w.Set(16,&xy);
283 m.Set(16,&mask);
284
8af755f4 285 w.SetBitNumber(strip);
286
287 w &= m;
288
289 UShort_t value;
290
291 w.Get(&value);
292
293 fXY[cathode][chamber] = value;
294}
295
ac809573 296
8af755f4 297//___________________________________________
57e2ad1a 298void AliMUONLocalTriggerBoard::Pattern(const Option_t *option) const
be972727 299{
71a2d3aa 300/// print bit pattern
301///
be972727 302 TString op = option;
303
304 if (op.Contains("X")) BP("X");
305
306 if (op.Contains("Y")) BP("Y");
307}
308
309
310//___________________________________________
57e2ad1a 311void AliMUONLocalTriggerBoard::BP(const Option_t *option) const
be972727 312{
b3a9876c 313/// Respect the old printout format
e44235aa 314
315 const Int_t kModuleId[126] =
316 {11,12,13,14,15,16,17, // right side of the chamber
317 21,22,23,24,25,26,27,
318 31,32,33,34,35,36,37,
319 41,42,43,44,45,46,47,
320 51,52,53,54,55,56,57,
321 61,62,63,64,65,66,67,
322 71,72,73,74,75,76,77,
323 81,82,83,84,85,86,87,
324 91,92,93,94,95,96,97,
325 -11,-12,-13,-14,-15,-16,-17, // right side of chamber
326 -21,-22,-23,-24,-25,-26,-27,
327 -31,-32,-33,-34,-35,-36,-37,
328 -41,-42,-43,-44,-45,-46,-47,
329 -51,-52,-53,-54,-55,-56,-57,
330 -61,-62,-63,-64,-65,-66,-67,
331 -71,-72,-73,-74,-75,-76,-77,
332 -81,-82,-83,-84,-85,-86,-87,
333 -91,-92,-93,-94,-95,-96,-97};
334
335 const Int_t kNstripY[126]=
336 { 8, 8, 8, 8, 8, 8,16, // right side of the chamber
337 8, 8, 8, 8, 8, 8,16,
338 16,16,16,16,16, 8,16,
339 16,16,16,16,16, 8,16,
340 0, 8,16,16,16, 8,16,
341 16,16,16,16,16, 8,16,
342 16,16,16,16,16, 8,16,
343 8, 8, 8, 8, 8, 8,16,
344 8, 8, 8, 8, 8, 8,16,
345 8, 8, 8, 8, 8, 8,16, // left side of the chamber
346 8, 8, 8, 8, 8, 8,16,
347 16,16,16,16,16, 8,16,
348 16,16,16,16,16, 8,16,
349 0, 8,16,16,16, 8,16,
350 16,16,16,16,16, 8,16,
351 16,16,16,16,16, 8,16,
352 8, 8, 8, 8, 8, 8,16,
353 8, 8, 8, 8, 8, 8,16};
be972727 354
355 TString op = option;
e44235aa 356
357 TString nn = GetName();
8af755f4 358
be972727 359 if (op.Contains("X"))
360 {
361 printf("-------- TRIGGER INPUT ---------\n");
75f27f52 362 printf("--- warning: switchs not activated at this level ---\n");
be972727 363 printf("===============================================================\n");
364 printf(" 5432109876543210");
365
a6e0ebfe 366 const char *x[4] = {"XMC11","XMC12","XMC21","XMC22"};
367 const char *s[4] = {" ",
be972727 368 " ",
369 " ",
370 " "};
371
372 for (Int_t ch=0; ch<4; ch++)
373 {
374 printf("\n %s%s", x[ch], s[ch]);
375
376 UShort_t xy = fXY[0][ch];
377
378 TBits w(16); w.Set(16,&xy);
379
380 if (ch<2) cout << w;
381 else
382 {
383 UShort_t xyd = fXYD[0][ch], xyu = fXYU[0][ch];
384 TBits dw(16), uw(16); dw.Set(16,&xyd); uw.Set(16,&xyu);
385
386 TBits ew(32);
387
388 for (Int_t i=0;i<16;i++) ew[i+8] = w[i];
389
390 for (Int_t i=0;i<8;i++)
391 {
392 ew[i] = dw[i+8]; // 8 MSB
393 ew[i+24] = uw[i]; //
394 }
395
396 cout << ew;
397 }
398 }
399
400 printf("\n ");
401 printf("10987654321098765432109876543210\n");
402 }
403
404 if (op.Contains("Y"))
405 {
406 printf("---------------------------------------------------------------\n");
407 printf(" ");
408
409/* OLD NUMBERING STYLE */
410/**/
8af755f4 411 Int_t idCircuit = 0, absidModule = 0;
412
e44235aa 413 if (!(nn.Contains("Int")))
414 {
415 idCircuit = fgkCircuitId[GetI()];
416 absidModule = TMath::Abs(Int_t(idCircuit/10));
417 }
8af755f4 418
be972727 419 Int_t iModule=0;
8af755f4 420
be972727 421 for (Int_t i=0; i<63; i++)
422 {
e44235aa 423 if (kModuleId[i]==absidModule)
be972727 424 {
425 iModule=i;
426 break;
427 }
428 }
429
e44235aa 430 Int_t nStrip = kNstripY[iModule];
be972727 431 for (Int_t istrip=nStrip-1; istrip>=0; istrip--) {
432 if (istrip>9) printf("%i",istrip-10*Int_t(istrip/10));
433 if (istrip<10) printf("%i",istrip);
434 }
435/**/
436/* */
437
438 UShort_t xyval = 0;
439
ac809573 440 if (GetSwitch(1))
be972727 441 {
442 xyval = fXY[1][0];
443 TBits v11(8); v11.Set(8,&xyval);
444 printf("\n YMC11 ");
445 cout << v11;
446
447 xyval = fXY[1][1];
448 TBits v12(8); v12.Set(8,&xyval);
449 printf("\n YMC12 ");
450 cout << v12;
451
452 xyval = fXY[1][2];
453 TBits v21(8); v21.Set(8,&xyval);
454 printf("\n YMC21 ");
455 cout << v21;
456
457 xyval = fXY[1][3];
458 TBits v22(8); v22.Set(8,&xyval);
459 printf("\n YMC22 ");
460 cout << v22 << endl;
461 }
462 else
463 {
464 xyval = fXY[1][0];
465 TBits v11(16); v11.Set(16,&xyval);
466 printf("\n YMC11 ");
467 cout << v11;
468
469 xyval = fXY[1][1];
470 TBits v12(16); v12.Set(16,&xyval);
471 printf("\n YMC12 ");
472 cout << v12;
473
474 xyval = fXY[1][2];
475 TBits v21(16); v21.Set(16,&xyval);
476 printf("\n YMC21 ");
477 cout << v21;
478
479 xyval = fXY[1][3];
480 TBits v22(16); v22.Set(16,&xyval);
481 printf("\n YMC22 ");
482 cout << v22 << endl;
483 }
484
be972727 485// tmp
486 printf("---------------------------------------------------------------");
487 printf("\n upper part of circuit %i",idCircuit);
488 printf("\n UMC21 ");
489 xyval = fXYU[1][2];
490 TBits wu21(16); wu21.Set(16,&xyval);
491 cout << wu21;
492 printf("\n UMC22 ");
493 xyval = fXYU[1][3];
494 TBits wu22(16); wu22.Set(16,&xyval);
495 cout << wu22;
496 printf("\n lower part of circuit %i",idCircuit);
497 printf("\n LMC21 ");
498 xyval = fXYD[1][2];
499 TBits wl21(16); wl21.Set(16,&xyval);
500 cout << wl21;
501 printf("\n LMC22 ");
502 xyval = fXYD[1][3];
503 TBits wl22(16); wl22.Set(16,&xyval);
504 cout << wl22;
505 printf("\n");
506 printf("===============================================================\n");
507 }
508}
509
510//___________________________________________
be972727 511void AliMUONLocalTriggerBoard::Module(char *mod)
512{
71a2d3aa 513/// get module from name
514///
f4678817 515 const Int_t kMaxfields = 2; char **fields = new char*[kMaxfields];
be972727 516
517 char s[100]; strcpy(s, GetName());
518
519 Int_t numlines = 0;
520
521 for (char *token = strtok(s, "B");
522 token != NULL;
523 token = strtok(NULL, " "))
524 {
525 fields[numlines] = new char[strlen(token)+1];
526 strcpy(fields[numlines++],token);
527 }
528
529 strcpy(mod,fields[0]);
530}
531
532//___________________________________________
e93b11d0 533void AliMUONLocalTriggerBoard::TrigX(Int_t ch1q[16], Int_t ch2q[16], Int_t ch3q[32], Int_t ch4q[32])
be972727 534{
71a2d3aa 535/// note : coinc44 = flag 0 or 1 (0 coincidence -> 3/4, 1 coincidence -> 4/4) \n
536///--------------------------------------------------------- \n
537/// step # 1 : declustering, reduction DS, calculate sgle & dble \n
538///---------------------------------------------------------
be972727 539 Int_t ch1e[19], ch2e[20], ch3e[35], ch4e[36];
540 Int_t sgleHit1[31], sgleHit2[63];
541 Int_t dbleHit1[31], dbleHit2[63];
542
543 Int_t i;
544 Int_t j;
545 Int_t istrip;
546
547 for (i=0; i<31; i++) {
548 sgleHit1[i]=0;
549 dbleHit1[i]=0;
550 }
551 for (i=0; i<63; i++) {
552 sgleHit2[i]=0;
553 dbleHit2[i]=0;
554 }
555
556//--- inititialize che using chq
75f27f52 557//--- switch zero_down, zero_middle & zero_up added 30/05/07
558//--- fSwitch[7/8/9] = zero_down/zero_middle/zero_up
be972727 559 for (i=0; i<19; i++) {
560 if (i<1||i>16) ch1e[i]=0;
561 else ch1e[i]=ch1q[i-1];
562 }
563 for (i=0; i<20; i++) {
564 if (i<2||i>17) ch2e[i]=0;
565 else ch2e[i]=ch2q[i-2];
566 }
567 for (i=0; i<35; i++) {
568 if (i<1||i>32) ch3e[i]=0;
ac809573 569 else if (i>=1 && i<=8) ch3e[i]=ch3q[i-1]&!GetSwitch(7);
570 else if (i>=9 && i<=24) ch3e[i]=ch3q[i-1]&!GetSwitch(8);
571 else if (i>=25 && i<=32) ch3e[i]=ch3q[i-1]&!GetSwitch(9);
be972727 572 }
573 for (i=0; i<36; i++) {
574 if (i<2||i>33) ch4e[i]=0;
ac809573 575 else if (i>=2 && i<=9) ch4e[i]=ch4q[i-2]&!GetSwitch(7);
576 else if (i>=10 && i<=25) ch4e[i]=ch4q[i-2]&!GetSwitch(8);
577 else if (i>=26 && i<=33) ch4e[i]=ch4q[i-2]&!GetSwitch(9);
be972727 578 }
579
be972727 580//--- calculate dble & sgle first station
581 for (i=0; i<=15; i++) {
582 sgleHit1[2*i] = (!ch1e[i+1]|(ch1e[i]^ch1e[i+2])) &
583 (!ch2e[i+2] | (ch2e[i+1]^ch2e[i+3]));
584
585 dbleHit1[2*i] = ch1e[i+1]&!(ch1e[i+2]^ch1e[i]) &
586 (ch2e[i+2] | (!ch2e[i]&ch2e[i+1]) | (ch2e[i+3]&!ch2e[i+4]));
587 }
588
589 for (i=0; i<=14; i++) {
590 sgleHit1[2*i+1] = (!ch1e[i+1]|!ch1e[i+2]|(ch1e[i]^ch1e[i+3])) &
591 (!ch2e[i+2] | !ch2e[i+3] | (ch2e[i+1]^ch2e[i+4]));
592 dbleHit1[2*i+1] = ch1e[i+1]&ch1e[i+2]&!(ch1e[i]^ch1e[i+3]) &
a6973cb8 593 ((ch2e[i+2]&(!ch2e[i+1]|!ch2e[i])) |
b88fed65 594 (ch2e[i+3]&(ch2e[i+2]|!ch2e[i+4]|!ch2e[i+5])));
be972727 595 }
596
597//--- calculate dble & sgle second station
598 for (i=0; i<=31; i++) {
599 sgleHit2[2*i] = (!ch3e[i+1]|(ch3e[i]^ch3e[i+2])) &
600 (!ch4e[i+2] | (ch4e[i+1]^ch4e[i+3]));
601 dbleHit2[2*i] = ch3e[i+1]&!(ch3e[i+2]^ch3e[i]) &
602 (ch4e[i+2] | (!ch4e[i]&ch4e[i+1]) | (ch4e[i+3]&!ch4e[i+4]));
603 }
604
605 for (i=0; i<=30; i++) {
606 sgleHit2[2*i+1] = (!ch3e[i+1]|!ch3e[i+2]|(ch3e[i]^ch3e[i+3])) &
607 (!ch4e[i+2] | !ch4e[i+3] | (ch4e[i+1]^ch4e[i+4]));
608 dbleHit2[2*i+1] = ch3e[i+1]&ch3e[i+2]&!(ch3e[i]^ch3e[i+3]) &
a6973cb8 609 ((ch4e[i+2]&(!ch4e[i+1]|!ch4e[i])) |
b88fed65 610 (ch4e[i+3]&((ch4e[i+2]|!ch4e[i+4])|!ch4e[i+5])));
be972727 611 }
612
613//---
614 if(AliDebugLevel()==3||AliDebugLevel()==5) {
615 printf("===============================================================\n");
616 printf(" X plane after sgle and dble \n");
617 printf(" 0987654321098765432109876543210");
618 printf("\n SGLE1 ");
619 for (istrip=30; istrip>=0; istrip--) printf("%i",(!sgleHit1[istrip]));
620 printf("\n DBLE1 ");
621 for (istrip=30; istrip>=0; istrip--) printf("%i",dbleHit1[istrip]);
622 printf("\n SGLE2 ");
623 for (istrip=62; istrip>=0; istrip--) printf("%i",(!sgleHit2[istrip]));
624 printf("\n DBLE2 ");
625 for (istrip=62; istrip>=0; istrip--) printf("%i",dbleHit2[istrip]);
626 printf("\n 210987654321098765432109876543210987654321098765432109876543210\n");
627 }
628
629//---------------------------------------------------------
630// step # 2 : coincidence 3/4
631//---------------------------------------------------------
632 Int_t rearImage[31][31];
633 for (i=0; i<31; i++) {
634 for (j=0; j<31; j++) {
635 rearImage[i][j]=0;
636 }
637 }
638
639 Int_t notOr1=!dbleHit1[30] & !dbleHit1[29] & !dbleHit1[28] & !dbleHit1[27] &
640 !dbleHit1[26] & !dbleHit1[25] & !dbleHit1[24] & !dbleHit1[23] &
641 !dbleHit1[22] & !dbleHit1[21] & !dbleHit1[20] & !dbleHit1[19] &
642 !dbleHit1[18] & !dbleHit1[17] & !dbleHit1[16] & !dbleHit1[15] &
643 !dbleHit1[14] & !dbleHit1[13] & !dbleHit1[12] & !dbleHit1[11] &
644 !dbleHit1[10] & !dbleHit1[9] & !dbleHit1[8] & !dbleHit1[7] &
645 !dbleHit1[6] & !dbleHit1[5] & !dbleHit1[4] & !dbleHit1[3] &
e93b11d0 646 !dbleHit1[2] & !dbleHit1[1] & !dbleHit1[0] & !fCoinc44;
be972727 647
648 Int_t notOr2= !dbleHit2[62] & !dbleHit2[61] & !dbleHit2[60] & !dbleHit2[59] &
649 !dbleHit2[58] & !dbleHit2[57] & !dbleHit2[56] & !dbleHit2[55] &
650 !dbleHit2[54] & !dbleHit2[53] & !dbleHit2[52] & !dbleHit2[51] &
651 !dbleHit2[50] & !dbleHit2[49] & !dbleHit2[48] & !dbleHit2[47] &
652 !dbleHit2[46] & !dbleHit2[45] & !dbleHit2[44] & !dbleHit2[43] &
653 !dbleHit2[42] & !dbleHit2[41] & !dbleHit2[40] & !dbleHit2[39] &
654 !dbleHit2[38] & !dbleHit2[37] & !dbleHit2[36] & !dbleHit2[35] &
655 !dbleHit2[34] & !dbleHit2[33] & !dbleHit2[32] & !dbleHit2[31] &
656 !dbleHit2[30] & !dbleHit2[29] & !dbleHit2[28] & !dbleHit2[27] &
657 !dbleHit2[26] & !dbleHit2[25] & !dbleHit2[24] & !dbleHit2[23] &
658 !dbleHit2[22] & !dbleHit2[21] & !dbleHit2[20] & !dbleHit2[19] &
659 !dbleHit2[18] & !dbleHit2[17] & !dbleHit2[16] & !dbleHit2[15] &
660 !dbleHit2[14] & !dbleHit2[13] & !dbleHit2[12] & !dbleHit2[11] &
661 !dbleHit2[10] & !dbleHit2[9] & !dbleHit2[8] & !dbleHit2[7] &
662 !dbleHit2[6] & !dbleHit2[5] & !dbleHit2[4] & !dbleHit2[3] &
e93b11d0 663 !dbleHit2[2] & !dbleHit2[1] & !dbleHit2[0] & !fCoinc44;
be972727 664
665// DS reduction
666 for (i=0; i<31; i++) {
667 sgleHit1[i] = !sgleHit1[i]&notOr1;
668 }
669 for (i=0; i<63; i++) {
670 sgleHit2[i] = !sgleHit2[i]&notOr2;
671 }
672
673// extract rearImage
674 for (i=0; i<31; i++){
675 Int_t tmpSgleHit2[31];
676 Int_t tmpDbleHit2[31];
677 for (j=0; j<31; j++){
678 tmpSgleHit2[j] = sgleHit2[i+j+1];
679 tmpDbleHit2[j] = dbleHit2[i+j+1];
680 }
681
682 for (Int_t k=0; k<31; k++) {
683 rearImage[i][k]=(sgleHit1[i]&tmpDbleHit2[k])|
684 (dbleHit1[i]&(tmpSgleHit2[k]|tmpDbleHit2[k]));
685 }
686 }
687
688 //-----------
689 if(AliDebugLevel()==3||AliDebugLevel()==5) {
690 printf("===============================================================\n");
691 for (i=30; i>=0; i--) {
692 printf("%i \t",i);
693 for (istrip=31; istrip>=0; istrip--) printf("%i",rearImage[i][istrip]);
694 printf("\n");
695 }
696 }
697
698//---------------------------------------------------------
699// step # 3 : calculate deviation
700//---------------------------------------------------------
701 Int_t dev[31][6];
702 for (i=0; i<31; i++) {
703 for (j=0; j<6; j++) {
704 dev[i][j]=0;
705 }
706 }
707
708 for (i=0; i<31; i++){
709 Int_t leftDev[5], rightDev[5];
710 Int_t orL1, andL1, andL2, orR1, orR2, andR1, andR2, andR3;
711
712// calculate Left deviation
713 orL1=rearImage[i][16]|rearImage[i][18]|rearImage[i][20]|rearImage[i][22];
714 andL1=!rearImage[i][17]&!rearImage[i][19]&!rearImage[i][21] & !orL1;
715 andL2=!rearImage[i][23]&!rearImage[i][24]&!rearImage[i][25]&!rearImage[i][26];
716
717 leftDev[0] = (rearImage[i][16]|!rearImage[i][17]) &
b88fed65 718 (rearImage[i][16]|rearImage[i][18]|(!rearImage[i][19]&
719 (rearImage[i][20]|!rearImage[i][21]))) &
720 (orL1|(!rearImage[i][23]&(rearImage[i][24]|!rearImage[i][25]))) &
721 (orL1|rearImage[i][24]|rearImage[i][26]|(!rearImage[i][27]&
722 (rearImage[i][28]|!rearImage[i][29])));
be972727 723
724 leftDev[1] = !rearImage[i][16] &
725 !(!rearImage[i][17]&!rearImage[i][18]&!rearImage[i][21]&!rearImage[i][22] &
b88fed65 726 (!rearImage[i][25]&!rearImage[i][26]&(rearImage[i][27]|rearImage[i][28]))) &
727 (rearImage[i][17]|rearImage[i][18] | (!rearImage[i][19]&!rearImage[i][20])) &
be972727 728 (rearImage[i][17]|rearImage[i][18]|rearImage[i][21]|rearImage[i][22] |
b88fed65 729 (!rearImage[i][23]&!rearImage[i][24]));
be972727 730
731 leftDev[2] = (!rearImage[i][16]&!rearImage[i][17]&!rearImage[i][18]) &
732 (rearImage[i][19]|rearImage[i][20]|rearImage[i][21]|rearImage[i][22] | andL2);
733
734 leftDev[3] = andL1;
735
736 leftDev[4] =
737 !rearImage[i][27]&!rearImage[i][28]&!rearImage[i][29]&!rearImage[i][30] &
738 andL1 & andL2;
739
740 // calculate Right deviation
741 orR1=rearImage[i][8]|rearImage[i][10]|rearImage[i][12]|rearImage[i][14];
742 orR2=rearImage[i][8]|rearImage[i][9]|rearImage[i][10]|rearImage[i][11];
743 andR1=!rearImage[i][12]&!rearImage[i][13]&!rearImage[i][14]&!rearImage[i][15];
744 andR2=
745 !rearImage[i][8]&!rearImage[i][9]&!rearImage[i][10]&!rearImage[i][11] & andR1;
746 andR3=!rearImage[i][4]&!rearImage[i][5]&!rearImage[i][6]&!rearImage[i][7];
747
748 rightDev[0] = !rearImage[i][15]&(rearImage[i][14]|!rearImage[i][13]) &
b88fed65 749 ((rearImage[i][12]|rearImage[i][14]|(!rearImage[i][11]&
750 (rearImage[i][10]|!rearImage[i][9]))) &
751 ((orR1|(!rearImage[i][7]&(rearImage[i][6]|!rearImage[i][5]))) &
752 (orR1|rearImage[i][4]|rearImage[i][6]|(!rearImage[i][3]&(rearImage[i][2]|
753 !rearImage[i][1])))));
be972727 754
755 rightDev[1] = !rearImage[i][15]&!rearImage[i][14] &
756 !(!rearImage[i][4]&!rearImage[i][5]&!rearImage[i][8]&!rearImage[i][9] &
757 (!rearImage[i][12]&!rearImage[i][13]&(rearImage[i][2]|rearImage[i][3]))) &
b88fed65 758 (rearImage[i][12]|rearImage[i][13] | (!rearImage[i][10]&!rearImage[i][11])) &
be972727 759 (rearImage[i][8]|rearImage[i][9]|rearImage[i][12]|rearImage[i][13] |
b88fed65 760 (!rearImage[i][6]&!rearImage[i][7]));
be972727 761
762 rightDev[2] = andR1 & (orR2 | andR3);
763 rightDev[3] = andR2;
764 rightDev[4] =
765 !rearImage[i][0]&!rearImage[i][1]&!rearImage[i][2]&!rearImage[i][3] &
766 andR2 & andR3 ;
767
768 // compare Left & Right deviations
769 Int_t tmpLeftDev=0, tmpRightDev=0;
770 for (j=0; j<5; j++){
771 tmpLeftDev = tmpLeftDev + Int_t(leftDev[j]<<j);
772 tmpRightDev = tmpRightDev + Int_t(rightDev[j]<<j);
773 }
774
775 // assign mimimum deviation do dev[][]
776 if (tmpLeftDev < tmpRightDev ){
777 for (j=0; j<5; j++){ dev[i][j]=leftDev[j];}
778 dev[i][5]=1;
779 } else {
780 for (j=0; j<5; j++){ dev[i][j]=rightDev[j];}
781 dev[i][5]=0;
782 }
783 }
784
785//---
786 if(AliDebugLevel()==3||AliDebugLevel()==5) {
787 printf("===============================================================\n");
788 for (i=30; i>=0; i--) {
789 printf("%i \t",i);
790 for (istrip=5; istrip>=0; istrip--) printf("%i",dev[i][istrip]);
791 printf(" \n");
792 }
793 }
794
795//---------------------------------------------------------
796// step # 4 : sort deviation
797//---------------------------------------------------------
798 Int_t bga1[16], bga2[8], bga3[4], bga4[2], bga5;
799 Int_t tmpbga1[16][6], tmpbga2[8][6], tmpbga3[4][6], tmpbga4[2][6], tmpbga5[6];
800 Int_t tmpMax[6]={1,1,1,1,1,0};
801
802 for (i=0; i<15; i++) {
803 Sort2x5(dev[2*i],dev[2*i+1],tmpbga1[i],bga1[i]);
804 }
805 Sort2x5(dev[30],tmpMax,tmpbga1[15],bga1[15]);
806
807//--
808 if(AliDebugLevel()==3||AliDebugLevel()==5) {
809 printf("===============================================================\n");
810 printf(" sorting : 1st level \n");
811 for (i=15; i>=0; i--) {
812 printf("\t %i \t",bga1[i]);
813 for (j=5; j>=0; j--) printf("%i",tmpbga1[i][j]);
814 printf(" \n");
815 }
816 }
817
818 for (i=0; i<8; i++) {
819 Sort2x5(tmpbga1[2*i],tmpbga1[2*i+1],tmpbga2[i],bga2[i]);
820 }
821
822//--
823 if(AliDebugLevel()==3||AliDebugLevel()==5) {
824 printf("===============================================================\n");
825 printf(" sorting : 2nd level \n");
826 for (i=7; i>=0; i--) {
827 printf("\t %i \t",bga2[i]);
828 for (j=5; j>=0; j--) printf("%i",tmpbga1[i][j]);
829 printf(" \n");
830 }
831 }
832
833 for (i=0; i<4; i++) {
834 Sort2x5(tmpbga2[2*i],tmpbga2[2*i+1],tmpbga3[i],bga3[i]);
835 }
836
837//--
838 if(AliDebugLevel()==3||AliDebugLevel()==5) {
839 printf("===============================================================\n");
840 printf(" sorting : 3rd level \n");
841 for (i=3; i>=0; i--) {
842 printf("\t %i \t",bga3[i]);
843 for (j=5; j>=0; j--) printf("%i",tmpbga3[i][j]);
844 printf(" \n");
845 }
846 }
847
848 for (i=0; i<2; i++) {
849 Sort2x5(tmpbga3[2*i],tmpbga3[2*i+1],tmpbga4[i],bga4[i]);
850 }
851
852//--
853 if(AliDebugLevel()==3||AliDebugLevel()==5) {
854 printf("===============================================================\n");
855 printf(" sorting : 4th level \n");
856 for (i=1; i>=0; i--) {
857 printf("\t %i \t",bga4[i]);
858 for (j=5; j>=0; j--) printf("%i",tmpbga4[i][j]);
859 printf(" \n");
860 }
861 }
862
863 Sort2x5(tmpbga4[0],tmpbga4[1],tmpbga5,bga5);
864
865 // coding from 6 to 5 bits
866 fMinDev[4] = tmpbga5[5] | tmpbga5[4];
867 for (i=0; i<4; i++) {
868 fMinDev[i]=tmpbga5[i] & !tmpbga5[4];
869 }
870
871 // find address of strip with minimum deviation
872 fMinDevStrip[4]=bga5;
873 if (bga5<=1) fMinDevStrip[3]=bga4[bga5];
874
875 Int_t tmpAd=fMinDevStrip[3]+fMinDevStrip[4]*2;
876 if (tmpAd<=3) fMinDevStrip[2]=bga3[tmpAd];
877
878 tmpAd=fMinDevStrip[2]+fMinDevStrip[3]*2+fMinDevStrip[4]*4;
879 if (tmpAd<=7) fMinDevStrip[1]=bga2[tmpAd];
880
881 tmpAd=fMinDevStrip[1]+fMinDevStrip[2]*2+fMinDevStrip[3]*4+fMinDevStrip[4]*8;
882 if (tmpAd<=15) fMinDevStrip[0]=bga1[tmpAd];
883
884 if(AliDebugLevel()==3||AliDebugLevel()==5) {
885 printf("===============================================================\n");
886 printf("minDevStrip = ");
887 for (i=4; i>=0; i--) printf("%i",fMinDevStrip[i]);
888 printf(" minDev = ");
889 for (i=4; i>=0; i--) printf("%i",fMinDev[i]);
890 printf(" \n");
891 printf("===============================================================\n");
892 }
893
894}
895
896//___________________________________________
897void AliMUONLocalTriggerBoard::Sort2x5(Int_t dev1[6], Int_t dev2[6],
898 Int_t minDev[6], Int_t &dev1GTdev2)
899{
71a2d3aa 900/// returns minimun between dev1 and dev2
be972727 901 Int_t tmpDev1=0, tmpDev2=0;
902
903 for (Int_t j=0; j<5; j++)
904 {
905 tmpDev1 += Int_t(dev1[j]<<j);
906 tmpDev2 += Int_t(dev2[j]<<j);
907 }
908
909 if (tmpDev1<=tmpDev2)
910 {
911 for (Int_t j=0; j<=5; j++) minDev[j]=dev1[j];
912 dev1GTdev2=0;
913 }
914 else
915 {
916 for (Int_t j=0; j<=5; j++) minDev[j]=dev2[j];
917 dev1GTdev2=1;
918 }
919}
920
921//___________________________________________
922void AliMUONLocalTriggerBoard::TrigY(Int_t y1[16], Int_t y2[16], Int_t y3[16], Int_t y4[16],
e93b11d0 923 Int_t y3u[16], Int_t y3d[16], Int_t y4u[16], Int_t y4d[16])
be972727 924{
71a2d3aa 925/// note : resMid = 1 -> cancel \n
926///---------------------------------------------------------\n
927/// step # 1 : prehandling Y \n
928///---------------------------------------------------------
be972727 929 Int_t i;
930 Int_t istrip;
931
932 for (i=0; i<16; i++)
933 {
ac809573 934 y3[i]=y3[i]&!GetSwitch(8);
935 y4[i]=y4[i]&!GetSwitch(8);
be972727 936 }
937
938// 10/29/04 fZeroAllYLSB added
939// for (i=0; i<8; i++)
940// {
941// y1[i] = y1[i]&!fSwitch[6];
942// y2[i] = y2[i]&!fSwitch[6];
943// y3[i] = y3[i]&!fSwitch[6];
944// y4[i] = y4[i]&!fSwitch[6];
945// }
946
947 Int_t ch1[16], ch2[16], ch3[16], ch4[16];
948
949 Int_t tmpy3to16[16], tmpy4to16[16];
950 Int_t tmpy3uto16[16], tmpy3dto16[16], tmpy4uto16[16], tmpy4dto16[16];
951 for (i=0; i<8; i++){
b88fed65 952 ch1[2*i] = (y1[i]&GetSwitch(1)) | (y1[2*i]&!GetSwitch(1));
953 ch1[2*i+1] = (y1[i]&GetSwitch(1)) | (y1[2*i+1]&!GetSwitch(1));
be972727 954
b88fed65 955 ch2[2*i] = (y2[i]&GetSwitch(1)) | (y2[2*i]&!GetSwitch(1));
956 ch2[2*i+1] = (y2[i]&GetSwitch(1)) | (y2[2*i+1]&!GetSwitch(1));
be972727 957
b88fed65 958 tmpy3to16[2*i ] = (y3[i]&GetSwitch(1)) | (y3[2*i ]&!GetSwitch(1));
959 tmpy3to16[2*i+1] = (y3[i]&GetSwitch(1)) | (y3[2*i+1]&!GetSwitch(1));
be972727 960
b88fed65 961 tmpy4to16[2*i ] = (y4[i]&GetSwitch(1)) | (y4[2*i ]&!GetSwitch(1));
962 tmpy4to16[2*i+1] = (y4[i]&GetSwitch(1)) | (y4[2*i+1]&!GetSwitch(1));
be972727 963
b88fed65 964 tmpy3uto16[2*i ] = (y3u[i]&GetSwitch(2)) | (y3u[2*i ]&!GetSwitch(2));
965 tmpy3uto16[2*i+1] = (y3u[i]&GetSwitch(2)) | (y3u[2*i+1]&!GetSwitch(2));
be972727 966
b88fed65 967 tmpy4uto16[2*i ] = (y4u[i]&GetSwitch(2)) | (y4u[2*i ]&!GetSwitch(2));
968 tmpy4uto16[2*i+1] = (y4u[i]&GetSwitch(2)) | (y4u[2*i+1]&!GetSwitch(2));
be972727 969
b88fed65 970 tmpy3dto16[2*i ] = (y3d[i]&GetSwitch(0)) | (y3d[2*i ]&!GetSwitch(0));
971 tmpy3dto16[2*i+1] = (y3d[i]&GetSwitch(0)) | (y3d[2*i+1]&!GetSwitch(0));
be972727 972
b88fed65 973 tmpy4dto16[2*i ] = (y4d[i]&GetSwitch(0)) | (y4d[2*i ]&!GetSwitch(0));
974 tmpy4dto16[2*i+1] = (y4d[i]&GetSwitch(0)) | (y4d[2*i+1]&!GetSwitch(0));
be972727 975 }
976
ac809573 977 if (GetSwitch(3)==0&&GetSwitch(4)==0){
be972727 978 for (i=0; i<16; i++){
979 ch3[i] = tmpy3to16[i];
980 ch4[i] = tmpy4to16[i];
981 }
982 }
ac809573 983 if (GetSwitch(3)==0&&GetSwitch(4)==1){
be972727 984 for (i=0; i<16; i++){
8af755f4 985 ch3[i] = tmpy3dto16[i]|tmpy3to16[i];
986 ch4[i] = tmpy4dto16[i]|tmpy4to16[i];
be972727 987 }
988 }
ac809573 989 if (GetSwitch(3)==1&&GetSwitch(4)==0){
be972727 990 for (i=0; i<16; i++){
8af755f4 991 ch3[i] = tmpy3uto16[i]|tmpy3to16[i];
992 ch4[i] = tmpy4uto16[i]|tmpy4to16[i];
be972727 993 }
994 }
ac809573 995 if (GetSwitch(3)==1&&GetSwitch(4)==1){
be972727 996 for (i=0; i<16; i++){
997 ch3[i] = tmpy3dto16[i]|tmpy3to16[i]|tmpy3uto16[i];
998 ch4[i] = tmpy4dto16[i]|tmpy4to16[i]|tmpy4uto16[i];
999 }
1000 }
1001
1002// debug
1003 if(AliDebugLevel()==4||AliDebugLevel()==5) {
1004 printf("===============================================================\n");
1005 printf(" Y plane after PreHandling x2m x2u x2d orMud %i %i %i %i %i \n",
ac809573 1006 GetSwitch(1),GetSwitch(2), GetSwitch(0),GetSwitch(3),GetSwitch(4));
be972727 1007 printf(" ");
1008 for (istrip=15; istrip>=0; istrip--) {
1009 if (istrip>9) printf("%i",istrip-10*Int_t(istrip/10));
1010 if (istrip<10) printf("%i",istrip);
1011 }
1012 printf("\n YMC11 ");
1013 for (istrip=15; istrip>=0; istrip--) printf("%i",ch1[istrip]);
1014 printf("\n YMC12 ");
1015 for (istrip=15; istrip>=0; istrip--) printf("%i",ch2[istrip]);
1016 printf("\n YMC21 ");
1017 for (istrip=15; istrip>=0; istrip--) printf("%i",ch3[istrip]);
1018 printf("\n YMC22 ");
1019 for (istrip=15; istrip>=0; istrip--) printf("%i",ch4[istrip]);
1020 printf(" \n");
1021 }
1022//debug
1023
1024//---------------------------------------------------------
1025// step # 2 : calculate sgle and dble, apply DS reduction
1026//---------------------------------------------------------
1027 Int_t sgle1[16], dble1[16];
1028 Int_t sgle2[16], dble2[16];
1029
1030 // Calculate simple and double hits
1031 for (i=0; i<16; i++) {
1032 dble1[i] = ch1[i] & ch2[i];
1033 dble2[i] = ch3[i] & ch4[i];
1034
1035 sgle1[i] = (ch1[i]|ch2[i]);
1036 sgle2[i] = (ch3[i]|ch4[i]);
1037 }
1038
1039 //debug
1040 if(AliDebugLevel()==4||AliDebugLevel()==5) {
1041 printf("===============================================================\n");
1042 printf(" Y plane after sgle dble \n");
1043 printf(" ");
1044 for (istrip=15; istrip>=0; istrip--) {
1045 if (istrip>9) printf("%i",istrip-10*Int_t(istrip/10));
1046 if (istrip<10) printf("%i",istrip);
1047 }
1048 printf("\n SGLE1 ");
1049 for (istrip=15; istrip>=0; istrip--) printf("%i",sgle1[istrip]);
1050 printf("\n DBLE1 ");
1051 for (istrip=15; istrip>=0; istrip--) printf("%i",dble1[istrip]);
1052 printf("\n SGLE2 ");
1053 for (istrip=15; istrip>=0; istrip--) printf("%i",sgle2[istrip]);
1054 printf("\n DBLE2 ");
1055 for (istrip=15; istrip>=0; istrip--) printf("%i",dble2[istrip]);
1056 printf(" \n");
1057 }
1058 //debug
1059
1060 // DS Reduction
1061 Int_t notOr1, notOr2;
1062
1063 notOr1=!dble1[15] & !dble1[14] & !dble1[13] & !dble1[12] &
1064 !dble1[11] & !dble1[10] & !dble1[9] & !dble1[8] &
1065 !dble1[7] & !dble1[6] & !dble1[5] & !dble1[4] &
1066 !dble1[3] & !dble1[2] & !dble1[1] & !dble1[0];
1067
1068 notOr2=!dble2[15] & !dble2[14] & !dble2[13] & !dble2[12] &
1069 !dble2[11] & !dble2[10] & !dble2[9] & !dble2[8] &
1070 !dble2[7] & !dble2[6] & !dble2[5] & !dble2[4] &
1071 !dble2[3] & !dble2[2] & !dble2[1] & !dble2[0];
1072
1073 for (i=0; i<16; i++) {
e93b11d0 1074 sgle1[i] = sgle1[i] & notOr1 & !fCoinc44;
1075 sgle2[i] = sgle2[i] & notOr2 & !fCoinc44;
be972727 1076 }
1077
1078//---------------------------------------------------------
1079// step # 3 : 3/4 coincidence
1080//---------------------------------------------------------
1081 Int_t frontImage[16];
1082
1083 for (i=1; i<15; i++) {
a6973cb8 1084 frontImage[i] = ((dble1[i] | sgle1[i]) &
1085 (dble2[i+1] | dble2[i] | dble2[i-1])) |
b88fed65 1086 (dble1[i] & (sgle2[i+1] | sgle2[i] | sgle2[i-1]));
be972727 1087 }
a6973cb8 1088 frontImage[0] = ((dble1[0] | sgle1[0]) &
1089 (dble2[1] | dble2[0])) | (dble1[0] & (sgle2[1] | sgle2[0]));
be972727 1090
a6973cb8 1091 frontImage[15] = ((dble1[15] | sgle1[15]) &
1092 (dble2[15] | dble2[14])) | (dble1[15] & (sgle2[15] | sgle2[14]));
be972727 1093
1094
1095//debug
1096 if(AliDebugLevel()==4||AliDebugLevel()==5) {
1097 printf("===============================================================\n");
1098 printf(" Y plane frontImage\n");
1099 printf(" ");
1100 for (istrip=15; istrip>=0; istrip--) {
1101 if (istrip>9) printf("%i",istrip-10*Int_t(istrip/10));
1102 if (istrip<10) printf("%i",istrip);
1103 }
1104 printf("\n ");
1105 for (istrip=15; istrip>=0; istrip--) printf("%i",frontImage[istrip]);
1106 printf("\n");
1107 }
1108//debug
1109
1110//---------------------------------------------------------
1111// step # 4 : Y position
1112//---------------------------------------------------------
1113 Int_t or1, or2, and1, and2, and3;
1114
1115 or1 = frontImage[7]|frontImage[5]|frontImage[3]|frontImage[1];
1116 or2 = frontImage[7]|frontImage[6]|frontImage[5]|frontImage[4];
1117 and1 = !frontImage[3]&!frontImage[2]&!frontImage[1]&!frontImage[0];
1118 and2 = !frontImage[7]&!frontImage[6]&!frontImage[5]&!frontImage[4] & and1;
1119 and3 = !frontImage[11]&!frontImage[10]&!frontImage[9]&!frontImage[8];
1120
1121 fCoordY[0] = !frontImage[0]&(frontImage[1]|!frontImage[2]) &
b88fed65 1122 (frontImage[3]|frontImage[1]|(!frontImage[4]&(frontImage[5]|!frontImage[6]))) &
1123 (or1|(!frontImage[8]&(frontImage[9]|!frontImage[10]))) &
1124 (or1|frontImage[11]|frontImage[9]|(!frontImage[12]&(frontImage[13]|!frontImage[14])));
be972727 1125
1126 fCoordY[1] = !frontImage[0]&!frontImage[1] &
1127 !(!frontImage[11]&!frontImage[10]&!frontImage[7]&!frontImage[6] &
1128 !frontImage[3]&!frontImage[2]&(frontImage[13]|frontImage[12])) &
b88fed65 1129 (frontImage[3]|frontImage[2] | (!frontImage[5]&!frontImage[4])) &
be972727 1130 (frontImage[7]|frontImage[6]|frontImage[3]|frontImage[2] |
b88fed65 1131 (!frontImage[9]&!frontImage[8]));
be972727 1132
1133 fCoordY[2] = and1 & (or2 | and3);
1134
1135 fCoordY[3] = and2;
1136
1137 fCoordY[4] = !frontImage[15]&!frontImage[14]&!frontImage[13]&!frontImage[12] &
1138 and2 & and3 ;
1139}
1140
1141//___________________________________________
1142void AliMUONLocalTriggerBoard::LocalTrigger()
1143{
71a2d3aa 1144/// L0 trigger after LUT
1145///
da81e4c2 1146 Int_t deviation=0;
1147 Int_t iStripY=0;
1148 Int_t iStripX=0;
5a8d50e9 1149 Bool_t xOutput = IsTrigX();
1150 Bool_t yOutput = IsTrigY();
be972727 1151
da81e4c2 1152 if (xOutput) {
1153 for (Int_t i=0; i<5; i++) iStripX += static_cast<int>( fMinDevStrip[i] << i );
5a8d50e9 1154 for (Int_t i=0; i<4; i++) deviation += static_cast<int>( fMinDev[i] << i );
be972727 1155 fDev = deviation;
fa8d035f 1156 fStripX11 = iStripX;
da81e4c2 1157 }
1158 if (yOutput) {
5a8d50e9 1159 for (Int_t i=0; i<4; i++) iStripY += static_cast<int>( fCoordY[i] << i );
da81e4c2 1160 fStripY11 = iStripY;
8aac6cce 1161 fTrigY = fCoordY[4];
da81e4c2 1162 }
1163
1164// cout << " Local Trigger " << " " << fNumber << " " <<
1165// xOutput << " " << yOutput << " " << fDev << " " << fStripX11 << " " <<
1166// fTrigY << " " << fStripY11 << "\n";
1167
abe5ad6f 1168 if (xOutput && yOutput){ // trigger in X and Y
da81e4c2 1169 fOutput =1;
1170
1171 Int_t sign = 0;
1172 if ( !fMinDev[4] && deviation ) sign=-1;
1173 if ( !fMinDev[4] && !deviation ) sign= 0;
1174 if ( fMinDev[4] == 1 ) sign=+1;
1175
1176 deviation *= sign;
1177
be972727 1178// calculate deviation in [0;+30]
da81e4c2 1179 deviation += 15;
1180
8af755f4 1181// GET LUT OUTPUT FOR icirc/istripX1/deviation/istripY
ac809573 1182 fLUT->GetLutOutput(GetNumber(), fStripX11, deviation, fStripY11, fLutLpt, fLutHpt);
1183 }
1184
1185 fResponse = fLutLpt[0] +
87280b64 1186 static_cast<int>(fLutLpt[1]<<1) +
1187 static_cast<int>(fLutHpt[0]<<2) +
1188 static_cast<int>(fLutHpt[1]<<3);
1189
be972727 1190}
1191
1192//___________________________________________
e44235aa 1193Int_t AliMUONLocalTriggerBoard::GetI() const
be972727 1194{
71a2d3aa 1195/// old numbering
1196///
f4678817 1197 const Int_t kMaxfields = 2; char **fields = new char*[kMaxfields];
be972727 1198
1199 char s[100]; strcpy(s, GetName());
1200
1201 Int_t numlines = 0;
1202
1203 for (char *token = strtok(s, "B");
1204 token != NULL;
1205 token = strtok(NULL, " "))
1206 {
1207 fields[numlines] = new char[strlen(token)+1];
1208 strcpy(fields[numlines++], token);
1209 }
1210
1211 TString l(fields[0]);
1212
1213 char copy = l[0];
1214
f4678817 1215 Int_t lL = atoi(&l[4]), cC = atoi(&l[2]), sS = (copy=='R') ? +1 : -1;
be972727 1216
a6e0ebfe 1217 const char *b[4] = {"12", "34", "56", "78"};
be972727 1218
1219 Int_t ib = 0;
1220
1221 for (Int_t i=0; i<4; i++) if (!strcmp(fields[1],b[i])) {ib = i; break;} ib++;
1222
f4678817 1223// lL=1 ON TOP
1224 lL -= 9; lL = abs(lL); lL++;
be972727 1225
f4678817 1226 Int_t code = 100 * lL + 10 * cC + ib;
be972727 1227
f4678817 1228 code *= sS;
be972727 1229
1230 Int_t ic = 0;
1231
e44235aa 1232 for (Int_t i=0; i<234; i++) if (fgkCircuitId[i] == code) {ic = i; break;}
be972727 1233
1234 return ic;
1235}
1236
1237//___________________________________________
41922e78 1238void AliMUONLocalTriggerBoard::Mask(Int_t index, UShort_t mask)
be972727 1239{
71a2d3aa 1240/// set mask
1241///
41922e78 1242 if ( index >= 0 && index < 2*4 )
1243 {
1244 Int_t i = index/4;
1245 Int_t j = index - i*4;
1246 fMask[i][j]=mask;
1247 }
1248 else
1249 {
1250 AliError(Form("Index %d out of bounds (max %d)",index,8));
1251 }
be972727 1252}
1253
1254//___________________________________________
e44235aa 1255void AliMUONLocalTriggerBoard::Scan(Option_t *option) const
be972727 1256{
71a2d3aa 1257/// full dump
1258///
ac809573 1259 TString op = option;
be972727 1260
1261 if (op.Contains("CONF")) Conf();
1262
8af755f4 1263 if (op.Contains("BITP")) Pattern();
be972727 1264
1265 if (op.Contains("RESPI")) Resp("I");
1266
1267 if (op.Contains("RESPF")) Resp("F");
1268
00d46f24 1269 if (op.Contains("RESPO")) Resp("O");
1270
be972727 1271 if (op.Contains("ALL"))
1272 {
1273 Conf();
1274 Pattern();
1275 Resp("I");
1276 Resp("F");
1277 }
1278}
1279
1280//___________________________________________
ac809573 1281void AliMUONLocalTriggerBoard::Conf() const
1282{
1283/// board switches
1284///
1285 cout << "Switch(" << GetName() << ")"
1286 << " x2d = " << GetSwitch(0)
1287 << " x2m = " << GetSwitch(1)
1288 << " x2u = " << GetSwitch(2)
1289 << " OR[0] = " << GetSwitch(3)
1290 << " OR[1] = " << GetSwitch(4)
1291 << " EN-Y = " << GetSwitch(5)
1292 << " ZERO-ALLY-LSB = " << GetSwitch(6)
1293 << " ZERO-down = " << GetSwitch(7)
1294 << " ZERO-middle = " << GetSwitch(8)
1295 << " ZERO-up = " << GetSwitch(9)
1296 << " trans. conn. " << GetTC()
1297 << " Slot = " << fSlot
1298 << endl;
1299}
1300
1301//___________________________________________
e44235aa 1302void AliMUONLocalTriggerBoard::Resp(Option_t *option) const
be972727 1303{
71a2d3aa 1304/// board I/O
1305///
be972727 1306 TString op = option;
1307
1308 if (op.Contains("I"))
1309 {
1310// print Local trigger output before the LuT step
1311 printf("===============================================================\n");
1312 printf("-------- TRIGGER OUTPUT --------\n");
1313 printf("minDevStrip = ");
1314 for (Int_t i=4; i>=0; i--) printf("%i",fMinDevStrip[i]);
1315 printf(" minDev = ");
1316 for (Int_t i=4; i>=0; i--) printf("%i",fMinDev[i]);
1317 printf(" coordY = ");
1318 for (Int_t i=4; i>=0; i--) printf("%i",fCoordY[i]);
1319 printf(" \n");
1320 }
00d46f24 1321 if (op.Contains("O")){
1322 printf("-- Output --\n");
1323 printf("Coinc44 %i\n", fCoinc44);
1324 printf("StripX11 %i StripY11 %i Dev %i TrigY %i\n", fStripX11, fStripY11, fDev, fTrigY);
1325 printf("LutLpt %i %i LutHpt %i %i\n", fLutLpt[0], fLutLpt[1], fLutHpt[0], fLutHpt[1]);
1326 printf("Response %i\n", fOutput);
1327 }
be972727 1328
be972727 1329}
1330
1331//___________________________________________
be972727 1332void AliMUONLocalTriggerBoard::Response()
1333{
71a2d3aa 1334/// algo
1335///
00d46f24 1336 if ( IsNull() ) {
1337 fMinDev[4] = 1;
1338 for (Int_t i=4; i>=0; i--) fCoordY[i] = 1;
1339 return; // Do nothing if strip pattern is null
1340 }
1341
f4678817 1342 Int_t xX1[16], xX2[16], xXX3[32], xXX4[32];
be972727 1343
1344 TBits x1(16), x2(16), x3(16), x4(16);
1345
1346 UShort_t xyv = 0;
1347
1348 xyv = fXY[0][0]; x1.Set(16,&xyv);
1349 xyv = fXY[0][1]; x2.Set(16,&xyv);
1350 xyv = fXY[0][2]; x3.Set(16,&xyv);
1351 xyv = fXY[0][3]; x4.Set(16,&xyv);
1352
1353 TBits x3u(16), x4u(16), x3d(16), x4d(16);
1354
1355 xyv = fXYU[0][2]; x3u.Set(16,&xyv);
1356 xyv = fXYU[0][3]; x4u.Set(16,&xyv);
1357
1358 xyv = fXYD[0][2]; x3d.Set(16,&xyv);
1359 xyv = fXYD[0][3]; x4d.Set(16,&xyv);
1360
1361 for (Int_t i=0;i<16;i++)
1362 {
f4678817 1363 xX1[i] = x1[i];
1364 xX2[i] = x2[i];
be972727 1365
f4678817 1366 xXX3[i+8] = x3[i];
1367 xXX4[i+8] = x4[i];
be972727 1368 }
1369
1370 for (Int_t i=0;i<8;i++)
1371 {
f4678817 1372 xXX3[i] = x3d[i+8];
1373 xXX4[i] = x4d[i+8];
be972727 1374
f4678817 1375 xXX3[i+24] = x3u[i];
1376 xXX4[i+24] = x4u[i];
be972727 1377 }
1378
e93b11d0 1379// Int_t coinc44 = 0;
be972727 1380
e93b11d0 1381 TrigX(xX1, xX2, xXX3, xXX4);
be972727 1382
f4678817 1383 Int_t yY1[16], yY2[16], yY3[16], yY4[16];
be972727 1384
f4678817 1385 Int_t yY3U[16], yY3D[16], yY4U[16], yY4D[16];
be972727 1386
1387 TBits y1(16), y2(16), y3(16), y4(16);
1388
1389 TBits y3u(16), y3d(16), y4u(16), y4d(16);
1390
1391 xyv = fXY[1][0]; y1.Set(16,&xyv);
1392 xyv = fXY[1][1]; y2.Set(16,&xyv);
1393 xyv = fXY[1][2]; y3.Set(16,&xyv);
1394 xyv = fXY[1][3]; y4.Set(16,&xyv);
1395
1396 xyv = fXYU[1][2]; y3u.Set(16,&xyv);
1397 xyv = fXYD[1][2]; y3d.Set(16,&xyv);
1398 xyv = fXYU[1][3]; y4u.Set(16,&xyv);
1399 xyv = fXYD[1][3]; y4d.Set(16,&xyv);
1400
1401 for (Int_t i=0;i<16;i++)
1402 {
f4678817 1403 yY1[i] = y1[i];
1404 yY2[i] = y2[i];
1405 yY3[i] = y3[i];
1406 yY4[i] = y4[i];
be972727 1407
f4678817 1408 yY3U[i] = y3u[i];
1409 yY3D[i] = y3d[i];
be972727 1410
f4678817 1411 yY4U[i] = y4u[i];
1412 yY4D[i] = y4d[i];
be972727 1413 }
1414
e93b11d0 1415 TrigY(yY1, yY2, yY3, yY4, yY3U, yY3D, yY4U, yY4D);
be972727 1416
ac809573 1417// ASIGN fLutLpt, fLutHpt
abe5ad6f 1418 LocalTrigger();
be972727 1419}
ac809573 1420
5a8d50e9 1421//___________________________________________
1422Bool_t AliMUONLocalTriggerBoard::IsTrigY() const
1423{
1424 /// Return the response of non-bending plane
1425 Int_t iStripY = 0;
1426 Bool_t output = kFALSE;
1427
1428 for (Int_t i=0; i<4; i++) iStripY += static_cast<int>( fCoordY[i] << i );
1429
1430 if (fCoordY[4]==1 && iStripY==15) output=kFALSE; // no trigger in Y
1431 else output=kTRUE; // trigger in Y
1432
1433 return output;
1434}
1435
1436//___________________________________________
1437Bool_t AliMUONLocalTriggerBoard::IsTrigX() const
1438{
1439 /// Return the response of bending plane
1440
1441 Int_t deviation = 0;
1442 Bool_t output = kFALSE;
1443
1444 for (Int_t i=0; i<4; i++) deviation += static_cast<int>( fMinDev[i] << i );
1445
1446 if (fMinDev[4]==1 && !deviation) output=kFALSE; // no trigger in X
1447 else output=kTRUE; // trigger in X
1448
1449 return output;
1450}
00d46f24 1451
1452//___________________________________________
1453Bool_t AliMUONLocalTriggerBoard::IsNull() const
1454{
1455 /// Check if board has fired strips in the first station
1456 for (Int_t icath=0; icath<2; icath++){
1457 for(Int_t ich=0; ich<2; ich++){
1458 if ( fXY[icath][ich] ) return kFALSE;
1459 }
1460 }
1461 return kTRUE;
1462}