New code from Piergiorgio added
[u/mrichter/AliRoot.git] / FMD / AliFMDv1.cxx
CommitLineData
4c039060 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/*
17$Log$
0d630091 18Revision 1.7 1999/09/29 09:24:14 fca
19Introduction of the Copyright and cvs Log
20
4c039060 21*/
22
fe4da5cc 23///////////////////////////////////////////////////////////////////////////////
24// //
0d630091 25// Forward Multiplicity Detector version 1 //
fe4da5cc 26// //
27//Begin_Html
28/*
1439f98e 29<img src="picts/AliFMDv1Class.gif">
fe4da5cc 30</pre>
31<br clear=left>
32<font size=+2 color=red>
33<p>The responsible person for this module is
34<a href="mailto:Valeri.Kondratiev@cern.ch">Valeri Kondratiev</a>.
35</font>
36<pre>
37*/
38//End_Html
39// //
40///////////////////////////////////////////////////////////////////////////////
41
42#include "AliRun.h"
43#include "AliFMDv1.h"
44#include "AliMC.h"
0d630091 45#include "stdlib.h"
fe4da5cc 46
47ClassImp(AliFMDv1)
48
49//_____________________________________________________________________________
99560b47 50AliFMDv1::AliFMDv1()
fe4da5cc 51{
52 //
0d630091 53 // Default constructor for FMD version 1
fe4da5cc 54 //
55}
56
57//_____________________________________________________________________________
58AliFMDv1::AliFMDv1(const char *name, const char *title)
59 : AliFMD(name,title)
60{
61 //
62 // Standard constructor for FMD version 1
63 //
99560b47 64 AliModule *start = gAlice->GetModule("START");
65 if(start) {
66 Error("ctor","This version of FMD is incompatible with START\n");
67 exit(1);
68 }
fe4da5cc 69}
70
0d630091 71//___________________________________________
fe4da5cc 72void AliFMDv1::CreateGeometry()
73{
74 //
75 // Creation of the geometry of the FMD version 1
76 //
77 //Begin_Html
78 /*
1439f98e 79 <img src="picts/AliFMDv1Tree.gif">
fe4da5cc 80 */
81 //End_Html
82 //Begin_Html
83 /*
1439f98e 84 <img src="picts/AliFMDv1.gif">
fe4da5cc 85 */
86 //End_Html
87
fe4da5cc 88
0d630091 89 Float_t rout, z;
fe4da5cc 90 Float_t par[3], rin;
0d630091 91
92 Int_t irotm[100];
fe4da5cc 93
0d630091 94 TArrayI &idtmed = *fIdtmed;
fe4da5cc 95
96 // ********************************************************
0d630091 97 // DEFINE DISK#3 OF FMD
fe4da5cc 98 // ********************************************************
99
0d630091 100 gMC->Gsvolu("BR3_", "BOX ", idtmed[4], par, 0);
101 gMC->Gsvolu("CB3_", "BOX ", idtmed[5], par, 0);
102 gMC->Gsvolu("BR5_", "BOX ", idtmed[4], par, 0);
103 gMC->Gsvolu("CB5_", "BOX ", idtmed[5], par, 0);
104
105 // Define parameters for the right disk
fe4da5cc 106
107 rin = 4.5;
108 rout = 10.5;
0d630091 109 z = 77.;
110
111 // Ring #1
112
fe4da5cc 113 par[0] = rin;
114 par[1] = rout;
0d630091 115 par[2] = 1.5;
116 gMC->Gsvolu("R1R3", "TUBE", idtmed[1], par, 3);
117 gMC->Gspos("R1R3", 1, "ALIC", 0., 0., z + 1.5, 0, "ONLY");
118
119 // Ring #2
120
fe4da5cc 121 par[0] = rout;
0d630091 122 par[1] = rout + .65;
123 par[2] = 1.5;
124 gMC->Gsvolu("R2R3", "TUBE", idtmed[2], par, 3);
125 gMC->Gspos("R2R3", 1, "ALIC", 0., 0., z + 1.5, 0, "ONLY");
126
127 // Ring #3
128
129 par[0] = rout + .65;
130 par[1] = rout + 5.65;
fe4da5cc 131 par[2] = .025;
0d630091 132 gMC->Gsvolu("R3R3", "TUBE", idtmed[3], par, 3);
133 gMC->Gspos("R3R3", 1, "ALIC", 0., 0., z + 1.525, 0, "ONLY");
134
135 // Bracket #1
136
137 par[0] = 1.5;
138 par[1] = 0.1;
139 par[2] = 15.35;
140 gMC->Matrix(irotm[11], 90., 0., 161.2, 90., 71.2, 90.);
141 gMC->Gsposp("BR3_",1,"ALIC", 0., 25.25, 85.0, irotm[11],"ONLY",par,3);
142 par[1] = 0.5;
143 gMC->Gsposp("CB3_",1,"ALIC", 0., 25.41, 84.53,irotm[11],"ONLY",par,3);
144
145
146 // Bracket #2
147
148 par[0] = 1.5;
149 par[1] = 0.1;
150 par[2] = 15.35;
151 gMC->Matrix(irotm[12], 90., 180., 161.2, 270., 71.2, 270.);
152 gMC->Gsposp("BR3_",2,"ALIC", 0.,-25.25, 85.0, irotm[12],"ONLY",par,3);
153 par[1] = 0.5;
154 gMC->Gsposp("CB3_",2,"ALIC", 0.,-25.41, 84.53,irotm[12],"ONLY",par,3);
155
156
157 // Bracket #3
158
159 par[0] = 1.5;
160 par[1] = 0.1;
161 par[2] = 15.35;
162 gMC->Matrix(irotm[13], 90., 270., 161.2, 0., 71.2, 0.);
163 gMC->Gsposp("BR3_",3,"ALIC", 25.25, 0., 85.0, irotm[13],"ONLY",par,3);
164 par[1] = 0.5;
165 gMC->Gsposp("CB3_",3,"ALIC", 25.41, 0., 84.53,irotm[13],"ONLY",par,3);
166
167
168 // Bracket #4
169
170 par[0] = 1.5;
171 par[1] = 0.1;
172 par[2] = 15.35;
173 gMC->Matrix(irotm[14], 90., 90., 161.2, 180., 71.2, 180.);
174 gMC->Gsposp("BR3_",4,"ALIC", -25.25, 0., 85.0, irotm[14],"ONLY",par,3);
175 par[1] = 0.5;
176 gMC->Gsposp("CB3_",4,"ALIC", -25.41, 0., 84.53,irotm[14],"ONLY",par,3);
177
178 // Right support ring
179
180 par[0] = 39.;
181 par[1] = 41.;
182 par[2] = .5;
183 gMC->Gsvolu("R1SP", "TUBE", idtmed[4], par, 3);
184 gMC->Gspos("R1SP", 1, "ALIC", 0., 0., 89.5, 0, "ONLY");
185
186 // Define parameters for the left disk
fe4da5cc 187
188 rin = 4.5;
189 rout = 10.5;
0d630091 190 z = -77.;
191
192 // Ring #1
193
fe4da5cc 194 par[0] = rin;
195 par[1] = rout;
0d630091 196 par[2] = 1.5;
197 gMC->Gsvolu("R1L3", "TUBE", idtmed[1], par, 3);
198 gMC->Gspos("R1L3", 1, "ALIC", 0., 0., z - 1.5, 0, "ONLY");
199
200 // Ring #2
201
fe4da5cc 202 par[0] = rout;
0d630091 203 par[1] = rout + .65;
204 par[2] = 1.5;
205 gMC->Gsvolu("R2L3", "TUBE", idtmed[2], par, 3);
206 gMC->Gspos("R2L3", 1, "ALIC", 0., 0., z - 1.5, 0, "ONLY");
207
208 // Ring #3
209
210 par[0] = rout + .65;
211 par[1] = rout + 5.65;
fe4da5cc 212 par[2] = .025;
0d630091 213 gMC->Gsvolu("R3L3", "TUBE", idtmed[3], par, 3);
214 gMC->Gspos("R3L3", 1, "ALIC", 0., 0., z - 1.525, 0, "ONLY");
215
216 // Bracket #1
217
218 par[0] = 1.5;
219 par[1] = 0.1;
220 par[2] = 15.35;
221 gMC->Matrix(irotm[61], 90., 0., 180.-161.2, 90., 180.-71.2, 90.);
222 gMC->Gsposp("BR3_",5,"ALIC", 0., 25.25, -85.0, irotm[61],"ONLY",par,3);
223 par[1] = 0.5;
224 gMC->Gsposp("CB3_",5,"ALIC", 0., 25.41, -84.53,irotm[61],"ONLY",par,3);
225
226 gMC->Matrix(irotm[91], 90., 0., 90.-10.6, 90., 180.-10.6, 90.);
227 par[1] = 0.1;
228 par[2] = 81.4;
229 gMC->Gsposp("BR5_",1,"ALIC", 0., 55.0, -170.0,irotm[91],"ONLY",par,3);
230 par[1] = 0.5;
231 gMC->Gsposp("CB5_",1,"ALIC", 0., 55.5, -169.5,irotm[91],"ONLY",par,3);
232
233 // Bracket #2
234
235 par[0] = 1.5;
236 par[1] = 0.1;
237 par[2] = 15.35;
238 gMC->Matrix(irotm[62], 90., 180., 180.-161.2, 270., 180.-71.2, 270.);
239 gMC->Gsposp("BR3_",6,"ALIC", 0.,-25.25, -85.0,irotm[62],"ONLY",par,3);
240 par[1] = 0.5;
241 gMC->Gsposp("CB3_",6,"ALIC", 0.,-25.41, -84.53,irotm[62],"ONLY",par,3);
242
243 gMC->Matrix(irotm[92], 90., 180., 90.-10.6, 270., 180.-10.6, 270.);
244 par[1] = 0.1;
245 par[2] = 81.4;
246 gMC->Gsposp("BR5_",2,"ALIC", 0., -55.0, -170.0,irotm[92],"ONLY",par,3);
247 par[1] = 0.5;
248 gMC->Gsposp("CB5_",2,"ALIC", 0., -55.5, -169.5,irotm[92],"ONLY",par,3);
249
250 // Bracket #3
251
252 par[0] = 1.5;
253 par[1] = 0.1;
254 par[2] = 15.35;
255 gMC->Matrix(irotm[63], 90., 270., 180.-161.2, 0., 180.-71.2, 0.);
256 gMC->Gsposp("BR3_",7,"ALIC", 25.25, 0., -85.0, irotm[63],"ONLY",par,3);
257 par[1] = 0.5;
258 gMC->Gsposp("CB3_",7,"ALIC", 25.41, 0., -84.53,irotm[63],"ONLY",par,3);
259
260 gMC->Matrix(irotm[93], 90., 270., 90.-10.6, 0., 180.-10.6, 0.);
261 par[1] = 0.1;
262 par[2] = 81.4;
263 gMC->Gsposp("BR5_",3,"ALIC", 55., 0., -170.,irotm[93],"ONLY",par,3);
264 par[1] = 0.5;
265 gMC->Gsposp("CB5_",3,"ALIC", 55.5, 0., -169.5,irotm[93],"ONLY",par,3);
266
267 // Bracket #4
268
269 par[0] = 1.5;
270 par[1] = 0.1;
271 par[2] = 15.35;
272 gMC->Matrix(irotm[64], 90., 90., 180.-161.2, 180., 180.-71.2, 180.);
273 gMC->Gsposp("BR3_",8,"ALIC", -25.25, 0., -85., irotm[64],"ONLY",par,3);
274 par[1] = 0.5;
275 gMC->Gsposp("CB3_",8,"ALIC", -25.41, 0., -84.53,irotm[64],"ONLY",par,3);
276
277 gMC->Matrix(irotm[94], 90., 90., 90.-10.6, 180., 180.-10.6, 180.);
278 par[1] = 0.1;
279 par[2] = 81.4;
280 gMC->Gsposp("BR5_",4,"ALIC", -55., 0., -170.,irotm[94],"ONLY",par,3);
281 par[1] = 0.5;
282 gMC->Gsposp("CB5_",4,"ALIC", -55.5, 0., -169.5,irotm[94],"ONLY",par,3);
283
284 // Central support ring
285
286 par[0] = 39.;
287 par[1] = 41.;
288 par[2] = .5;
289 gMC->Gsvolu("R2SP", "TUBE", idtmed[4], par, 3);
290 gMC->Gspos("R2SP", 1, "ALIC", 0., 0., -89.5, 0, "ONLY");
291
292 // Left support ring
293
294 par[0] = 69.;
295 par[1] = 71.;
296 par[2] = .5;
297 gMC->Gsvolu("R3SP", "TUBE", idtmed[4], par, 3);
298 gMC->Gspos("R3SP", 1, "ALIC", 0., 0., -249.5, 0, "ONLY");
299
fe4da5cc 300 // ********************************************************
0d630091 301 // DEFINE DISK#2 OF FMD
fe4da5cc 302 // ********************************************************
303
0d630091 304 gMC->Gsvolu("BR2_", "BOX ", idtmed[4], par, 0);
305 gMC->Gsvolu("CB2_", "BOX ", idtmed[5], par, 0);
306
307 // Define parameters for right disk #2
fe4da5cc 308
0d630091 309 rin = 7.7;
310 rout = 13.7;
311 z = 64.4;
312
313 // Ring #1
314
fe4da5cc 315 par[0] = rin;
316 par[1] = rout;
0d630091 317 par[2] = 1.5;
318 gMC->Gsvolu("R1R2", "TUBE", idtmed[1], par, 3);
319 gMC->Gspos("R1R2", 1, "ALIC", 0., 0., z + 1.5, 0, "ONLY");
320
321 // Ring #2
322
fe4da5cc 323 par[0] = rout;
0d630091 324 par[1] = rout + .65;
325 par[2] = 1.5;
326 gMC->Gsvolu("R2R2", "TUBE", idtmed[2], par, 3);
327 gMC->Gspos("R2R2", 1, "ALIC", 0., 0., z + 1.5, 0, "ONLY");
328
329 // Ring #3
330
331 par[0] = rout + .65;
332 par[1] = rout + 5.65;
fe4da5cc 333 par[2] = .025;
0d630091 334 gMC->Gsvolu("R3R2", "TUBE", idtmed[3], par, 3);
335 gMC->Gspos("R3R2", 1, "ALIC", 0., 0., z + 1.525, 0, "ONLY");
336
337 // Bracket #1
338
339 par[0] = 1.5;
340 par[1] = 0.1;
341 par[2] = 17.35;
342 gMC->Matrix(irotm[21], 90., 30., 139.3, 120., 49.3, 120.);
343 gMC->Gsposp("BR2_",1,"ALIC", -13.4, 23.3, 78.7, irotm[21],"ONLY",par,3);
344 par[1] = 0.5;
345 gMC->Gsposp("CB2_",1,"ALIC", -13.56, 23.6, 78.4, irotm[21],"ONLY",par,3);
346
347 // Bracket #2
348
349 par[0] = 1.5;
350 par[1] = 0.1;
351 par[2] = 17.35;
352 gMC->Matrix(irotm[22], 90., 210., 139.3, 300., 49.3, 300.);
353 gMC->Gsposp("BR2_",2,"ALIC", 13.4,-23.3, 78.7, irotm[22],"ONLY",par,3);
354 par[1] = 0.5;
355 gMC->Gsposp("CB2_",2,"ALIC", 13.5,-23.6, 78.4,irotm[22],"ONLY",par,3);
356
357
358 // Bracket #3
359
360 par[0] = 1.5;
361 par[1] = 0.1;
362 par[2] = 17.35;
363 gMC->Matrix(irotm[23], 90., 300., 139.3, 30., 49.3, 30.);
364 gMC->Gsposp("BR2_",3,"ALIC", 23.3, 13.4, 78.7, irotm[23],"ONLY",par,3);
365 par[1] = 0.5;
366 gMC->Gsposp("CB2_",3,"ALIC", 23.6, 13.56, 78.4,irotm[23],"ONLY",par,3);
367
368 // Bracket #4
369
370 par[0] = 1.5;
371 par[1] = 0.1;
372 par[2] = 17.35;
373 gMC->Matrix(irotm[24], 90., 120., 139.3, 210., 49.3, 210.);
374 gMC->Gsposp("BR2_",4,"ALIC", -23.3, -13.4, 78.7,irotm[24],"ONLY",par,3);
375 par[1] = 0.5;
376 gMC->Gsposp("CB2_",4,"ALIC", -23.6, -13.56, 78.4,irotm[24],"ONLY",par,3);
377
378 // Define parameters for left disk
fe4da5cc 379
0d630091 380 rin = 7.7;
381 rout = 13.7;
382 z = -64.4;
383
384 // Ring #1
385
fe4da5cc 386 par[0] = rin;
387 par[1] = rout;
0d630091 388 par[2] = 1.5;
389 gMC->Gsvolu("R1L2", "TUBE", idtmed[1], par, 3);
390 gMC->Gspos("R1L2", 1, "ALIC", 0., 0., z - 1.5, 0, "ONLY");
391
392 // Ring #2
393
fe4da5cc 394 par[0] = rout;
0d630091 395 par[1] = rout + .65;
396 par[2] = 1.5;
397 gMC->Gsvolu("R2L2", "TUBE", idtmed[2], par, 3);
398 gMC->Gspos("R2L2", 1, "ALIC", 0., 0., z - 1.5, 0, "ONLY");
399
400 // Ring #3
401
402 par[0] = rout + .65;
403 par[1] = rout + 5.65;
fe4da5cc 404 par[2] = .025;
0d630091 405 gMC->Gsvolu("R3L2", "TUBE", idtmed[3], par, 3);
406 gMC->Gspos("R3L2", 1, "ALIC", 0., 0., z - 1.525, 0, "ONLY");
407
408 // Bracket #1
409
410 par[0] = 1.5;
411 par[1] = 0.1;
412 par[2] = 17.35;
413 gMC->Matrix(irotm[51], 90., 30., 180.-139.3, 120., 180.-49.3, 120.);
414 gMC->Gsposp("BR2_",5,"ALIC", -13.4, 23.3, -78.7, irotm[51],"ONLY",par,3);
415 par[1] = 0.5;
416 gMC->Gsposp("CB2_",5,"ALIC", -13.56, 23.6, -78.4,irotm[51],"ONLY",par,3);
417
418 gMC->Matrix(irotm[81], 90., 30., 90.-10.6, 120., 180.-10.6, 120.);
419 par[1] = 0.1;
420 par[2] = 81.4;
421 gMC->Gsposp("BR5_",5,"ALIC", -27.5, 47.6, -170.0,irotm[81],"ONLY",par,3);
422 par[1] = 0.5;
423 gMC->Gsposp("CB5_",5,"ALIC", -27.85, 48., -169.5,irotm[81],"ONLY",par,3);
424
425 // Bracket #2
426
427 par[0] = 1.5;
428 par[1] = 0.1;
429 par[2] = 17.35;
430 gMC->Matrix(irotm[52], 90., 210., 180.-139.3, 300., 180.-49.3, 300.);
431 gMC->Gsposp("BR2_",6,"ALIC", 13.4, -23.3, -78.7,irotm[52],"ONLY",par,3);
432 par[1] = 0.5;
433 gMC->Gsposp("CB2_",6,"ALIC", 13.56, -23.6, -78.4,irotm[52],"ONLY",par,3);
434
435 gMC->Matrix(irotm[82], 90., 210., 90.-10.6, 300., 180.-10.6, 300.);
436 par[1] = 0.1;
437 par[2] = 81.4;
438 gMC->Gsposp("BR5_",6,"ALIC", 27.5, -47.6, -170.0,irotm[82],"ONLY",par,3);
439 par[1] = 0.5;
440 gMC->Gsposp("CB5_",6,"ALIC", 27.85, -48., -169.5,irotm[82],"ONLY",par,3);
441
442 // Bracket #3
443
444 par[0] = 1.5;
445 par[1] = 0.1;
446 par[2] = 17.35;
447 gMC->Matrix(irotm[53], 90., 300., 180.-139.3, 30., 180.-49.3, 30.);
448 gMC->Gsposp("BR2_",7,"ALIC", 23.3, 13.4, -78.7, irotm[53],"ONLY",par,3);
449 par[1] = 0.5;
450 gMC->Gsposp("CB2_",7,"ALIC", 23.6, 13.56, -78.4,irotm[53],"ONLY",par,3);
451
452 gMC->Matrix(irotm[83], 90., 300., 90.-10.6, 30., 180.-10.6, 30.);
453 par[1] = 0.1;
454 par[2] = 81.4;
455 gMC->Gsposp("BR5_",7,"ALIC", 47.6, 27.5, -170.,irotm[83],"ONLY",par,3);
456 par[1] = 0.5;
457 gMC->Gsposp("CB5_",7,"ALIC", 48., 27.85, -169.5,irotm[83],"ONLY",par,3);
458
459 // Bracket #4
460
461 par[0] = 1.5;
462 par[1] = 0.1;
463 par[2] = 17.35;
464 gMC->Matrix(irotm[54], 90., 120., 180.-139.3, 210., 180.-49.3, 210.);
465 gMC->Gsposp("BR2_",8,"ALIC", -23.3, -13.4, -78.7, irotm[54],"ONLY",par,3);
466 par[1] = 0.5;
467 gMC->Gsposp("CB2_",8,"ALIC", -23.6, -13.56, -78.4,irotm[54],"ONLY",par,3);
468
469 gMC->Matrix(irotm[84], 90., 120., 90.-10.6, 210., 180.-10.6, 210.);
470 par[1] = 0.1;
471 par[2] = 81.4;
472 gMC->Gsposp("BR5_",8,"ALIC", -47.6, -27.5, -170.,irotm[84],"ONLY",par,3);
473 par[1] = 0.5;
474 gMC->Gsposp("CB5_",8,"ALIC", -48., -27.85, -169.5,irotm[84],"ONLY",par,3);
475
fe4da5cc 476 // ********************************************************
0d630091 477 // DEFINE DISK#1 OF FMD
fe4da5cc 478 // ********************************************************
479
0d630091 480 gMC->Gsvolu("BR1_", "BOX ", idtmed[4], par, 0);
481 gMC->Gsvolu("CB1_", "BOX ", idtmed[5], par, 0);
fe4da5cc 482
0d630091 483 // Define parameters for right disk #1
484
485 rin = 12.4;
486 rout = 18.4;
487 z = 59.4;
488
489 // Ring #1
490
fe4da5cc 491 par[0] = rin;
492 par[1] = rout;
0d630091 493 par[2] = 1.5;
494 gMC->Gsvolu("R1R1", "TUBE", idtmed[1], par, 3);
495 gMC->Gspos("R1R1", 1, "ALIC", 0., 0., z + 1.5, 0, "ONLY");
496
497 // Ring #2
498
fe4da5cc 499 par[0] = rout;
0d630091 500 par[1] = rout + .65;
501 par[2] = 1.5;
502 gMC->Gsvolu("R2R1", "TUBE", idtmed[2], par, 3);
503 gMC->Gspos("R2R1", 1, "ALIC", 0., 0., z + 1.5, 0, "ONLY");
504
505 // Ring #3
506
507 par[0] = rout + .65;
508 par[1] = rout + 5.65;
fe4da5cc 509 par[2] = .025;
0d630091 510 gMC->Gsvolu("R3R1", "TUBE", idtmed[3], par, 3);
511 gMC->Gspos("R3R1", 1, "ALIC", 0., 0., z + 1.525, 0, "ONLY");
512
513 // Bracket #1
514
515 par[0] = 1.5;
516 par[1] = 0.1;
517 par[2] = 17.5;
518 gMC->Matrix(irotm[31], 90., 60., 128., 150., 38., 150.);
519 gMC->Gsposp("BR1_",1,"ALIC", -25.3, 14.6, 76.2, irotm[31],"ONLY",par,3);
520 par[1] = 0.5;
521 gMC->Gsposp("CB1_",1,"ALIC", -25.35, 14.8, 75.9, irotm[31],"ONLY",par,3);
522
523 // Bracket #2
524
525 par[0] = 1.5;
526 par[1] = 0.1;
527 par[2] = 17.5;
528 gMC->Matrix(irotm[32], 90., 240., 128., 330., 38., 330.);
529 gMC->Gsposp("BR1_",2,"ALIC", 25.3, -14.6, 76.2, irotm[32],"ONLY",par,3);
530 par[1] = 0.5;
531 gMC->Gsposp("CB1_",2,"ALIC", 25.35, -14.8, 75.9,irotm[32],"ONLY",par,3);
532
533 // Bracket #3
534
535 par[0] = 1.5;
536 par[1] = 0.1;
537 par[2] = 17.5;
538 gMC->Matrix(irotm[33], 90., 330., 128., 60., 38., 60.);
539 gMC->Gsposp("BR1_",3,"ALIC", 14.6, 25.3, 76.2, irotm[33],"ONLY",par,3);
540 par[1] = 0.5;
541 gMC->Gsposp("CB1_",3,"ALIC", 14.8, 25.35, 75.9,irotm[33],"ONLY",par,3);
542
543 // Bracket #4
544
545 par[0] = 1.5;
546 par[1] = 0.1;
547 par[2] = 17.5;
548 gMC->Matrix(irotm[34], 90., 150., 128., 240., 38., 240.);
549 gMC->Gsposp("BR1_",4,"ALIC", -14.6, -25.3, 76.2,irotm[34],"ONLY",par,3);
550 par[1] = 0.5;
551 gMC->Gsposp("CB1_",4,"ALIC", -14.8, -25.35, 75.9,irotm[34],"ONLY",par,3);
552
553 // Define parameters for left disk #1
fe4da5cc 554
0d630091 555 rin = 12.4;
556 rout = 18.4;
557 z = -59.4;
558
559 // Ring #1
560
fe4da5cc 561 par[0] = rin;
562 par[1] = rout;
0d630091 563 par[2] = 1.5;
564 gMC->Gsvolu("R1L1", "TUBE", idtmed[1], par, 3);
565 gMC->Gspos("R1L1", 1, "ALIC", 0., 0., z - 1.5, 0, "ONLY");
566
567 // Ring #2
568
fe4da5cc 569 par[0] = rout;
0d630091 570 par[1] = rout + .65;
571 par[2] = 1.5;
572 gMC->Gsvolu("R2L1", "TUBE", idtmed[2], par, 3);
573 gMC->Gspos("R2L1", 1, "ALIC", 0., 0., z - 1.5, 0, "ONLY");
574
575 // Ring #3
576
577 par[0] = rout + .65;
578 par[1] = rout + 5.65;
fe4da5cc 579 par[2] = .025;
0d630091 580 gMC->Gsvolu("R3L1", "TUBE", idtmed[3], par, 3);
581 gMC->Gspos("R3L1", 1, "ALIC", 0., 0., z - 1.525, 0, "ONLY");
582
583 // Bracket #1
584
585 par[0] = 1.5;
586 par[1] = 0.1;
587 par[2] = 17.5;
588 gMC->Matrix(irotm[41], 90., 60., 180.-128., 150., 180.-38.0, 150.);
589 gMC->Gsposp("BR1_",5,"ALIC", -25.3, 14.6, -76.2, irotm[41],"ONLY",par,3);
590 par[1] = 0.5;
591 gMC->Gsposp("CB1_",5,"ALIC", -25.35, 14.8, -75.9,irotm[41],"ONLY",par,3);
592
593 gMC->Matrix(irotm[1], 90., 60., 90.-10.6, 150., 180.-10.6, 150.);
594 par[1] = 0.1;
595 par[2] = 81.4;
596 gMC->Gsposp("BR5_",9,"ALIC", -47.6, 27.5, -170.,irotm[1],"ONLY",par,3);
597 par[1] = 0.5;
598 gMC->Gsposp("CB5_",9,"ALIC", -48., 27.85, -169.5,irotm[1],"ONLY",par,3);
599
600 // Bracket #2
601
602 par[0] = 1.5;
603 par[1] = 0.1;
604 par[2] = 17.5;
605 gMC->Matrix(irotm[42], 90., 240., 180.-128., 330., 180.-38., 330.);
606 gMC->Gsposp("BR1_",6,"ALIC", 25.3,-14.6, -76.2,irotm[42],"ONLY",par,3);
607 par[1] = 0.5;
608 gMC->Gsposp("CB1_",6,"ALIC", 25.35,-14.8, -75.9,irotm[42],"ONLY",par,3);
609
610 gMC->Matrix(irotm[2], 90., 240., 90.-10.6, 330., 180.-10.6, 330.);
611 par[1] = 0.1;
612 par[2] = 81.4;
613 gMC->Gsposp("BR5_",10,"ALIC", 47.6, -27.5, -170.0,irotm[2],"ONLY",par,3);
614 par[1] = 0.5;
615 gMC->Gsposp("CB5_",10,"ALIC", 48., -27.85, -169.5,irotm[2],"ONLY",par,3);
616
617 // Bracket #3
618
619 par[0] = 1.5;
620 par[1] = 0.1;
621 par[2] = 17.5;
622 gMC->Matrix(irotm[43], 90., 330., 180.-128., 60., 180.-38., 60.);
623 gMC->Gsposp("BR1_",7,"ALIC", 14.6, 25.3, -76.2, irotm[43],"ONLY",par,3);
624 par[1] = 0.5;
625 gMC->Gsposp("CB1_",7,"ALIC", 14.8, 25.35, -75.9,irotm[43],"ONLY",par,3);
626
627 gMC->Matrix(irotm[3], 90., 330., 90.-10.6, 60., 180.-10.6, 60.);
628 par[1] = 0.1;
629 par[2] = 81.4;
630 gMC->Gsposp("BR5_",11,"ALIC", 27.5, 47.6, -170.,irotm[3],"ONLY",par,3);
631 par[1] = 0.5;
632 gMC->Gsposp("CB5_",11,"ALIC", 27.85, 48., -169.5,irotm[3],"ONLY",par,3);
633
634 // Bracket #4
635
636 par[0] = 1.5;
637 par[1] = 0.1;
638 par[2] = 17.5;
639 gMC->Matrix(irotm[44], 90., 150., 180.-128., 240., 180.-38., 240.);
640 gMC->Gsposp("BR1_",8,"ALIC", -14.6, -25.3, -76.2, irotm[44],"ONLY",par,3);
641 par[1] = 0.5;
642 gMC->Gsposp("CB1_",8,"ALIC", -14.8, -25.35, -75.9,irotm[44],"ONLY",par,3);
643
644 gMC->Matrix(irotm[4], 90., 150., 90.-10.6, 240., 180.-10.6, 240.);
645 par[1] = 0.1;
646 par[2] = 81.4;
647 gMC->Gsposp("BR5_",12,"ALIC", -27.5, -47.6, -170.,irotm[4],"ONLY",par,3);
648 par[1] = 0.5;
649 gMC->Gsposp("CB5_",12,"ALIC", -27.85, -48., -169.5,irotm[4],"ONLY",par,3);
650
fe4da5cc 651 // ***********************************************************
0d630091 652 // DEFINE LEFT DISK#4 OF FMD
fe4da5cc 653 // ***********************************************************
0d630091 654
655 gMC->Gsvolu("BR4_", "BOX ", idtmed[4], par, 0);
656 gMC->Gsvolu("CB4_", "BOX ", idtmed[5], par, 0);
657
658 // Define parameters
659
660 rin = 4.5;
661 rout = 10.5;
fe4da5cc 662 z = -229.5;
0d630091 663
664 // Ring #1
665
fe4da5cc 666 par[0] = rin;
667 par[1] = rout;
0d630091 668 par[2] = 1.5;
669 gMC->Gsvolu("R1L4", "TUBE", idtmed[1], par, 3);
670 gMC->Gspos("R1L4", 1, "ALIC", 0., 0., z - 1.5, 0, "ONLY");
671
672 // Ring #2
673
fe4da5cc 674 par[0] = rout;
0d630091 675 par[1] = rout + .65;
676 par[2] = 1.5;
677 gMC->Gsvolu("R2L4", "TUBE", idtmed[2], par, 3);
678 gMC->Gspos("R2L4", 1, "ALIC", 0., 0., z - 1.5, 0, "ONLY");
679
680 // Ring #3
681
682 par[0] = rout + .65;
683 par[1] = rout + 5.65;
fe4da5cc 684 par[2] = .025;
0d630091 685 gMC->Gsvolu("R3L4", "TUBE", idtmed[3], par, 3);
686 gMC->Gspos("R3L4", 1, "ALIC", 0., 0., z - 1.525, 0, "ONLY");
687
688 // Bracket #1
689
690 par[0] = 1.5;
691 par[1] = 0.1;
692 par[2] = 31.25;
693 gMC->Matrix(irotm[71], 90., 0., 90.-71., 90., 180.-71., 90.);
694 gMC->Gsposp("BR4_",1,"ALIC", 0., 40.25, -240., irotm[71],"ONLY",par,3);
695 par[1] = 0.5;
696 gMC->Gsposp("CB4_",1,"ALIC", 0., 40.75, -239.5, irotm[71],"ONLY",par,3);
697
698 // Bracket #2
699
700 par[0] = 1.5;
701 par[1] = 0.1;
702 par[2] = 31.25;
703 gMC->Matrix(irotm[72], 90., 180., 90.-71., 270., 180.-71., 270.);
704 gMC->Gsposp("BR4_",2,"ALIC", 0., -40.25, -240., irotm[72],"ONLY",par,3);
705 par[1] = 0.5;
706 gMC->Gsposp("CB4_",2,"ALIC", 0., -40.75, -239.5,irotm[72],"ONLY",par,3);
707
708 // Bracket #3
709
710 par[0] = 1.5;
711 par[1] = 0.1;
712 par[2] = 31.25;
713 gMC->Matrix(irotm[73], 90., 270., 90.-71., 0., 180.-71., 0.);
714 gMC->Gsposp("BR4_",3,"ALIC", 40.25, 0., -240., irotm[73],"ONLY",par,3);
715 par[1] = 0.5;
716 gMC->Gsposp("CB4_",3,"ALIC", 40.75, 0., -239.5,irotm[73],"ONLY",par,3);
717
718 // Bracket #4
719
720 par[0] = 1.5;
721 par[1] = 0.1;
722 par[2] = 31.25;
723 gMC->Matrix(irotm[74], 90., 90., 90.-71., 180., 180.-71., 180.);
724 gMC->Gsposp("BR4_",4,"ALIC", -40.25, 0., -240.,irotm[74],"ONLY",par,3);
725 par[1] = 0.5;
726 gMC->Gsposp("CB4_",4,"ALIC", -40.75, 0., -239.5,irotm[74],"ONLY",par,3);
fe4da5cc 727}
0d630091 728
fe4da5cc 729//_____________________________________________________________________________
b316044f 730void AliFMDv1::DrawModule()
fe4da5cc 731{
732 //
733 // Draw a shaded view of the FMD version 1
734 //
735
fe4da5cc 736 // Set everything unseen
cfce8870 737 gMC->Gsatt("*", "seen", -1);
fe4da5cc 738 //
739 // Set ALIC mother transparent
cfce8870 740 gMC->Gsatt("ALIC","SEEN",0);
fe4da5cc 741 //
742 // Set the volumes visible
0d630091 743 gMC->Gsatt("R1R3","SEEN",1);
744 gMC->Gsatt("R2R3","SEEN",1);
745 gMC->Gsatt("R3R3","SEEN",1);
746 gMC->Gsatt("R1L3","SEEN",1);
747 gMC->Gsatt("R2L3","SEEN",1);
748 gMC->Gsatt("R3L3","SEEN",1);
749 gMC->Gsatt("R1R2","SEEN",1);
750 gMC->Gsatt("R2R2","SEEN",1);
751 gMC->Gsatt("R3R2","SEEN",1);
752 gMC->Gsatt("R1L2","SEEN",1);
753 gMC->Gsatt("R2L2","SEEN",1);
754 gMC->Gsatt("R3L2","SEEN",1);
755 gMC->Gsatt("R1R1","SEEN",1);
756 gMC->Gsatt("R2R1","SEEN",1);
757 gMC->Gsatt("R3R1","SEEN",1);
758 gMC->Gsatt("R1L1","SEEN",1);
759 gMC->Gsatt("R2L1","SEEN",1);
760 gMC->Gsatt("R3L1","SEEN",1);
761 gMC->Gsatt("R1L4","SEEN",1);
762 gMC->Gsatt("R2L4","SEEN",1);
763 gMC->Gsatt("R3L4","SEEN",1);
764 gMC->Gsatt("BR1_","SEEN",1);
765 gMC->Gsatt("BR2_","SEEN",1);
766 gMC->Gsatt("BR3_","SEEN",1);
767 gMC->Gsatt("BR4_","SEEN",1);
768 gMC->Gsatt("BR5_","SEEN",1);
769 gMC->Gsatt("CB1_","SEEN",1);
770 gMC->Gsatt("CB2_","SEEN",1);
771 gMC->Gsatt("CB3_","SEEN",1);
772 gMC->Gsatt("CB4_","SEEN",1);
773 gMC->Gsatt("CB5_","SEEN",1);
774 gMC->Gsatt("R1SP","SEEN",1);
775 gMC->Gsatt("R2SP","SEEN",1);
776 gMC->Gsatt("R3SP","SEEN",1);
fe4da5cc 777 //
cfce8870 778 gMC->Gdopt("hide", "on");
779 gMC->Gdopt("shad", "on");
780 gMC->Gsatt("*", "fill", 7);
781 gMC->SetClipBox(".");
782 gMC->SetClipBox("*", 0, 1000, -1000, 1000, -1000, 1000);
783 gMC->DefaultRange();
784 gMC->Gdraw("alic", 40, 30, 0, 6, 9, .08, .08);
785 gMC->Gdhead(1111, "Forward Multiplicity Detector version 1");
786 gMC->Gdman(13, 9, "MAN");
fe4da5cc 787}
788
789//_____________________________________________________________________________
790void AliFMDv1::CreateMaterials()
791{
792 //
0d630091 793 // Create Materials for version 1 of FMD
fe4da5cc 794 //
795
0d630091 796 // Material for ring #1
797 Float_t ar1[8] = { 55.8,58.7,52.,47.9,16.,28.,207.2,27. };
798 Float_t zr1[8] = { 26.,28.,24.,22.,8.,14.,82.,13. };
799 Float_t wr1[8] = { .27,.081,.054,.045,.18,.25,.06,.06 };
800 // Material for ring #2
801 Float_t ar2[3] = { 55.8,27.,16. };
802 Float_t zr2[3] = { 26.,13.,8. };
803 Float_t wr2[3] = { .35,.34,.31 };
804 // Material for ring #3
805 Float_t ar3[3] = { 28.,27.,16. };
806 Float_t zr3[3] = { 14.,13.,8. };
807 Float_t wr3[3] = { .37,.33,.3 };
808 // Material for brackets
809 Float_t abr[2] = { 1.,12. };
810 Float_t zbr[2] = { 1.,6. };
811 Float_t wbr[2] = { .1,.9 };
812 // Material for cables
813 Float_t acb[4] = { 1.,12.,37.,63.54 };
814 Float_t zcb[4] = { 1.,6.,17.,2.9 };
815 Float_t wcb[4] = { .02,.14,.12,.72 };
fe4da5cc 816
817 Float_t epsil, stmin, deemax, tmaxfd, stemax;
818
819 Int_t ISXFLD = gAlice->Field()->Integ();
820 Float_t SXMGMX = gAlice->Field()->Max();
821
0d630091 822 // Ring #1
fe4da5cc 823
0d630091 824 AliMixture(1, "FMD_R1$", ar1, zr1, 2.69, 8, wr1);
fe4da5cc 825
0d630091 826 // Ring #2
fe4da5cc 827
0d630091 828 AliMixture(2, "FMD_R2$", ar2, zr2, 2.63, 3, wr2);
fe4da5cc 829
0d630091 830 // Ring #3
fe4da5cc 831
0d630091 832 AliMixture(3, "FMD_R3$", ar3, zr3, 3.15, 3, wr3);
833
834 // Brackets
fe4da5cc 835
0d630091 836 AliMixture(4, "FMD_BR$", abr, zbr, 1.8, 2, wbr);
837
838 // Cables
fe4da5cc 839
0d630091 840 AliMixture(5, "FMD_CB$", acb, zcb, 3.11, 4, wcb);
841
fe4da5cc 842 // *******************************************************
843 // Defines tracking media parameters.
844 // *******************************************************
845 epsil = .001; // Tracking precision, DLS
846 stemax = -1.; // Maximum displacement for multiple scattering
847 tmaxfd = -20.; // Maximum angle due to field deflection
848 deemax = -.3; // Maximum fractional energy loss, DLS
849 stmin = -.8;
850 // ********************************************************
0d630091 851 AliMedium(1, "FMD_R1_L3 ", 1, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
852 AliMedium(2, "FMD_R2_L3 ", 2, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
853 AliMedium(3, "FMD_R3_L3 ", 3, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
854 AliMedium(4, "FMD_BR_L3 ", 4, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
855 AliMedium(5, "FMD_CB_L3 ", 5, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
fe4da5cc 856}
857
0d630091 858
859