]>
Commit | Line | Data |
---|---|---|
d810d0de | 1 | // $Id$ |
2 | // Main authors: Matevz Tadel & Alja Mrak-Tadel: 2006, 2007 | |
3 | ||
4 | /************************************************************************** | |
5 | * Copyright(c) 1998-2008, ALICE Experiment at CERN, all rights reserved. * | |
6 | * See http://aliceinfo.cern.ch/Offline/AliRoot/License.html for * | |
51346b82 | 7 | * full copyright notice. * |
d810d0de | 8 | **************************************************************************/ |
9 | #include "AliEveTOFSector.h" | |
8616d353 | 10 | |
8616d353 | 11 | #include <AliTOFdigit.h> |
28a21457 | 12 | #include <AliTOFGeometry.h> |
8616d353 | 13 | |
a15e6d7d | 14 | #include <TEveTrans.h> |
15 | #include <TEveManager.h> | |
8616d353 | 16 | #include <TStyle.h> |
d810d0de | 17 | |
fd31e9de | 18 | Bool_t AliEveTOFSector::fgStaticInitDone = kFALSE; |
d810d0de | 19 | TEveFrameBox* AliEveTOFSector::fgTOFsectorFrameBox = 0; |
20 | TEveRGBAPalette* AliEveTOFSector::fgTOFsectorPalette = 0; | |
8616d353 | 21 | |
22 | //_______________________________________________________ | |
d810d0de | 23 | ClassImp(AliEveTOFSector) |
8616d353 | 24 | |
25 | /* ************************************************************************ */ | |
26 | ||
d810d0de | 27 | AliEveTOFSector::AliEveTOFSector(const Text_t* n, const Text_t* t) : |
84aff7a4 | 28 | TEveQuadSet(n, t), |
28a21457 | 29 | fTOFgeometry(new AliTOFGeometry()), |
8616d353 | 30 | fTOFarray(0x0), |
31 | fTOFtree(0x0), | |
32 | fSector(-1), | |
33 | fDx(0), fDy(0), fDz(0), | |
34 | fAutoTrans (kTRUE), | |
ee9d8eca | 35 | //fMinTime (0), fMaxTime (0), |
36 | fThreshold (5), fMaxVal (80), | |
8616d353 | 37 | fSectorID (0), |
ee9d8eca | 38 | //fPlateFlag0(kTRUE), fPlateFlag1(kTRUE), fPlateFlag2(kTRUE), fPlateFlag3(kTRUE), fPlateFlag4(kTRUE), |
39 | fPlateFlag (0x0), | |
40 | //fFrameColor(4), | |
41 | //fRnrFrame (kFALSE), | |
8616d353 | 42 | fGeoManager(0) |
43 | { | |
44 | ||
45 | fPlateFlag = new Bool_t[5]; | |
46 | for (Int_t ii=0; ii<5; ii++) fPlateFlag[ii]=kTRUE; | |
47 | ||
ee9d8eca | 48 | |
81515106 | 49 | //fGeoManager = AliEveEventManager::AssertGeometry(); |
8616d353 | 50 | if (!fGeoManager) { |
51 | printf("ERROR: no TGeo\n"); | |
52 | } | |
fd31e9de | 53 | |
ee9d8eca | 54 | } |
8616d353 | 55 | /* ************************************************************************ */ |
56 | ||
ee9d8eca | 57 | AliEveTOFSector::AliEveTOFSector(TGeoManager *localGeoManager, |
58 | Int_t nSector) | |
59 | : | |
60 | TEveQuadSet(Form("Sector%i",nSector)), | |
28a21457 | 61 | fTOFgeometry(new AliTOFGeometry()), |
8616d353 | 62 | fTOFarray(0x0), |
63 | fTOFtree(0x0), | |
64 | fSector(nSector), | |
65 | fDx(0), fDy(0), fDz(0), | |
66 | fAutoTrans (kTRUE), | |
ee9d8eca | 67 | //fMinTime (0), fMaxTime (0), |
fd31e9de | 68 | fThreshold (5), fMaxVal (80), |
8616d353 | 69 | fSectorID (nSector), |
ee9d8eca | 70 | //fPlateFlag0(kTRUE), fPlateFlag1(kTRUE), fPlateFlag2(kTRUE), fPlateFlag3(kTRUE), fPlateFlag4(kTRUE), |
71 | fPlateFlag (0x0), | |
72 | //fFrameColor(4), | |
73 | //fRnrFrame (kFALSE), | |
8616d353 | 74 | fGeoManager(localGeoManager) |
75 | { | |
76 | ||
77 | fPlateFlag = new Bool_t[5]; | |
78 | for (Int_t ii=0; ii<5; ii++) fPlateFlag[ii]=kTRUE; | |
79 | ||
80 | /* | |
81 | if (!fGeoManager) { | |
82 | printf("ERROR: no TGeo\n"); | |
83 | } | |
84 | */ | |
85 | ||
86 | InitModule(); | |
fd31e9de | 87 | |
ee9d8eca | 88 | } |
8616d353 | 89 | /* ************************************************************************ */ |
90 | ||
d810d0de | 91 | AliEveTOFSector::AliEveTOFSector(TGeoManager *localGeoManager, |
8616d353 | 92 | Int_t nSector, |
93 | TClonesArray *tofArray) | |
94 | : | |
84aff7a4 | 95 | TEveQuadSet(Form("Sector%i",nSector)), |
28a21457 | 96 | fTOFgeometry(new AliTOFGeometry()), |
8616d353 | 97 | fTOFarray(tofArray), |
98 | fTOFtree(0x0), | |
99 | fSector(nSector), | |
100 | fDx(0), fDy(0), fDz(0), | |
101 | fAutoTrans (kTRUE), | |
ee9d8eca | 102 | //fMinTime (0), fMaxTime (0), |
103 | fThreshold (5), fMaxVal (80), | |
8616d353 | 104 | fSectorID (nSector), |
ee9d8eca | 105 | //fPlateFlag0(kTRUE), fPlateFlag1(kTRUE), fPlateFlag2(kTRUE), fPlateFlag3(kTRUE), fPlateFlag4(kTRUE), |
106 | fPlateFlag (0x0), | |
107 | //fFrameColor(4), | |
108 | //fRnrFrame (kFALSE), | |
8616d353 | 109 | fGeoManager(localGeoManager) |
110 | { | |
111 | ||
112 | fPlateFlag = new Bool_t[5]; | |
113 | for (Int_t ii=0; ii<5; ii++) fPlateFlag[ii]=kTRUE; | |
114 | ||
115 | InitModule(); | |
116 | ||
117 | } | |
118 | /* ************************************************************************ */ | |
119 | ||
ee9d8eca | 120 | AliEveTOFSector::AliEveTOFSector(TGeoManager *localGeoManager, |
121 | Int_t nSector, | |
122 | TTree *tofTree) | |
123 | : | |
124 | TEveQuadSet(Form("Sector%i",nSector)), | |
28a21457 | 125 | fTOFgeometry(new AliTOFGeometry()), |
8616d353 | 126 | fTOFarray(0x0), |
127 | fTOFtree(tofTree), | |
128 | fSector(nSector), | |
129 | fDx(0), fDy(0), fDz(0), | |
130 | fAutoTrans (kTRUE), | |
ee9d8eca | 131 | //fMinTime (0), fMaxTime (0), |
132 | fThreshold (5), fMaxVal (80), | |
8616d353 | 133 | fSectorID (nSector), |
ee9d8eca | 134 | //fPlateFlag0(kTRUE), fPlateFlag1(kTRUE), fPlateFlag2(kTRUE), fPlateFlag3(kTRUE), fPlateFlag4(kTRUE), |
135 | fPlateFlag (0x0), | |
136 | //fFrameColor(4), | |
137 | //fRnrFrame (kFALSE), | |
8616d353 | 138 | fGeoManager(localGeoManager) |
139 | { | |
140 | ||
141 | fPlateFlag = new Bool_t[5]; | |
142 | for (Int_t ii=0; ii<5; ii++) fPlateFlag[ii]=kTRUE; | |
143 | ||
144 | InitModule(); | |
145 | ||
146 | } | |
147 | /* ************************************************************************ */ | |
148 | ||
d810d0de | 149 | AliEveTOFSector::~AliEveTOFSector() |
8616d353 | 150 | { |
151 | /* | |
152 | fGeoManager = 0x0; | |
153 | delete fGeoManager; | |
154 | ||
155 | delete fTOFarray; | |
156 | fTOFarray = 0x0; | |
157 | */ | |
158 | delete fPlateFlag; | |
159 | ||
160 | } | |
161 | ||
162 | /* ************************************************************************ */ | |
163 | /* | |
d810d0de | 164 | void AliEveTOFSector::SetDigitsInfo(AliEveTOFDigitsInfo* info) |
8616d353 | 165 | { |
166 | if(fInfo) fInfo->DecRefCount(); | |
167 | fInfo = info; | |
168 | if(fInfo) fInfo->IncRefCount(); | |
169 | ||
170 | } | |
171 | */ | |
172 | /* ************************************************************************ */ | |
d810d0de | 173 | void AliEveTOFSector::InitStatics() |
8616d353 | 174 | { |
175 | if (fgStaticInitDone) return; | |
176 | ||
177 | Float_t dx = 124.5; | |
178 | Float_t dz = 29.; | |
179 | Float_t dy = 370.6*2.; | |
84aff7a4 | 180 | fgTOFsectorFrameBox = new TEveFrameBox(); |
8616d353 | 181 | fgTOFsectorFrameBox->SetAABox(-dx*0.5, -dy*0.5, -dz*0.5, dx, dy, dz); |
ab73f3eb | 182 | fgTOFsectorFrameBox->SetFrameColor(Color_t(32)); |
55ee7431 | 183 | fgTOFsectorFrameBox->IncRefCount(); |
8616d353 | 184 | |
84aff7a4 | 185 | //fgTOFsectorPalette = new TEveRGBAPalette(0, 2048); // TOT |
ee9d8eca | 186 | //fgTOFsectorPalette->SetLimits(0, 2048); |
82dabae7 | 187 | //fgTOFsectorPalette = new TEveRGBAPalette(0, 8192); // TDC |
188 | fgTOFsectorPalette = new TEveRGBAPalette(0, 100000); // TDC | |
189 | fgTOFsectorPalette->SetOverflowAction(2); | |
190 | fgTOFsectorPalette->SetUnderflowAction(2); | |
55ee7431 | 191 | fgTOFsectorPalette->IncRefCount(); |
8616d353 | 192 | |
193 | fgStaticInitDone = kTRUE; | |
194 | } | |
195 | ||
196 | /* ************************************************************************ */ | |
d810d0de | 197 | void AliEveTOFSector::InitModule() |
8616d353 | 198 | { |
199 | ||
200 | fDx = fTOFgeometry->XPad()*fTOFgeometry->NpadX(); | |
201 | //fDy = fTOFgeometry->XPad()*fTOFgeometry->NpadX(); | |
202 | fDz = fTOFgeometry->ZPad()*fTOFgeometry->NpadZ(); | |
203 | ||
204 | if (!fgStaticInitDone) InitStatics(); | |
205 | ||
206 | SetFrame(fgTOFsectorFrameBox); | |
207 | SetPalette(fgTOFsectorPalette); | |
208 | //fFrame = fgTOFsectorFrameBox; | |
209 | //fPalette = fgTOFsectorPalette; | |
210 | ||
51346b82 | 211 | LoadQuads(); |
8616d353 | 212 | ComputeBBox(); |
213 | SetTrans(); | |
214 | ||
215 | } | |
216 | ||
217 | /* ************************************************************************ */ | |
d810d0de | 218 | void AliEveTOFSector::LoadQuads() |
8616d353 | 219 | { |
220 | ||
8616d353 | 221 | //Int_t n_col = gStyle->GetNumberOfColors(); |
222 | ||
223 | Int_t vol[5] = {fSectorID, -1, -1, -1, -1}; | |
224 | Int_t informations[4] = {-1, -1, -1, -1}; | |
225 | Int_t dummy[3] = {-1, -1, -1}; | |
226 | Int_t tdc = -1; | |
227 | Int_t tot = -1; | |
228 | ||
229 | Double_t **coord = new Double_t*[4]; | |
230 | for (Int_t ii=0; ii<4; ii++) coord[ii] = new Double_t[3]; | |
231 | ||
ee9d8eca | 232 | Reset(kQT_FreeQuad, kFALSE, 32); |
233 | ||
8616d353 | 234 | //printf(" fTOFarray->GetEntries() = %4i \n",fTOFarray->GetEntries()); |
235 | ||
236 | if (fTOFtree && !fTOFarray) { | |
237 | //printf("Hello world\n"); | |
238 | TClonesArray* digitsTOFnew = new TClonesArray("AliTOFdigit", 300); | |
239 | ||
240 | fTOFarray = new TClonesArray("AliTOFdigit", 300); | |
241 | TClonesArray &ldigits = *fTOFarray; | |
242 | Int_t newCounter = 0; | |
243 | ||
244 | AliTOFdigit *digs; | |
245 | ||
246 | fTOFtree->SetBranchAddress("TOF",&digitsTOFnew); | |
247 | fTOFtree->GetEntry(0); | |
248 | for (Int_t digitNumber=0; digitNumber<digitsTOFnew->GetEntries(); digitNumber++) { | |
249 | ||
250 | //if (digitNumber==digitsTOF->GetEntries()-1) printf(" Hello 4 -> %3i digit of %i \n", digitNumber+1, digitsTOF->GetEntries()); | |
51346b82 | 251 | |
8616d353 | 252 | digs = (AliTOFdigit*)digitsTOFnew->UncheckedAt(digitNumber); |
253 | ||
254 | if (digs->GetSector()!=fSectorID) continue; | |
255 | ||
256 | vol[1] = digs->GetPlate(); // Plate Number (0-4) | |
257 | vol[2] = digs->GetStrip(); // Strip Number (0-14/18) | |
ee9d8eca | 258 | vol[3] = digs->GetPadx(); // Pad Number in x direction (0-47) |
259 | vol[4] = digs->GetPadz(); // Pad Number in z direction (0-1) | |
8616d353 | 260 | |
261 | informations[0] = digs->GetTdc(); | |
262 | informations[1] = digs->GetAdc(); | |
263 | informations[2] = digs->GetToT(); | |
264 | informations[3] = digs->GetTdcND(); | |
265 | new (ldigits[newCounter++]) AliTOFdigit(dummy, vol, informations); | |
266 | } | |
267 | } | |
268 | ||
269 | AliTOFdigit *tofDigit; | |
270 | //printf(" 0x%lx\n",fTOFarray); | |
271 | for (Int_t ii=0; ii<fTOFarray->GetEntries(); ii++) { | |
272 | ||
273 | tofDigit = (AliTOFdigit*)fTOFarray->UncheckedAt(ii); | |
51346b82 | 274 | |
8616d353 | 275 | if (fPlateFlag[tofDigit->GetPlate()]) { |
276 | ||
277 | vol[1] = tofDigit->GetPlate(); | |
278 | vol[2] = tofDigit->GetStrip(); | |
279 | vol[3] = tofDigit->GetPadz(); | |
280 | vol[4] = tofDigit->GetPadx(); | |
281 | ||
282 | tot = tofDigit->GetToT(); | |
283 | tdc = tofDigit->GetTdc(); | |
284 | ||
285 | //if (fSector==4 && fPlate==2 && fStrip==0) printf(" %2i %1i\n", iPadX, iPadZ); | |
286 | //if (iPadX==23 || iPadX==24) printf(" %2i %1i %2i \n", fSector, fPlate, fStrip); | |
287 | //if (vol[0]==4 && vol[1]==2 && vol[2]==7) { | |
288 | ||
289 | for (Int_t kk=0; kk<4; kk++) for (Int_t jj=0; jj<3; jj++) coord[kk][jj]=0.; | |
290 | ||
291 | fTOFgeometry->DetToSectorRF(vol, coord); | |
292 | /* | |
293 | printf("\n"); | |
294 | printf(" %1i %2i, %f %f %f \n", vol[3], vol[4], coord[0][0], coord[0][1], coord[0][2]); | |
295 | printf(" %1i %2i, %f %f %f \n", vol[3], vol[4], coord[1][0], coord[1][1], coord[1][2]); | |
296 | printf(" %1i %2i, %f %f %f \n", vol[3], vol[4], coord[2][0], coord[2][1], coord[2][2]); | |
297 | printf(" %1i %2i, %f %f %f \n", vol[3], vol[4], coord[3][0], coord[3][1], coord[3][2]); | |
298 | */ | |
299 | Float_t vertices[12]={(Float_t)coord[0][0], (Float_t)coord[0][1], (Float_t)coord[0][2], | |
300 | (Float_t)coord[1][0], (Float_t)coord[1][1], (Float_t)coord[1][2], | |
301 | (Float_t)coord[2][0], (Float_t)coord[2][1], (Float_t)coord[2][2], | |
302 | (Float_t)coord[3][0], (Float_t)coord[3][1], (Float_t)coord[3][2]}; | |
51346b82 | 303 | |
8616d353 | 304 | AddQuad(vertices); |
305 | //AddQuad((Float_t*)coord); | |
306 | //AddQuad(coord[0], coord[1], coord[2], 2.5, 3.5); | |
307 | //AddQuad(-2.5*0.5, -3.5*0.5, 2.5, 3.5); | |
308 | ||
309 | // In principle could have color based on number of neigbours. We | |
310 | // can insert the time-of-flight value for each pad | |
ee9d8eca | 311 | |
8616d353 | 312 | //QuadValue((Int_t)tot); |
313 | QuadValue((Int_t)tdc); | |
314 | QuadId(tofDigit); | |
51346b82 | 315 | |
8616d353 | 316 | //} |
317 | } // closed if control on plates switched on | |
318 | } // closed loop on TOF sector digits | |
319 | ||
320 | RefitPlex(); | |
321 | ||
322 | fTOFarray = 0x0; | |
51346b82 | 323 | |
8616d353 | 324 | } |
325 | ||
326 | /* ************************************************************ */ | |
d810d0de | 327 | void AliEveTOFSector::SetTrans() |
8616d353 | 328 | { |
a15e6d7d | 329 | InitMainTrans(); |
330 | TEveTrans& t = RefMainTrans(); | |
8616d353 | 331 | |
332 | //Int_t det[5] = {fSector, -1, -1, -1, -1}; | |
333 | Char_t path[100]; | |
334 | ||
335 | Int_t localSector = fSector; | |
336 | if (!fAutoTrans) localSector = 4; | |
337 | ||
338 | //fTOFgeometry->GetVolumePath(det,path); | |
339 | fTOFgeometry->GetVolumePath(localSector,path); | |
340 | fGeoManager->cd(path); | |
341 | TGeoHMatrix global = *fGeoManager->GetCurrentMatrix(); | |
342 | Double_t *rotMat = global.GetRotationMatrix(); | |
343 | Double_t *tr = global.GetTranslation(); | |
344 | ||
a15e6d7d | 345 | t.SetBaseVec(1, rotMat[0], rotMat[3], rotMat[6]); |
346 | t.SetBaseVec(2, rotMat[1], rotMat[4], rotMat[7]); | |
347 | t.SetBaseVec(3, rotMat[2], rotMat[5], rotMat[8]); | |
348 | t.SetBaseVec(4, tr); | |
8616d353 | 349 | } |
350 | ||
351 | //----------------------------------------------------- | |
352 | ||
d810d0de | 353 | void AliEveTOFSector::SetSectorID(Int_t id) |
8616d353 | 354 | { |
355 | fSectorID = id; | |
356 | fSector = id; | |
357 | if (fAutoTrans) | |
358 | SetTrans(); // Force repositioning. | |
359 | ||
360 | LoadQuads(); | |
361 | } | |
362 | ||
363 | //----------------------------------------------------- | |
364 | ||
d810d0de | 365 | void AliEveTOFSector::SetPlate(Int_t nPlate, Bool_t r) |
8616d353 | 366 | { |
367 | ||
368 | fPlateFlag[nPlate] = r; | |
369 | ||
370 | //printf(" HELLO World ! %i %i %i \n", nPlate, r, fPlateFlag[nPlate]); | |
371 | } | |
372 | ||
57ffa5fb | 373 | /******************************************************************************/ |
8616d353 | 374 | |
d810d0de | 375 | void AliEveTOFSector::SetThreshold(Short_t t) |
8616d353 | 376 | { |
377 | fThreshold = TMath::Min(t, (Short_t)(fMaxVal - 1)); | |
378 | // ClearColorArray(); | |
379 | } | |
380 | ||
57ffa5fb | 381 | /******************************************************************************/ |
8616d353 | 382 | |
d810d0de | 383 | void AliEveTOFSector::SetMaxVal(Int_t mv) |
8616d353 | 384 | { |
385 | fMaxVal = TMath::Max(mv, (Int_t)(fThreshold + 1)); | |
386 | //ClearColorArray(); | |
387 | } | |
388 | ||
57ffa5fb | 389 | /******************************************************************************/ |
8616d353 | 390 | |
d810d0de | 391 | void AliEveTOFSector::DigitSelected(Int_t idx) |
8616d353 | 392 | { |
ee9d8eca | 393 | // Override control-click from TEveQuadSet |
8616d353 | 394 | |
84aff7a4 | 395 | DigitBase_t* qb = GetDigit(idx); |
8616d353 | 396 | TObject* obj = qb->fId.GetObject(); |
397 | AliTOFdigit* digs = dynamic_cast<AliTOFdigit*>(obj); | |
d810d0de | 398 | // printf("AliEveTOFSector::QuadSelected "); Print(); |
ee9d8eca | 399 | /* |
8616d353 | 400 | printf(" idx = %5i, value = %5d, obj = 0x%lx, digit = 0x%lx ", |
401 | idx, qb->fValue, (ULong_t)obj, (ULong_t)digs); | |
ee9d8eca | 402 | */ |
8616d353 | 403 | if (digs) |
ee9d8eca | 404 | printf("\n Sector = %2i Plate = %1i Strip = %2i PadZ = %1i PadX = %2i ToT = %3i Tof = %5i\n", |
405 | fSector , digs->GetPlate(), digs->GetStrip(), digs->GetPadz(), digs->GetPadx(), digs->GetToT(), digs->GetTdc()); | |
8616d353 | 406 | else printf("\n"); |
407 | ||
408 | } | |
409 | ||
57ffa5fb | 410 | /******************************************************************************/ |