AliRoot Core  3dc7879 (3dc7879)
AliMFTPlane.cxx
Go to the documentation of this file.
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 "TAxis.h"
28 #include "TPave.h"
29 #include "TCanvas.h"
30 #include "TH2D.h"
31 #include "TEllipse.h"
32 #include "TMath.h"
33 #include "AliLog.h"
34 #include "AliMFTConstants.h"
35 #include "AliMFTPlane.h"
36 
37 const Double_t AliMFTPlane::fActiveSuperposition = AliMFTConstants::fActiveSuperposition;
38 const Double_t AliMFTPlane::fHeightActive = AliMFTConstants::fHeightActive;
39 const Double_t AliMFTPlane::fHeightReadout = AliMFTConstants::fHeightReadout;
40 const Double_t AliMFTPlane::fSupportExtMargin = AliMFTConstants::fSupportExtMargin;
41 
42 ClassImp(AliMFTPlane)
43 
44 //====================================================================================================================================================
45 
46 AliMFTPlane::AliMFTPlane():
47  TNamed(),
48  fPlaneNumber(-1),
49  fZCenter(0),
50  fRMinSupport(0),
51  fRMax(0),
52  fRMaxSupport(0),
53  fPixelSizeX(0),
54  fPixelSizeY(0),
55  fThicknessActive(0),
56  fThicknessSupport(0),
57  fThicknessReadout(0),
58  fZCenterActiveFront(0),
59  fZCenterActiveBack(0),
60  fEquivalentSilicon(0),
61  fEquivalentSiliconBeforeFront(0),
62  fEquivalentSiliconBeforeBack(0),
63  fActiveElements(0),
64  fReadoutElements(0),
65  fSupportElements(0),
66  fHasPixelRectangularPatternAlongY(kFALSE),
67  fPlaneIsOdd(kFALSE)
68 {
69 
70  // default constructor
71 
72 }
73 
74 //====================================================================================================================================================
75 
76 AliMFTPlane::AliMFTPlane(const Char_t *name, const Char_t *title):
77  TNamed(name, title),
78  fPlaneNumber(-1),
79  fZCenter(0),
80  fRMinSupport(0),
81  fRMax(0),
82  fRMaxSupport(0),
83  fPixelSizeX(0),
84  fPixelSizeY(0),
85  fThicknessActive(0),
86  fThicknessSupport(0),
87  fThicknessReadout(0),
88  fZCenterActiveFront(0),
89  fZCenterActiveBack(0),
90  fEquivalentSilicon(0),
91  fEquivalentSiliconBeforeFront(0),
92  fEquivalentSiliconBeforeBack(0),
93  fActiveElements(0),
94  fReadoutElements(0),
95  fSupportElements(0),
96  fHasPixelRectangularPatternAlongY(kFALSE),
97  fPlaneIsOdd(kFALSE)
98 {
99 
100  // constructor
101  fActiveElements = new TClonesArray("THnSparseC");
102  fReadoutElements = new TClonesArray("THnSparseC");
103  fSupportElements = new TClonesArray("THnSparseC");
104  fActiveElements->SetOwner(kTRUE);
105  fReadoutElements->SetOwner(kTRUE);
106  fSupportElements->SetOwner(kTRUE);
107 
108 }
109 
110 //====================================================================================================================================================
111 
112 AliMFTPlane::AliMFTPlane(const AliMFTPlane& plane):
113  TNamed(plane),
114  fPlaneNumber(plane.fPlaneNumber),
115  fZCenter(plane.fZCenter),
116  fRMinSupport(plane.fRMinSupport),
117  fRMax(plane.fRMax),
118  fRMaxSupport(plane.fRMaxSupport),
119  fPixelSizeX(plane.fPixelSizeX),
120  fPixelSizeY(plane.fPixelSizeY),
121  fThicknessActive(plane.fThicknessActive),
122  fThicknessSupport(plane.fThicknessSupport),
123  fThicknessReadout(plane.fThicknessReadout),
124  fZCenterActiveFront(plane.fZCenterActiveFront),
125  fZCenterActiveBack(plane.fZCenterActiveBack),
126  fEquivalentSilicon(plane.fEquivalentSilicon),
127  fEquivalentSiliconBeforeFront(plane.fEquivalentSiliconBeforeFront),
128  fEquivalentSiliconBeforeBack(plane.fEquivalentSiliconBeforeBack),
129  fActiveElements(0),
130  fReadoutElements(0),
131  fSupportElements(0),
132  fHasPixelRectangularPatternAlongY(plane.fHasPixelRectangularPatternAlongY),
133  fPlaneIsOdd(plane.fPlaneIsOdd)
134 {
135 
136  // copy constructor
137  fActiveElements = new TClonesArray(*(plane.fActiveElements));
138  fActiveElements -> SetOwner(kTRUE);
139  fReadoutElements = new TClonesArray(*(plane.fReadoutElements));
140  fReadoutElements -> SetOwner(kTRUE);
141  fSupportElements = new TClonesArray(*(plane.fSupportElements));
142  fSupportElements -> SetOwner(kTRUE);
143 
144 
145 }
146 
147 //====================================================================================================================================================
148 
149 AliMFTPlane::~AliMFTPlane() {
150 
151  AliInfo("Delete AliMFTPlane");
152  if(fActiveElements) fActiveElements->Delete();
153  delete fActiveElements;
154  if(fReadoutElements) fReadoutElements->Delete();
155  delete fReadoutElements;
156  if(fSupportElements) fSupportElements->Delete();
157  delete fSupportElements;
158 
159 }
160 
161 //====================================================================================================================================================
162 
163 void AliMFTPlane::Clear(const Option_t* /*opt*/) {
164 
165  AliInfo("Clear AliMFTPlane");
166  if(fActiveElements) fActiveElements->Delete();
167  delete fActiveElements; fActiveElements=NULL;
168  if(fReadoutElements) fReadoutElements->Delete();
169  delete fReadoutElements; fReadoutElements=NULL;
170  if(fSupportElements) fSupportElements->Delete();
171  delete fSupportElements; fSupportElements=NULL;
172 
173 }
174 
175 //====================================================================================================================================================
176 
177 AliMFTPlane& AliMFTPlane::operator=(const AliMFTPlane& plane) {
178 
179  // Assignment operator
180 
181  // check assignement to self
182  if (this != &plane) {
183 
184  // base class assignement
185  TNamed::operator=(plane);
186 
187  // clear memory
188  Clear("");
189 
190  fPlaneNumber = plane.fPlaneNumber;
191  fZCenter = plane.fZCenter;
192  fRMinSupport = plane.fRMinSupport;
193  fRMax = plane.fRMax;
194  fRMaxSupport = plane.fRMaxSupport;
195  fPixelSizeX = plane.fPixelSizeX;
196  fPixelSizeY = plane.fPixelSizeY;
197  fThicknessActive = plane.fThicknessActive;
198  fThicknessSupport = plane.fThicknessSupport;
199  fThicknessReadout = plane.fThicknessReadout;
200  fZCenterActiveFront = plane.fZCenterActiveFront;
201  fZCenterActiveBack = plane.fZCenterActiveBack;
202  fEquivalentSilicon = plane.fEquivalentSilicon;
203  fEquivalentSiliconBeforeFront = plane.fEquivalentSiliconBeforeFront;
204  fEquivalentSiliconBeforeBack = plane.fEquivalentSiliconBeforeBack;
205  fActiveElements = new TClonesArray(*(plane.fActiveElements));
206  fActiveElements -> SetOwner(kTRUE);
207  fReadoutElements = new TClonesArray(*(plane.fReadoutElements));
208  fReadoutElements -> SetOwner(kTRUE);
209  fSupportElements = new TClonesArray(*(plane.fSupportElements));
210  fSupportElements -> SetOwner(kTRUE);
211  fHasPixelRectangularPatternAlongY = plane.fHasPixelRectangularPatternAlongY;
212  fPlaneIsOdd = plane.fPlaneIsOdd;
213 
214  }
215 
216  return *this;
217 
218 }
219 
220 //====================================================================================================================================================
221 
222 Bool_t AliMFTPlane::Init(Int_t planeNumber,
223  Double_t zCenter,
224  Double_t rMin,
225  Double_t rMax,
226  Double_t pixelSizeX,
227  Double_t pixelSizeY,
228  Double_t thicknessActive,
229  Double_t thicknessSupport,
230  Double_t thicknessReadout,
231  Bool_t hasPixelRectangularPatternAlongY) {
232 
233  AliDebug(1, Form("Initializing Plane Structure for Plane %s", GetName()));
234 
235  fPlaneNumber = planeNumber;
236  fZCenter = zCenter;
237  fRMinSupport = rMin;
238  fRMax = rMax;
239  fPixelSizeX = pixelSizeX;
240  fPixelSizeY = pixelSizeY;
241  fThicknessActive = thicknessActive;
242  fThicknessSupport = thicknessSupport;
243  fThicknessReadout = thicknessReadout;
244 
245  fHasPixelRectangularPatternAlongY = hasPixelRectangularPatternAlongY;
246 
247  fZCenterActiveFront = fZCenter - 0.5*fThicknessSupport - 0.5*fThicknessActive;
248  fZCenterActiveBack = fZCenter + 0.5*fThicknessSupport + 0.5*fThicknessActive;
249 
250  if (fRMax < fRMinSupport+fHeightActive) fRMax = fRMinSupport + fHeightActive;
251 
252  Int_t nLaddersWithinPipe = Int_t(fRMinSupport/(fHeightActive-fActiveSuperposition));
253  if (fRMinSupport-nLaddersWithinPipe*(fHeightActive-fActiveSuperposition) > 0.5*(fHeightActive-2*fActiveSuperposition)) fPlaneIsOdd = kTRUE;
254  else fPlaneIsOdd = kFALSE;
255 
256  fRMax = fRMinSupport + (fHeightActive-fActiveSuperposition) *
257  (Int_t((fRMax-fRMinSupport-fHeightActive)/(fHeightActive-fActiveSuperposition))+1) + fHeightActive;
258 
259  fRMaxSupport = TMath::Sqrt(fHeightActive*(2.*rMax-fHeightActive) + fRMax*fRMax) + fSupportExtMargin;
260 
261  return kTRUE;
262 
263 }
264 
265 //====================================================================================================================================================
266 
267 Bool_t AliMFTPlane::CreateStructure() {
268 
269  Int_t nBins[3]={0};
270  Double_t minPosition[3]={0}, maxPosition[3]={0};
271 
272  // ------------------- det elements: active + readout ----------------------------------
273 
274  Double_t lowEdgeActive = -1.*fRMax;
275  Double_t supEdgeActive = lowEdgeActive + fHeightActive;
276  Double_t zMinFront = fZCenter - 0.5*fThicknessSupport - fThicknessActive;
277  Double_t zMinBack = fZCenter + 0.5*fThicknessSupport;
278  Double_t zMin = 0.;
279  Bool_t isFront = kTRUE;
280 
281  while (lowEdgeActive < 0) {
282 
283  Double_t extLimitAtLowEdgeActive = TMath::Sqrt((fRMax-TMath::Abs(lowEdgeActive)) * TMath::Abs(2*fRMax - (fRMax-TMath::Abs(lowEdgeActive))));
284  Double_t extLimitAtSupEdgeActive = TMath::Sqrt((fRMax-TMath::Abs(supEdgeActive)) * TMath::Abs(2*fRMax - (fRMax-TMath::Abs(supEdgeActive))));
285 
286  // creating new det element: active + readout
287 
288  Double_t extLimitDetElem = TMath::Max(extLimitAtLowEdgeActive, extLimitAtSupEdgeActive);
289 
290  if (supEdgeActive<-1.*fRMinSupport+0.01 || lowEdgeActive>1.*fRMinSupport-0.01) { // single element covering the row
291 
292  nBins[0] = TMath::Nint(2.*extLimitDetElem/fPixelSizeX);
293  nBins[1] = TMath::Nint(fHeightActive/fPixelSizeY);
294  nBins[2] = 1;
295 
296  // element below the pipe
297 
298  if (isFront) zMin = zMinFront;
299  else zMin = zMinBack;
300 
301  minPosition[0] = -1.*extLimitDetElem;
302  minPosition[1] = lowEdgeActive;
303  minPosition[2] = zMin;
304 
305  maxPosition[0] = +1.*extLimitDetElem;
306  maxPosition[1] = supEdgeActive;
307  maxPosition[2] = zMin+fThicknessActive;
308 
309  new ((*fActiveElements)[fActiveElements->GetEntries()]) THnSparseC(Form("MFTActiveElemHist_%02d%03d", fPlaneNumber, fActiveElements->GetEntries()),
310  Form("MFTActiveElemHist_%02d%03d", fPlaneNumber, fActiveElements->GetEntries()),
311  3, nBins, minPosition, maxPosition);
312 
313  minPosition[1] = lowEdgeActive-fHeightReadout;
314  maxPosition[1] = lowEdgeActive;
315 
316  new ((*fReadoutElements)[fReadoutElements->GetEntries()]) THnSparseC(Form("MFTReadoutElemHist_%02d%03d", fPlaneNumber, fReadoutElements->GetEntries()),
317  Form("MFTReadoutElemHist_%02d%03d", fPlaneNumber, fReadoutElements->GetEntries()),
318  3, nBins, minPosition, maxPosition);
319 
320  // specular element above the pipe
321 
322  if (fPlaneIsOdd) {
323  if (isFront) zMin = zMinBack;
324  else zMin = zMinFront;
325  }
326 
327  minPosition[0] = -1.*extLimitDetElem;
328  minPosition[1] = -1.*supEdgeActive;
329  minPosition[2] = zMin;
330 
331  maxPosition[0] = +1.*extLimitDetElem;
332  maxPosition[1] = -1.*lowEdgeActive;
333  maxPosition[2] = zMin+fThicknessActive;
334 
335  new ((*fActiveElements)[fActiveElements->GetEntries()]) THnSparseC(Form("MFTActiveElemHist_%02d%03d", fPlaneNumber, fActiveElements->GetEntries()),
336  Form("MFTActiveElemHist_%02d%03d", fPlaneNumber, fActiveElements->GetEntries()),
337  3, nBins, minPosition, maxPosition);
338 
339  minPosition[1] = -1.*lowEdgeActive;
340  maxPosition[1] = -1.*(lowEdgeActive-fHeightReadout);
341 
342  new ((*fReadoutElements)[fReadoutElements->GetEntries()]) THnSparseC(Form("MFTReadoutElemHist_%02d%03d", fPlaneNumber, fReadoutElements->GetEntries()),
343  Form("MFTReadoutElemHist_%02d%03d", fPlaneNumber, fReadoutElements->GetEntries()),
344  3, nBins, minPosition, maxPosition);
345 
346  }
347 
348  else { // two elements covering the row
349 
350  Double_t intLimitAtLowEdge = 0., intLimitAtSupEdge = 0.;
351  if (fRMinSupport-TMath::Abs(lowEdgeActive)>0.) intLimitAtLowEdge = TMath::Sqrt((fRMinSupport-TMath::Abs(lowEdgeActive)) * TMath::Abs(2*fRMinSupport - (fRMinSupport-TMath::Abs(lowEdgeActive))));
352  if (fRMinSupport-TMath::Abs(supEdgeActive)>0.) intLimitAtSupEdge = TMath::Sqrt((fRMinSupport-TMath::Abs(supEdgeActive)) * TMath::Abs(2*fRMinSupport - (fRMinSupport-TMath::Abs(supEdgeActive))));
353  Double_t intLimitDetElem = TMath::Max(intLimitAtLowEdge, intLimitAtSupEdge);
354 
355  nBins[0] = TMath::Nint((extLimitDetElem-intLimitDetElem)/fPixelSizeX);
356  nBins[1] = TMath::Nint(fHeightActive/fPixelSizeY);
357  nBins[2] = 1;
358 
359  // left element: y < 0
360 
361  if (isFront) zMin = zMinFront;
362  else zMin = zMinBack;
363 
364  minPosition[0] = -1.*extLimitDetElem;
365  minPosition[1] = lowEdgeActive;
366  minPosition[2] = zMin;
367 
368  maxPosition[0] = -1.*intLimitDetElem;
369  maxPosition[1] = supEdgeActive;
370  maxPosition[2] = zMin+fThicknessActive;
371 
372  new ((*fActiveElements)[fActiveElements->GetEntries()]) THnSparseC(Form("MFTActiveElemHist_%02d%03d", fPlaneNumber, fActiveElements->GetEntries()),
373  Form("MFTActiveElemHist_%02d%03d", fPlaneNumber, fActiveElements->GetEntries()),
374  3, nBins, minPosition, maxPosition);
375 
376  minPosition[1] = lowEdgeActive-fHeightReadout;
377  maxPosition[1] = lowEdgeActive;
378 
379  new ((*fReadoutElements)[fReadoutElements->GetEntries()]) THnSparseC(Form("MFTReadoutElemHist_%02d%03d", fPlaneNumber, fReadoutElements->GetEntries()),
380  Form("MFTReadoutElemHist_%02d%03d", fPlaneNumber, fReadoutElements->GetEntries()),
381  3, nBins, minPosition, maxPosition);
382 
383  // left element: y > 0
384 
385  if (supEdgeActive < 0.5*fHeightActive) {
386 
387  if (fPlaneIsOdd) {
388  if (isFront) zMin = zMinBack;
389  else zMin = zMinFront;
390  }
391 
392  minPosition[0] = -1.*extLimitDetElem;
393  minPosition[1] = -1.*supEdgeActive;
394  minPosition[2] = zMin;
395 
396  maxPosition[0] = -1.*intLimitDetElem;
397  maxPosition[1] = -1.*lowEdgeActive;
398  maxPosition[2] = zMin+fThicknessActive;
399 
400  new ((*fActiveElements)[fActiveElements->GetEntries()]) THnSparseC(Form("MFTActiveElemHist_%02d%03d", fPlaneNumber, fActiveElements->GetEntries()),
401  Form("MFTActiveElemHist_%02d%03d", fPlaneNumber, fActiveElements->GetEntries()),
402  3, nBins, minPosition, maxPosition);
403 
404  minPosition[1] = -1.*lowEdgeActive;
405  maxPosition[1] = -1.*(lowEdgeActive-fHeightReadout);
406 
407  new ((*fReadoutElements)[fReadoutElements->GetEntries()]) THnSparseC(Form("MFTReadoutElemHist_%02d%03d", fPlaneNumber, fReadoutElements->GetEntries()),
408  Form("MFTReadoutElemHist_%02d%03d", fPlaneNumber, fReadoutElements->GetEntries()),
409  3, nBins, minPosition, maxPosition);
410 
411  }
412 
413  // right element: y < 0
414 
415  if (isFront) zMin = zMinFront;
416  else zMin = zMinBack;
417 
418  minPosition[0] = +1.*intLimitDetElem;
419  minPosition[1] = lowEdgeActive;
420  minPosition[2] = zMin;
421 
422  maxPosition[0] = +1.*extLimitDetElem;
423  maxPosition[1] = supEdgeActive;
424  maxPosition[2] = zMin+fThicknessActive;
425 
426  new ((*fActiveElements)[fActiveElements->GetEntries()]) THnSparseC(Form("MFTActiveElemHist_%02d%03d", fPlaneNumber, fActiveElements->GetEntries()),
427  Form("MFTActiveElemHist_%02d%03d", fPlaneNumber, fActiveElements->GetEntries()),
428  3, nBins, minPosition, maxPosition);
429 
430  minPosition[1] = lowEdgeActive-fHeightReadout;
431  maxPosition[1] = lowEdgeActive;
432 
433  new ((*fReadoutElements)[fReadoutElements->GetEntries()]) THnSparseC(Form("MFTReadoutElemHist_%02d%03d", fPlaneNumber, fReadoutElements->GetEntries()),
434  Form("MFTReadoutElemHist_%02d%03d", fPlaneNumber, fReadoutElements->GetEntries()),
435  3, nBins, minPosition, maxPosition);
436 
437  // right element: y > 0
438 
439  if (supEdgeActive < 0.5*fHeightActive) {
440 
441  if (fPlaneIsOdd) {
442  if (isFront) zMin = zMinBack;
443  else zMin = zMinFront;
444  }
445 
446  minPosition[0] = +1.*intLimitDetElem;
447  minPosition[1] = -1.*supEdgeActive;
448  minPosition[2] = zMin;
449 
450  maxPosition[0] = +1.*extLimitDetElem;
451  maxPosition[1] = -1.*lowEdgeActive;
452  maxPosition[2] = zMin+fThicknessActive;
453 
454  new ((*fActiveElements)[fActiveElements->GetEntries()]) THnSparseC(Form("MFTActiveElemHist_%02d%03d", fPlaneNumber, fActiveElements->GetEntries()),
455  Form("MFTActiveElemHist_%02d%03d", fPlaneNumber, fActiveElements->GetEntries()),
456  3, nBins, minPosition, maxPosition);
457 
458  minPosition[1] = -1.*lowEdgeActive;
459  maxPosition[1] = -1.*(lowEdgeActive-fHeightReadout);
460 
461  new ((*fReadoutElements)[fReadoutElements->GetEntries()]) THnSparseC(Form("MFTReadoutElemHist_%02d%03d", fPlaneNumber, fReadoutElements->GetEntries()),
462  Form("MFTReadoutElemHist_%02d%03d", fPlaneNumber, fReadoutElements->GetEntries()),
463  3, nBins, minPosition, maxPosition);
464 
465  }
466 
467  }
468 
469  lowEdgeActive += fHeightActive - fActiveSuperposition;
470  supEdgeActive = lowEdgeActive + fHeightActive;
471  isFront = !isFront;
472 
473  }
474 
475  // ------------------- support element -------------------------------------------------
476 
477  nBins[0] = 1;
478  nBins[1] = 1;
479  nBins[2] = 1;
480 
481  minPosition[0] = -1.*fRMaxSupport;
482  minPosition[1] = -1.*fRMaxSupport;
483  minPosition[2] = fZCenter - 0.5*fThicknessSupport;
484 
485  maxPosition[0] = +1.*fRMaxSupport;
486  maxPosition[1] = +1.*fRMaxSupport;
487  maxPosition[2] = fZCenter + 0.5*fThicknessSupport;
488 
489  new ((*fSupportElements)[fSupportElements->GetEntries()]) THnSparseC(Form("MFTSupportElemHist_%02d%03d", fPlaneNumber, fSupportElements->GetEntries()),
490  Form("MFTSupportElemHist_%02d%03d", fPlaneNumber, fSupportElements->GetEntries()),
491  3, nBins, minPosition, maxPosition);
492 
493  // --------------------------------------------------------------------------------------
494 
495  AliDebug(1, Form("Structure completed for MFT plane %s", GetName()));
496 
497  return kTRUE;
498 
499 }
500 
501 //====================================================================================================================================================
502 
503 THnSparseC* AliMFTPlane::GetActiveElement(Int_t id) {
504 
505  if (id<0 || id>=GetNActiveElements()) return NULL;
506  else return (THnSparseC*) fActiveElements->At(id);
507 
508 }
509 
510 //====================================================================================================================================================
511 
512 THnSparseC* AliMFTPlane::GetReadoutElement(Int_t id) {
513 
514  if (id<0 || id>=GetNReadoutElements()) return NULL;
515  else return (THnSparseC*) fReadoutElements->At(id);
516 
517 }
518 
519 //====================================================================================================================================================
520 
521 THnSparseC* AliMFTPlane::GetSupportElement(Int_t id) {
522 
523  if (id<0 || id>=GetNSupportElements()) return NULL;
524  else return (THnSparseC*) fSupportElements->At(id);
525 
526 }
527 
528 //====================================================================================================================================================
529 
530 void AliMFTPlane::DrawPlane(Option_t *opt) {
531 
532  // ------------------- "FRONT" option ------------------
533 
534  if (!strcmp(opt, "front")) {
535 
536  TCanvas *cnv = new TCanvas("cnv", GetName(), 900, 900);
537  cnv->Draw();
538 
539  TH2D *h = new TH2D("tmp", GetName(),
540  1, 1.1*GetSupportElement(0)->GetAxis(0)->GetXmin(), 1.1*GetSupportElement(0)->GetAxis(0)->GetXmax(),
541  1, 1.1*GetSupportElement(0)->GetAxis(1)->GetXmin(), 1.1*GetSupportElement(0)->GetAxis(1)->GetXmax());
542  h->SetXTitle("x [cm]");
543  h->SetYTitle("y [cm]");
544  h->Draw();
545 
546  AliInfo("Created hist");
547 
548  TEllipse *supportExt = new TEllipse(0.0, 0.0, fRMaxSupport, fRMaxSupport);
549  TEllipse *supportInt = new TEllipse(0.0, 0.0, fRMinSupport, fRMinSupport);
550  supportExt->SetFillColor(kCyan-10);
551  supportExt -> Draw("same");
552  supportInt -> Draw("same");
553 
554  for (Int_t iEl=0; iEl<GetNActiveElements(); iEl++) {
555  if (!IsFront(GetActiveElement(iEl))) continue;
556  TPave *pave = new TPave(GetActiveElement(iEl)->GetAxis(0)->GetXmin(),
557  GetActiveElement(iEl)->GetAxis(1)->GetXmin(),
558  GetActiveElement(iEl)->GetAxis(0)->GetXmax(),
559  GetActiveElement(iEl)->GetAxis(1)->GetXmax(), 1);
560  pave -> SetFillColor(kGreen);
561  pave -> Draw("same");
562  }
563 
564  for (Int_t iEl=0; iEl<GetNReadoutElements(); iEl++) {
565  if (!IsFront(GetReadoutElement(iEl))) continue;
566  TPave *pave = new TPave(GetReadoutElement(iEl)->GetAxis(0)->GetXmin(),
567  GetReadoutElement(iEl)->GetAxis(1)->GetXmin(),
568  GetReadoutElement(iEl)->GetAxis(0)->GetXmax(),
569  GetReadoutElement(iEl)->GetAxis(1)->GetXmax(), 1);
570  pave -> SetFillColor(kRed);
571  pave -> Draw("same");
572  }
573 
574  }
575 
576  // ------------------- "BACK" option ------------------
577 
578  else if (!strcmp(opt, "back")) {
579 
580  TCanvas *cnv = new TCanvas("cnv", GetName(), 900, 900);
581  cnv->Draw();
582 
583  TH2D *h = new TH2D("tmp", GetName(),
584  1, 1.1*GetSupportElement(0)->GetAxis(0)->GetXmin(), 1.1*GetSupportElement(0)->GetAxis(0)->GetXmax(),
585  1, 1.1*GetSupportElement(0)->GetAxis(1)->GetXmin(), 1.1*GetSupportElement(0)->GetAxis(1)->GetXmax());
586  h->SetXTitle("x [cm]");
587  h->SetYTitle("y [cm]");
588  h->Draw();
589 
590  TEllipse *supportExt = new TEllipse(0.0, 0.0, fRMaxSupport, fRMaxSupport);
591  TEllipse *supportInt = new TEllipse(0.0, 0.0, fRMinSupport, fRMinSupport);
592  supportExt -> SetFillColor(kCyan-10);
593  supportExt -> Draw("same");
594  supportInt -> Draw("same");
595 
596  for (Int_t iEl=0; iEl<GetNActiveElements(); iEl++) {
597  if (IsFront(GetActiveElement(iEl))) continue;
598  TPave *pave = new TPave(GetActiveElement(iEl)->GetAxis(0)->GetXmin(),
599  GetActiveElement(iEl)->GetAxis(1)->GetXmin(),
600  GetActiveElement(iEl)->GetAxis(0)->GetXmax(),
601  GetActiveElement(iEl)->GetAxis(1)->GetXmax(), 1);
602  pave -> SetFillColor(kGreen);
603  pave -> Draw("same");
604  }
605 
606  for (Int_t iEl=0; iEl<GetNReadoutElements(); iEl++) {
607  if (IsFront(GetReadoutElement(iEl))) continue;
608  TPave *pave = new TPave(GetReadoutElement(iEl)->GetAxis(0)->GetXmin(),
609  GetReadoutElement(iEl)->GetAxis(1)->GetXmin(),
610  GetReadoutElement(iEl)->GetAxis(0)->GetXmax(),
611  GetReadoutElement(iEl)->GetAxis(1)->GetXmax(), 1);
612  pave -> SetFillColor(kRed);
613  pave -> Draw("same");
614  }
615 
616  }
617 
618  // ------------------- "BOTH" option ------------------
619 
620  else if (!strcmp(opt, "both")) {
621 
622  TCanvas *cnv = new TCanvas("cnv", GetName(), 900, 900);
623  cnv->Draw();
624 
625  TH2D *h = new TH2D("tmp", GetName(),
626  1, 1.1*GetSupportElement(0)->GetAxis(0)->GetXmin(), 1.1*GetSupportElement(0)->GetAxis(0)->GetXmax(),
627  1, 1.1*GetSupportElement(0)->GetAxis(1)->GetXmin(), 1.1*GetSupportElement(0)->GetAxis(1)->GetXmax());
628  h->SetXTitle("x [cm]");
629  h->SetYTitle("y [cm]");
630  h->Draw();
631 
632  TEllipse *supportExt = new TEllipse(0.0, 0.0, fRMaxSupport, fRMaxSupport);
633  TEllipse *supportInt = new TEllipse(0.0, 0.0, fRMinSupport, fRMinSupport);
634  supportExt -> SetFillColor(kCyan-10);
635  supportExt -> Draw("same");
636  supportInt -> Draw("same");
637 
638  for (Int_t iEl=0; iEl<GetNActiveElements(); iEl++) {
639  if (IsFront(GetActiveElement(iEl)) && GetActiveElement(iEl)->GetAxis(0)->GetXmin()<0.) {
640  TPave *pave = new TPave(GetActiveElement(iEl)->GetAxis(0)->GetXmin(),
641  GetActiveElement(iEl)->GetAxis(1)->GetXmin(),
642  TMath::Min(GetActiveElement(iEl)->GetAxis(0)->GetXmax(), 0.),
643  GetActiveElement(iEl)->GetAxis(1)->GetXmax(), 1);
644  pave -> SetFillColor(kGreen);
645  pave -> Draw("same");
646  }
647  else if (!IsFront(GetActiveElement(iEl)) && GetActiveElement(iEl)->GetAxis(0)->GetXmax()>0.) {
648  TPave *pave = new TPave(TMath::Max(GetActiveElement(iEl)->GetAxis(0)->GetXmin(), 0.),
649  GetActiveElement(iEl)->GetAxis(1)->GetXmin(),
650  GetActiveElement(iEl)->GetAxis(0)->GetXmax(),
651  GetActiveElement(iEl)->GetAxis(1)->GetXmax(), 1);
652  pave -> SetFillColor(kGreen);
653  pave -> Draw("same");
654  }
655  }
656 
657  for (Int_t iEl=0; iEl<GetNReadoutElements(); iEl++) {
658  if (IsFront(GetReadoutElement(iEl)) && GetReadoutElement(iEl)->GetAxis(0)->GetXmin()<0.) {
659  TPave *pave = new TPave(GetReadoutElement(iEl)->GetAxis(0)->GetXmin(),
660  GetReadoutElement(iEl)->GetAxis(1)->GetXmin(),
661  TMath::Min(GetReadoutElement(iEl)->GetAxis(0)->GetXmax(), 0.),
662  GetReadoutElement(iEl)->GetAxis(1)->GetXmax(), 1);
663  pave -> SetFillColor(kRed);
664  pave -> Draw("same");
665  }
666  else if (!IsFront(GetReadoutElement(iEl)) && GetReadoutElement(iEl)->GetAxis(0)->GetXmax()>0.) {
667  TPave *pave = new TPave(TMath::Max(GetReadoutElement(iEl)->GetAxis(0)->GetXmin(), 0.),
668  GetReadoutElement(iEl)->GetAxis(1)->GetXmin(),
669  GetReadoutElement(iEl)->GetAxis(0)->GetXmax(),
670  GetReadoutElement(iEl)->GetAxis(1)->GetXmax(), 1);
671  pave -> SetFillColor(kRed);
672  pave -> Draw("same");
673  }
674  }
675 
676  }
677 
678  // ------------------- "PROFILE" option ------------------
679 
680  else if (!strcmp(opt, "profile")) {
681 
682  TCanvas *cnv = new TCanvas("cnv", GetName(), 300, 900);
683  cnv->Draw();
684 
685  TH2D *h = new TH2D("tmp", GetName(),
686  1, fZCenter-0.5, fZCenter+0.5,
687  1, 1.1*GetSupportElement(0)->GetAxis(1)->GetXmin(), 1.1*GetSupportElement(0)->GetAxis(1)->GetXmax());
688  h->SetXTitle("z [cm]");
689  h->SetYTitle("y [cm]");
690  h->Draw();
691 
692  TPave *supportExt = new TPave(GetSupportElement(0)->GetAxis(2)->GetXmin(), -fRMaxSupport,
693  GetSupportElement(0)->GetAxis(2)->GetXmax(), fRMaxSupport);
694  TPave *supportInt = new TPave(GetSupportElement(0)->GetAxis(2)->GetXmin(), -fRMinSupport,
695  GetSupportElement(0)->GetAxis(2)->GetXmax(), fRMinSupport);
696  supportExt -> SetFillColor(kCyan-10);
697  supportInt -> SetFillColor(kCyan-10);
698  supportExt -> SetBorderSize(1);
699  supportInt -> SetBorderSize(1);
700  supportExt -> Draw("same");
701  supportInt -> Draw("same");
702 
703  for (Int_t iEl=0; iEl<GetNActiveElements(); iEl++) {
704  TPave * pave = 0;
705  if (IsFront(GetActiveElement(iEl))) {
706  pave = new TPave(GetActiveElement(iEl)->GetAxis(2)->GetXmax() -
707  5*(GetActiveElement(iEl)->GetAxis(2)->GetXmax()-GetActiveElement(iEl)->GetAxis(2)->GetXmin()),
708  GetActiveElement(iEl)->GetAxis(1)->GetXmin(),
709  GetActiveElement(iEl)->GetAxis(2)->GetXmax(),
710  GetActiveElement(iEl)->GetAxis(1)->GetXmax(), 1);
711  }
712  else {
713  pave = new TPave(GetActiveElement(iEl)->GetAxis(2)->GetXmin(),
714  GetActiveElement(iEl)->GetAxis(1)->GetXmin(),
715  GetActiveElement(iEl)->GetAxis(2)->GetXmin() +
716  5*(GetActiveElement(iEl)->GetAxis(2)->GetXmax()-GetActiveElement(iEl)->GetAxis(2)->GetXmin()),
717  GetActiveElement(iEl)->GetAxis(1)->GetXmax(), 1);
718  }
719  pave -> SetFillColor(kGreen);
720  pave -> Draw("same");
721  }
722 
723  for (Int_t iEl=0; iEl<GetNReadoutElements(); iEl++) {
724  TPave *pave = 0;
725  if (IsFront(GetReadoutElement(iEl))) {
726  pave = new TPave(GetReadoutElement(iEl)->GetAxis(2)->GetXmax() -
727  5*(GetReadoutElement(iEl)->GetAxis(2)->GetXmax()-GetReadoutElement(iEl)->GetAxis(2)->GetXmin()),
728  GetReadoutElement(iEl)->GetAxis(1)->GetXmin(),
729  GetReadoutElement(iEl)->GetAxis(2)->GetXmax(),
730  GetReadoutElement(iEl)->GetAxis(1)->GetXmax(), 1);
731  }
732  else {
733  pave = new TPave(GetReadoutElement(iEl)->GetAxis(2)->GetXmin(),
734  GetReadoutElement(iEl)->GetAxis(1)->GetXmin(),
735  GetReadoutElement(iEl)->GetAxis(2)->GetXmin() +
736  5*(GetReadoutElement(iEl)->GetAxis(2)->GetXmax()-GetReadoutElement(iEl)->GetAxis(2)->GetXmin()),
737  GetReadoutElement(iEl)->GetAxis(1)->GetXmax(), 1);
738  }
739  pave -> SetFillColor(kRed);
740  pave -> Draw("same");
741  }
742 
743  }
744 
745 }
746 
747 //====================================================================================================================================================
748 
749 Int_t AliMFTPlane::GetNumberOfChips(Option_t *opt) {
750 
751  Int_t nChips = 0;
752 
753  if (!strcmp(opt, "front")) {
754  for (Int_t iEl=0; iEl<GetNActiveElements(); iEl++) {
755  if (!IsFront(GetActiveElement(iEl))) continue;
756  Double_t length = GetActiveElement(iEl)->GetAxis(0)->GetXmax() - GetActiveElement(iEl)->GetAxis(0)->GetXmin();
757  nChips += Int_t (length/AliMFTConstants::fWidthChip) + 1;
758  }
759  }
760 
761  else if (!strcmp(opt, "back")) {
762  for (Int_t iEl=0; iEl<GetNActiveElements(); iEl++) {
763  if (IsFront(GetActiveElement(iEl))) continue;
764  Double_t length = GetActiveElement(iEl)->GetAxis(0)->GetXmax() - GetActiveElement(iEl)->GetAxis(0)->GetXmin();
765  nChips += Int_t (length/AliMFTConstants::fWidthChip) + 1;
766  }
767  }
768 
769  return nChips;
770 
771 }
772 
773 //====================================================================================================================================================
AliMFTPlane::fZCenter
Double_t fZCenter
Definition: AliMFTPlane.h:101
AliMFTPlane::GetNReadoutElements
Int_t GetNReadoutElements() const
Definition: AliMFTPlane.h:55
AliMFTPlane::fEquivalentSiliconBeforeBack
Double_t fEquivalentSiliconBeforeBack
Definition: AliMFTPlane.h:102
AliMFTPlane::fRMaxSupport
Double_t fRMaxSupport
Definition: AliMFTPlane.h:101
AliMFTPlane::fRMax
Double_t fRMax
Definition: AliMFTPlane.h:101
AliMFTPlane::IsFront
Bool_t IsFront(THnSparseC *element) const
Definition: AliMFTPlane.h:66
AliMFTPlane::fThicknessActive
Double_t fThicknessActive
Definition: AliMFTPlane.h:101
AliMFTPlane::fEquivalentSilicon
Double_t fEquivalentSilicon
Definition: AliMFTPlane.h:102
AliMFTPlane::fReadoutElements
TClonesArray * fReadoutElements
Definition: AliMFTPlane.h:104
AliMFTPlane::GetNSupportElements
Int_t GetNSupportElements() const
Definition: AliMFTPlane.h:56
AliMFTPlane::DrawPlane
void DrawPlane(Option_t *opt="")
Definition: AliMFTPlane.cxx:530
AliMFTPlane::fPlaneIsOdd
Bool_t fPlaneIsOdd
Definition: AliMFTPlane.h:106
AliMFTPlane::fThicknessSupport
Double_t fThicknessSupport
Definition: AliMFTPlane.h:101
AliMFTPlane::Init
Bool_t Init(Int_t planeNumber, Double_t zCenter, Double_t rMin, Double_t rMax, Double_t pixelSizeX, Double_t pixelSizeY, Double_t thicknessActive, Double_t thicknessSupport, Double_t thicknessReadout, Bool_t hasPixelRectangularPatternAlongY)
Definition: AliMFTPlane.cxx:222
AliMFTPlane::fZCenterActiveFront
Double_t fZCenterActiveFront
Definition: AliMFTPlane.h:102
AliMFTConstants::fHeightActive
static const Double_t fHeightActive
height of the active elements
Definition: AliMFTConstants.h:69
AliLog.h
AliInfo
#define AliInfo(message)
Definition: AliLog.h:484
AliMFTConstants::fActiveSuperposition
static const Double_t fActiveSuperposition
Definition: AliMFTConstants.h:67
AliMFTPlane::GetNActiveElements
Int_t GetNActiveElements() const
Definition: AliMFTPlane.h:54
AliMFTPlane::fPixelSizeY
Double_t fPixelSizeY
Definition: AliMFTPlane.h:101
AliMFTPlane::GetSupportElement
THnSparseC * GetSupportElement(Int_t id)
Definition: AliMFTPlane.cxx:521
AliMFTPlane::CreateStructure
Bool_t CreateStructure()
Definition: AliMFTPlane.cxx:267
AliMFTPlane::fSupportExtMargin
static const Double_t fSupportExtMargin
Definition: AliMFTPlane.h:97
AliMFTPlane::GetReadoutElement
THnSparseC * GetReadoutElement(Int_t id)
Definition: AliMFTPlane.cxx:512
AliMFTPlane::Clear
virtual void Clear(const Option_t *)
Definition: AliMFTPlane.cxx:163
AliMFTPlane::operator=
AliMFTPlane & operator=(const AliMFTPlane &source)
Definition: AliMFTPlane.cxx:177
AliMFTPlane::fPixelSizeX
Double_t fPixelSizeX
Definition: AliMFTPlane.h:101
AliMFTPlane::fThicknessReadout
Double_t fThicknessReadout
Definition: AliMFTPlane.h:101
nBins
Int_t nBins
Definition: AliNDFormulaBrowser.cxx:60
AliMFTConstants::fHeightReadout
static const Double_t fHeightReadout
height of the readout elements attached to the active ones
Definition: AliMFTConstants.h:70
AliMFTConstants.h
AliMFTPlane::fActiveElements
TClonesArray * fActiveElements
Definition: AliMFTPlane.h:104
AliDebug
#define AliDebug(logLevel, message)
Definition: AliLog.h:300
AliMFTPlane::fHasPixelRectangularPatternAlongY
Bool_t fHasPixelRectangularPatternAlongY
Definition: AliMFTPlane.h:106
Draw
gr Draw("same*")
AliMFTPlane::GetNumberOfChips
Int_t GetNumberOfChips(Option_t *opt)
Definition: AliMFTPlane.cxx:749
AliMFTPlane::fEquivalentSiliconBeforeFront
Double_t fEquivalentSiliconBeforeFront
Definition: AliMFTPlane.h:102
AliMFTPlane::fPlaneNumber
Int_t fPlaneNumber
Definition: AliMFTPlane.h:99
AliMFTPlane::fActiveSuperposition
static const Double_t fActiveSuperposition
Definition: AliMFTPlane.h:92
AliMFTPlane::fRMinSupport
Double_t fRMinSupport
Definition: AliMFTPlane.h:101
AliMFTPlane::fHeightActive
static const Double_t fHeightActive
Definition: AliMFTPlane.h:94
AliMFTPlane.h
AliMFTPlane::fHeightReadout
static const Double_t fHeightReadout
Definition: AliMFTPlane.h:95
AliMFTConstants::fSupportExtMargin
static const Double_t fSupportExtMargin
Definition: AliMFTConstants.h:73
AliMFTPlane::~AliMFTPlane
virtual ~AliMFTPlane()
Definition: AliMFTPlane.cxx:149
AliMFTPlane::fSupportElements
TClonesArray * fSupportElements
Definition: AliMFTPlane.h:104
AliMFTConstants::fWidthChip
static const Double_t fWidthChip
expressed in cm
Definition: AliMFTConstants.h:79
AliMFTPlane::fZCenterActiveBack
Double_t fZCenterActiveBack
Definition: AliMFTPlane.h:102
AliMFTPlane::GetActiveElement
THnSparseC * GetActiveElement(Int_t id)
Definition: AliMFTPlane.cxx:503