STRUCT/AliPIPEv3.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.10  2002/10/14 14:57:39  hristov
  19 Merging the VirtualMC branch to the main development branch (HEAD)
  20
  21 Revision 1.8.12.1  2002/07/24 10:08:27  alibrary
  22 Updating VirtualMC
  23
  24 Revision 1.9  2002/07/23 13:08:46  morsch
  25 All volume names start with "Q".
  26
  27 Revision 1.8  2000/10/02 21:28:15  fca
  28 Removal of useless dependecies via forward declarations
  29
  30 Revision 1.7  2000/06/11 12:38:00  morsch
  31 Coding rule violations corrected
  32
  33 Revision 1.6  1999/10/06 19:57:07  fca
  34 Correct materials in pipe
  35
  36 Revision 1.5  1999/09/29 09:24:30  fca
  37 Introduction of the Copyright and cvs Log
  38
  39 */
  40
  41 ///////////////////////////////////////////////////////////////////////////////
  42 //                                                                           //
  43 //  Beam pipe class                                                          //
  44 //                                                                           //
  45 //Begin_Html
  46 /*
  47 <img src="picts/AliPIPEClass.gif">
  48 */
  49 //End_Html
  50 //                                                                           //
  51 ///////////////////////////////////////////////////////////////////////////////
  52
  53 #include "AliPIPEv3.h"
  54 #include "AliRun.h"
  55 #include "AliMagF.h"
  56
  57 ClassImp(AliPIPEv3)
  58
  59 //_____________________________________________________________________________
  60 AliPIPEv3::AliPIPEv3()
  61 {
  62   //
  63   // Default constructor for beam pipe
  64   //
  65 }
  66
  67 //_____________________________________________________________________________
  68 AliPIPEv3::AliPIPEv3(const char *name, const char *title)
  69        : AliPIPE(name,title)
  70 {
  71   //
  72   // Standard constructor for beam pipe
  73   //
  74 }
  75
  76 //_____________________________________________________________________________
  77 void AliPIPEv3::CreateGeometry()
  78 {
  79   //
  80   // Create Beam Pipe geometry
  81   //
  82   //Begin_Html
  83   /*
  84     <img src="picts/AliPIPE.gif">
  85   */
  86   //End_Html
  87   //Begin_Html
  88   /*
  89     <img src="picts/AliPIPETree.gif">
  90   */
  91   //End_Html
  92
  93   Float_t tpar[3], dzmo, zpos;
  94   Float_t bepar[3], alpar[3],sspar[3],flange[3],vacpar[3];
  95   Float_t bellow[3];
  96 //  Float_t undul[3];
  97 //  const Double_t z_flange = 150;
  98 //for undulated structure
  99   char cn18[][5]={"CN01","CN02","CN03","CN04","CN05","CN06","CN07","CN08"};
 100   char cn48[][5]={"CN21","CN22","CN23","CN24","CN25","CN26","CN27","CN28"};
 101 //  char undul[][5]={'BELO','UNDL'};
 102   Float_t zundul;
 103   Float_t rundul;
 104   Float_t pitch;
 105   Float_t thick;
 106
 107
 108   Int_t *idtmed = fIdtmed->GetArray()-1999;
 109 //     the mother of all beam pipes
 110
 111   tpar[0] = 0.;
 112   tpar[1] = 10.;
 113   tpar[2] = 1400. / 2;
 114   dzmo = tpar[2];
 115   gMC->Gsvolu("QQMO", "TUBE", idtmed[2015], tpar, 3);
 116
 117 //      All beam pipe details as per the provisonal drawings given by Lars
 118 //      Leistam on 31.5.99
 119
 120 //     Beryllium  beam pipe, length 56.6 cm, centered at vertex
 121
 122   bepar[0]=0.0;
 123   bepar[1]=3.0;
 124   bepar[2]=28.3;
 125   zpos=0.0;
 126   vacpar[0]=0.0;
 127   vacpar[1]=2.9;
 128   vacpar[2]=bepar[2];
 129   //
 130   gMC->Gsvolu("QQBE", "TUBE", idtmed[2004], bepar, 3);
 131   gMC->Gsvolu("QAC1", "TUBE", idtmed[2015], vacpar, 3);
 132   gMC->Gspos("QAC1", 1, "QQBE", 0., 0., 0., 0, "ONLY");
 133   gMC->Gspos("QQBE", 1, "QQMO", 0., 0., zpos, 0, "ONLY");
 134
 135   // now beam pipes only in negative z-part for use in PMD.
 136
 137   // SS Flange 4 cm thick, 5.8 cm ID, 6.3 cm OD
 138   flange[0]=0.0;
 139   flange[1]=3.15;
 140   flange[2]=2.0;
 141   zpos = zpos -bepar[2] - flange[2];
 142   vacpar[0]=0.0;
 143   vacpar[1]=2.9;
 144   vacpar[2]=flange[2];
 145   //
 146   gMC->Gsvolu("QFL1", "TUBE", idtmed[2018], flange, 3);
 147   gMC->Gsvolu("QAC2", "TUBE", idtmed[2015], vacpar, 3);
 148   gMC->Gspos("QAC2", 1, "QFL1", 0., 0., 0., 0, "ONLY");
 149   gMC->Gspos("QFL1", 1, "QQMO", 0., 0., zpos, 0, "ONLY");
 150
 151   // Aluminium alloy beam pipe, 1mm thick, 230 cm long
 152   alpar[0]=0.0;
 153   alpar[1]=3.0;
 154   alpar[2]=115.;
 155   zpos = zpos - flange[2] - alpar[2];
 156
 157   vacpar[0]=0.0;
 158   vacpar[1]=2.9;
 159   vacpar[2]=alpar[2];
 160   gMC->Gsvolu("QQAL", "TUBE", idtmed[2003], alpar, 3);
 161   gMC->Gsvolu("QAC3", "TUBE", idtmed[2015], vacpar, 3);
 162   gMC->Gspos("QAC3", 1, "QQAL", 0., 0., 0., 0, "ONLY");
 163   gMC->Gspos("QQAL", 1, "QQMO", 0., 0., zpos, 0, "ONLY");
 164
 165
 166   // SS tube 2.0 cm long, 0.8 mm thick, 5.96 cm OD
 167
 168   sspar[0]=0.0;
 169   sspar[1]=2.98;
 170   sspar[2]=1.0;
 171   zpos = zpos - alpar[2] - sspar[2];
 172
 173   vacpar[0]=0.0;
 174   vacpar[1]=2.9;
 175   vacpar[2]=sspar[2];
 176   gMC->Gsvolu("QSS1", "TUBE", idtmed[2018], sspar, 3);
 177   gMC->Gsvolu("QAC4", "TUBE", idtmed[2015], vacpar, 3);
 178   gMC->Gspos("QAC4", 1, "QSS1", 0., 0., 0., 0, "ONLY");
 179   gMC->Gspos("QSS1", 1, "QQMO", 0., 0., zpos, 0, "ONLY");
 180
 181
 182  // SS Flange 3 cm thick 7.4 cm OD, 5.8 cm ID
 183
 184   flange[0]=0.0;
 185   flange[1]=3.7;
 186   flange[2]=1.5;
 187   zpos = zpos - sspar[2] - flange[2];
 188
 189   vacpar[0]=0.0;
 190   vacpar[1]=2.9;
 191   vacpar[2]=flange[2];
 192   gMC->Gsvolu("QFL2", "TUBE", idtmed[2018], flange, 3);
 193   gMC->Gsvolu("QAC5", "TUBE", idtmed[2015], vacpar, 3);
 194   gMC->Gspos("QAC5", 1, "QFL2", 0., 0., 0., 0, "ONLY");
 195   gMC->Gspos("QFL2", 1, "QQMO", 0., 0., zpos, 0, "ONLY");
 196
 197
 198   // SS tube 4.0 cm long, 0.8 mm thick, 5.96 cm OD
 199
 200   sspar[0]=0.0;
 201   sspar[1]=2.98;
 202   sspar[2]=2.0;
 203   zpos = zpos - flange[2] - sspar[2];
 204
 205   vacpar[0]=0.0;
 206   vacpar[1]=2.9;
 207   vacpar[2]=sspar[2];
 208   gMC->Gsvolu("QSS2", "TUBE", idtmed[2018], sspar, 3);
 209   gMC->Gsvolu("QAC6", "TUBE", idtmed[2015], vacpar, 3);
 210   gMC->Gspos("QAC6", 1, "QSS2", 0., 0., 0., 0, "ONLY");
 211   gMC->Gspos("QSS2", 1, "QQMO", 0., 0., zpos, 0, "ONLY");
 212
 213
 214   // *************
 215   // SS Bellow 8.4 cm long, 6.5 cm ID, 7.5 cm OD
 216   // 0.8 mm thick material, 0.3 cm pitch.
 217   // zundul=4.2, rundul=6.5, thick=0.08
 218   // **************
 219   pitch=0.3;
 220   thick=0.08;
 221   zundul=4.2;
 222   rundul=6.5;
 223   Undulation("BELO",pitch,thick,zundul,rundul,cn18);
 224 //
 225   bellow[2]=zundul;
 226   zpos = zpos - sspar[2] - bellow[2];
 227   gMC->Gspos("BELO", 1, "QQMO", 0., 0., zpos, 0, "ONLY");
 228
 229   // SS tube 20.0 cm long, 0.8 mm thick, 5.96 cm OD
 230
 231   sspar[0]=0.0;
 232   sspar[1]=2.98;
 233   sspar[2]=10.0;
 234   zpos = zpos - bellow[2] - sspar[2];
 235
 236   vacpar[0]=0.0;
 237   vacpar[1]=2.9;
 238   vacpar[2]=sspar[2];
 239   gMC->Gsvolu("QSS3", "TUBE", idtmed[2018], sspar, 3);
 240   gMC->Gsvolu("QAC7", "TUBE", idtmed[2015], vacpar, 3);
 241   gMC->Gspos("QAC7", 1, "QSS3", 0., 0., 0., 0, "ONLY");
 242   gMC->Gspos("QSS3", 1, "QQMO", 0., 0., zpos, 0, "ONLY");
 243
 244   // *************
 245   // SS Bellow 8.4 cm long, 6.5 cm ID, 7.5 cm OD
 246   // 0.8 mm thick material, 0.3 cm pitch.
 247   // **************
 248 //
 249   zpos = zpos - sspar[2] - bellow[2];
 250   gMC->Gspos("BELO", 2, "QQMO", 0., 0., zpos, 0, "ONLY");
 251
 252   // SS tube 4.7 cm long, 0.8 mm thick,
 253
 254   sspar[0]=0.0;
 255   sspar[1]=2.98;
 256   sspar[2]=4.7/2.;
 257   zpos = zpos - bellow[2] - sspar[2];
 258
 259   vacpar[0]=0.0;
 260   vacpar[1]=2.9;
 261   vacpar[2]=sspar[2];
 262   gMC->Gsvolu("QSS4", "TUBE", idtmed[2018], sspar, 3);
 263   gMC->Gsvolu("QAC8", "TUBE", idtmed[2015], vacpar, 3);
 264   gMC->Gspos("QAC8", 1, "QSS4", 0., 0., 0., 0, "ONLY");
 265   gMC->Gspos("QSS4", 1, "QQMO", 0., 0., zpos, 0, "ONLY");
 266
 267   // SS Flange 2.2 cm thick, ID=5.8 cm, OD=9.8 cm
 268
 269   flange[0]=0.0;
 270   flange[1]=4.9;
 271   flange[2]=1.1;
 272   zpos = zpos - sspar[2] - flange[2];
 273
 274   vacpar[0]=0.0;
 275   vacpar[1]=2.9;
 276   vacpar[2]=flange[2];
 277   gMC->Gsvolu("QFL3", "TUBE", idtmed[2018], flange, 3);
 278   gMC->Gsvolu("QAC9", "TUBE", idtmed[2015], vacpar, 3);
 279   gMC->Gspos("QAC9", 1, "QFL3", 0., 0., 0., 0, "ONLY");
 280   gMC->Gspos("QFL3", 1, "QQMO", 0., 0., zpos, 0, "ONLY");
 281
 282 //Total of 3150 mm from vertex on the negative side upto this point.
 283
 284 // SS tube 20.0 cm long, 0.15 cm thick, 5.8 cm ID, to support vac. pump
 285
 286   sspar[0]=0.0;
 287   sspar[1]=3.05;
 288   sspar[2]=10.0;
 289   zpos = zpos - flange[2] - sspar[2];
 290
 291   vacpar[0]=0.0;
 292   vacpar[1]=2.9;
 293   vacpar[2]=sspar[2];
 294   gMC->Gsvolu("QSS5", "TUBE", idtmed[2018], sspar, 3);
 295   gMC->Gsvolu("QA10", "TUBE", idtmed[2015], vacpar, 3);
 296   gMC->Gspos("QA10", 1, "QSS5", 0., 0., 0., 0, "ONLY");
 297   gMC->Gspos("QSS5", 1, "QQMO", 0., 0., zpos, 0, "ONLY");
 298
 299 //
 300   // last item, undulated SS beam pipe, pitch=0.25, length= 342.0 cm
 301   // material thickness 0.015 cm, ID=6.0 cm,
 302   // zundul=171.0, thick=0.015, rundul=3.0
 303   pitch=0.25;
 304   thick=0.015;
 305   zundul=171;
 306   rundul=3.0;
 307   Undulation("UNDL",pitch,thick,zundul,rundul,cn48);
 308   //
 309   zpos = zpos - sspar[2] - zundul;
 310   gMC->Gspos("UNDL", 1, "QQMO", 0., 0., zpos, 0, "ONLY");
 311 //
 312   gMC->Gspos("QQMO", 1, "ALIC", 0., 0., 0.1, 0, "ONLY");
 313
 314 //      total of 6770 mm length upto this point, end of undulated beam
 315 //      pipe section.
 316
 317 //      SS flange 22*2 mm thick
 318
 319
 320   flange[0]=0.0;
 321   flange[1]=6.3;
 322   flange[2]=2.2;
 323   zpos = zpos  - zundul - flange[2];
 324
 325   vacpar[0]=0.0;
 326   vacpar[1]=2.9;
 327   vacpar[2]=flange[2];
 328   gMC->Gsvolu("QFL4", "TUBE", idtmed[2018], flange, 3);
 329   gMC->Gsvolu("QC11", "TUBE", idtmed[2015], vacpar, 3);
 330   gMC->Gspos("QC11", 1, "QFL4", 0., 0., 0., 0, "ONLY");
 331   gMC->Gspos("QFL4", 1, "QQMO", 0., 0., zpos, 0, "ONLY");
 332
 333 }
 334
 335 //_____________________________________________________________________________
 336 void AliPIPEv3::DrawModule()
 337 {
 338   //
 339   // Draw a shaded view of the Beam Pipe
 340   //
 341
 342   // Set everything unseen
 343   gMC->Gsatt("*", "seen", -1);
 344   //
 345   // Set ALIC mother transparent
 346   gMC->Gsatt("ALIC","SEEN",0);
 347   //
 348   // Set the volumes visible
 349   gMC->Gsatt("QQMO","seen",1);
 350   gMC->Gsatt("QQBE","seen",1);
 351   gMC->Gsatt("QFL1","seen",1);
 352   gMC->Gsatt("QQAL","seen",1);
 353   gMC->Gsatt("QSS1","seen",1);
 354   gMC->Gsatt("QFL2","seen",1);
 355   gMC->Gsatt("QSS2","seen",1);
 356   gMC->Gsatt("QSS3","seen",1);
 357   gMC->Gsatt("QSS4","seen",1);
 358   gMC->Gsatt("QFL3","seen",1);
 359   gMC->Gsatt("QSS5","seen",1);
 360   gMC->Gsatt("BELO","seen",1);
 361   gMC->Gsatt("UNDL","seen",1);
 362   //
 363   gMC->Gdopt("hide", "on");
 364   gMC->Gdopt("shad", "on");
 365   gMC->Gsatt("*", "fill", 7);
 366   gMC->SetClipBox(".");
 367   gMC->SetClipBox("*", 0, 3000, -3000, 3000, -6000, 6000);
 368   gMC->DefaultRange();
 369   gMC->Gdraw("alic", 40, 30, 0, 3, 5, .04, .04);
 370   gMC->Gdhead(1111, "Beam Pipe");
 371   gMC->Gdman(16, 6, "MAN");
 372   gMC->Gdopt("hide","off");
 373 }
 374
 375 //_____________________________________________________________________________
 376 void AliPIPEv3::CreateMaterials()
 377 {
 378   //
 379   // Create materials for beam pipe
 380   //
 381
 382   Int_t   isxfld = gAlice->Field()->Integ();
 383   Float_t sxmgmx = gAlice->Field()->Max();
 384
 385   Float_t asteel[4] = { 55.847,51.9961,58.6934,28.0855 };
 386   Float_t zsteel[4] = { 26.,24.,28.,14. };
 387   Float_t wsteel[4] = { .715,.18,.1,.005 };
 388
 389   Float_t epsil, stmin, tmaxfd, deemax, stemax;
 390
 391   //     STEEL
 392
 393
 394   // --- Define the various materials for GEANT ---
 395   AliMaterial(5, "BERILLIUM$", 9.01, 4., 1.848, 35.3, 36.7);
 396   AliMaterial(4, "ALUMINIUM$", 26.98, 13., 2.7, 8.9, 18.5);
 397   AliMaterial(16, "VACUUM$ ", 1e-16, 1e-16, 1e-16, 1e16, 1e16);
 398   AliMaterial(15, "AIR$      ", 14.61, 7.3, .001205, 30423.24, 67500);
 399   AliMixture(19, "STAINLESS STEEL$", asteel, zsteel, 7.88, 4, wsteel);
 400
 401   // ****************
 402   //     Defines tracking media parameters.
 403   //     Les valeurs sont commentees pour laisser le defaut
 404   //     a GEANT (version 3-21, page CONS200), f.m.
 405   epsil  = .001;  // Tracking precision,
 406   stemax = -1.;   // Maximum displacement for multiple scat
 407   tmaxfd = -20.;  // Maximum angle due to field deflection
 408   deemax = -.3;   // Maximum fractional energy loss, DLS
 409   stmin  = -.8;
 410
 411   //    Air
 412
 413   AliMedium(15, "AIR_L3_US", 15, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
 414
 415   //    Beryllium
 416
 417   AliMedium(5, "BE_L3_US", 5, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
 418
 419
 420     //    Aluminium
 421
 422   AliMedium(4, "AL_L3_US", 4, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
 423
 424   //   Vacuum
 425
 426   AliMedium(16, "VA_L3_US", 16, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
 427
 428   //    Steel
 429
 430   AliMedium(19, "ST_L3_US", 19, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
 431 }
 432 //
 433 void AliPIPEv3::Undulation(char *undul, Float_t pitch, Float_t thick,
 434                         Float_t zundul, Float_t rundul, char (*cone)[5])
 435 {
 436   //
 437   // RUNDUL   : Internal radius of the undulated chamber
 438   // THICK    : material thickness
 439   // PITCH    : one-QUARTER wave of undulation (cm)
 440   // ZUNDUL   : half length (cm)
 441   //
 442   // The undulated structure is desgned as a superposition of eight CONES
 443   // of suitable sizes, where the inner/outer radius of the cone increases,
 444   // then decreases, each half of the wave is assumed to be a semicircle,
 445   // which allows to calculate the thickness and the radii of the cone, by
 446   // dividing the semicircle into 4 parts of equal arc length.
 447   // Thus apear the constants 0.293 and 0.707.
 448   //
 449
 450   const Float_t kConst1 = .293;
 451   const Float_t kConst2 = .707;
 452
 453   // Local variables
 454   Int_t j, nwave;
 455   Float_t dcone1[5], dcone2[5], dcone3[5], dcone4[5], dcone5[5],
 456     dcone6[5], dcone7[5], dcone8[5];
 457   Float_t xc, yc, zc, dundul[3];
 458   Int_t *idtmed = fIdtmed->GetArray()-1999;
 459
 460   // Function Body
 461
 462   dcone1[0] = kConst1 * pitch / 2;
 463   dcone1[1] = rundul;
 464   dcone1[2] = dcone1[1] + thick;
 465   dcone1[3] = dcone1[1] + kConst2 * pitch;
 466   dcone1[4] = dcone1[3] + thick;
 467
 468   dcone2[0] = kConst2 * pitch / 2;
 469   dcone2[1] = dcone1[3];
 470   dcone2[2] = dcone1[4];
 471   dcone2[3] = dcone2[1] + kConst1 * pitch;
 472   dcone2[4] = dcone2[3] + thick;
 473
 474   dcone3[0] = dcone2[0];
 475   dcone3[1] = dcone2[3];
 476   dcone3[2] = dcone2[4];
 477   dcone3[3] = dcone2[1];
 478   dcone3[4] = dcone2[2];
 479
 480   dcone4[0] = dcone1[0];
 481   dcone4[1] = dcone1[3];
 482   dcone4[2] = dcone1[4];
 483   dcone4[3] = dcone1[1];
 484   dcone4[4] = dcone1[2];
 485
 486   dcone5[0] = dcone1[0];
 487   dcone5[1] = dcone1[1] - thick;
 488   dcone5[2] = dcone1[1];
 489   dcone5[3] = dcone5[1] - kConst2 * pitch;
 490   dcone5[4] = dcone5[3] + thick;
 491
 492   dcone6[0] = dcone2[0];
 493   dcone6[1] = dcone5[3];
 494   dcone6[2] = dcone5[4];
 495   dcone6[3] = dcone6[1] - kConst1 * pitch;
 496   dcone6[4] = dcone6[3] + thick;
 497   dcone7[0] = dcone6[0];
 498   dcone7[1] = dcone6[3];
 499   dcone7[2] = dcone6[4];
 500   dcone7[3] = dcone5[3];
 501   dcone7[4] = dcone5[4];
 502
 503   dcone8[0] = dcone5[0];
 504   dcone8[1] = dcone7[3];
 505   dcone8[2] = dcone7[4];
 506   dcone8[3] = dcone5[1];
 507   dcone8[4] = dcone5[2];
 508
 509   gMC->Gsvolu(cone[0], "CONE", idtmed[2018], dcone1, 5);
 510   gMC->Gsvolu(cone[1], "CONE", idtmed[2018], dcone2, 5);
 511   gMC->Gsvolu(cone[2], "CONE", idtmed[2018], dcone3, 5);
 512   gMC->Gsvolu(cone[3], "CONE", idtmed[2018], dcone4, 5);
 513   gMC->Gsvolu(cone[4], "CONE", idtmed[2018], dcone5, 5);
 514   gMC->Gsvolu(cone[5], "CONE", idtmed[2018], dcone6, 5);
 515   gMC->Gsvolu(cone[6], "CONE", idtmed[2018], dcone7, 5);
 516   gMC->Gsvolu(cone[7], "CONE", idtmed[2018], dcone8, 5);
 517   gMC->Gsatt(cone[0], "SEEN", 0);
 518   gMC->Gsatt(cone[1], "SEEN", 0);
 519   gMC->Gsatt(cone[2], "SEEN", 0);
 520   gMC->Gsatt(cone[3], "SEEN", 0);
 521   gMC->Gsatt(cone[4], "SEEN", 0);
 522   gMC->Gsatt(cone[5], "SEEN", 0);
 523   gMC->Gsatt(cone[6], "SEEN", 0);
 524   gMC->Gsatt(cone[7], "SEEN", 0);
 525
 526   // DEFINE AN IMAGINARY TUBE VOLUME FOR UNDULATED CHAMBER, FILL WITH VACUUM
 527
 528   nwave = Int_t (zundul / (pitch * 2) + .1);
 529   dundul[2] = pitch * 2 * nwave;
 530   dundul[1] = rundul + pitch + thick * 2;
 531   //
 532   dundul[0] = 1e-4;
 533   gMC->Gsvolu(undul, "TUBE", idtmed[2015], dundul, 3);
 534
 535   xc = 0;
 536   yc = 0;
 537   zc = -dundul[2] + dcone1[0];
 538   for (j = 1; j <= nwave; ++j) {
 539     gMC->Gspos(cone[0], j, undul, xc, yc, zc, 0, "ONLY");
 540     zc = zc + dcone1[0] + dcone2[0];
 541     gMC->Gspos(cone[1], j, undul, xc, yc, zc, 0, "ONLY");
 542     zc = zc + dcone2[0] + dcone3[0];
 543     gMC->Gspos(cone[2], j, undul, xc, yc, zc, 0, "ONLY");
 544     zc = zc + dcone3[0] + dcone4[0];
 545     gMC->Gspos(cone[3], j, undul, xc, yc, zc, 0, "ONLY");
 546     zc = zc + dcone4[0] + dcone5[0];
 547     gMC->Gspos(cone[4], j, undul, xc, yc, zc, 0, "ONLY");
 548     zc = zc + dcone5[0] + dcone6[0];
 549     gMC->Gspos(cone[5], j, undul, xc, yc, zc, 0, "ONLY");
 550     zc = zc + dcone6[0] + dcone7[0];
 551     gMC->Gspos(cone[6], j, undul, xc, yc, zc, 0, "ONLY");
 552     zc = zc + dcone7[0] + dcone8[0];
 553     gMC->Gspos(cone[7], j, undul, xc, yc, zc, 0, "ONLY");
 554     zc = zc + dcone8[0] + dcone1[0];
 555   }
 556 }