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 /*
  17 $Log$
  18 Revision 1.18  2001/06/21 14:16:43  morsch
  19 Use absolute coordinates for mother volume PCON.
  20
  21 Revision 1.17  2001/06/21 11:59:25  morsch
  22 Some more details in compensator geometry.
  23
  24 Revision 1.16  2001/06/20 16:07:08  morsch
  25 Compensator dipole MBWMD (MCB@SPS) added.
  26
  27 Revision 1.15  2001/03/16 15:34:37  morsch
  28 Mothervolume defined MANY because overlap with station 3 mothervolume not avoidable (A. de Falco)
  29
  30 Revision 1.14  2000/12/21 16:37:23  morsch
  31 Use Al for coil and cable material. The materials used before cause the dipole to
  32 have hydrogene on the outer surface leading to unrealistic gamma rates due to
  33 n-capture.
  34
  35 Revision 1.13  2000/10/02 21:28:15  fca
  36 Removal of useless dependecies via forward declarations
  37
  38 Revision 1.12  2000/06/20 10:53:01  morsch
  39 Volume placed outside mother volume (DDIP) corrected (Galina Chabratova)
  40
  41 Revision 1.11  2000/06/11 12:33:46  morsch
  42 Coding rule violations corrected
  43
  44 Revision 1.10  2000/06/09 19:32:56  morsch
  45 New detailed and corrected version from Galina Chabratova
  46
  47 Revision 1.9  2000/04/27 09:29:53  fca
  48 Reverting to version 1.6.2
  49
  50 Revision 1.6.2.1  1999/12/03 16:38:51  fca
  51 Correct overlap in magnet
  52
  53 Revision 1.6  1999/09/29 09:24:30  fca
  54 Introduction of the Copyright and cvs Log
  55
  56 */
  57
  58 ///////////////////////////////////////////////////////////////////////////////
  59 //                                                                           //
  60 //  Magnetic Dipole version 1                                                //
  61 //                                                                           //
  62 //Begin_Html
  63 /*
  64 <img src="picts/AliDIPOv2Class.gif">
  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 //                                                                           //
  76
  77 #include "AliDIPOv2.h"
  78 #include "AliRun.h"
  79 #include "AliMagF.h"
  80 #include "AliConst.h"
  81
  82 ClassImp(AliDIPOv2)
  83
  84 //_____________________________________________________________________________
  85 AliDIPOv2::AliDIPOv2()
  86 {
  87   //
  88   // Last design of magnetic dipole version 2
  89   //
  90 }
  91
  92 //_____________________________________________________________________________
  93 AliDIPOv2::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 }
 102
 103 void AliDIPOv2::CreateGeometry()
 104 {
 105   //
 106   // Creation of the geometry of the magnetic DIPOLE version 2
 107   //
 108
 109     CreateSpectrometerDipole();
 110     CreateCompensatorDipole();
 111 }
 112
 113
 114 //_____________________________________________________________________________
 115 void AliDIPOv2::CreateSpectrometerDipole()
 116 {
 117   //
 118   // Creation of the geometry of the magnetic DIPOLE version 2
 119   //
 120
 121   Float_t cpar[5], tpar[18], ypar[12];
 122   Float_t dz, dx, dy;
 123   Int_t idrotm[1899];
 124   Float_t accMax, the1, phi1, the2, phi2, the3, phi3;
 125
 126   Int_t *idtmed = fIdtmed->GetArray()-1799;
 127
 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
 134   accMax = 9.;   // ANGLE POLAIRE MAXIMUM
 135
 136   //       DIPOLE MAGNET
 137   const Float_t kZDipole = 975;
 138
 139   tpar[0] = 0.;
 140   tpar[1] = 360.;
 141   tpar[2] = 5.;
 142   //
 143   tpar[3] = -250.55+kZDipole;
 144   tpar[4] = 30.5;
 145   tpar[5] = 527.34;
 146   //
 147   tpar[6] = -160.7+kZDipole;
 148   tpar[7] = 30.5;
 149   tpar[8] = 527.34;
 150   //
 151   tpar[9] = 30.+kZDipole;
 152   tpar[10] = 30.5;
 153   tpar[11] = 527.34;
 154
 155   tpar[12] = 150.8+kZDipole;
 156   tpar[13] = 193.3;
 157   tpar[14] = 527.34;
 158   //
 159   tpar[15] = 250.55+kZDipole;
 160   tpar[16] = 193.3;
 161   tpar[17] = 527.34;
 162
 163
 164   gMC->Gsvolu("DDIP", "PCON", idtmed[1814], tpar, 18);
 165   //       COILS
 166   // air - m.f.
 167   cpar[0] = 207.;
 168   cpar[1] = 274.;
 169   cpar[2] = 37.65;
 170   cpar[3] = 119.;
 171   cpar[4] = 241.;
 172   //   coil - high cuts
 173   gMC->Gsvolu("DC1 ", "TUBS", idtmed[kCoil+40], cpar, 5);
 174   cpar[3] = -61.;
 175   cpar[4] = 61.;
 176   gMC->Gsvolu("DC2 ", "TUBS", idtmed[kCoil+40], cpar, 5);
 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
 185   gMC->Gsvolu("DC3 ", "TUBS", idtmed[kCoil], cpar, 5);
 186   cpar[0] = 207.;
 187   cpar[1] = 217;
 188   cpar[3] = -61.;
 189   cpar[4] = 61.;
 190   gMC->Gsvolu("DC4 ", "TUBS", idtmed[kCoil], cpar, 5);
 191
 192   gMC->Gspos("DC3 ", 1, "DC1 ", 0., 0., 0., 0, "ONLY");
 193   gMC->Gspos("DC4 ", 1, "DC2 ", 0., 0., 0., 0, "ONLY");
 194
 195 //  dz =  37.65 - 243.55
 196   dz = -205.9-2.45;
 197   dx = 5.;
 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");
 202   the1 = 180.;
 203   phi1 = 0.;
 204   the2 = 90.;
 205   phi2 = 151.;
 206   the3 = 90.;
 207   phi3 = 61.;
 208   AliMatrix(idrotm[1800], the1, phi1, the2, phi2, the3, phi3);
 209   phi2 = 29.;  //90-61
 210   the3 = -90.;
 211   phi3 = -61.;
 212   AliMatrix(idrotm[1801], the1, phi1, the2, phi2, the3, phi3);
 213   the1 = 0.;
 214   phi1 = 0.;
 215   the2 = 90.;
 216   phi2 = 151.;
 217   the3 = 90.;
 218   phi3 = 61.;
 219   AliMatrix(idrotm[1802], the1, phi1, the2, phi2, the3, phi3);
 220   phi2 = 29.;
 221   the3 = -90.;
 222   phi3 = -61.;
 223   AliMatrix(idrotm[1803], the1, phi1, the2, phi2, the3, phi3);
 224
 225   cpar[0] = 25.;
 226   cpar[1] = 100.3; //25+75.3
 227   cpar[2] = 33.5;
 228   cpar[3] = 270.;
 229   cpar[4] = 360.;
 230 //*  coil high cuts
 231   gMC->Gsvolu("DC11", "TUBS", idtmed[kCoil+40], cpar, 5);
 232
 233   dx = TMath::Sin(30.5*kDegrad) * -(207.+33.5)+5./TMath::Sin(30.5*kDegrad);
 234   dy = TMath::Cos(30.5*kDegrad) * -(207.+33.5);
 235   dz = cpar[1] - 243.55-2.45;
 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");
 240
 241
 242
 243 //* ... higher cuts
 244   cpar[0] = 25.;
 245   cpar[1] = 100.3; //25+75.3
 246   cpar[2] = 33.5;
 247   cpar[3] = 0.;
 248   cpar[4] = 90.;
 249 //*  coil high cuts
 250   gMC->Gsvolu("DC12", "TUBS", idtmed[kCoil+40], cpar, 5);
 251
 252   dx = TMath::Sin(30.5*kDegrad) * -(207.+33.5)+5./TMath::Sin(30.5*kDegrad);
 253   dy = TMath::Cos(30.5*kDegrad) *(207.+33.5);
 254   dz = cpar[1] - 243.55-2.45;
 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");
 259
 260   the1 = 90.;
 261   phi1 = 61.;
 262   the2 = 90.;
 263   phi2 = 151.;
 264   the3 = 0.;
 265   phi3 = 0.;
 266   AliMatrix(idrotm[1804], the1, phi1, the2, phi2, the3, phi3);
 267   the1 = 90.;
 268   phi1 = -61.;
 269   the2 = 90.;
 270   phi2 = -151.;
 271   AliMatrix(idrotm[1805], the1, phi1, the2, phi2, the3, phi3);
 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;
 288   gMC->Gsvolu("DL1 ", "BOX ", idtmed[kCoil+40], tpar, 3);
 289
 290 // coil - low cuts
 291
 292   tpar[0] = 5.;
 293   dx = 37.65  - 5.;
 294   gMC->Gsvolu("DL2 ", "BOX ", idtmed[kCoil], tpar, 3);
 295   gMC->Gspos("DL2 ", 1, "DL1 ", dx, 0., 0., 0, "ONLY");
 296
 297   dx =-53.62;
 298   dy =-241.26819;
 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");
 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.;
 320   gMC->Gspos("DCO1", 1, "DDIP", dx, 0, dz+kZDipole, 0, "ONLY");
 321   dz = 243.55+4.5+1.5+1.;
 322   gMC->Gspos("DCO1", 2, "DDIP", dx, 0, dz+kZDipole, 0, "ONLY");
 323
 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.;
 335   gMC->Gspos("DCO2", 1, "DDIP", dx, 0, dz+kZDipole, 0, "ONLY");
 336   dz = 243.55+4.5+1.5+1.;
 337   gMC->Gspos("DCO2", 2, "DDIP", dx, 0, dz+kZDipole, 0, "ONLY");
 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;
 352   gMC->Gspos("DCO3", 1, "DDIP", dx,  0, dz+kZDipole, 0, "ONLY");
 353   dz = 243.55+4.5+0.75;
 354   gMC->Gspos("DCO3", 2, "DDIP", dx,  0, dz+kZDipole, 0, "ONLY");
 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;
 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");
 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;
 380   gMC->Gspos("DCO5", 1, "DDIP", dx,  0, dz+kZDipole, 0, "ONLY");
 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;
 390   dz = 243.55+2.25;
 391   gMC->Gspos("DCO6", 1, "DDIP", dx,  0, dz+kZDipole, 0, "ONLY");
 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.;
 405   gMC->Gspos("DCO7", 1, "DDIP", dx, 0, dz+kZDipole, 0, "ONLY");
 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.;
 419   gMC->Gspos("DCO8", 1, "DDIP", dx, 0, dz+kZDipole, 0, "ONLY");
 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.;
 433   gMC->Gspos("DCO9", 1, "DDIP", dx, 0, dz+kZDipole, 0, "ONLY");
 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.;
 447   gMC->Gspos("DCOA", 1, "DDIP", dx, 0, dz+kZDipole, 0, "ONLY");
 448
 449
 450   // Sides steel planes
 451
 452   cpar[0] = 207. - 1.5 -2.;
 453   cpar[1] = 207. - 1.5;
 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;
 467   gMC->Gspos("DCOB", 1, "DDIP", dx,  0, dz+kZDipole, 0, "ONLY");
 468   gMC->Gspos("DCOC", 1, "DDIP", dx,  0, dz+kZDipole, 0, "ONLY");
 469
 470   // 9.06.2000
 471
 472   cpar[0] = 207. - 1.5 -2.;
 473   cpar[1] = 207. - 1.5;
 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;
 487   gMC->Gspos("DCOD", 1, "DDIP", dx,  0, dz+kZDipole, 0, "ONLY");
 488   gMC->Gspos("DCOE", 1, "DDIP", dx,  0, dz+kZDipole, 0, "ONLY");
 489
 490
 491   // Top and bottom resin  planes
 492
 493   cpar[0] = 207. - 1.5;
 494   cpar[1] = 207.;
 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;
 509   gMC->Gspos("DCOF", 1, "DDIP", dx,  0, dz+kZDipole, 0, "ONLY");
 510   gMC->Gspos("DCOG", 1, "DDIP", dx,  0, dz+kZDipole, 0, "ONLY");
 511
 512   // 9.06.2000
 513   cpar[0] = 207. - 1.5;
 514   cpar[1] = 207.;
 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;
 530   gMC->Gspos("DCOH", 1, "DDIP", dx,  0, dz+kZDipole, 0, "ONLY");
 531   gMC->Gspos("DCOI", 1, "DDIP", dx,  0, dz+kZDipole, 0, "ONLY");
 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
 542   gMC->Gsvolu("DCOJ", "TUBS", idtmed[kCable], cpar, 5);
 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;
 549   gMC->Gspos("DCOJ", 1, "DDIP", dx,  0, dz+kZDipole, 0, "ONLY");
 550
 551   dz = 243.55 - 5.05/2;
 552   gMC->Gspos("DCOJ", 2, "DDIP", dx,  0, dz+kZDipole, 0, "ONLY");
 553
 554   // 9.06.2000
 555
 556   cpar[3] = 180.-24.;
 557   cpar[4] = 180.+24.;
 558
 559    gMC->Gsvolu("DCOK", "TUBS", idtmed[kCable], cpar, 5);
 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;
 566   gMC->Gspos("DCOK", 1, "DDIP", dx,  0, dz+kZDipole, 0, "ONLY");
 567
 568   dz = 243.55 - 5.05/2;
 569   gMC->Gspos("DCOK", 2, "DDIP", dx,  0, dz+kZDipole, 0, "ONLY");
 570
 571
 572   //   YOKE
 573
 574 // Top and bottom blocks
 575   ypar[0] = 298.1;
 576   ypar[1] = 69.5;
 577   ypar[2] = 155.75;
 578
 579 // iron- high cuts
 580   gMC->Gsvolu("DY1 ", "BOX ", idtmed[1858], ypar, 3);
 581   ypar[0] = 144.+10.;
 582   ypar[1] = 193.3+10.;
 583   ypar[2] = 5.;
 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;
 592   gMC->Gspos("DY1 ", 1, "DDIP", 0.,  dy, -dz+kZDipole, 0, "ONLY");
 593
 594   the1 = 270.;
 595   phi1 = 0.;
 596   the2 = 270.;
 597   phi2 = 90.;
 598   the3 = 0.;
 599   phi3 = 0.;
 600   AliMatrix(idrotm[1808], the1, phi1, the2, phi2, the3, phi3);
 601   gMC->Gspos("DY1 ", 2, "DDIP", 0., -dy, -dz+kZDipole, idrotm[1808] , "ONLY");
 602
 603 // side walls
 604   //  ypar[0] = 579./2.;
 605   ypar[0] = 296.;
 606   ypar[1] = 0.;
 607   ypar[2] = 0.;
 608   ypar[3] = 155.75;
 609   ypar[4] = 47.9;
 610   ypar[5] = 72.55;
 611   ypar[6] = 4.3058039629;
 612   // z+
 613   ypar[7] = 155.75;
 614   ypar[8] = 47.9;
 615   ypar[9] = 72.55;
 616   ypar[10] = 4.3058039629;
 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.;
 647   phi1 = 0.;
 648   the2 = 180.;
 649   phi2 = 0.;
 650   the3 = 90.;
 651   phi3 = 90.;
 652   AliMatrix(idrotm[1810], the1, phi1, the2, phi2, the3, phi3);
 653
 654   dx = 228.875;
 655   dz = - 4.95;
 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");
 659
 660   gMC->Gspos("DDIP", 1, "ALIC", 0., 0., 0., 0, "ONLY");
 661
 662   gMC->Gsatt("DDIP", "SEEN", 0);
 663 //  gMC->Gsatt("DC21", "SEEN", 0);
 664 //  gMC->Gsatt("DC22", "SEEN", 0);
 665 //  gMC->Gsatt("DC3 ", "SEEN", 0);
 666 //  gMC->Gsatt("DC4 ", "SEEN", 0);
 667 }
 668
 669
 670 void 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);
 682
 683 //
 684 //  Mother volume containing lower coil
 685     pbox[0] = 58.5/2.;
 686     pbox[1] = 30.0;
 687     pbox[2] = 152.5;
 688
 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);
 695 //
 696 // Coil: straight sections, horizontal
 697     pbox[0] =   6.;
 698     pbox[1] =  11.;
 699     pbox[2] = 135.;
 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
 717     pbox[0] =  6.0;
 718     pbox[1] =  9.5;
 719     pbox[2] = 11.0;
 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.;
 731 //    gMC->Gsvolu("DCC1", "TUBS", idtmed[1809], ptubs, 5);
 732     ptubs[0] = 13.;
 733     ptubs[1] = 35.;
 734     ptubs[2] =  6.;
 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);
 752 //
 753 // Assemble all
 754 //
 755     AliMatrix(idrotm[1811], -90., 0., 90., 90.,   0., 0.);
 756     AliMatrix(idrotm[1812],   0., 0., 90., 90.,  90., 0.);
 757     AliMatrix(idrotm[1813], 180., 0., 90., 90.,  90., 0.);
 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");
 765
 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");
 768
 769
 770     gMC->Gspos("DCH2", 1, "DCMU", 2., 6.5, 0., 0, "ONLY");
 771     gMC->Gspos("DCMU", 1, "DCM0", -12., 45., 0., 0, "ONLY");
 772     gMC->Gspos("DCMU", 2, "DCM0",  12., 45., 0., idrotm[1811], "ONLY");
 773
 774 //    gMC->Gspos("DCC2", 1, "DCC1", 0., 0., 0., 0, "ONLY");
 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
 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");
 795
 796
 797     gMC->Gspos("DCBA", 1, "DCM0",  0., -47.5 , 17.5, 0, "ONLY");
 798     gMC->Gspos("DCM0", 1, "ALIC",  0., -6.75, -975., 0, "ONLY");
 799
 800
 801 }
 802
 803 //_____________________________________________________________________________
 804 void AliDIPOv2::DrawModule()
 805 {
 806   //
 807   // Draw a shaded view of the muon absorber
 808   //
 809
 810   // Set everything unseen
 811   gMC->Gsatt("*", "seen", -1);
 812   //
 813   // Set ALIC mother transparent
 814   gMC->Gsatt("ALIC","SEEN",0);
 815   //
 816   // Set the volumes visible
 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);
 832  // gMC->Gsatt("DY4 ","seen",1);
 833  // gMC->Gsatt("DY5 ","seen",1);
 834  // gMC->Gsatt("DY6 ","seen",1);
 835 //  gMC->Gsatt("DY7 ","seen",1);
 836   //
 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");
 846 }
 847
 848
 849
 850
 851 //_____________________________________________________________________________
 852 void AliDIPOv2::CreateMaterials()
 853 {
 854   //
 855   // Create Materials for Magnetic Dipole version 2
 856   //
 857
 858   Int_t isxfld   = gAlice->Field()->Integ();
 859   Float_t sxmgmx = gAlice->Field()->Max();
 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 };
 864
 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 };
 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
 881
 882   Float_t epsil, stmin, deemax, tmaxfd, stemax;
 883
 884   // --- Define the various materials for GEANT ---
 885   //     Aluminum
 886   AliMaterial( 9, "ALUMINIUM$", 26.98, 13., 2.7, 8.9, 37.2);
 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);
 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
 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);
 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
 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
 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);
 948
 949   //    Iron
 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);
 953
 954   //    Air
 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);
 958
 959   //    Vacuum
 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);
 963
 964   //    Steel
 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);
 968
 969   //    Coil
 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);
 973
 974   //    Resin
 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);
 978
 979   //    G10
 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);
 983
 984   //Aluminium Contactor
 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);
 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
 995 }
 996
 997
 998
 999
1000
1001