]>
Commit | Line | Data |
---|---|---|
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 | } | |
46 | ||
47 | //___________________________________________ | |
48 | AliMUONGlobalTriggerBoard::AliMUONGlobalTriggerBoard(const char *name, Int_t a) : AliMUONTriggerBoard(name, a) | |
49 | { | |
50 | /// Standard constructor | |
51 | ||
52 | for (Int_t i=0;i<16;i++) fRegionalResponse[i] = 0; | |
53 | } | |
54 | ||
55 | //___________________________________________ | |
56 | AliMUONGlobalTriggerBoard::~AliMUONGlobalTriggerBoard() | |
57 | { | |
58 | /// Destructor | |
59 | } | |
60 | ||
61 | //___________________________________________ | |
62 | void AliMUONGlobalTriggerBoard::Mask(Int_t index, UShort_t mask) | |
63 | { | |
64 | /// mask global trigger board input index with value mask | |
65 | if ( index>=0 && index < 2 ) | |
66 | { | |
67 | fMask[index]=mask; | |
68 | } | |
69 | else | |
70 | { | |
71 | AliError(Form("Index %d out of bounds (max %d)",index,2)); | |
72 | } | |
73 | } | |
74 | ||
75 | //___________________________________________ | |
76 | void AliMUONGlobalTriggerBoard::Response() | |
77 | { | |
78 | /// compute the global trigger board | |
79 | /// response according to the algo() method | |
80 | // output from global trigger algorithm | |
81 | // [+, -, US, LS] * [Hpt, Lpt] | |
82 | // transformed to [usHpt, usLpt, lsHpt, lsLpt, sHpt, sLpt] according | |
83 | // to Global Trigger Unit user manual | |
84 | ||
85 | Int_t t[16]; | |
86 | ||
87 | for (Int_t i = 0; i < 16; ++i) | |
88 | { | |
89 | Int_t index = i/8; | |
90 | Int_t shift = i % 8; | |
91 | UShort_t enable = !((fMask[index] >> shift) & 0x1); | |
92 | if (enable) | |
93 | t[i] = fRegionalResponse[i]; | |
94 | else | |
95 | t[i] = 0; | |
96 | } | |
97 | ||
98 | ||
99 | Int_t rank = 8; | |
100 | ||
101 | for (Int_t i=0;i<4;i++) | |
102 | { | |
103 | Int_t ip = 0; | |
104 | ||
105 | for (Int_t j=0;j<rank;j++) | |
106 | { | |
107 | UShort_t lthres = Algo(t[2*j],t[2*j+1],"LPT"); | |
108 | ||
109 | UShort_t hthres = Algo(t[2*j],t[2*j+1],"HPT"); hthres <<= 4; | |
110 | ||
111 | t[ip] = lthres | hthres; | |
112 | ||
113 | ip++; | |
114 | } | |
115 | ||
116 | rank /= 2; | |
117 | } | |
118 | UChar_t sLpt, sHpt, lsLpt, lsHpt, usLpt, usHpt; | |
119 | sLpt = ((t[0] & 0xC) != 0); | |
120 | sHpt = ((t[0] & 0xC0) != 0); | |
121 | //lsLpt = ((t[0] & 0x2) != 0); | |
122 | lsLpt = ((t[0] & 0x1) != 0); | |
123 | //lsHpt = ((t[0] & 0x20) != 0); | |
124 | lsHpt = ((t[0] & 0x10) != 0); | |
125 | //usLpt = ((t[0] & 0x1 ) != 0); | |
126 | usLpt = ((t[0] & 0x2 ) != 0); | |
127 | //usHpt = ((t[0] & 0x10) != 0); | |
128 | usHpt = ((t[0] & 0x20) != 0); | |
129 | ||
130 | sHpt <<= 1; | |
131 | lsLpt <<= 2; | |
132 | lsHpt <<= 3; | |
133 | usLpt <<= 4; | |
134 | usHpt <<= 5; | |
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, 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[1]; trgUS1[0] = a[0]; | |
167 | //trgLS2[0] = b[1]; trgUS2[0] = b[0]; | |
168 | trgLS1[0] = a[0]; trgUS1[0] = a[1]; | |
169 | trgLS2[0] = b[0]; trgUS2[0] = b[1]; | |
170 | } | |
171 | else | |
172 | { | |
173 | //trgLS1[0] = a[5]; trgUS1[0] = a[4]; | |
174 | //trgLS2[0] = b[5]; trgUS2[0] = b[4]; | |
175 | trgLS1[0] = a[4]; trgUS1[0] = a[5]; | |
176 | trgLS2[0] = b[4]; trgUS2[0] = b[5]; | |
177 | } | |
178 | ||
179 | trgLS[0] = ( trg1[0] & trg2[0] ) | ( trg1[1] & trg2[1] ) | trgLS1[0] | trgLS2[0]; | |
180 | trgUS[0] = ( trg1[0] & trg2[1] ) | ( trg1[1] & trg2[0] ) | trgUS1[0] | trgUS2[0]; | |
181 | ||
182 | trg[0] = trg1[0] | trg2[0]; | |
183 | trg[1] = trg1[1] | trg2[1]; | |
184 | ||
185 | TBits v(4); | |
186 | ||
187 | //v[0] = trgUS[0]; | |
188 | //v[1] = trgLS[0]; | |
189 | v[0] = trgLS[0]; | |
190 | v[1] = trgUS[0]; | |
191 | v[2] = trg[0]; | |
192 | v[3] = trg[1]; | |
193 | ||
194 | UShort_t rv = 0; | |
195 | v.Get(&rv); | |
196 | ||
197 | return rv; | |
198 | } | |
199 | ||
200 | //___________________________________________ | |
201 | void AliMUONGlobalTriggerBoard::Scan(Option_t*) const | |
202 | { | |
203 | /// print global trigger output | |
204 | TBits w(6); w.Set(6,&fResponse); | |
205 | ||
206 | // TRG[1:0] | |
207 | // 00 noth | |
208 | // 01 negative track | |
209 | // 10 positive track | |
210 | // 11 undef | |
211 | ||
212 | Int_t iS[2] = {0,0}; | |
213 | ||
214 | iS[0] = (Int_t)w.TestBitNumber(0); | |
215 | iS[1] = (Int_t)w.TestBitNumber(1); | |
216 | ||
217 | Int_t iPU[2] = {w[4],w[5]}; | |
218 | Int_t iPL[2] = {w[2],w[3]}; | |
219 | ||
220 | printf("============================================\n"); | |
221 | printf(" Global Trigger output Low pt High pt\n"); | |
222 | printf(" number of Single :\t"); | |
223 | for (Int_t i=0; i<2; i++) printf("%i\t",iS[i]); | |
224 | printf("\n"); | |
225 | printf(" number of UnlikeSign pair :\t"); | |
226 | for (Int_t i=0; i<2; i++) printf("%i\t",iPU[i]); | |
227 | printf("\n"); | |
228 | printf(" number of LikeSign pair :\t"); | |
229 | for (Int_t i=0; i<2; i++) printf("%i\t",iPL[i]); | |
230 | printf("\n"); | |
231 | printf("===================================================\n"); | |
232 | printf("\n"); | |
233 | } | |
234 |