2 // Main authors: Matevz Tadel & Alja Mrak-Tadel: 2006, 2007
4 /**************************************************************************
5 * Copyright(c) 1998-2008, ALICE Experiment at CERN, all rights reserved. *
6 * See http://aliceinfo.cern.ch/Offline/AliRoot/License.html for *
7 * full copyright notice. *
8 **************************************************************************/
10 #include "AliEveITSModule.h"
12 #include <AliITSgeomTGeo.h>
13 #include <AliITSsegmentationSPD.h>
14 #include <AliITSsegmentationSDD.h>
15 #include <AliITSsegmentationSSD.h>
17 #include <AliITSdigitSPD.h>
18 #include <AliITSdigitSDD.h>
19 #include <AliITSdigitSSD.h>
21 #include <TEveTrans.h>
22 #include <TClonesArray.h>
26 //______________________________________________________________________________
28 // Visualization of an ITS module.
30 ClassImp(AliEveITSModule)
32 Bool_t AliEveITSModule::fgStaticInitDone = kFALSE;
34 TEveFrameBox* AliEveITSModule::fgSPDFrameBox = 0;
35 TEveFrameBox* AliEveITSModule::fgSDDFrameBox = 0;
36 TEveFrameBox* AliEveITSModule::fgSSDFrameBox = 0;
38 TEveRGBAPalette* AliEveITSModule::fgSPDPalette = 0;
39 TEveRGBAPalette* AliEveITSModule::fgSDDPalette = 0;
40 TEveRGBAPalette* AliEveITSModule::fgSSDPalette = 0;
42 /******************************************************************************/
44 AliEveITSModule::AliEveITSModule(const Text_t* n, const Text_t* t) :
48 fLayer(-1), fLadder(-1), fDet(-1),
49 fDx(0), fDz(0), fDy(0)
54 AliEveITSModule::AliEveITSModule(Int_t gid, AliEveITSDigitsInfo* info) :
55 TEveQuadSet(Form("ITS module %d", gid)),
58 fLayer(-1), fLadder(-1), fDet(-1),
59 fDx(0), fDz(0), fDy(0)
61 // Constructor with module id and data-source.
67 AliEveITSModule::~AliEveITSModule()
71 if(fInfo) fInfo->DecRefCount();
74 /******************************************************************************/
76 void AliEveITSModule::InitStatics(AliEveITSDigitsInfo* info)
78 // Initialize static variables.
80 // Warning all sensor sizes are in microns, here we transform them
81 // to cm. In Eve half-lengths/widths are used, hence another 1/2.
83 if (fgStaticInitDone) return;
84 fgStaticInitDone = kTRUE;
87 Float_t dx = info->fSegSPD->Dx()*0.00005;
88 Float_t dz = info->fSegSPD->Dz()*0.00005;
90 fgSPDFrameBox = new TEveFrameBox();
91 fgSPDFrameBox->SetAAQuadXZ(-dx, 0, -dz, 2*dx, 2*dz);
92 fgSPDFrameBox->SetFrameColor((Color_t) 31);
93 fgSPDFrameBox->SetFrameFill(kTRUE);
94 fgSPDFrameBox->IncRefCount();
95 fgSPDPalette = new TEveRGBAPalette(info->fSPDMinVal,info->fSPDMaxVal);
96 fgSPDPalette->IncRefCount();
100 Float_t dx = info->fSegSDD->Dx()*0.0001;
101 Float_t dz = info->fSegSDD->Dz()*0.00005;
103 fgSDDFrameBox = new TEveFrameBox();
104 fgSDDFrameBox->SetAAQuadXZ(-dx, 0, -dz, 2*dx, 2*dz);
105 fgSDDFrameBox->SetFrameColor((Color_t) 32);
106 fgSDDFrameBox->SetFrameFill(kTRUE);
107 fgSDDFrameBox->IncRefCount();
108 fgSDDPalette = new TEveRGBAPalette(info->fSDDMinVal,info->fSDDMaxVal);
109 fgSDDPalette->SetLimits(0, info->fSDDHighLim); // Set proper ADC range.
110 fgSDDPalette->IncRefCount();
114 Float_t dx = info->fSegSSD->Dx()*0.00005;
115 Float_t dz = info->fSegSSD->Dz()*0.00005;
117 fgSSDFrameBox = new TEveFrameBox();
118 fgSSDFrameBox->SetAAQuadXZ(-dx, 0, -dz, 2*dx, 2*dz);
119 fgSSDFrameBox->SetFrameColor((Color_t) 33);
120 fgSSDFrameBox->SetFrameFill(kTRUE);
121 fgSSDFrameBox->IncRefCount();
122 fgSSDPalette = new TEveRGBAPalette(info->fSSDMinVal,info->fSSDMaxVal);
123 fgSSDPalette->SetLimits(0, info->fSSDHighLim); // Set proper ADC range.
124 fgSSDPalette->IncRefCount();
129 /******************************************************************************/
131 void AliEveITSModule::SetDigitsInfo(AliEveITSDigitsInfo* info)
133 // Set data and geometry source.
135 if (fInfo == info) return;
136 if (fInfo) fInfo->DecRefCount();
138 if (fInfo) fInfo->IncRefCount();
141 /******************************************************************************/
143 void AliEveITSModule::SetID(Int_t gid, Bool_t trans)
147 static const TEveException kEH("AliEveITSModule::SetID ");
150 throw(kEH + "AliEveITSDigitsInfo not set.");
152 Int_t firstSPD = AliITSgeomTGeo::GetModuleIndex(1,1,1);
153 Int_t lastSSD = AliITSgeomTGeo::GetNModules() - 1;
154 if (gid < firstSPD || gid > lastSSD)
156 throw(kEH + Form("%d is not valid. ID range from %d to %d", gid,
157 firstSPD, lastSSD ));
162 if (!fgStaticInitDone) {
165 fgSPDFrameBox->IncRefCount(this);
166 fgSPDPalette->IncRefCount();
168 fgSDDFrameBox->IncRefCount(this);
169 fgSDDPalette->IncRefCount();
171 fgSSDFrameBox->IncRefCount(this);
172 fgSSDPalette->IncRefCount();
175 AliITSgeomTGeo::GetModuleId(fID, fLayer, fLadder, fDet);
176 TString strLadder = "Ladder";
177 TString strSensor = "Sensor";
179 Int_t id, nsector, nstave, nladder, rest;
181 if (fID <= (AliITSgeomTGeo::GetModuleIndex(3,1,1) - 1))
185 SetFrame(fgSPDFrameBox);
186 SetPalette(fgSPDPalette);
188 symname += strLadder;
192 rest = fID - 8*nsector;
199 rest = id - 8*nsector;
202 if (rest < 4) nstave = 0;
207 fDx = fInfo->fSegSPD->Dx()*0.00005;
209 fDy = fInfo->fSegSPD->Dy()*0.00005;
212 else if (fID <= (AliITSgeomTGeo::GetModuleIndex(5,1,1) - 1))
216 SetFrame(fgSDDFrameBox);
217 SetPalette(fgSDDPalette);
219 symname += strSensor;
224 rest = id - 6*nladder;
230 rest = id - 8*nladder;
235 fDx = fInfo->fSegSDD->Dx()*0.0001;
236 fDz = fInfo->fSegSDD->Dz()*0.00005;
237 fDy = fInfo->fSegSDD->Dy()*0.00005;
244 SetFrame(fgSSDFrameBox);
245 SetPalette(fgSSDPalette);
247 symname += strSensor;
252 rest = id - nladder*22;
258 rest = id - nladder*25;
263 fInfo->fSegSSD->SetLayer(fLayer);
264 fDx = fInfo->fSegSSD->Dx()*0.00005;
265 fDz = fInfo->fSegSSD->Dz()*0.00005;
266 fDy = fInfo->fSegSSD->Dy()*0.00005;
277 void AliEveITSModule::LoadQuads()
279 // Read module data from source and create low-level objects for
280 // visualization - called quads.
282 TClonesArray *digits = fInfo->GetDigits(fID, fDetID);
283 Int_t ndigits = digits ? digits->GetEntriesFast() : 0;
285 Float_t x, z, dpx, dpz;
292 AliITSsegmentationSPD* seg = fInfo->fSegSPD;
294 Reset(kQT_RectangleXZFixedY, kFALSE, 32);
296 for (Int_t k=0; k<ndigits; ++k)
298 AliITSdigit *d = (AliITSdigit*) digits->UncheckedAt(k);
301 x = -0.5*seg->Dx() + i*seg->Dpx(0);
303 fInfo->GetSPDLocalZ(j, z);
304 dpx = seg->Dpx(i)*0.0001;
305 dpz = seg->Dpz(j)*0.0001;
307 AddQuad(x, z, dpx, dpz);
308 QuadValue(1); // In principle could have color based on number of neigbours
315 AliITSsegmentationSDD *seg = fInfo->fSegSDD;
317 Reset(kQT_RectangleXZFixedY, kFALSE, 32);
319 for (Int_t k=0; k<ndigits; ++k)
321 AliITSdigit* d = (AliITSdigit*) digits->UncheckedAt(k);
323 // if (d->GetSignal() > fgSDDThreshold)
327 seg->DetToLocal(i, j, x, z);
328 dpx = seg->Dpx(i)*0.0001;
329 dpz = seg->Dpz(j)*0.0001;
331 AddQuad(x-2*dpx, z - dpz*0.5, 4*dpx, dpz);
332 QuadValue(d->GetSignal());
340 AliITSsegmentationSSD* seg = fInfo->fSegSSD;
342 Reset(kQT_LineXZFixedY, kFALSE, 32);
344 Float_t ap, an; // positive/negative angles -> offsets
346 ap = TMath::Tan(ap) * fDz;
347 an = - TMath::Tan(an) * fDz;
349 for (Int_t k=0; k<ndigits; ++k)
351 AliITSdigit *d = (AliITSdigit*) digits->UncheckedAt(k);
352 // if(d->GetSignal() > fgSSDThreshold)
356 seg->DetToLocal(i,j,x,z);
358 Float_t a = ( d->GetCoord1() == 1) ? ap : an;
360 AddLine(x-a, -fDz, 2*a, 2*fDz);
361 QuadValue(d->GetSignal());
363 // printf("%3d -> %3d -> %8x\n", d->GetSignal(), ci, fQuads.back().color);
374 /******************************************************************************/
376 void AliEveITSModule::SetTrans()
378 // Set transformation matrix based on module id (use geometry to
379 // retrieve this information).
381 fMainTrans->SetFrom(*AliITSgeomTGeo::GetMatrix(fID));
384 /******************************************************************************/
386 void AliEveITSModule::DigitSelected(Int_t idx)
388 // Override secondary select (alt-click) from TEveQuadSet.
390 DigitBase_t *qb = GetDigit(idx);
391 TObject *obj = qb->fId.GetObject();
392 AliITSdigit *d = dynamic_cast<AliITSdigit*>(obj);
393 printf("AliEveITSModule::QuadSelected "); Print();
394 printf(" idx=%d, value=%d, obj=0x%lx, digit=0x%lx\n",
395 idx, qb->fValue, (ULong_t)obj, (ULong_t)d);
397 printf(" coord1=%3d coord2=%3d signal=%d\n",
398 d->GetCoord1(), d->GetCoord2(), d->GetSignal());
402 /******************************************************************************/
404 void AliEveITSModule::Print(Option_t* ) const
406 // Print object summary information.
408 printf("AliEveITSModule: ID %d, layer %d, ladder %d, det %d\n",
409 fID, fLayer, fLadder, fDetID);