]>
Commit | Line | Data |
---|---|---|
820b4d9e | 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 | //==================================================================================================================================================== | |
17 | // | |
18 | // Class for the description of the structure for the planes of the ALICE Muon Forward Tracker | |
19 | // | |
20 | // Contact author: antonio.uras@cern.ch | |
21 | // | |
22 | //==================================================================================================================================================== | |
23 | ||
24 | #include "TNamed.h" | |
25 | #include "THnSparse.h" | |
26 | #include "TClonesArray.h" | |
27 | #include "AliMFTPlane.h" | |
28 | #include "TAxis.h" | |
29 | #include "TPave.h" | |
30 | #include "TCanvas.h" | |
31 | #include "TH2D.h" | |
32 | #include "TEllipse.h" | |
33 | #include "TMath.h" | |
34 | #include "AliLog.h" | |
35 | ||
36 | ClassImp(AliMFTPlane) | |
37 | ||
38 | //==================================================================================================================================================== | |
39 | ||
40 | AliMFTPlane::AliMFTPlane(): | |
41 | TNamed(), | |
42 | fPlaneNumber(0), | |
43 | fZCenter(0), | |
44 | fRMinSupport(0), | |
45 | fRMax(0), | |
46 | fRMaxSupport(0), | |
47 | fPixelSizeX(0), | |
48 | fPixelSizeY(0), | |
49 | fThicknessActive(0), | |
50 | fThicknessSupport(0), | |
51 | fThicknessReadout(0), | |
52 | fZCenterActiveFront(0), | |
53 | fZCenterActiveBack(0), | |
54 | fEquivalentSilicon(0), | |
55 | fEquivalentSiliconBeforeFront(0), | |
56 | fEquivalentSiliconBeforeBack(0), | |
57 | fActiveElements(0), | |
58 | fReadoutElements(0), | |
59 | fSupportElements(0) | |
60 | { | |
61 | ||
62 | fPlaneNumber = -1; | |
63 | ||
64 | fActiveElements = new TClonesArray("THnSparseC"); | |
65 | fReadoutElements = new TClonesArray("THnSparseC"); | |
66 | fSupportElements = new TClonesArray("THnSparseC"); | |
67 | ||
68 | // default constructor | |
69 | ||
70 | } | |
71 | ||
72 | //==================================================================================================================================================== | |
73 | ||
74 | AliMFTPlane::AliMFTPlane(const Char_t *name, const Char_t *title): | |
75 | TNamed(name, title), | |
76 | fPlaneNumber(0), | |
77 | fZCenter(0), | |
78 | fRMinSupport(0), | |
79 | fRMax(0), | |
80 | fRMaxSupport(0), | |
81 | fPixelSizeX(0), | |
82 | fPixelSizeY(0), | |
83 | fThicknessActive(0), | |
84 | fThicknessSupport(0), | |
85 | fThicknessReadout(0), | |
86 | fZCenterActiveFront(0), | |
87 | fZCenterActiveBack(0), | |
88 | fEquivalentSilicon(0), | |
89 | fEquivalentSiliconBeforeFront(0), | |
90 | fEquivalentSiliconBeforeBack(0), | |
91 | fActiveElements(0), | |
92 | fReadoutElements(0), | |
93 | fSupportElements(0) | |
94 | { | |
95 | ||
96 | fPlaneNumber = -1; | |
97 | ||
98 | fActiveElements = new TClonesArray("THnSparseC"); | |
99 | fReadoutElements = new TClonesArray("THnSparseC"); | |
100 | fSupportElements = new TClonesArray("THnSparseC"); | |
101 | ||
102 | // default constructor | |
103 | ||
104 | } | |
105 | ||
106 | //==================================================================================================================================================== | |
107 | ||
108 | AliMFTPlane::AliMFTPlane(const AliMFTPlane& plane): | |
109 | TNamed(plane), | |
110 | fPlaneNumber(plane.fPlaneNumber), | |
111 | fZCenter(plane.fZCenter), | |
112 | fRMinSupport(plane.fRMinSupport), | |
113 | fRMax(plane.fRMax), | |
114 | fRMaxSupport(plane.fRMaxSupport), | |
115 | fPixelSizeX(plane.fPixelSizeX), | |
116 | fPixelSizeY(plane.fPixelSizeY), | |
117 | fThicknessActive(plane.fThicknessActive), | |
118 | fThicknessSupport(plane.fThicknessSupport), | |
119 | fThicknessReadout(plane.fThicknessReadout), | |
120 | fZCenterActiveFront(plane.fZCenterActiveFront), | |
121 | fZCenterActiveBack(plane.fZCenterActiveBack), | |
122 | fEquivalentSilicon(plane.fEquivalentSilicon), | |
123 | fEquivalentSiliconBeforeFront(plane.fEquivalentSiliconBeforeFront), | |
124 | fEquivalentSiliconBeforeBack(plane.fEquivalentSiliconBeforeBack), | |
125 | fActiveElements(plane.fActiveElements), | |
126 | fReadoutElements(plane.fReadoutElements), | |
127 | fSupportElements(plane.fSupportElements) | |
128 | { | |
129 | ||
130 | // copy constructor | |
131 | ||
132 | } | |
133 | ||
134 | //==================================================================================================================================================== | |
135 | ||
136 | AliMFTPlane& AliMFTPlane::operator=(const AliMFTPlane& plane) { | |
137 | ||
138 | // Assignment operator | |
139 | ||
140 | // check assignement to self | |
141 | if (this == &plane) return *this; | |
142 | ||
143 | // base class assignement | |
144 | TNamed::operator=(plane); | |
145 | ||
146 | // clear memory | |
147 | Clear(); | |
148 | ||
149 | fPlaneNumber = plane.fPlaneNumber; | |
150 | fZCenter = plane.fZCenter; | |
151 | fRMinSupport = plane.fRMinSupport; | |
152 | fRMax = plane.fRMax; | |
153 | fRMaxSupport = plane.fRMaxSupport; | |
154 | fPixelSizeX = plane.fPixelSizeX; | |
155 | fPixelSizeY = plane.fPixelSizeY; | |
156 | fThicknessActive = plane.fThicknessActive; | |
157 | fThicknessSupport = plane.fThicknessSupport; | |
158 | fThicknessReadout = plane.fThicknessReadout; | |
159 | fZCenterActiveFront = plane.fZCenterActiveFront; | |
160 | fZCenterActiveBack = plane.fZCenterActiveBack; | |
161 | fEquivalentSilicon = plane.fEquivalentSilicon; | |
162 | fEquivalentSiliconBeforeFront = plane.fEquivalentSiliconBeforeFront; | |
163 | fEquivalentSiliconBeforeBack = plane.fEquivalentSiliconBeforeBack; | |
164 | fActiveElements = plane.fActiveElements; | |
165 | fReadoutElements = plane.fReadoutElements; | |
166 | fSupportElements = plane.fSupportElements; | |
167 | ||
168 | return *this; | |
169 | ||
170 | } | |
171 | ||
172 | //==================================================================================================================================================== | |
173 | ||
174 | Bool_t AliMFTPlane::Init(Int_t planeNumber, | |
175 | Double_t zCenter, | |
176 | Double_t rMin, | |
177 | Double_t rMax, | |
178 | Double_t pixelSizeX, | |
179 | Double_t pixelSizeY, | |
180 | Double_t thicknessActive, | |
181 | Double_t thicknessSupport, | |
182 | Double_t thicknessReadout) { | |
183 | ||
184 | AliDebug(1, Form("Initializing Plane Structure for Plane %s", GetName())); | |
185 | ||
186 | fPlaneNumber = planeNumber; | |
187 | fZCenter = zCenter; | |
188 | fRMinSupport = rMin; | |
189 | fRMax = rMax; | |
190 | fPixelSizeX = pixelSizeX; | |
191 | fPixelSizeY = pixelSizeY; | |
192 | fThicknessActive = thicknessActive; | |
193 | fThicknessSupport = thicknessSupport; | |
194 | fThicknessReadout = thicknessReadout; | |
195 | ||
196 | fZCenterActiveFront = fZCenter - 0.5*fThicknessSupport - 0.5*fThicknessActive; | |
197 | fZCenterActiveBack = fZCenter + 0.5*fThicknessSupport + 0.5*fThicknessActive; | |
198 | ||
199 | if (fRMinSupport <= fRadiusMin) fRMinSupport = fRadiusMin; | |
200 | else { | |
201 | fRMinSupport = fRadiusMin + (fHeightActive-fActiveSuperposition) * Int_t((fRMinSupport-fRadiusMin)/(fHeightActive-fActiveSuperposition)); | |
202 | } | |
203 | ||
204 | if (fRMax < fRMinSupport+fHeightActive) fRMax = fRMinSupport + fHeightActive; | |
205 | ||
206 | fRMax = fRMinSupport + (fHeightActive-fActiveSuperposition) * | |
207 | (Int_t((fRMax-fRMinSupport-fHeightActive)/(fHeightActive-fActiveSuperposition))+1) + fHeightActive; | |
208 | ||
209 | fRMaxSupport = TMath::Sqrt(fHeightActive*(rMax-fHeightActive) + fRMax*fRMax) + fSupportExtMargin; | |
210 | ||
211 | return kTRUE; | |
212 | ||
213 | } | |
214 | ||
215 | //==================================================================================================================================================== | |
216 | ||
217 | Bool_t AliMFTPlane::CreateStructure() { | |
218 | ||
219 | Int_t nBins[3]={0}; | |
220 | Double_t minPosition[3]={0}, maxPosition[3]={0}; | |
221 | ||
222 | // ------------------- support element ------------------------------------------------- | |
223 | ||
224 | nBins[0] = 1; | |
225 | nBins[1] = 1; | |
226 | nBins[2] = 1; | |
227 | ||
228 | minPosition[0] = -1.*fRMaxSupport; | |
229 | minPosition[1] = -1.*fRMaxSupport; | |
230 | minPosition[2] = fZCenter - 0.5*fThicknessSupport; | |
231 | ||
232 | maxPosition[0] = +1.*fRMaxSupport; | |
233 | maxPosition[1] = +1.*fRMaxSupport; | |
234 | maxPosition[2] = fZCenter + 0.5*fThicknessSupport; | |
235 | ||
236 | new ((*fSupportElements)[fSupportElements->GetEntries()]) THnSparseC(Form("MFTSupportElemHist_%02d%03d", fPlaneNumber, fSupportElements->GetEntries()), | |
237 | Form("MFTSupportElemHist_%02d%03d", fPlaneNumber, fSupportElements->GetEntries()), | |
238 | 3, nBins, minPosition, maxPosition); | |
239 | ||
240 | // ------------------- det elements: active + readout ---------------------------------- | |
241 | ||
242 | Double_t lowEdgeActive = -1.*fRMax; | |
243 | Double_t supEdgeActive = lowEdgeActive + fHeightActive; | |
244 | Double_t zMin = 0.; | |
245 | Bool_t isFront = kTRUE; | |
246 | ||
247 | while (supEdgeActive < fRMax+0.01) { | |
248 | ||
249 | if (isFront) zMin = fZCenter - 0.5*fThicknessSupport - fThicknessActive; | |
250 | else zMin = fZCenter + 0.5*fThicknessSupport; | |
251 | ||
252 | Double_t extLimitAtLowEdgeActive = TMath::Sqrt((fRMax-TMath::Abs(lowEdgeActive)) * TMath::Abs(2*fRMax - (fRMax-TMath::Abs(lowEdgeActive)))); | |
253 | Double_t extLimitAtSupEdgeActive = TMath::Sqrt((fRMax-TMath::Abs(supEdgeActive)) * TMath::Abs(2*fRMax - (fRMax-TMath::Abs(supEdgeActive)))); | |
254 | ||
255 | // creating new det element: active + readout | |
256 | ||
257 | Double_t extLimitDetElem = TMath::Max(extLimitAtLowEdgeActive, extLimitAtSupEdgeActive); | |
258 | ||
259 | if (supEdgeActive<-1.*fRMinSupport+0.01 || lowEdgeActive>1.*fRMinSupport-0.01) { // single element covering the row | |
260 | ||
261 | nBins[0] = TMath::Nint(2.*extLimitDetElem/fPixelSizeX); | |
262 | nBins[1] = TMath::Nint(fHeightActive/fPixelSizeY); | |
263 | nBins[2] = 1; | |
264 | ||
265 | minPosition[0] = -1.*extLimitDetElem; | |
266 | minPosition[1] = lowEdgeActive; | |
267 | minPosition[2] = zMin; | |
268 | ||
269 | maxPosition[0] = +1.*extLimitDetElem; | |
270 | maxPosition[1] = supEdgeActive; | |
271 | maxPosition[2] = zMin+fThicknessActive; | |
272 | ||
273 | new ((*fActiveElements)[fActiveElements->GetEntries()]) THnSparseC(Form("MFTActiveElemHist_%02d%03d", fPlaneNumber, fActiveElements->GetEntries()), | |
274 | Form("MFTActiveElemHist_%02d%03d", fPlaneNumber, fActiveElements->GetEntries()), | |
275 | 3, nBins, minPosition, maxPosition); | |
276 | ||
277 | if (supEdgeActive>0.) { | |
278 | minPosition[1] = supEdgeActive; | |
279 | maxPosition[1] = supEdgeActive+fHeightReadout; | |
280 | } | |
281 | else { | |
282 | minPosition[1] = lowEdgeActive-fHeightReadout; | |
283 | maxPosition[1] = lowEdgeActive; | |
284 | } | |
285 | ||
286 | new ((*fReadoutElements)[fReadoutElements->GetEntries()]) THnSparseC(Form("MFTReadoutElemHist_%02d%03d", fPlaneNumber, fReadoutElements->GetEntries()), | |
287 | Form("MFTReadoutElemHist_%02d%03d", fPlaneNumber, fReadoutElements->GetEntries()), | |
288 | 3, nBins, minPosition, maxPosition); | |
289 | ||
290 | } | |
291 | ||
292 | else { // two elements covering the row | |
293 | ||
294 | Double_t intLimitAtLowEdge = 0., intLimitAtSupEdge = 0.; | |
295 | if (fRMinSupport-TMath::Abs(lowEdgeActive)>0.) intLimitAtLowEdge = TMath::Sqrt((fRMinSupport-TMath::Abs(lowEdgeActive)) * TMath::Abs(2*fRMinSupport - (fRMinSupport-TMath::Abs(lowEdgeActive)))); | |
296 | if (fRMinSupport-TMath::Abs(supEdgeActive)>0.) intLimitAtSupEdge = TMath::Sqrt((fRMinSupport-TMath::Abs(supEdgeActive)) * TMath::Abs(2*fRMinSupport - (fRMinSupport-TMath::Abs(supEdgeActive)))); | |
297 | Double_t intLimitDetElem = TMath::Max(intLimitAtLowEdge, intLimitAtSupEdge); | |
298 | ||
299 | nBins[0] = TMath::Nint((extLimitDetElem-intLimitDetElem)/fPixelSizeX); | |
300 | nBins[1] = TMath::Nint(fHeightActive/fPixelSizeY); | |
301 | nBins[2] = 1; | |
302 | ||
303 | // left element | |
304 | ||
305 | minPosition[0] = -1.*extLimitDetElem; | |
306 | minPosition[1] = lowEdgeActive; | |
307 | minPosition[2] = zMin; | |
308 | ||
309 | maxPosition[0] = -1.*intLimitDetElem; | |
310 | maxPosition[1] = supEdgeActive; | |
311 | maxPosition[2] = zMin+fThicknessActive; | |
312 | ||
313 | new ((*fActiveElements)[fActiveElements->GetEntries()]) THnSparseC(Form("MFTActiveElemHist_%02d%03d", fPlaneNumber, fActiveElements->GetEntries()), | |
314 | Form("MFTActiveElemHist_%02d%03d", fPlaneNumber, fActiveElements->GetEntries()), | |
315 | 3, nBins, minPosition, maxPosition); | |
316 | ||
317 | if (supEdgeActive>0.) { | |
318 | minPosition[1] = supEdgeActive; | |
319 | maxPosition[1] = supEdgeActive+fHeightReadout; | |
320 | } | |
321 | else { | |
322 | minPosition[1] = lowEdgeActive-fHeightReadout; | |
323 | maxPosition[1] = lowEdgeActive; | |
324 | } | |
325 | ||
326 | new ((*fReadoutElements)[fReadoutElements->GetEntries()]) THnSparseC(Form("MFTReadoutElemHist_%02d%03d", fPlaneNumber, fReadoutElements->GetEntries()), | |
327 | Form("MFTReadoutElemHist_%02d%03d", fPlaneNumber, fReadoutElements->GetEntries()), | |
328 | 3, nBins, minPosition, maxPosition); | |
329 | ||
330 | // right element | |
331 | ||
332 | minPosition[0] = +1.*intLimitDetElem; | |
333 | minPosition[1] = lowEdgeActive; | |
334 | minPosition[2] = zMin; | |
335 | ||
336 | maxPosition[0] = +1.*extLimitDetElem; | |
337 | maxPosition[1] = supEdgeActive; | |
338 | maxPosition[2] = zMin+fThicknessActive; | |
339 | ||
340 | new ((*fActiveElements)[fActiveElements->GetEntries()]) THnSparseC(Form("MFTActiveElemHist_%02d%03d", fPlaneNumber, fActiveElements->GetEntries()), | |
341 | Form("MFTActiveElemHist_%02d%03d", fPlaneNumber, fActiveElements->GetEntries()), | |
342 | 3, nBins, minPosition, maxPosition); | |
343 | ||
344 | if (supEdgeActive>0.) { | |
345 | minPosition[1] = supEdgeActive; | |
346 | maxPosition[1] = supEdgeActive+fHeightReadout; | |
347 | } | |
348 | else { | |
349 | minPosition[1] = lowEdgeActive-fHeightReadout; | |
350 | maxPosition[1] = lowEdgeActive; | |
351 | } | |
352 | ||
353 | new ((*fReadoutElements)[fReadoutElements->GetEntries()]) THnSparseC(Form("MFTReadoutElemHist_%02d%03d", fPlaneNumber, fReadoutElements->GetEntries()), | |
354 | Form("MFTReadoutElemHist_%02d%03d", fPlaneNumber, fReadoutElements->GetEntries()), | |
355 | 3, nBins, minPosition, maxPosition); | |
356 | ||
357 | } | |
358 | ||
359 | lowEdgeActive += fHeightActive - fActiveSuperposition; | |
360 | supEdgeActive = lowEdgeActive + fHeightActive; | |
361 | isFront = !isFront; | |
362 | ||
363 | } | |
364 | ||
365 | AliDebug(1, Form("Structure completed for MFT plane %s", GetName())); | |
366 | ||
367 | return kTRUE; | |
368 | ||
369 | } | |
370 | ||
371 | //==================================================================================================================================================== | |
372 | ||
373 | THnSparseC* AliMFTPlane::GetActiveElement(Int_t id) { | |
374 | ||
375 | if (id<0 || id>=GetNActiveElements()) return NULL; | |
376 | else return (THnSparseC*) fActiveElements->At(id); | |
377 | ||
378 | } | |
379 | ||
380 | //==================================================================================================================================================== | |
381 | ||
382 | THnSparseC* AliMFTPlane::GetReadoutElement(Int_t id) { | |
383 | ||
384 | if (id<0 || id>=GetNReadoutElements()) return NULL; | |
385 | else return (THnSparseC*) fReadoutElements->At(id); | |
386 | ||
387 | } | |
388 | ||
389 | //==================================================================================================================================================== | |
390 | ||
391 | THnSparseC* AliMFTPlane::GetSupportElement(Int_t id) { | |
392 | ||
393 | if (id<0 || id>=GetNSupportElements()) return NULL; | |
394 | else return (THnSparseC*) fSupportElements->At(id); | |
395 | ||
396 | } | |
397 | ||
398 | //==================================================================================================================================================== | |
399 | ||
400 | void AliMFTPlane::DrawPlane(Char_t *opt) { | |
401 | ||
402 | // ------------------- "FRONT" option ------------------ | |
403 | ||
404 | if (!strcmp(opt, "front")) { | |
405 | ||
406 | TCanvas *cnv = new TCanvas("cnv", GetName(), 900, 900); | |
407 | cnv->Draw(); | |
408 | ||
409 | // printf("Created Canvas\n"); | |
410 | ||
411 | TH2D *h = new TH2D("tmp", GetName(), | |
412 | 1, 1.1*GetSupportElement(0)->GetAxis(0)->GetXmin(), 1.1*GetSupportElement(0)->GetAxis(0)->GetXmax(), | |
413 | 1, 1.1*GetSupportElement(0)->GetAxis(1)->GetXmin(), 1.1*GetSupportElement(0)->GetAxis(1)->GetXmax()); | |
414 | h->SetXTitle("x [cm]"); | |
415 | h->SetYTitle("y [cm]"); | |
416 | h->Draw(); | |
417 | ||
418 | printf("Created hist\n"); | |
419 | ||
420 | TEllipse *supportExt = new TEllipse(0.0, 0.0, fRMaxSupport, fRMaxSupport); | |
421 | TEllipse *supportInt = new TEllipse(0.0, 0.0, fRMinSupport, fRMinSupport); | |
422 | supportExt->SetFillColor(kCyan-10); | |
423 | supportExt -> Draw("same"); | |
424 | supportInt -> Draw("same"); | |
425 | ||
426 | // printf("Created Ellipses\n"); | |
427 | ||
428 | for (Int_t iEl=0; iEl<GetNActiveElements(); iEl++) { | |
429 | // printf("Active element %d\n", iEl); | |
430 | if (!IsFront(GetActiveElement(iEl))) continue; | |
431 | TPave *pave = new TPave(GetActiveElement(iEl)->GetAxis(0)->GetXmin(), | |
432 | GetActiveElement(iEl)->GetAxis(1)->GetXmin(), | |
433 | GetActiveElement(iEl)->GetAxis(0)->GetXmax(), | |
434 | GetActiveElement(iEl)->GetAxis(1)->GetXmax(), 1); | |
435 | pave -> SetFillColor(kGreen); | |
436 | pave -> Draw("same"); | |
437 | } | |
438 | ||
439 | for (Int_t iEl=0; iEl<GetNReadoutElements(); iEl++) { | |
440 | // printf("Readout element %d\n", iEl); | |
441 | if (!IsFront(GetReadoutElement(iEl))) continue; | |
442 | TPave *pave = new TPave(GetReadoutElement(iEl)->GetAxis(0)->GetXmin(), | |
443 | GetReadoutElement(iEl)->GetAxis(1)->GetXmin(), | |
444 | GetReadoutElement(iEl)->GetAxis(0)->GetXmax(), | |
445 | GetReadoutElement(iEl)->GetAxis(1)->GetXmax(), 1); | |
446 | pave -> SetFillColor(kRed); | |
447 | pave -> Draw("same"); | |
448 | } | |
449 | ||
450 | } | |
451 | ||
452 | // ------------------- "BACK" option ------------------ | |
453 | ||
454 | else if (!strcmp(opt, "back")) { | |
455 | ||
456 | TCanvas *cnv = new TCanvas("cnv", GetName(), 900, 900); | |
457 | cnv->Draw(); | |
458 | ||
459 | TH2D *h = new TH2D("tmp", GetName(), | |
460 | 1, 1.1*GetSupportElement(0)->GetAxis(0)->GetXmin(), 1.1*GetSupportElement(0)->GetAxis(0)->GetXmax(), | |
461 | 1, 1.1*GetSupportElement(0)->GetAxis(1)->GetXmin(), 1.1*GetSupportElement(0)->GetAxis(1)->GetXmax()); | |
462 | h->SetXTitle("x [cm]"); | |
463 | h->SetYTitle("y [cm]"); | |
464 | h->Draw(); | |
465 | ||
466 | TEllipse *supportExt = new TEllipse(0.0, 0.0, fRMaxSupport, fRMaxSupport); | |
467 | TEllipse *supportInt = new TEllipse(0.0, 0.0, fRMinSupport, fRMinSupport); | |
468 | supportExt -> SetFillColor(kCyan-10); | |
469 | supportExt -> Draw("same"); | |
470 | supportInt -> Draw("same"); | |
471 | ||
472 | for (Int_t iEl=0; iEl<GetNActiveElements(); iEl++) { | |
473 | if (IsFront(GetActiveElement(iEl))) continue; | |
474 | TPave *pave = new TPave(GetActiveElement(iEl)->GetAxis(0)->GetXmin(), | |
475 | GetActiveElement(iEl)->GetAxis(1)->GetXmin(), | |
476 | GetActiveElement(iEl)->GetAxis(0)->GetXmax(), | |
477 | GetActiveElement(iEl)->GetAxis(1)->GetXmax(), 1); | |
478 | pave -> SetFillColor(kGreen); | |
479 | pave -> Draw("same"); | |
480 | } | |
481 | ||
482 | for (Int_t iEl=0; iEl<GetNReadoutElements(); iEl++) { | |
483 | if (IsFront(GetReadoutElement(iEl))) continue; | |
484 | TPave *pave = new TPave(GetReadoutElement(iEl)->GetAxis(0)->GetXmin(), | |
485 | GetReadoutElement(iEl)->GetAxis(1)->GetXmin(), | |
486 | GetReadoutElement(iEl)->GetAxis(0)->GetXmax(), | |
487 | GetReadoutElement(iEl)->GetAxis(1)->GetXmax(), 1); | |
488 | pave -> SetFillColor(kRed); | |
489 | pave -> Draw("same"); | |
490 | } | |
491 | ||
492 | } | |
493 | ||
494 | // ------------------- "BOTH" option ------------------ | |
495 | ||
496 | else if (!strcmp(opt, "both")) { | |
497 | ||
498 | TCanvas *cnv = new TCanvas("cnv", GetName(), 900, 900); | |
499 | cnv->Draw(); | |
500 | ||
501 | TH2D *h = new TH2D("tmp", GetName(), | |
502 | 1, 1.1*GetSupportElement(0)->GetAxis(0)->GetXmin(), 1.1*GetSupportElement(0)->GetAxis(0)->GetXmax(), | |
503 | 1, 1.1*GetSupportElement(0)->GetAxis(1)->GetXmin(), 1.1*GetSupportElement(0)->GetAxis(1)->GetXmax()); | |
504 | h->SetXTitle("x [cm]"); | |
505 | h->SetYTitle("y [cm]"); | |
506 | h->Draw(); | |
507 | ||
508 | TEllipse *supportExt = new TEllipse(0.0, 0.0, fRMaxSupport, fRMaxSupport); | |
509 | TEllipse *supportInt = new TEllipse(0.0, 0.0, fRMinSupport, fRMinSupport); | |
510 | supportExt -> SetFillColor(kCyan-10); | |
511 | supportExt -> Draw("same"); | |
512 | supportInt -> Draw("same"); | |
513 | ||
514 | for (Int_t iEl=0; iEl<GetNActiveElements(); iEl++) { | |
515 | if (IsFront(GetActiveElement(iEl)) && GetActiveElement(iEl)->GetAxis(0)->GetXmin()<0.) { | |
516 | TPave *pave = new TPave(GetActiveElement(iEl)->GetAxis(0)->GetXmin(), | |
517 | GetActiveElement(iEl)->GetAxis(1)->GetXmin(), | |
518 | TMath::Min(GetActiveElement(iEl)->GetAxis(0)->GetXmax(), 0.), | |
519 | GetActiveElement(iEl)->GetAxis(1)->GetXmax(), 1); | |
520 | pave -> SetFillColor(kGreen); | |
521 | pave -> Draw("same"); | |
522 | } | |
523 | else if (!IsFront(GetActiveElement(iEl)) && GetActiveElement(iEl)->GetAxis(0)->GetXmax()>0.) { | |
524 | TPave *pave = new TPave(TMath::Max(GetActiveElement(iEl)->GetAxis(0)->GetXmin(), 0.), | |
525 | GetActiveElement(iEl)->GetAxis(1)->GetXmin(), | |
526 | GetActiveElement(iEl)->GetAxis(0)->GetXmax(), | |
527 | GetActiveElement(iEl)->GetAxis(1)->GetXmax(), 1); | |
528 | pave -> SetFillColor(kGreen); | |
529 | pave -> Draw("same"); | |
530 | } | |
531 | } | |
532 | ||
533 | for (Int_t iEl=0; iEl<GetNReadoutElements(); iEl++) { | |
534 | if (IsFront(GetReadoutElement(iEl)) && GetReadoutElement(iEl)->GetAxis(0)->GetXmin()<0.) { | |
535 | TPave *pave = new TPave(GetReadoutElement(iEl)->GetAxis(0)->GetXmin(), | |
536 | GetReadoutElement(iEl)->GetAxis(1)->GetXmin(), | |
537 | TMath::Min(GetReadoutElement(iEl)->GetAxis(0)->GetXmax(), 0.), | |
538 | GetReadoutElement(iEl)->GetAxis(1)->GetXmax(), 1); | |
539 | pave -> SetFillColor(kRed); | |
540 | pave -> Draw("same"); | |
541 | } | |
542 | else if (!IsFront(GetReadoutElement(iEl)) && GetReadoutElement(iEl)->GetAxis(0)->GetXmax()>0.) { | |
543 | TPave *pave = new TPave(TMath::Max(GetReadoutElement(iEl)->GetAxis(0)->GetXmin(), 0.), | |
544 | GetReadoutElement(iEl)->GetAxis(1)->GetXmin(), | |
545 | GetReadoutElement(iEl)->GetAxis(0)->GetXmax(), | |
546 | GetReadoutElement(iEl)->GetAxis(1)->GetXmax(), 1); | |
547 | pave -> SetFillColor(kRed); | |
548 | pave -> Draw("same"); | |
549 | } | |
550 | } | |
551 | ||
552 | } | |
553 | ||
554 | // ------------------- "PROFILE" option ------------------ | |
555 | ||
556 | else if (!strcmp(opt, "profile")) { | |
557 | ||
558 | TCanvas *cnv = new TCanvas("cnv", GetName(), 300, 900); | |
559 | cnv->Draw(); | |
560 | ||
561 | TH2D *h = new TH2D("tmp", GetName(), | |
562 | 1, fZCenter-0.5, fZCenter+0.5, | |
563 | 1, 1.1*GetSupportElement(0)->GetAxis(1)->GetXmin(), 1.1*GetSupportElement(0)->GetAxis(1)->GetXmax()); | |
564 | h->SetXTitle("z [cm]"); | |
565 | h->SetYTitle("y [cm]"); | |
566 | h->Draw(); | |
567 | ||
568 | TPave *supportExt = new TPave(GetSupportElement(0)->GetAxis(2)->GetXmin(), -fRMaxSupport, | |
569 | GetSupportElement(0)->GetAxis(2)->GetXmax(), fRMaxSupport); | |
570 | TPave *supportInt = new TPave(GetSupportElement(0)->GetAxis(2)->GetXmin(), -fRMinSupport, | |
571 | GetSupportElement(0)->GetAxis(2)->GetXmax(), fRMinSupport); | |
572 | supportExt -> SetFillColor(kCyan-10); | |
573 | supportInt -> SetFillColor(kCyan-10); | |
574 | supportExt -> SetBorderSize(1); | |
575 | supportInt -> SetBorderSize(1); | |
576 | supportExt -> Draw("same"); | |
577 | supportInt -> Draw("same"); | |
578 | ||
579 | for (Int_t iEl=0; iEl<GetNActiveElements(); iEl++) { | |
580 | TPave * pave = 0; | |
581 | if (IsFront(GetActiveElement(iEl))) { | |
582 | pave = new TPave(GetActiveElement(iEl)->GetAxis(2)->GetXmax() - | |
583 | 5*(GetActiveElement(iEl)->GetAxis(2)->GetXmax()-GetActiveElement(iEl)->GetAxis(2)->GetXmin()), | |
584 | GetActiveElement(iEl)->GetAxis(1)->GetXmin(), | |
585 | GetActiveElement(iEl)->GetAxis(2)->GetXmax(), | |
586 | GetActiveElement(iEl)->GetAxis(1)->GetXmax(), 1); | |
587 | } | |
588 | else { | |
589 | pave = new TPave(GetActiveElement(iEl)->GetAxis(2)->GetXmin(), | |
590 | GetActiveElement(iEl)->GetAxis(1)->GetXmin(), | |
591 | GetActiveElement(iEl)->GetAxis(2)->GetXmin() + | |
592 | 5*(GetActiveElement(iEl)->GetAxis(2)->GetXmax()-GetActiveElement(iEl)->GetAxis(2)->GetXmin()), | |
593 | GetActiveElement(iEl)->GetAxis(1)->GetXmax(), 1); | |
594 | } | |
595 | pave -> SetFillColor(kGreen); | |
596 | pave -> Draw("same"); | |
597 | } | |
598 | ||
599 | for (Int_t iEl=0; iEl<GetNReadoutElements(); iEl++) { | |
600 | TPave *pave = 0; | |
601 | if (IsFront(GetReadoutElement(iEl))) { | |
602 | pave = new TPave(GetReadoutElement(iEl)->GetAxis(2)->GetXmax() - | |
603 | 5*(GetReadoutElement(iEl)->GetAxis(2)->GetXmax()-GetReadoutElement(iEl)->GetAxis(2)->GetXmin()), | |
604 | GetReadoutElement(iEl)->GetAxis(1)->GetXmin(), | |
605 | GetReadoutElement(iEl)->GetAxis(2)->GetXmax(), | |
606 | GetReadoutElement(iEl)->GetAxis(1)->GetXmax(), 1); | |
607 | } | |
608 | else { | |
609 | pave = new TPave(GetReadoutElement(iEl)->GetAxis(2)->GetXmin(), | |
610 | GetReadoutElement(iEl)->GetAxis(1)->GetXmin(), | |
611 | GetReadoutElement(iEl)->GetAxis(2)->GetXmin() + | |
612 | 5*(GetReadoutElement(iEl)->GetAxis(2)->GetXmax()-GetReadoutElement(iEl)->GetAxis(2)->GetXmin()), | |
613 | GetReadoutElement(iEl)->GetAxis(1)->GetXmax(), 1); | |
614 | } | |
615 | pave -> SetFillColor(kRed); | |
616 | pave -> Draw("same"); | |
617 | } | |
618 | ||
619 | } | |
620 | ||
621 | } | |
622 | ||
623 | //==================================================================================================================================================== | |
624 |