remove sign info on trigger single (Christian, Philippe)
[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 //*-- Author: Rachid Guernane (LPCCFd)
19 //    GLOBAL TRIGGER IMPLEMENTATION
20 //    INPUTS ARE REGIONAL RESPONSES
21 //    OUTPUT IS A 12-BIT WORD
22 //    4 BITS PER TRIGGER LEVEL
23
24 #include "AliMUONGlobalTriggerBoard.h"
25 #include "AliLog.h"
26 #include "TBits.h"
27
28 #include <Riostream.h>
29
30 ClassImp(AliMUONGlobalTriggerBoard)
31
32 //___________________________________________
33 AliMUONGlobalTriggerBoard::AliMUONGlobalTriggerBoard()
34 {
35    for (Int_t i=0;i<16;i++) fRegionalResponse[i] = 0;
36 }
37
38 //___________________________________________
39 AliMUONGlobalTriggerBoard::AliMUONGlobalTriggerBoard(const char *name, Int_t a) : AliMUONTriggerBoard(name, a)
40 {
41    for (Int_t i=0;i<16;i++) fRegionalResponse[i] = 0;
42 }
43
44 //___________________________________________
45 void AliMUONGlobalTriggerBoard::Mask(Int_t index, UShort_t mask)
46 {
47   // MASK GLOBAL TRIGGER BOARD INPUT index WITH VALUE mask
48   if ( index>=0 && index < 16 ) 
49   {
50     fMask[index]=mask;
51   }
52   else
53   {
54     AliError(Form("Index %d out of bounds (max %d)",index,16));
55   }  
56 }
57
58 //___________________________________________
59 void AliMUONGlobalTriggerBoard::Response()
60 {
61    // COMPUTE THE GLOBAL TRIGGER BOARD
62    // RESPONSE ACCORDING TO THE Algo() METHOD
63 /// output from global trigger algorithm
64 /// [+, -, LS, US] * [Hpt, Lpt]
65 /// transformed to [usHpt, usLpt, lsHpt, lsLpt, sHpt, sLpt] according
66 /// to Global Trigger Unit user manual
67
68    Int_t t[16];
69    for (Int_t i=0;i<16;i++) t[i] = fRegionalResponse[i] & fMask[i];
70
71    Int_t rank = 8;
72
73    for (Int_t i=0;i<4;i++)
74    {
75       Int_t ip = 0;
76       
77       for (Int_t j=0;j<rank;j++)
78       {
79          UShort_t lthres = Algo(t[2*j],t[2*j+1],"LPT");
80
81          UShort_t hthres = Algo(t[2*j],t[2*j+1],"HPT"); hthres <<= 4;
82
83          t[ip] = lthres | hthres;
84
85          ip++;
86       }
87       
88       rank /= 2; 
89    }
90    cout << " ===============================================" << "\n";
91   cout << " here ----------------- " << t[0] << "\n";
92   cout << " ===============================================" << "\n";
93    UChar_t sLpt, sHpt, lsLpt, lsHpt, usLpt, usHpt;
94    sLpt  = ((t[0] & 0xC)  != 0);
95    sHpt  = ((t[0] & 0xC0) != 0);
96    lsLpt = ((t[0] & 0x2)  != 0);
97    lsHpt = ((t[0] & 0x20) != 0);
98    usLpt = ((t[0] & 0x1 ) != 0);
99    usHpt = ((t[0] & 0x10) != 0);
100
101    sHpt  <<= 1;
102    lsLpt <<= 2;
103    lsHpt <<= 3;
104    usLpt <<= 4;
105    usHpt <<= 5;
106
107    fResponse = sLpt | sHpt | lsLpt | lsHpt | usLpt |usHpt;
108
109    cout << " ===============================================" << "\n";
110    cout << sLpt << " " << sHpt << " " << lsLpt << " "<<  lsHpt << " " 
111         << " " << usLpt  << " " << usHpt <<"\n";
112   cout << " here +++++++++++++++++++ " << fResponse << "\n";
113   cout << " ===============================================" << "\n";
114 }
115
116 //___________________________________________
117 UShort_t AliMUONGlobalTriggerBoard::Algo(UShort_t i, UShort_t j, char *thres)
118 {
119    // GLOBAL TRIGGER ALGORITHM
120    TBits a(8), b(8); a.Set(8,&i); b.Set(8,&j);
121
122    TBits trg1(2), trg2(2), trg(2);
123
124    if (!strcmp(thres,"LPT"))
125    {
126       trg1[0] = a[2]; trg1[1] = a[3]; 
127       trg2[0] = b[2]; trg2[1] = b[3];
128    }
129    else
130    {
131       trg1[0] = a[6]; trg1[1] = a[7]; 
132       trg2[0] = b[6]; trg2[1] = b[7];         
133    }
134        
135    TBits trgLS1(1), trgUS1(1), trgLS2(1), trgUS2(1), trgLS(1), trgUS(1);
136
137    if (!strcmp(thres,"LPT"))
138    {
139       trgLS1[0] = a[1]; trgUS1[0] = a[0]; 
140       trgLS2[0] = b[1]; trgUS2[0] = b[0];
141    }
142    else
143    {
144       trgLS1[0] = a[5]; trgUS1[0] = a[4]; 
145       trgLS2[0] = b[5]; trgUS2[0] = b[4];         
146    }
147
148    trgLS[0] = ( trg1[0] & trg2[0] ) | ( trg1[1] & trg2[1] ) | trgLS1[0] | trgLS2[0];
149    trgUS[0] = ( trg1[0] & trg2[1] ) | ( trg1[1] & trg2[0] ) | trgUS1[0] | trgUS2[0];
150    
151    trg[0] = trg1[0] | trg2[0];
152    trg[1] = trg1[1] | trg2[1];
153    
154    TBits v(4);
155    
156    v[0] = trgUS[0];
157    v[1] = trgLS[0];
158    v[2] = trg[0];
159    v[3] = trg[1];
160
161    UShort_t rv = 0;
162    v.Get(&rv);
163    
164    return rv;
165 }
166
167 //___________________________________________
168 void AliMUONGlobalTriggerBoard::Scan(Option_t*) const
169 {
170   // PRINT GLOBAL TRIGGER OUTPUT 
171   TBits w(8); w.Set(8,&fResponse);
172
173 // TRG[1:0]
174 // 00 noth
175 // 01 negative track
176 // 10 positive track
177 // 11 undef
178
179    Int_t iSP[2] = {0,0}, iSM[2] = {0,0}, iSU[2] = {0,0};
180
181    TBits a(2), n(2), p(2), u(2);
182    
183    UShort_t val;
184
185    val = 1; n.Set(2,&val);
186    val = 2; p.Set(2,&val);
187    val = 3; u.Set(2,&val);
188    
189    a[0] = w[2];
190    a[1] = w[3];
191    
192    if      (a==p) iSP[0] = 1;
193    else if (a==n) iSM[0] = 1;
194    else if (a==u) iSU[0] = 1;   
195
196    a[0] = w[6];
197    a[1] = w[7];
198
199    if      (a==p) iSP[1] = 1;
200    else if (a==n) iSM[1] = 1;
201    else if (a==u) iSU[1] = 1;
202    
203    Int_t iPU[2] = {w[0],w[4]};
204    Int_t iPL[2] = {w[1],w[5]};
205
206    printf("============================================\n");
207    printf(" Global Trigger output       Low pt  High pt\n");
208    printf(" number of Single Plus      :\t");
209    for (Int_t i=0; i<2; i++) printf("%i\t",iSP[i]);
210    printf("\n");
211    printf(" number of Single Minus     :\t");
212    for (Int_t i=0; i<2; i++) printf("%i\t",iSM[i]);
213    printf("\n");
214    printf(" number of Single Undefined :\t"); 
215    for (Int_t i=0; i<2; i++) printf("%i\t",iSU[i]);
216    printf("\n");
217    printf(" number of UnlikeSign pair  :\t"); 
218    for (Int_t i=0; i<2; i++) printf("%i\t",iPU[i]);
219    printf("\n");
220    printf(" number of LikeSign pair    :\t");  
221    for (Int_t i=0; i<2; i++) printf("%i\t",iPL[i]);
222    printf("\n");
223    printf("===================================================\n");
224    printf("\n");
225 }
226
227 ClassImp(AliMUONGlobalTriggerBoard)