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"
17 class AliTRDparameter;
19 class AliTRDgeometry : public AliGeometry {
23 enum { kNplan = 6, kNcham = 5, kNsect = 18, kNdet = 540, kNdets = 30 };
26 virtual ~AliTRDgeometry();
28 virtual void CreateGeometry(Int_t *idtmed);
29 virtual Int_t IsVersion() const = 0;
31 virtual Bool_t Impact(const TParticle* ) const { return kTRUE; };
32 virtual Bool_t Local2Global(Int_t d, Double_t *local, Double_t *global) const;
33 virtual Bool_t Local2Global(Int_t p, Int_t c, Int_t s
34 , Double_t *local, Double_t *global) const;
36 virtual Bool_t Global2Local(Int_t mode, Double_t *local, Double_t *global
37 , Int_t* index, AliTRDparameter *par) const;
38 virtual Bool_t Global2Detector(Double_t global[3], Int_t index[3]);
40 virtual Bool_t Rotate(Int_t d, Double_t *pos, Double_t *rot) const;
41 virtual Bool_t RotateBack(Int_t d, Double_t *rot, Double_t *pos) const;
43 static Int_t Nsect() { return fgkNsect; };
44 static Int_t Nplan() { return fgkNplan; };
45 static Int_t Ncham() { return fgkNcham; };
46 static Int_t Ndet() { return fgkNdet; };
48 static Float_t Rmin() { return fgkRmin; };
49 static Float_t Rmax() { return fgkRmax; };
50 static Float_t Zmax1() { return fgkZmax1; };
51 static Float_t Zmax2() { return fgkZmax2; };
53 static Float_t Cwidcha() { return (fgkSwidth2 - fgkSwidth1)
54 / fgkSheight * (fgkCH + fgkVspace); };
55 static Float_t Cheight() { return fgkCH; };
56 static Float_t Cspace() { return fgkVspace; };
57 static Float_t CraHght() { return fgkCraH; };
58 static Float_t CdrHght() { return fgkCdrH; };
59 static Float_t CamHght() { return fgkCamH; };
60 static Float_t CroHght() { return fgkCroH; };
61 static Float_t CroWid() { return fgkCroW; };
62 static Float_t MyThick() { return fgkMyThick; };
63 static Float_t DrThick() { return fgkDrThick; };
64 static Float_t AmThick() { return fgkAmThick; };
65 static Float_t DrZpos() { return fgkDrZpos; };
66 static Float_t RpadW() { return fgkRpadW; };
67 static Float_t CpadW() { return fgkCpadW; };
69 virtual void SetPHOShole() = 0;
70 virtual void SetRICHhole() = 0;
72 virtual Bool_t GetPHOShole() const = 0;
73 virtual Bool_t GetRICHhole() const = 0;
74 virtual Bool_t IsHole(Int_t /*iplan*/, Int_t /*icham*/, Int_t /*isect*/) const {return kFALSE;}
75 static Int_t GetDetectorSec(Int_t p, Int_t c);
76 static Int_t GetDetector(Int_t p, Int_t c, Int_t s);
77 virtual Int_t GetPlane(Int_t d) const;
78 virtual Int_t GetChamber(Int_t d) const;
79 virtual Int_t GetSector(Int_t d) const;
81 Float_t GetChamberWidth(Int_t p) const { return fCwidth[p]; };
82 Float_t GetChamberLength(Int_t p, Int_t c) const { return fClength[p][c]; };
84 virtual void GetGlobal(const AliRecPoint* , TVector3& , TMatrixF& ) const { };
85 virtual void GetGlobal(const AliRecPoint* , TVector3& ) const { };
87 static Double_t GetAlpha() { return 2 * 3.14159265358979323846 / fgkNsect; };
89 static AliTRDgeometry* GetGeometry(AliRunLoader* runLoader = NULL);
91 static Float_t GetTime0(Int_t p) { return fgkTime0[p]; };
95 static const Int_t fgkNsect; // Number of sectors in the full detector (18)
96 static const Int_t fgkNplan; // Number of planes of the TRD (6)
97 static const Int_t fgkNcham; // Number of chambers in z-direction (5)
98 static const Int_t fgkNdet; // Total number of detectors (18 * 6 * 5 = 540)
100 static const Float_t fgkRmin; // Minimal radius of the TRD
101 static const Float_t fgkRmax; // Maximal radius of the TRD
103 static const Float_t fgkZmax1; // Half-length of the TRD at outer radius
104 static const Float_t fgkZmax2; // Half-length of the TRD at inner radius
106 static const Float_t fgkSheight; // Height of the TRD-volume in spaceframe (BTR1-3)
107 static const Float_t fgkSwidth1; // Lower width of the TRD-volume in spaceframe (BTR1-3)
108 static const Float_t fgkSwidth2; // Upper width of the TRD-volume in spaceframe (BTR1-3)
109 static const Float_t fgkSlenTR1; // Length of the TRD-volume in spaceframe (BTR1)
110 static const Float_t fgkSlenTR2; // Length of the TRD-volume in spaceframe (BTR2)
111 static const Float_t fgkSlenTR3; // Length of the TRD-volume in spaceframe (BTR3)
113 static const Float_t fgkSMpltT; // Thickness of the super module side plates
114 static const Float_t fgkSMgapT; // Thickness of the gap between side plates and space frame
116 static const Float_t fgkCraH; // Height of the radiator part of the chambers
117 static const Float_t fgkCdrH; // Height of the drift region of the chambers
118 static const Float_t fgkCamH; // Height of the amplification region of the chambers
119 static const Float_t fgkCroH; // Height of the readout of the chambers
120 static const Float_t fgkCH; // Total height of the chambers
122 static const Float_t fgkVspace; // Vertical spacing of the chambers
123 static const Float_t fgkHspace; // Horizontal spacing of the chambers
125 static const Float_t fgkCalT; // Thickness of the lower aluminum frame
126 static const Float_t fgkCclsT; // Thickness of the lower G10 frame sides
127 static const Float_t fgkCclfT; // Thickness of the lower G10 frame front
128 static const Float_t fgkCcuT; // Thickness of the upper G10 frame
129 static const Float_t fgkCauT; // Thickness of the upper aluminum frame
131 static const Float_t fgkCroW; // Additional width of the readout chamber frames
133 static const Float_t fgkCpadW; // Difference of outer chamber width and pad plane width
134 static const Float_t fgkRpadW; // Difference of outer chamber width and pad plane width
136 static const Float_t fgkRaThick; // Thickness of the radiator
137 static const Float_t fgkMyThick; // Thickness of the mylar-layer
138 static const Float_t fgkXeThick; // Thickness of the gas volume
139 static const Float_t fgkDrThick; // Thickness of the drift region
140 static const Float_t fgkAmThick; // Thickness of the amplification region
141 static const Float_t fgkCuThick; // Thickness of the pad plane
142 static const Float_t fgkSuThick; // Thickness of the HEXCEL+G10 support structure
143 static const Float_t fgkFeThick; // Thickness of the FEE + signal lines
144 static const Float_t fgkCoThick; // Thickness of the PE of the cooling device
145 static const Float_t fgkWaThick; // Thickness of the cooling water
147 static const Float_t fgkRaZpos; // Position of the radiator
148 static const Float_t fgkMyZpos; // Position of the mylar-layer
149 static const Float_t fgkDrZpos; // Position of the drift region
150 static const Float_t fgkAmZpos; // Position of the amplification region
151 static const Float_t fgkCuZpos; // Position of the pad plane
152 static const Float_t fgkSuZpos; // Position of the HEXCEL+G10 support structure
153 static const Float_t fgkFeZpos; // Position of the FEE + signal lines
154 static const Float_t fgkCoZpos; // Position of the PE of the cooling device
155 static const Float_t fgkWaZpos; // Position of the colling water
157 Float_t fCwidth[kNplan]; // Outer widths of the chambers
158 Float_t fClength[kNplan][kNcham]; // Outer lengths of the chambers
159 Float_t fClengthPH[kNplan][kNcham]; // For sectors with holes for the PHOS
160 Float_t fClengthRH[kNplan][kNcham]; // For sectors with holes for the RICH
162 Float_t fRotA11[kNsect]; // Matrix elements for the rotation
163 Float_t fRotA12[kNsect]; // Matrix elements for the rotation
164 Float_t fRotA21[kNsect]; // Matrix elements for the rotation
165 Float_t fRotA22[kNsect]; // Matrix elements for the rotation
167 Float_t fRotB11[kNsect]; // Matrix elements for the backward rotation
168 Float_t fRotB12[kNsect]; // Matrix elements for the backward rotation
169 Float_t fRotB21[kNsect]; // Matrix elements for the backward rotation
170 Float_t fRotB22[kNsect]; // Matrix elements for the backward rotation
172 static const Double_t fgkTime0Base; // Base value for calculation of Time-position of pad 0
173 static const Float_t fgkTime0[kNplan]; // Time-position of pad 0
175 Float_t fChamberUAorig[3*kNdets][3]; // Volumes origin in
176 Float_t fChamberUDorig[3*kNdets][3]; // the chamber
177 Float_t fChamberUForig[3*kNdets][3]; // [3] = x, y, z
178 Float_t fChamberUUorig[3*kNdets][3];
180 Float_t fChamberUAboxd[3*kNdets][3]; // Volumes box
181 Float_t fChamberUDboxd[3*kNdets][3]; // dimensions (half)
182 Float_t fChamberUFboxd[3*kNdets][3]; // [3] = x, y, z
183 Float_t fChamberUUboxd[3*kNdets][3];
185 ClassDef(AliTRDgeometry,7) // TRD geometry base class