]>
Commit | Line | Data |
---|---|---|
d1cd2474 | 1 | // $Id$ |
2 | // | |
3 | // Class AliMUONSt1GeometryBuilder | |
4 | // ------------------------------- | |
5 | // MUON Station1 coarse geometry construction class. | |
d1cd2474 | 6 | // Extracted from AliMUONv1 |
7 | // by Ivana Hrivnacova, IPN Orsay | |
5f1df83a | 8 | // Included in AliRoot 2004/01/23 |
d1cd2474 | 9 | |
10 | #include <TVirtualMC.h> | |
11 | #include <TGeoMatrix.h> | |
12 | ||
13 | #include "AliMUONSt1GeometryBuilder.h" | |
14 | #include "AliMUON.h" | |
15 | #include "AliMUONChamber.h" | |
16 | #include "AliMUONChamberGeometry.h" | |
a432117a | 17 | #include "AliMUONGeometryEnvelopeStore.h" |
8c343c7c | 18 | #include "AliLog.h" |
d1cd2474 | 19 | |
20 | ClassImp(AliMUONSt1GeometryBuilder) | |
21 | ||
22 | //______________________________________________________________________________ | |
23 | AliMUONSt1GeometryBuilder::AliMUONSt1GeometryBuilder(AliMUON* muon) | |
a432117a | 24 | : AliMUONVGeometryBuilder("st1.dat", |
25 | &muon->Chamber(0), &muon->Chamber(1)), | |
d1cd2474 | 26 | fMUON(muon) |
27 | { | |
28 | // Standard constructor | |
29 | ||
30 | } | |
31 | ||
32 | //______________________________________________________________________________ | |
33 | AliMUONSt1GeometryBuilder::AliMUONSt1GeometryBuilder() | |
34 | : AliMUONVGeometryBuilder(), | |
35 | fMUON(0) | |
36 | { | |
37 | // Default constructor | |
38 | } | |
39 | ||
40 | ||
41 | //______________________________________________________________________________ | |
42 | AliMUONSt1GeometryBuilder::AliMUONSt1GeometryBuilder(const AliMUONSt1GeometryBuilder& rhs) | |
43 | : AliMUONVGeometryBuilder(rhs) | |
44 | { | |
30178c30 | 45 | // Protected copy constructor |
46 | ||
8c343c7c | 47 | AliFatal("Copy constructor is not implemented."); |
d1cd2474 | 48 | } |
49 | ||
50 | //______________________________________________________________________________ | |
51 | AliMUONSt1GeometryBuilder::~AliMUONSt1GeometryBuilder() { | |
52 | // | |
53 | } | |
54 | ||
55 | //______________________________________________________________________________ | |
56 | AliMUONSt1GeometryBuilder& | |
57 | AliMUONSt1GeometryBuilder::operator = (const AliMUONSt1GeometryBuilder& rhs) | |
58 | { | |
30178c30 | 59 | // Protected assignement operator |
60 | ||
d1cd2474 | 61 | // check assignement to self |
62 | if (this == &rhs) return *this; | |
63 | ||
8c343c7c | 64 | AliFatal("Assignment operator is not implemented."); |
d1cd2474 | 65 | |
66 | return *this; | |
67 | } | |
68 | ||
69 | // | |
70 | // public methods | |
71 | // | |
72 | ||
73 | //______________________________________________________________________________ | |
74 | void AliMUONSt1GeometryBuilder::CreateGeometry() | |
75 | { | |
76 | // From AliMUONv1::CreateGeometry() | |
77 | // --- | |
78 | ||
79 | //******************************************************************** | |
80 | // Station 1 ** | |
81 | //******************************************************************** | |
82 | // CONCENTRIC | |
83 | // indices 1 and 2 for first and second chambers in the station | |
84 | // iChamber (first chamber) kept for other quanties than Z, | |
85 | // assumed to be the same in both chambers | |
86 | ||
87 | // Get tracking medias Ids | |
88 | Int_t *idtmed = fMUON->GetIdtmed()->GetArray()-1099; | |
89 | Int_t idAir= idtmed[1100]; // medium 1 | |
90 | Int_t idAlu1=idtmed[1103]; // medium 4 | |
91 | Int_t idAlu2=idtmed[1104]; // medium 5 | |
92 | Int_t idGas=idtmed[1108]; // medium 9 = Ar-CO2 gas (80%+20%) | |
93 | Bool_t frameCrosses=kTRUE; | |
94 | ||
95 | // Rotation matrices in the x-y plane | |
96 | // phi= 0 deg | |
97 | Int_t irot1; | |
98 | fMUON->AliMatrix(irot1, 90., 0., 90., 90., 0., 0.); | |
99 | // phi= 90 deg | |
100 | Int_t irot2; | |
101 | fMUON->AliMatrix(irot2, 90., 90., 90., 180., 0., 0.); | |
102 | ||
103 | AliMUONChamber* iChamber1 = GetChamber(0); | |
104 | AliMUONChamber* iChamber2 = GetChamber(1); | |
105 | AliMUONChamber* iChamber = iChamber1; | |
106 | ||
107 | // DGas decreased from standard one (0.5) | |
108 | iChamber->SetDGas(0.4); | |
109 | iChamber2->SetDGas(0.4); | |
110 | ||
111 | // DAlu increased from standard one (3% of X0), | |
112 | // because more electronics with smaller pads | |
113 | iChamber->SetDAlu(3.5 * 8.9 / 100.); | |
114 | iChamber2->SetDAlu(3.5 * 8.9 / 100.); | |
115 | ||
116 | // Half of the total thickness of frame crosses (including DAlu) | |
117 | // for each chamber in stations 1 and 2: | |
118 | // 3% of X0 of composite material, | |
119 | // but taken as Aluminium here, with same thickness in number of X0 | |
120 | Float_t dframez = 3. * 8.9 / 100; | |
121 | Float_t zfpos=-(iChamber->DGas()+dframez+iChamber->DAlu())/2; | |
122 | // The same parameters are defined in builder for station 2 | |
123 | ||
124 | // Mother volume | |
125 | // Outer excess and inner recess for mother volume radius | |
126 | // with respect to ROuter and RInner | |
127 | Float_t dframep=.001; // Value for station 3 should be 6 ... | |
128 | // Width (RdPhi) of the frame crosses for stations 1 and 2 (cm) | |
129 | // Float_t dframep1=.001; | |
130 | Float_t dframep1 = 11.0; | |
131 | Float_t phi=2*TMath::Pi()/12/2; | |
132 | // The same parameters are defined in builder for station 2 | |
133 | ||
134 | Float_t tpar[3]; | |
135 | Float_t dstation = (-iChamber2->Z()) - (-iChamber1->Z()); | |
136 | tpar[0] = iChamber->RInner()-dframep; | |
137 | tpar[1] = (iChamber->ROuter()+dframep)/TMath::Cos(phi); | |
138 | tpar[2] = dstation/5; | |
139 | ||
140 | gMC->Gsvolu("S01M", "TUBE", idAir, tpar, 3); | |
141 | gMC->Gsvolu("S02M", "TUBE", idAir, tpar, 3); | |
142 | ||
143 | // CHANGED | |
144 | //gMC->Gspos("S01M", 1, "ALIC", 0., 0., zpos1 , 0, "ONLY"); | |
a432117a | 145 | //gMC->Gspos("S02M", 1, "ALIC", 0., 0., zpos2 , 0, "ONLY"); |
146 | ||
147 | GetEnvelopes(0)->AddEnvelope("S01M", 100, false); | |
148 | GetEnvelopes(1)->AddEnvelope("S02M", 200, false); | |
149 | ||
d1cd2474 | 150 | |
151 | // // Aluminium frames | |
152 | // // Outer frames | |
153 | // pgpar[0] = 360/12/2; | |
154 | // pgpar[1] = 360.; | |
155 | // pgpar[2] = 12.; | |
156 | // pgpar[3] = 2; | |
157 | // pgpar[4] = -dframez/2; | |
158 | // pgpar[5] = iChamber->ROuter(); | |
159 | // pgpar[6] = pgpar[5]+dframep1; | |
160 | // pgpar[7] = +dframez/2; | |
161 | // pgpar[8] = pgpar[5]; | |
162 | // pgpar[9] = pgpar[6]; | |
163 | // gMC->Gsvolu("S01O", "PGON", idAlu1, pgpar, 10); | |
164 | // gMC->Gsvolu("S02O", "PGON", idAlu1, pgpar, 10); | |
165 | // gMC->Gspos("S01O",1,"S01M", 0.,0.,-zfpos, 0,"ONLY"); | |
166 | // gMC->Gspos("S01O",2,"S01M", 0.,0.,+zfpos, 0,"ONLY"); | |
167 | // gMC->Gspos("S02O",1,"S02M", 0.,0.,-zfpos, 0,"ONLY"); | |
168 | // gMC->Gspos("S02O",2,"S02M", 0.,0.,+zfpos, 0,"ONLY"); | |
169 | // // | |
170 | // // Inner frame | |
171 | // tpar[0]= iChamber->RInner()-dframep1; | |
172 | // tpar[1]= iChamber->RInner(); | |
173 | // tpar[2]= dframez/2; | |
174 | // gMC->Gsvolu("S01I", "TUBE", idAlu1, tpar, 3); | |
175 | // gMC->Gsvolu("S02I", "TUBE", idAlu1, tpar, 3); | |
176 | ||
177 | // gMC->Gspos("S01I",1,"S01M", 0.,0.,-zfpos, 0,"ONLY"); | |
178 | // gMC->Gspos("S01I",2,"S01M", 0.,0.,+zfpos, 0,"ONLY"); | |
179 | // gMC->Gspos("S02I",1,"S02M", 0.,0.,-zfpos, 0,"ONLY"); | |
180 | // gMC->Gspos("S02I",2,"S02M", 0.,0.,+zfpos, 0,"ONLY"); | |
181 | // | |
182 | // Frame Crosses | |
183 | if (frameCrosses) { | |
184 | // outside gas | |
185 | // security for inside mother volume | |
186 | Float_t bpar[3]; | |
187 | bpar[0] = (iChamber->ROuter() - iChamber->RInner()) | |
188 | * TMath::Cos(TMath::ASin(dframep1 / | |
189 | (iChamber->ROuter() - iChamber->RInner()))) | |
190 | / 2.0; | |
191 | bpar[1] = dframep1/2; | |
192 | // total thickness will be (4 * bpar[2]) for each chamber, | |
193 | // which has to be equal to (2 * dframez) - DAlu | |
194 | bpar[2] = (2.0 * dframez - iChamber->DAlu()) / 4.0; | |
195 | gMC->Gsvolu("S01B", "BOX", idAlu1, bpar, 3); | |
196 | gMC->Gsvolu("S02B", "BOX", idAlu1, bpar, 3); | |
197 | ||
198 | gMC->Gspos("S01B",1,"S01M", +iChamber->RInner()+bpar[0] , 0,-zfpos, | |
199 | irot1,"ONLY"); | |
200 | gMC->Gspos("S01B",2,"S01M", -iChamber->RInner()-bpar[0] , 0,-zfpos, | |
201 | irot1,"ONLY"); | |
202 | gMC->Gspos("S01B",3,"S01M", 0, +iChamber->RInner()+bpar[0] ,-zfpos, | |
203 | irot2,"ONLY"); | |
204 | gMC->Gspos("S01B",4,"S01M", 0, -iChamber->RInner()-bpar[0] ,-zfpos, | |
205 | irot2,"ONLY"); | |
206 | gMC->Gspos("S01B",5,"S01M", +iChamber->RInner()+bpar[0] , 0,+zfpos, | |
207 | irot1,"ONLY"); | |
208 | gMC->Gspos("S01B",6,"S01M", -iChamber->RInner()-bpar[0] , 0,+zfpos, | |
209 | irot1,"ONLY"); | |
210 | gMC->Gspos("S01B",7,"S01M", 0, +iChamber->RInner()+bpar[0] ,+zfpos, | |
211 | irot2,"ONLY"); | |
212 | gMC->Gspos("S01B",8,"S01M", 0, -iChamber->RInner()-bpar[0] ,+zfpos, | |
213 | irot2,"ONLY"); | |
214 | ||
215 | gMC->Gspos("S02B",1,"S02M", +iChamber->RInner()+bpar[0] , 0,-zfpos, | |
216 | irot1,"ONLY"); | |
217 | gMC->Gspos("S02B",2,"S02M", -iChamber->RInner()-bpar[0] , 0,-zfpos, | |
218 | irot1,"ONLY"); | |
219 | gMC->Gspos("S02B",3,"S02M", 0, +iChamber->RInner()+bpar[0] ,-zfpos, | |
220 | irot2,"ONLY"); | |
221 | gMC->Gspos("S02B",4,"S02M", 0, -iChamber->RInner()-bpar[0] ,-zfpos, | |
222 | irot2,"ONLY"); | |
223 | gMC->Gspos("S02B",5,"S02M", +iChamber->RInner()+bpar[0] , 0,+zfpos, | |
224 | irot1,"ONLY"); | |
225 | gMC->Gspos("S02B",6,"S02M", -iChamber->RInner()-bpar[0] , 0,+zfpos, | |
226 | irot1,"ONLY"); | |
227 | gMC->Gspos("S02B",7,"S02M", 0, +iChamber->RInner()+bpar[0] ,+zfpos, | |
228 | irot2,"ONLY"); | |
229 | gMC->Gspos("S02B",8,"S02M", 0, -iChamber->RInner()-bpar[0] ,+zfpos, | |
230 | irot2,"ONLY"); | |
231 | } | |
232 | // | |
233 | // Chamber Material represented by Alu sheet | |
234 | tpar[0]= iChamber->RInner(); | |
235 | tpar[1]= iChamber->ROuter(); | |
236 | tpar[2] = (iChamber->DGas()+iChamber->DAlu())/2; | |
237 | gMC->Gsvolu("S01A", "TUBE", idAlu2, tpar, 3); | |
238 | gMC->Gsvolu("S02A", "TUBE",idAlu2, tpar, 3); | |
239 | gMC->Gspos("S01A", 1, "S01M", 0., 0., 0., 0, "ONLY"); | |
240 | gMC->Gspos("S02A", 1, "S02M", 0., 0., 0., 0, "ONLY"); | |
241 | // | |
242 | // Sensitive volumes | |
243 | // tpar[2] = iChamber->DGas(); | |
244 | tpar[2] = iChamber->DGas()/2; | |
245 | gMC->Gsvolu("S01G", "TUBE", idGas, tpar, 3); | |
246 | gMC->Gsvolu("S02G", "TUBE", idGas, tpar, 3); | |
247 | gMC->Gspos("S01G", 1, "S01A", 0., 0., 0., 0, "ONLY"); | |
248 | gMC->Gspos("S02G", 1, "S02A", 0., 0., 0., 0, "ONLY"); | |
249 | // | |
250 | // Frame Crosses to be placed inside gas | |
251 | // NONE: chambers are sensitive everywhere | |
252 | // if (frameCrosses) { | |
253 | ||
254 | // dr = (iChamber->ROuter() - iChamber->RInner()); | |
255 | // bpar[0] = TMath::Sqrt(dr*dr-dframep1*dframep1/4)/2; | |
256 | // bpar[1] = dframep1/2; | |
257 | // bpar[2] = iChamber->DGas()/2; | |
258 | // gMC->Gsvolu("S01F", "BOX", idAlu1, bpar, 3); | |
259 | // gMC->Gsvolu("S02F", "BOX", idAlu1, bpar, 3); | |
260 | ||
261 | // gMC->Gspos("S01F",1,"S01G", +iChamber->RInner()+bpar[0] , 0, 0, | |
262 | // irot1,"ONLY"); | |
263 | // gMC->Gspos("S01F",2,"S01G", -iChamber->RInner()-bpar[0] , 0, 0, | |
264 | // irot1,"ONLY"); | |
265 | // gMC->Gspos("S01F",3,"S01G", 0, +iChamber->RInner()+bpar[0] , 0, | |
266 | // irot2,"ONLY"); | |
267 | // gMC->Gspos("S01F",4,"S01G", 0, -iChamber->RInner()-bpar[0] , 0, | |
268 | // irot2,"ONLY"); | |
269 | ||
270 | // gMC->Gspos("S02F",1,"S02G", +iChamber->RInner()+bpar[0] , 0, 0, | |
271 | // irot1,"ONLY"); | |
272 | // gMC->Gspos("S02F",2,"S02G", -iChamber->RInner()-bpar[0] , 0, 0, | |
273 | // irot1,"ONLY"); | |
274 | // gMC->Gspos("S02F",3,"S02G", 0, +iChamber->RInner()+bpar[0] , 0, | |
275 | // irot2,"ONLY"); | |
276 | // gMC->Gspos("S02F",4,"S02G", 0, -iChamber->RInner()-bpar[0] , 0, | |
277 | // irot2,"ONLY"); | |
278 | // } | |
279 | } | |
280 | ||
281 | //______________________________________________________________________________ | |
282 | void AliMUONSt1GeometryBuilder::SetTransformations() | |
283 | { | |
284 | // Defines the transformations for the station2 chambers. | |
285 | // --- | |
286 | ||
287 | AliMUONChamber* iChamber1 = GetChamber(0); | |
288 | Double_t zpos1= - iChamber1->Z(); | |
289 | iChamber1->GetGeometry() | |
290 | ->SetTranslation(TGeoTranslation(0., 0., zpos1)); | |
291 | ||
292 | AliMUONChamber* iChamber2 = GetChamber(1); | |
293 | Double_t zpos2 = - iChamber2->Z(); | |
294 | iChamber2->GetGeometry() | |
295 | ->SetTranslation(TGeoTranslation(0., 0., zpos2)); | |
296 | } | |
297 | ||
298 | //______________________________________________________________________________ | |
299 | void AliMUONSt1GeometryBuilder::SetSensitiveVolumes() | |
300 | { | |
301 | // Defines the sensitive volumes for station1 chambers. | |
302 | // --- | |
303 | ||
304 | GetChamber(0)->GetGeometry()->SetSensitiveVolume("S01G"); | |
305 | GetChamber(1)->GetGeometry()->SetSensitiveVolume("S02G"); | |
306 | } |