[u/mrichter/AliRoot.git] / MFT / SetMFTGeometry.C

//==========================================================================================================================================

// Macro to create the config file (AliMFTGeometry.root) for the geometry of the ALICE Muon Forward Tracker
//
// Contact author: antonio.uras@cern.ch

//==========================================================================================================================================

void SetMFTGeometry() {

  const Int_t nPlanes = 6;
  
  const Float_t zCenter[nPlanes]          = {  -45.5,   -49.3,   -53.1,   -56.9,   -68.8,    -76.8 };   // expressed in cm
				          
  const Float_t rMin[nPlanes]             = {   2.50,    2.50,    2.50,    3.00,    3.50,     3.90 };   // expressed in cm  
  const Float_t rMax[nPlanes]             = {  10.80,   11.50,   12.30,   13.00,   15.30,    16.90 };   // expressed in cm

  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
  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

  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
  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
  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

  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
  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
  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

  const Float_t hasPixelRectangularPatternAlongY[nPlanes] = {0., 0., 0., 0., 0., 0.};
				         
  TNtuple *geomMFT = new TNtuple("AliMFTGeometry", "ALICE MFT Geometry", "zCenter:rMin:rMax:pixelSizeX:pixelSizeY:thicknessActive:thicknessSupport:thicknessReadout:equivalentSilicon:equivalentSiliconBeforeFront:equivalentSiliconBeforeBack:hasPixelRectangularPatternAlongY");

  for (Int_t iPlane=0; iPlane<nPlanes; iPlane++) geomMFT -> Fill(zCenter[iPlane],
								 rMin[iPlane],
								 rMax[iPlane],
								 pixelSizeX[iPlane],
								 pixelSizeY[iPlane],
								 thicknessActive[iPlane],
								 thicknessSupport[iPlane],
								 thicknessReadout[iPlane],
								 equivalentSilicon[iPlane],
								 equivalentSiliconBeforeFront[iPlane],
								 equivalentSiliconBeforeBack[iPlane],
								 hasPixelRectangularPatternAlongY[iPlane]);

  TFile *fileGeomMFT = new TFile("AliMFTGeometry.root", "recreate");
  geomMFT -> Write();
  fileGeomMFT -> Close();

}

//==========================================================================================================================================
Commit	Line	Data
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
	9	void 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
190a8304	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
190a8304	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],
bcaf50eb	44	hasPixelRectangularPatternAlongY[iPlane]);
820b4d9e	45
	46	TFile *fileGeomMFT = new TFile("AliMFTGeometry.root", "recreate");
	47	geomMFT -> Write();
	48	fileGeomMFT -> Close();
	49
	50	}
	51
	52	//==========================================================================================================================================