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 = 5;
  12
  13   const Float_t zCenter[nPlanes]          = {  -50.0,   -58.0,   -66.0,   -72.0,   -76.0 };   // expressed in cm
  14
  15   const Float_t rMin[nPlanes]             = {   2.50,    2.50,    3.00,    3.50,    3.50 };   // expressed in cm
  16   const Float_t rMax[nPlanes]             = {   9.70,   11.00,   12.40,   13.40,   14.00 };   // expressed in cm
  17
  18   const Float_t pixelSizeX[nPlanes]       = { 25.e-4,  25.e-4,  25.e-4,  25.e-4,  25.e-4 };   // expressed in cm
  19   const Float_t pixelSizeY[nPlanes]       = { 25.e-4,  25.e-4,  25.e-4,  25.e-4,  25.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 };   // expressed in cm
  22   const Float_t thicknessSupport[nPlanes] = {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 };   // 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 };    // expressed in cm
  26   const Float_t equivalentSiliconBeforeFront[nPlanes] = {   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 };    // expressed in cm
  28
  29   const Float_t hasPixelRectangularPatternAlongY[nPlanes] = {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 //==========================================================================================================================================