]>
Commit | Line | Data |
---|---|---|
1a1fdef7 | 1 | #ifndef ALIFMDGEOMETRY_H |
2 | #define ALIFMDGEOMETRY_H | |
088f8e79 | 3 | /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights |
4 | * reserved. | |
5 | * | |
6 | * Latest changes by Christian Holm Christensen <cholm@nbi.dk> | |
7 | * | |
8 | * See cxx source for full Copyright notice | |
9 | */ | |
c2fc1258 | 10 | /** @file AliFMDGeometry.h |
11 | @author Christian Holm Christensen <cholm@nbi.dk> | |
12 | @date Mon Mar 27 12:40:37 2006 | |
13 | @brief Geometry mananger for the FMD | |
14 | */ | |
02a27b50 | 15 | //____________________________________________________________________ |
16 | // | |
69b696b9 | 17 | // Forward Multiplicity Detector based on Silicon wafers. |
18 | // | |
19 | // This class is a singleton that handles the geometry parameters of | |
20 | // the FMD detectors. | |
6169f936 | 21 | // The actual code is done by various separate classes. |
69b696b9 | 22 | // |
1a1fdef7 | 23 | #ifndef ALIGEOMETRY_H |
24 | # include <AliGeometry.h> | |
25 | #endif | |
57c3c593 | 26 | #ifndef ROOT_TArrayI |
27 | # include <TArrayI.h> | |
28 | #endif | |
29 | #ifndef ROOT_TMatrixFfwd | |
30 | # include <TMatrixFfwd.h> | |
31 | #endif | |
32 | class TVector3; | |
1a1fdef7 | 33 | class TParticle; |
34 | class AliRecPoint; | |
35 | class AliFMDRing; | |
36 | class AliFMDDetector; | |
37 | class AliFMD1; | |
38 | class AliFMD2; | |
39 | class AliFMD3; | |
54e415a8 | 40 | class AliFMDGeometryBuilder; |
1a1fdef7 | 41 | |
42 | ||
43 | //__________________________________________________________________ | |
c2fc1258 | 44 | /** @brief Singleton object of FMD geometry descriptions and parameters. |
9f662337 | 45 | This class is a singleton that handles the geometry parameters of |
46 | the FMD detectors. | |
47 | ||
48 | The actual code is done by various separate classes. Below is | |
49 | diagram showing the relationship between the various FMD classes | |
50 | that handles the geometry | |
51 | @verbatim | |
52 | +------------+ | |
53 | +- | AliFMDRing | | |
54 | 2 | +------------+ | |
55 | +----------------+<>--+ | | |
56 | | AliFMDGeometry | ^ | |
57 | +----------------+<>--+ V 1..2 | |
58 | 3 | +----------------+ | |
59 | +-| AliFMDDetector | | |
60 | +----------------+ | |
61 | ^ | |
62 | | | |
63 | +-------------+-------------+ | |
64 | | | | | |
65 | +---------+ +---------+ +---------+ | |
66 | | AliFMD1 | | AliFMD2 | | AliFMD3 | | |
67 | +---------+ +---------+ +---------+ | |
68 | ||
69 | @endverbatim | |
70 | - AliFMDRing | |
71 | This class contains all stuff needed to do with a ring. It's | |
72 | used by the AliFMDDetector objects to instantise inner and | |
73 | outer rings. The AliFMDRing objects are shared by the | |
74 | AliFMDDetector objects, and owned by the AliFMDv1 object. | |
75 | - AliFMD1, AliFMD2, and AliFMD3 | |
76 | These are specialisation of AliFMDDetector, that contains the | |
77 | particularities of each of the sub-detector system. It is | |
78 | envisioned that the classes should also define the support | |
79 | volumes and material for each of the detectors. | |
80 | ||
81 | @ingroup FMD_base | |
1a1fdef7 | 82 | */ |
83 | class AliFMDGeometry : public AliGeometry | |
84 | { | |
85 | public: | |
09b6c804 | 86 | /** |
87 | * singleton access | |
88 | * | |
89 | * @return Singleton | |
90 | */ | |
1a1fdef7 | 91 | static AliFMDGeometry* Instance(); |
09b6c804 | 92 | /** |
93 | * Initialize the the singleton if not done so already | |
94 | */ | |
1a1fdef7 | 95 | virtual void Init(); |
09b6c804 | 96 | /** |
97 | * Find all local <-> global transforms | |
98 | */ | |
352482dc | 99 | virtual void InitTransformations(Bool_t force=kFALSE); |
09b6c804 | 100 | /** |
101 | * @return Get inner description | |
102 | */ | |
1a1fdef7 | 103 | AliFMDRing* GetInner() const { return fInner; } |
09b6c804 | 104 | /** |
105 | * @return Get outer description | |
106 | */ | |
1a1fdef7 | 107 | AliFMDRing* GetOuter() const { return fOuter; } |
09b6c804 | 108 | /** |
109 | * @return Get FMD1 description | |
110 | */ | |
1a1fdef7 | 111 | AliFMD1* GetFMD1() const { return (fUseFMD1 ? fFMD1 : 0); } |
09b6c804 | 112 | /** |
113 | * @return Get FMD2 description | |
114 | */ | |
1a1fdef7 | 115 | AliFMD2* GetFMD2() const { return (fUseFMD2 ? fFMD2 : 0); } |
09b6c804 | 116 | /** |
117 | * @return Get FMD3 description | |
118 | */ | |
1a1fdef7 | 119 | AliFMD3* GetFMD3() const { return (fUseFMD3 ? fFMD3 : 0); } |
09b6c804 | 120 | /** |
121 | * Get description of a sub-detector | |
122 | * | |
123 | * @param i Sub-detector # | |
124 | * @return Description of sub-detector, or 0 | |
125 | */ | |
1a1fdef7 | 126 | AliFMDDetector* GetDetector(Int_t i) const; |
09b6c804 | 127 | /** |
128 | * Get description of a ring, i should be one of 'I' or 'O' (case | |
129 | * insensitive). If an invalid parameter is passed, 0 (NULL) is | |
130 | * returned. | |
131 | * | |
132 | * @param i Ring id | |
133 | * @return Description of ring, or 0 | |
134 | */ | |
1a1fdef7 | 135 | AliFMDRing* GetRing(Char_t i) const; |
09b6c804 | 136 | /** |
137 | * Disable the ith detector | |
138 | * | |
139 | * @param i IF true, disable sub-detector @a i | |
140 | */ | |
1a1fdef7 | 141 | void Disable(Int_t i); |
09b6c804 | 142 | /** |
143 | * Enable the ith detector | |
144 | * | |
145 | * @param i IF true, enable sub-detector @a i | |
146 | */ | |
1a1fdef7 | 147 | void Enable(Int_t i); |
09b6c804 | 148 | /** |
149 | * @return Density @f$ \rho@f$ of silicon | |
150 | */ | |
1a1fdef7 | 151 | Double_t GetSiDensity() const { return 2.33; } |
09b6c804 | 152 | /** |
153 | * Translate detector coordinates (detector, ring, sector, strip) | |
154 | * to spatial coordinates (x, y, z) in the master reference frame | |
155 | * of ALICE. The member function uses the transformations | |
156 | * previously obtained from the TGeoManager. | |
157 | * | |
158 | * @param detector Detector number | |
159 | * @param ring Ring id | |
160 | * @param sector Sector number | |
161 | * @param strip Strip number | |
162 | * @param x On return, X coordinate | |
163 | * @param y On return, Y coordinate | |
164 | * @param z On return, Z coordinate | |
165 | */ | |
54e415a8 | 166 | void Detector2XYZ(UShort_t detector, Char_t ring, |
167 | UShort_t sector, UShort_t strip, | |
168 | Double_t& x, Double_t& y, Double_t& z) const; | |
09b6c804 | 169 | /** |
170 | * Translate spatial coordinates (x,y,z) in the master reference | |
171 | * frame of ALICE to the detector coordinates (detector, ring, | |
172 | * sector, strip). Note, that if this method is to be used in | |
173 | * reconstruction or the like, then the input z-coordinate should | |
174 | * be corrected for the events interactions points z-coordinate, | |
175 | * like | |
176 | * @code | |
177 | * geom->XYZ2Detector(x,y,z-ipz,d,r,s,t); | |
178 | * @endcode | |
179 | * | |
180 | * @param x X coordinate | |
181 | * @param y Y coordinate | |
182 | * @param z Z coordinate | |
183 | * @param detector On return, Detector number | |
184 | * @param ring On return, Ring id | |
185 | * @param sector On return, Sector number | |
186 | * @param strip On return, Strip number | |
187 | * @return @c false of (@a x, @a y, @a z) is not within this | |
188 | * detector. | |
189 | */ | |
54e415a8 | 190 | Bool_t XYZ2Detector(Double_t x, Double_t y, Double_t z, |
191 | UShort_t& detector, Char_t& ring, | |
192 | UShort_t& sector, UShort_t& strip) const; | |
09b6c804 | 193 | /** |
194 | * Make the geometry. This delegates to AliFMDGeometryBuilder | |
195 | */ | |
54e415a8 | 196 | void Build(); |
09b6c804 | 197 | /** |
198 | * @return Get detector offset in paths | |
199 | */ | |
54e415a8 | 200 | Int_t GetDetectorOff() const { return fDetectorOff; } |
09b6c804 | 201 | /** |
202 | * @return Get sensor offset in paths | |
203 | */ | |
54e415a8 | 204 | Int_t GetModuleOff() const { return fModuleOff; } |
09b6c804 | 205 | /** |
206 | * @return Get ring offset in paths | |
207 | */ | |
54e415a8 | 208 | Int_t GetRingOff() const { return fRingOff; } |
09b6c804 | 209 | /** |
210 | * @return Get ring sector in paths | |
211 | */ | |
54e415a8 | 212 | Int_t GetSectorOff() const { return fSectorOff; } |
09b6c804 | 213 | /** |
214 | * @param off Detector off-set set in geometry path | |
215 | */ | |
54e415a8 | 216 | void SetDetectorOff(Int_t off) { fDetectorOff = off; } |
09b6c804 | 217 | /** |
218 | * @param off Module off-set set in geometry path | |
219 | */ | |
54e415a8 | 220 | void SetModuleOff(Int_t off) { fModuleOff = off; } |
09b6c804 | 221 | /** |
222 | * @param off Ring off-set set in geometry path | |
223 | */ | |
54e415a8 | 224 | void SetRingOff(Int_t off) { fRingOff = off; } |
09b6c804 | 225 | /** |
226 | * @param off Sectord off-set set in geometry path | |
227 | */ | |
54e415a8 | 228 | void SetSectorOff(Int_t off) { fSectorOff = off; } |
09b6c804 | 229 | /** |
230 | * Check if volume @a vol is marked as active | |
231 | * | |
232 | * @param vol Volume ID | |
233 | * @return @c true if @a vol is declared active | |
234 | */ | |
54e415a8 | 235 | Bool_t IsActive(Int_t vol) const; |
09b6c804 | 236 | /** |
237 | * Set active volumes | |
238 | * | |
239 | * @param active Active volume id array | |
240 | * @param n elements of @a active | |
241 | */ | |
54e415a8 | 242 | void SetActive(Int_t* active, Int_t n); |
09b6c804 | 243 | /** |
244 | * Add an active volume | |
245 | * | |
246 | * @param id Register volume @a id to be active | |
247 | */ | |
54e415a8 | 248 | void AddActive(Int_t id); |
09b6c804 | 249 | /** |
250 | * Get Array of active volume numbers | |
251 | * | |
252 | * @return constant reference to active volume numbers | |
253 | */ | |
45855cd7 | 254 | const TArrayI& ActiveIds() const { return fActive; } |
09b6c804 | 255 | /** |
256 | * Set an external geometry builder | |
257 | * | |
258 | * @param b Geometry builder | |
259 | */ | |
54e415a8 | 260 | void SetBuilder(AliFMDGeometryBuilder* b) { fBuilder = b; } |
09b6c804 | 261 | /** |
262 | * Check the volume depth of some nodes, get the active volume | |
263 | * numbers, and so forth. | |
264 | * | |
265 | * @todo Here, we should actually also get the parameters of the | |
266 | * shapes, like the verticies of the polygon shape that makes up the | |
267 | * silicon sensor, the strip pitch, the ring radii, the z-positions, | |
268 | * and so on - that is, all the geometric information we need for | |
269 | * futher processing, such as simulation, digitization, | |
270 | * reconstruction, etc. | |
271 | */ | |
54e415a8 | 272 | void ExtractGeomInfo(); |
09b6c804 | 273 | /** |
274 | * Whether we are to use a detailed geometry or not | |
275 | * | |
276 | * @param det if @c true, make a detailed geometry. | |
277 | */ | |
54e415a8 | 278 | void SetDetailed(Bool_t det) { fDetailed = det; } |
09b6c804 | 279 | /** |
280 | * @return @c true if geometry is detailed | |
281 | */ | |
54e415a8 | 282 | Bool_t IsDetailed() const { return fDetailed; } |
09b6c804 | 283 | /** |
284 | * @param ass Whether to use assemblies or not | |
285 | */ | |
54e415a8 | 286 | void UseAssembly(Bool_t ass) { fUseAssembly = ass; } |
1a1fdef7 | 287 | |
288 | // AliGeometry member functions | |
09b6c804 | 289 | /** |
290 | * Get global coordinates cooresponding to a rec point. | |
291 | * | |
292 | * @param p Reconstructed point. | |
293 | * @param pos On return, the position | |
294 | * @param mat On return, the material at @a post | |
295 | */ | |
1a1fdef7 | 296 | virtual void GetGlobal(const AliRecPoint* p, TVector3& pos, |
e8d02863 | 297 | TMatrixF& mat) const; |
09b6c804 | 298 | /** |
299 | * Get global coordinates cooresponding to a rec point. | |
300 | * | |
301 | * @param p Reconstructed point. | |
302 | * @param pos On return, the position | |
303 | */ | |
1a1fdef7 | 304 | virtual void GetGlobal(const AliRecPoint* p, TVector3& pos) const; |
09b6c804 | 305 | /** |
306 | * Check if particle will hit an active detector element. | |
307 | * | |
308 | * @todo implement this function | |
309 | * | |
310 | * @param particle Track | |
311 | * @return @c true if @a particle will hit this detector | |
312 | */ | |
1a1fdef7 | 313 | virtual Bool_t Impact(const TParticle* particle) const; |
09b6c804 | 314 | /** |
315 | * Declare alignable volumes | |
316 | */ | |
9de78b35 | 317 | virtual void SetAlignableVolumes() const; |
75609cab | 318 | |
319 | ||
320 | /** | |
321 | * Service function to convert Cartisean XYZ to r, eta, phi, and theta. | |
322 | * | |
323 | * Note, that the z input should be corrected for the vertex location | |
324 | * if needed. | |
325 | * | |
326 | * @param x Cartisean X coordinate | |
327 | * @param y Cartisean Y coordinate | |
328 | * @param z Cartisean Z coordinate | |
329 | * @param r On return, the radius | |
330 | * @param eta On return, the pseudo-rapidity | |
331 | * @param phi On return, the azimuthal angle | |
332 | * @param theta On return, the polar angle; | |
333 | * | |
334 | * @return kTRUE on success, kFALSE in case of problems | |
335 | */ | |
336 | static Bool_t XYZ2REtaPhiTheta(Double_t x, Double_t y, | |
337 | Double_t z, | |
338 | Double_t& r, Double_t& eta, | |
339 | Double_t& phi, Double_t& theta); | |
340 | ||
1a1fdef7 | 341 | protected: |
69b696b9 | 342 | Bool_t fIsInitialized; // Whether singleton is initalized |
1a1fdef7 | 343 | AliFMDRing* fInner; // Inner ring geometry information |
344 | AliFMDRing* fOuter; // Outer ring geometry information | |
345 | AliFMD1* fFMD1; // FMD1 geometry information | |
346 | AliFMD2* fFMD2; // FMD2 geometry information | |
347 | AliFMD3* fFMD3; // FMD3 geometry information | |
348 | Bool_t fUseFMD1; // Wheter to Use FMD1 or not | |
349 | Bool_t fUseFMD2; // Wheter to Use FMD2 or not | |
350 | Bool_t fUseFMD3; // Wheter to Use FMD3 or not | |
352482dc | 351 | Bool_t fIsInitTrans; // Transforms initialised? |
69b696b9 | 352 | static AliFMDGeometry* fgInstance; // Singleton instance |
09b6c804 | 353 | /** |
354 | * CTOR | |
355 | */ | |
1a1fdef7 | 356 | AliFMDGeometry(); |
09b6c804 | 357 | /** |
358 | * CTOR | |
359 | * | |
360 | * @param name Not used | |
361 | */ | |
362 | AliFMDGeometry(const char* name); | |
363 | /** | |
364 | * Copy CTOR | |
365 | * | |
366 | * @param other To copy from | |
367 | */ | |
69b696b9 | 368 | AliFMDGeometry(const AliFMDGeometry& other); |
09b6c804 | 369 | /** |
370 | * Assignment operator | |
371 | * | |
372 | * @param other To assig from | |
373 | * @return reference to this. | |
374 | */ | |
69b696b9 | 375 | AliFMDGeometry& operator=(const AliFMDGeometry& other); |
09b6c804 | 376 | /** |
377 | * Destructor | |
378 | */ | |
1a1fdef7 | 379 | virtual ~AliFMDGeometry() {} |
69b696b9 | 380 | |
02a27b50 | 381 | AliFMDGeometryBuilder* fBuilder; // Geometry builder |
382 | Int_t fDetectorOff; // Detector off-set | |
383 | Int_t fModuleOff; // Module off-set | |
384 | Int_t fRingOff; // ring offset | |
385 | Int_t fSectorOff; // Sector offset | |
386 | TArrayI fActive; // Active volumes | |
387 | Bool_t fDetailed; // Whether to make detailed geom | |
388 | Bool_t fUseAssembly; // Whther to use assemblies | |
54e415a8 | 389 | |
02a27b50 | 390 | ClassDef(AliFMDGeometry,1); // Geometry parameters and manager |
1a1fdef7 | 391 | }; |
392 | ||
393 | ||
394 | #endif | |
395 | //____________________________________________________________________ | |
396 | // | |
397 | // Local Variables: | |
398 | // mode: C++ | |
399 | // End: | |
400 | // | |
401 | // EOF | |
402 | // |