1 #ifndef ALITRDGEOMETRY_H
2 #define ALITRDGEOMETRY_H
3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
8 ///////////////////////////////////////////////////////////////////////////////
10 // TRD geometry class //
12 ///////////////////////////////////////////////////////////////////////////////
14 #include "AliGeometry.h"
16 #include "TObjArray.h"
21 class AliTRDgeometry : public AliGeometry {
32 AliTRDgeometry(const AliTRDgeometry &g);
33 virtual ~AliTRDgeometry();
34 AliTRDgeometry &operator=(const AliTRDgeometry &g);
37 virtual void CreateGeometry(Int_t *idtmed);
38 virtual Int_t IsVersion() { return 1; }
39 virtual Bool_t Impact(const TParticle* ) const { return kTRUE; }
40 virtual Bool_t IsHole(Int_t /*p*/, Int_t /*c*/, Int_t /*s*/) const { return kFALSE; }
42 virtual Bool_t Local2Global(Int_t d, Double_t *local, Double_t *global) const;
43 virtual Bool_t Local2Global(Int_t p, Int_t c, Int_t s
44 , Double_t *local, Double_t *global) const;
45 virtual Bool_t Global2Local(Int_t mode, Double_t *local, Double_t *global
46 , Int_t *index) const;
47 virtual Bool_t Global2Detector(Double_t global[3], Int_t index[3]);
48 virtual Bool_t Rotate(Int_t d, Double_t *pos, Double_t *rot) const;
49 virtual Bool_t RotateBack(Int_t d, Double_t *rot, Double_t *pos) const;
51 void GroupChamber(Int_t iplan, Int_t icham, Int_t *idtmed);
52 void CreateFrame(Int_t *idtmed);
53 void CreateServices(Int_t *idtmed);
55 Bool_t ReadGeoMatrices();
57 void SetSMstatus(Int_t sm, Char_t status) { fSMstatus[sm] = status; }
59 static AliTRDgeometry* GetGeometry(AliRunLoader *runLoader = NULL);
61 static Int_t GetDetectorSec(Int_t p, Int_t c);
62 static Int_t GetDetector(Int_t p, Int_t c, Int_t s);
63 virtual Int_t GetPlane(Int_t d) const;
64 virtual Int_t GetChamber(Int_t d) const;
65 virtual Int_t GetSector(Int_t d) const;
67 static Float_t GetTime0(Int_t p) { return fgkTime0[p]; }
69 Char_t GetSMstatus(Int_t sm) const { return fSMstatus[sm]; }
70 Float_t GetChamberWidth(Int_t p) const { return fCwidth[p]; }
71 Float_t GetChamberLength(Int_t p, Int_t c) const { return fClength[p][c]; }
73 virtual void GetGlobal(const AliRecPoint*, TVector3&, TMatrixF& ) const { };
74 virtual void GetGlobal(const AliRecPoint*, TVector3& ) const { };
76 static Double_t GetAlpha() { return 2.0
80 static Int_t Nsect() { return fgkNsect; }
81 static Int_t Nplan() { return fgkNplan; }
82 static Int_t Ncham() { return fgkNcham; }
83 static Int_t Ndet() { return fgkNdet; }
85 static Float_t Cheight() { return fgkCH; }
86 static Float_t Cspace() { return fgkVspace; }
87 static Float_t CraHght() { return fgkCraH; }
88 static Float_t CdrHght() { return fgkCdrH; }
89 static Float_t CamHght() { return fgkCamH; }
90 static Float_t CroHght() { return fgkCroH; }
91 static Float_t CroWid() { return fgkCroW; }
92 static Float_t MyThick() { return fgkMyThick; }
93 static Float_t DrThick() { return fgkDrThick; }
94 static Float_t AmThick() { return fgkAmThick; }
95 static Float_t DrZpos() { return fgkDrZpos; }
96 static Float_t RpadW() { return fgkRpadW; }
97 static Float_t CpadW() { return fgkCpadW; }
99 static Float_t Cwidcha() { return (fgkSwidth2 - fgkSwidth1)
101 * (fgkCH + fgkVspace); }
103 TGeoHMatrix *GetGeoMatrix(Int_t det) { return (TGeoHMatrix *)
104 fMatrixGeo->At(det); }
105 TGeoHMatrix *GetMatrix(Int_t det) { return (TGeoHMatrix *)
106 fMatrixArray->At(det); }
107 TGeoHMatrix *GetCorrectionMatrix(Int_t det) { return (TGeoHMatrix *)
108 fMatrixCorrectionArray->At(det); }
112 static const Int_t fgkNsect; // Number of sectors in the full detector (18)
113 static const Int_t fgkNplan; // Number of planes of the TRD (6)
114 static const Int_t fgkNcham; // Number of chambers in z-direction (5)
115 static const Int_t fgkNdet; // Total number of detectors (18 * 6 * 5 = 540)
117 static const Float_t fgkSheight; // Height of the TRD-volume in spaceframe (BTRD)
118 static const Float_t fgkSwidth1; // Lower width of the TRD-volume in spaceframe (BTRD)
119 static const Float_t fgkSwidth2; // Upper width of the TRD-volume in spaceframe (BTRD)
120 static const Float_t fgkSlength; // Length of the TRD-volume in spaceframe (BTRD)
122 static const Float_t fgkSMpltT; // Thickness of the super module side plates
124 static const Float_t fgkCraH; // Height of the radiator part of the chambers
125 static const Float_t fgkCdrH; // Height of the drift region of the chambers
126 static const Float_t fgkCamH; // Height of the amplification region of the chambers
127 static const Float_t fgkCroH; // Height of the readout of the chambers
128 static const Float_t fgkCH; // Total height of the chambers
130 static const Float_t fgkVspace; // Vertical spacing of the chambers
131 static const Float_t fgkHspace; // Horizontal spacing of the chambers
132 static const Float_t fgkVrocsm; // Radial distance of the first ROC to the outer SM plates
133 static const Float_t fgkCalT; // Thickness of the lower aluminum frame
134 static const Float_t fgkCalW; // Width of additional aluminum on lower frame
135 static const Float_t fgkCclsT; // Thickness of the lower Wacosit frame sides
136 static const Float_t fgkCclfT; // Thickness of the lower Wacosit frame front
137 static const Float_t fgkCglT; // Thichness of the glue around the radiator
138 static const Float_t fgkCcuT; // Thickness of the upper Wacosit frame
139 static const Float_t fgkCauT; // Thickness of the aluminum frame of the back panel
141 static const Float_t fgkCroW; // Additional width of the readout chamber frames
143 static const Float_t fgkCpadW; // Difference of outer chamber width and pad plane width
144 static const Float_t fgkRpadW; // Difference of outer chamber width and pad plane width
146 static const Float_t fgkMyThick; // Thickness of the mylar-layer
147 static const Float_t fgkRaThick; // Thickness of the radiator
148 static const Float_t fgkXeThick; // Thickness of the gas volume
149 static const Float_t fgkDrThick; // Thickness of the drift region
150 static const Float_t fgkAmThick; // Thickness of the amplification region
151 static const Float_t fgkWrThick; // Thickness of the wire planes
152 static const Float_t fgkCuThick; // Thickness of the pad plane
153 static const Float_t fgkGlThick; // Thickness of the glue layer
154 static const Float_t fgkSuThick; // Thickness of the NOMEX support structure
155 static const Float_t fgkRpThick; // Thickness of the PCB readout boards
156 static const Float_t fgkRcThick; // Thickness of the PCB copper layers
157 static const Float_t fgkRoThick; // Thickness of all other ROB componentes (caps, etc.)
159 static const Float_t fgkRaZpos; // Position of the radiator
160 static const Float_t fgkDrZpos; // Position of the drift region
161 static const Float_t fgkAmZpos; // Position of the amplification region
162 static const Float_t fgkWrZpos; // Position of the wire planes
163 static const Float_t fgkCuZpos; // Position of the pad plane
164 static const Float_t fgkGlZpos; // Position of the glue layer
165 static const Float_t fgkSuZpos; // Position of the HEXCEL+G10 support structure
166 static const Float_t fgkRpZpos; // Position of the PCB readout boards
167 static const Float_t fgkRcZpos; // Position of the PCB copper layers
168 static const Float_t fgkRoZpos; // Position of all other ROB componentes (caps, etc.)
170 Char_t fSMstatus[kNsect]; // Super module status byte
172 Float_t fCwidth[kNplan]; // Outer widths of the chambers
173 Float_t fClength[kNplan][kNcham]; // Outer lengths of the chambers
175 Float_t fRotA11[kNsect]; // Matrix elements for the rotation
176 Float_t fRotA12[kNsect]; // Matrix elements for the rotation
177 Float_t fRotA21[kNsect]; // Matrix elements for the rotation
178 Float_t fRotA22[kNsect]; // Matrix elements for the rotation
180 Float_t fRotB11[kNsect]; // Matrix elements for the backward rotation
181 Float_t fRotB12[kNsect]; // Matrix elements for the backward rotation
182 Float_t fRotB21[kNsect]; // Matrix elements for the backward rotation
183 Float_t fRotB22[kNsect]; // Matrix elements for the backward rotation
185 static const Double_t fgkTime0Base; // Base value for calculation of Time-position of pad 0
186 static const Float_t fgkTime0[kNplan]; // Time-position of pad 0
188 Float_t fChamberUAorig[3*kNdets][3]; // Volumes origin in
189 Float_t fChamberUDorig[3*kNdets][3]; // the chamber
190 Float_t fChamberUForig[3*kNdets][3]; // [3] = x, y, z
191 Float_t fChamberUUorig[3*kNdets][3]; //
193 Float_t fChamberUAboxd[3*kNdets][3]; // Volumes box
194 Float_t fChamberUDboxd[3*kNdets][3]; // dimensions (half)
195 Float_t fChamberUFboxd[3*kNdets][3]; // [3] = x, y, z
196 Float_t fChamberUUboxd[3*kNdets][3]; //
198 TObjArray * fMatrixArray; //! Transformation Global to Local
199 TObjArray * fMatrixCorrectionArray; //! Transformation Cluster to Tracking systerm
200 TObjArray * fMatrixGeo; //! Geo matrices
202 ClassDef(AliTRDgeometry,11) // TRD geometry class