]>
Commit | Line | Data |
---|---|---|
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 | // $Id$ | |
17 | ||
18 | //----------------------------------------------------------------------------- | |
19 | // Class AliMUONSt2GeometryBuilderV2 | |
20 | // --------------------------------- | |
21 | // MUON Station2 coarse geometry construction class. | |
22 | // Author: SANJOY PAL ,Prof. SUKALYAN CHATTOPADHAYAY [SINP, KOLKATA] | |
23 | // & Dr.SHAKEEL AHMAD (AMU), INDIA | |
24 | //----------------------------------------------------------------------------- | |
25 | ||
26 | ||
27 | #include <TVirtualMC.h> | |
28 | #include <TGeoMatrix.h> | |
29 | #include <Riostream.h> | |
30 | #include <TArrayI.h> | |
31 | ||
32 | #include "AliRun.h" | |
33 | #include "AliLog.h" | |
34 | ||
35 | #include "AliMUONSt2GeometryBuilderV2.h" | |
36 | #include "AliMUON.h" | |
37 | #include "AliMUONConstants.h" | |
38 | #include "AliMUONGeometryModule.h" | |
39 | #include "AliMUONGeometryEnvelopeStore.h" | |
40 | #include "AliMUONConstants.h" | |
41 | ||
42 | #define PI 3.14159 | |
43 | ||
44 | /// \cond CLASSIMP | |
45 | ClassImp(AliMUONSt2GeometryBuilderV2) | |
46 | /// \endcond | |
47 | ||
48 | //______________________________________________________________________________ | |
49 | AliMUONSt2GeometryBuilderV2::AliMUONSt2GeometryBuilderV2(AliMUON* muon) | |
50 | : AliMUONVGeometryBuilder(2, 2), | |
51 | fMUON(muon) | |
52 | { | |
53 | /// Standard constructor | |
54 | ||
55 | } | |
56 | ||
57 | //______________________________________________________________________________ | |
58 | AliMUONSt2GeometryBuilderV2::AliMUONSt2GeometryBuilderV2() | |
59 | : AliMUONVGeometryBuilder(), | |
60 | fMUON(0) | |
61 | { | |
62 | /// Default constructor | |
63 | } | |
64 | ||
65 | //______________________________________________________________________________ | |
66 | AliMUONSt2GeometryBuilderV2::~AliMUONSt2GeometryBuilderV2() | |
67 | { | |
68 | /// Destructor | |
69 | } | |
70 | ||
71 | // | |
72 | // public methods | |
73 | // | |
74 | ||
75 | //______________________________________________________________________________ | |
76 | void AliMUONSt2GeometryBuilderV2::CreateGeometry() | |
77 | { | |
78 | /// Geometry construction | |
79 | ||
80 | // | |
81 | //******************************************************************** | |
82 | // Station 2 ** | |
83 | //******************************************************************** | |
84 | // indices 1 and 2 for first and second chambers in the station | |
85 | // iChamber (first chamber) kept for other quanties than Z, | |
86 | // assumed to be the same in both chambers | |
87 | ||
88 | // Get tracking medias Ids | |
89 | Int_t *idtmed = fMUON->GetIdtmed()->GetArray()-1099; | |
90 | Int_t idAir = idtmed[1100]; // medium 1 | |
91 | Int_t idGas = idtmed[1108]; // medium Ar-CO2 gas (80%+20%) | |
92 | Int_t idPCB = idtmed[1122]; // medium FR4 | |
93 | Int_t idCU = idtmed[1110]; // medium copper | |
94 | Int_t idRoha = idtmed[1113]; // medium roha cell | |
95 | Int_t idPGF30= idtmed[1123]; // medium for Frame Eq.to Bakelite | |
96 | //Int_t idScru = idtmed[1128]; // screw material - Stainless Steel(18%Cr,9%Ni,Fe) | |
97 | ||
98 | ||
99 | ||
100 | /*######################################################################################## | |
101 | Create volume for one Quadrant | |
102 | ##########################################################################################*/ | |
103 | Float_t tpar1[5]; | |
104 | tpar1[0] = 20.55; | |
105 | tpar1[1] = 123.5; | |
106 | tpar1[2] = 6.8/2; | |
107 | tpar1[3] = -12.0; | |
108 | tpar1[4] = 102.0; | |
109 | ||
110 | ||
111 | gMC->Gsvolu("SQM3","TUBS", idAir, tpar1, 5); | |
112 | gMC->Gsvolu("SQM4","TUBS", idAir, tpar1, 5); | |
113 | ||
114 | ||
115 | //================================================================================== | |
116 | // Plane | |
117 | //================================================================================== | |
118 | ||
119 | //Thickness of variour parts | |
120 | Float_t zCu = 0.005; // eff. cu in cathode pcb | |
121 | Float_t zCbb = 0.04; //cathode pcb | |
122 | Float_t zRoha = 2.5; // Rhocell | |
123 | Float_t zMeb = 0.04; //Mech. exit board //0.08 | |
124 | Float_t zEeb = 0.04; //Effective electronic readout board //0.02 | |
125 | ||
126 | //starting Z-positions of various parts--- in Plane-1 | |
127 | Float_t zposCu = 0.25; // 2.5 mm => gap between anode & chatode plane | |
128 | Float_t zposCbb = zposCu + zCu; | |
129 | Float_t zposRoha = zposCbb + zCbb; | |
130 | Float_t zposMeb = zposRoha + zRoha; | |
131 | Float_t zposEeb = zposMeb + zMeb; | |
132 | ||
133 | Float_t zposCuBar = zposCu + zCu/2.; //for segment 0 & 6 | |
134 | Float_t zposCbbBar = zposCbb + zCbb/2.; | |
135 | Float_t zposRohaBar = zposRoha + zRoha/2.; | |
136 | Float_t zposMebBar = zposMeb + zMeb/2.; | |
137 | Float_t zposEebBar = zposEeb + zEeb/2.; | |
138 | ||
139 | ||
140 | //Cathode PCB + Copper sheet + Rohacell + mech exit board + eff. electronic exit board | |
141 | ||
142 | //Segment-0 ~~~Horizantal box ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ | |
143 | Float_t bparH[3]; | |
144 | bparH[0] = 95.5/2.; // extension beyond 0 deg in x direction // 94.5 | |
145 | bparH[1] = 1.2/2.; // extension beyond 0 deg in y direction (3.7[total extn] - 2.5[frame dim]) | |
146 | bparH[2] = zCu/2.; //Thickness of Copper sheet | |
147 | gMC->Gsvolu("SCU0L", "BOX", idCU, bparH, 3); | |
148 | ||
149 | bparH[2] = zCbb/2.; // thickness of cathode sheet in z direction | |
150 | gMC->Gsvolu("SCB0L", "BOX", idPCB, bparH, 3); | |
151 | ||
152 | bparH[2] = zRoha/2.; //Thickness of Roha cell | |
153 | gMC->Gsvolu("SRH0L", "BOX", idRoha, bparH, 3); | |
154 | ||
155 | bparH[2] = zMeb/2; //Thickness of mechanical exit board | |
156 | gMC->Gsvolu("SMB0L", "BOX", idPCB, bparH, 3); | |
157 | ||
158 | bparH[2] = zEeb/2; //Thickness of effective electronic exit board | |
159 | gMC->Gsvolu("SEB0L", "BOX", idCU, bparH, 3); | |
160 | ||
161 | //Segment-1 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ | |
162 | Float_t pgpar[10]; // polygon | |
163 | pgpar[0] = 0.; // initial angle | |
164 | pgpar[1] = 90.; // increment in angle starting from initial angle | |
165 | pgpar[2] = 5; // number of side | |
166 | pgpar[3] = 2.; // number of plane | |
167 | pgpar[4] = -zCu/2.; // z-position of the first plane | |
168 | pgpar[5] = 23.1; // innner radius first plane | |
169 | pgpar[6] = 117.6; // outer radious first plane | |
170 | pgpar[7] = zCu/2.; // z-position of the second plane | |
171 | pgpar[8] = pgpar[5]; // innner radius of second plane | |
172 | pgpar[9] = pgpar[6]; // outer radious of second plane | |
173 | gMC->Gsvolu("SCU1L", "PGON", idCU, pgpar, 10); | |
174 | ||
175 | pgpar[4] = -zCbb/2.; // z-position of the first plane | |
176 | pgpar[7] = zCbb/2.; // Thickness of copper-sheet | |
177 | gMC->Gsvolu("SCB1L", "PGON", idPCB, pgpar, 10); | |
178 | ||
179 | pgpar[4] = -zRoha/2.; // Thickness of Roha cell | |
180 | pgpar[7] = zRoha/2.; // Thickness of Roha cell | |
181 | gMC->Gsvolu("SRH1L", "PGON", idRoha, pgpar, 10); | |
182 | ||
183 | pgpar[4] = -zMeb/2.; // Thickness of mechanical exit board | |
184 | pgpar[7] = zMeb/2.; // Thickness of mechanical exit board | |
185 | gMC->Gsvolu("SMB1L", "PGON", idPCB, pgpar, 10); | |
186 | ||
187 | pgpar[4] = -zEeb/2.; // Thickness of electronic readout board | |
188 | pgpar[7] = zEeb/2.; // Thickness of electronic readout board | |
189 | gMC->Gsvolu("SEB1L", "PGON", idCU, pgpar, 10); | |
190 | ||
191 | ||
192 | //Segment-2 - vertical box (simalar to horizontal bar as in Segment 0)~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ | |
193 | ||
194 | Float_t bparV[3]; | |
195 | ||
196 | bparV[0] = 1.0/2.; | |
197 | bparV[1] = 95.5/2.; // 94.5 | |
198 | bparV[2] = zCu/2.; | |
199 | gMC->Gsvolu("SCU2L", "BOX", idCU, bparV, 3); | |
200 | ||
201 | bparV[2] = zCbb/2.; | |
202 | gMC->Gsvolu("SCB2L", "BOX", idPCB, bparV, 3); | |
203 | ||
204 | ||
205 | bparV[2] = zRoha/2.; | |
206 | gMC->Gsvolu("SRH2L", "BOX", idRoha, bparV, 3); | |
207 | ||
208 | bparV[2] = zMeb/2; | |
209 | gMC->Gsvolu("SMB2L", "BOX", idPCB, bparV, 3); | |
210 | ||
211 | bparV[2] = zEeb/2; | |
212 | gMC->Gsvolu("SEB2L", "BOX", idCU, bparV, 3); | |
213 | ||
214 | //....(Setting posion of Segment 0,1,2).................................................................. | |
215 | ||
216 | Float_t xposHorBox = bparH[0] + 23.1; // 23.1 = 20.6(inner radius of qrd) + 2.5 (width of frame) | |
217 | Float_t yposHorBox = -bparH[1]; | |
218 | ||
219 | Float_t xposVerBox = -bparV[0]; | |
220 | Float_t yposVerBox = bparV[1] + 23.1; | |
221 | ||
222 | ||
223 | //Positioning the PCB | |
224 | ||
225 | ||
226 | // chamber 3 | |
227 | gMC->Gspos("SCB0L",1, "SQM3",xposHorBox,yposHorBox,zposCbbBar,0, "ONLY"); | |
228 | gMC->Gspos("SCB1L",1, "SQM3", 0.0,0.0,zposCbbBar,0, "ONLY"); | |
229 | gMC->Gspos("SCB2L",1, "SQM3",xposVerBox,yposVerBox,zposCbbBar,0, "ONLY"); | |
230 | ||
231 | gMC->Gspos("SCB0L",2, "SQM3",xposHorBox,yposHorBox,-zposCbbBar,0, "ONLY"); | |
232 | gMC->Gspos("SCB1L",2, "SQM3", 0.0,0.0,-zposCbbBar,0, "ONLY"); | |
233 | gMC->Gspos("SCB2L",2, "SQM3",xposVerBox,yposVerBox,-zposCbbBar,0, "ONLY"); | |
234 | ||
235 | // chamber 4 | |
236 | gMC->Gspos("SCB0L",3, "SQM4",xposHorBox,yposHorBox,zposCbbBar,0, "ONLY"); | |
237 | gMC->Gspos("SCB1L",3, "SQM4", 0.0,0.0,zposCbbBar,0, "ONLY"); | |
238 | gMC->Gspos("SCB2L",3, "SQM4",xposVerBox,yposVerBox,zposCbbBar,0, "ONLY"); | |
239 | ||
240 | gMC->Gspos("SCB0L",4, "SQM4",xposHorBox,yposHorBox,-zposCbbBar,0, "ONLY"); | |
241 | gMC->Gspos("SCB1L",4, "SQM4", 0.0,0.0,-zposCbbBar,0, "ONLY"); | |
242 | gMC->Gspos("SCB2L",4, "SQM4",xposVerBox,yposVerBox,-zposCbbBar,0, "ONLY"); | |
243 | ||
244 | ||
245 | //Positioning Copper sheet | |
246 | ||
247 | // chamber 3 | |
248 | ||
249 | gMC->Gspos("SCU0L",1, "SQM3",xposHorBox,yposHorBox,zposCuBar,0, "ONLY"); | |
250 | gMC->Gspos("SCU1L",1, "SQM3", 0.0,0.0,zposCuBar,0, "ONLY"); | |
251 | gMC->Gspos("SCU2L",1, "SQM3",xposVerBox,yposVerBox,zposCuBar,0, "ONLY"); | |
252 | ||
253 | gMC->Gspos("SCU0L",2, "SQM3",xposHorBox,yposHorBox,-zposCuBar,0, "ONLY"); | |
254 | gMC->Gspos("SCU1L",2, "SQM3", 0.0,0.0,-zposCuBar,0, "ONLY"); | |
255 | gMC->Gspos("SCU2L",2, "SQM3",xposVerBox,yposVerBox,-zposCuBar,0, "ONLY"); | |
256 | ||
257 | // chamber 4 | |
258 | gMC->Gspos("SCU0L",3, "SQM4",xposHorBox,yposHorBox,zposCuBar,0, "ONLY"); | |
259 | gMC->Gspos("SCU1L",3, "SQM4", 0.0,0.0,zposCuBar,0, "ONLY"); | |
260 | gMC->Gspos("SCU2L",3, "SQM4",xposVerBox,yposVerBox,zposCuBar,0, "ONLY"); | |
261 | ||
262 | gMC->Gspos("SCU0L",4, "SQM4",xposHorBox,yposHorBox,-zposCuBar,0, "ONLY"); | |
263 | gMC->Gspos("SCU1L",4, "SQM4", 0.0,0.0,-zposCuBar,0, "ONLY"); | |
264 | gMC->Gspos("SCU2L",4, "SQM4",xposVerBox,yposVerBox,-zposCuBar,0, "ONLY"); | |
265 | ||
266 | //Positioning Roha cell | |
267 | ||
268 | // chamber 3 | |
269 | ||
270 | gMC->Gspos("SRH0L",1, "SQM3",xposHorBox,yposHorBox,zposRohaBar,0, "ONLY"); | |
271 | gMC->Gspos("SRH1L",1, "SQM3", 0.0,0.0,zposRohaBar,0, "ONLY"); | |
272 | gMC->Gspos("SRH2L",1, "SQM3",xposVerBox,yposVerBox,zposRohaBar,0, "ONLY"); | |
273 | ||
274 | gMC->Gspos("SRH0L",2, "SQM3",xposHorBox,yposHorBox,-zposRohaBar,0, "ONLY"); | |
275 | gMC->Gspos("SRH1L",2, "SQM3", 0.0,0.0,-zposRohaBar,0, "ONLY"); | |
276 | gMC->Gspos("SRH2L",2, "SQM3",xposVerBox,yposVerBox,-zposRohaBar,0, "ONLY"); | |
277 | ||
278 | // chamber 4 | |
279 | ||
280 | gMC->Gspos("SRH0L",3, "SQM4",xposHorBox,yposHorBox,zposRohaBar,0, "ONLY"); | |
281 | gMC->Gspos("SRH1L",3, "SQM4", 0.0,0.0,zposRohaBar,0, "ONLY"); | |
282 | gMC->Gspos("SRH2L",3, "SQM4",xposVerBox,yposVerBox,zposRohaBar,0, "ONLY"); | |
283 | ||
284 | gMC->Gspos("SRH0L",4, "SQM4",xposHorBox,yposHorBox,-zposRohaBar,0, "ONLY"); | |
285 | gMC->Gspos("SRH1L",4, "SQM4", 0.0,0.0,-zposRohaBar,0, "ONLY"); | |
286 | gMC->Gspos("SRH2L",4, "SQM4",xposVerBox,yposVerBox,-zposRohaBar,0, "ONLY"); | |
287 | ||
288 | //Positioning of Mech. exit board | |
289 | ||
290 | // chamber 3 | |
291 | ||
292 | gMC->Gspos("SMB0L",1, "SQM3",xposHorBox,yposHorBox,zposMebBar,0, "ONLY"); | |
293 | gMC->Gspos("SMB1L",1, "SQM3", 0.0,0.0,zposMebBar,0, "ONLY"); | |
294 | gMC->Gspos("SMB2L",1, "SQM3",xposVerBox,yposVerBox,zposMebBar,0, "ONLY"); | |
295 | ||
296 | gMC->Gspos("SMB0L",2, "SQM3",xposHorBox,yposHorBox,-zposMebBar,0, "ONLY"); | |
297 | gMC->Gspos("SMB1L",2, "SQM3", 0.0,0.0,-zposMebBar,0, "ONLY"); | |
298 | gMC->Gspos("SMB2L",2, "SQM3",xposVerBox,yposVerBox,-zposMebBar,0, "ONLY"); | |
299 | // chamber 4 | |
300 | ||
301 | gMC->Gspos("SMB0L",3, "SQM4",xposHorBox,yposHorBox,zposMebBar,0, "ONLY"); | |
302 | gMC->Gspos("SMB1L",3, "SQM4", 0.0,0.0,zposMebBar,0, "ONLY"); | |
303 | gMC->Gspos("SMB2L",3, "SQM4",xposVerBox,yposVerBox,zposMebBar,0, "ONLY"); | |
304 | ||
305 | gMC->Gspos("SMB0L",4, "SQM4",xposHorBox,yposHorBox,-zposMebBar,0, "ONLY"); | |
306 | gMC->Gspos("SMB1L",4, "SQM4", 0.0,0.0,-zposMebBar,0, "ONLY"); | |
307 | gMC->Gspos("SMB2L",4, "SQM4",xposVerBox,yposVerBox,-zposMebBar,0, "ONLY"); | |
308 | ||
309 | //Positioning of Electronic exit board | |
310 | ||
311 | ||
312 | // chamber 3 | |
313 | ||
314 | gMC->Gspos("SEB0L",1, "SQM3",xposHorBox,yposHorBox,zposEebBar,0, "ONLY"); | |
315 | gMC->Gspos("SEB1L",1, "SQM3", 0.0,0.0,zposEebBar,0, "ONLY"); | |
316 | gMC->Gspos("SEB2L",1, "SQM3",xposVerBox,yposVerBox,zposEebBar,0, "ONLY"); | |
317 | ||
318 | gMC->Gspos("SEB0L",2, "SQM3",xposHorBox,yposHorBox,-zposEebBar,0, "ONLY"); | |
319 | gMC->Gspos("SEB1L",2, "SQM3", 0.0,0.0,-zposEebBar,0, "ONLY"); | |
320 | gMC->Gspos("SEB2L",2, "SQM3",xposVerBox,yposVerBox,-zposEebBar,0, "ONLY"); | |
321 | // chamber 4 | |
322 | ||
323 | gMC->Gspos("SEB0L",3, "SQM4",xposHorBox,yposHorBox,zposEebBar,0, "ONLY"); | |
324 | gMC->Gspos("SEB1L",3, "SQM4", 0.0,0.0,zposEebBar,0, "ONLY"); | |
325 | gMC->Gspos("SEB2L",3, "SQM4",xposVerBox,yposVerBox,zposEebBar,0, "ONLY"); | |
326 | ||
327 | gMC->Gspos("SEB0L",4, "SQM4",xposHorBox,yposHorBox,-zposEebBar,0, "ONLY"); | |
328 | gMC->Gspos("SEB1L",4, "SQM4", 0.0,0.0,-zposEebBar,0, "ONLY"); | |
329 | gMC->Gspos("SEB2L",4, "SQM4",xposVerBox,yposVerBox,-zposEebBar,0, "ONLY"); | |
330 | ||
331 | ||
332 | //---------------------------------------------------------------------- | |
333 | // Frames | |
334 | //---------------------------------------------------------------------- | |
335 | //Frame-1 | |
336 | ||
337 | Float_t frame1[3] ; | |
338 | frame1[0] = 101.0/2.; //100.6 = 94.5 + 2.5 + 3.6 | |
339 | frame1[1] = 2.5/2.; | |
340 | frame1[2] = 5.0/2.; | |
341 | ||
342 | gMC->Gsvolu("SFRM1", "BOX", idPGF30, frame1, 3); //Frame - 1 // fill with pkk GF30 | |
343 | ||
344 | Float_t arib1[3]; | |
345 | arib1[0] = frame1[0]; | |
346 | arib1[1] = 0.9/2.; | |
347 | arib1[2] =(frame1[2]-0.95)/2.0; | |
348 | ||
349 | gMC->Gsvolu("SFRA1", "BOX", idAir, arib1, 3); // fill with air | |
350 | ||
351 | Float_t xposarib1 = 0; | |
352 | Float_t yposarib1 = -frame1[1] + arib1[1]; | |
353 | Float_t zposarib1 = frame1[2] - arib1[2]; | |
354 | ||
355 | gMC->Gspos("SFRA1",1, "SFRM1", xposarib1, yposarib1, zposarib1,0, "ONLY"); //replace pkk GF30 with air(b) | |
356 | gMC->Gspos("SFRA1",2, "SFRM1", xposarib1, yposarib1, -zposarib1,0, "ONLY"); //replace pkk GF30 with air(nb) | |
357 | ||
358 | Float_t rrib1[3]; | |
359 | rrib1[0] = frame1[0]; | |
360 | rrib1[1] = 0.6/2.; | |
361 | rrib1[2] =(frame1[2]-0.95)/2.0; | |
362 | ||
363 | gMC->Gsvolu("SFRR1", "BOX", idRoha, rrib1, 3); // fill with rohacell | |
364 | ||
365 | Float_t xposrrib1 = 0.0; | |
366 | Float_t yposrrib1 = frame1[1] - rrib1[1]; | |
367 | Float_t zposrrib1 = frame1[2] - rrib1[2]; | |
368 | ||
369 | gMC->Gspos("SFRR1",1, "SFRM1", xposrrib1, yposrrib1, zposrrib1,0, "ONLY");//replace pkk GF30 with rohacell | |
370 | gMC->Gspos("SFRR1",2, "SFRM1", xposrrib1, yposrrib1, -zposrrib1,0, "ONLY");//replace pkk GF30 with rohacell | |
371 | ||
372 | ||
373 | Float_t xposFr1 = frame1[0] + 20.6; | |
374 | Float_t yposFr1 = -3.7 + frame1[1] ; | |
375 | Float_t zposFr1 = 0.0; | |
376 | ||
377 | gMC->Gspos("SFRM1",1, "SQM3", xposFr1, yposFr1, zposFr1,0, "ONLY");// frame -1 | |
378 | gMC->Gspos("SFRM1",2, "SQM4", xposFr1, yposFr1, zposFr1,0, "ONLY");// frame -1 | |
379 | ||
380 | ||
381 | //...................................................................................... | |
382 | //Frame-2 | |
383 | ||
384 | Float_t frame2[3]; | |
385 | frame2[0] = 4.0/2.; | |
386 | frame2[1] = 1.2/2.; | |
387 | frame2[2] = 5.0/2; | |
388 | ||
389 | gMC->Gsvolu("SFRM2", "BOX", idPGF30, frame2, 3); //Frame - 2 | |
390 | ||
391 | Float_t rrib2[3]; | |
392 | rrib2[0] = frame2[0]-1.0/2.0; | |
393 | rrib2[1] = frame2[1]; | |
394 | rrib2[2] =(frame2[2]-0.95)/2.0; | |
395 | ||
396 | gMC->Gsvolu("SFRR2", "BOX", idRoha, rrib2, 3); | |
397 | ||
398 | Float_t xposrrib2 = -1.0/2.0; | |
399 | Float_t yposrrib2 = 0.0; | |
400 | Float_t zposrrib2 = frame2[2] - rrib2[2]; | |
401 | ||
402 | gMC->Gspos("SFRR2",1, "SFRM2", xposrrib2, yposrrib2, zposrrib2,0, "ONLY");//replace pkk GF30 with rohacell | |
403 | gMC->Gspos("SFRR2",2, "SFRM2", xposrrib2, yposrrib2, -zposrrib2,0, "ONLY");//replace pkk GF30 with roha | |
404 | ||
405 | ||
406 | ||
407 | Float_t xposFr2 = frame2[0] + 117.6; | |
408 | Float_t yposFr2 = -frame2[1]; | |
409 | Float_t zposFr2 = 0.0; | |
410 | ||
411 | gMC->Gspos("SFRM2",1, "SQM3", xposFr2, yposFr2, zposFr2,0, "MANY");//global positing of frame in SQM3 | |
412 | gMC->Gspos("SFRM2",2, "SQM4", xposFr2, yposFr2, zposFr2,0, "MANY");//global positing of frame in SQM4 | |
413 | ||
414 | ||
415 | //...................................................................................... | |
416 | ||
417 | //Frame-3 | |
418 | ||
419 | Float_t pgparFr3[10]; | |
420 | pgparFr3[0] = 0.; | |
421 | pgparFr3[1] = 90.; | |
422 | pgparFr3[2] = 5; | |
423 | pgparFr3[3] = 2.; | |
424 | pgparFr3[4] = 0.; | |
425 | pgparFr3[5] = 117.6; | |
426 | pgparFr3[6] = 121.6; | |
427 | pgparFr3[7] = pgparFr3[4] + 5.0; | |
428 | pgparFr3[8] = pgparFr3[5]; | |
429 | pgparFr3[9] = pgparFr3[6]; | |
430 | ||
431 | gMC->Gsvolu("SFRM3", "PGON", idPGF30, pgparFr3, 10); | |
432 | ||
433 | Float_t pgparRrib3[10]; | |
434 | pgparRrib3[0] = 0.; | |
435 | pgparRrib3[1] = 90.; | |
436 | pgparRrib3[2] = 5; | |
437 | pgparRrib3[3] = 2.; | |
438 | pgparRrib3[4] = 0.; | |
439 | pgparRrib3[5] = 117.6; | |
440 | pgparRrib3[6] = 120.6; | |
441 | pgparRrib3[7] = pgparRrib3[4] +1.55 ; | |
442 | pgparRrib3[8] = pgparRrib3[5]; | |
443 | pgparRrib3[9] = pgparRrib3[6]; | |
444 | ||
445 | gMC->Gsvolu("SFRR3", "PGON", idRoha, pgparRrib3, 10); | |
446 | ||
447 | Float_t xposrrib3 = 0.0; | |
448 | Float_t yposrrib3 = 0.0; | |
449 | Float_t zposrrib3 = 0.0; | |
450 | ||
451 | gMC->Gspos("SFRR3",1, "SFRM3", xposrrib3, yposrrib3, zposrrib3,0, "ONLY"); | |
452 | ||
453 | zposrrib3 = 3.45; | |
454 | ||
455 | gMC->Gspos("SFRR3",2, "SFRM3", xposrrib3, yposrrib3, zposrrib3,0, "ONLY"); | |
456 | ||
457 | ||
458 | ||
459 | Float_t xposFr3 = 0.0; | |
460 | Float_t yposFr3 = 0.0; | |
461 | Float_t zposFr3 = -frame1[2]; | |
462 | ||
463 | gMC->Gspos("SFRM3",1, "SQM3", xposFr3, yposFr3, zposFr3,0, "ONLY");// frame -1 | |
464 | gMC->Gspos("SFRM3",2, "SQM4", xposFr3, yposFr3, zposFr3,0, "ONLY");// frame -1 | |
465 | ||
466 | ||
467 | //...................................................................................... | |
468 | //Frame-4 | |
469 | ||
470 | Float_t frame4[3]; | |
471 | frame4[0] = 1.0/2.; | |
472 | frame4[1] = 4.0/2.; | |
473 | frame4[2] = frame1[2]; | |
474 | ||
475 | gMC->Gsvolu("SFRM4", "BOX", idPGF30, frame4, 3); | |
476 | ||
477 | Float_t rrib4[3]; | |
478 | rrib4[0] = frame4[0]; | |
479 | rrib4[1] = frame4[1]-1.0/2; | |
480 | rrib4[2] =(frame4[2]-0.95)/2.0; | |
481 | ||
482 | gMC->Gsvolu("SFRR4", "BOX", idRoha, rrib4, 3); | |
483 | ||
484 | Float_t xposrrib4 = 0.0; | |
485 | Float_t yposrrib4 = -1.0/2; | |
486 | Float_t zposrrib4 = frame4[2] - rrib4[2]; | |
487 | ||
488 | gMC->Gspos("SFRR4",1, "SFRM4", xposrrib4, yposrrib4, zposrrib4,0, "ONLY"); | |
489 | gMC->Gspos("SFRR4",2, "SFRM4", xposrrib4, yposrrib4, -zposrrib4,0, "ONLY"); | |
490 | ||
491 | ||
492 | ||
493 | Float_t xposFr4 = -frame4[0]; | |
494 | Float_t yposFr4 = -frame4[1] + 117.6; | |
495 | Float_t zposFr4 = 0.0; | |
496 | ||
497 | gMC->Gspos("SFRM4",1, "SQM3", xposFr4, yposFr4, zposFr4,0, "MANY"); | |
498 | gMC->Gspos("SFRM4",2, "SQM4", xposFr4, yposFr4, zposFr4,0, "MANY"); | |
499 | ||
500 | ||
501 | //...................................................................................... | |
502 | //Frame-5 | |
503 | Float_t frame5[3] ; | |
504 | frame5[0] = 2.7/2.; | |
505 | frame5[1] = 101.0/2.; | |
506 | frame5[2] = 5.0/2.; | |
507 | ||
508 | gMC->Gsvolu("SFRM5", "BOX", idPGF30, frame5, 3); //Frame - 1 | |
509 | ||
510 | Float_t arib5[3]; | |
511 | arib5[0] = 0.9/2.0; | |
512 | arib5[1] = frame5[1]; | |
513 | arib5[2] = (frame5[2]-0.95)/2.0; | |
514 | ||
515 | gMC->Gsvolu("SFRA5", "BOX", idAir, arib5, 3); | |
516 | ||
517 | Float_t xposarib5 = -frame5[0] + arib5[0]; | |
518 | Float_t yposarib5 = 0.0; | |
519 | Float_t zposarib5 = frame5[2] - arib5[2]; | |
520 | ||
521 | gMC->Gspos("SFRA5",1, "SFRM5", xposarib5, yposarib5, zposarib5,0, "ONLY"); | |
522 | gMC->Gspos("SFRA5",2, "SFRM5", xposarib5, yposarib5, -zposarib5,0, "ONLY"); | |
523 | ||
524 | Float_t rrib5[3]; | |
525 | rrib5[0] = 0.8/2.0; | |
526 | rrib5[1] = frame5[1]; | |
527 | rrib5[2] = (frame5[2]-0.95)/2.0; | |
528 | ||
529 | gMC->Gsvolu("SFRR5", "BOX", idRoha, rrib5, 3); | |
530 | ||
531 | Float_t xposrrib5 = frame5[0] - rrib5[0]; | |
532 | Float_t yposrrib5 = 0.0; | |
533 | Float_t zposrrib5 = frame5[2] - rrib5[2]; | |
534 | ||
535 | gMC->Gspos("SFRR5",1, "SFRM5", xposrrib5, yposrrib5, zposrrib5,0, "ONLY"); | |
536 | gMC->Gspos("SFRR5",2, "SFRM5", xposrrib5, yposrrib5, -zposrrib5,0, "ONLY"); | |
537 | ||
538 | Float_t xposFr5 = -3.7 + frame5[0]; | |
539 | Float_t yposFr5 = frame5[1] + 20.6; | |
540 | Float_t zposFr5 = 0.0; | |
541 | ||
542 | gMC->Gspos("SFRM5",1, "SQM3", xposFr5, yposFr5, zposFr5,0, "ONLY"); | |
543 | gMC->Gspos("SFRM5",2, "SQM4", xposFr5, yposFr5, zposFr5,0, "ONLY"); | |
544 | ||
545 | //...................................................................................... | |
546 | //Frame -6 | |
547 | ||
548 | Float_t frame6[3]; | |
549 | frame6[0] = 1.0/2.; | |
550 | frame6[1] = 2.5/2.; | |
551 | frame6[2] = frame1[2]; | |
552 | ||
553 | gMC->Gsvolu("SFRM6", "BOX", idPGF30, frame6, 3); | |
554 | ||
555 | Float_t rrib6[3]; | |
556 | rrib6[0] = frame6[0]; | |
557 | rrib6[1] = 1.5/2.; | |
558 | rrib6[2] =(frame2[2]-0.95)/2.0; | |
559 | ||
560 | gMC->Gsvolu("SFRR6", "BOX", idRoha, rrib6, 3); | |
561 | ||
562 | Float_t xposrrib6 = 0.0; | |
563 | Float_t yposrrib6 = 1.0/2.0; | |
564 | Float_t zposrrib6 = frame6[2] - rrib6[2]; | |
565 | ||
566 | gMC->Gspos("SFRR6",1, "SFRM6", xposrrib6, yposrrib6, zposrrib6,0, "ONLY"); | |
567 | gMC->Gspos("SFRR6",2, "SFRM6", xposrrib6, yposrrib6, -zposrrib6,0, "ONLY"); | |
568 | ||
569 | ||
570 | ||
571 | Float_t xposFr6 = -frame6[0]; | |
572 | Float_t yposFr6 = frame6[1] + 20.6; | |
573 | Float_t zposFr6 = 0.0; | |
574 | ||
575 | gMC->Gspos("SFRM6",1, "SQM3", xposFr6, yposFr6, zposFr6,0, "ONLY"); | |
576 | gMC->Gspos("SFRM6",2, "SQM4", xposFr6, yposFr6, zposFr6,0, "ONLY"); | |
577 | ||
578 | ||
579 | //...................................................................................... | |
580 | //Frame - 7 inner pgon | |
581 | ||
582 | Float_t pgparFr7[10]; | |
583 | pgparFr7[0] = 0.; | |
584 | pgparFr7[1] = 90.; | |
585 | pgparFr7[2] = 5; | |
586 | pgparFr7[3] = 2.; | |
587 | pgparFr7[4] = 0.; | |
588 | pgparFr7[5] = 20.6; | |
589 | pgparFr7[6] = 23.1; | |
590 | pgparFr7[7] = pgparFr7[4] + 5.0; | |
591 | pgparFr7[8] = pgparFr7[5]; | |
592 | pgparFr7[9] = pgparFr7[6]; | |
593 | ||
594 | gMC->Gsvolu("SFRM7", "PGON", idPGF30, pgparFr7, 10); | |
595 | ||
596 | Float_t pgparRrib7[10]; | |
597 | pgparRrib7[0] = 0.; | |
598 | pgparRrib7[1] = 90.; | |
599 | pgparRrib7[2] = 5; | |
600 | pgparRrib7[3] = 2.; | |
601 | pgparRrib7[4] = 0.; | |
602 | pgparRrib7[5] = 21.6; | |
603 | pgparRrib7[6] = 23.1; | |
604 | pgparRrib7[7] = pgparRrib7[4] + 1.55; | |
605 | pgparRrib7[8] = pgparRrib7[5]; | |
606 | pgparRrib7[9] = pgparRrib7[6]; | |
607 | ||
608 | gMC->Gsvolu("SFRR7", "PGON", idRoha, pgparRrib7, 10); | |
609 | ||
610 | Float_t xposrrib7 = 0.0; | |
611 | Float_t yposrrib7 = 0.0; | |
612 | Float_t zposrrib7 = 0.0; | |
613 | ||
614 | gMC->Gspos("SFRR7",1, "SFRM7", xposrrib7, yposrrib7, zposrrib7,0, "ONLY"); | |
615 | ||
616 | zposrrib7 = 3.45; | |
617 | ||
618 | gMC->Gspos("SFRR7",2, "SFRM7", xposrrib7, yposrrib7, zposrrib7,0, "ONLY"); | |
619 | ||
620 | Float_t xposFr7 = 0.0; | |
621 | Float_t yposFr7 = 0.0; | |
622 | Float_t zposFr7 = -frame1[2]; | |
623 | ||
624 | gMC->Gspos("SFRM7",1, "SQM3", xposFr7, yposFr7, zposFr7,0, "ONLY"); | |
625 | gMC->Gspos("SFRM7",2, "SQM4", xposFr7, yposFr7, zposFr7,0, "ONLY"); | |
626 | ||
627 | ||
628 | //...................................................................................... | |
629 | //Frame - 8 | |
630 | ||
631 | ||
632 | Float_t frame8[3] ; | |
633 | frame8[0] = 2.5/2.0; | |
634 | frame8[1] = 1.2/2.0; | |
635 | frame8[2] = frame1[2]; | |
636 | ||
637 | gMC->Gsvolu("SFRM8", "BOX", idPGF30, frame8, 3); //Frame - 2 | |
638 | ||
639 | Float_t rrib8[3]; | |
640 | rrib8[0] = frame8[0]-1.0/2; | |
641 | rrib8[1] = frame8[1]; | |
642 | rrib8[2] =(frame8[2]-0.95)/2.0; | |
643 | ||
644 | gMC->Gsvolu("SFRR8", "BOX", idRoha, rrib8, 3); | |
645 | ||
646 | Float_t xposrrib8 = -1.0/2; | |
647 | Float_t yposrrib8 = 0.0; | |
648 | Float_t zposrrib8 = frame8[2] - rrib8[2]; | |
649 | ||
650 | gMC->Gspos("SFRR8",1, "SFRM8", xposrrib8, yposrrib8, zposrrib8,0, "ONLY"); | |
651 | gMC->Gspos("SFRR8",2, "SFRM8", xposrrib8, yposrrib8, -zposrrib8,0, "ONLY"); | |
652 | ||
653 | ||
654 | ||
655 | Float_t xposFr8 = frame8[0] + 20.6; | |
656 | Float_t yposFr8 = -frame8[1]; | |
657 | Float_t zposFr8 = 0.0; | |
658 | ||
659 | gMC->Gspos("SFRM8",1, "SQM3", xposFr8, yposFr8, zposFr8,0, "ONLY"); | |
660 | gMC->Gspos("SFRM8",2, "SQM4", xposFr8, yposFr8, zposFr8,0, "ONLY"); | |
661 | ||
662 | ||
663 | ||
664 | ||
665 | //^^^^^^^^^^^^^^^^^^^^^^^^^ Sensitive volumes ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ | |
666 | ||
667 | Float_t zsenv = 0.5; // distance between two cathode plane | |
668 | ||
669 | //Segment-0 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ | |
670 | bparH[0] = 94.5/2.; | |
671 | bparH[1] = 1.2/2.; | |
672 | bparH[2] = zsenv/2.; | |
673 | gMC->Gsvolu("SC3G0", "BOX", idGas, bparH, 3); | |
674 | gMC->Gsvolu("SC4G0", "BOX", idGas, bparH, 3); | |
675 | ||
676 | //Segment-1 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ | |
677 | pgpar[0] = 0.; | |
678 | pgpar[1] = 90.; | |
679 | pgpar[2] = 5; | |
680 | pgpar[3] = 2.; | |
681 | pgpar[4] = -zsenv/2.; | |
682 | pgpar[5] = 23.1; | |
683 | pgpar[6] = 117.6; | |
684 | pgpar[7] = zsenv/2.; | |
685 | pgpar[8] = pgpar[5]; | |
686 | pgpar[9] = pgpar[6]; | |
687 | gMC->Gsvolu("SC3G1", "PGON", idGas, pgpar, 10); | |
688 | gMC->Gsvolu("SC4G1", "PGON", idGas, pgpar, 10); | |
689 | ||
690 | ||
691 | //Segment-2 - vertical box ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ | |
692 | ||
693 | bparV[0] = 1.0/2.; | |
694 | bparV[1] = 95.5/2.; | |
695 | bparV[2] = zsenv/2.; | |
696 | gMC->Gsvolu("SC3G2", "BOX", idGas, bparV, 3); | |
697 | gMC->Gsvolu("SC4G2", "BOX", idGas, bparV, 3); | |
698 | ||
699 | //........................................................................................... | |
700 | ||
701 | xposHorBox = bparH[0] + 23.1; | |
702 | yposHorBox = -bparH[1]; | |
703 | ||
704 | xposVerBox = -bparV[0]; | |
705 | yposVerBox = bparV[1] + 23.1; | |
706 | ||
707 | gMC->Gspos("SC3G0",1, "SQM3", xposHorBox,yposHorBox,0.,0, "ONLY"); | |
708 | gMC->Gspos("SC3G1",1, "SQM3", 0.,0.,0.,0, "ONLY"); | |
709 | gMC->Gspos("SC3G2",1, "SQM3", xposVerBox,yposVerBox,0.,0, "ONLY"); | |
710 | ||
711 | ||
712 | gMC->Gspos("SC4G0",1, "SQM4", xposHorBox,yposHorBox,0.,0, "ONLY"); | |
713 | gMC->Gspos("SC4G1",1, "SQM4", 0.,0.,0.,0, "ONLY"); | |
714 | gMC->Gspos("SC4G2",1, "SQM4", xposVerBox,yposVerBox,0.,0, "ONLY"); | |
715 | ||
716 | ||
717 | //^^^^^^^^^^^^^^^^^^^^^^^^^ Sensitive volumes ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ | |
718 | ||
719 | ||
720 | ||
721 | //################################################################################################## | |
722 | // Positioning Quadrant in chamber#3 and chamber#4 | |
723 | //################################################################################################## | |
724 | /******Transformations for Quadrant********************************************** | |
725 | || I => Quadrant I: no rotation | |
726 | || | |
727 | II. || I. II => Quadrant II: Reflaction of Quadrant I in XZ plane | |
728 | || => TGeoRotation("Qrot2",90.,180.,90.,90.,180.,0.); | |
729 | ============= | |
730 | || III => Quadrant III: 180 degree rotation of Quadrant I in XY plane | |
731 | III. || IV. => TGeoRotation("Qrot3",90.,180.,90.,90.,180.,0.); | |
732 | || IV => Quadrant IV:-180 degree rotation of Quadrant II in XY plane | |
733 | => TGeoRotation("Qrot4",90.,0.,90.,-90.,180.,0.); | |
734 | **********************************************************************************************/ | |
735 | ||
736 | Int_t detElemId1 = 1; // quadrant I | |
737 | Int_t detElemId2 = 0; // quadrant II | |
738 | Int_t detElemId3 = 3; // quadrant III | |
739 | Int_t detElemId4 = 2; // quadrant IV | |
740 | ||
741 | //Float_t halfChamber = zCbb + zCu + zRoha + zMeb + zEeb + zsenv/2; | |
742 | // cout<< "\n half_chamber \t" << half_chamber << endl; | |
743 | ||
744 | Float_t halfChamber = tpar1[2];// tpar1[2] = 6.8/2; | |
745 | ||
746 | // ------------------------------St2 Chamber3------------------------------------------------ | |
747 | ||
748 | // GetEnvelopes(2)->AddEnvelope("S3M0", 300, true,TGeoTranslation(0.,0.,0.)); | |
749 | GetEnvelopes(2)->AddEnvelope("SQM3", 300+detElemId1, 1, TGeoTranslation( 0., 0., - halfChamber)); | |
750 | GetEnvelopes(2)->AddEnvelope("SQM3", 300+detElemId2, 2, TGeoTranslation( 0., 0., + halfChamber), | |
751 | TGeoRotation("Qrot3",90.,180.,90.,90.,180.,0.)); | |
752 | GetEnvelopes(2)->AddEnvelope("SQM3", 300+detElemId3, 3, TGeoTranslation( 0., 0., - halfChamber), | |
753 | TGeoRotation("Qrot3",90.,180.,90.,270.,0.,0.)); | |
754 | GetEnvelopes(2)->AddEnvelope("SQM3", 300+detElemId4, 4, TGeoTranslation( 0., 0., + halfChamber), | |
755 | TGeoRotation("Qrot4",90.,0.,90.,-90.,180.,0.)); | |
756 | ||
757 | //--------------------------------St2 Chamber4------------------------------------------------- | |
758 | ||
759 | GetEnvelopes(3)->AddEnvelope("SQM4", 400+detElemId1, 1, TGeoTranslation( 0., 0., - halfChamber)); | |
760 | GetEnvelopes(3)->AddEnvelope("SQM4", 400+detElemId2, 2, TGeoTranslation( 0., 0., + halfChamber), | |
761 | TGeoRotation("Qrot2",90.,180.,90.,90.,180.,0.)); | |
762 | GetEnvelopes(3)->AddEnvelope("SQM4", 400+detElemId3, 3, TGeoTranslation( 0., 0., - halfChamber), | |
763 | TGeoRotation("Qrot3",90.,180.,90.,270.,0.,0.)); | |
764 | GetEnvelopes(3)->AddEnvelope("SQM4", 400+detElemId4, 4, TGeoTranslation( 0., 0., + halfChamber), | |
765 | TGeoRotation("Qrot4",90.,0.,90.,-90.,180.,0.)); | |
766 | ||
767 | //********************************************************************************************** | |
768 | ||
769 | } | |
770 | //______________________________________________________________________________ | |
771 | void AliMUONSt2GeometryBuilderV2::SetVolumes() | |
772 | { | |
773 | /// Defines the volumes for the station2 chambers. | |
774 | ||
775 | if (gAlice->GetModule("SHIL")) { | |
776 | SetMotherVolume(2, "YOUT1"); | |
777 | SetMotherVolume(3, "YOUT1"); | |
778 | } | |
779 | ||
780 | // Define chamber volumes as virtual | |
781 | SetVolume(2, "SC03", true); | |
782 | SetVolume(3, "SC04", true); | |
783 | } | |
784 | ||
785 | //______________________________________________________________________________ | |
786 | void AliMUONSt2GeometryBuilderV2::SetTransformations() | |
787 | { | |
788 | /// Defines the transformations for the station2 chambers. | |
789 | ||
790 | Double_t zpos1 = - AliMUONConstants::DefaultChamberZ(2); | |
791 | SetTranslation(2, TGeoTranslation(0., 0., zpos1)); | |
792 | ||
793 | Double_t zpos2 = - AliMUONConstants::DefaultChamberZ(3); | |
794 | SetTranslation(3, TGeoTranslation(0., 0., zpos2)); | |
795 | } | |
796 | ||
797 | //______________________________________________________________________________ | |
798 | void AliMUONSt2GeometryBuilderV2::SetSensitiveVolumes() | |
799 | { | |
800 | /// Defines the sensitive volumes for station2 chambers. | |
801 | ||
802 | GetGeometry(2)->SetSensitiveVolume("SC3G0"); | |
803 | GetGeometry(2)->SetSensitiveVolume("SC3G1"); | |
804 | GetGeometry(2)->SetSensitiveVolume("SC3G2"); | |
805 | ||
806 | GetGeometry(3)->SetSensitiveVolume("SC4G0"); | |
807 | GetGeometry(3)->SetSensitiveVolume("SC4G1"); | |
808 | GetGeometry(3)->SetSensitiveVolume("SC4G2"); | |
809 | } | |
810 |