fix AliHLTGlobalEsdConverterComponent
[u/mrichter/AliRoot.git] / MFT / AliMFTPlane.cxx
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 (supEdgeActive < 0.5*(fHeightActive+fHeightReadout)) {
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-1.1*(0.5*fThicknessSupport+fThicknessActive), fZCenter+1.1*(0.5*fThicknessSupport+fThicknessActive),
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 //====================================================================================================================================================