Geometry and simulation file modified
[u/mrichter/AliRoot.git] / MFT / SetMFTGeometry.C
CommitLineData
820b4d9e 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
9void SetMFTGeometry() {
10
190a8304 11 const Int_t nPlanes = 6;
820b4d9e 12
190a8304 13 const Float_t zCenter[nPlanes] = { -45.5, -49.3, -53.1, -56.9, -68.8, -76.8 }; // expressed in cm
820b4d9e 14
190a8304 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
58b93f36 17
190a8304 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
820b4d9e 20
190a8304 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
820b4d9e 24
190a8304 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
820b4d9e 28
190a8304 29 const Float_t hasPixelRectangularPatternAlongY[nPlanes] = {0., 0., 0., 0., 0., 0.};
58b93f36 30
bcaf50eb 31 TNtuple *geomMFT = new TNtuple("AliMFTGeometry", "ALICE MFT Geometry", "zCenter:rMin:rMax:pixelSizeX:pixelSizeY:thicknessActive:thicknessSupport:thicknessReadout:equivalentSilicon:equivalentSiliconBeforeFront:equivalentSiliconBeforeBack:hasPixelRectangularPatternAlongY");
820b4d9e 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],
bcaf50eb 43 equivalentSiliconBeforeBack[iPlane],
44 hasPixelRectangularPatternAlongY[iPlane]);
820b4d9e 45
46 TFile *fileGeomMFT = new TFile("AliMFTGeometry.root", "recreate");
47 geomMFT -> Write();
48 fileGeomMFT -> Close();
49
50}
51
52//==========================================================================================================================================