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,   -74.0,   -82.0 };   // expressed in cm
  14
  15   const Float_t rMin[nPlanes]             = {   2.00,    2.31,    2.66,    3.01,    3.36 };   // expressed in cm
  16   const Float_t rMax[nPlanes]             = {   9.70,   11.11,   12.52,   13.93,   15.34 };   // 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 };   // expressed in cm
  19   const Float_t pixelSizeY[nPlanes]       = { 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 };   // 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]            = { 300.e-4, 300.e-4, 300.e-4, 300.e-4, 300.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]  = { 250.e-4, 250.e-4, 250.e-4, 250.e-4, 250.e-4};    // expressed in cm
  28
  29   TNtuple *geomMFT = new TNtuple("AliMFTGeometry", "ALICE MFT Geometry", "zCenter:rMin:rMax:pixelSizeX:pixelSizeY:thicknessActive:thicknessSupport:thicknessReadout:equivalentSilicon:equivalentSiliconBeforeFront:equivalentSiliconBeforeBack");
  30
  31   for (Int_t iPlane=0; iPlane<nPlanes; iPlane++) geomMFT -> Fill(zCenter[iPlane],
  32                                                                  rMin[iPlane],
  33                                                                  rMax[iPlane],
  34                                                                  pixelSizeX[iPlane],
  35                                                                  pixelSizeY[iPlane],
  36                                                                  thicknessActive[iPlane],
  37                                                                  thicknessSupport[iPlane],
  38                                                                  thicknessReadout[iPlane],
  39                                                                  equivalentSilicon[iPlane],
  40                                                                  equivalentSiliconBeforeFront[iPlane],
  41                                                                  equivalentSiliconBeforeBack[iPlane]);
  42
  43   TFile *fileGeomMFT = new TFile("AliMFTGeometry.root", "recreate");
  44   geomMFT -> Write();
  45   fileGeomMFT -> Close();
  46
  47 }
  48
  49 //==========================================================================================================================================