Fix coverity defect
[u/mrichter/AliRoot.git] / MUON / AliMUONGlobalTriggerBoard.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$ */
17
18 //-----------------------------------------------------------------------------
19 /// \class AliMUONGlobalTriggerBoard
20 /// Global trigger implementation:
21 /// - inputs are regional responses
22 /// - output is a 12-bit word
23 /// - 4 bits per trigger level
24 ///
25 /// \author Rachid Guernane (LPCCFd), 
26 /// Corrected by Christian Finck (Subatech)
27 //-----------------------------------------------------------------------------
28
29 #include "AliMUONGlobalTriggerBoard.h"
30 #include "AliLog.h"
31 #include "TBits.h"
32
33 #include <Riostream.h>
34
35 /// \cond CLASSIMP
36 ClassImp(AliMUONGlobalTriggerBoard)
37 /// \endcond
38
39 //___________________________________________
40 AliMUONGlobalTriggerBoard::AliMUONGlobalTriggerBoard(): AliMUONTriggerBoard()
41 {
42 /// Default constructor
43
44    for (Int_t i=0;i<16;i++) fRegionalResponse[i] = 0;
45    for (Int_t i=0;i< 4;i++) fGlobalInput[i] = 0;
46    for (Int_t i=0;i< 4;i++) fMask[i] = 0xffffffff;
47 }
48
49 //___________________________________________
50 AliMUONGlobalTriggerBoard::AliMUONGlobalTriggerBoard(const char *name, Int_t a) : AliMUONTriggerBoard(name, a)
51 {
52 /// Standard constructor
53
54    for (Int_t i=0;i<16;i++) fRegionalResponse[i] = 0;
55    for (Int_t i=0;i< 4;i++) fGlobalInput[i] = 0;
56    for (Int_t i=0;i< 4;i++) fMask[i] = 0xffffffff;
57 }
58
59 //___________________________________________
60 AliMUONGlobalTriggerBoard::~AliMUONGlobalTriggerBoard()
61 {
62 /// Destructor
63 }
64
65 //___________________________________________
66 void AliMUONGlobalTriggerBoard::Mask(Int_t index, UInt_t mask)
67 {
68   /// mask global trigger board input index with value mask
69   if ( index >= 0 && index < 4 ) 
70   {
71     fMask[index]=mask;
72   }
73   else
74   {
75     AliError(Form("Index %d out of bounds (max %d)",index,3));
76   }  
77 }
78
79 //___________________________________________
80 void AliMUONGlobalTriggerBoard::Response()
81 {
82    /// compute the global trigger board
83    /// response according to the algo() method
84 // output from global trigger algorithm
85 // [+, -, US, LS] * [Hpt, Lpt]
86 // transformed to [usHpt, usLpt, lsHpt, lsLpt, sHpt, sLpt] according
87 // to Global Trigger Unit user manual
88
89    Int_t t[16];
90
91    BuildGlobalInput();
92    MaskGlobalInput();
93
94    for (Int_t i = 0; i < 16; ++i) 
95    {
96      t[i] = fRegionalResponse[i];
97    }
98    
99    
100    Int_t rank = 8;
101
102    for (Int_t i=0;i<4;i++)
103    {
104       Int_t ip = 0;
105       
106       for (Int_t j=0;j<rank;j++)
107       {
108          UShort_t lthres = Algo(t[2*j],t[2*j+1],"LPT");
109
110          UShort_t hthres = Algo(t[2*j],t[2*j+1],"HPT"); hthres <<= 4;
111
112          t[ip] = lthres | hthres;
113
114          ip++;
115       }
116       
117       rank /= 2; 
118    }
119    UChar_t sLpt, sHpt, lsLpt, lsHpt, usLpt, usHpt;
120    sLpt  = ((t[0] & 0xC)  != 0);
121    sHpt  = ((t[0] & 0xC0) != 0);
122    lsLpt = ((t[0] & 0x1)  != 0);
123    lsHpt = ((t[0] & 0x10) != 0);
124    usLpt = ((t[0] & 0x2 ) != 0);
125    usHpt = ((t[0] & 0x20) != 0);
126
127    // LSB is zero (trigger choice to send to CTP: sLpt or sHpt)
128
129    sLpt  <<= 1;
130    sHpt  <<= 2;
131    lsLpt <<= 3;
132    lsHpt <<= 4;
133    usLpt <<= 5;
134    usHpt <<= 6;
135
136    fResponse = sLpt | sHpt | lsLpt | lsHpt | usLpt |usHpt;
137    
138
139 }
140
141 //___________________________________________
142 UShort_t AliMUONGlobalTriggerBoard::Algo(UShort_t i, UShort_t j, const char *thres)
143 {
144 /// global trigger algorithm
145 ///   a ,b = reg  response  =  Hpt (+|-|us|ls) |  Lpt (+|-|us|ls)  
146                            
147    TBits a(8), b(8); a.Set(8,&i); b.Set(8,&j);
148
149    TBits trg1(2), trg2(2), trg(2);
150
151    if (!strcmp(thres,"LPT"))
152    {
153       trg1[0] = a[2]; trg1[1] = a[3]; 
154       trg2[0] = b[2]; trg2[1] = b[3];
155    }
156    else
157    {
158       trg1[0] = a[6]; trg1[1] = a[7]; 
159       trg2[0] = b[6]; trg2[1] = b[7];         
160    }
161        
162    TBits trgLS1(1), trgUS1(1), trgLS2(1), trgUS2(1), trgLS(1), trgUS(1);
163
164    if (!strcmp(thres,"LPT"))
165    {
166       trgLS1[0] = a[0]; trgUS1[0] = a[1]; 
167       trgLS2[0] = b[0]; trgUS2[0] = b[1];
168    }
169    else
170    {
171       trgLS1[0] = a[4]; trgUS1[0] = a[5]; 
172       trgLS2[0] = b[4]; trgUS2[0] = b[5];         
173    }
174
175    trgLS[0] = ( trg1[0] & trg2[0] ) | ( trg1[1] & trg2[1] ) | trgLS1[0] | trgLS2[0];
176    trgUS[0] = ( trg1[0] & trg2[1] ) | ( trg1[1] & trg2[0] ) | trgUS1[0] | trgUS2[0];
177    
178    trg[0] = trg1[0] | trg2[0];
179    trg[1] = trg1[1] | trg2[1];
180    
181    TBits v(4);
182    
183    v[0] = trgLS[0];
184    v[1] = trgUS[0];
185    v[2] = trg[0];
186    v[3] = trg[1];
187
188    UShort_t rv = 0;
189    v.Get(&rv);
190    
191    return rv;
192 }
193
194 //___________________________________________
195 void AliMUONGlobalTriggerBoard::BuildGlobalInput()
196 {
197   /// build the 4 words (32bits) global input from the regional responses
198   /// the order of regional responses is:
199   /// 1R, 2R, 2-3R, 3R, 4R, 5R, 6R, 7R, 1L, 2L, 2-3L, 3L, 4L, 5L, 6L, 7L
200
201   for (Int_t i=0;i< 4;i++) fGlobalInput[i] = 0;
202
203   UShort_t regRespInv;
204   for (Int_t iReg = 0; iReg < 16; iReg++) {
205
206     regRespInv = InvertPairBits(iReg);
207
208     if (iReg < 8) {    // right
209       // Lpt word
210       fGlobalInput[0] |=  (regRespInv & 0x0F)       << (4*iReg);
211       // Hpt word
212       fGlobalInput[2] |= ((regRespInv & 0xF0) >> 4) << (4*iReg);
213     } else {           // left
214       // Lpt word
215       fGlobalInput[1] |=  (regRespInv & 0x0F)       << (4*(iReg-8));
216       // Hpt word
217       fGlobalInput[3] |= ((regRespInv & 0xF0) >> 4) << (4*(iReg-8));
218     }
219
220   }
221
222 }
223
224 //___________________________________________
225 void AliMUONGlobalTriggerBoard::MaskGlobalInput()
226 {
227   /// Apply masks to global input and recalculate regional inputs before
228   /// applying the global response
229
230   UInt_t gitmp[4];
231   for (Int_t i = 0; i < 4; i++) {
232     fGlobalInput[i] &= fMask[i];
233     gitmp[i] = fGlobalInput[i];
234   }
235
236   RecomputeRegional(gitmp);
237 }
238
239
240 //___________________________________________
241 void AliMUONGlobalTriggerBoard::RecomputeRegional(UInt_t gitmp[4])
242 {
243   //
244   /// Recomput regional response from global input
245   //
246   for (Int_t iReg = 0; iReg < 16; iReg++) {
247     fRegionalResponse[iReg] = 0;
248     if (iReg < 8) {    // right
249       // Lpt
250       fRegionalResponse[iReg] |=  (gitmp[0] >> (4*iReg))     & 0xF;
251       // Hpt
252       fRegionalResponse[iReg] |= ((gitmp[2] >> (4*iReg))     & 0xF) << 4;
253     } else {           // left
254       // Lpt
255       fRegionalResponse[iReg] |=  (gitmp[1] >> (4*(iReg-8))) & 0xF;
256       // Hpt
257       fRegionalResponse[iReg] |= ((gitmp[3] >> (4*(iReg-8))) & 0xF) << 4;
258     }
259
260     fRegionalResponse[iReg] = InvertPairBits(iReg);
261   }
262 }
263
264 //___________________________________________
265 void AliMUONGlobalTriggerBoard::Scan(Option_t*) const
266 {
267   /// print global trigger output 
268   TBits w(7); w.Set(7,&fResponse);
269
270 // TRG[1:0]
271 // 00 noth
272 // 01 negative track
273 // 10 positive track
274 // 11 undef
275
276    Int_t iS[2] = {0,0};
277
278    iS[0] = (Int_t)w.TestBitNumber(1);
279    iS[1] = (Int_t)w.TestBitNumber(2);
280
281    Int_t iPU[2] = {w[5],w[6]};
282    Int_t iPL[2] = {w[3],w[4]};
283
284    printf("============================================\n");
285    printf(" Global Trigger output       Low pt  High pt\n");
286    printf(" number of Single           :\t");
287    for (Int_t i=0; i<2; i++) printf("%i\t",iS[i]);
288    printf("\n");
289    printf(" number of UnlikeSign pair  :\t"); 
290    for (Int_t i=0; i<2; i++) printf("%i\t",iPU[i]);
291    printf("\n");
292    printf(" number of LikeSign pair    :\t");  
293    for (Int_t i=0; i<2; i++) printf("%i\t",iPL[i]);
294    printf("\n");
295    printf("===================================================\n");
296    printf("\n");
297 }
298
299
300 //___________________________________________
301 UShort_t AliMUONGlobalTriggerBoard::InvertPairBits(Int_t iReg)
302 {
303   //
304   /// invert "pair" bits in regional response
305   /// [+, -, US, LS] becomes [+, -, LS, US]
306   //
307   TBits rs(8), rsi(8);
308   rs.Set(8,&fRegionalResponse[iReg]);
309     for (Int_t i = 0; i < 4; i++) {
310       if (i%2 == 0) {
311         rsi[2*i]   = rs[2*i+1];
312         rsi[2*i+1] = rs[2*i];
313       } else {
314         rsi[2*i]   = rs[2*i];
315         rsi[2*i+1] = rs[2*i+1];
316       }
317     }
318     UShort_t regRespInv = 0;
319     rsi.Get(&regRespInv);
320     return regRespInv;
321 }