STRUCT/AliDIPOv2.cxx

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