[u/mrichter/AliRoot.git] / STRUCT / AliPIPEv3.cxx

/**************************************************************************
 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/

/*
$Log$
Revision 1.7  2000/06/11 12:38:00  morsch
Coding rule violations corrected

Revision 1.6  1999/10/06 19:57:07  fca
Correct materials in pipe

Revision 1.5  1999/09/29 09:24:30  fca
Introduction of the Copyright and cvs Log

*/

///////////////////////////////////////////////////////////////////////////////
//                                                                           //
//  Beam pipe class                                                          //
//                                                                           //
//Begin_Html
/*
<img src="picts/AliPIPEClass.gif">
*/
//End_Html
//                                                                           //
///////////////////////////////////////////////////////////////////////////////

#include "AliPIPEv3.h"
#include "AliRun.h"
#include "AliMC.h"
#include "AliMagF.h"
 
ClassImp(AliPIPEv3)
 
//_____________________________________________________________________________
AliPIPEv3::AliPIPEv3()
{
  //
  // Default constructor for beam pipe
  //
}
 
//_____________________________________________________________________________
AliPIPEv3::AliPIPEv3(const char *name, const char *title)
       : AliPIPE(name,title)
{
  //
  // Standard constructor for beam pipe
  //
}
 
//_____________________________________________________________________________
void AliPIPEv3::CreateGeometry()
{
  //
  // Create Beam Pipe geometry
  //
  //Begin_Html
  /*
    <img src="picts/AliPIPE.gif">
  */
  //End_Html
  //Begin_Html
  /*
    <img src="picts/AliPIPETree.gif">
  */
  //End_Html

  Float_t tpar[3], dzmo, zpos;
  Float_t bepar[3], alpar[3],sspar[3],flange[3],vacpar[3];
  Float_t bellow[3];
//  Float_t undul[3];
//  const Double_t z_flange = 150;
//for undulated structure
  char cn18[][5]={"CN01","CN02","CN03","CN04","CN05","CN06","CN07","CN08"};
  char cn48[][5]={"CN21","CN22","CN23","CN24","CN25","CN26","CN27","CN28"};
//  char undul[][5]={'BELO','UNDL'};
  Float_t zundul;
  Float_t rundul;
  Float_t pitch;
  Float_t thick;

  
  Int_t *idtmed = fIdtmed->GetArray()-1999;
//     the mother of all beam pipes

  tpar[0] = 0.;
  tpar[1] = 10.;
  tpar[2] = 1400. / 2;
  dzmo = tpar[2];
  gMC->Gsvolu("QQMO", "TUBE", idtmed[2015], tpar, 3);

//	All beam pipe details as per the provisonal drawings given by Lars
//	Leistam on 31.5.99 
    
//     Beryllium  beam pipe, length 56.6 cm, centered at vertex 
  
  bepar[0]=0.0;
  bepar[1]=3.0;
  bepar[2]=28.3;
  zpos=0.0;
  vacpar[0]=0.0;
  vacpar[1]=2.9;
  vacpar[2]=bepar[2];
  //
  gMC->Gsvolu("QQBE", "TUBE", idtmed[2004], bepar, 3);
  gMC->Gsvolu("VAC1", "TUBE", idtmed[2015], vacpar, 3);
  gMC->Gspos("VAC1", 1, "QQBE", 0., 0., 0., 0, "ONLY");
  gMC->Gspos("QQBE", 1, "QQMO", 0., 0., zpos, 0, "ONLY");
  
  // now beam pipes only in negative z-part for use in PMD.
 
  // SS Flange 4 cm thick, 5.8 cm ID, 6.3 cm OD
  flange[0]=0.0;
  flange[1]=3.15;
  flange[2]=2.0;
  zpos = zpos -bepar[2] - flange[2];
  vacpar[0]=0.0;
  vacpar[1]=2.9;
  vacpar[2]=flange[2];
  //
  gMC->Gsvolu("QFL1", "TUBE", idtmed[2018], flange, 3);
  gMC->Gsvolu("VAC2", "TUBE", idtmed[2015], vacpar, 3);
  gMC->Gspos("VAC2", 1, "QFL1", 0., 0., 0., 0, "ONLY");
  gMC->Gspos("QFL1", 1, "QQMO", 0., 0., zpos, 0, "ONLY");
  
  // Aluminium alloy beam pipe, 1mm thick, 230 cm long
  alpar[0]=0.0;
  alpar[1]=3.0;
  alpar[2]=115.;
  zpos = zpos - flange[2] - alpar[2];

  vacpar[0]=0.0;
  vacpar[1]=2.9;
  vacpar[2]=alpar[2];
  gMC->Gsvolu("QQAL", "TUBE", idtmed[2003], alpar, 3);
  gMC->Gsvolu("VAC3", "TUBE", idtmed[2015], vacpar, 3);
  gMC->Gspos("VAC3", 1, "QQAL", 0., 0., 0., 0, "ONLY");
  gMC->Gspos("QQAL", 1, "QQMO", 0., 0., zpos, 0, "ONLY");

 
  // SS tube 2.0 cm long, 0.8 mm thick, 5.96 cm OD

  sspar[0]=0.0;
  sspar[1]=2.98;
  sspar[2]=1.0;
  zpos = zpos - alpar[2] - sspar[2];

  vacpar[0]=0.0;
  vacpar[1]=2.9;
  vacpar[2]=sspar[2];
  gMC->Gsvolu("QSS1", "TUBE", idtmed[2018], sspar, 3);
  gMC->Gsvolu("VAC4", "TUBE", idtmed[2015], vacpar, 3);
  gMC->Gspos("VAC4", 1, "QSS1", 0., 0., 0., 0, "ONLY");
  gMC->Gspos("QSS1", 1, "QQMO", 0., 0., zpos, 0, "ONLY");


 // SS Flange 3 cm thick 7.4 cm OD, 5.8 cm ID
   
  flange[0]=0.0;
  flange[1]=3.7;
  flange[2]=1.5;
  zpos = zpos - sspar[2] - flange[2];

  vacpar[0]=0.0;
  vacpar[1]=2.9;
  vacpar[2]=flange[2];
  gMC->Gsvolu("QFL2", "TUBE", idtmed[2018], flange, 3);
  gMC->Gsvolu("VAC5", "TUBE", idtmed[2015], vacpar, 3);
  gMC->Gspos("VAC5", 1, "QFL2", 0., 0., 0., 0, "ONLY");
  gMC->Gspos("QFL2", 1, "QQMO", 0., 0., zpos, 0, "ONLY");


  // SS tube 4.0 cm long, 0.8 mm thick, 5.96 cm OD

  sspar[0]=0.0;
  sspar[1]=2.98;
  sspar[2]=2.0;
  zpos = zpos - flange[2] - sspar[2];

  vacpar[0]=0.0;
  vacpar[1]=2.9;
  vacpar[2]=sspar[2];
  gMC->Gsvolu("QSS2", "TUBE", idtmed[2018], sspar, 3);
  gMC->Gsvolu("VAC6", "TUBE", idtmed[2015], vacpar, 3);
  gMC->Gspos("VAC6", 1, "QSS2", 0., 0., 0., 0, "ONLY");
  gMC->Gspos("QSS2", 1, "QQMO", 0., 0., zpos, 0, "ONLY");


  // *************
  // SS Bellow 8.4 cm long, 6.5 cm ID, 7.5 cm OD
  // 0.8 mm thick material, 0.3 cm pitch.
  // zundul=4.2, rundul=6.5, thick=0.08
  // **************
  pitch=0.3;
  thick=0.08;
  zundul=4.2;
  rundul=6.5;
  Undulation("BELO",pitch,thick,zundul,rundul,cn18);
//
  bellow[2]=zundul;
  zpos = zpos - sspar[2] - bellow[2];
  gMC->Gspos("BELO", 1, "QQMO", 0., 0., zpos, 0, "ONLY");

  // SS tube 20.0 cm long, 0.8 mm thick, 5.96 cm OD

  sspar[0]=0.0;
  sspar[1]=2.98;
  sspar[2]=10.0;
  zpos = zpos - bellow[2] - sspar[2];

  vacpar[0]=0.0;
  vacpar[1]=2.9;
  vacpar[2]=sspar[2];
  gMC->Gsvolu("QSS3", "TUBE", idtmed[2018], sspar, 3);
  gMC->Gsvolu("VAC7", "TUBE", idtmed[2015], vacpar, 3);
  gMC->Gspos("VAC7", 1, "QSS3", 0., 0., 0., 0, "ONLY");
  gMC->Gspos("QSS3", 1, "QQMO", 0., 0., zpos, 0, "ONLY");

  // *************
  // SS Bellow 8.4 cm long, 6.5 cm ID, 7.5 cm OD
  // 0.8 mm thick material, 0.3 cm pitch.
  // **************
//  
  zpos = zpos - sspar[2] - bellow[2];
  gMC->Gspos("BELO", 2, "QQMO", 0., 0., zpos, 0, "ONLY");

  // SS tube 4.7 cm long, 0.8 mm thick, 

  sspar[0]=0.0;
  sspar[1]=2.98;
  sspar[2]=4.7/2.;
  zpos = zpos - bellow[2] - sspar[2];

  vacpar[0]=0.0;
  vacpar[1]=2.9;
  vacpar[2]=sspar[2];
  gMC->Gsvolu("QSS4", "TUBE", idtmed[2018], sspar, 3);
  gMC->Gsvolu("VAC8", "TUBE", idtmed[2015], vacpar, 3);
  gMC->Gspos("VAC8", 1, "QSS4", 0., 0., 0., 0, "ONLY");
  gMC->Gspos("QSS4", 1, "QQMO", 0., 0., zpos, 0, "ONLY");

  // SS Flange 2.2 cm thick, ID=5.8 cm, OD=9.8 cm

  flange[0]=0.0;
  flange[1]=4.9;
  flange[2]=1.1;
  zpos = zpos - sspar[2] - flange[2];

  vacpar[0]=0.0;
  vacpar[1]=2.9;
  vacpar[2]=flange[2];
  gMC->Gsvolu("QFL3", "TUBE", idtmed[2018], flange, 3);
  gMC->Gsvolu("VAC9", "TUBE", idtmed[2015], vacpar, 3);
  gMC->Gspos("VAC9", 1, "QFL3", 0., 0., 0., 0, "ONLY");
  gMC->Gspos("QFL3", 1, "QQMO", 0., 0., zpos, 0, "ONLY");

//Total of 3150 mm from vertex on the negative side upto this point.

// SS tube 20.0 cm long, 0.15 cm thick, 5.8 cm ID, to support vac. pump

  sspar[0]=0.0;
  sspar[1]=3.05;
  sspar[2]=10.0;
  zpos = zpos - flange[2] - sspar[2];

  vacpar[0]=0.0;
  vacpar[1]=2.9;
  vacpar[2]=sspar[2];
  gMC->Gsvolu("QSS5", "TUBE", idtmed[2018], sspar, 3);
  gMC->Gsvolu("VA10", "TUBE", idtmed[2015], vacpar, 3);
  gMC->Gspos("VA10", 1, "QSS5", 0., 0., 0., 0, "ONLY");
  gMC->Gspos("QSS5", 1, "QQMO", 0., 0., zpos, 0, "ONLY");

// 
  // last item, undulated SS beam pipe, pitch=0.25, length= 342.0 cm
  // material thickness 0.015 cm, ID=6.0 cm,
  // zundul=171.0, thick=0.015, rundul=3.0
  pitch=0.25;
  thick=0.015;
  zundul=171;
  rundul=3.0;
  Undulation("UNDL",pitch,thick,zundul,rundul,cn48);
  //
  zpos = zpos - sspar[2] - zundul;
  gMC->Gspos("UNDL", 1, "QQMO", 0., 0., zpos, 0, "ONLY");
//
  gMC->Gspos("QQMO", 1, "ALIC", 0., 0., 0.1, 0, "ONLY");

// 	total of 6770 mm length upto this point, end of undulated beam
//	pipe section.

//	SS flange 22*2 mm thick


  flange[0]=0.0;
  flange[1]=6.3;
  flange[2]=2.2;
  zpos = zpos  - zundul - flange[2];

  vacpar[0]=0.0;
  vacpar[1]=2.9;
  vacpar[2]=flange[2];
  gMC->Gsvolu("QFL4", "TUBE", idtmed[2018], flange, 3);
  gMC->Gsvolu("VC11", "TUBE", idtmed[2015], vacpar, 3);
  gMC->Gspos("VC11", 1, "QFL4", 0., 0., 0., 0, "ONLY");
  gMC->Gspos("QFL4", 1, "QQMO", 0., 0., zpos, 0, "ONLY");

}

//_____________________________________________________________________________
void AliPIPEv3::DrawModule()
{  
  //
  // Draw a shaded view of the Beam Pipe
  //

  // Set everything unseen
  gMC->Gsatt("*", "seen", -1);
  // 
  // Set ALIC mother transparent
  gMC->Gsatt("ALIC","SEEN",0);
  //
  // Set the volumes visible
  gMC->Gsatt("QQMO","seen",1);
  gMC->Gsatt("QQBE","seen",1);
  gMC->Gsatt("QFL1","seen",1);
  gMC->Gsatt("QQAL","seen",1);
  gMC->Gsatt("QSS1","seen",1);
  gMC->Gsatt("QFL2","seen",1);
  gMC->Gsatt("QSS2","seen",1);
  gMC->Gsatt("QSS3","seen",1);
  gMC->Gsatt("QSS4","seen",1);
  gMC->Gsatt("QFL3","seen",1);
  gMC->Gsatt("QSS5","seen",1);
  gMC->Gsatt("BELO","seen",1);
  gMC->Gsatt("UNDL","seen",1);
  //
  gMC->Gdopt("hide", "on");
  gMC->Gdopt("shad", "on");
  gMC->Gsatt("*", "fill", 7);
  gMC->SetClipBox(".");
  gMC->SetClipBox("*", 0, 3000, -3000, 3000, -6000, 6000);
  gMC->DefaultRange();
  gMC->Gdraw("alic", 40, 30, 0, 3, 5, .04, .04);
  gMC->Gdhead(1111, "Beam Pipe");
  gMC->Gdman(16, 6, "MAN");
  gMC->Gdopt("hide","off");
}

//_____________________________________________________________________________
void AliPIPEv3::CreateMaterials()
{
  //
  // Create materials for beam pipe
  //

  Int_t   isxfld = gAlice->Field()->Integ();
  Float_t sxmgmx = gAlice->Field()->Max();
  
  Float_t asteel[4] = { 55.847,51.9961,58.6934,28.0855 };
  Float_t zsteel[4] = { 26.,24.,28.,14. };
  Float_t wsteel[4] = { .715,.18,.1,.005 };
  
  Float_t epsil, stmin, tmaxfd, deemax, stemax;
  
  //     STEEL 
  
  
  // --- Define the various materials for GEANT --- 
  AliMaterial(5, "BERILLIUM$", 9.01, 4., 1.848, 35.3, 36.7);
  AliMaterial(4, "ALUMINIUM$", 26.98, 13., 2.7, 8.9, 18.5);
  AliMaterial(16, "VACUUM$ ", 1e-16, 1e-16, 1e-16, 1e16, 1e16);
  AliMaterial(15, "AIR$      ", 14.61, 7.3, .001205, 30423.24, 67500);
  AliMixture(19, "STAINLESS STEEL$", asteel, zsteel, 7.88, 4, wsteel);
  
  // **************** 
  //     Defines tracking media parameters. 
  //     Les valeurs sont commentees pour laisser le defaut 
  //     a GEANT (version 3-21, page CONS200), f.m. 
  epsil  = .001;  // Tracking precision, 
  stemax = -1.;   // Maximum displacement for multiple scat 
  tmaxfd = -20.;  // Maximum angle due to field deflection 
  deemax = -.3;   // Maximum fractional energy loss, DLS 
  stmin  = -.8;
  
  //    Air 
  
  AliMedium(15, "AIR_L3_US", 15, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
  
  //    Beryllium 
  
  AliMedium(5, "BE_L3_US", 5, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);

  
    //    Aluminium 
  
  AliMedium(4, "AL_L3_US", 4, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);

  //   Vacuum

  AliMedium(16, "VA_L3_US", 16, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
  
  //    Steel 
  
  AliMedium(19, "ST_L3_US", 19, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
}
//
void AliPIPEv3::Undulation(char *undul, Float_t pitch, Float_t thick,
                        Float_t zundul, Float_t rundul, char (*cone)[5])
{
  //
  // RUNDUL   : Internal radius of the undulated chamber
  // THICK    : material thickness
  // PITCH    : one-QUARTER wave of undulation (cm)
  // ZUNDUL   : half length (cm)
  //
  // The undulated structure is desgned as a superposition of eight CONES
  // of suitable sizes, where the inner/outer radius of the cone increases,
  // then decreases, each half of the wave is assumed to be a semicircle,
  // which allows to calculate the thickness and the radii of the cone, by
  // dividing the semicircle into 4 parts of equal arc length.
  // Thus apear the constants 0.293 and 0.707.
  //

  const Float_t kConst1 = .293;
  const Float_t kConst2 = .707;

  // Local variables
  Int_t j, nwave;
  Float_t dcone1[5], dcone2[5], dcone3[5], dcone4[5], dcone5[5],
    dcone6[5], dcone7[5], dcone8[5];
  Float_t xc, yc, zc, dundul[3];
  Int_t *idtmed = fIdtmed->GetArray()-1999;

  // Function Body

  dcone1[0] = kConst1 * pitch / 2;
  dcone1[1] = rundul;
  dcone1[2] = dcone1[1] + thick;
  dcone1[3] = dcone1[1] + kConst2 * pitch;
  dcone1[4] = dcone1[3] + thick;

  dcone2[0] = kConst2 * pitch / 2;
  dcone2[1] = dcone1[3];
  dcone2[2] = dcone1[4];
  dcone2[3] = dcone2[1] + kConst1 * pitch;
  dcone2[4] = dcone2[3] + thick;

  dcone3[0] = dcone2[0];
  dcone3[1] = dcone2[3];
  dcone3[2] = dcone2[4];
  dcone3[3] = dcone2[1];
  dcone3[4] = dcone2[2];

  dcone4[0] = dcone1[0];
  dcone4[1] = dcone1[3];
  dcone4[2] = dcone1[4];
  dcone4[3] = dcone1[1];
  dcone4[4] = dcone1[2];

  dcone5[0] = dcone1[0];
  dcone5[1] = dcone1[1] - thick;
  dcone5[2] = dcone1[1];
  dcone5[3] = dcone5[1] - kConst2 * pitch;
  dcone5[4] = dcone5[3] + thick;

  dcone6[0] = dcone2[0];
  dcone6[1] = dcone5[3];
  dcone6[2] = dcone5[4];
  dcone6[3] = dcone6[1] - kConst1 * pitch;
  dcone6[4] = dcone6[3] + thick;
  dcone7[0] = dcone6[0];
  dcone7[1] = dcone6[3];
  dcone7[2] = dcone6[4];
  dcone7[3] = dcone5[3];
  dcone7[4] = dcone5[4];

  dcone8[0] = dcone5[0];
  dcone8[1] = dcone7[3];
  dcone8[2] = dcone7[4];
  dcone8[3] = dcone5[1];
  dcone8[4] = dcone5[2];

  gMC->Gsvolu(cone[0], "CONE", idtmed[2018], dcone1, 5);
  gMC->Gsvolu(cone[1], "CONE", idtmed[2018], dcone2, 5);
  gMC->Gsvolu(cone[2], "CONE", idtmed[2018], dcone3, 5);
  gMC->Gsvolu(cone[3], "CONE", idtmed[2018], dcone4, 5);
  gMC->Gsvolu(cone[4], "CONE", idtmed[2018], dcone5, 5);
  gMC->Gsvolu(cone[5], "CONE", idtmed[2018], dcone6, 5);
  gMC->Gsvolu(cone[6], "CONE", idtmed[2018], dcone7, 5);
  gMC->Gsvolu(cone[7], "CONE", idtmed[2018], dcone8, 5);
  gMC->Gsatt(cone[0], "SEEN", 0);
  gMC->Gsatt(cone[1], "SEEN", 0);
  gMC->Gsatt(cone[2], "SEEN", 0);
  gMC->Gsatt(cone[3], "SEEN", 0);
  gMC->Gsatt(cone[4], "SEEN", 0);
  gMC->Gsatt(cone[5], "SEEN", 0);
  gMC->Gsatt(cone[6], "SEEN", 0);
  gMC->Gsatt(cone[7], "SEEN", 0);

  // DEFINE AN IMAGINARY TUBE VOLUME FOR UNDULATED CHAMBER, FILL WITH VACUUM

  nwave = Int_t (zundul / (pitch * 2) + .1);
  dundul[2] = pitch * 2 * nwave;
  dundul[1] = rundul + pitch + thick * 2;
  //
  dundul[0] = 1e-4;
  gMC->Gsvolu(undul, "TUBE", idtmed[2015], dundul, 3);

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