Changes to compile with Root6
[u/mrichter/AliRoot.git] / MFT / SetMFTGeometry.C
1 //==========================================================================================================================================
2
3 // Macro to create the config file (AliMFTGeometry.root) for the geometry of the ALICE Muon Forward Tracker
4 //
5 // Contact author: antonio.uras@cern.ch
6
7 //==========================================================================================================================================
8
9 void SetMFTGeometry() {
10
11   const Int_t nPlanes = 6;
12   
13   const Float_t zCenter[nPlanes]          = {  -45.5,   -49.3,   -53.1,   -56.9,   -68.8,    -76.8 };   // expressed in cm
14                                           
15   const Float_t rMin[nPlanes]             = {   2.50,    2.50,    2.50,    3.00,    3.50,     3.90 };   // expressed in cm  
16   const Float_t rMax[nPlanes]             = {  10.80,   11.50,   12.30,   13.00,   15.30,    16.90 };   // expressed in cm
17
18   const Float_t pixelSizeX[nPlanes]       = { 20.e-4,  20.e-4,  20.e-4,  20.e-4,  20.e-4,   20.e-4 };   // expressed in cm
19   const Float_t pixelSizeY[nPlanes]       = { 20.e-4,  20.e-4,  20.e-4,  20.e-4,  20.e-4,   20.e-4 };   // expressed in cm
20
21   const Float_t thicknessActive[nPlanes]  = {  50.e-4,   50.e-4,   50.e-4,   50.e-4,   50.e-4,   50.e-4 };   // expressed in cm
22   const Float_t thicknessSupport[nPlanes] = {2000.e-4, 2000.e-4, 2000.e-4, 2000.e-4, 2000.e-4, 2000.e-4 };   // expressed in cm
23   const Float_t thicknessReadout[nPlanes] = {  50.e-4,   50.e-4,   50.e-4,   50.e-4,   50.e-4,   50.e-4 };   // expressed in cm
24
25   const Float_t equivalentSilicon[nPlanes]            = { 400.e-4, 400.e-4, 400.e-4, 400.e-4, 400.e-4, 400.e-4 };    // expressed in cm
26   const Float_t equivalentSiliconBeforeFront[nPlanes] = {   0.e-4,   0.e-4,   0.e-4,   0.e-4,   0.e-4,   0.e-4 };    // expressed in cm
27   const Float_t equivalentSiliconBeforeBack[nPlanes]  = { 350.e-4, 350.e-4, 350.e-4, 350.e-4, 350.e-4, 350.e-4 };    // expressed in cm
28
29   const Float_t hasPixelRectangularPatternAlongY[nPlanes] = {0., 0., 0., 0., 0., 0.};
30                                          
31   TNtuple *geomMFT = new TNtuple("AliMFTGeometry", "ALICE MFT Geometry", "zCenter:rMin:rMax:pixelSizeX:pixelSizeY:thicknessActive:thicknessSupport:thicknessReadout:equivalentSilicon:equivalentSiliconBeforeFront:equivalentSiliconBeforeBack:hasPixelRectangularPatternAlongY");
32
33   for (Int_t iPlane=0; iPlane<nPlanes; iPlane++) geomMFT -> Fill(zCenter[iPlane],
34                                                                  rMin[iPlane],
35                                                                  rMax[iPlane],
36                                                                  pixelSizeX[iPlane],
37                                                                  pixelSizeY[iPlane],
38                                                                  thicknessActive[iPlane],
39                                                                  thicknessSupport[iPlane],
40                                                                  thicknessReadout[iPlane],
41                                                                  equivalentSilicon[iPlane],
42                                                                  equivalentSiliconBeforeFront[iPlane],
43                                                                  equivalentSiliconBeforeBack[iPlane],
44                                                                  hasPixelRectangularPatternAlongY[iPlane]);
45
46   TFile *fileGeomMFT = new TFile("AliMFTGeometry.root", "recreate");
47   geomMFT -> Write();
48   fileGeomMFT -> Close();
49
50 }
51
52 //==========================================================================================================================================