]> git.uio.no Git - u/mrichter/AliRoot.git/blame - STRUCT/AliDIPOv2.cxx
Further modifications in OpenOutput and WriteCluster
[u/mrichter/AliRoot.git] / STRUCT / AliDIPOv2.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$
43ea1292 18Revision 1.18 2001/06/21 14:16:43 morsch
19Use absolute coordinates for mother volume PCON.
20
395de7f7 21Revision 1.17 2001/06/21 11:59:25 morsch
22Some more details in compensator geometry.
23
34d8746c 24Revision 1.16 2001/06/20 16:07:08 morsch
25Compensator dipole MBWMD (MCB@SPS) added.
26
c557d782 27Revision 1.15 2001/03/16 15:34:37 morsch
28Mothervolume defined MANY because overlap with station 3 mothervolume not avoidable (A. de Falco)
29
a37078b5 30Revision 1.14 2000/12/21 16:37:23 morsch
31Use Al for coil and cable material. The materials used before cause the dipole to
32have hydrogene on the outer surface leading to unrealistic gamma rates due to
33n-capture.
34
ee7389e1 35Revision 1.13 2000/10/02 21:28:15 fca
36Removal of useless dependecies via forward declarations
37
94de3818 38Revision 1.12 2000/06/20 10:53:01 morsch
39Volume placed outside mother volume (DDIP) corrected (Galina Chabratova)
40
e094dc1b 41Revision 1.11 2000/06/11 12:33:46 morsch
42Coding rule violations corrected
43
b43eb0dc 44Revision 1.10 2000/06/09 19:32:56 morsch
45New detailed and corrected version from Galina Chabratova
46
1b206195 47Revision 1.9 2000/04/27 09:29:53 fca
48Reverting to version 1.6.2
49
d811c12b 50Revision 1.6.2.1 1999/12/03 16:38:51 fca
51Correct overlap in magnet
9e8c9316 52
e821f75f 53Revision 1.6 1999/09/29 09:24:30 fca
54Introduction of the Copyright and cvs Log
55
4c039060 56*/
57
fe4da5cc 58///////////////////////////////////////////////////////////////////////////////
59// //
60// Magnetic Dipole version 1 //
61// //
62//Begin_Html
63/*
d811c12b 64<img src="picts/AliDIPOv2Class.gif">
fe4da5cc 65</pre>
66<br clear=left>
67<font size=+2 color=red>
68<p>The responsible person for this module is
69<a href="mailto:andreas.morsch@cern.ch">Andreas Morsch</a>.
70</font>
71<pre>
72*/
73//End_Html
74// //
75// //
fe4da5cc 76
77#include "AliDIPOv2.h"
78#include "AliRun.h"
94de3818 79#include "AliMagF.h"
fe4da5cc 80#include "AliConst.h"
81
82ClassImp(AliDIPOv2)
83
84//_____________________________________________________________________________
1b206195 85AliDIPOv2::AliDIPOv2()
fe4da5cc 86{
87 //
1b206195 88 // Last design of magnetic dipole version 2
fe4da5cc 89 //
90}
91
92//_____________________________________________________________________________
93AliDIPOv2::AliDIPOv2(const char *name, const char *title)
94 : AliDIPO(name,title)
95{
96 //
97 // Standard constructor for the magnetic dipole version 2
98 SetMarkerColor(7);
99 SetMarkerStyle(2);
100 SetMarkerSize(0.4);
101}
c557d782 102
fe4da5cc 103void AliDIPOv2::CreateGeometry()
c557d782 104{
105 //
106 // Creation of the geometry of the magnetic DIPOLE version 2
107 //
108
109 CreateSpectrometerDipole();
110 CreateCompensatorDipole();
111}
112
113
114//_____________________________________________________________________________
115void AliDIPOv2::CreateSpectrometerDipole()
fe4da5cc 116{
117 //
118 // Creation of the geometry of the magnetic DIPOLE version 2
119 //
1b206195 120
395de7f7 121 Float_t cpar[5], tpar[18], ypar[12];
fe4da5cc 122 Float_t dz, dx, dy;
123 Int_t idrotm[1899];
b43eb0dc 124 Float_t accMax, the1, phi1, the2, phi2, the3, phi3;
d811c12b 125
1b206195 126 Int_t *idtmed = fIdtmed->GetArray()-1799;
127
ee7389e1 128// const Int_t kCoil = 1813;
129// const Int_t kCable= 1811;
130
131 const Int_t kCoil = 1808;
132 const Int_t kCable= 1808;
133
b43eb0dc 134 accMax = 9.; // ANGLE POLAIRE MAXIMUM
1b206195 135
fe4da5cc 136 // DIPOLE MAGNET
395de7f7 137 const Float_t kZDipole = 975;
1b206195 138
e094dc1b 139 tpar[0] = 0.;
140 tpar[1] = 360.;
395de7f7 141 tpar[2] = 5.;
e094dc1b 142 //
395de7f7 143 tpar[3] = -250.55+kZDipole;
144 tpar[4] = 30.5;
e094dc1b 145 tpar[5] = 527.34;
146 //
395de7f7 147 tpar[6] = -160.7+kZDipole;
148 tpar[7] = 30.5;
e094dc1b 149 tpar[8] = 527.34;
150 //
395de7f7 151 tpar[9] = 30.+kZDipole;
152 tpar[10] = 30.5;
e094dc1b 153 tpar[11] = 527.34;
395de7f7 154
155 tpar[12] = 150.8+kZDipole;
e094dc1b 156 tpar[13] = 193.3;
157 tpar[14] = 527.34;
395de7f7 158 //
159 tpar[15] = 250.55+kZDipole;
160 tpar[16] = 193.3;
161 tpar[17] = 527.34;
e094dc1b 162
1b206195 163
395de7f7 164 gMC->Gsvolu("DDIP", "PCON", idtmed[1814], tpar, 18);
fe4da5cc 165 // COILS
fe4da5cc 166 // air - m.f.
1b206195 167 cpar[0] = 207.;
168 cpar[1] = 274.;
169 cpar[2] = 37.65;
170 cpar[3] = 119.;
395de7f7 171 cpar[4] = 241.;
1b206195 172 // coil - high cuts
ee7389e1 173 gMC->Gsvolu("DC1 ", "TUBS", idtmed[kCoil+40], cpar, 5);
1b206195 174 cpar[3] = -61.;
175 cpar[4] = 61.;
ee7389e1 176 gMC->Gsvolu("DC2 ", "TUBS", idtmed[kCoil+40], cpar, 5);
1b206195 177
178 // coil - low cuts cuts
179 cpar[0] = 207.;
180// cpar[1] = cpar[0] + 10.;
181 cpar[1] = 217;
182 cpar[3] = 119.;
183 cpar[4] = 241.;
184
ee7389e1 185 gMC->Gsvolu("DC3 ", "TUBS", idtmed[kCoil], cpar, 5);
1b206195 186 cpar[0] = 207.;
187 cpar[1] = 217;
188 cpar[3] = -61.;
189 cpar[4] = 61.;
ee7389e1 190 gMC->Gsvolu("DC4 ", "TUBS", idtmed[kCoil], cpar, 5);
1b206195 191
cfce8870 192 gMC->Gspos("DC3 ", 1, "DC1 ", 0., 0., 0., 0, "ONLY");
cfce8870 193 gMC->Gspos("DC4 ", 1, "DC2 ", 0., 0., 0., 0, "ONLY");
1b206195 194
195// dz = 37.65 - 243.55
196 dz = -205.9-2.45;
197 dx = 5.;
395de7f7 198 gMC->Gspos("DC1 ", 1, "DDIP", dx, 0., dz+kZDipole, 0, "ONLY");
199 gMC->Gspos("DC1 ", 2, "DDIP", dx, 0., -dz+kZDipole, 0, "ONLY");
200 gMC->Gspos("DC2 ", 1, "DDIP", -dx, 0., dz+kZDipole, 0, "ONLY");
201 gMC->Gspos("DC2 ", 2, "DDIP", -dx, 0., -dz+kZDipole, 0, "ONLY");
fe4da5cc 202 the1 = 180.;
203 phi1 = 0.;
204 the2 = 90.;
1b206195 205 phi2 = 151.;
fe4da5cc 206 the3 = 90.;
1b206195 207 phi3 = 61.;
fe4da5cc 208 AliMatrix(idrotm[1800], the1, phi1, the2, phi2, the3, phi3);
1b206195 209 phi2 = 29.; //90-61
fe4da5cc 210 the3 = -90.;
1b206195 211 phi3 = -61.;
fe4da5cc 212 AliMatrix(idrotm[1801], the1, phi1, the2, phi2, the3, phi3);
213 the1 = 0.;
214 phi1 = 0.;
215 the2 = 90.;
1b206195 216 phi2 = 151.;
fe4da5cc 217 the3 = 90.;
1b206195 218 phi3 = 61.;
fe4da5cc 219 AliMatrix(idrotm[1802], the1, phi1, the2, phi2, the3, phi3);
1b206195 220 phi2 = 29.;
fe4da5cc 221 the3 = -90.;
1b206195 222 phi3 = -61.;
fe4da5cc 223 AliMatrix(idrotm[1803], the1, phi1, the2, phi2, the3, phi3);
1b206195 224
fe4da5cc 225 cpar[0] = 25.;
1b206195 226 cpar[1] = 100.3; //25+75.3
227 cpar[2] = 33.5;
fe4da5cc 228 cpar[3] = 270.;
229 cpar[4] = 360.;
1b206195 230//* coil high cuts
ee7389e1 231 gMC->Gsvolu("DC11", "TUBS", idtmed[kCoil+40], cpar, 5);
1b206195 232
395de7f7 233 dx = TMath::Sin(30.5*kDegrad) * -(207.+33.5)+5./TMath::Sin(30.5*kDegrad);
1b206195 234 dy = TMath::Cos(30.5*kDegrad) * -(207.+33.5);
235 dz = cpar[1] - 243.55-2.45;
395de7f7 236 gMC->Gspos("DC11", 1, "DDIP", dx, dy, dz+kZDipole, idrotm[1800], "ONLY");
237 gMC->Gspos("DC11", 2, "DDIP", dx, dy, -dz+kZDipole, idrotm[1802], "ONLY");
238 gMC->Gspos("DC11", 3, "DDIP", -dx, dy, dz+kZDipole, idrotm[1801], "ONLY");
239 gMC->Gspos("DC11", 4, "DDIP", -dx, dy, -dz+kZDipole, idrotm[1803], "ONLY");
1b206195 240
241
242
243//* ... higher cuts
fe4da5cc 244 cpar[0] = 25.;
1b206195 245 cpar[1] = 100.3; //25+75.3
246 cpar[2] = 33.5;
fe4da5cc 247 cpar[3] = 0.;
248 cpar[4] = 90.;
1b206195 249//* coil high cuts
ee7389e1 250 gMC->Gsvolu("DC12", "TUBS", idtmed[kCoil+40], cpar, 5);
1b206195 251
395de7f7 252 dx = TMath::Sin(30.5*kDegrad) * -(207.+33.5)+5./TMath::Sin(30.5*kDegrad);
1b206195 253 dy = TMath::Cos(30.5*kDegrad) *(207.+33.5);
254 dz = cpar[1] - 243.55-2.45;
395de7f7 255 gMC->Gspos("DC12", 1, "DDIP", dx, dy, dz+kZDipole, idrotm[1801], "ONLY");
256 gMC->Gspos("DC12", 2, "DDIP", dx, dy, -dz+kZDipole, idrotm[1803], "ONLY");
257 gMC->Gspos("DC12", 3, "DDIP", -dx, dy, dz+kZDipole, idrotm[1800], "ONLY");
258 gMC->Gspos("DC12", 4, "DDIP", -dx, dy, -dz+kZDipole, idrotm[1802], "ONLY");
1b206195 259
fe4da5cc 260 the1 = 90.;
1b206195 261 phi1 = 61.;
fe4da5cc 262 the2 = 90.;
1b206195 263 phi2 = 151.;
fe4da5cc 264 the3 = 0.;
265 phi3 = 0.;
266 AliMatrix(idrotm[1804], the1, phi1, the2, phi2, the3, phi3);
267 the1 = 90.;
1b206195 268 phi1 = -61.;
fe4da5cc 269 the2 = 90.;
1b206195 270 phi2 = -151.;
fe4da5cc 271 AliMatrix(idrotm[1805], the1, phi1, the2, phi2, the3, phi3);
1b206195 272 the1 = 90.;
273 phi1 = 119.; //180 -61
274 the2 = 90.;
275 phi2 = 209.; //270-61
276 AliMatrix(idrotm[1806], the1, phi1, the2, phi2, the3, phi3);
277 the1 = 90.;
278 phi1 = -119.;
279 the2 = 90.;
280 phi2 = -209.;
281 AliMatrix(idrotm[1807], the1, phi1, the2, phi2, the3, phi3);
282
283//* coil - high cuts
284
285 tpar[0] = 37.65;
286 tpar[1] = 33.5;
287 tpar[2] = 145.5;
ee7389e1 288 gMC->Gsvolu("DL1 ", "BOX ", idtmed[kCoil+40], tpar, 3);
1b206195 289
290// coil - low cuts
291
292 tpar[0] = 5.;
293 dx = 37.65 - 5.;
ee7389e1 294 gMC->Gsvolu("DL2 ", "BOX ", idtmed[kCoil], tpar, 3);
1b206195 295 gMC->Gspos("DL2 ", 1, "DL1 ", dx, 0., 0., 0, "ONLY");
296
297 dx =-53.62;
298 dy =-241.26819;
395de7f7 299 dz = 0.;
300 gMC->Gspos("DL1 ", 1, "DDIP", dx, dy, dz+kZDipole, idrotm[1804], "ONLY");
301 gMC->Gspos("DL1 ", 2, "DDIP", dx, -dy, dz+kZDipole, idrotm[1805], "ONLY");
302 gMC->Gspos("DL1 ", 3, "DDIP",-dx, dy, dz+kZDipole, idrotm[1806], "ONLY");
303 gMC->Gspos("DL1 ", 4, "DDIP",-dx, -dy, dz+kZDipole, idrotm[1807], "ONLY");
1b206195 304
305 // Contactor
306
307 // high cuts
308
309 //Steel outer face planes
310
311 cpar[0] = 207.-18.6;
312 cpar[1] = 274.+18.6;
313 cpar[2] = 1.;
314 cpar[3] = -50.;
315 cpar[4] = 50.;
316
317 gMC->Gsvolu("DCO1", "TUBS", idtmed[1818], cpar, 5);
318 dx = -5.;
319 dz = 168.25-1.5-1.;
395de7f7 320 gMC->Gspos("DCO1", 1, "DDIP", dx, 0, dz+kZDipole, 0, "ONLY");
1b206195 321 dz = 243.55+4.5+1.5+1.;
395de7f7 322 gMC->Gspos("DCO1", 2, "DDIP", dx, 0, dz+kZDipole, 0, "ONLY");
d811c12b 323
1b206195 324 // 9.06.2000
325
326 // cpar[0] = 207.-18.6;
327 // cpar[1] = 274.+18.6;
328 // cpar[2] = 1.;
329 cpar[3] = 180.-50.;
330 cpar[4] = 180.+50.;
331
332 gMC->Gsvolu("DCO2", "TUBS", idtmed[1818], cpar, 5);
333 dx = +5.;
334 dz = 168.25-1.5-1.;
395de7f7 335 gMC->Gspos("DCO2", 1, "DDIP", dx, 0, dz+kZDipole, 0, "ONLY");
1b206195 336 dz = 243.55+4.5+1.5+1.;
395de7f7 337 gMC->Gspos("DCO2", 2, "DDIP", dx, 0, dz+kZDipole, 0, "ONLY");
1b206195 338
339
340
341 // Resin face planes
342
343 cpar[0] = 207.;
344 cpar[1] = 274.;
345 cpar[2] = .75;
346 cpar[3] = -50.;
347 cpar[4] = 50.;
348
349 gMC->Gsvolu("DCO3", "TUBS", idtmed[1812], cpar, 5);
350 dx = -5;
351 dz = 168.25-0.75;
395de7f7 352 gMC->Gspos("DCO3", 1, "DDIP", dx, 0, dz+kZDipole, 0, "ONLY");
1b206195 353 dz = 243.55+4.5+0.75;
395de7f7 354 gMC->Gspos("DCO3", 2, "DDIP", dx, 0, dz+kZDipole, 0, "ONLY");
1b206195 355
356 // 9.06.2000
357
358 cpar[3] = 180.-50.;
359 cpar[4] = 180.+50.;
360 gMC->Gsvolu("DCO4", "TUBS", idtmed[1812], cpar, 5);
361 dx = +5;
362 dz = 168.25-0.75;
395de7f7 363 gMC->Gspos("DCO4", 1, "DDIP", dx, 0, dz+kZDipole, 0, "ONLY");
364 dz = 243.55+4.5+0.75 ;
365 gMC->Gspos("DCO4", 2, "DDIP", dx, 0, dz+kZDipole, 0, "ONLY");
1b206195 366
367
368 // G10 face plane
369
370 cpar[0] = 207.;
371 cpar[1] = 274.;
372 cpar[2] = 2.25;
373 cpar[3] = -50.;
374 cpar[4] = 50.;
375
376 gMC->Gsvolu("DCO5", "TUBS", idtmed[1810], cpar, 5);
377
378 dx = -5;
379 dz = 243.55+2.25;
395de7f7 380 gMC->Gspos("DCO5", 1, "DDIP", dx, 0, dz+kZDipole, 0, "ONLY");
1b206195 381
382 // 9.06.2000
383
384 cpar[3] = 180.-50.;
385 cpar[4] = 180.+50.;
386
387 gMC->Gsvolu("DCO6", "TUBS", idtmed[1810], cpar, 5);
388
389 dx = +5;
395de7f7 390 dz = 243.55+2.25;
391 gMC->Gspos("DCO6", 1, "DDIP", dx, 0, dz+kZDipole, 0, "ONLY");
1b206195 392
393 //Steel supported planes
394
395 cpar[0] = 274.+1.5+2.;
396 cpar[1] = 274.+18.6;
397 cpar[2] = 1.;
398 cpar[3] = -50.;
399 cpar[4] = 50.;
400
401 gMC->Gsvolu("DCO7", "TUBS", idtmed[1818], cpar, 5);
402
403 dx = -5;
404 dz = 168.25+1.;
395de7f7 405 gMC->Gspos("DCO7", 1, "DDIP", dx, 0, dz+kZDipole, 0, "ONLY");
1b206195 406
407 // 9.06.2000
408 cpar[0] = 274.+1.5+2.;
409 cpar[1] = 274.+18.6;
410 cpar[2] = 1.;
411 cpar[3] = 180.-50.;
412 cpar[4] = 180.+50.;
413
414
415 gMC->Gsvolu("DCO8", "TUBS", idtmed[1818], cpar, 5);
416
417 dx = +5;
418 dz = 168.25+1.;
395de7f7 419 gMC->Gspos("DCO8", 1, "DDIP", dx, 0, dz+kZDipole, 0, "ONLY");
1b206195 420
421 //
422
423 cpar[0] = 207.- 18.6;
424 cpar[1] = 207.- 2.- 1.5;
425 cpar[2] = 1.;
426 cpar[3] = -50.;
427 cpar[4] = 50.;
428
429 gMC->Gsvolu("DCO9", "TUBS", idtmed[1818], cpar, 5);
430
431 dx = -5;
432 dz = 168.25+1.;
395de7f7 433 gMC->Gspos("DCO9", 1, "DDIP", dx, 0, dz+kZDipole, 0, "ONLY");
1b206195 434
435 // 9.06.2000
436
437 cpar[0] = 207.- 18.6;
438 cpar[1] = 207.- 2.- 1.5;
439 cpar[2] = 1.;
440 cpar[3] = 180.-50.;
441 cpar[4] = 180.+50.;
442
443 gMC->Gsvolu("DCOA", "TUBS", idtmed[1818], cpar, 5);
444
445 dx = +5;
446 dz = 168.25+1.;
395de7f7 447 gMC->Gspos("DCOA", 1, "DDIP", dx, 0, dz+kZDipole, 0, "ONLY");
1b206195 448
449
450 // Sides steel planes
451
452 cpar[0] = 207. - 1.5 -2.;
395de7f7 453 cpar[1] = 207. - 1.5;
1b206195 454 cpar[2] = ((243.55+4.5+1.5)-168.25)/2;
455 cpar[3] = -50.;
456 cpar[4] = 50.;
457
458 gMC->Gsvolu("DCOB", "TUBS", idtmed[1818], cpar, 5);
459
460 cpar[0] = 274. + 1.5;
461 cpar[1] = 274. + 1.5 +2.;
462
463 gMC->Gsvolu("DCOC", "TUBS", idtmed[1818], cpar, 5);
464
465 dx=-5.;
466 dz = ((243.55+4.5+1.5)+168.25)/2;
395de7f7 467 gMC->Gspos("DCOB", 1, "DDIP", dx, 0, dz+kZDipole, 0, "ONLY");
468 gMC->Gspos("DCOC", 1, "DDIP", dx, 0, dz+kZDipole, 0, "ONLY");
1b206195 469
470 // 9.06.2000
471
472 cpar[0] = 207. - 1.5 -2.;
395de7f7 473 cpar[1] = 207. - 1.5;
1b206195 474 cpar[2] = ((243.55+4.5+1.5)-168.25)/2;
475 cpar[3] = 180.-50.;
476 cpar[4] = 180.+50.;
477
478 gMC->Gsvolu("DCOD", "TUBS", idtmed[1818], cpar, 5);
479
480 cpar[0] = 274. + 1.5;
481 cpar[1] = 274. + 1.5 +2.;
482
483 gMC->Gsvolu("DCOE", "TUBS", idtmed[1818], cpar, 5);
484
485 dx=+5.;
486 dz = ((243.55+4.5+1.5)+168.25)/2;
395de7f7 487 gMC->Gspos("DCOD", 1, "DDIP", dx, 0, dz+kZDipole, 0, "ONLY");
488 gMC->Gspos("DCOE", 1, "DDIP", dx, 0, dz+kZDipole, 0, "ONLY");
1b206195 489
490
491 // Top and bottom resin planes
492
395de7f7 493 cpar[0] = 207. - 1.5;
494 cpar[1] = 207.;
1b206195 495 cpar[2] = ((243.55+4.5+1.5)-168.25)/2;
496 cpar[3] = -50.;
497 cpar[4] = 50.;
498
499 gMC->Gsvolu("DCOF", "TUBS", idtmed[1812], cpar, 5);
500
501 cpar[0] = 274.;
502 cpar[1] = 274. + 1.5;
503
504 gMC->Gsvolu("DCOG", "TUBS", idtmed[1812], cpar, 5);
505
506
507 dx=-5.;
508 dz = ((243.55+4.5+1.5)+168.25)/2;
395de7f7 509 gMC->Gspos("DCOF", 1, "DDIP", dx, 0, dz+kZDipole, 0, "ONLY");
510 gMC->Gspos("DCOG", 1, "DDIP", dx, 0, dz+kZDipole, 0, "ONLY");
1b206195 511
512 // 9.06.2000
395de7f7 513 cpar[0] = 207. - 1.5;
514 cpar[1] = 207.;
1b206195 515 cpar[2] = ((243.55+4.5+1.5)-168.25)/2;
516
517 cpar[3] = 180.-50.;
518 cpar[4] = 180.+50.;
519
520 gMC->Gsvolu("DCOH", "TUBS", idtmed[1812], cpar, 5);
521
522 cpar[0] = 274.;
523 cpar[1] = 274. + 1.5;
524
525 gMC->Gsvolu("DCOI", "TUBS", idtmed[1812], cpar, 5);
526
527
528 dx=+5.;
529 dz = ((243.55+4.5+1.5)+168.25)/2;
395de7f7 530 gMC->Gspos("DCOH", 1, "DDIP", dx, 0, dz+kZDipole, 0, "ONLY");
531 gMC->Gspos("DCOI", 1, "DDIP", dx, 0, dz+kZDipole, 0, "ONLY");
1b206195 532
533
534 // Aluminum cabels
535
536 cpar[0] = 274. + 1.5 +2.;
537 cpar[1] = 274. + 1.5 +2. + 80.;
538 cpar[2] = 5.05/2;
539 cpar[3] = -24.;
540 cpar[4] = 24.;
541
ee7389e1 542 gMC->Gsvolu("DCOJ", "TUBS", idtmed[kCable], cpar, 5);
1b206195 543
544 // dx = 274. + 1.5 +2. +40.;
545 // dx = 5. + 1.5 +2. +40.;
546 // dx = 5. + 1.5 +2.;
547 dx=-5.;
548 dz = 168.25 + 5.05 + 5.05/2;
395de7f7 549 gMC->Gspos("DCOJ", 1, "DDIP", dx, 0, dz+kZDipole, 0, "ONLY");
1b206195 550
551 dz = 243.55 - 5.05/2;
395de7f7 552 gMC->Gspos("DCOJ", 2, "DDIP", dx, 0, dz+kZDipole, 0, "ONLY");
1b206195 553
554 // 9.06.2000
555
556 cpar[3] = 180.-24.;
557 cpar[4] = 180.+24.;
558
ee7389e1 559 gMC->Gsvolu("DCOK", "TUBS", idtmed[kCable], cpar, 5);
1b206195 560
561 // dx = 274. + 1.5 +2. +40.;
562 // dx = 5. + 1.5 +2. +40.;
563 // dx = 5. + 1.5 +2.;
564 dx=+5.;
565 dz = 168.25 + 5.05 + 5.05/2;
395de7f7 566 gMC->Gspos("DCOK", 1, "DDIP", dx, 0, dz+kZDipole, 0, "ONLY");
1b206195 567
568 dz = 243.55 - 5.05/2;
395de7f7 569 gMC->Gspos("DCOK", 2, "DDIP", dx, 0, dz+kZDipole, 0, "ONLY");
1b206195 570
571
fe4da5cc 572 // YOKE
1b206195 573
574// Top and bottom blocks
395de7f7 575 ypar[0] = 298.1;
1b206195 576 ypar[1] = 69.5;
577 ypar[2] = 155.75;
578
579// iron- high cuts
580 gMC->Gsvolu("DY1 ", "BOX ", idtmed[1858], ypar, 3);
395de7f7 581 ypar[0] = 144.+10.;
1b206195 582 ypar[1] = 193.3+10.;
fe4da5cc 583 ypar[2] = 5.;
1b206195 584 ypar[3] = 155.75;
585 dy = -69.5 + 5.;
586// iron- low cuts
587 gMC->Gsvolu("DY11", "TRD1", idtmed[1818], ypar, 4);
588 gMC->Gspos("DY11", 1, "DY1 ", 0., dy, 0., 0, "ONLY");
589
590 dy = 365.5;
591 dz = 4.95;
395de7f7 592 gMC->Gspos("DY1 ", 1, "DDIP", 0., dy, -dz+kZDipole, 0, "ONLY");
1b206195 593
594 the1 = 270.;
fe4da5cc 595 phi1 = 0.;
1b206195 596 the2 = 270.;
fe4da5cc 597 phi2 = 90.;
1b206195 598 the3 = 0.;
fe4da5cc 599 phi3 = 0.;
1b206195 600 AliMatrix(idrotm[1808], the1, phi1, the2, phi2, the3, phi3);
395de7f7 601 gMC->Gspos("DY1 ", 2, "DDIP", 0., -dy, -dz+kZDipole, idrotm[1808] , "ONLY");
1b206195 602
603// side walls
395de7f7 604 // ypar[0] = 579./2.;
605 ypar[0] = 296.;
1b206195 606 ypar[1] = 0.;
607 ypar[2] = 0.;
608 ypar[3] = 155.75;
395de7f7 609 ypar[4] = 47.9;
1b206195 610 ypar[5] = 72.55;
395de7f7 611 ypar[6] = 4.3058039629;
1b206195 612 // z+
613 ypar[7] = 155.75;
395de7f7 614 ypar[8] = 47.9;
1b206195 615 ypar[9] = 72.55;
395de7f7 616 ypar[10] = 4.3058039629;
1b206195 617
618// iron - high cuts
619
620 gMC->Gsvolu("DY2 ", "TRAP", idtmed[1858], ypar,11);
621
622 ypar[4] = 47.9 -5.;
623 ypar[5] = 72.55 -5.;
624
625 ypar[8] = 47.9 -5.;
626 ypar[9] = 72.55 -5.;
627
628
629// iron - low cuts
630
631 gMC->Gsvolu("DY22", "TRAP", idtmed[1818], ypar,11);
632
633 dy = 0.;
634 dx = -5.;
635
636 gMC->Gspos("DY22", 1, "DY2 ", dx, dy, 0., 0, "ONLY");
637
638 the1 = 90.;
639 phi1 = 180.;
640 the2 = 180.;
641 phi2 = 180.;
642 the3 = 90.;
643 phi3 = 90.;
644 AliMatrix(idrotm[1809], the1, phi1, the2, phi2, the3, phi3);
645
646 the1 = 90.;
fe4da5cc 647 phi1 = 0.;
1b206195 648 the2 = 180.;
649 phi2 = 0.;
395de7f7 650 the3 = 90.;
1b206195 651 phi3 = 90.;
652 AliMatrix(idrotm[1810], the1, phi1, the2, phi2, the3, phi3);
653
654 dx = 228.875;
655 dz = - 4.95;
395de7f7 656
657 gMC->Gspos("DY2 ", 1, "DDIP", dx, 0.0, dz+kZDipole, idrotm[1809], "ONLY");
658 gMC->Gspos("DY2 ", 2, "DDIP", -dx, 0.0, dz+kZDipole, idrotm[1810], "ONLY");
1b206195 659
395de7f7 660 gMC->Gspos("DDIP", 1, "ALIC", 0., 0., 0., 0, "ONLY");
1b206195 661
cfce8870 662 gMC->Gsatt("DDIP", "SEEN", 0);
1b206195 663// gMC->Gsatt("DC21", "SEEN", 0);
664// gMC->Gsatt("DC22", "SEEN", 0);
665// gMC->Gsatt("DC3 ", "SEEN", 0);
666// gMC->Gsatt("DC4 ", "SEEN", 0);
fe4da5cc 667}
668
c557d782 669
670void AliDIPOv2::CreateCompensatorDipole()
671{
672 //
673 // Geometry of the compensator Dipole MBWMD (was MCB @ SPS)
674 //
675 Int_t *idtmed = fIdtmed->GetArray()-1799;
676 Int_t idrotm[1899];
677//
678 Float_t pbox[3] = {62.5, 62.5, 170.};
679
680// Mother volumes
681 gMC->Gsvolu("DCM0", "BOX", idtmed[1814], pbox, 3);
34d8746c 682
c557d782 683//
34d8746c 684// Mother volume containing lower coil
c557d782 685 pbox[0] = 58.5/2.;
686 pbox[1] = 30.0;
34d8746c 687 pbox[2] = 152.5;
688
c557d782 689 gMC->Gsvolu("DCML", "BOX", idtmed[1809], pbox, 3);
690//
691// Base
692 pbox[0] = 62.5;
693 pbox[1] = 15.0;
694 gMC->Gsvolu("DCBA", "BOX", idtmed[1809], pbox, 3);
c557d782 695//
696// Coil: straight sections, horizontal
697 pbox[0] = 6.;
698 pbox[1] = 11.;
34d8746c 699 pbox[2] = 135.;
c557d782 700 gMC->Gsvolu("DCH1", "BOX", idtmed[1816], pbox, 3);
701//
702// Coil: straight sections, horizontal
703 pbox[0] = 6.;
704 pbox[1] = 11.;
705 pbox[2] = 135.;
706 gMC->Gsvolu("DCH2", "BOX", idtmed[1816], pbox, 3);
707
708//
709// Mother volume containing upper coil
710 pbox[0] = 8.0;
711 pbox[1] = 17.5;
712 pbox[2] = 135.0;
713 gMC->Gsvolu("DCMU", "BOX", idtmed[1809], pbox, 3);
714
715//
716// Coil: straight sections, vertical
34d8746c 717 pbox[0] = 6.0;
718 pbox[1] = 9.5;
719 pbox[2] = 11.0;
c557d782 720
721 gMC->Gsvolu("DCCV", "BOX", idtmed[1816], pbox, 3);
722//
723// Coil: circular section
724
725 Float_t ptubs[5];
726 ptubs[0] = 0.;
727 ptubs[1] = 35.;
728 ptubs[2] = 8.;
729 ptubs[3] = 0.;
730 ptubs[4] = 90.;
34d8746c 731// gMC->Gsvolu("DCC1", "TUBS", idtmed[1809], ptubs, 5);
c557d782 732 ptubs[0] = 13.;
733 ptubs[1] = 35.;
734 ptubs[2] = 6.;
34d8746c 735 ptubs[3] = 0.;
736 ptubs[4] = 90.;
737 gMC->Gsvolu("DCC1", "TUBS", idtmed[1816], ptubs, 5);
738//
739// Clamps
740 Float_t ppgon[10];
741 ppgon[0] = 0.;
742 ppgon[1] = 90.;
743 ppgon[2] = 1.;
744 ppgon[3] = 2.;
745 ppgon[4] = -1.;
746 ppgon[5] = 0.;
747 ppgon[6] = 24.75;
748 ppgon[7] = 1.;
749 ppgon[8] = 0.;
750 ppgon[9] = 24.75;
751 gMC->Gsvolu("DCLA", "PGON", idtmed[1809], ppgon, 10);
c557d782 752//
753// Assemble all
754//
755 AliMatrix(idrotm[1811], -90., 0., 90., 90., 0., 0.);
34d8746c 756 AliMatrix(idrotm[1812], 0., 0., 90., 90., 90., 0.);
c557d782 757 AliMatrix(idrotm[1813], 180., 0., 90., 90., 90., 0.);
34d8746c 758 AliMatrix(idrotm[1814], 0., 180., 90., 270., 90., 0.);
759 AliMatrix(idrotm[1815], 180., 180., 90., 270., 90., 0.);
760
761 gMC->Gspos("DCH1", 1, "DCML", 23.25, -13., -17.5, 0, "ONLY");
762 gMC->Gspos("DCCV", 1, "DCM0", 12., 19., -159., 0, "ONLY");
763 gMC->Gspos("DCCV", 2, "DCM0", -12., 19., -159., 0, "ONLY");
764 gMC->Gspos("DCCV", 3, "DCML", 23.25, 20.5, 141.5, 0, "ONLY");
c557d782 765
34d8746c 766 gMC->Gspos("DCML", 1, "DCM0", -33.25, -2.5, 17.5, 0, "ONLY");
767 gMC->Gspos("DCML", 2, "DCM0", 33.25, -2.5, 17.5, idrotm[1811], "ONLY");
c557d782 768
769
34d8746c 770 gMC->Gspos("DCH2", 1, "DCMU", 2., 6.5, 0., 0, "ONLY");
c557d782 771 gMC->Gspos("DCMU", 1, "DCM0", -12., 45., 0., 0, "ONLY");
772 gMC->Gspos("DCMU", 2, "DCM0", 12., 45., 0., idrotm[1811], "ONLY");
773
34d8746c 774// gMC->Gspos("DCC2", 1, "DCC1", 0., 0., 0., 0, "ONLY");
c557d782 775
776 gMC->Gspos("DCC1", 1, "DCM0", -12., 27.5, 135., idrotm[1812], "ONLY");
777 gMC->Gspos("DCC1", 2, "DCM0", 12., 27.5, 135., idrotm[1812], "ONLY");
778 gMC->Gspos("DCC1", 3, "DCM0", -12., 27.5, -135., idrotm[1813], "ONLY");
779 gMC->Gspos("DCC1", 4, "DCM0", 12., 27.5, -135., idrotm[1813], "ONLY");
780
34d8746c 781 gMC->Gspos("DCC1", 5, "DCM0", 12., 27.5-32.+13., -135., idrotm[1815], "ONLY");
782 gMC->Gspos("DCC1", 6, "DCM0", -12., 27.5-32.+13., -135., idrotm[1815], "ONLY");
783
784 gMC->Gspos("DCC1", 7, "DCML", 23.25, -13+13.+11., 117.5, idrotm[1814], "ONLY");
785
786 gMC->Gspos("DCLA", 1, "DCM0", 20., 27.5, -134., 0, "ONLY");
787 gMC->Gspos("DCLA", 2, "DCM0", 20., 27.5, -44., 0, "ONLY");
788 gMC->Gspos("DCLA", 3, "DCM0", 20., 27.5, 46., 0, "ONLY");
789 gMC->Gspos("DCLA", 4, "DCM0", 20., 27.5, 134., 0, "ONLY");
790
791 gMC->Gspos("DCLA", 5, "DCM0", -20., 27.5, -134., idrotm[1811], "ONLY");
792 gMC->Gspos("DCLA", 6, "DCM0", -20., 27.5, -44., idrotm[1811], "ONLY");
793 gMC->Gspos("DCLA", 7, "DCM0", -20., 27.5, 46., idrotm[1811], "ONLY");
794 gMC->Gspos("DCLA", 8, "DCM0", -20., 27.5, 134., idrotm[1811], "ONLY");
c557d782 795
c557d782 796
34d8746c 797 gMC->Gspos("DCBA", 1, "DCM0", 0., -47.5 , 17.5, 0, "ONLY");
c557d782 798 gMC->Gspos("DCM0", 1, "ALIC", 0., -6.75, -975., 0, "ONLY");
799
800
801}
802
fe4da5cc 803//_____________________________________________________________________________
b8032157 804void AliDIPOv2::DrawModule()
fe4da5cc 805{
806 //
807 // Draw a shaded view of the muon absorber
808 //
809
fe4da5cc 810 // Set everything unseen
cfce8870 811 gMC->Gsatt("*", "seen", -1);
fe4da5cc 812 //
813 // Set ALIC mother transparent
cfce8870 814 gMC->Gsatt("ALIC","SEEN",0);
fe4da5cc 815 //
816 // Set the volumes visible
cfce8870 817 gMC->Gsatt("DDIP","seen",0);
818 gMC->Gsatt("DC1 ","seen",1);
819 gMC->Gsatt("DC2 ","seen",1);
820 gMC->Gsatt("DC3 ","seen",1);
821 gMC->Gsatt("DC4 ","seen",1);
822 gMC->Gsatt("DC11","seen",1);
823 gMC->Gsatt("DC21","seen",1);
824 gMC->Gsatt("DC12","seen",1);
825 gMC->Gsatt("DC22","seen",1);
826 gMC->Gsatt("DL1 ","seen",1);
827 gMC->Gsatt("DL2 ","seen",1);
828 gMC->Gsatt("DY1 ","seen",1);
829 gMC->Gsatt("DY2 ","seen",1);
830 gMC->Gsatt("DYL ","seen",1);
831 gMC->Gsatt("DY3 ","seen",1);
1b206195 832 // gMC->Gsatt("DY4 ","seen",1);
833 // gMC->Gsatt("DY5 ","seen",1);
834 // gMC->Gsatt("DY6 ","seen",1);
835// gMC->Gsatt("DY7 ","seen",1);
fe4da5cc 836 //
cfce8870 837 gMC->Gdopt("hide", "on");
838 gMC->Gdopt("shad", "on");
839 gMC->Gsatt("*", "fill", 7);
840 gMC->SetClipBox(".");
841 gMC->SetClipBox(".");
842 gMC->DefaultRange();
843 gMC->Gdraw("alic", 30, 30, 0, 17, 13.5, .019, .019);
844 gMC->Gdhead(1111, "Magnetic Dipole Version 2");
845 gMC->Gdman(16, 4, "MAN");
fe4da5cc 846}
847
c557d782 848
849
850
fe4da5cc 851//_____________________________________________________________________________
852void AliDIPOv2::CreateMaterials()
853{
854 //
855 // Create Materials for Magnetic Dipole version 2
856 //
857
b43eb0dc 858 Int_t isxfld = gAlice->Field()->Integ();
859 Float_t sxmgmx = gAlice->Field()->Max();
fe4da5cc 860
861 Float_t asteel[4] = { 55.847,51.9961,58.6934,28.0855 };
862 Float_t zsteel[4] = { 26.,24.,28.,14. };
863 Float_t wsteel[4] = { .715,.18,.1,.005 };
1b206195 864
fe4da5cc 865 Float_t acoil[3] = { 26.98,1.01,16. };
866 Float_t zcoil[3] = { 13.,1.,8. };
867 Float_t wcoil[3] = { .66,.226,.114 };
1b206195 868
869 Float_t aresi[3] = { 1.01,12.011,16.};
870 Float_t zresi[3] = { 1.,6.,8. };
871 Float_t wresi[3] = { .0644,.7655,.1701 };
872
873 Float_t aG10[5] = { 1.01,12.011,16.,28.085 ,79.904 };
874 Float_t zG10[5] = { 1.,6.,8.,14.,35. };
875 Float_t wG10[5] = { .02089,.22338,.28493,.41342,.05738 };
876
877 Float_t aAlCon[2] = { 14.61, 26.98};
878 Float_t zAlCon[2] = { 7.3, 13.};
879 Float_t wAlCon[2] = { .0004,.9996};
880
fe4da5cc 881
882 Float_t epsil, stmin, deemax, tmaxfd, stemax;
883
884 // --- Define the various materials for GEANT ---
885 // Aluminum
c557d782 886 AliMaterial( 9, "ALUMINIUM$", 26.98, 13., 2.7, 8.9, 37.2);
fe4da5cc 887 AliMaterial(29, "ALUMINIUM$", 26.98, 13., 2.7, 8.9, 37.2);
888 AliMaterial(49, "ALUMINIUM$", 26.98, 13., 2.7, 8.9, 37.2);
889
890 // Iron
891 AliMaterial(10, "IRON$ ", 55.85, 26., 7.87, 1.76, 17.1);
892 AliMaterial(30, "IRON$ ", 55.85, 26., 7.87, 1.76, 17.1);
893 AliMaterial(50, "IRON$ ", 55.85, 26., 7.87, 1.76, 17.1);
c557d782 894 // Copper
895 AliMaterial(17, "COPPER$ ", 63.55, 29., 8.96, 1.43, 15.1);
896 AliMaterial(37, "COPPER$ ", 63.55, 29., 8.96, 1.43, 15.1);
897 AliMaterial(57, "COPPER$ ", 63.55, 29., 8.96, 1.43, 15.1);
898 // Air
fe4da5cc 899 AliMaterial(15, "AIR$ ", 14.61, 7.3, .001205, 30423.24, 67500);
900 AliMaterial(35, "AIR$ ", 14.61, 7.3, .001205, 30423.24, 67500);
901 AliMaterial(55, "AIR$ ", 14.61, 7.3, .001205, 30423.24, 67500);
902
903 // Vacuum
904 AliMaterial(16, "VACUUM$ ", 1e-16, 1e-16, 1e-16, 1e16, 1e16);
905 AliMaterial(36, "VACUUM$ ", 1e-16, 1e-16, 1e-16, 1e16, 1e16);
906 AliMaterial(56, "VACUUM$ ", 1e-16, 1e-16, 1e-16, 1e16, 1e16);
907
908 // stainless Steel
909 AliMixture(19, "STAINLESS STEEL$", asteel, zsteel, 7.88, 4, wsteel);
910 AliMixture(39, "STAINLESS STEEL$", asteel, zsteel, 7.88, 4, wsteel);
911 AliMixture(59, "STAINLESS STEEL$", asteel, zsteel, 7.88, 4, wsteel);
912
913 // Coil
914 AliMixture(14, "Al$", acoil, zcoil, 2.122, 3, wcoil);
915 AliMixture(34, "Al$", acoil, zcoil, 2.122, 3, wcoil);
916 AliMixture(54, "Al$", acoil, zcoil, 2.122, 3, wcoil);
1b206195 917
918 //RESIN
919 AliMixture(13, "RESIN$", aresi, zresi, 1.05, 3, wresi);
920 AliMixture(33, "RESIN$", aresi, zresi, 1.05, 3, wresi);
921 AliMixture(53, "RESIN$", aresi, zresi, 1.05, 3, wresi);
922
923 //G10
924 AliMixture(11, "G10$", aG10, zG10, 1.7, 5, wG10);
925 AliMixture(31, "G10$", aG10, zG10, 1.7, 5, wG10);
926 AliMixture(51, "G10$", aG10, zG10, 1.7, 5, wG10);
927
928 //Aluminium Conductor
929 AliMixture(12, "AlCond$", aAlCon, zAlCon, 1.3506, 2, wAlCon);
930 AliMixture(32, "AlCond$", aAlCon, zAlCon, 1.3506, 2, wAlCon);
931 AliMixture(52, "AlCond$", aAlCon, zAlCon, 1.3506, 2, wAlCon);
932
fe4da5cc 933 // ****************
934 // Defines tracking media parameters.
935 // Les valeurs sont commentees pour laisser le defaut
936 // a GEANT (version 3-21, page CONS200), f.m.
937 epsil = .001; // Tracking precision,
938 stemax = -1.; // Maximum displacement for multiple scat
939 tmaxfd = -20.; // Maximum angle due to field deflection
940 deemax = -.3; // Maximum fractional energy loss, DLS
941 stmin = -.8;
942 // ***************
943
944 // Aluminum
b43eb0dc 945 AliMedium(9, "ALU_C0 ", 9, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
946 AliMedium(29, "ALU_C1 ", 29, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
947 AliMedium(49, "ALU_C2 ", 49, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
fe4da5cc 948
949 // Iron
b43eb0dc 950 AliMedium(10, "FE_C0 ", 10, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
951 AliMedium(30, "FE_C1 ", 30, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
952 AliMedium(50, "FE_C2 ", 50, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
fe4da5cc 953
954 // Air
b43eb0dc 955 AliMedium(15, "AIR_C0 ", 15, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
956 AliMedium(35, "AIR_C1 ", 35, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
957 AliMedium(55, "AIR_C2 ", 55, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
fe4da5cc 958
959 // Vacuum
b43eb0dc 960 AliMedium(16, "VA_C0 ", 16, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
961 AliMedium(36, "VA_C1 ", 36, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
962 AliMedium(56, "VA_C2 ", 56, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
fe4da5cc 963
964 // Steel
b43eb0dc 965 AliMedium(19, "ST_C0 ", 19, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
966 AliMedium(39, "ST_C1 ", 39, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
967 AliMedium(59, "ST_C3 ", 59, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
fe4da5cc 968
969 // Coil
b43eb0dc 970 AliMedium(14, "Coil_C1 ", 14, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
971 AliMedium(34, "Coil_C2 ", 34, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
972 AliMedium(54, "Coil_C3 ", 54, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
1b206195 973
974 // Resin
b43eb0dc 975 AliMedium(13, "RESIN_C0 ", 13, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
976 AliMedium(33, "RESIN_C1 ", 33, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
977 AliMedium(53, "RESIN_C2 ", 53, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
1b206195 978
979 // G10
b43eb0dc 980 AliMedium(11, "G10_C0 ", 11, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
981 AliMedium(31, "G10_C1 ", 31, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
982 AliMedium(51, "G10_C2 ", 51, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
1b206195 983
984 //Aluminium Contactor
b43eb0dc 985 AliMedium(12, "AlCond_C0 ", 12, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
986 AliMedium(32, "AlCond_C1 ", 32, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
987 AliMedium(52, "AlCond_C2 ", 52, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
c557d782 988
989 //
990 // Copper
991 AliMedium(17, "Cu_C0 ", 17, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
992 AliMedium(37, "Cu_C1 ", 37, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
993 AliMedium(57, "Cu_C2 ", 57, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
994
fe4da5cc 995}
996
1b206195 997
998
999
1000
1001