Support files added
[u/mrichter/AliRoot.git] / MFT / AliMFTSegmentation.cxx
1 /**************************************************************************
2  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3  *                                                                        *
4  * Author: The ALICE Off-line Project.                                    *
5  * Contributors are mentioned in the code where appropriate.              *
6  *                                                                        *
7  * Permission to use, copy, modify and distribute this software and its   *
8  * documentation strictly for non-commercial purposes is hereby granted   *
9  * without fee, provided that the above copyright notice appears in all   *
10  * copies and that both the copyright notice and this permission notice   *
11  * appear in the supporting documentation. The authors make no claims     *
12  * about the suitability of this software for any purpose. It is          *
13  * provided "as is" without express or implied warranty.                  *
14  **************************************************************************/
15
16 //====================================================================================================================================================
17 //
18 //      Segmentation class for the planes of the ALICE Muon Forward Tracker
19 //
20 //      Contact author: antonio.uras@cern.ch
21 //
22 //====================================================================================================================================================
23
24 #include "TFile.h"
25 #include "TNtuple.h"
26 #include "TClonesArray.h"
27 #include "TMath.h"
28 #include "AliMFTPlane.h"
29 #include "AliMFTSegmentation.h"
30
31 ClassImp(AliMFTSegmentation)
32
33 //====================================================================================================================================================
34
35 AliMFTSegmentation::AliMFTSegmentation(): 
36   TObject(),
37   fMFTPlanes(0)
38
39
40   // default constructor
41
42 }
43
44 //====================================================================================================================================================
45
46 AliMFTSegmentation::AliMFTSegmentation(const Char_t *nameGeomFile): 
47   TObject(),
48   fMFTPlanes(new TClonesArray("AliMFTPlane", fNMaxPlanes))
49
50
51   // constructor
52
53   Float_t zCenter, rMin, rMax, pixelSizeX, pixelSizeY, thicknessActive, thicknessSupport, thicknessReadout;
54   Float_t equivalentSilicon, equivalentSiliconBeforeFront, equivalentSiliconBeforeBack;
55
56   TFile *geomFile = new TFile(nameGeomFile);
57   TNtuple *geomNtuple = (TNtuple*) geomFile->Get("AliMFTGeometry");
58
59   geomNtuple -> SetBranchAddress("zCenter", &zCenter);
60   geomNtuple -> SetBranchAddress("rMin",    &rMin);
61   geomNtuple -> SetBranchAddress("rMax",    &rMax);
62   geomNtuple -> SetBranchAddress("pixelSizeX", &pixelSizeX);
63   geomNtuple -> SetBranchAddress("pixelSizeY", &pixelSizeY);
64   geomNtuple -> SetBranchAddress("thicknessActive",  &thicknessActive);
65   geomNtuple -> SetBranchAddress("thicknessSupport", &thicknessSupport);
66   geomNtuple -> SetBranchAddress("thicknessReadout", &thicknessReadout);
67   geomNtuple -> SetBranchAddress("equivalentSilicon",            &equivalentSilicon);
68   geomNtuple -> SetBranchAddress("equivalentSiliconBeforeFront", &equivalentSiliconBeforeFront);
69   geomNtuple -> SetBranchAddress("equivalentSiliconBeforeBack",  &equivalentSiliconBeforeBack);
70   
71   Int_t nPlanes = geomNtuple->GetEntries();
72
73   for (Int_t iPlane=0; iPlane<nPlanes; iPlane++) {
74
75     // Create new plane
76
77     printf("Setting segmentation for MFT plane #%02d\n", iPlane);
78
79     geomNtuple -> GetEntry(iPlane);
80     zCenter = TMath::Abs(zCenter);
81
82     AliMFTPlane *plane = new AliMFTPlane(Form("MFTPlane_%02d", iPlane), Form("MFTPlane_%02d", iPlane));
83     plane -> Init(iPlane, zCenter, rMin, rMax, pixelSizeX, pixelSizeY, thicknessActive, thicknessSupport, thicknessReadout);
84     plane -> SetEquivalentSilicon(equivalentSilicon);
85     plane -> SetEquivalentSiliconBeforeFront(equivalentSiliconBeforeFront);
86     plane -> SetEquivalentSiliconBeforeBack(equivalentSiliconBeforeBack);
87     plane -> CreateStructure();
88     
89     new ((*fMFTPlanes)[fMFTPlanes->GetEntries()]) AliMFTPlane(*plane);
90
91   }
92   
93   delete geomFile;
94
95   printf("MFT segmentation set!\n");
96
97 }
98
99 //====================================================================================================================================================
100
101 THnSparseC* AliMFTSegmentation::GetDetElem(Int_t detElemID) const {
102       
103   // Find det elem
104
105   Int_t planeNb = detElemID/fNMaxDetElemPerPlane;
106   Int_t detElemNb = detElemID - planeNb*fNMaxDetElemPerPlane;
107   
108   THnSparseC *detElem = GetPlane(planeNb)->GetActiveElement(detElemNb);
109
110   return detElem;
111
112 }
113
114 //====================================================================================================================================================
115
116 Bool_t AliMFTSegmentation::Hit2PixelID(Double_t xHit, Double_t yHit, Int_t detElemID, Int_t &xPixel, Int_t &yPixel) {
117
118   THnSparseC *detElem = GetDetElem(detElemID);
119
120   if ( xHit<detElem->GetAxis(0)->GetXmin() ||
121        xHit>detElem->GetAxis(0)->GetXmax() ||
122        yHit<detElem->GetAxis(1)->GetXmin() ||
123        yHit>detElem->GetAxis(1)->GetXmax() ) return kFALSE;
124
125   xPixel = detElem->GetAxis(0)->FindBin(xHit) - 1;
126   yPixel = detElem->GetAxis(1)->FindBin(yHit) - 1;
127
128   return kTRUE;
129
130 }
131
132 //====================================================================================================================================================
133