made const array arguments of GetPredictedChi2 methods
[u/mrichter/AliRoot.git] / STRUCT / AliPIPEv4.cxx
CommitLineData
820b4d9e 1
2/**************************************************************************
3 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * *
5 * Author: The ALICE Off-line Project. *
6 * Contributors are mentioned in the code where appropriate. *
7 * *
8 * Permission to use, copy, modify and distribute this software and its *
9 * documentation strictly for non-commercial purposes is hereby granted *
10 * without fee, provided that the above copyright notice appears in all *
11 * copies and that both the copyright notice and this permission notice *
12 * appear in the supporting documentation. The authors make no claims *
13 * about the suitability of this software for any purpose. It is *
14 * provided "as is" without express or implied warranty. *
15 **************************************************************************/
16
820b4d9e 17
18//-------------------------------------------------------------------------
af5f3976 19// Beam pipe class for ALICE MFT upgrade
820b4d9e 20// This version uses TGeo
af5f3976 21// Authors:
22// F. Manso
23// A. Morsch
04a81bb0 24// R. Tieulent
820b4d9e 25//-------------------------------------------------------------------------
26
27
28#include <Riostream.h>
29
30#include <TSystem.h>
31#include <TVirtualMC.h>
32#include <TGeoManager.h>
33#include <TGeoMatrix.h>
34#include <TGeoVolume.h>
35#include <TGeoTorus.h>
36#include <TGeoTube.h>
37#include <TGeoCone.h>
38#include <TGeoPcon.h>
39#include <TGeoBBox.h>
40#include <TGeoXtru.h>
41#include <TGeoCompositeShape.h>
42#include <TGeoGlobalMagField.h>
43
44#include "AliConst.h"
45#include "AliMagF.h"
46#include "AliPIPEv4.h"
47#include "AliRun.h"
48#include "AliLog.h"
af5f3976 49
820b4d9e 50ClassImp(AliPIPEv4)
af5f3976 51
820b4d9e 52//_____________________________________________________________________________
236bcb7a 53AliPIPEv4::AliPIPEv4(): AliPIPE()
54
820b4d9e 55{
56 // Constructor
57}
58
59//_____________________________________________________________________________
60AliPIPEv4::AliPIPEv4(const char *name, const char *title)
236bcb7a 61 : AliPIPE(name,title)
820b4d9e 62{
63 // Constructor
64}
65
66
820b4d9e 67
68//___________________________________________
69void AliPIPEv4::CreateGeometry()
70{
71 AliDebug(1,"Create PIPEv4 geometry");
72 //
73 // Class describing the beam pipe geometry
74 //
af5f3976 75 Float_t z, zsh, z0;
820b4d9e 76 //
77 // Rotation Matrices
78 //
79 const Float_t kDegRad = TMath::Pi() / 180.;
80 // Rotation by 180 deg
81 TGeoRotation* rot180 = new TGeoRotation("rot180", 90., 180., 90., 90., 180., 0.);
82 TGeoRotation* rotyz = new TGeoRotation("rotyz", 90., 180., 0., 180., 90., 90.);
83 TGeoRotation* rotxz = new TGeoRotation("rotxz", 0., 0., 90., 90., 90., 180.);
af5f3976 84 //TGeoRotation* rot045 = new TGeoRotation("rot045", 90., 45., 90., 135., 0., 0.);
85 //TGeoRotation* rot135 = new TGeoRotation("rot135", 90. ,135., 90., 225., 0., 0.);
86 //TGeoRotation* rot225 = new TGeoRotation("rot225", 90. ,225., 90., 315., 0., 0.);
87 //TGeoRotation* rot315 = new TGeoRotation("rot315", 90. ,315., 90., 45., 0., 0.);
820b4d9e 88 //
89 // Media
af5f3976 90 //const TGeoMedium* kMedSi = gGeoManager->GetMedium("PIPE_SILICON"); //FM
820b4d9e 91 const TGeoMedium* kMedAir = gGeoManager->GetMedium("PIPE_AIR");
92 const TGeoMedium* kMedAirHigh = gGeoManager->GetMedium("PIPE_AIR_HIGH");
93 const TGeoMedium* kMedVac = gGeoManager->GetMedium("PIPE_VACUUM");
94 const TGeoMedium* kMedInsu = gGeoManager->GetMedium("PIPE_INS_C0");
95 const TGeoMedium* kMedSteel = gGeoManager->GetMedium("PIPE_INOX");
96 const TGeoMedium* kMedBe = gGeoManager->GetMedium("PIPE_BE");
97 const TGeoMedium* kMedCu = gGeoManager->GetMedium("PIPE_CU");
98 //const TGeoMedium* kMedKapton = gGeoManager->GetMedium("PIPE_KAPTON");
af5f3976 99 //const TGeoMedium* kMedAco = gGeoManager->GetMedium("PIPE_ANTICORODAL");
820b4d9e 100 //const TGeoMedium* kMedNEG = gGeoManager->GetMedium("PIPE_NEG COATING");
24679816 101 //const TGeoMedium* kMedAlu = gGeoManager->GetMedium("PIPE_ALU"); // fm
102 const TGeoMedium* kMedAlu2219 = gGeoManager->GetMedium("PIPE_AA2219"); // fm
236bcb7a 103 const TGeoMedium* kMedRohacell = gGeoManager->GetMedium("PIPE_ROHACELL");
104 const TGeoMedium* kMedPolyimide = gGeoManager->GetMedium("PIPE_POLYIMIDE");
105 const TGeoMedium* kMedCarbonFiber = gGeoManager->GetMedium("PIPE_M55J6K");
106
820b4d9e 107 // Top volume
108 TGeoVolume* top = gGeoManager->GetVolume("ALIC");
109 //
110 //
111 ////////////////////////////////////////////////////////////////////////////////
112 // //
113 // The Central Vacuum system //
114 // //
115 ////////////////////////////////////////////////////////////////////////////////
116 //
117 //
118 // The ALICE central beam-pipe according to drawing LHCVC2C_0001
119 // Drawings of sub-elements:
120 //
121 // Pos 7 - Minimised Flange: LHCVFX_P0025
122 // Pos 6 - Standard Flange: STDVFUHV0009
123 // Pos 8 - Bellow: LHCVBX__0001
124 //
125 // Absolute z-coordinates -82.0 - 400.0 cm
126 // Total length: 482.0 cm
127 // It consists of 3 main parts:
128 // CP/2 The flange on the non-absorber side: 36.5 cm
129 // CP/1 The central Be pipe: 405.0 cm
130 // CP/3 The double-bellow and flange on the absorber side: 40.5 cm
820b4d9e 131
132 //
af5f3976 133 /*
820b4d9e 134 // Starting position in z
135 const Float_t kCPz0 = -400.0;
136 // Length of the CP/1 section
137 const Float_t kCP1Length = 405.0;
138 // Length of the CP/2 section
139 const Float_t kCP2Length = 36.5;
140 // Length of the CP/3 section
141 const Float_t kCP3Length = 40.5;
142 // Position of the CP/2 section
143 // const Float_t kCP2pos = kCPz0 + kCP2Length / 2.;
144 // Position of the CP/3 section
145 const Float_t kCP3pos = kCPz0 + kCP2Length + kCP1Length + kCP3Length/2.;
af5f3976 146 */
820b4d9e 147
148
236bcb7a 149 //////////////////// NEW BEAM PIPE GEOMETRY FOR MuonForwardTracker ,
150 // Authors: F. Manso, R. Tieulent
151 // Drawings from C. Gargiulo :
152 // \\cern.ch\dfs\Workspaces\c\cgargiul\EXPERIMENT\ALICE\ALICE_MECHANICS\ALICE_DATA_PACKAGE\IN\DETECTORS\ITS_UPGRADE\1-DESIGN\3D_cad_model\R14_20140311_ALI\
153 //
154 //------------------- Pipe version 4.7 March 2014 -----------------------------
04a81bb0 155
236bcb7a 156 TGeoVolumeAssembly * beamPipeCsideSection = new TGeoVolumeAssembly("BeamPipeCsideSection");
157
158 Float_t fBeryliumSectionOuterRadius = 1.9;
159 Float_t fBeryliumSectionZmax = 44.4;
160 Float_t fBeryliumSectionZmin = -44.4;
161 Float_t fBeryliumSectionThickness = 0.08;
162
163 Float_t fBellowSectionOuterRadius = 2.15;
164 Float_t fCSideBPSOuterRadius = 2.22;
165 Float_t fCSideBPSWallThickness = 0.15;
166 Float_t fBellowSectionZmax = -55.35;
167 Float_t fBellowOuterRadius = 2.8;
168 Float_t fFirstConeAngle = 15. * TMath::DegToRad();
169 Float_t fChangeThicknessAngle = 45. * TMath::DegToRad();
170 // Float_t fCSideBPSLength = 3.53;
171 Float_t fCSideBPSLength = 3.53+1.52;
172 Float_t fDzFirstCone = (fCSideBPSOuterRadius - fBeryliumSectionOuterRadius) / TMath::Tan(fFirstConeAngle);
173 // Float_t fReduceThicknessPartAfterBPSLength = 1.52;
174 Float_t fReduceThicknessPartAfterBPSLength = 0.;
175 Float_t fThinPartBeforeBellowLength = 1.025;
04a81bb0 176
236bcb7a 177
178 Float_t fDistanceBetweenBellows = 2.5;
179
180
181 Float_t fAdaptConeZmax = -77.43;
182 Float_t fAdaptConeZmin = -80.6;
183 Float_t fAdaptConeRmax = 3.0;
184 Float_t fFlangeRmax = 4.3;
185 Float_t fFlangeLength = 1.4;
186
187 Float_t fBellowPlieRadius = 0.17; // radius of bellow plies
188 Float_t fBellowPlieThickness = 0.03; // Thickness of bellow plies 300 microns
189 Int_t fNBellowConvolutions = 7;
190
191
820b4d9e 192
236bcb7a 193 Float_t fZ1 = fBeryliumSectionZmin; // z of Be - Al jonction on the C-side
194 Float_t fZ2 = fBellowSectionZmax +fDzFirstCone ; // z of end of small diameter part (beginning of first cone before the bellow
195 Float_t fZ3 = fBellowSectionZmax +(fCSideBPSOuterRadius - fBellowSectionOuterRadius) / TMath::Tan(fFirstConeAngle); // z of End of first cone part with 0.8mm thickness
196 Float_t fZ4 = fBellowSectionZmax; // z of End of first Cone
197 Float_t fZ5 = fBellowSectionZmax - fCSideBPSLength; // z of End of Beam Pipe support section
198 Float_t fZ6 = fBellowSectionZmax - fCSideBPSLength - (fCSideBPSOuterRadius-fBellowSectionOuterRadius) / TMath::Tan(fChangeThicknessAngle); // z of End of Beam Pipe support section after reduction of thickness
199 Float_t fZ7 = fZ6 - fReduceThicknessPartAfterBPSLength ; // Z of end of 800 microns section after Beam Pipe Support
200 Float_t fZ8 = fZ7 - (fBeryliumSectionThickness-fBellowPlieThickness) / TMath::Tan(fChangeThicknessAngle);
201 Float_t fZ9 = fZ7 - fThinPartBeforeBellowLength; // Z of the start of first bellow
202 Float_t fFirstBellowZmax = fZ9;
203
24679816 204 //---------------- Be pipe around the IP ----------
236bcb7a 205 TGeoPcon* berylliumTube = new TGeoPcon(0., 360., 2);
206 berylliumTube->DefineSection(0,fBeryliumSectionZmax,fBeryliumSectionOuterRadius-fBeryliumSectionThickness,fBeryliumSectionOuterRadius);
207 berylliumTube->DefineSection(1,fBeryliumSectionZmin,fBeryliumSectionOuterRadius-fBeryliumSectionThickness,fBeryliumSectionOuterRadius);
208 TGeoVolume* voberylliumTube = new TGeoVolume("berylliumTube",berylliumTube,kMedBe);
209 voberylliumTube->SetLineColor(kRed);
210 beamPipeCsideSection->AddNode(voberylliumTube,1,new TGeoTranslation(0., 0., 0.));
211
212 TGeoPcon* berylliumTubeVacuum = new TGeoPcon(0., 360., 2);
213 berylliumTubeVacuum->DefineSection(0,fBeryliumSectionZmax, 0.,fBeryliumSectionOuterRadius-fBeryliumSectionThickness);
214 berylliumTubeVacuum->DefineSection(1,fBeryliumSectionZmin, 0.,fBeryliumSectionOuterRadius-fBeryliumSectionThickness);
215 TGeoVolume* voberylliumTubeVacuum = new TGeoVolume("berylliumTubeVacuum",berylliumTubeVacuum,kMedVac);
216 voberylliumTubeVacuum->SetVisibility(0);voberylliumTubeVacuum->SetLineColor(kGreen);
217 beamPipeCsideSection->AddNode(voberylliumTubeVacuum,1,new TGeoTranslation(0., 0., 0.));
24679816 218 //-------------------------------------------------
af5f3976 219
af5f3976 220
236bcb7a 221 //---------------- Al tube ------------------
222 TGeoPcon* aluBeforeBellows = new TGeoPcon(0., 360., 9);
223 aluBeforeBellows->DefineSection(0,fZ1, fBeryliumSectionOuterRadius-fBeryliumSectionThickness,fBeryliumSectionOuterRadius);
224 aluBeforeBellows->DefineSection(1,fZ2,fBeryliumSectionOuterRadius-fBeryliumSectionThickness,fBeryliumSectionOuterRadius);
225 aluBeforeBellows->DefineSection(2,fZ3,fBellowSectionOuterRadius-fBeryliumSectionThickness,fBellowSectionOuterRadius);
226 aluBeforeBellows->DefineSection(3,fZ4,fCSideBPSOuterRadius-fCSideBPSWallThickness,fCSideBPSOuterRadius);
227 aluBeforeBellows->DefineSection(4,fZ5,fCSideBPSOuterRadius-fCSideBPSWallThickness,fCSideBPSOuterRadius);
228 aluBeforeBellows->DefineSection(5,fZ6,fBellowSectionOuterRadius-fBeryliumSectionThickness,fBellowSectionOuterRadius);
229 aluBeforeBellows->DefineSection(6,fZ7,fBellowSectionOuterRadius-fBeryliumSectionThickness,fBellowSectionOuterRadius);
230 aluBeforeBellows->DefineSection(7,fZ8,fBellowSectionOuterRadius-fBeryliumSectionThickness,fBellowSectionOuterRadius-fBeryliumSectionThickness+fBellowPlieThickness);
231 aluBeforeBellows->DefineSection(8,fZ9,fBellowSectionOuterRadius-fBeryliumSectionThickness,fBellowSectionOuterRadius-fBeryliumSectionThickness+fBellowPlieThickness);
232 TGeoVolume* voaluBeforeBellows = new TGeoVolume("aluBeforeBellows",aluBeforeBellows,kMedAlu2219);
233 voaluBeforeBellows->SetLineColor(kBlue);
234 beamPipeCsideSection->AddNode(voaluBeforeBellows,1,new TGeoTranslation(0., 0., 0.));
235
236 TGeoPcon* aluBeforeBellowsVacuum = new TGeoPcon(0., 360., 7);
237 aluBeforeBellowsVacuum->DefineSection(0,fZ1,0.,fBeryliumSectionOuterRadius-fBeryliumSectionThickness);
238 aluBeforeBellowsVacuum->DefineSection(1,fZ2,0.,fBeryliumSectionOuterRadius-fBeryliumSectionThickness);
239 aluBeforeBellowsVacuum->DefineSection(2,fZ3,0.,fBellowSectionOuterRadius-fBeryliumSectionThickness);
240 aluBeforeBellowsVacuum->DefineSection(3,fZ4,0.,fCSideBPSOuterRadius-fCSideBPSWallThickness);
241 aluBeforeBellowsVacuum->DefineSection(4,fZ5,0.,fCSideBPSOuterRadius-fCSideBPSWallThickness);
242 aluBeforeBellowsVacuum->DefineSection(5,fZ6,0.,fBellowSectionOuterRadius-fBeryliumSectionThickness);
243 aluBeforeBellowsVacuum->DefineSection(6,fZ9,0.,fBellowSectionOuterRadius-fBeryliumSectionThickness);
244 TGeoVolume* voaluBeforeBellowsVacuum = new TGeoVolume("aluBeforeBellowsVacuum",aluBeforeBellowsVacuum,kMedVac);
245 voaluBeforeBellowsVacuum->SetVisibility(0);voaluBeforeBellowsVacuum->SetLineColor(kGreen);
246 beamPipeCsideSection->AddNode(voaluBeforeBellowsVacuum,1,new TGeoTranslation(0., 0., 0.));
af5f3976 247 //-------------------------------------------------
248
236bcb7a 249
250 Float_t fBellowLength = fNBellowConvolutions * (4.*fBellowPlieRadius - 2. *fBellowPlieThickness);
251 // ------------------ First Bellow --------------------
252 TGeoVolume* vobellows1 = MakeBellowCside("bellows1", fNBellowConvolutions, fBellowSectionOuterRadius-fBeryliumSectionThickness, fBellowOuterRadius, fBellowPlieRadius ,fBellowPlieThickness);
253 beamPipeCsideSection->AddNode(vobellows1, 1, new TGeoTranslation(0., 0., fFirstBellowZmax-fBellowLength/2. - 2.*fBellowPlieRadius));
254 //------------------------------------------------------
fb4b4c5e 255
236bcb7a 256 Float_t fZ10 = fFirstBellowZmax - fBellowLength; // End of First bellow
257 Float_t fZ12 = fZ10 - fThinPartBeforeBellowLength;
258 Float_t fZ11 = fZ12 + (fBeryliumSectionThickness-fBellowPlieThickness) / TMath::Tan(fChangeThicknessAngle); // End of 300 microns thickness part after first bellow
259 Float_t fZ13 = fZ12 - fDistanceBetweenBellows;
260 Float_t fZ14 = fZ13 -(fBeryliumSectionThickness-fBellowPlieThickness) / TMath::Tan(fChangeThicknessAngle);
261 Float_t fZ15 = fZ14 -fThinPartBeforeBellowLength;
262 Float_t fSecondBellowZmax = fZ15;
af5f3976 263
236bcb7a 264
24679816 265 //---------- Al tube between the bellows ----------
236bcb7a 266 TGeoPcon* tube4 = new TGeoPcon(0., 360., 6);
267 tube4->DefineSection(0,fZ10, fBellowSectionOuterRadius-fBeryliumSectionThickness,fBellowSectionOuterRadius-fBeryliumSectionThickness+fBellowPlieThickness);
268 tube4->DefineSection(1,fZ11,fBellowSectionOuterRadius-fBeryliumSectionThickness,fBellowSectionOuterRadius-fBeryliumSectionThickness+fBellowPlieThickness);
269 tube4->DefineSection(2,fZ12,fBellowSectionOuterRadius-fBeryliumSectionThickness,fBellowSectionOuterRadius);
270 tube4->DefineSection(3,fZ13,fBellowSectionOuterRadius-fBeryliumSectionThickness,fBellowSectionOuterRadius);
271 tube4->DefineSection(4,fZ14,fBellowSectionOuterRadius-fBeryliumSectionThickness,fBellowSectionOuterRadius-fBeryliumSectionThickness+fBellowPlieThickness);
272 tube4->DefineSection(5,fZ15,fBellowSectionOuterRadius-fBeryliumSectionThickness,fBellowSectionOuterRadius-fBeryliumSectionThickness+fBellowPlieThickness);
24679816 273 TGeoVolume* votube4 = new TGeoVolume("votube4",tube4,kMedAlu2219);
274 votube4->SetLineColor(kBlue);
236bcb7a 275 beamPipeCsideSection->AddNode(votube4,1,new TGeoTranslation(0., 0., 0.));
276
277 TGeoPcon* tube4Vacuum = new TGeoPcon(0., 360., 2);
278 tube4Vacuum->DefineSection(0,fZ10,0., fBellowSectionOuterRadius-fBeryliumSectionThickness);
279 tube4Vacuum->DefineSection(1,fZ15,0.,fBellowSectionOuterRadius-fBeryliumSectionThickness);
280 TGeoVolume* votube4Vacuum = new TGeoVolume("tube4Vacuum",tube4Vacuum,kMedVac);
281 votube4Vacuum->SetVisibility(0);
282
283 beamPipeCsideSection->AddNode(votube4Vacuum,1,new TGeoTranslation(0., 0., 0.));
284
285
286 // ------------------ Second Bellow --------------------
287 TGeoVolume* vobellows2 = MakeBellowCside("bellows2", fNBellowConvolutions, fBellowSectionOuterRadius-fBeryliumSectionThickness, fBellowOuterRadius, fBellowPlieRadius ,fBellowPlieThickness);
288 beamPipeCsideSection->AddNode(vobellows2, 1, new TGeoTranslation(0., 0., fSecondBellowZmax-fBellowLength/2. - 2.*fBellowPlieRadius));
289 // -----------------------------------------------------
290
291 Float_t fZ16 = fSecondBellowZmax - fBellowLength; // End of Second bellow
292 Float_t fZ18 = fZ16 - fThinPartBeforeBellowLength;
293 Float_t fZ17 = fZ18 + (fBeryliumSectionThickness-fBellowPlieThickness) / TMath::Tan(fChangeThicknessAngle); // End of 300 microns thickness part after first bellow
294 Float_t fZ19 = fAdaptConeZmax; // Start of the Adpation Cone
295 Float_t fZ20 = fAdaptConeZmin; // End of the Adpation Cone
296 Float_t fZ21 = fAdaptConeZmin - fFlangeLength; // End of the Flange
297
298
299 //----------- 15 deg Conical adaptator + flange ----------
300 TGeoPcon* adaptator = new TGeoPcon(0., 360., 7);
301 adaptator->DefineSection(0,fZ16, fBellowSectionOuterRadius-fBeryliumSectionThickness,fBellowSectionOuterRadius-fBeryliumSectionThickness+fBellowPlieThickness);
302 adaptator->DefineSection(1,fZ17, fBellowSectionOuterRadius-fBeryliumSectionThickness,fBellowSectionOuterRadius-fBeryliumSectionThickness+fBellowPlieThickness);
303 adaptator->DefineSection(2,fZ18, fBellowSectionOuterRadius-fBeryliumSectionThickness,fBellowSectionOuterRadius);
304 adaptator->DefineSection(3,fZ19, fBellowSectionOuterRadius-fBeryliumSectionThickness,fBellowSectionOuterRadius);
305 adaptator->DefineSection(4,fZ20, fAdaptConeRmax-fBeryliumSectionThickness,fAdaptConeRmax);
306 adaptator->DefineSection(5,fZ20, fAdaptConeRmax-fBeryliumSectionThickness,fFlangeRmax);
307 adaptator->DefineSection(6,fZ21, fAdaptConeRmax-fBeryliumSectionThickness,fFlangeRmax);
af5f3976 308 TGeoVolume* voadaptator = new TGeoVolume("voadaptator",adaptator,kMedAlu2219);
236bcb7a 309 voadaptator->SetLineColor(kBlue);
310 beamPipeCsideSection->AddNode(voadaptator,1,new TGeoTranslation(0., 0., 0.));
fb4b4c5e 311
24679816 312 TGeoPcon* adaptatorvide = new TGeoPcon(0., 360., 4);
236bcb7a 313 adaptatorvide->DefineSection(0,fZ16, 0., fBellowSectionOuterRadius-fBeryliumSectionThickness);
314 adaptatorvide->DefineSection(1,fZ19, 0., fBellowSectionOuterRadius-fBeryliumSectionThickness);
315 adaptatorvide->DefineSection(2,fZ20, 0., fAdaptConeRmax-fBeryliumSectionThickness);
316 adaptatorvide->DefineSection(3,fZ21, 0., fAdaptConeRmax-fBeryliumSectionThickness);
24679816 317 TGeoVolume* voadaptatorvide = new TGeoVolume("voadaptatorvide",adaptatorvide,kMedVac);
236bcb7a 318 voadaptatorvide->SetVisibility(0);
319// voadaptatorvide->SetLineColor(kGreen);
320 beamPipeCsideSection->AddNode(voadaptatorvide,1,new TGeoTranslation(0., 0., 0.));
24679816 321 //------------------------------------------------------
af5f3976 322
236bcb7a 323 top->AddNode(beamPipeCsideSection,1);
04a81bb0
AM
324
325 ///////////////////////////////////
326 // Beam Pipe support //
04a81bb0
AM
327 ///////////////////////////////////
328
236bcb7a 329
330 // Beam Pipe Support
331 TGeoVolume *beamPipeSupport = new TGeoVolumeAssembly("BeamPipeSupport");
332 Float_t beamPipesupportZpos = fZ5;
333
334 // Dimensions :
335
336 Float_t supportXdim= 20.67;
337 Float_t beamPipeRingZdim = 4.0;
338 Float_t vespelRmax = 2.3;
339 Float_t vespelRmin = 2.22;
340 Float_t beampipeCarbonCollarRmin = 2.4;
341 Float_t beampipeCarbonCollarRmax = 2.7;
342
343 Float_t fixationCarbonCollarRmin = 1.5;
344 Float_t fixationCarbonCollarRmax = 1.7;
345 Float_t fixationCarbonCollarDZ = 2.5;
346
347
348 Float_t skinThickness = 0.1;
349 Float_t skinXdim = 14.25;
350 Float_t skinYdim = 1.;
351 Float_t skinZdim = fixationCarbonCollarDZ;
352 Float_t carbonEarsXdim = 1.01;
353 Float_t carbonEarsYdim = 0.2;
354 Float_t carbonEarsZdim = fixationCarbonCollarDZ;
355
356 // Support Bar
357 TGeoVolumeAssembly *supportBar = new TGeoVolumeAssembly("BPS_SupportBar");
358
359 TGeoBBox * carbonSkinBPS = new TGeoBBox(skinXdim/2.,skinYdim/2.,skinZdim/2.);
360 carbonSkinBPS->SetName("carbonSkinBPS");
361
362 TGeoBBox * foambarBPS = new TGeoBBox("foambarBPS",skinXdim/2.-skinThickness,skinYdim/2.-skinThickness,skinZdim/2.-skinThickness/2.);
363 TGeoBBox * carbonEarsBPS = new TGeoBBox(carbonEarsXdim/2.,carbonEarsYdim/2.,carbonEarsZdim/2.);
364 carbonEarsBPS->SetName("carbonEarsBPS");
365
366 TGeoTranslation * transBP1 = new TGeoTranslation("transBP1",(skinXdim+carbonEarsXdim)/2.,0.,0.);
367 transBP1->RegisterYourself();
368 TGeoTranslation * transBP2 = new TGeoTranslation("transBP2",-(skinXdim+carbonEarsXdim)/2.,0.,0.);
369 transBP2->RegisterYourself();
370 TGeoCompositeShape *supportBarCarbon = new TGeoCompositeShape("BPS_supportBarCarbon", "(carbonSkinBPS-foambarBPS)+carbonEarsBPS:transBP1+carbonEarsBPS:transBP2");
371
372 TGeoVolume *supportBarCarbonVol = new TGeoVolume("BPS_supportBarCarbon",supportBarCarbon,kMedCarbonFiber);
373 supportBarCarbonVol->SetLineColor(kGray+3);
374
375 supportBar->AddNode(supportBarCarbonVol, 1, new TGeoTranslation(skinXdim/2.+carbonEarsXdim+beampipeCarbonCollarRmax,0,0));
376 supportBar->AddNode(supportBarCarbonVol, 2, new TGeoTranslation(-(skinXdim/2.+carbonEarsXdim+beampipeCarbonCollarRmax),0,0));
377
378 TGeoVolume *foamVol = new TGeoVolume("supportBarFoam",foambarBPS,kMedRohacell);
379 foamVol->SetLineColor(kGray);
380 supportBar->AddNode(foamVol, 1, new TGeoTranslation(skinXdim/2.+carbonEarsXdim+beampipeCarbonCollarRmax,0,0));
381 supportBar->AddNode(foamVol, 2, new TGeoTranslation(-(skinXdim/2.+carbonEarsXdim+beampipeCarbonCollarRmax),0,0));
382
383 beamPipeSupport->AddNode(supportBar,1);
384
385
386 // Fixation to wings
387
388 TGeoVolumeAssembly *fixationToWings = new TGeoVolumeAssembly("BPS_fixationToWings");
389
390 Float_t delatX = 0.1;
391
392 TGeoTubeSeg * fixationTube = new TGeoTubeSeg(fixationCarbonCollarRmin,fixationCarbonCollarRmax,fixationCarbonCollarDZ/2.,-90.,90.);
393 fixationTube->SetName("fixationTube");
394 TGeoBBox * fixationToBar = new TGeoBBox(carbonEarsXdim/2.+delatX,carbonEarsYdim/2.,carbonEarsZdim/2.);
395 fixationToBar->SetName("fixationToBar");
396
397 TGeoTranslation * transBP3 = new TGeoTranslation("transBP3",fixationCarbonCollarRmax+carbonEarsXdim/2.-delatX,carbonEarsYdim,0.);
398 transBP3->RegisterYourself();
399 TGeoTranslation * transBP4 = new TGeoTranslation("transBP4",fixationCarbonCollarRmax+carbonEarsXdim/2.-delatX,-carbonEarsYdim,0.);
400 transBP4->RegisterYourself();
401 TGeoCompositeShape *fixationToWing = new TGeoCompositeShape("fixationToWing", "fixationTube+fixationToBar:transBP3+fixationToBar:transBP4");
402
403 TGeoVolume *fixationToWingVol = new TGeoVolume("fixationToWing",fixationToWing,kMedCarbonFiber);
404 fixationToWingVol->SetLineColor(kGray+2);
405
406
407 fixationToWings->AddNode(fixationToWingVol,1, new TGeoTranslation(-supportXdim,0,0));
408 fixationToWings->AddNode(fixationToWingVol,2, new TGeoCombiTrans(+supportXdim,0,0,new TGeoRotation("rot",0.,0.,180.)));
409
410
411 beamPipeSupport->AddNode(fixationToWings,1);
412
413
414 // Fixation to pipe
415
416 TGeoVolumeAssembly *fixationToPipe = new TGeoVolumeAssembly("fixationToPipe");
417
418 TGeoTubeSeg * pipeSupportTubeCarbon = new TGeoTubeSeg(beampipeCarbonCollarRmin,beampipeCarbonCollarRmax,fixationCarbonCollarDZ/2.,0.,180.);
419 pipeSupportTubeCarbon->SetName("pipeSupportTubeCarbon");
420
421 TGeoBBox * fixationTubeToBar = new TGeoBBox(carbonEarsXdim/2.+delatX,carbonEarsYdim/2.,carbonEarsZdim/2.);
422 fixationTubeToBar->SetName("fixationTubeToBar");
423 TGeoBBox * hole = new TGeoBBox((beampipeCarbonCollarRmax-vespelRmin)/2.,carbonEarsYdim/2.,carbonEarsZdim/2.+1e-3);
424 hole->SetName("hole");
425
426 TGeoTranslation * transBP5 = new TGeoTranslation("transBP5",beampipeCarbonCollarRmax+carbonEarsXdim/2.-delatX,carbonEarsYdim,0.);
427 transBP5->RegisterYourself();
428 TGeoTranslation * transBP6 = new TGeoTranslation("transBP6",-(beampipeCarbonCollarRmax+carbonEarsXdim/2.-delatX),carbonEarsYdim,0.);
429 transBP6->RegisterYourself();
430 TGeoTranslation * transBP7 = new TGeoTranslation("transBP7",(beampipeCarbonCollarRmax+vespelRmin)/2.,0.,0.);
431 transBP7->RegisterYourself();
432 TGeoTranslation * transBP8 = new TGeoTranslation("transBP8",-((beampipeCarbonCollarRmax+vespelRmin)/2.),0.,0.);
433 transBP8->RegisterYourself();
434 TGeoCompositeShape *halfFixationToPipe = new TGeoCompositeShape("halfFixationToPipe", "(pipeSupportTubeCarbon-hole:transBP7-hole:transBP8)+fixationTubeToBar:transBP5+fixationTubeToBar:transBP6");
435
436 TGeoVolume *halfFixationToPipeVol = new TGeoVolume("halfFixationToPipe",halfFixationToPipe,kMedCarbonFiber);
437 halfFixationToPipeVol->SetLineColor(kRed+2);
438
439 fixationToPipe->AddNode(halfFixationToPipeVol,1);
440 fixationToPipe->AddNode(halfFixationToPipeVol,2, new TGeoCombiTrans(0,0,0,new TGeoRotation("rot",0.,0.,180.)));
441
442 beamPipeSupport->AddNode(fixationToPipe,1);
443
444
445 // Beam Pipe Ring
446
447 TGeoVolumeAssembly *beamPipeRing = new TGeoVolumeAssembly("beamPipeRing");
448
449 TGeoTube * beamPipeRingCarbon = new TGeoTube(vespelRmax,beampipeCarbonCollarRmin,beamPipeRingZdim/2.);
450 TGeoVolume *beamPipeRingCarbonVol = new TGeoVolume("beamPipeRingCarbon",beamPipeRingCarbon,kMedCarbonFiber);
451 beamPipeRingCarbonVol->SetLineColor(kGreen+2);
452 beamPipeRing->AddNode(beamPipeRingCarbonVol,1, new TGeoTranslation(0.,0,(beamPipeRingZdim-fixationCarbonCollarDZ)/2.));
453
454 TGeoTube * beamPipeRingVespel = new TGeoTube(vespelRmin,vespelRmax,beamPipeRingZdim/2.);
455 TGeoVolume *beamPipeRingVespelVol = new TGeoVolume("beamPipeRingVespel",beamPipeRingVespel,kMedPolyimide);
456 beamPipeRingVespelVol->SetLineColor(kGreen+4);
457 beamPipeRing->AddNode(beamPipeRingVespelVol,1, new TGeoTranslation(0.,0,(beamPipeRingZdim-fixationCarbonCollarDZ)/2.));
458
459 beamPipeSupport->AddNode(beamPipeRing,1);
460 beamPipeSupport->SetVisibility(0);
461
462 top->AddNode(beamPipeSupport,1,new TGeoTranslation(0.,0,beamPipesupportZpos+fixationCarbonCollarDZ/2.));
04a81bb0
AM
463
464
465
466 ///////////// END NEW BEAM PIPE GEOMETRY fOR MFT ////////////////////
820b4d9e 467
236bcb7a 468 ////////////////////////////////////////////////////////////////////////////////
469 // Side A section after Beryllium
470 // !!!!!! THIS PART NEED TO BE WORKED OUT !!!!!
471 ////////////////////////////////////////////////////////////////////////////////
820b4d9e 472
236bcb7a 473 //---------------- Al tube ------------------
474 Float_t fSmallRadiusZmax =fBeryliumSectionZmax + 20.43;
475 Float_t fLargeRadiusZmin =fSmallRadiusZmax + 2.61;
476 Float_t fAdaptConeSideAZmin =fLargeRadiusZmin + 200.; // THIS PART NEED TO BE WORKED OUT
477 Float_t fAluSideARmax = 2.5;
478
479 TGeoPcon* aluSideA = new TGeoPcon(0., 360., 4);
480 aluSideA->DefineSection(0,fBeryliumSectionZmax, fBeryliumSectionOuterRadius-fBeryliumSectionThickness,fBeryliumSectionOuterRadius);
481 aluSideA->DefineSection(1,fSmallRadiusZmax, fBeryliumSectionOuterRadius-fBeryliumSectionThickness,fBeryliumSectionOuterRadius);
482 aluSideA->DefineSection(2,fLargeRadiusZmin, fAluSideARmax-fBeryliumSectionThickness,fAluSideARmax);
483 aluSideA->DefineSection(3,fAdaptConeSideAZmin, fAluSideARmax-fBeryliumSectionThickness,fAluSideARmax);
484 TGeoVolume* voaluSideA = new TGeoVolume("aluSideA",aluSideA,kMedAlu2219);
485 voaluSideA->SetLineColor(kBlue);
486 top->AddNode(voaluSideA,1,new TGeoTranslation(0., 0., 0.));
487
488 TGeoPcon* aluSideAVac = new TGeoPcon(0., 360., 4);
489 aluSideAVac->DefineSection(0,fBeryliumSectionZmax, 0., fBeryliumSectionOuterRadius-fBeryliumSectionThickness);
490 aluSideAVac->DefineSection(1,fSmallRadiusZmax, 0., fBeryliumSectionOuterRadius-fBeryliumSectionThickness);
491 aluSideAVac->DefineSection(2,fLargeRadiusZmin, 0., fAluSideARmax-fBeryliumSectionThickness);
492 aluSideAVac->DefineSection(3,fAdaptConeSideAZmin, 0., fAluSideARmax-fBeryliumSectionThickness);
493 TGeoVolume* voaluSideAVac = new TGeoVolume("aluSideAVac",aluSideAVac,kMedVac);
494 voaluSideAVac->SetLineColor(kGreen);
495 voaluSideAVac->SetVisibility(0);
496 top->AddNode(voaluSideAVac,1,new TGeoTranslation(0., 0., 0.));
497 //-------------------------------------------------
820b4d9e 498
236bcb7a 499
820b4d9e 500 ////////////////////////////////////////////////////////////////////////////////
501 // //
502 // RB24/1 //
503 // //
504 ////////////////////////////////////////////////////////////////////////////////
505 //
506 //
507 // Drawing LHCVC2U_0001
508 // Copper Tube RB24/1 393.5 cm
509 // Warm module VMACA 18.0 cm
510 // Annular Ion Pump 35.0 cm
511 // Valve 7.5 cm
512 // Warm module VMABC 28.0 cm
513 // ================================
514 // 462.0 cm
515 //
516
517
518 // Copper Tube RB24/1
519 const Float_t kRB24CuTubeL = 393.5;
520 const Float_t kRB24CuTubeRi = 8.0/2.;
521 const Float_t kRB24CuTubeRo = 8.4/2.;
522 const Float_t kRB24CuTubeFRo = 7.6;
523 const Float_t kRB24CuTubeFL = 1.86;
524
525 TGeoVolume* voRB24CuTubeM = new TGeoVolume("voRB24CuTubeM",
526 new TGeoTube(0., kRB24CuTubeRo, kRB24CuTubeL/2.), kMedVac);
527 voRB24CuTubeM->SetVisibility(0);
528 TGeoVolume* voRB24CuTube = new TGeoVolume("voRB24CuTube",
529 new TGeoTube(kRB24CuTubeRi, kRB24CuTubeRo, kRB24CuTubeL/2.), kMedCu);
530 voRB24CuTubeM->AddNode(voRB24CuTube, 1, gGeoIdentity);
531 // Air outside tube with higher transport cuts
532 TGeoVolume* voRB24CuTubeA = new TGeoVolume("voRB24CuTubeA",
533 new TGeoTube(25., 100., kRB24CuTubeL/2.), kMedAirHigh);
534 voRB24CuTubeA->SetVisibility(0);
535 // Simplified DN 100 Flange
536 TGeoVolume* voRB24CuTubeF = new TGeoVolume("voRB24CuTubeF",
537 new TGeoTube(kRB24CuTubeRo, kRB24CuTubeFRo, kRB24CuTubeFL/2.), kMedSteel);
538
539 // Warm Module Type VMACA
540 // LHCVMACA_0002
541 //
542 // Pos 1 Warm Bellows DN100 LHCVBU__0012
543 // Pos 2 RF Contact D80 LHCVSR__0005
544 // Pos 3 Trans. Tube Flange LHCVSR__0065
545 // [Pos 4 Hex. Countersunk Screw Bossard BN4719]
546 // [Pos 5 Tension spring LHCVSR__0011]
547 //
548 //
549 //
550 // Pos1 Warm Bellows DN100
551 // Pos1.1 Bellows LHCVBU__0006
552 //
553 //
554 // Connection Tubes
555 // Connection tube inner r
556 const Float_t kRB24B1ConTubeRin = 10.0/2.;
557 // Connection tube outer r
558 const Float_t kRB24B1ConTubeRou = 10.3/2.;
559 // Connection tube length
560 const Float_t kRB24B1ConTubeL = 2.5;
561 //
562 const Float_t kRB24B1CompL = 16.00; // Length of the compensator
563 const Float_t kRB24B1BellowRi = 10.25/2.; // Bellow inner radius
564 const Float_t kRB24B1BellowRo = 11.40/2.; // Bellow outer radius
565 const Int_t kRB24B1NumberOfPlies = 27; // Number of plies
566 const Float_t kRB24B1BellowUndL = 11.00; // Length of undulated region
567 const Float_t kRB24B1PlieThickness = 0.015; // Plie thickness
568
569 const Float_t kRB24B1PlieRadius =
570 (kRB24B1BellowUndL + (2. * kRB24B1NumberOfPlies - 2.) * kRB24B1PlieThickness) / (4. * kRB24B1NumberOfPlies);
af5f3976 571
820b4d9e 572 const Float_t kRB24B1ProtTubeThickness = 0.02; // Thickness of the protection tube
573 const Float_t kRB24B1ProtTubeLength = 4.2; // Length of the protection tube
574
575 const Float_t kRB24B1RFlangeL = 1.86; // Length of the flanges
576 const Float_t kRB24B1RFlangeLO = 0.26; // Flange overlap
577 const Float_t kRB24B1RFlangeRO = 11.18/2; // Inner radius at Flange overlap
578 const Float_t kRB24B1RFlangeRou = 15.20/2.; // Outer radius of flange
579 const Float_t kRB24B1RFlangeRecess = 0.98; // Flange recess
580 const Float_t kRB24B1L = kRB24B1CompL + 2. * (kRB24B1RFlangeL - kRB24B1RFlangeRecess);
581
582 ///
583 //
584 // Bellow mother volume
585 TGeoPcon* shRB24B1BellowM = new TGeoPcon(0., 360., 14);
586 // Connection Tube and Flange
587 z = 0.;
588 shRB24B1BellowM->DefineSection( 0, z, 0., kRB24B1RFlangeRou);
589 z += kRB24B1RFlangeLO;
590 shRB24B1BellowM->DefineSection( 1, z, 0., kRB24B1RFlangeRou);
591 shRB24B1BellowM->DefineSection( 2, z, 0., kRB24B1RFlangeRou);
592 z = kRB24B1RFlangeL;
593 shRB24B1BellowM->DefineSection( 3, z, 0., kRB24B1RFlangeRou);
594 shRB24B1BellowM->DefineSection( 4, z, 0., kRB24B1ConTubeRou);
595 z = kRB24B1ConTubeL + kRB24B1RFlangeL - kRB24B1RFlangeRecess;
596 shRB24B1BellowM->DefineSection( 5, z, 0., kRB24B1ConTubeRou);
597 // Plie
598 shRB24B1BellowM->DefineSection( 6, z, 0., kRB24B1BellowRo + kRB24B1ProtTubeThickness);
599 z += kRB24B1BellowUndL;
600 shRB24B1BellowM->DefineSection( 7, z, 0., kRB24B1BellowRo + kRB24B1ProtTubeThickness);
601 shRB24B1BellowM->DefineSection( 8, z, 0., kRB24B1ConTubeRou);
602 // Connection Tube and Flange
603 z = kRB24B1L - shRB24B1BellowM->GetZ(3);
604 shRB24B1BellowM->DefineSection( 9, z, 0., kRB24B1ConTubeRou);
605 shRB24B1BellowM->DefineSection(10, z, 0., kRB24B1RFlangeRou);
606 z = kRB24B1L - shRB24B1BellowM->GetZ(1);
607 shRB24B1BellowM->DefineSection(11, z, 0., kRB24B1RFlangeRou);
608 shRB24B1BellowM->DefineSection(12, z, 0., kRB24B1RFlangeRou);
609 z = kRB24B1L - shRB24B1BellowM->GetZ(0);
610 shRB24B1BellowM->DefineSection(13, z, 0., kRB24B1RFlangeRou);
611
612 TGeoVolume* voRB24B1BellowM = new TGeoVolume("RB24B1BellowM", shRB24B1BellowM, kMedVac);
613 voRB24B1BellowM->SetVisibility(0);
614 //
615 // Bellow Section
616 TGeoVolume* voRB24B1Bellow
617 = MakeBellow("RB24B1", kRB24B1NumberOfPlies, kRB24B1BellowRi, kRB24B1BellowRo,
618 kRB24B1BellowUndL, kRB24B1PlieRadius ,kRB24B1PlieThickness);
619 voRB24B1Bellow->SetVisibility(0);
620
621 //
622 // End Parts (connection tube)
623 TGeoVolume* voRB24B1CT = new TGeoVolume("RB24B1CT", new TGeoTube(kRB24B1ConTubeRin, kRB24B1ConTubeRou, kRB24B1ConTubeL/2.), kMedSteel);
624 //
625 // Protection Tube
626 TGeoVolume* voRB24B1PT = new TGeoVolume("RB24B1PT", new TGeoTube(kRB24B1BellowRo, kRB24B1BellowRo + kRB24B1ProtTubeThickness,
627 kRB24B1ProtTubeLength / 2.), kMedSteel);
628
629 z = kRB24B1ConTubeL/2. + (kRB24B1RFlangeL - kRB24B1RFlangeRecess);
630
631 voRB24B1BellowM->AddNode(voRB24B1CT, 1, new TGeoTranslation(0., 0., z));
632 z += (kRB24B1ConTubeL/2.+ kRB24B1BellowUndL/2.);
633 voRB24B1BellowM->AddNode(voRB24B1Bellow, 1, new TGeoTranslation(0., 0., z));
634 z += (kRB24B1BellowUndL/2. + kRB24B1ConTubeL/2);
635 voRB24B1BellowM->AddNode(voRB24B1CT, 2, new TGeoTranslation(0., 0., z));
636 z = kRB24B1ConTubeL + kRB24B1ProtTubeLength / 2. + 1. + kRB24B1RFlangeLO;
637 voRB24B1BellowM->AddNode(voRB24B1PT, 1, new TGeoTranslation(0., 0., z));
638 z += kRB24B1ProtTubeLength + 0.6;
639 voRB24B1BellowM->AddNode(voRB24B1PT, 2, new TGeoTranslation(0., 0., z));
640
641
642
643 // Pos 1/2 Rotatable Flange LHCVBU__0013
644 // Pos 1/3 Flange DN100/103 LHCVBU__0018
645 // The two flanges can be represented by the same volume
646 // Outer Radius (including the outer movable ring).
647 // The inner ring has a diameter of 12.04 cm
648
649
650 TGeoPcon* shRB24B1RFlange = new TGeoPcon(0., 360., 10);
651 z = 0.;
652 shRB24B1RFlange->DefineSection(0, z, 10.30/2., kRB24B1RFlangeRou);
653 z += 0.55; // 5.5 mm added for outer ring
654 z += 0.43;
655 shRB24B1RFlange->DefineSection(1, z, 10.30/2., kRB24B1RFlangeRou);
656 shRB24B1RFlange->DefineSection(2, z, 10.06/2., kRB24B1RFlangeRou);
657 z += 0.15;
658 shRB24B1RFlange->DefineSection(3, z, 10.06/2., kRB24B1RFlangeRou);
659 // In reality this part is rounded
660 shRB24B1RFlange->DefineSection(4, z, 10.91/2., kRB24B1RFlangeRou);
661 z += 0.15;
662 shRB24B1RFlange->DefineSection(5, z, 10.91/2., kRB24B1RFlangeRou);
663 shRB24B1RFlange->DefineSection(6, z, 10.06/2., kRB24B1RFlangeRou);
664 z += 0.32;
665 shRB24B1RFlange->DefineSection(7, z, 10.06/2., kRB24B1RFlangeRou);
666 shRB24B1RFlange->DefineSection(8, z, kRB24B1RFlangeRO, kRB24B1RFlangeRou);
667 z += kRB24B1RFlangeLO;
668 shRB24B1RFlange->DefineSection(9, z, kRB24B1RFlangeRO, kRB24B1RFlangeRou);
669
670 TGeoVolume* voRB24B1RFlange = new TGeoVolume("RB24B1RFlange", shRB24B1RFlange, kMedSteel);
671
672
673 z = kRB24B1L - kRB24B1RFlangeL;
674 voRB24B1BellowM->AddNode(voRB24B1RFlange, 1, new TGeoTranslation(0., 0., z));
675 z = kRB24B1RFlangeL;
676 voRB24B1BellowM->AddNode(voRB24B1RFlange, 2, new TGeoCombiTrans(0., 0., z, rot180));
677 //
678 // Pos 2 RF Contact D80 LHCVSR__0005
679 //
680 // Pos 2.1 RF Contact Flange LHCVSR__0003
681 //
682 TGeoPcon* shRB24B1RCTFlange = new TGeoPcon(0., 360., 6);
683 const Float_t kRB24B1RCTFlangeRin = 8.06/2. + 0.05; // Inner radius
684 const Float_t kRB24B1RCTFlangeL = 1.45; // Length
685
686 z = 0.;
687 shRB24B1RCTFlange->DefineSection(0, z, kRB24B1RCTFlangeRin, 8.20/2.);
688 z += 0.15;
689 shRB24B1RCTFlange->DefineSection(1, z, kRB24B1RCTFlangeRin, 8.20/2.);
690 shRB24B1RCTFlange->DefineSection(2, z, kRB24B1RCTFlangeRin, 8.60/2.);
691 z += 1.05;
692 shRB24B1RCTFlange->DefineSection(3, z, kRB24B1RCTFlangeRin, 8.60/2.);
693 shRB24B1RCTFlange->DefineSection(4, z, kRB24B1RCTFlangeRin, 11.16/2.);
694 z += 0.25;
695 shRB24B1RCTFlange->DefineSection(5, z, kRB24B1RCTFlangeRin, 11.16/2.);
696 TGeoVolume* voRB24B1RCTFlange = new TGeoVolume("RB24B1RCTFlange", shRB24B1RCTFlange, kMedCu);
697 z = kRB24B1L - kRB24B1RCTFlangeL;
698
699 voRB24B1BellowM->AddNode(voRB24B1RCTFlange, 1, new TGeoTranslation(0., 0., z));
700 //
701 // Pos 2.2 RF-Contact LHCVSR__0004
702 //
703 TGeoPcon* shRB24B1RCT = new TGeoPcon(0., 360., 3);
704 const Float_t kRB24B1RCTRin = 8.00/2.; // Inner radius
705 const Float_t kRB24B1RCTCRin = 8.99/2.; // Max. inner radius conical section
706 const Float_t kRB24B1RCTL = 11.78; // Length
707 const Float_t kRB24B1RCTSL = 10.48; // Length of straight section
708 const Float_t kRB24B1RCTd = 0.03; // Thickness
709
710 z = 0;
711 shRB24B1RCT->DefineSection(0, z, kRB24B1RCTCRin, kRB24B1RCTCRin + kRB24B1RCTd);
712 z = kRB24B1RCTL - kRB24B1RCTSL;
713 // In the (VSR0004) this section is straight in (LHCVC2U_0001) it is conical ????
714 shRB24B1RCT->DefineSection(1, z, kRB24B1RCTRin + 0.35, kRB24B1RCTRin + 0.35 + kRB24B1RCTd);
715 z = kRB24B1RCTL - 0.03;
716 shRB24B1RCT->DefineSection(2, z, kRB24B1RCTRin, kRB24B1RCTRin + kRB24B1RCTd);
717
718 TGeoVolume* voRB24B1RCT = new TGeoVolume("RB24B1RCT", shRB24B1RCT, kMedCu);
719 z = kRB24B1L - kRB24B1RCTL - 0.45;
720 voRB24B1BellowM->AddNode(voRB24B1RCT, 1, new TGeoTranslation(0., 0., z));
721
722 //
723 // Pos 3 Trans. Tube Flange LHCVSR__0065
724 //
725 // Pos 3.1 Transition Tube D53 LHCVSR__0064
726 // Pos 3.2 Transition Flange LHCVSR__0060
727 // Pos 3.3 Transition Tube LHCVSR__0058
728 TGeoPcon* shRB24B1TTF = new TGeoPcon(0., 360., 7);
729 // Flange
730 z = 0.;
731 shRB24B1TTF->DefineSection(0, z, 6.30/2., 11.16/2.);
732 z += 0.25;
733 shRB24B1TTF->DefineSection(1, z, 6.30/2., 11.16/2.);
734 shRB24B1TTF->DefineSection(2, z, 6.30/2., 9.3/2.);
735 z += 0.55;
736 shRB24B1TTF->DefineSection(3, z, 6.30/2., 9.3/2.);
737 // Tube
738 shRB24B1TTF->DefineSection(4, z, 6.30/2., 6.7/2.);
739 z += 5.80;
740 shRB24B1TTF->DefineSection(5, z, 6.30/2., 6.7/2.);
741 // Transition Tube
742 z += 3.75;
743 shRB24B1TTF->DefineSection(6, z, 8.05/2., 8.45/2.);
744 TGeoVolume* voRB24B1TTF = new TGeoVolume("RB24B1TTF", shRB24B1TTF, kMedSteel);
745 z = 0.;
746 voRB24B1BellowM->AddNode(voRB24B1TTF, 1, new TGeoTranslation(0., 0., z));
747
748 // Annular Ion Pump
749 // LHCVC2U_0003
750 //
751 // Pos 1 Rotable Flange LHCVFX__0031
752 // Pos 2 RF Screen Tube LHCVC2U_0005
753 // Pos 3 Shell LHCVC2U_0007
754 // Pos 4 Extruded Shell LHCVC2U_0006
755 // Pos 5 Feedthrough Tube LHCVC2U_0004
756 // Pos 6 Tubulated Flange STDVFUHV0021
757 // Pos 7 Fixed Flange LHCVFX__0032
758 // Pos 8 Pumping Elements
759
760 //
761 // Pos 1 Rotable Flange LHCVFX__0031
762 // pos 7 Fixed Flange LHCVFX__0032
763 //
764 // Mother volume
765 const Float_t kRB24AIpML = 35.;
766
767 TGeoVolume* voRB24AIpM = new TGeoVolume("voRB24AIpM", new TGeoTube(0., 10., kRB24AIpML/2.), kMedAir);
768 voRB24AIpM->SetVisibility(0);
769
770 //
771 // Length 35 cm
772 // Flange 2 x 1.98 = 3.96
773 // Tube = 32.84
774 //==========================
775 // 36.80
776 // Overlap 2 * 0.90 = 1.80
777
778 const Float_t kRB24IpRFD1 = 0.68; // Length of section 1
779 const Float_t kRB24IpRFD2 = 0.30; // Length of section 2
780 const Float_t kRB24IpRFD3 = 0.10; // Length of section 3
781 const Float_t kRB24IpRFD4 = 0.35; // Length of section 4
782 const Float_t kRB24IpRFD5 = 0.55; // Length of section 5
783
784 const Float_t kRB24IpRFRo = 15.20/2.; // Flange outer radius
785 const Float_t kRB24IpRFRi1 = 6.30/2.; // Flange inner radius section 1
786 const Float_t kRB24IpRFRi2 = 6.00/2.; // Flange inner radius section 2
787 const Float_t kRB24IpRFRi3 = 5.84/2.; // Flange inner radius section 3
788 const Float_t kRB24IpRFRi4 = 6.00/2.; // Flange inner radius section 1
789 const Float_t kRB24IpRFRi5 = 10.50/2.; // Flange inner radius section 2
790
791 TGeoPcon* shRB24IpRF = new TGeoPcon(0., 360., 9);
792 z0 = 0.;
793 shRB24IpRF->DefineSection(0, z0, kRB24IpRFRi1, kRB24IpRFRo);
794 z0 += kRB24IpRFD1;
795 shRB24IpRF->DefineSection(1, z0, kRB24IpRFRi2, kRB24IpRFRo);
796 z0 += kRB24IpRFD2;
797 shRB24IpRF->DefineSection(2, z0, kRB24IpRFRi2, kRB24IpRFRo);
798 shRB24IpRF->DefineSection(3, z0, kRB24IpRFRi3, kRB24IpRFRo);
799 z0 += kRB24IpRFD3;
800 shRB24IpRF->DefineSection(4, z0, kRB24IpRFRi3, kRB24IpRFRo);
801 shRB24IpRF->DefineSection(5, z0, kRB24IpRFRi4, kRB24IpRFRo);
802 z0 += kRB24IpRFD4;
803 shRB24IpRF->DefineSection(6, z0, kRB24IpRFRi4, kRB24IpRFRo);
804 shRB24IpRF->DefineSection(7, z0, kRB24IpRFRi5, kRB24IpRFRo);
805 z0 += kRB24IpRFD5;
806 shRB24IpRF->DefineSection(8, z0, kRB24IpRFRi5, kRB24IpRFRo);
807
808 TGeoVolume* voRB24IpRF = new TGeoVolume("RB24IpRF", shRB24IpRF, kMedSteel);
809
810 //
811 // Pos 2 RF Screen Tube LHCVC2U_0005
812 //
813
814 //
815 // Tube
816 Float_t kRB24IpSTTL = 32.84; // Total length of the tube
817 Float_t kRB24IpSTTRi = 5.80/2.; // Inner Radius
818 Float_t kRB24IpSTTRo = 6.00/2.; // Outer Radius
819 TGeoVolume* voRB24IpSTT = new TGeoVolume("RB24IpSTT", new TGeoTube(kRB24IpSTTRi, kRB24IpSTTRo, kRB24IpSTTL/2.), kMedSteel);
820 // Screen
821 Float_t kRB24IpSTCL = 0.4; // Lenth of the crochet detail
822 // Length of the screen
823 Float_t kRB24IpSTSL = 9.00 - 2. * kRB24IpSTCL;
824 // Rel. position of the screen
825 Float_t kRB24IpSTSZ = 7.00 + kRB24IpSTCL;
826 TGeoVolume* voRB24IpSTS = new TGeoVolume("RB24IpSTS", new TGeoTube(kRB24IpSTTRi, kRB24IpSTTRo, kRB24IpSTSL/2.), kMedSteel);
827 // Vacuum
828 TGeoVolume* voRB24IpSTV = new TGeoVolume("RB24IpSTV", new TGeoTube(0., kRB24IpSTTRi, kRB24AIpML/2.), kMedVac);
829 //
830 voRB24IpSTT->AddNode(voRB24IpSTS, 1, new TGeoTranslation(0., 0., kRB24IpSTSZ - kRB24IpSTTL/2. + kRB24IpSTSL/2.));
831
832 // Crochets
833 // Inner radius
834 Float_t kRB24IpSTCRi = kRB24IpSTTRo + 0.25;
835 // Outer radius
836 Float_t kRB24IpSTCRo = kRB24IpSTTRo + 0.35;
837 // Length of 1stsection
838 Float_t kRB24IpSTCL1 = 0.15;
839 // Length of 2nd section
840 Float_t kRB24IpSTCL2 = 0.15;
841 // Length of 3rd section
842 Float_t kRB24IpSTCL3 = 0.10;
843 // Rel. position of 1st Crochet
844
845
846 TGeoPcon* shRB24IpSTC = new TGeoPcon(0., 360., 5);
847 z0 = 0;
848 shRB24IpSTC->DefineSection(0, z0, kRB24IpSTCRi, kRB24IpSTCRo);
849 z0 += kRB24IpSTCL1;
850 shRB24IpSTC->DefineSection(1, z0, kRB24IpSTCRi, kRB24IpSTCRo);
851 shRB24IpSTC->DefineSection(2, z0, kRB24IpSTTRo, kRB24IpSTCRo);
852 z0 += kRB24IpSTCL2;
853 shRB24IpSTC->DefineSection(3, z0, kRB24IpSTTRo, kRB24IpSTCRo);
854 z0 += kRB24IpSTCL3;
855 shRB24IpSTC->DefineSection(4, z0, kRB24IpSTTRo, kRB24IpSTTRo + 0.001);
856 TGeoVolume* voRB24IpSTC = new TGeoVolume("RB24IpSTC", shRB24IpSTC, kMedSteel);
857
858 // Pos 3 Shell LHCVC2U_0007
859 // Pos 4 Extruded Shell LHCVC2U_0006
860 Float_t kRB24IpShellL = 4.45; // Length of the Shell
861 Float_t kRB24IpShellD = 0.10; // Wall thickness of the shell
862 Float_t kRB24IpShellCTRi = 6.70/2.; // Inner radius of the connection tube
863 Float_t kRB24IpShellCTL = 1.56; // Length of the connection tube
864 Float_t kRB24IpShellCARi = 17.80/2.; // Inner radius of the cavity
865 Float_t kRB24IpShellCCRo = 18.20/2.; // Inner radius at the centre
866
867 TGeoPcon* shRB24IpShell = new TGeoPcon(0., 360., 7);
868 z0 = 0;
869 shRB24IpShell->DefineSection(0, z0, kRB24IpShellCTRi, kRB24IpShellCTRi + kRB24IpShellD);
870 z0 += kRB24IpShellCTL;
871 shRB24IpShell->DefineSection(1, z0, kRB24IpShellCTRi, kRB24IpShellCTRi + kRB24IpShellD);
872 shRB24IpShell->DefineSection(2, z0, kRB24IpShellCTRi, kRB24IpShellCARi + kRB24IpShellD);
873 z0 += kRB24IpShellD;
874 shRB24IpShell->DefineSection(3, z0, kRB24IpShellCARi, kRB24IpShellCARi + kRB24IpShellD);
875 z0 = kRB24IpShellL - kRB24IpShellD;
876 shRB24IpShell->DefineSection(4, z0, kRB24IpShellCARi, kRB24IpShellCARi + kRB24IpShellD);
877 shRB24IpShell->DefineSection(5, z0, kRB24IpShellCARi, kRB24IpShellCCRo);
878 z0 = kRB24IpShellL;
879 shRB24IpShell->DefineSection(6, z0, kRB24IpShellCARi, kRB24IpShellCCRo);
880 TGeoVolume* voRB24IpShell = new TGeoVolume("RB24IpShell", shRB24IpShell, kMedSteel);
881
882 TGeoPcon* shRB24IpShellM = MakeMotherFromTemplate(shRB24IpShell, 0, 6, kRB24IpShellCTRi , 13);
883
884
885 for (Int_t i = 0; i < 6; i++) {
886 z = 2. * kRB24IpShellL - shRB24IpShellM->GetZ(5-i);
887 Float_t rmin = shRB24IpShellM->GetRmin(5-i);
888 Float_t rmax = shRB24IpShellM->GetRmax(5-i);
889 shRB24IpShellM->DefineSection(7+i, z, rmin, rmax);
890 }
891
892 TGeoVolume* voRB24IpShellM = new TGeoVolume("RB24IpShellM", shRB24IpShellM, kMedVac);
893 voRB24IpShellM->SetVisibility(0);
894 voRB24IpShellM->AddNode(voRB24IpShell, 1, gGeoIdentity);
895 voRB24IpShellM->AddNode(voRB24IpShell, 2, new TGeoCombiTrans(0., 0., 2. * kRB24IpShellL, rot180));
896 //
897 // Pos 8 Pumping Elements
898 //
899 // Anode array
900 TGeoVolume* voRB24IpPE = new TGeoVolume("voRB24IpPE", new TGeoTube(0.9, 1., 2.54/2.), kMedSteel);
901 Float_t kRB24IpPEAR = 5.5;
902
903 for (Int_t i = 0; i < 15; i++) {
904 Float_t phi = Float_t(i) * 24.;
905 Float_t x = kRB24IpPEAR * TMath::Cos(kDegRad * phi);
906 Float_t y = kRB24IpPEAR * TMath::Sin(kDegRad * phi);
907 voRB24IpShellM->AddNode(voRB24IpPE, i+1, new TGeoTranslation(x, y, kRB24IpShellL));
908 }
909
910
911 //
912 // Cathodes
913 //
914 // Here we could add some Ti strips
915
916 // Postioning of elements
917 voRB24AIpM->AddNode(voRB24IpRF, 1, new TGeoTranslation(0., 0., -kRB24AIpML/2.));
918 voRB24AIpM->AddNode(voRB24IpRF, 2, new TGeoCombiTrans (0., 0., +kRB24AIpML/2., rot180));
919 voRB24AIpM->AddNode(voRB24IpSTT, 1, new TGeoTranslation(0., 0., 0.));
920 voRB24AIpM->AddNode(voRB24IpSTV, 1, new TGeoTranslation(0., 0., 0.));
921 voRB24AIpM->AddNode(voRB24IpShellM, 1, new TGeoTranslation(0., 0., -kRB24AIpML/2. + 8.13));
922 voRB24AIpM->AddNode(voRB24IpSTC, 1, new TGeoTranslation(0., 0., 8.13 - kRB24AIpML/2.));
923 voRB24AIpM->AddNode(voRB24IpSTC, 2, new TGeoCombiTrans (0., 0., 8.14 + 8.9 - kRB24AIpML/2., rot180));
924
925 //
926 // Valve
927 // VAC Series 47 DN 63 with manual actuator
928 //
929 const Float_t kRB24ValveWz = 7.5;
930 const Float_t kRB24ValveDN = 10.0/2.;
931 //
932 // Body containing the valve plate
933 //
934 const Float_t kRB24ValveBoWx = 15.6;
935 const Float_t kRB24ValveBoWy = (21.5 + 23.1 - 5.);
936 const Float_t kRB24ValveBoWz = 4.6;
937 const Float_t kRB24ValveBoD = 0.5;
938
939 TGeoVolume* voRB24ValveBoM =
940 new TGeoVolume("RB24ValveBoM",
941 new TGeoBBox( kRB24ValveBoWx/2., kRB24ValveBoWy/2., kRB24ValveBoWz/2.), kMedAir);
942 voRB24ValveBoM->SetVisibility(0);
943 TGeoVolume* voRB24ValveBo =
944 new TGeoVolume("RB24ValveBo",
945 new TGeoBBox( kRB24ValveBoWx/2., kRB24ValveBoWy/2., kRB24ValveBoWz/2.), kMedSteel);
946 voRB24ValveBoM->AddNode(voRB24ValveBo, 1, gGeoIdentity);
947 //
948 // Inner volume
949 //
950 TGeoVolume* voRB24ValveBoI = new TGeoVolume("RB24ValveBoI",
951 new TGeoBBox( kRB24ValveBoWx/2. - kRB24ValveBoD,
952 kRB24ValveBoWy/2. - kRB24ValveBoD/2.,
953 kRB24ValveBoWz/2. - kRB24ValveBoD),
954 kMedVac);
955 voRB24ValveBo->AddNode(voRB24ValveBoI, 1, new TGeoTranslation(0., kRB24ValveBoD/2., 0.));
956 //
957 // Opening and Flanges
958 const Float_t kRB24ValveFlRo = 18./2.;
959 const Float_t kRB24ValveFlD = 1.45;
960 TGeoVolume* voRB24ValveBoA = new TGeoVolume("RB24ValveBoA",
961 new TGeoTube(0., kRB24ValveDN/2., kRB24ValveBoD/2.), kMedVac);
962 voRB24ValveBo->AddNode(voRB24ValveBoA, 1, new TGeoTranslation(0., - kRB24ValveBoWy/2. + 21.5, -kRB24ValveBoWz/2. + kRB24ValveBoD/2.));
963 voRB24ValveBo->AddNode(voRB24ValveBoA, 2, new TGeoTranslation(0., - kRB24ValveBoWy/2. + 21.5, +kRB24ValveBoWz/2. - kRB24ValveBoD/2.));
964
965 TGeoVolume* voRB24ValveFl = new TGeoVolume("RB24ValveFl", new TGeoTube(kRB24ValveDN/2., kRB24ValveFlRo, kRB24ValveFlD/2.), kMedSteel);
966 TGeoVolume* voRB24ValveFlI = new TGeoVolume("RB24ValveFlI", new TGeoTube(0., kRB24ValveFlRo, kRB24ValveFlD/2.), kMedVac);
967 voRB24ValveFlI->AddNode(voRB24ValveFl, 1, gGeoIdentity);
968
969 //
970 // Actuator Flange
971 const Float_t kRB24ValveAFlWx = 18.9;
972 const Float_t kRB24ValveAFlWy = 5.0;
973 const Float_t kRB24ValveAFlWz = 7.7;
974 TGeoVolume* voRB24ValveAFl = new TGeoVolume("RB24ValveAFl", new TGeoBBox(kRB24ValveAFlWx/2., kRB24ValveAFlWy/2., kRB24ValveAFlWz/2.), kMedSteel);
975 //
976 // Actuator Tube
977 const Float_t kRB24ValveATRo = 9.7/2.;
978 const Float_t kRB24ValveATH = 16.6;
979 TGeoVolume* voRB24ValveAT = new TGeoVolume("RB24ValveAT", new TGeoTube(kRB24ValveATRo - 2. * kRB24ValveBoD,kRB24ValveATRo, kRB24ValveATH/2.),
980 kMedSteel);
981 //
982 // Manual Actuator (my best guess)
983 TGeoVolume* voRB24ValveMA1 = new TGeoVolume("RB24ValveMA1", new TGeoCone(2.5/2., 0., 0.5, 4.5, 5.), kMedSteel);
984 TGeoVolume* voRB24ValveMA2 = new TGeoVolume("RB24ValveMA2", new TGeoTorus(5., 0., 1.25), kMedSteel);
985 TGeoVolume* voRB24ValveMA3 = new TGeoVolume("RB24ValveMA3", new TGeoTube (0., 1.25, 2.5), kMedSteel);
986
987
988 //
989 // Position all volumes
990 Float_t y0;
991 TGeoVolumeAssembly* voRB24ValveMo = new TGeoVolumeAssembly("RB24ValveMo");
992 voRB24ValveMo->AddNode(voRB24ValveFl, 1, new TGeoTranslation(0., 0., - 7.5/2. + kRB24ValveFlD/2.));
993 voRB24ValveMo->AddNode(voRB24ValveFl, 2, new TGeoTranslation(0., 0., + 7.5/2. - kRB24ValveFlD/2.));
994 y0 = -21.5;
995 voRB24ValveMo->AddNode(voRB24ValveBoM, 1, new TGeoTranslation(0., y0 + kRB24ValveBoWy/2., 0.));
996 y0 += kRB24ValveBoWy;
997 voRB24ValveMo->AddNode(voRB24ValveAFl, 1, new TGeoTranslation(0., y0 + kRB24ValveAFlWy/2., 0.));
998 y0 += kRB24ValveAFlWy;
999 voRB24ValveMo->AddNode(voRB24ValveAT, 1, new TGeoCombiTrans(0., y0 + kRB24ValveATH/2., 0., rotyz));
1000 y0 += kRB24ValveATH;
1001 voRB24ValveMo->AddNode(voRB24ValveMA1, 1, new TGeoCombiTrans(0., y0 + 2.5/2., 0., rotyz));
1002 y0 += 2.5;
1003 voRB24ValveMo->AddNode(voRB24ValveMA2, 1, new TGeoCombiTrans(0., y0 + 2.5/2., 0., rotyz));
1004 y0 += 2.5;
1005 voRB24ValveMo->AddNode(voRB24ValveMA3, 1, new TGeoCombiTrans(5./TMath::Sqrt(2.), y0 + 5.0/2., 5./TMath::Sqrt(2.), rotyz));
1006 //
1007 // Warm Module Type VMABC
1008 // LHCVMABC_0002
1009 //
1010 //
1011 //
1012 // Flange 1.00
1013 // Central Piece 11.50
1014 // Bellow 14.50
1015 // End Flange 1.00
1016 //===================================
1017 // Total 28.00
1018 //
1019 // Pos 1 Warm Bellows DN100 LHCVBU__0016
1020 // Pos 2 Trans. Tube Flange LHCVSR__0062
1021 // Pos 3 RF Contact D63 LHCVSR__0057
1022 // [Pos 4 Hex. Countersunk Screw Bossard BN4719]
1023 // [Pos 5 Tension spring LHCVSR__00239]
1024 //
1025
1026 // Pos 1 Warm Bellows DN100 LHCVBU__0016
1027 // Pos 1.1 Right Body 2 Ports with Support LHCVBU__0014
1028 //
1029 // Tube 1
1030 const Float_t kRB24VMABCRBT1Ri = 10.0/2.;
1031 const Float_t kRB24VMABCRBT1Ro = 10.3/2.;
1032 const Float_t kRB24VMABCRBT1L = 11.5;
1033 const Float_t kRB24VMABCRBT1L2 = 8.;
1034 const Float_t kRB24VMABCL = 28.;
1035
1036 TGeoTube* shRB24VMABCRBT1 = new TGeoTube(kRB24VMABCRBT1Ri, kRB24VMABCRBT1Ro, kRB24VMABCRBT1L/2.);
1037 shRB24VMABCRBT1->SetName("RB24VMABCRBT1");
1038 TGeoTube* shRB24VMABCRBT1o = new TGeoTube(0., kRB24VMABCRBT1Ro, kRB24VMABCRBT1L/2.);
1039 shRB24VMABCRBT1o->SetName("RB24VMABCRBT1o");
1040 TGeoTube* shRB24VMABCRBT1o2 = new TGeoTube(0., kRB24VMABCRBT1Ro + 0.3, kRB24VMABCRBT1L/2.);
1041 shRB24VMABCRBT1o2->SetName("RB24VMABCRBT1o2");
1042 // Lower inforcement
1043 TGeoVolume* voRB24VMABCRBT12 = new TGeoVolume("RB24VMABCRBT12",
1044 new TGeoTubeSeg(kRB24VMABCRBT1Ro, kRB24VMABCRBT1Ro + 0.3, kRB24VMABCRBT1L2/2., 220., 320.)
1045 , kMedSteel);
1046 //
1047 // Tube 2
1048 const Float_t kRB24VMABCRBT2Ri = 6.0/2.;
1049 const Float_t kRB24VMABCRBT2Ro = 6.3/2.;
1050 const Float_t kRB24VMABCRBF2Ro = 11.4/2.;
1051 const Float_t kRB24VMABCRBT2L = 5.95 + 2.; // 2. cm added for welding
1052 const Float_t kRB24VMABCRBF2L = 1.75;
1053 TGeoTube* shRB24VMABCRBT2 = new TGeoTube(kRB24VMABCRBT2Ri, kRB24VMABCRBT2Ro, kRB24VMABCRBT2L/2.);
1054 shRB24VMABCRBT2->SetName("RB24VMABCRBT2");
1055 TGeoTube* shRB24VMABCRBT2i = new TGeoTube(0., kRB24VMABCRBT2Ri, kRB24VMABCRBT2L/2. + 2.);
1056 shRB24VMABCRBT2i->SetName("RB24VMABCRBT2i");
1057 TGeoCombiTrans* tRBT2 = new TGeoCombiTrans(-11.5 + kRB24VMABCRBT2L/2., 0., 7.2 - kRB24VMABCRBT1L/2. , rotxz);
1058 tRBT2->SetName("tRBT2");
1059 tRBT2->RegisterYourself();
1060 TGeoCompositeShape* shRB24VMABCRBT2c = new TGeoCompositeShape("shRB24VMABCRBT2c","RB24VMABCRBT2:tRBT2-RB24VMABCRBT1o");
1061 TGeoVolume* voRB24VMABCRBT2 = new TGeoVolume("shRB24VMABCRBT2", shRB24VMABCRBT2c, kMedSteel);
1062 // Flange
1063 // Pos 1.4 Flange DN63 LHCVBU__0008
1064 TGeoVolume* voRB24VMABCRBF2 = new TGeoVolume("RB24VMABCRBF2",
1065 new TGeoTube(kRB24VMABCRBT2Ro, kRB24VMABCRBF2Ro, kRB24VMABCRBF2L/2.), kMedSteel);
1066 // DN63 Blank Flange (my best guess)
1067 TGeoVolume* voRB24VMABCRBF2B = new TGeoVolume("RB24VMABCRBF2B",
1068 new TGeoTube(0., kRB24VMABCRBF2Ro, kRB24VMABCRBF2L/2.), kMedSteel);
1069 //
1070 // Tube 3
1071 const Float_t kRB24VMABCRBT3Ri = 3.5/2.;
1072 const Float_t kRB24VMABCRBT3Ro = 3.8/2.;
1073 const Float_t kRB24VMABCRBF3Ro = 7.0/2.;
1074 const Float_t kRB24VMABCRBT3L = 4.95 + 2.; // 2. cm added for welding
1075 const Float_t kRB24VMABCRBF3L = 1.27;
1076 TGeoTube* shRB24VMABCRBT3 = new TGeoTube(kRB24VMABCRBT3Ri, kRB24VMABCRBT3Ro, kRB24VMABCRBT3L/2);
1077 shRB24VMABCRBT3->SetName("RB24VMABCRBT3");
1078 TGeoTube* shRB24VMABCRBT3i = new TGeoTube(0., kRB24VMABCRBT3Ri, kRB24VMABCRBT3L/2. + 2.);
1079 shRB24VMABCRBT3i->SetName("RB24VMABCRBT3i");
1080 TGeoCombiTrans* tRBT3 = new TGeoCombiTrans(0., 10.5 - kRB24VMABCRBT3L/2., 7.2 - kRB24VMABCRBT1L/2. , rotyz);
1081 tRBT3->SetName("tRBT3");
1082 tRBT3->RegisterYourself();
1083 TGeoCompositeShape* shRB24VMABCRBT3c = new TGeoCompositeShape("shRB24VMABCRBT3c","RB24VMABCRBT3:tRBT3-RB24VMABCRBT1o");
1084 TGeoVolume* voRB24VMABCRBT3 = new TGeoVolume("shRB24VMABCRBT3", shRB24VMABCRBT3c, kMedSteel);
1085 // Flange
1086 // Pos 1.4 Flange DN35 LHCVBU__0007
1087 TGeoVolume* voRB24VMABCRBF3 = new TGeoVolume("RB24VMABCRBF3",
1088 new TGeoTube(kRB24VMABCRBT3Ro, kRB24VMABCRBF3Ro, kRB24VMABCRBF3L/2.), kMedSteel);
1089 //
1090 // Tube 4
1091 const Float_t kRB24VMABCRBT4Ri = 6.0/2.;
1092 const Float_t kRB24VMABCRBT4Ro = 6.4/2.;
1093 const Float_t kRB24VMABCRBT4L = 6.6;
1094 TGeoTube* shRB24VMABCRBT4 = new TGeoTube(kRB24VMABCRBT4Ri, kRB24VMABCRBT4Ro, kRB24VMABCRBT4L/2.);
1095 shRB24VMABCRBT4->SetName("RB24VMABCRBT4");
1096 TGeoCombiTrans* tRBT4 = new TGeoCombiTrans(0.,-11.+kRB24VMABCRBT4L/2., 7.2 - kRB24VMABCRBT1L/2. , rotyz);
1097 tRBT4->SetName("tRBT4");
1098 tRBT4->RegisterYourself();
1099 TGeoCompositeShape* shRB24VMABCRBT4c = new TGeoCompositeShape("shRB24VMABCRBT4c","RB24VMABCRBT4:tRBT4-RB24VMABCRBT1o2");
1100 TGeoVolume* voRB24VMABCRBT4 = new TGeoVolume("shRB24VMABCRBT4", shRB24VMABCRBT4c, kMedSteel);
1101 TGeoCompositeShape* shRB24VMABCRB = new TGeoCompositeShape("shRB24VMABCRB", "RB24VMABCRBT1-(RB24VMABCRBT2i:tRBT2+RB24VMABCRBT3i:tRBT3)");
1102 TGeoVolume* voRB24VMABCRBI = new TGeoVolume("RB24VMABCRBI", shRB24VMABCRB, kMedSteel);
1103 //
1104 // Plate
1105 const Float_t kRB24VMABCRBBx = 16.0;
1106 const Float_t kRB24VMABCRBBy = 1.5;
1107 const Float_t kRB24VMABCRBBz = 15.0;
1108
1109 // Relative position of tubes
1110 const Float_t kRB24VMABCTz = 7.2;
1111 // Relative position of plate
1112 const Float_t kRB24VMABCPz = 3.6;
1113 const Float_t kRB24VMABCPy = -12.5;
1114
1115 TGeoVolume* voRB24VMABCRBP = new TGeoVolume("RB24VMABCRBP", new TGeoBBox(kRB24VMABCRBBx/2., kRB24VMABCRBBy/2., kRB24VMABCRBBz/2.), kMedSteel);
1116 //
1117 // Pirani Gauge (my best guess)
1118 //
1119 TGeoPcon* shRB24VMABCPirani = new TGeoPcon(0., 360., 15);
1120 // DN35/16 Coupling
1121 z = 0;
1122 shRB24VMABCPirani->DefineSection( 0, z, 0.8 , kRB24VMABCRBF3Ro);
1123 z += kRB24VMABCRBF3L; // 1.3
1124 shRB24VMABCPirani->DefineSection( 1, z, 0.8 , kRB24VMABCRBF3Ro);
1125 shRB24VMABCPirani->DefineSection( 2, z, 0.8 , 1.0);
1126 // Pipe
1127 z += 2.8;
1128 shRB24VMABCPirani->DefineSection( 3, z, 0.8 , 1.0);
1129 // Flange
1130 shRB24VMABCPirani->DefineSection( 4, z, 0.8 , 1.75);
1131 z += 1.6;
1132 shRB24VMABCPirani->DefineSection( 5, z, 0.8 , 1.75);
1133 shRB24VMABCPirani->DefineSection( 6, z, 0.8 , 1.0);
1134 z += 5.2;
1135 shRB24VMABCPirani->DefineSection( 7, z, 0.8 , 1.0);
1136 shRB24VMABCPirani->DefineSection( 8, z, 0.8 , 2.5);
1137 z += 2.0;
1138 shRB24VMABCPirani->DefineSection( 9, z, 0.80, 2.50);
1139 shRB24VMABCPirani->DefineSection(10, z, 1.55, 1.75);
1140 z += 5.7;
1141 shRB24VMABCPirani->DefineSection(11, z, 1.55, 1.75);
1142 shRB24VMABCPirani->DefineSection(11, z, 0.00, 1.75);
1143 z += 0.2;
1144 shRB24VMABCPirani->DefineSection(12, z, 0.00, 1.75);
1145 shRB24VMABCPirani->DefineSection(13, z, 0.00, 0.75);
1146 z += 0.5;
1147 shRB24VMABCPirani->DefineSection(14, z, 0.00, 0.75);
1148 TGeoVolume* voRB24VMABCPirani = new TGeoVolume("RB24VMABCPirani", shRB24VMABCPirani, kMedSteel);
1149 //
1150 //
1151 //
1152
1153
1154 //
1155 // Positioning of elements
1156 TGeoVolumeAssembly* voRB24VMABCRB = new TGeoVolumeAssembly("RB24VMABCRB");
1157 //
1158 voRB24VMABCRB->AddNode(voRB24VMABCRBI, 1, gGeoIdentity);
1159 // Plate
1160 voRB24VMABCRB->AddNode(voRB24VMABCRBP, 1, new TGeoTranslation(0., kRB24VMABCPy + kRB24VMABCRBBy /2.,
1161 kRB24VMABCRBBz/2. - kRB24VMABCRBT1L/2. + kRB24VMABCPz));
1162 // Tube 2
1163 voRB24VMABCRB->AddNode(voRB24VMABCRBT2, 1, gGeoIdentity);
1164 // Flange Tube 2
1165 voRB24VMABCRB->AddNode(voRB24VMABCRBF2, 1, new TGeoCombiTrans(kRB24VMABCPy + kRB24VMABCRBF2L/2., 0., kRB24VMABCTz - kRB24VMABCRBT1L/2., rotxz));
1166 // Blank Flange Tube 2
1167 voRB24VMABCRB->AddNode(voRB24VMABCRBF2B, 1, new TGeoCombiTrans(kRB24VMABCPy- kRB24VMABCRBF2L/2., 0., kRB24VMABCTz - kRB24VMABCRBT1L/2., rotxz));
1168 // Tube 3
1169 voRB24VMABCRB->AddNode(voRB24VMABCRBT3, 1, gGeoIdentity);
1170 // Flange Tube 3
1171 voRB24VMABCRB->AddNode(voRB24VMABCRBF3, 1, new TGeoCombiTrans(0., 11.2 - kRB24VMABCRBF3L/2., kRB24VMABCTz - kRB24VMABCRBT1L/2., rotyz));
1172 // Pirani Gauge
1173 voRB24VMABCRB->AddNode(voRB24VMABCPirani, 1, new TGeoCombiTrans(0., 11.2, kRB24VMABCTz - kRB24VMABCRBT1L/2., rotyz));
1174 // Tube 4
1175 voRB24VMABCRB->AddNode(voRB24VMABCRBT4, 1, gGeoIdentity);
1176 // Inforcement
1177 voRB24VMABCRB->AddNode(voRB24VMABCRBT12, 1, new TGeoTranslation(0., 0., kRB24VMABCRBT1L2/2. - kRB24VMABCRBT1L/2. + 2.8));
1178
1179
1180 // Pos 1.3 Bellows with end part LHCVBU__0002
1181 //
1182 // Connection Tube
1183 // Connection tube inner r
1184 const Float_t kRB24VMABBEConTubeRin = 10.0/2.;
1185 // Connection tube outer r
1186 const Float_t kRB24VMABBEConTubeRou = 10.3/2.;
1187 // Connection tube length
1188 const Float_t kRB24VMABBEConTubeL1 = 0.9;
1189 const Float_t kRB24VMABBEConTubeL2 = 2.6;
1190 // const Float_t RB24VMABBEBellowL = kRB24VMABBEConTubeL1 + kRB24VMABBEConTubeL2 + kRB24B1BellowUndL;
1191
1192 // Mother volume
1193 TGeoPcon* shRB24VMABBEBellowM = new TGeoPcon(0., 360., 6);
1194 // Connection Tube and Flange
1195 z = 0.;
1196 shRB24VMABBEBellowM->DefineSection( 0, z, kRB24VMABBEConTubeRin, kRB24VMABBEConTubeRou);
1197 z += kRB24VMABBEConTubeL1;
1198 shRB24VMABBEBellowM->DefineSection( 1, z, kRB24VMABBEConTubeRin, kRB24VMABBEConTubeRou);
1199 shRB24VMABBEBellowM->DefineSection( 2, z, kRB24B1BellowRi, kRB24B1BellowRo + kRB24B1ProtTubeThickness);
1200 z += kRB24B1BellowUndL;
1201 shRB24VMABBEBellowM->DefineSection( 3, z, kRB24B1BellowRi, kRB24B1BellowRo + kRB24B1ProtTubeThickness);
1202 shRB24VMABBEBellowM->DefineSection( 4, z, kRB24VMABBEConTubeRin, kRB24VMABBEConTubeRou);
1203 z += kRB24VMABBEConTubeL2;
1204 shRB24VMABBEBellowM->DefineSection( 5, z, kRB24VMABBEConTubeRin, kRB24VMABBEConTubeRou);
1205 TGeoVolume* voRB24VMABBEBellowM = new TGeoVolume("RB24VMABBEBellowM", shRB24VMABBEBellowM, kMedVac);
1206 voRB24VMABBEBellowM->SetVisibility(0);
1207
1208 // Connection tube left
1209 TGeoVolume* voRB24VMABBECT1 = new TGeoVolume("RB24VMABBECT1",
1210 new TGeoTube(kRB24VMABBEConTubeRin, kRB24VMABBEConTubeRou,kRB24VMABBEConTubeL1/2.),
1211 kMedSteel);
1212 // Connection tube right
1213 TGeoVolume* voRB24VMABBECT2 = new TGeoVolume("RB24VMABBECT2",
1214 new TGeoTube(kRB24VMABBEConTubeRin, kRB24VMABBEConTubeRou,kRB24VMABBEConTubeL2/2.),
1215 kMedSteel);
1216 z = kRB24VMABBEConTubeL1/2.;
1217 voRB24VMABBEBellowM->AddNode(voRB24VMABBECT1, 1, new TGeoTranslation(0., 0., z));
1218 z += kRB24VMABBEConTubeL1/2.;
1219 z += kRB24B1BellowUndL/2.;
1220 voRB24VMABBEBellowM->AddNode(voRB24B1Bellow, 2, new TGeoTranslation(0., 0., z));
1221 z += kRB24B1BellowUndL/2.;
1222 z += kRB24VMABBEConTubeL2/2.;
1223 voRB24VMABBEBellowM->AddNode(voRB24VMABBECT2, 1, new TGeoTranslation(0., 0., z));
1224 z += kRB24VMABBEConTubeL2/2.;
1225
1226 voRB24VMABCRB->AddNode(voRB24VMABBEBellowM, 1, new TGeoTranslation(0., 0., kRB24VMABCRBT1L/2.));
1227
1228 // Pos 1.2 Rotable flange LHCVBU__0013[*]
1229 // Front
1230 voRB24VMABCRB->AddNode(voRB24B1RFlange, 3, new TGeoCombiTrans(0., 0., - kRB24VMABCRBT1L/2. + 0.86, rot180));
1231 // End
1232 z = kRB24VMABCRBT1L/2. + kRB24B1BellowUndL +kRB24VMABBEConTubeL1 + kRB24VMABBEConTubeL2;
1233 voRB24VMABCRB->AddNode(voRB24B1RFlange, 4, new TGeoTranslation(0., 0., z - 0.86));
1234
1235
1236 // Pos 2 Trans. Tube Flange LHCVSR__0062
1237 // Pos 2.1 Transition Tube LHCVSR__0063
1238 // Pos 2.2 Transition Flange LHCVSR__0060
1239 //
1240 // Transition Tube with Flange
1241 TGeoPcon* shRB24VMABCTT = new TGeoPcon(0., 360., 7);
1242 z = 0.;
1243 shRB24VMABCTT->DefineSection(0, z, 6.3/2., 11.16/2.);
1244 z += 0.25;
1245 shRB24VMABCTT->DefineSection(1, z, 6.3/2., 11.16/2.);
1246 shRB24VMABCTT->DefineSection(2, z, 6.3/2., 9.30/2.);
1247 z += 0.25;
1248 shRB24VMABCTT->DefineSection(3, z, 6.3/2., 9.30/2.);
1249 shRB24VMABCTT->DefineSection(4, z, 6.3/2., 6.70/2.);
1250 z += (20.35 - 0.63);
1251 shRB24VMABCTT->DefineSection(5, z, 6.3/2., 6.7/2.);
1252 z += 0.63;
1253 shRB24VMABCTT->DefineSection(6, z, 6.3/2., 6.7/2.);
1254 TGeoVolume* voRB24VMABCTT = new TGeoVolume("RB24VMABCTT", shRB24VMABCTT, kMedSteel);
1255 voRB24VMABCRB->AddNode(voRB24VMABCTT, 1, new TGeoTranslation(0., 0., - kRB24VMABCRBT1L/2.-1.));
1256
1257 // Pos 3 RF Contact D63 LHCVSR__0057
1258 // Pos 3.1 RF Contact Flange LHCVSR__0017
1259 //
1260 TGeoPcon* shRB24VMABCCTFlange = new TGeoPcon(0., 360., 6);
1261 const Float_t kRB24VMABCCTFlangeRin = 6.36/2.; // Inner radius
1262 const Float_t kRB24VMABCCTFlangeL = 1.30; // Length
1263
1264 z = 0.;
1265 shRB24VMABCCTFlange->DefineSection(0, z, kRB24VMABCCTFlangeRin, 6.5/2.);
1266 z += 0.15;
1267 shRB24VMABCCTFlange->DefineSection(1, z, kRB24VMABCCTFlangeRin, 6.5/2.);
1268 shRB24VMABCCTFlange->DefineSection(2, z, kRB24VMABCCTFlangeRin, 6.9/2.);
1269 z += 0.9;
1270 shRB24VMABCCTFlange->DefineSection(3, z, kRB24VMABCCTFlangeRin, 6.9/2.);
1271 shRB24VMABCCTFlange->DefineSection(4, z, kRB24VMABCCTFlangeRin, 11.16/2.);
1272 z += 0.25;
1273 shRB24VMABCCTFlange->DefineSection(5, z, kRB24VMABCCTFlangeRin, 11.16/2.);
1274 TGeoVolume* voRB24VMABCCTFlange = new TGeoVolume("RB24VMABCCTFlange", shRB24VMABCCTFlange, kMedCu);
1275 //
1276 // Pos 3.2 RF-Contact LHCVSR__0056
1277 //
1278 TGeoPcon* shRB24VMABCCT = new TGeoPcon(0., 360., 4);
1279 const Float_t kRB24VMABCCTRin = 6.30/2.; // Inner radius
1280 const Float_t kRB24VMABCCTCRin = 7.29/2.; // Max. inner radius conical section
1281 const Float_t kRB24VMABCCTL = 11.88; // Length
1282 const Float_t kRB24VMABCCTSL = 10.48; // Length of straight section
1283 const Float_t kRB24VMABCCTd = 0.03; // Thickness
1284 z = 0;
1285 shRB24VMABCCT->DefineSection(0, z, kRB24VMABCCTCRin, kRB24VMABCCTCRin + kRB24VMABCCTd);
1286 z = kRB24VMABCCTL - kRB24VMABCCTSL;
1287 shRB24VMABCCT->DefineSection(1, z, kRB24VMABCCTRin + 0.35, kRB24VMABCCTRin + 0.35 + kRB24VMABCCTd);
1288 z = kRB24VMABCCTL - kRB24VMABCCTFlangeL;
1289 shRB24VMABCCT->DefineSection(2, z, kRB24VMABCCTRin, kRB24VMABCCTRin + kRB24VMABCCTd);
1290 z = kRB24VMABCCTL;
1291 shRB24VMABCCT->DefineSection(3, z, kRB24VMABCCTRin, kRB24VMABCCTRin + kRB24VMABCCTd);
1292
1293 TGeoVolume* voRB24VMABCCT = new TGeoVolume("RB24VMABCCT", shRB24VMABCCT, kMedCu);
1294
1295 TGeoVolumeAssembly* voRB24VMABRFCT = new TGeoVolumeAssembly("RB24VMABRFCT");
1296 voRB24VMABRFCT->AddNode(voRB24VMABCCT, 1, gGeoIdentity);
1297 voRB24VMABRFCT->AddNode( voRB24VMABCCTFlange, 1, new TGeoTranslation(0., 0., kRB24VMABCCTL - kRB24VMABCCTFlangeL));
1298
1299 z = kRB24VMABCRBT1L/2. + kRB24B1BellowUndL + kRB24VMABBEConTubeL1 + kRB24VMABBEConTubeL2 - kRB24VMABCCTL + 1.;
1300 voRB24VMABCRB->AddNode(voRB24VMABRFCT, 1, new TGeoTranslation(0., 0., z));
1301
1302
1303 //
1304 // Assembling RB24/1
1305 //
1306 TGeoVolumeAssembly* voRB24 = new TGeoVolumeAssembly("RB24");
1307 // Cu Tube with two simplified flanges
1308 voRB24->AddNode(voRB24CuTubeM, 1, gGeoIdentity);
1309 voRB24->AddNode(voRB24CuTubeA, 1, gGeoIdentity);
1310 z = - kRB24CuTubeL/2 + kRB24CuTubeFL/2.;
1311 voRB24->AddNode(voRB24CuTubeF, 1, new TGeoTranslation(0., 0., z));
1312 z = + kRB24CuTubeL/2 - kRB24CuTubeFL/2.;
1313 voRB24->AddNode(voRB24CuTubeF, 2, new TGeoTranslation(0., 0., z));
1314 // VMABC close to compensator magnet
1315 z = - kRB24CuTubeL/2. - (kRB24VMABCL - kRB24VMABCRBT1L/2) + 1.;
1316
1317 voRB24->AddNode(voRB24VMABCRB, 2, new TGeoTranslation(0., 0., z));
1318 // Bellow
1319 z = kRB24CuTubeL/2;
1320 voRB24->AddNode(voRB24B1BellowM, 1, new TGeoTranslation(0., 0., z));
1321 z += (kRB24B1L + kRB24AIpML/2.);
1322 // Annular ion pump
1323 voRB24->AddNode(voRB24AIpM, 1, new TGeoTranslation(0., 0., z));
1324 z += (kRB24AIpML/2. + kRB24ValveWz/2.);
1325 // Valve
1326 voRB24->AddNode(voRB24ValveMo, 1, new TGeoTranslation(0., 0., z));
1327 z += (kRB24ValveWz/2.+ kRB24VMABCRBT1L/2. + 1.);
1328 // VMABC close to forward detectors
1329 voRB24->AddNode(voRB24VMABCRB, 3, new TGeoTranslation(0., 0., z));
1330 //
1331 // RB24/2
1332 //
1333 // Copper Tube RB24/2
1334 const Float_t kRB242CuTubeL = 330.0;
1335
1336 TGeoVolume* voRB242CuTubeM = new TGeoVolume("voRB242CuTubeM",
1337 new TGeoTube(0., kRB24CuTubeRo, kRB242CuTubeL/2.), kMedVac);
1338 voRB24CuTubeM->SetVisibility(0);
1339 TGeoVolume* voRB242CuTube = new TGeoVolume("voRB242CuTube",
1340 new TGeoTube(kRB24CuTubeRi, kRB24CuTubeRo, kRB242CuTubeL/2.), kMedCu);
1341 voRB242CuTubeM->AddNode(voRB242CuTube, 1, gGeoIdentity);
1342
1343
1344 TGeoVolumeAssembly* voRB242 = new TGeoVolumeAssembly("RB242");
1345 voRB242->AddNode(voRB242CuTube, 1, gGeoIdentity);
1346 z = - kRB242CuTubeL/2 + kRB24CuTubeFL/2.;
1347 voRB242->AddNode(voRB24CuTubeF, 3, new TGeoTranslation(0., 0., z));
1348 z = + kRB242CuTubeL/2 - kRB24CuTubeFL/2.;
1349 voRB242->AddNode(voRB24CuTubeF, 4, new TGeoTranslation(0., 0., z));
1350 z = - kRB24CuTubeL/2 - kRB24VMABCL - kRB242CuTubeL/2.;
1351 voRB24->AddNode(voRB242, 1, new TGeoTranslation(0., 0., z));
1352 //
1353 // RB24/3
1354 //
1355 // Copper Tube RB24/3
1356 const Float_t kRB243CuTubeL = 303.35;
1357
1358 TGeoVolume* voRB243CuTubeM = new TGeoVolume("voRB243CuTubeM",
1359 new TGeoTube(0., kRB24CuTubeRo, kRB243CuTubeL/2.), kMedVac);
1360 voRB24CuTubeM->SetVisibility(0);
1361 TGeoVolume* voRB243CuTube = new TGeoVolume("voRB243CuTube",
1362 new TGeoTube(kRB24CuTubeRi, kRB24CuTubeRo, kRB243CuTubeL/2.), kMedCu);
1363 voRB243CuTubeM->AddNode(voRB243CuTube, 1, gGeoIdentity);
1364
1365
1366 TGeoVolumeAssembly* voRB243 = new TGeoVolumeAssembly("RB243");
1367 TGeoVolumeAssembly* voRB243A = new TGeoVolumeAssembly("RB243A");
1368
1369 voRB243A->AddNode(voRB243CuTube, 1, gGeoIdentity);
1370 z = - kRB243CuTubeL/2 + kRB24CuTubeFL/2.;
1371 voRB243A->AddNode(voRB24CuTubeF, 5, new TGeoTranslation(0., 0., z));
1372 z = + kRB243CuTubeL/2 - kRB24CuTubeFL/2.;
1373 voRB243A->AddNode(voRB24CuTubeF, 6, new TGeoTranslation(0., 0., z));
1374 z = + kRB243CuTubeL/2;
1375 voRB243A->AddNode(voRB24B1BellowM, 2, new TGeoTranslation(0., 0., z));
1376
1377 z = - kRB243CuTubeL/2. - kRB24B1L;
1378 voRB243->AddNode(voRB243A, 1, new TGeoTranslation(0., 0., z));
1379 z = - (1.5 * kRB243CuTubeL + 2. * kRB24B1L);
1380 voRB243->AddNode(voRB243A, 2, new TGeoTranslation(0., 0., z));
1381
1382 z = - 2. * (kRB243CuTubeL + kRB24B1L) - (kRB24VMABCL - kRB24VMABCRBT1L/2) + 1.;
1383 voRB243->AddNode(voRB24VMABCRB, 3, new TGeoTranslation(0., 0., z));
1384
1385 z = - kRB24CuTubeL/2 - kRB24VMABCL - kRB242CuTubeL;
1386 voRB24->AddNode(voRB243, 1, new TGeoTranslation(0., 0., z));
1387
1388
1389 //
1390 //
1391 top->AddNode(voRB24, 1, new TGeoCombiTrans(0., 0., kRB24CuTubeL/2 + 88.5 + 400., rot180));
1392
1393
1394 //
1395 ////////////////////////////////////////////////////////////////////////////////
1396 // //
1397 // The Absorber Vacuum system //
1398 // //
1399 ////////////////////////////////////////////////////////////////////////////////
1400 //
1401 // Rotable Flange starts at: 82.00 cm from IP
1402 // Length of rotable flange section: 10.68 cm
1403 // Weld 0.08 cm
1404 // Length of straight section 207.21 cm
1405 // =======================================================================
1406 // 299.97 cm [0.03 cm missing ?]
1407 // Length of opening cone 252.09 cm
1408 // Weld 0.15 cm
1409 // Length of compensator 30.54 cm
1410 // Weld 0.15 cm
1411 // Length of fixed flange 2.13 - 0.97 1.16 cm
1412 // =======================================================================
1413 // 584.06 cm [584.80 installed] [0.74 cm missing]
1414 // RB26/3
1415 // Length of split flange 2.13 - 1.2 0.93 cm
1416 // Weld 0.15 cm
1417 // Length of fixed point section 16.07 cm
1418 // Weld 0.15 cm
1419 // Length of opening cone 629.20 cm
1420 // Weld 0.30 cm
1421 // Kength of the compensator 41.70 cm
1422 // Weld 0.30 cm
1423 // Length of fixed flange 2.99 - 1.72 1.27 cm
1424 // =================================================
1425 // Length of RB26/3 690.07 cm [689.20 installed] [0.87 cm too much]
1426 //
1427 // RB26/4-5
1428 // Length of split flange 2.13 - 1.2 0.93 cm
1429 // Weld 0.15 cm
1430 // Length of fixed point section 16.07 cm
1431 // Weld 0.15 cm
1432 // Length of opening cone 629.20 cm
1433 // Weld 0.30 cm
1434 // Length of closing cone
1435 // Weld
1436 // Lenth of straight section
1437 // Kength of the compensator 41.70 cm
1438 // Weld 0.30 cm
1439 // Length of fixed flange 2.99 - 1.72 1.27 cm
1440 // =================================================
1441 // Length of RB26/3 690.07 cm [689.20 installed] [0.87 cm too much]
1442
1443 ///////////////////////////////////////////
1444 // //
1445 // RB26/1-2 //
1446 // Drawing LHCV2a_0050 [as installed] //
1447 // Drawing LHCV2a_0008 //
1448 // Drawing LHCV2a_0001 //
1449 ///////////////////////////////////////////
1450 // Pos1 Vacuum Tubes LHCVC2A__0010
1451 // Pos2 Compensator LHCVC2A__0064
1452 // Pos3 Rotable Flange LHCVFX___0016
1453 // Pos4 Fixed Flange LHCVFX___0006
1454 // Pos5 Bellow Tooling LHCVFX___0003
1455 //
1456 //
1457 //
1458 ///////////////////////////////////
1459 // RB26/1-2 Vacuum Tubes //
1460 // Drawing LHCVC2a_0010 //
1461 ///////////////////////////////////
1462 const Float_t kRB26s12TubeL = 459.45; // 0.15 cm added for welding
1463 //
1464 // Add 1 cm on outer diameter for insulation
1465 //
1466 TGeoPcon* shRB26s12Tube = new TGeoPcon(0., 360., 5);
1467 // Section 1: straight section
1468 shRB26s12Tube->DefineSection(0, 0.00, 5.84/2., 6.00/2.);
1469 shRB26s12Tube->DefineSection(1, 207.21, 5.84/2., 6.00/2.);
1470 // Section 2: 0.72 deg opening cone
1471 shRB26s12Tube->DefineSection(2, 207.21, 5.84/2., 6.14/2.);
1472 shRB26s12Tube->DefineSection(3, 452.30, 12.00/2., 12.30/2.);
1473 shRB26s12Tube->DefineSection(4, kRB26s12TubeL, 12.00/2., 12.30/2.);
1474 TGeoVolume* voRB26s12Tube = new TGeoVolume("RB26s12Tube", shRB26s12Tube, kMedSteel);
1475 // Add the insulation layer
1476 TGeoVolume* voRB26s12TubeIns = new TGeoVolume("RB26s12TubeIns", MakeInsulationFromTemplate(shRB26s12Tube), kMedInsu);
1477 voRB26s12Tube->AddNode(voRB26s12TubeIns, 1, gGeoIdentity);
1478
1479
1480 TGeoVolume* voRB26s12TubeM = new TGeoVolume("RB26s12TubeM", MakeMotherFromTemplate(shRB26s12Tube), kMedVac);
1481 voRB26s12TubeM->AddNode(voRB26s12Tube, 1, gGeoIdentity);
1482
1483
1484
1485 ///////////////////////////////////
1486 // RB26/2 Axial Compensator //
1487 // Drawing LHCVC2a_0064 //
1488 ///////////////////////////////////
1489 const Float_t kRB26s2CompL = 30.65; // Length of the compensator
1490 const Float_t kRB26s2BellowRo = 14.38/2.; // Bellow outer radius [Pos 1]
1491 const Float_t kRB26s2BellowRi = 12.12/2.; // Bellow inner radius [Pos 1]
1492 const Int_t kRB26s2NumberOfPlies = 14; // Number of plies [Pos 1]
1493 const Float_t kRB26s2BellowUndL = 10.00; // Length of undulated region [Pos 1] [+10 mm installed including pretension ?]
1494 const Float_t kRB26s2PlieThickness = 0.025; // Plie thickness [Pos 1]
1495 const Float_t kRB26s2ConnectionPlieR = 0.21; // Connection plie radius [Pos 1]
1496 // Plie radius
1497 const Float_t kRB26s2PlieR =
1498 (kRB26s2BellowUndL - 4. * kRB26s2ConnectionPlieR + 2. * kRB26s2PlieThickness +
1499 (2. * kRB26s2NumberOfPlies - 2.) * kRB26s2PlieThickness) / (4. * kRB26s2NumberOfPlies - 2.);
1500 const Float_t kRB26s2CompTubeInnerR = 12.00/2.; // Connection tubes inner radius [Pos 2 + 3]
1501 const Float_t kRB26s2CompTubeOuterR = 12.30/2.; // Connection tubes outer radius [Pos 2 + 3]
1502 const Float_t kRB26s2WeldingTubeLeftL = 9.00/2.; // Left connection tube half length [Pos 2]
1503 const Float_t kRB26s2WeldingTubeRightL = 11.65/2.; // Right connection tube half length [Pos 3] [+ 0.15 cm for welding]
1504 const Float_t kRB26s2RingOuterR = 18.10/2.; // Ring inner radius [Pos 4]
1505 const Float_t kRB26s2RingL = 0.40/2.; // Ring half length [Pos 4]
1506 const Float_t kRB26s2RingZ = 6.50 ; // Ring z-position [Pos 4]
1507 const Float_t kRB26s2ProtOuterR = 18.20/2.; // Protection tube outer radius [Pos 5]
1508 const Float_t kRB26s2ProtL = 15.00/2.; // Protection tube half length [Pos 5]
1509 const Float_t kRB26s2ProtZ = 6.70 ; // Protection tube z-position [Pos 5]
1510
1511
1512 // Mother volume
1513 //
1514 TGeoPcon* shRB26s2Compensator = new TGeoPcon(0., 360., 6);
1515 shRB26s2Compensator->DefineSection( 0, 0.0, 0., kRB26s2CompTubeOuterR);
1516 shRB26s2Compensator->DefineSection( 1, kRB26s2RingZ, 0., kRB26s2CompTubeOuterR);
1517 shRB26s2Compensator->DefineSection( 2, kRB26s2RingZ, 0., kRB26s2ProtOuterR);
1518 shRB26s2Compensator->DefineSection( 3, kRB26s2ProtZ + 2. * kRB26s2ProtL, 0., kRB26s2ProtOuterR);
1519 shRB26s2Compensator->DefineSection( 4, kRB26s2ProtZ + 2. * kRB26s2ProtL, 0., kRB26s2CompTubeOuterR);
1520 shRB26s2Compensator->DefineSection( 5, kRB26s2CompL , 0., kRB26s2CompTubeOuterR);
1521 TGeoVolume* voRB26s2Compensator = new TGeoVolume("RB26s2Compensator", shRB26s2Compensator, kMedVac);
1522
1523 //
1524 // [Pos 1] Bellow
1525 //
1526 //
1527 TGeoVolume* voRB26s2Bellow = new TGeoVolume("RB26s2Bellow", new TGeoTube(kRB26s2BellowRi, kRB26s2BellowRo, kRB26s2BellowUndL/2.), kMedVac);
1528 //
1529 // Upper part of the undulation
1530 //
1531 TGeoTorus* shRB26s2PlieTorusU = new TGeoTorus(kRB26s2BellowRo - kRB26s2PlieR, kRB26s2PlieR - kRB26s2PlieThickness, kRB26s2PlieR);
1532 shRB26s2PlieTorusU->SetName("RB26s2TorusU");
1533 TGeoTube* shRB26s2PlieTubeU = new TGeoTube (kRB26s2BellowRo - kRB26s2PlieR, kRB26s2BellowRo, kRB26s2PlieR);
1534 shRB26s2PlieTubeU->SetName("RB26s2TubeU");
1535 TGeoCompositeShape* shRB26s2UpperPlie = new TGeoCompositeShape("RB26s2UpperPlie", "RB26s2TorusU*RB26s2TubeU");
1536
1537 TGeoVolume* voRB26s2WiggleU = new TGeoVolume("RB26s2UpperPlie", shRB26s2UpperPlie, kMedSteel);
1538 //
1539 // Lower part of the undulation
1540 TGeoTorus* shRB26s2PlieTorusL = new TGeoTorus(kRB26s2BellowRi + kRB26s2PlieR, kRB26s2PlieR - kRB26s2PlieThickness, kRB26s2PlieR);
1541 shRB26s2PlieTorusL->SetName("RB26s2TorusL");
1542 TGeoTube* shRB26s2PlieTubeL = new TGeoTube (kRB26s2BellowRi, kRB26s2BellowRi + kRB26s2PlieR, kRB26s2PlieR);
1543 shRB26s2PlieTubeL->SetName("RB26s2TubeL");
1544 TGeoCompositeShape* shRB26s2LowerPlie = new TGeoCompositeShape("RB26s2LowerPlie", "RB26s2TorusL*RB26s2TubeL");
1545
1546 TGeoVolume* voRB26s2WiggleL = new TGeoVolume("RB26s2LowerPlie", shRB26s2LowerPlie, kMedSteel);
1547
1548 //
1549 // Connection between upper and lower part of undulation
1550 TGeoVolume* voRB26s2WiggleC1 = new TGeoVolume("RB26s2PlieConn1",
1551 new TGeoTube(kRB26s2BellowRi + kRB26s2PlieR,
1552 kRB26s2BellowRo - kRB26s2PlieR, kRB26s2PlieThickness / 2.), kMedSteel);
1553 //
1554 // One wiggle
1555 TGeoVolumeAssembly* voRB26s2Wiggle = new TGeoVolumeAssembly("RB26s2Wiggle");
1556 z0 = - kRB26s2PlieThickness / 2.;
1557 voRB26s2Wiggle->AddNode(voRB26s2WiggleC1, 1 , new TGeoTranslation(0., 0., z0));
1558 z0 += kRB26s2PlieR - kRB26s2PlieThickness / 2.;
1559 voRB26s2Wiggle->AddNode(voRB26s2WiggleU, 1 , new TGeoTranslation(0., 0., z0));
1560 z0 += kRB26s2PlieR - kRB26s2PlieThickness / 2.;
1561 voRB26s2Wiggle->AddNode(voRB26s2WiggleC1, 2 , new TGeoTranslation(0., 0., z0));
1562 z0 += kRB26s2PlieR - kRB26s2PlieThickness;
1563 voRB26s2Wiggle->AddNode(voRB26s2WiggleL , 1 , new TGeoTranslation(0., 0., z0));
1564 // Positioning of the volumes
1565 z0 = - kRB26s2BellowUndL/2.+ kRB26s2ConnectionPlieR;
1566 voRB26s2Bellow->AddNode(voRB26s2WiggleL, 1, new TGeoTranslation(0., 0., z0));
1567 z0 += kRB26s2ConnectionPlieR;
1568 zsh = 4. * kRB26s2PlieR - 2. * kRB26s2PlieThickness;
1569 for (Int_t iw = 0; iw < kRB26s2NumberOfPlies; iw++) {
1570 Float_t zpos = z0 + iw * zsh;
1571 voRB26s2Bellow->AddNode(voRB26s2Wiggle, iw + 1, new TGeoTranslation(0., 0., zpos - kRB26s2PlieThickness));
1572 }
1573
1574 voRB26s2Compensator->AddNode(voRB26s2Bellow, 1, new TGeoTranslation(0., 0., 2. * kRB26s2WeldingTubeLeftL + kRB26s2BellowUndL/2.));
1575
1576 //
1577 // [Pos 2] Left Welding Tube
1578 //
1579 TGeoTube* shRB26s2CompLeftTube = new TGeoTube(kRB26s2CompTubeInnerR, kRB26s2CompTubeOuterR, kRB26s2WeldingTubeLeftL);
1580 TGeoVolume* voRB26s2CompLeftTube = new TGeoVolume("RB26s2CompLeftTube", shRB26s2CompLeftTube, kMedSteel);
1581 voRB26s2Compensator->AddNode(voRB26s2CompLeftTube, 1, new TGeoTranslation(0., 0., kRB26s2WeldingTubeLeftL));
1582 //
1583 // [Pos 3] Right Welding Tube
1584 //
1585 TGeoTube* shRB26s2CompRightTube = new TGeoTube(kRB26s2CompTubeInnerR, kRB26s2CompTubeOuterR, kRB26s2WeldingTubeRightL);
1586 TGeoVolume* voRB26s2CompRightTube = new TGeoVolume("RB26s2CompRightTube", shRB26s2CompRightTube, kMedSteel);
1587 voRB26s2Compensator->AddNode(voRB26s2CompRightTube, 1, new TGeoTranslation(0., 0., kRB26s2CompL - kRB26s2WeldingTubeRightL));
1588 //
1589 // [Pos 4] Ring
1590 //
1591 TGeoTube* shRB26s2CompRing = new TGeoTube(kRB26s2CompTubeOuterR, kRB26s2RingOuterR, kRB26s2RingL);
1592 TGeoVolume* voRB26s2CompRing = new TGeoVolume("RB26s2CompRing", shRB26s2CompRing, kMedSteel);
1593 voRB26s2Compensator->AddNode(voRB26s2CompRing, 1, new TGeoTranslation(0., 0., kRB26s2RingZ + kRB26s2RingL));
1594
1595 //
1596 // [Pos 5] Outer Protecting Tube
1597 //
1598 TGeoTube* shRB26s2CompProtTube = new TGeoTube(kRB26s2RingOuterR, kRB26s2ProtOuterR, kRB26s2ProtL);
1599 TGeoVolume* voRB26s2CompProtTube = new TGeoVolume("RB26s2CompProtTube", shRB26s2CompProtTube, kMedSteel);
1600 voRB26s2Compensator->AddNode(voRB26s2CompProtTube, 1, new TGeoTranslation(0., 0., kRB26s2ProtZ + kRB26s2ProtL));
1601
1602 ///////////////////////////////////
1603 // Rotable Flange //
1604 // Drawing LHCVFX_0016 //
1605 ///////////////////////////////////
af5f3976 1606 const Float_t kRB26s1RFlangeTubeRi = 5.84/2.; // Tube inner radius
1607 const Float_t kRB26s1RFlangeTubeRo = 6.00/2.; // Tube outer radius
820b4d9e 1608
1609 // Pos 1 Clamp Ring LHCVFX__0015
1610 const Float_t kRB26s1RFlangeCrL = 1.40 ; // Lenth of the clamp ring
af5f3976 1611 const Float_t kRB26s1RFlangeCrRi1 = 6.72/2.; // Ring inner radius section 1
1612 const Float_t kRB26s1RFlangeCrRi2 = 6.06/2.; // Ring inner radius section 2
1613 const Float_t kRB26s1RFlangeCrRo = 8.60/2.;// Ring outer radius
820b4d9e 1614 const Float_t kRB26s1RFlangeCrD = 0.800 ; // Width section 1
1615
1616 TGeoPcon* shRB26s1RFlangeCr = new TGeoPcon(0., 360., 4);
1617 z0 = 0.;
1618 shRB26s1RFlangeCr->DefineSection(0, z0, kRB26s1RFlangeCrRi1, kRB26s1RFlangeCrRo);
1619 z0 += kRB26s1RFlangeCrD;
1620 shRB26s1RFlangeCr->DefineSection(1, z0, kRB26s1RFlangeCrRi1, kRB26s1RFlangeCrRo);
1621 shRB26s1RFlangeCr->DefineSection(2, z0, kRB26s1RFlangeCrRi2, kRB26s1RFlangeCrRo);
1622 z0 = kRB26s1RFlangeCrL;
1623 shRB26s1RFlangeCr->DefineSection(3, z0, kRB26s1RFlangeCrRi2, kRB26s1RFlangeCrRo);
1624 TGeoVolume* voRB26s1RFlangeCr =
1625 new TGeoVolume("RB26s1RFlangeCr", shRB26s1RFlangeCr, kMedSteel);
1626
1627 // Pos 2 Insert LHCVFX__0015
1628 const Float_t kRB26s1RFlangeIsL = 4.88 ; // Lenth of the insert
af5f3976 1629 const Float_t kRB26s1RFlangeIsR = 6.70/2. ; // Ring radius
820b4d9e 1630 const Float_t kRB26s1RFlangeIsD = 0.80 ; // Ring Width
1631
1632 TGeoPcon* shRB26s1RFlangeIs = new TGeoPcon(0., 360., 4);
1633 z0 = 0.;
1634 shRB26s1RFlangeIs->DefineSection(0, z0, kRB26s1RFlangeTubeRi, kRB26s1RFlangeIsR);
1635 z0 += kRB26s1RFlangeIsD;
1636 shRB26s1RFlangeIs->DefineSection(1, z0, kRB26s1RFlangeTubeRi, kRB26s1RFlangeIsR);
1637 shRB26s1RFlangeIs->DefineSection(2, z0, kRB26s1RFlangeTubeRi, kRB26s1RFlangeTubeRo);
1638 z0 = kRB26s1RFlangeIsL;
1639 shRB26s1RFlangeIs->DefineSection(3, z0, kRB26s1RFlangeTubeRi, kRB26s1RFlangeTubeRo);
1640 TGeoVolume* voRB26s1RFlangeIs =
1641 new TGeoVolume("RB26s1RFlangeIs", shRB26s1RFlangeIs, kMedSteel);
1642 // 4.88 + 3.7 = 8.58 (8.7 to avoid overlap)
1643 // Pos 3 Fixed Point Section LHCVC2A_0021
1644 const Float_t kRB26s1RFlangeFpL = 5.88 ; // Length of the fixed point section (0.08 cm added for welding)
1645 const Float_t kRB26s1RFlangeFpZ = 3.82 ; // Position of the ring
1646 const Float_t kRB26s1RFlangeFpD = 0.59 ; // Width of the ring
af5f3976 1647 const Float_t kRB26s1RFlangeFpR = 7.00/2. ; // Radius of the ring
820b4d9e 1648
1649 TGeoPcon* shRB26s1RFlangeFp = new TGeoPcon(0., 360., 6);
1650 z0 = 0.;
1651 shRB26s1RFlangeFp->DefineSection(0, z0, kRB26s1RFlangeTubeRi, kRB26s1RFlangeTubeRo);
1652 z0 += kRB26s1RFlangeFpZ;
1653 shRB26s1RFlangeFp->DefineSection(1, z0, kRB26s1RFlangeTubeRi, kRB26s1RFlangeTubeRo);
1654 shRB26s1RFlangeFp->DefineSection(2, z0, kRB26s1RFlangeTubeRi, kRB26s1RFlangeFpR);
1655 z0 += kRB26s1RFlangeFpD;
1656 shRB26s1RFlangeFp->DefineSection(3, z0, kRB26s1RFlangeTubeRi, kRB26s1RFlangeFpR);
1657 shRB26s1RFlangeFp->DefineSection(4, z0, kRB26s1RFlangeTubeRi, kRB26s1RFlangeTubeRo);
1658 z0 = kRB26s1RFlangeFpL;
1659 shRB26s1RFlangeFp->DefineSection(5, z0, kRB26s1RFlangeTubeRi, kRB26s1RFlangeTubeRo);
1660 TGeoVolume* voRB26s1RFlangeFp = new TGeoVolume("RB26s1RFlangeFp", shRB26s1RFlangeFp, kMedSteel);
1661
1662 // Put everything in a mother volume
1663 TGeoPcon* shRB26s1RFlange = new TGeoPcon(0., 360., 8);
1664 z0 = 0.;
1665 shRB26s1RFlange->DefineSection(0, z0, 0., kRB26s1RFlangeCrRo);
1666 z0 += kRB26s1RFlangeCrL;
1667 shRB26s1RFlange->DefineSection(1, z0, 0., kRB26s1RFlangeCrRo);
1668 shRB26s1RFlange->DefineSection(2, z0, 0., kRB26s1RFlangeTubeRo);
1669 z0 = kRB26s1RFlangeIsL + kRB26s1RFlangeFpZ;
1670 shRB26s1RFlange->DefineSection(3, z0, 0., kRB26s1RFlangeTubeRo);
1671 shRB26s1RFlange->DefineSection(4, z0, 0., kRB26s1RFlangeFpR);
1672 z0 += kRB26s1RFlangeFpD;
1673 shRB26s1RFlange->DefineSection(5, z0, 0., kRB26s1RFlangeFpR);
1674 shRB26s1RFlange->DefineSection(6, z0, 0., kRB26s1RFlangeTubeRo);
1675 z0 = kRB26s1RFlangeIsL + kRB26s1RFlangeFpL;
1676 shRB26s1RFlange->DefineSection(7, z0, 0., kRB26s1RFlangeTubeRo);
1677 TGeoVolume* voRB26s1RFlange = new TGeoVolume("RB26s1RFlange", shRB26s1RFlange, kMedVac);
1678
1679 voRB26s1RFlange->AddNode(voRB26s1RFlangeIs, 1, gGeoIdentity);
1680 voRB26s1RFlange->AddNode(voRB26s1RFlangeCr, 1, gGeoIdentity);
1681 voRB26s1RFlange->AddNode(voRB26s1RFlangeFp, 1, new TGeoTranslation(0., 0., kRB26s1RFlangeIsL));
1682
1683 ///////////////////////////////////
1684 // Fixed Flange //
1685 // Drawing LHCVFX_0006 //
1686 ///////////////////////////////////
1687 const Float_t kRB26s2FFlangeL = 2.13; // Length of the flange
1688 const Float_t kRB26s2FFlangeD1 = 0.97; // Length of section 1
1689 const Float_t kRB26s2FFlangeD2 = 0.29; // Length of section 2
1690 const Float_t kRB26s2FFlangeD3 = 0.87; // Length of section 3
1691 const Float_t kRB26s2FFlangeRo = 17.15/2.; // Flange outer radius
1692 const Float_t kRB26s2FFlangeRi1 = 12.30/2.; // Flange inner radius section 1
1693 const Float_t kRB26s2FFlangeRi2 = 12.00/2.; // Flange inner radius section 2
1694 const Float_t kRB26s2FFlangeRi3 = 12.30/2.; // Flange inner radius section 3
1695 z0 = 0;
1696 TGeoPcon* shRB26s2FFlange = new TGeoPcon(0., 360., 6);
1697 z0 = 0.;
1698 shRB26s2FFlange->DefineSection(0, z0, kRB26s2FFlangeRi1, kRB26s2FFlangeRo);
1699 z0 += kRB26s2FFlangeD1;
1700 shRB26s2FFlange->DefineSection(1, z0, kRB26s2FFlangeRi1, kRB26s2FFlangeRo);
1701 shRB26s2FFlange->DefineSection(2, z0, kRB26s2FFlangeRi2, kRB26s2FFlangeRo);
1702 z0 += kRB26s2FFlangeD2;
1703 shRB26s2FFlange->DefineSection(3, z0, kRB26s2FFlangeRi2, kRB26s2FFlangeRo);
1704 shRB26s2FFlange->DefineSection(4, z0, kRB26s2FFlangeRi3, kRB26s2FFlangeRo);
1705 z0 += kRB26s2FFlangeD3;
1706 shRB26s2FFlange->DefineSection(5, z0, kRB26s2FFlangeRi3, kRB26s2FFlangeRo);
1707 TGeoVolume* voRB26s2FFlange = new TGeoVolume("RB26s2FFlange", shRB26s2FFlange, kMedSteel);
1708
1709 TGeoVolume* voRB26s2FFlangeM = new TGeoVolume("RB26s2FFlangeM", MakeMotherFromTemplate(shRB26s2FFlange, 2, 5), kMedVac);
1710 voRB26s2FFlangeM->AddNode(voRB26s2FFlange, 1, gGeoIdentity);
1711
1712
1713
1714 ////////////////////////////////////////
1715 // //
1716 // RB26/3 //
1717 // Drawing LHCV2a_0048 //
1718 // Drawing LHCV2a_0002 //
1719 ////////////////////////////////////////
1720 //
1721 // Pos 1 Vacuum Tubes LHCVC2A__0003
1722 // Pos 2 Fixed Point LHCVFX___0005
1723 // Pos 3 Split Flange LHCVFX___0007
1724 // Pos 4 Fixed Flange LHCVFX___0004
1725 // Pos 5 Axial Compensator LHCVC2A__0065
1726 //
1727 //
1728 //
1729 //
1730 ///////////////////////////////////
1731 // Vacuum Tube //
1732 // Drawing LHCVC2A_0003 //
1733 ///////////////////////////////////
1734 const Float_t kRB26s3TubeL = 629.35 + 0.3; // 0.3 cm added for welding
1735 const Float_t kRB26s3TubeR1 = 12./2.;
1736 const Float_t kRB26s3TubeR2 = kRB26s3TubeR1 + 215.8 * TMath::Tan(0.829 / 180. * TMath::Pi());
1737
1738
1739 TGeoPcon* shRB26s3Tube = new TGeoPcon(0., 360., 7);
1740 // Section 1: straight section
1741 shRB26s3Tube->DefineSection(0, 0.00, kRB26s3TubeR1, kRB26s3TubeR1 + 0.15);
1742 shRB26s3Tube->DefineSection(1, 2.00, kRB26s3TubeR1, kRB26s3TubeR1 + 0.15);
1743 // Section 2: 0.829 deg opening cone
1744 shRB26s3Tube->DefineSection(2, 2.00, kRB26s3TubeR1, kRB26s3TubeR1 + 0.20);
1745
1746 shRB26s3Tube->DefineSection(3, 217.80, kRB26s3TubeR2, kRB26s3TubeR2 + 0.20);
1747 shRB26s3Tube->DefineSection(4, 217.80, kRB26s3TubeR2, kRB26s3TubeR2 + 0.30);
1748
1749 shRB26s3Tube->DefineSection(5, 622.20, 30.00/2., 30.60/2.);
1750 shRB26s3Tube->DefineSection(6, kRB26s3TubeL, 30.00/2., 30.60/2.);
1751
1752 TGeoVolume* voRB26s3Tube = new TGeoVolume("RB26s3Tube", shRB26s3Tube, kMedSteel);
1753 // Add the insulation layer
1754 TGeoVolume* voRB26s3TubeIns = new TGeoVolume("RB26s3TubeIns", MakeInsulationFromTemplate(shRB26s3Tube), kMedInsu);
1755 voRB26s3Tube->AddNode(voRB26s3TubeIns, 1, gGeoIdentity);
1756
1757 TGeoVolume* voRB26s3TubeM = new TGeoVolume("RB26s3TubeM", MakeMotherFromTemplate(shRB26s3Tube), kMedVac);
1758 voRB26s3TubeM->AddNode(voRB26s3Tube, 1, gGeoIdentity);
1759
1760
1761
1762 ///////////////////////////////////
1763 // Fixed Point //
1764 // Drawing LHCVFX_0005 //
1765 ///////////////////////////////////
1766 const Float_t kRB26s3FixedPointL = 16.37 ; // Length of the fixed point section (0.3 cm added for welding)
1767 const Float_t kRB26s3FixedPointZ = 9.72 ; // Position of the ring (0.15 cm added for welding)
1768 const Float_t kRB26s3FixedPointD = 0.595 ; // Width of the ring
1769 const Float_t kRB26s3FixedPointR = 13.30/2. ; // Radius of the ring
1770 const Float_t kRB26s3FixedPointRi = 12.00/2. ; // Inner radius of the tube
1771 const Float_t kRB26s3FixedPointRo1 = 12.30/2. ; // Outer radius of the tube (in)
1772 const Float_t kRB26s3FixedPointRo2 = 12.40/2. ; // Outer radius of the tube (out)
1773 const Float_t kRB26s3FixedPointDs = 1.5 ; // Width of straight section behind ring
1774 const Float_t kRB26s3FixedPointDc = 3.15 ; // Width of conical section behind ring (0.15 cm added for welding)
1775
1776 TGeoPcon* shRB26s3FixedPoint = new TGeoPcon(0., 360., 8);
1777 z0 = 0.;
1778 shRB26s3FixedPoint->DefineSection(0, z0, kRB26s3FixedPointRi, kRB26s3FixedPointRo1);
1779 z0 += kRB26s3FixedPointZ;
1780 shRB26s3FixedPoint->DefineSection(1, z0, kRB26s3FixedPointRi, kRB26s3FixedPointRo1);
1781 shRB26s3FixedPoint->DefineSection(2, z0, kRB26s3FixedPointRi, kRB26s3FixedPointR);
1782 z0 += kRB26s3FixedPointD;
1783 shRB26s3FixedPoint->DefineSection(3, z0, kRB26s3FixedPointRi, kRB26s3FixedPointR);
1784 shRB26s3FixedPoint->DefineSection(4, z0, kRB26s3FixedPointRi, kRB26s3FixedPointRo1);
1785 z0 += kRB26s3FixedPointDs;
1786 shRB26s3FixedPoint->DefineSection(5, z0, kRB26s3FixedPointRi, kRB26s3FixedPointRo1);
1787 z0 += kRB26s3FixedPointDc;
1788 shRB26s3FixedPoint->DefineSection(6, z0, kRB26s3FixedPointRi, kRB26s3FixedPointRo2);
1789 z0 = kRB26s3FixedPointL;
1790 shRB26s3FixedPoint->DefineSection(7, z0, kRB26s3FixedPointRi, kRB26s3FixedPointRo2);
1791 TGeoVolume* voRB26s3FixedPoint = new TGeoVolume("RB26s3FixedPoint", shRB26s3FixedPoint, kMedSteel);
1792
1793 TGeoVolume* voRB26s3FixedPointM = new TGeoVolume("RB26s3FixedPointM", MakeMotherFromTemplate(shRB26s3FixedPoint), kMedVac);
1794 voRB26s3FixedPointM->AddNode(voRB26s3FixedPoint, 1, gGeoIdentity);
1795
1796 ///////////////////////////////////
1797 // Split Flange //
1798 // Drawing LHCVFX_0005 //
1799 ///////////////////////////////////
1800 const Float_t kRB26s3SFlangeL = 2.13; // Length of the flange
1801 const Float_t kRB26s3SFlangeD1 = 0.57; // Length of section 1
1802 const Float_t kRB26s3SFlangeD2 = 0.36; // Length of section 2
1803 const Float_t kRB26s3SFlangeD3 = 0.50 + 0.70; // Length of section 3
1804 const Float_t kRB26s3SFlangeRo = 17.15/2.; // Flange outer radius
1805 const Float_t kRB26s3SFlangeRi1 = 12.30/2.; // Flange inner radius section 1
1806 const Float_t kRB26s3SFlangeRi2 = 12.00/2.; // Flange inner radius section 2
1807 const Float_t kRB26s3SFlangeRi3 = 12.30/2.; // Flange inner radius section 3
1808 z0 = 0;
1809 TGeoPcon* shRB26s3SFlange = new TGeoPcon(0., 360., 6);
1810 z0 = 0.;
1811 shRB26s3SFlange->DefineSection(0, z0, kRB26s3SFlangeRi1, kRB26s3SFlangeRo);
1812 z0 += kRB26s3SFlangeD1;
1813 shRB26s3SFlange->DefineSection(1, z0, kRB26s3SFlangeRi1, kRB26s3SFlangeRo);
1814 shRB26s3SFlange->DefineSection(2, z0, kRB26s3SFlangeRi2, kRB26s3SFlangeRo);
1815 z0 += kRB26s3SFlangeD2;
1816 shRB26s3SFlange->DefineSection(3, z0, kRB26s3SFlangeRi2, kRB26s3SFlangeRo);
1817 shRB26s3SFlange->DefineSection(4, z0, kRB26s3SFlangeRi3, kRB26s3SFlangeRo);
1818 z0 += kRB26s3SFlangeD3;
1819 shRB26s3SFlange->DefineSection(5, z0, kRB26s3SFlangeRi3, kRB26s3SFlangeRo);
1820 TGeoVolume* voRB26s3SFlange = new TGeoVolume("RB26s3SFlange", shRB26s3SFlange, kMedSteel);
1821
1822 TGeoVolume* voRB26s3SFlangeM = new TGeoVolume("RB26s3SFlangeM", MakeMotherFromTemplate(shRB26s3SFlange, 0, 3), kMedVac);
1823 voRB26s3SFlangeM->AddNode(voRB26s3SFlange, 1, gGeoIdentity);
1824
1825 ///////////////////////////////////
1826 // RB26/3 Fixed Flange //
1827 // Drawing LHCVFX___0004 //
1828 ///////////////////////////////////
1829 const Float_t kRB26s3FFlangeL = 2.99; // Length of the flange
1830 const Float_t kRB26s3FFlangeD1 = 1.72; // Length of section 1
1831 const Float_t kRB26s3FFlangeD2 = 0.30; // Length of section 2
1832 const Float_t kRB26s3FFlangeD3 = 0.97; // Length of section 3
1833 const Float_t kRB26s3FFlangeRo = 36.20/2.; // Flange outer radius
1834 const Float_t kRB26s3FFlangeRi1 = 30.60/2.; // Flange inner radius section 1
1835 const Float_t kRB26s3FFlangeRi2 = 30.00/2.; // Flange inner radius section 2
1836 const Float_t kRB26s3FFlangeRi3 = 30.60/2.; // Flange inner radius section 3
1837 z0 = 0;
1838 TGeoPcon* shRB26s3FFlange = new TGeoPcon(0., 360., 6);
1839 z0 = 0.;
1840 shRB26s3FFlange->DefineSection(0, z0, kRB26s3FFlangeRi1, kRB26s3FFlangeRo);
1841 z0 += kRB26s3FFlangeD1;
1842 shRB26s3FFlange->DefineSection(1, z0, kRB26s3FFlangeRi1, kRB26s3FFlangeRo);
1843 shRB26s3FFlange->DefineSection(2, z0, kRB26s3FFlangeRi2, kRB26s3FFlangeRo);
1844 z0 += kRB26s3FFlangeD2;
1845 shRB26s3FFlange->DefineSection(3, z0, kRB26s3FFlangeRi2, kRB26s3FFlangeRo);
1846 shRB26s3FFlange->DefineSection(4, z0, kRB26s3FFlangeRi3, kRB26s3FFlangeRo);
1847 z0 += kRB26s3FFlangeD3;
1848 shRB26s3FFlange->DefineSection(5, z0, kRB26s3FFlangeRi3, kRB26s3FFlangeRo);
1849 TGeoVolume* voRB26s3FFlange = new TGeoVolume("RB26s3FFlange", shRB26s3FFlange, kMedSteel);
1850
1851 TGeoVolume* voRB26s3FFlangeM = new TGeoVolume("RB26s3FFlangeM", MakeMotherFromTemplate(shRB26s3FFlange, 2, 5), kMedVac);
1852 voRB26s3FFlangeM->AddNode(voRB26s3FFlange, 1, gGeoIdentity);
1853
1854
1855
1856 ///////////////////////////////////
1857 // RB26/3 Axial Compensator //
1858 // Drawing LHCVC2a_0065 //
1859 ///////////////////////////////////
1860 const Float_t kRB26s3CompL = 42.0; // Length of the compensator (0.3 cm added for welding)
1861 const Float_t kRB26s3BellowRo = 34.00/2.; // Bellow outer radius [Pos 1]
1862 const Float_t kRB26s3BellowRi = 30.10/2.; // Bellow inner radius [Pos 1]
1863 const Int_t kRB26s3NumberOfPlies = 13; // Number of plies [Pos 1]
1864 const Float_t kRB26s3BellowUndL = 17.70; // Length of undulated region [Pos 1]
1865 const Float_t kRB26s3PlieThickness = 0.06; // Plie thickness [Pos 1]
1866 const Float_t kRB26s3ConnectionPlieR = 0.21; // Connection plie radius [Pos 1]
1867 // Plie radius
1868 const Float_t kRB26s3PlieR =
1869 (kRB26s3BellowUndL - 4. * kRB26s3ConnectionPlieR + 2. * kRB26s3PlieThickness +
1870 (2. * kRB26s3NumberOfPlies - 2.) * kRB26s3PlieThickness) / (4. * kRB26s3NumberOfPlies - 2.);
1871
1872 //
1873 // The welding tubes have 3 sections with different radii and 2 transition regions.
1874 // Section 1: connection to the outside
1875 // Section 2: commection to the bellow
1876 // Section 3: between 1 and 2
1877 const Float_t kRB26s3CompTubeInnerR1 = 30.0/2.; // Outer Connection tubes inner radius [Pos 4 + 3]
1878 const Float_t kRB26s3CompTubeOuterR1 = 30.6/2.; // Outer Connection tubes outer radius [Pos 4 + 3]
1879 const Float_t kRB26s3CompTubeInnerR2 = 29.4/2.; // Connection tubes inner radius [Pos 4 + 3]
1880 const Float_t kRB26s3CompTubeOuterR2 = 30.0/2.; // Connection tubes outer radius [Pos 4 + 3]
1881 const Float_t kRB26s3CompTubeInnerR3 = 30.6/2.; // Connection tubes inner radius at bellow [Pos 4 + 3]
1882 const Float_t kRB26s3CompTubeOuterR3 = 32.2/2.; // Connection tubes outer radius at bellow [Pos 4 + 3]
1883
1884 const Float_t kRB26s3WeldingTubeLeftL1 = 2.0; // Left connection tube length [Pos 4]
1885 const Float_t kRB26s3WeldingTubeLeftL2 = 3.4; // Left connection tube length [Pos 4]
1886 const Float_t kRB26s3WeldingTubeLeftL = 7.0; // Left connection tube total length [Pos 4]
1887 const Float_t kRB26s3WeldingTubeRightL1 = 2.3; // Right connection tube length [Pos 3] (0.3 cm added for welding)
1888 const Float_t kRB26s3WeldingTubeRightL2 = 13.4; // Right connection tube length [Pos 3]
1889
1890 const Float_t kRB26s3WeldingTubeT1 = 0.6; // Length of first r-transition [Pos 4 + 3]
1891 const Float_t kRB26s3WeldingTubeT2 = 1.0; // Length of 2nd r-transition [Pos 4 + 3]
1892
1893
1894
1895 const Float_t kRB26s3RingOuterR = 36.1/2.; // Ring inner radius [Pos 4]
1896 const Float_t kRB26s3RingL = 0.8/2.; // Ring half length [Pos 4]
1897 const Float_t kRB26s3RingZ = 3.7 ; // Ring z-position [Pos 4]
1898 const Float_t kRB26s3ProtOuterR = 36.2/2.; // Protection tube outer radius [Pos 2]
1899 const Float_t kRB26s3ProtL = 27.0/2.; // Protection tube half length [Pos 2]
1900 const Float_t kRB26s3ProtZ = 4.0 ; // Protection tube z-position [Pos 2]
1901
1902
1903 // Mother volume
1904 //
1905 TGeoPcon* shRB26s3Compensator = new TGeoPcon(0., 360., 6);
1906 shRB26s3Compensator->DefineSection( 0, 0.0, 0., kRB26s3CompTubeOuterR1);
1907 shRB26s3Compensator->DefineSection( 1, kRB26s3RingZ, 0., kRB26s3CompTubeOuterR1);
1908 shRB26s3Compensator->DefineSection( 2, kRB26s3RingZ, 0., kRB26s3ProtOuterR);
1909 shRB26s3Compensator->DefineSection( 3, kRB26s3ProtZ + 2. * kRB26s3ProtL, 0., kRB26s3ProtOuterR);
1910 shRB26s3Compensator->DefineSection( 4, kRB26s3ProtZ + 2. * kRB26s3ProtL, 0., kRB26s3CompTubeOuterR1);
1911 shRB26s3Compensator->DefineSection( 5, kRB26s3CompL , 0., kRB26s3CompTubeOuterR1);
1912 TGeoVolume* voRB26s3Compensator =
1913 new TGeoVolume("RB26s3Compensator", shRB26s3Compensator, kMedVac);
1914
1915 //
1916 // [Pos 1] Bellow
1917 //
1918 //
1919 TGeoVolume* voRB26s3Bellow = new TGeoVolume("RB26s3Bellow",
1920 new TGeoTube(kRB26s3BellowRi, kRB26s3BellowRo, kRB26s3BellowUndL/2.), kMedVac);
1921 //
1922 // Upper part of the undulation
1923 //
1924 TGeoTorus* shRB26s3PlieTorusU = new TGeoTorus(kRB26s3BellowRo - kRB26s3PlieR, kRB26s3PlieR - kRB26s3PlieThickness, kRB26s3PlieR);
1925 shRB26s3PlieTorusU->SetName("RB26s3TorusU");
1926 TGeoTube* shRB26s3PlieTubeU = new TGeoTube (kRB26s3BellowRo - kRB26s3PlieR, kRB26s3BellowRo, kRB26s3PlieR);
1927 shRB26s3PlieTubeU->SetName("RB26s3TubeU");
1928 TGeoCompositeShape* shRB26s3UpperPlie = new TGeoCompositeShape("RB26s3UpperPlie", "RB26s3TorusU*RB26s3TubeU");
1929
1930 TGeoVolume* voRB26s3WiggleU = new TGeoVolume("RB26s3UpperPlie", shRB26s3UpperPlie, kMedSteel);
1931 //
1932 // Lower part of the undulation
1933 TGeoTorus* shRB26s3PlieTorusL = new TGeoTorus(kRB26s3BellowRi + kRB26s3PlieR, kRB26s3PlieR - kRB26s3PlieThickness, kRB26s3PlieR);
1934 shRB26s3PlieTorusL->SetName("RB26s3TorusL");
1935 TGeoTube* shRB26s3PlieTubeL = new TGeoTube (kRB26s3BellowRi, kRB26s3BellowRi + kRB26s3PlieR, kRB26s3PlieR);
1936 shRB26s3PlieTubeL->SetName("RB26s3TubeL");
1937 TGeoCompositeShape* shRB26s3LowerPlie = new TGeoCompositeShape("RB26s3LowerPlie", "RB26s3TorusL*RB26s3TubeL");
1938
1939 TGeoVolume* voRB26s3WiggleL = new TGeoVolume("RB26s3LowerPlie", shRB26s3LowerPlie, kMedSteel);
1940
1941 //
1942 // Connection between upper and lower part of undulation
1943 TGeoVolume* voRB26s3WiggleC1 = new TGeoVolume("RB26s3PlieConn1",
1944 new TGeoTube(kRB26s3BellowRi + kRB26s3PlieR,
1945 kRB26s3BellowRo - kRB26s3PlieR, kRB26s3PlieThickness / 2.), kMedSteel);
1946 //
1947 // One wiggle
1948 TGeoVolumeAssembly* voRB26s3Wiggle = new TGeoVolumeAssembly("RB26s3Wiggle");
1949 z0 = - kRB26s3PlieThickness / 2.;
1950 voRB26s3Wiggle->AddNode(voRB26s3WiggleC1, 1 , new TGeoTranslation(0., 0., z0));
1951 z0 += kRB26s3PlieR - kRB26s3PlieThickness / 2.;
1952 voRB26s3Wiggle->AddNode(voRB26s3WiggleU, 1 , new TGeoTranslation(0., 0., z0));
1953 z0 += kRB26s3PlieR - kRB26s3PlieThickness / 2.;
1954 voRB26s3Wiggle->AddNode(voRB26s3WiggleC1, 2 , new TGeoTranslation(0., 0., z0));
1955 z0 += kRB26s3PlieR - kRB26s3PlieThickness;
1956 voRB26s3Wiggle->AddNode(voRB26s3WiggleL, 1 , new TGeoTranslation(0., 0., z0));
1957 // Positioning of the volumes
1958 z0 = - kRB26s3BellowUndL/2.+ kRB26s3ConnectionPlieR;
1959 voRB26s3Bellow->AddNode(voRB26s3WiggleL, 1, new TGeoTranslation(0., 0., z0));
1960 z0 += kRB26s3ConnectionPlieR;
1961 zsh = 4. * kRB26s3PlieR - 2. * kRB26s3PlieThickness;
1962 for (Int_t iw = 0; iw < kRB26s3NumberOfPlies; iw++) {
1963 Float_t zpos = z0 + iw * zsh;
1964 voRB26s3Bellow->AddNode(voRB26s3Wiggle, iw + 1, new TGeoTranslation(0., 0., zpos - kRB26s3PlieThickness));
1965 }
1966
1967 voRB26s3Compensator->AddNode(voRB26s3Bellow, 1, new TGeoTranslation(0., 0., kRB26s3WeldingTubeLeftL + kRB26s3BellowUndL/2.));
1968
1969
1970 //
1971 // [Pos 2] Outer Protecting Tube
1972 //
1973 TGeoTube* shRB26s3CompProtTube = new TGeoTube(kRB26s3RingOuterR, kRB26s3ProtOuterR, kRB26s3ProtL);
1974 TGeoVolume* voRB26s3CompProtTube =
1975 new TGeoVolume("RB26s3CompProtTube", shRB26s3CompProtTube, kMedSteel);
1976 voRB26s3Compensator->AddNode(voRB26s3CompProtTube, 1, new TGeoTranslation(0., 0., kRB26s3ProtZ + kRB26s3ProtL));
1977
1978
1979 //
1980 // [Pos 3] Right Welding Tube
1981 //
1982 TGeoPcon* shRB26s3CompRightTube = new TGeoPcon(0., 360., 5);
1983 z0 = 0.;
1984 shRB26s3CompRightTube->DefineSection(0, z0, kRB26s3CompTubeInnerR3, kRB26s3CompTubeOuterR3);
1985 z0 += kRB26s3WeldingTubeT2;
1986 shRB26s3CompRightTube->DefineSection(1, z0, kRB26s3CompTubeInnerR2, kRB26s3CompTubeOuterR2);
1987 z0 += kRB26s3WeldingTubeRightL2;
1988 shRB26s3CompRightTube->DefineSection(2, z0, kRB26s3CompTubeInnerR2, kRB26s3CompTubeOuterR2);
1989 z0 += kRB26s3WeldingTubeT1;
1990 shRB26s3CompRightTube->DefineSection(3, z0, kRB26s3CompTubeInnerR1, kRB26s3CompTubeOuterR1);
1991 z0 += kRB26s3WeldingTubeRightL1;
1992 shRB26s3CompRightTube->DefineSection(4, z0, kRB26s3CompTubeInnerR1, kRB26s3CompTubeOuterR1);
1993
1994 TGeoVolume* voRB26s3CompRightTube =
1995 new TGeoVolume("RB26s3CompRightTube", shRB26s3CompRightTube, kMedSteel);
1996 voRB26s3Compensator->AddNode(voRB26s3CompRightTube, 1, new TGeoTranslation(0., 0., kRB26s3CompL - z0));
1997
1998 //
1999 // [Pos 4] Left Welding Tube
2000 //
2001 TGeoPcon* shRB26s3CompLeftTube = new TGeoPcon(0., 360., 5);
2002 z0 = 0.;
2003 shRB26s3CompLeftTube->DefineSection(0, z0, kRB26s3CompTubeInnerR1, kRB26s3CompTubeOuterR1);
2004 z0 += kRB26s3WeldingTubeLeftL1;
2005 shRB26s3CompLeftTube->DefineSection(1, z0, kRB26s3CompTubeInnerR1, kRB26s3CompTubeOuterR1);
2006 z0 += kRB26s3WeldingTubeT1;
2007 shRB26s3CompLeftTube->DefineSection(2, z0, kRB26s3CompTubeInnerR2, kRB26s3CompTubeOuterR2);
2008 z0 += kRB26s3WeldingTubeLeftL2;
2009 shRB26s3CompLeftTube->DefineSection(3, z0, kRB26s3CompTubeInnerR2, kRB26s3CompTubeOuterR2);
2010 z0 += kRB26s3WeldingTubeT2;
2011 shRB26s3CompLeftTube->DefineSection(4, z0, kRB26s3CompTubeInnerR3, kRB26s3CompTubeOuterR3);
2012
2013 TGeoVolume* voRB26s3CompLeftTube =
2014 new TGeoVolume("RB26s3CompLeftTube", shRB26s3CompLeftTube, kMedSteel);
2015 voRB26s3Compensator->AddNode(voRB26s3CompLeftTube, 1, gGeoIdentity);
2016 //
2017 // [Pos 5] Ring
2018 //
2019 TGeoTube* shRB26s3CompRing = new TGeoTube(kRB26s3CompTubeOuterR2, kRB26s3RingOuterR, kRB26s3RingL);
2020 TGeoVolume* voRB26s3CompRing =
2021 new TGeoVolume("RB26s3CompRing", shRB26s3CompRing, kMedSteel);
2022 voRB26s3Compensator->AddNode(voRB26s3CompRing, 1, new TGeoTranslation(0., 0., kRB26s3RingZ + kRB26s3RingL));
2023
2024
2025
2026 ///////////////////////////////////////////
2027 // //
2028 // RB26/4-5 //
2029 // Drawing LHCV2a_0012 [as installed] //
2030 ////////////////////////////////////////////
2031 // Pos1 Vacuum Tubes LHCVC2A__0014
2032 // Pos2 Compensator LHCVC2A__0066
2033 // Pos3 Fixed Point Section LHCVC2A__0016
2034 // Pos4 Split Flange LHCVFX___0005
2035 // Pos5 RotableFlange LHCVFX___0009
2036 ////////////////////////////////////////////
2037
2038 ///////////////////////////////////
2039 // RB26/4-5 Vacuum Tubes //
2040 // Drawing LHCVC2a_0014 //
2041 ///////////////////////////////////
2042 const Float_t kRB26s45TubeL = 593.12 + 0.3; // 0.3 cm added for welding
2043
2044 TGeoPcon* shRB26s45Tube = new TGeoPcon(0., 360., 11);
2045 // Section 1: straight section
2046 shRB26s45Tube->DefineSection( 0, 0.00, 30.00/2., 30.60/2.);
2047 shRB26s45Tube->DefineSection( 1, 1.20, 30.00/2., 30.60/2.);
2048 shRB26s45Tube->DefineSection( 2, 1.20, 30.00/2., 30.80/2.);
2049 shRB26s45Tube->DefineSection( 3, 25.10, 30.00/2., 30.80/2.);
2050 // Section 2: 0.932 deg opening cone
2051 shRB26s45Tube->DefineSection( 4, 486.10, 45.00/2., 45.80/2.);
2052 // Section 3: straight section 4 mm
2053 shRB26s45Tube->DefineSection( 5, 512.10, 45.00/2., 45.80/2.);
2054 // Section 4: straight section 3 mm
2055 shRB26s45Tube->DefineSection( 6, 512.10, 45.00/2., 45.60/2.);
2056 shRB26s45Tube->DefineSection( 7, 527.70, 45.00/2., 45.60/2.);
2057 // Section 4: closing cone
2058 shRB26s45Tube->DefineSection( 8, 591.30, 10.00/2., 10.60/2.);
2059 shRB26s45Tube->DefineSection( 9, 591.89, 10.00/2., 10.30/2.);
2060
2061 shRB26s45Tube->DefineSection(10, kRB26s45TubeL, 10.00/2., 10.30/2.);
2062 TGeoVolume* voRB26s45Tube =
2063 new TGeoVolume("RB26s45Tube", shRB26s45Tube, kMedSteel);
2064
2065 TGeoVolume* voRB26s45TubeM = new TGeoVolume("RB26s45TubeM", MakeMotherFromTemplate(shRB26s45Tube), kMedVac);
2066 voRB26s45TubeM->AddNode(voRB26s45Tube, 1, gGeoIdentity);
2067
2068
2069
2070 ///////////////////////////////////
2071 // RB26/5 Axial Compensator //
2072 // Drawing LHCVC2a_0066 //
2073 ///////////////////////////////////
2074 const Float_t kRB26s5CompL = 27.60; // Length of the compensator (0.30 cm added for welding)
2075 const Float_t kRB26s5BellowRo = 12.48/2.; // Bellow outer radius [Pos 1]
2076 const Float_t kRB26s5BellowRi = 10.32/2.; // Bellow inner radius [Pos 1]
2077 const Int_t kRB26s5NumberOfPlies = 15; // Number of plies [Pos 1]
2078 const Float_t kRB26s5BellowUndL = 10.50; // Length of undulated region [Pos 1]
2079 const Float_t kRB26s5PlieThickness = 0.025; // Plie thickness [Pos 1]
2080 const Float_t kRB26s5ConnectionPlieR = 0.21; // Connection plie radius [Pos 1]
2081 const Float_t kRB26s5ConnectionR = 11.2/2.; // Bellow connection radius [Pos 1]
2082 // Plie radius
2083 const Float_t kRB26s5PlieR =
2084 (kRB26s5BellowUndL - 4. * kRB26s5ConnectionPlieR + 2. * kRB26s5PlieThickness +
2085 (2. * kRB26s5NumberOfPlies - 2.) * kRB26s5PlieThickness) / (4. * kRB26s5NumberOfPlies - 2.);
2086 const Float_t kRB26s5CompTubeInnerR = 10.00/2.; // Connection tubes inner radius [Pos 2 + 3]
2087 const Float_t kRB26s5CompTubeOuterR = 10.30/2.; // Connection tubes outer radius [Pos 2 + 3]
2088 const Float_t kRB26s5WeldingTubeLeftL = 3.70/2.; // Left connection tube half length [Pos 2]
2089 const Float_t kRB26s5WeldingTubeRightL = 13.40/2.; // Right connection tube half length [Pos 3] (0.3 cm added for welding)
2090 const Float_t kRB26s5RingInnerR = 11.2/2.; // Ring inner radius [Pos 4]
2091 const Float_t kRB26s5RingOuterR = 16.0/2.; // Ring inner radius [Pos 4]
2092 const Float_t kRB26s5RingL = 0.4/2.; // Ring half length [Pos 4]
2093 const Float_t kRB26s5RingZ = 14.97; // Ring z-position [Pos 4]
2094 const Float_t kRB26s5ProtOuterR = 16.2/2.; // Protection tube outer radius [Pos 5]
2095 const Float_t kRB26s5ProtL = 13.0/2.; // Protection tube half length [Pos 5]
2096 const Float_t kRB26s5ProtZ = 2.17; // Protection tube z-position [Pos 5]
2097 const Float_t kRB26s5DetailZR = 11.3/2.; // Detail Z max radius
2098
2099
2100 // Mother volume
2101 //
2102 TGeoPcon* shRB26s5Compensator = new TGeoPcon(0., 360., 8);
2103 shRB26s5Compensator->DefineSection( 0, 0.0, 0., kRB26s5CompTubeOuterR);
2104 shRB26s5Compensator->DefineSection( 1, kRB26s5ProtZ, 0., kRB26s5CompTubeOuterR);
2105 shRB26s5Compensator->DefineSection( 2, kRB26s5ProtZ, 0., kRB26s5ProtOuterR);
2106 shRB26s5Compensator->DefineSection( 3, kRB26s5ProtZ + 2. * kRB26s5ProtL + 2. * kRB26s5RingL, 0., kRB26s5ProtOuterR);
2107 shRB26s5Compensator->DefineSection( 4, kRB26s5ProtZ + 2. * kRB26s5ProtL + 2. * kRB26s5RingL, 0., kRB26s5DetailZR);
2108 shRB26s5Compensator->DefineSection( 5, kRB26s5CompL - 8., 0., kRB26s5DetailZR);
2109 shRB26s5Compensator->DefineSection( 6, kRB26s5CompL - 8., 0., kRB26s5CompTubeOuterR);
2110 shRB26s5Compensator->DefineSection( 7, kRB26s5CompL, 0., kRB26s5CompTubeOuterR);
2111 TGeoVolume* voRB26s5Compensator = new TGeoVolume("RB26s5Compensator", shRB26s5Compensator, kMedVac);
2112
2113 //
2114 // [Pos 1] Bellow
2115 //
2116 //
2117 TGeoVolume* voRB26s5Bellow = new TGeoVolume("RB26s5Bellow",
2118 new TGeoTube(kRB26s5BellowRi, kRB26s5BellowRo, kRB26s5BellowUndL/2.), kMedVac);
2119 //
2120 // Upper part of the undulation
2121 //
2122 TGeoTorus* shRB26s5PlieTorusU = new TGeoTorus(kRB26s5BellowRo - kRB26s5PlieR, kRB26s5PlieR - kRB26s5PlieThickness, kRB26s5PlieR);
2123 shRB26s5PlieTorusU->SetName("RB26s5TorusU");
2124 TGeoTube* shRB26s5PlieTubeU = new TGeoTube (kRB26s5BellowRo - kRB26s5PlieR, kRB26s5BellowRo, kRB26s5PlieR);
2125 shRB26s5PlieTubeU->SetName("RB26s5TubeU");
2126 TGeoCompositeShape* shRB26s5UpperPlie = new TGeoCompositeShape("RB26s5UpperPlie", "RB26s5TorusU*RB26s5TubeU");
2127
2128 TGeoVolume* voRB26s5WiggleU = new TGeoVolume("RB26s5UpperPlie", shRB26s5UpperPlie, kMedSteel);
2129 //
2130 // Lower part of the undulation
2131 TGeoTorus* shRB26s5PlieTorusL = new TGeoTorus(kRB26s5BellowRi + kRB26s5PlieR, kRB26s5PlieR - kRB26s5PlieThickness, kRB26s5PlieR);
2132 shRB26s5PlieTorusL->SetName("RB26s5TorusL");
2133 TGeoTube* shRB26s5PlieTubeL = new TGeoTube (kRB26s5BellowRi, kRB26s5BellowRi + kRB26s5PlieR, kRB26s5PlieR);
2134 shRB26s5PlieTubeL->SetName("RB26s5TubeL");
2135 TGeoCompositeShape* shRB26s5LowerPlie = new TGeoCompositeShape("RB26s5LowerPlie", "RB26s5TorusL*RB26s5TubeL");
2136
2137 TGeoVolume* voRB26s5WiggleL = new TGeoVolume("RB26s5LowerPlie", shRB26s5LowerPlie, kMedSteel);
2138
2139 //
2140 // Connection between upper and lower part of undulation
2141 TGeoVolume* voRB26s5WiggleC1 = new TGeoVolume("RB26s5PlieConn1",
2142 new TGeoTube(kRB26s5BellowRi + kRB26s5PlieR,
2143 kRB26s5BellowRo - kRB26s5PlieR, kRB26s5PlieThickness / 2.), kMedSteel);
2144 //
2145 // One wiggle
2146 TGeoVolumeAssembly* voRB26s5Wiggle = new TGeoVolumeAssembly("RB26s5Wiggle");
2147 z0 = - kRB26s5PlieThickness / 2.;
2148 voRB26s5Wiggle->AddNode(voRB26s5WiggleC1, 1 , new TGeoTranslation(0., 0., z0));
2149 z0 += kRB26s5PlieR - kRB26s5PlieThickness / 2.;
2150 voRB26s5Wiggle->AddNode(voRB26s5WiggleU, 1 , new TGeoTranslation(0., 0., z0));
2151 z0 += kRB26s5PlieR - kRB26s5PlieThickness / 2.;
2152 voRB26s5Wiggle->AddNode(voRB26s5WiggleC1, 2 , new TGeoTranslation(0., 0., z0));
2153 z0 += kRB26s5PlieR - kRB26s5PlieThickness;
2154 voRB26s5Wiggle->AddNode(voRB26s5WiggleL , 1 , new TGeoTranslation(0., 0., z0));
2155 // Positioning of the volumes
2156 z0 = - kRB26s5BellowUndL/2.+ kRB26s5ConnectionPlieR;
2157 voRB26s5Bellow->AddNode(voRB26s5WiggleL, 1, new TGeoTranslation(0., 0., z0));
2158 z0 += kRB26s5ConnectionPlieR;
2159 zsh = 4. * kRB26s5PlieR - 2. * kRB26s5PlieThickness;
2160 for (Int_t iw = 0; iw < kRB26s5NumberOfPlies; iw++) {
2161 Float_t zpos = z0 + iw * zsh;
2162 voRB26s5Bellow->AddNode(voRB26s5Wiggle, iw + 1, new TGeoTranslation(0., 0., zpos - kRB26s5PlieThickness));
2163 }
2164
2165 voRB26s5Compensator->AddNode(voRB26s5Bellow, 1, new TGeoTranslation(0., 0., 2. * kRB26s5WeldingTubeLeftL + kRB26s5BellowUndL/2.));
2166
2167 //
2168 // [Pos 2] Left Welding Tube
2169 //
2170 TGeoPcon* shRB26s5CompLeftTube = new TGeoPcon(0., 360., 3);
2171 z0 = 0;
2172 shRB26s5CompLeftTube->DefineSection(0, z0, kRB26s5CompTubeInnerR, kRB26s5CompTubeOuterR);
2173 z0 += 2 * kRB26s5WeldingTubeLeftL - ( kRB26s5ConnectionR - kRB26s5CompTubeOuterR);
2174 shRB26s5CompLeftTube->DefineSection(1, z0, kRB26s5CompTubeInnerR, kRB26s5CompTubeOuterR);
2175 z0 += ( kRB26s5ConnectionR - kRB26s5CompTubeOuterR);
2176 shRB26s5CompLeftTube->DefineSection(2, z0, kRB26s5ConnectionR - 0.15, kRB26s5ConnectionR);
2177 TGeoVolume* voRB26s5CompLeftTube = new TGeoVolume("RB26s5CompLeftTube", shRB26s5CompLeftTube, kMedSteel);
2178 voRB26s5Compensator->AddNode(voRB26s5CompLeftTube, 1, gGeoIdentity);
2179 //
2180 // [Pos 3] Right Welding Tube
2181 //
2182 TGeoPcon* shRB26s5CompRightTube = new TGeoPcon(0., 360., 11);
2183 // Detail Z
2184 shRB26s5CompRightTube->DefineSection( 0, 0. , kRB26s5CompTubeInnerR + 0.22, 11.2/2.);
2185 shRB26s5CompRightTube->DefineSection( 1, 0.05, kRB26s5CompTubeInnerR + 0.18, 11.2/2.);
2186 shRB26s5CompRightTube->DefineSection( 2, 0.22, kRB26s5CompTubeInnerR , 11.2/2. - 0.22);
2187 shRB26s5CompRightTube->DefineSection( 3, 0.44, kRB26s5CompTubeInnerR , 11.2/2.);
2188 shRB26s5CompRightTube->DefineSection( 4, 1.70, kRB26s5CompTubeInnerR , 11.2/2.);
2189 shRB26s5CompRightTube->DefineSection( 5, 2.10, kRB26s5CompTubeInnerR , kRB26s5CompTubeOuterR);
2190 shRB26s5CompRightTube->DefineSection( 6, 2.80, kRB26s5CompTubeInnerR , kRB26s5CompTubeOuterR);
2191 shRB26s5CompRightTube->DefineSection( 7, 2.80, kRB26s5CompTubeInnerR , 11.3/2.);
2192 shRB26s5CompRightTube->DefineSection( 8, 3.40, kRB26s5CompTubeInnerR , 11.3/2.);
2193 // Normal pipe
2194 shRB26s5CompRightTube->DefineSection( 9, 3.50, kRB26s5CompTubeInnerR , kRB26s5CompTubeOuterR);
2195 shRB26s5CompRightTube->DefineSection(10, 2. * kRB26s5WeldingTubeRightL, kRB26s5CompTubeInnerR, kRB26s5CompTubeOuterR);
2196
2197 TGeoVolume* voRB26s5CompRightTube =
2198 new TGeoVolume("RB26s5CompRightTube", shRB26s5CompRightTube, kMedSteel);
2199 voRB26s5Compensator->AddNode(voRB26s5CompRightTube, 1,
2200 new TGeoTranslation(0., 0., kRB26s5CompL - 2. * kRB26s5WeldingTubeRightL));
2201 //
2202 // [Pos 4] Ring
2203 //
2204 TGeoTube* shRB26s5CompRing = new TGeoTube(kRB26s5RingInnerR, kRB26s5RingOuterR, kRB26s5RingL);
2205 TGeoVolume* voRB26s5CompRing =
2206 new TGeoVolume("RB26s5CompRing", shRB26s5CompRing, kMedSteel);
2207 voRB26s5Compensator->AddNode(voRB26s5CompRing, 1, new TGeoTranslation(0., 0., kRB26s5RingZ + kRB26s5RingL));
2208
2209 //
2210 // [Pos 5] Outer Protecting Tube
2211 //
2212 TGeoTube* shRB26s5CompProtTube = new TGeoTube(kRB26s5RingOuterR, kRB26s5ProtOuterR, kRB26s5ProtL);
2213 TGeoVolume* voRB26s5CompProtTube =
2214 new TGeoVolume("RB26s5CompProtTube", shRB26s5CompProtTube, kMedSteel);
2215 voRB26s5Compensator->AddNode(voRB26s5CompProtTube, 1, new TGeoTranslation(0., 0., kRB26s5ProtZ + kRB26s5ProtL));
2216
2217 ///////////////////////////////////////
2218 // RB26/4 Fixed Point Section //
2219 // Drawing LHCVC2a_0016 //
2220 ///////////////////////////////////////
2221 const Float_t kRB26s4TubeRi = 30.30/2. ; // Tube inner radius (0.3 cm added for welding)
2222 const Float_t kRB26s4TubeRo = 30.60/2. ; // Tube outer radius
2223 const Float_t kRB26s4FixedPointL = 12.63 ; // Length of the fixed point section
2224 const Float_t kRB26s4FixedPointZ = 10.53 ; // Position of the ring (0.15 added for welding)
2225 const Float_t kRB26s4FixedPointD = 0.595 ; // Width of the ring
2226 const Float_t kRB26s4FixedPointR = 31.60/2. ; // Radius of the ring
2227
2228 TGeoPcon* shRB26s4FixedPoint = new TGeoPcon(0., 360., 6);
2229 z0 = 0.;
2230 shRB26s4FixedPoint->DefineSection(0, z0, kRB26s4TubeRi, kRB26s4TubeRo);
2231 z0 += kRB26s4FixedPointZ;
2232 shRB26s4FixedPoint->DefineSection(1, z0, kRB26s4TubeRi, kRB26s4TubeRo);
2233 shRB26s4FixedPoint->DefineSection(2, z0, kRB26s4TubeRi, kRB26s4FixedPointR);
2234 z0 += kRB26s4FixedPointD;
2235 shRB26s4FixedPoint->DefineSection(3, z0, kRB26s4TubeRi, kRB26s4FixedPointR);
2236 shRB26s4FixedPoint->DefineSection(4, z0, kRB26s4TubeRi, kRB26s4TubeRo);
2237 z0 = kRB26s4FixedPointL;
2238 shRB26s4FixedPoint->DefineSection(5, z0, kRB26s4TubeRi, kRB26s4TubeRo);
2239 TGeoVolume* voRB26s4FixedPoint = new TGeoVolume("RB26s4FixedPoint", shRB26s4FixedPoint, kMedSteel);
2240
2241 TGeoVolume* voRB26s4FixedPointM = new TGeoVolume("RB26s4FixedPointM", MakeMotherFromTemplate(shRB26s4FixedPoint), kMedVac);
2242 voRB26s4FixedPointM->AddNode(voRB26s4FixedPoint, 1, gGeoIdentity);
2243
2244
2245 ///////////////////////////////////////
2246 // RB26/4 Split Flange //
2247 // Drawing LHCVFX__0005 //
2248 ///////////////////////////////////////
2249 const Float_t kRB26s4SFlangeL = 2.99; // Length of the flange
2250 const Float_t kRB26s4SFlangeD1 = 0.85; // Length of section 1
2251 const Float_t kRB26s4SFlangeD2 = 0.36; // Length of section 2
2252 const Float_t kRB26s4SFlangeD3 = 0.73 + 1.05; // Length of section 3
2253 const Float_t kRB26s4SFlangeRo = 36.20/2.; // Flange outer radius
2254 const Float_t kRB26s4SFlangeRi1 = 30.60/2.; // Flange inner radius section 1
2255 const Float_t kRB26s4SFlangeRi2 = 30.00/2.; // Flange inner radius section 2
2256 const Float_t kRB26s4SFlangeRi3 = 30.60/2.; // Flange inner radius section 3
2257 z0 = 0;
2258 TGeoPcon* shRB26s4SFlange = new TGeoPcon(0., 360., 6);
2259 z0 = 0.;
2260 shRB26s4SFlange->DefineSection(0, z0, kRB26s4SFlangeRi1, kRB26s4SFlangeRo);
2261 z0 += kRB26s4SFlangeD1;
2262 shRB26s4SFlange->DefineSection(1, z0, kRB26s4SFlangeRi1, kRB26s4SFlangeRo);
2263 shRB26s4SFlange->DefineSection(2, z0, kRB26s4SFlangeRi2, kRB26s4SFlangeRo);
2264 z0 += kRB26s4SFlangeD2;
2265 shRB26s4SFlange->DefineSection(3, z0, kRB26s4SFlangeRi2, kRB26s4SFlangeRo);
2266 shRB26s4SFlange->DefineSection(4, z0, kRB26s4SFlangeRi3, kRB26s4SFlangeRo);
2267 z0 += kRB26s4SFlangeD3;
2268 shRB26s4SFlange->DefineSection(5, z0, kRB26s4SFlangeRi3, kRB26s4SFlangeRo);
2269 TGeoVolume* voRB26s4SFlange = new TGeoVolume("RB26s4SFlange", shRB26s4SFlange, kMedSteel);
2270
2271 TGeoVolume* voRB26s4SFlangeM = new TGeoVolume("RB26s4SFlangeM", MakeMotherFromTemplate(shRB26s4SFlange, 0, 3), kMedVac);
2272 voRB26s4SFlangeM->AddNode(voRB26s4SFlange, 1, gGeoIdentity);
2273
2274 ///////////////////////////////////////
2275 // RB26/5 Rotable Flange //
2276 // Drawing LHCVFX__0009 //
2277 ///////////////////////////////////////
2278 const Float_t kRB26s5RFlangeL = 1.86; // Length of the flange
2279 const Float_t kRB26s5RFlangeD1 = 0.61; // Length of section 1
2280 const Float_t kRB26s5RFlangeD2 = 0.15; // Length of section 2
2281 const Float_t kRB26s5RFlangeD3 = 0.60; // Length of section 3
2282 const Float_t kRB26s5RFlangeD4 = 0.50; // Length of section 4
2283 const Float_t kRB26s5RFlangeRo = 15.20/2.; // Flange outer radius
2284 const Float_t kRB26s5RFlangeRi1 = 10.30/2.; // Flange inner radius section 1
2285 const Float_t kRB26s5RFlangeRi2 = 10.00/2.; // Flange inner radius section 2
2286 const Float_t kRB26s5RFlangeRi3 = 10.30/2.; // Flange inner radius section 3
2287 const Float_t kRB26s5RFlangeRi4 = 10.50/2.; // Flange inner radius section 4
2288
2289 z0 = 0;
2290 TGeoPcon* shRB26s5RFlange = new TGeoPcon(0., 360., 8);
2291 z0 = 0.;
2292 shRB26s5RFlange->DefineSection(0, z0, kRB26s5RFlangeRi4, kRB26s5RFlangeRo);
2293 z0 += kRB26s5RFlangeD4;
2294 shRB26s5RFlange->DefineSection(1, z0, kRB26s5RFlangeRi4, kRB26s5RFlangeRo);
2295 shRB26s5RFlange->DefineSection(2, z0, kRB26s5RFlangeRi3, kRB26s5RFlangeRo);
2296 z0 += kRB26s5RFlangeD3;
2297 shRB26s5RFlange->DefineSection(3, z0, kRB26s5RFlangeRi3, kRB26s5RFlangeRo);
2298 shRB26s5RFlange->DefineSection(4, z0, kRB26s5RFlangeRi2, kRB26s5RFlangeRo);
2299 z0 += kRB26s5RFlangeD2;
2300 shRB26s5RFlange->DefineSection(5, z0, kRB26s5RFlangeRi2, kRB26s5RFlangeRo);
2301 shRB26s5RFlange->DefineSection(6, z0, kRB26s5RFlangeRi1, kRB26s5RFlangeRo);
2302 z0 += kRB26s5RFlangeD1;
2303 shRB26s5RFlange->DefineSection(7, z0, kRB26s5RFlangeRi1, kRB26s5RFlangeRo);
2304 TGeoVolume* voRB26s5RFlange = new TGeoVolume("RB26s5RFlange", shRB26s5RFlange, kMedSteel);
2305
2306 TGeoVolume* voRB26s5RFlangeM = new TGeoVolume("RB26s5RFlangeM", MakeMotherFromTemplate(shRB26s5RFlange, 4, 7), kMedVac);
2307 voRB26s5RFlangeM->AddNode(voRB26s5RFlange, 1, gGeoIdentity);
2308
2309 //
2310 // Assemble RB26/1-2
2311 //
2312 TGeoVolumeAssembly* asRB26s12 = new TGeoVolumeAssembly("RB26s12");
2313 z0 = 0.;
2314 asRB26s12->AddNode(voRB26s1RFlange, 1, gGeoIdentity);
2315 z0 += kRB26s1RFlangeIsL + kRB26s1RFlangeFpL;
2316 asRB26s12->AddNode(voRB26s12TubeM, 1, new TGeoTranslation(0., 0., z0));
2317 z0 += kRB26s12TubeL;
2318 asRB26s12->AddNode(voRB26s2Compensator, 1, new TGeoTranslation(0., 0., z0));
2319 z0 += kRB26s2CompL;
2320 z0 -= kRB26s2FFlangeD1;
2321 asRB26s12->AddNode(voRB26s2FFlangeM, 1, new TGeoTranslation(0., 0., z0));
2322 z0 += kRB26s2FFlangeL;
2323 const Float_t kRB26s12L = z0;
2324
2325 //
2326 // Assemble RB26/3
2327 //
2328 TGeoVolumeAssembly* asRB26s3 = new TGeoVolumeAssembly("RB26s3");
2329 z0 = 0.;
2330 asRB26s3->AddNode(voRB26s3SFlangeM, 1, gGeoIdentity);
2331 z0 += kRB26s3SFlangeL;
2332 z0 -= kRB26s3SFlangeD3;
2333 asRB26s3->AddNode(voRB26s3FixedPointM, 1, new TGeoTranslation(0., 0., z0));
2334 z0 += kRB26s3FixedPointL;
2335 asRB26s3->AddNode(voRB26s3TubeM, 1, new TGeoTranslation(0., 0., z0));
2336 z0 += kRB26s3TubeL;
2337 asRB26s3->AddNode(voRB26s3Compensator, 1, new TGeoTranslation(0., 0., z0));
2338 z0 += kRB26s3CompL;
2339 z0 -= kRB26s3FFlangeD1;
2340 asRB26s3->AddNode(voRB26s3FFlangeM, 1, new TGeoTranslation(0., 0., z0));
2341 z0 += kRB26s3FFlangeL;
2342 const Float_t kRB26s3L = z0;
2343
2344
2345 //
2346 // Assemble RB26/4-5
2347 //
2348 TGeoVolumeAssembly* asRB26s45 = new TGeoVolumeAssembly("RB26s45");
2349 z0 = 0.;
2350 asRB26s45->AddNode(voRB26s4SFlangeM, 1, gGeoIdentity);
2351 z0 += kRB26s4SFlangeL;
2352 z0 -= kRB26s4SFlangeD3;
2353 asRB26s45->AddNode(voRB26s4FixedPointM, 1, new TGeoTranslation(0., 0., z0));
2354 z0 += kRB26s4FixedPointL;
2355 asRB26s45->AddNode(voRB26s45TubeM, 1, new TGeoTranslation(0., 0., z0));
2356 z0 += kRB26s45TubeL;
2357 asRB26s45->AddNode(voRB26s5Compensator, 1, new TGeoTranslation(0., 0., z0));
2358 z0 += kRB26s5CompL;
2359 z0 -= kRB26s5RFlangeD3;
2360 z0 -= kRB26s5RFlangeD4;
2361 asRB26s45->AddNode(voRB26s5RFlangeM, 1, new TGeoTranslation(0., 0., z0));
2362 z0 += kRB26s5RFlangeL;
2363 const Float_t kRB26s45L = z0;
2364
2365 //
2366 // Assemble RB26
2367 //
2368 TGeoVolumeAssembly* asRB26Pipe = new TGeoVolumeAssembly("RB26Pipe");
2369 z0 = 0.;
2370 asRB26Pipe->AddNode(asRB26s12, 1, new TGeoTranslation(0., 0., z0));
2371 z0 += kRB26s12L;
2372 asRB26Pipe->AddNode(asRB26s3, 1, new TGeoTranslation(0., 0., z0));
2373 z0 += kRB26s3L;
2374 asRB26Pipe->AddNode(asRB26s45, 1, new TGeoTranslation(0., 0., z0));
2375 z0 += kRB26s45L;
2376 top->AddNode(asRB26Pipe, 1, new TGeoCombiTrans(0., 0., -82., rot180));
2377}
2378
2379
2380
2381//___________________________________________
2382void AliPIPEv4::CreateMaterials()
2383{
2384 //
2385 // Define materials for beam pipe
2386 //
2387
2388 AliDebugClass(1,"Create PIPEv4 materials");
2389 Int_t isxfld = ((AliMagF*)TGeoGlobalMagField::Instance()->GetField())->Integ();
2390 Float_t sxmgmx = ((AliMagF*)TGeoGlobalMagField::Instance()->GetField())->Max();
2391 // Steel (Inox)
2392 Float_t asteel[4] = { 55.847,51.9961,58.6934,28.0855 };
2393 Float_t zsteel[4] = { 26.,24.,28.,14. };
2394 Float_t wsteel[4] = { .715,.18,.1,.005 };
2395 // AlBe - alloy
af5f3976 2396 Float_t aAlBe[2] = { 26.98, 9.01}; // al=2.702 be=1.8477
820b4d9e 2397 Float_t zAlBe[2] = { 13.00, 4.00};
2398 Float_t wAlBe[2] = { 0.4, 0.6};
2399 //
2400 // Polyamid
2401 Float_t aPA[4] = {16., 14., 12., 1.};
2402 Float_t zPA[4] = { 8., 7., 6., 1.};
2403 Float_t wPA[4] = { 1., 1., 6., 11.};
2404 //
236bcb7a 2405 // Polyimide film
2406 Float_t aPI[4] = {16., 14., 12., 1.};
2407 Float_t zPI[4] = { 8., 7., 6., 1.};
2408 Float_t wPI[4] = { 5., 2., 22., 10.};
2409 // Rohacell
2410 Float_t aRohacell[4] = {16., 14., 12., 1.};
2411 Float_t zRohacell[4] = { 8., 7., 6., 1.};
2412 Float_t wRohacell[4] = { 2., 1., 9., 13.};
2413 //
820b4d9e 2414 // Air
2415 //
2416 Float_t aAir[4]={12.0107,14.0067,15.9994,39.948};
2417 Float_t zAir[4]={6.,7.,8.,18.};
2418 Float_t wAir[4]={0.000124,0.755267,0.231781,0.012827};
2419 Float_t dAir = 1.20479E-3;
2420 Float_t dAir1 = 1.20479E-10;
2421 //
2422 // Insulation powder
2423 // Si O Ti Al
2424 Float_t ains[4] ={28.0855, 15.9994, 47.867, 26.982};
2425 Float_t zins[4] ={14., 8. , 22. , 13. };
2426 Float_t wins[4] ={ 0.3019, 0.4887, 0.1914, 0.018};
2427 //
2428 //
2429 // Anticorodal
2430 //
2431 // Al Si7 Mg 0.6
2432 //
2433 Float_t aaco[3] ={26.982, 28.0855, 24.035};
2434 Float_t zaco[3] ={13., 14. , 12. };
2435 Float_t waco[3] ={ 0.924, 0.07, 0.006};
2436 // Kapton
2437 //
2438 Float_t aKapton[4]={1.00794,12.0107, 14.010,15.9994};
2439 Float_t zKapton[4]={1.,6.,7.,8.};
2440 Float_t wKapton[4]={0.026362,0.69113,0.07327,0.209235};
2441 Float_t dKapton = 1.42;
2442 // NEG coating
2443 // Ti V Zr
2444 Float_t aNEG[4] = {47.87, 50.94, 91.24};
2445 Float_t zNEG[4] = {22.00, 23.00, 40.00};
2446 Float_t wNEG[4] = {1./3., 1./3., 1./3.};
2447 Float_t dNEG = 5.6; // ?
2448
af5f3976 2449 //---------------------------------
2450 // Aluminium AA 5083 for MFT: Al Manganese(Mn) Magnesium(Mg) Chrome(Cr)
2451 Float_t aALU5083[4]={26.982, 54.938, 24.305, 51.996}; // Mg pas meme a que la ligne Anticorodal!
2452 Float_t zALU5083[4] ={13., 25., 12., 24.};
2453 Float_t wALU5083[4] ={0.947, 0.007, 0.044, 0.0015};
2454 // Aluminium AA 2219 for MFT: Al Cu Mn Ti V Zr
2455 Float_t aALU2219[6]={26.982, 63.546, 54.938, 47.867, 50.941, 91.224};
2456 Float_t zALU2219[6] ={13., 29., 25., 22., 23., 40.};
2457 Float_t wALU2219[6] ={0.93, 0.063, 0.003, 0.0006, 0.001, 0.0018};
2458 //---------------------------------
2459
820b4d9e 2460 //
af5f3976 2461 // Silicon for ITS UPGRADE
2462 AliMaterial(2, "SILICON$",28.09 , 14.00 , 2.33 , 9.36 , 45.);
2463
820b4d9e 2464 //
2465 // Berillium
2466 AliMaterial(5, "BERILLIUM$", 9.01, 4., 1.848, 35.3, 36.7);
2467 //
2468 // Carbon
2469 AliMaterial(6, "CARBON$ ", 12.01, 6., 2.265, 18.8, 49.9);
2470 //
2471 // Aluminum
2472 AliMaterial(9, "ALUMINIUM$", 26.98, 13., 2.7, 8.9, 37.2);
2473 //
2474 // Copper
2475 AliMaterial(10, "COPPER", 63.55, 29, 8.96, 1.43, 85.6/8.96);
2476 //
2477 // Air
2478 AliMixture(15, "AIR$ ", aAir, zAir, dAir, 4, wAir);
2479 AliMixture(35, "AIR_HIGH$ ", aAir, zAir, dAir, 4, wAir);
2480 //
2481 // Vacuum
2482 AliMixture(16, "VACUUM$ ", aAir, zAir, dAir1, 4, wAir);
2483 //
2484 // stainless Steel
2485 AliMixture(19, "STAINLESS STEEL$", asteel, zsteel, 7.88, 4, wsteel);
2486 //
2487 // reduced density steel to approximate pump getter material
2488 AliMixture(20, "GETTER$", asteel, zsteel, 1.00, 4, wsteel);
2489 // Al-Be alloy
2490 //
af5f3976 2491 AliMixture(21, "AlBe$", aAlBe, zAlBe, 2.07, 2, wAlBe);
820b4d9e 2492 // Polyamid
2493 //
2494 AliMixture(22, "PA$", aPA, zPA, 1.14, -4, wPA);
2495 //
2496 // Kapton
2497 AliMixture(23, "KAPTON", aKapton, zKapton, dKapton, 4, wKapton);
2498 // Anticorodal
2499 AliMixture(24, "ANTICORODAL", aaco, zaco, 2.66, 3, waco);
2500
2501 //
2502 // Insulation powder
2503 AliMixture(14, "INSULATION0$", ains, zins, 0.41, 4, wins);
2504 AliMixture(34, "INSULATION1$", ains, zins, 0.41, 4, wins);
2505 AliMixture(54, "INSULATION2$", ains, zins, 0.41, 4, wins);
2506
2507 // NEG
2508 AliMixture(25, "NEG COATING", aNEG, zNEG, dNEG, -3, wNEG);
af5f3976 2509
2510 //---------------------------------
2511 // Aluminium AA5083 for MFT
2512 AliMixture(63, "ALUMINIUM5083$",aALU5083,zALU5083, 2.66 ,4,wALU5083); // from aubertduval.fr
2513 // Aluminium AA2219 for MFT
2514 AliMixture(64, "ALUMINIUM2219$",aALU2219,zALU2219, 2.84 ,6,wALU2219); // from aubertduval.fr
2515 //---------------------------------
236bcb7a 2516 // Polyimide Film
2517 //
2518 AliMixture(65, "PI$", aPI, zPI, 1.42, -4, wPI);
2519 //---------------------------------
2520 // Carbon Fiber M55J
2521 AliMaterial(66,"M55J6K$",12.0107,6,1.92,999,999);
2522 // Rohacell C9 H13 N1 O2 0.03 g/cm^3
2523 AliMixture(67,"Rohacell$", aRohacell, zRohacell, 0.03, -4, wRohacell);
2524
2525 // ****************
820b4d9e 2526 // Defines tracking media parameters.
2527 //
2528 Float_t epsil = .001; // Tracking precision,
2529 Float_t stemax = -0.01; // Maximum displacement for multiple scat
2530 Float_t tmaxfd = -20.; // Maximum angle due to field deflection
2531 Float_t deemax = -.3; // Maximum fractional energy loss, DLS
2532 Float_t stmin = -.8;
2533 // ***************
2534 //
af5f3976 2535 // Silicon for ITS UPGRADE
2536 AliMedium(2, "SILICON", 2, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
2537
2538
820b4d9e 2539 // Beryllium
2540
2541 AliMedium(5, "BE", 5, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
2542
2543 // Carbon
2544 AliMedium(6, "C", 6, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
2545 //
2546 // Aluminum
2547 AliMedium(9, "ALU", 9, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
2548 // Copper
2549 AliMedium(10, "CU", 10, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
2550 //
2551 // Air
2552 AliMedium(15, "AIR", 15, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
2553 AliMedium(35, "AIR_HIGH",35, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
2554 //
2555 // Vacuum
2556 AliMedium(16, "VACUUM", 16, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
2557 //
2558 // Steel
2559 AliMedium(19, "INOX", 19, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
2560 //
2561 // Getter
2562 AliMedium(20, "GETTER", 20, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
2563 //
2564 // AlBe - Aloy
2565 AliMedium(21, "AlBe" , 21, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
2566 //
2567 // Polyamid
2568 AliMedium(22, "PA" , 22, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
2569 // Antocorodal
2570 AliMedium(24, "ANTICORODAL", 24, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
2571 // Insulation Powder
2572 AliMedium(14, "INS_C0 ", 14, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
2573 AliMedium(34, "INS_C1 ", 34, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
2574 AliMedium(54, "INS_C2 ", 54, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
2575 //
2576 // KAPTON
2577 AliMedium(23, "KAPTON", 23, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
2578
2579 //
2580 // NEG
2581 AliMedium(25, "NEG COATING", 25, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
af5f3976 2582
2583 //----------------- for the MFT ----------------------
2584 AliMedium(63,"AA5083", 63, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
2585 AliMedium(64,"AA2219", 64, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
2586 //----------------------------------------------------
236bcb7a 2587 AliMedium(65,"POLYIMIDE", 65, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
2588 //---------------------------------
2589 // Carbon Fiber M55J
2590 AliMedium(66, "M55J6K",66,0,isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
2591 // Rohacell
2592 AliMedium(67,"ROHACELL",67,0,isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
af5f3976 2593
2594
820b4d9e 2595}
2596
2597
2598TGeoPcon* AliPIPEv4::MakeMotherFromTemplate(TGeoPcon* shape, Int_t imin, Int_t imax, Float_t r0, Int_t nz)
2599{
2600 //
2601 // Create a mother shape from a template setting some min radii to 0
2602 //
2603 Int_t nz0 = shape->GetNz();
2604 // if nz > -1 the number of planes is given by nz
2605 if (nz != -1) nz0 = nz;
2606 TGeoPcon* mother = new TGeoPcon(0., 360., nz0);
2607
2608 if (imin == -1 || imax == -1) {
2609 imin = 0;
2610 imax = shape->GetNz();
2611 } else if (imax >= nz0) {
2612 imax = nz0 - 1;
2613 printf("Warning: imax reset to nz-1 %5d %5d %5d %5d\n", imin, imax, nz, nz0);
2614 }
2615
2616
2617
2618 for (Int_t i = 0; i < shape->GetNz(); i++) {
2619 Double_t rmin = shape->GetRmin(i);
2620 if ((i >= imin) && (i <= imax) ) rmin = r0;
2621 Double_t rmax = shape->GetRmax(i);
2622 Double_t z = shape->GetZ(i);
2623 mother->DefineSection(i, z, rmin, rmax);
2624 }
2625 return mother;
2626
2627}
2628
2629TGeoPcon* AliPIPEv4::MakeInsulationFromTemplate(TGeoPcon* shape)
2630{
2631 //
2632 // Create an beam pipe insulation layer shape from a template
2633 //
2634 Int_t nz = shape->GetNz();
2635 TGeoPcon* insu = new TGeoPcon(0., 360., nz);
2636
2637 for (Int_t i = 0; i < nz; i++) {
2638 Double_t z = shape->GetZ(i);
2639 Double_t rmin = shape->GetRmin(i);
2640 Double_t rmax = shape->GetRmax(i);
2641 rmax += 0.5;
2642 shape->DefineSection(i, z, rmin, rmax);
2643 rmin = rmax - 0.5;
2644 insu->DefineSection(i, z, rmin, rmax);
2645 }
2646 return insu;
2647
2648}
2649
2650
2651TGeoVolume* AliPIPEv4::MakeBellow(const char* ext, Int_t nc, Float_t rMin, Float_t rMax, Float_t dU, Float_t rPlie, Float_t dPlie)
2652{
2653 // nc Number of convolution
2654 // rMin Inner radius of the bellow
2655 // rMax Outer radius of the bellow
2656 // dU Undulation length
2657 // rPlie Plie radius
2658 // dPlie Plie thickness
2659 const TGeoMedium* kMedVac = gGeoManager->GetMedium("PIPE_VACUUM");
af5f3976 2660 //const TGeoMedium* kMedSteel = gGeoManager->GetMedium("PIPE_INOX");
24679816 2661 const TGeoMedium* kMedAlu5083 = gGeoManager->GetMedium("PIPE_AA5083"); //fm
820b4d9e 2662
2663 char name[64], nameA[64], nameB[64], bools[64];
267ba80c 2664 snprintf(name, 64, "%sBellowUS", ext);
820b4d9e 2665 TGeoVolume* voBellow = new TGeoVolume(name, new TGeoTube(rMin, rMax, dU/2.), kMedVac);
2666 //
2667 // Upper part of the undulation
2668 //
2669 TGeoTorus* shPlieTorusU = new TGeoTorus(rMax - rPlie, rPlie - dPlie, rPlie);
267ba80c 2670 snprintf(nameA, 64, "%sTorusU", ext);
820b4d9e 2671 shPlieTorusU->SetName(nameA);
2672 TGeoTube* shPlieTubeU = new TGeoTube (rMax - rPlie, rMax, rPlie);
267ba80c 2673 snprintf(nameB, 64, "%sTubeU", ext);
820b4d9e 2674 shPlieTubeU->SetName(nameB);
267ba80c 2675 snprintf(name, 64, "%sUpperPlie", ext);
2676 snprintf(bools, 64, "%s*%s", nameA, nameB);
820b4d9e 2677 TGeoCompositeShape* shUpperPlie = new TGeoCompositeShape(name, bools);
2678
af5f3976 2679 TGeoVolume* voWiggleU = new TGeoVolume(name, shUpperPlie, kMedAlu5083);
24679816 2680 voWiggleU->SetLineColor(kOrange); // fm
2681
820b4d9e 2682 // Lower part of the undulation
2683 TGeoTorus* shPlieTorusL = new TGeoTorus(rMin + rPlie, rPlie - dPlie, rPlie);
267ba80c 2684 snprintf(nameA, 64, "%sTorusL", ext);
820b4d9e 2685 shPlieTorusL->SetName(nameA);
2686 TGeoTube* shPlieTubeL = new TGeoTube (rMin, rMin + rPlie, rPlie);
267ba80c 2687 snprintf(nameB, 64, "%sTubeL", ext);
820b4d9e 2688 shPlieTubeL->SetName(nameB);
267ba80c 2689 snprintf(name, 64, "%sLowerPlie", ext);
2690 snprintf(bools, 64, "%s*%s", nameA, nameB);
820b4d9e 2691 TGeoCompositeShape* shLowerPlie = new TGeoCompositeShape(name, bools);
2692
af5f3976 2693 TGeoVolume* voWiggleL = new TGeoVolume(name, shLowerPlie, kMedAlu5083);
24679816 2694 voWiggleL->SetLineColor(kOrange); // fm
2695
820b4d9e 2696 // Connection between upper and lower part of undulation
267ba80c 2697 snprintf(name, 64, "%sPlieConn1", ext);
af5f3976 2698 TGeoVolume* voWiggleC1 = new TGeoVolume(name, new TGeoTube(rMin + rPlie, rMax - rPlie, dPlie/2.), kMedAlu5083);
24679816 2699 voWiggleC1->SetLineColor(kOrange); // fm
2700
820b4d9e 2701 // One wiggle
2702 Float_t dz = rPlie - dPlie / 2.;
2703 Float_t z0 = - dPlie / 2.;
267ba80c 2704 snprintf(name, 64, "%sWiggle", ext);
820b4d9e 2705 TGeoVolumeAssembly* asWiggle = new TGeoVolumeAssembly(name);
af5f3976 2706
820b4d9e 2707 asWiggle->AddNode(voWiggleC1, 1 , new TGeoTranslation(0., 0., z0));
2708 z0 += dz;
2709 asWiggle->AddNode(voWiggleU, 1 , new TGeoTranslation(0., 0., z0));
2710 z0 += dz;
2711 asWiggle->AddNode(voWiggleC1, 2 , new TGeoTranslation(0., 0., z0));
2712 z0 += dz;
2713 asWiggle->AddNode(voWiggleL , 1 , new TGeoTranslation(0., 0., z0));
2714 // Positioning of the volumes
2715 z0 = - dU / 2.+ rPlie;
24679816 2716 ////////////voBellow->AddNode(voWiggleL, 2, new TGeoTranslation(0., 0., z0)); removing the first 1/2 plie, fm
820b4d9e 2717 z0 += rPlie;
2718 Float_t zsh = 4. * rPlie - 2. * dPlie;
2719 for (Int_t iw = 0; iw < nc; iw++) {
2720 Float_t zpos = z0 + iw * zsh;
af5f3976 2721 voBellow->AddNode(asWiggle, iw + 1, new TGeoTranslation(0., 0., zpos - dPlie));
2722
820b4d9e 2723 }
2724 return voBellow;
2725}
2726
236bcb7a 2727TGeoVolume* AliPIPEv4::MakeBellowCside(const char* ext, Int_t nc, Float_t rMin, Float_t rMax, Float_t rPlie, Float_t dPlie)
04a81bb0
AM
2728{
2729 // nc Number of convolution
2730 // rMin Inner radius of the bellow
2731 // rMax Outer radius of the bellow
2732 // dU Undulation length
2733 // rPlie Plie radius
2734 // dPlie Plie thickness
2735 const TGeoMedium* kMedVac = gGeoManager->GetMedium("PIPE_VACUUM");
2736 //const TGeoMedium* kMedSteel = gGeoManager->GetMedium("PIPE_INOX");
2737 const TGeoMedium* kMedAlu5083 = gGeoManager->GetMedium("PIPE_AA5083"); //fm
2738
236bcb7a 2739 Float_t dU = nc * (4.*rPlie - 2. *dPlie);
2740
04a81bb0
AM
2741 char name[64], nameA[64], nameB[64], bools[64];
2742 snprintf(name, 64, "%sBellowUS", ext);
236bcb7a 2743// TGeoVolume* voBellow = new TGeoVolume(name, new TGeoTube(rMin, rMax, dU/2.), kMedVac);
2744 TGeoVolumeAssembly *voBellow = new TGeoVolumeAssembly(name);
04a81bb0
AM
2745 //
2746 // Upper part of the undulation
2747 //
236bcb7a 2748
04a81bb0
AM
2749 TGeoTorus* shPlieTorusU = new TGeoTorus(rMax - rPlie, rPlie - dPlie, rPlie);
2750 snprintf(nameA, 64, "%sTorusU", ext);
2751 shPlieTorusU->SetName(nameA);
2752 TGeoTube* shPlieTubeU = new TGeoTube (rMax - rPlie, rMax, rPlie);
2753 snprintf(nameB, 64, "%sTubeU", ext);
2754 shPlieTubeU->SetName(nameB);
2755 snprintf(name, 64, "%sUpperPlie", ext);
2756 snprintf(bools, 64, "%s*%s", nameA, nameB);
2757 TGeoCompositeShape* shUpperPlie = new TGeoCompositeShape(name, bools);
2758
2759 TGeoVolume* voWiggleU = new TGeoVolume(name, shUpperPlie, kMedAlu5083);
2760 voWiggleU->SetLineColor(kOrange); // fm
2761
236bcb7a 2762 // First Lower part of the ondulation
04a81bb0
AM
2763 TGeoTorus* shPlieTorusL = new TGeoTorus(rMin + rPlie, rPlie - dPlie, rPlie);
2764 snprintf(nameA, 64, "%sTorusL", ext);
2765 shPlieTorusL->SetName(nameA);
2766 TGeoTranslation *t1 = new TGeoTranslation("t1",0,0,-rPlie/2.);
2767 t1->RegisterYourself();
2768
2769 TGeoTube* shPlieTubeL = new TGeoTube (rMin, rMin + rPlie, rPlie/2.);
2770 snprintf(nameB, 64, "%sTubeL", ext);
2771 shPlieTubeL->SetName(nameB);
2772 snprintf(name, 64, "%sLowerPlie", ext);
2773 snprintf(bools, 64, "%s*%s:t1", nameA, nameB);
2774 TGeoCompositeShape* shLowerPlie1 = new TGeoCompositeShape(name, bools);
2775
2776 TGeoVolume* voWiggleL1 = new TGeoVolume(name, shLowerPlie1, kMedAlu5083);
2777 voWiggleL1->SetLineColor(kOrange); // fm
2778
2779 // Second Lower part of the undulation
2780 TGeoTranslation *t2 = new TGeoTranslation("t2",0,0,rPlie/2.);
2781 t2->RegisterYourself();
2782
2783 snprintf(bools, 64, "%s*%s:t2", nameA, nameB);
2784 TGeoCompositeShape* shLowerPlie2 = new TGeoCompositeShape(name, bools);
2785
2786 TGeoVolume* voWiggleL2 = new TGeoVolume(name, shLowerPlie2, kMedAlu5083);
2787 voWiggleL2->SetLineColor(kOrange); // fm
2788
2789 // Connection between upper and lower part of undulation
2790 snprintf(name, 64, "%sPlieConn1", ext);
2791 TGeoVolume* voWiggleC1 = new TGeoVolume(name, new TGeoTube(rMin + rPlie, rMax - rPlie, dPlie/2.), kMedAlu5083);
2792 voWiggleC1->SetLineColor(kOrange); // fm
236bcb7a 2793
2794 //
2795 // Vacuum Part
2796 //
2797
2798 //--Upper part of the ondulation
2799
2800 TGeoTorus* vacPlieTorusU = new TGeoTorus(rMax - rPlie, 0., rPlie- dPlie);
2801 snprintf(nameA, 64, "%svacTorusU", ext);
2802 vacPlieTorusU->SetName(nameA);
2803 TGeoTube* vacPlieTubeU = new TGeoTube (0., rMax- rPlie, rPlie-dPlie);
2804 snprintf(nameB, 64, "%svacTubeU", ext);
2805 vacPlieTubeU->SetName(nameB);
2806 snprintf(name, 64, "%svacUpperPlie", ext);
2807 snprintf(bools, 64, "%s+%s", nameA, nameB);
2808 TGeoCompositeShape* vacUpperPlie = new TGeoCompositeShape(name, bools);
2809
2810 TGeoVolume* voVacWiggleU = new TGeoVolume(name, vacUpperPlie, kMedVac);
2811 voVacWiggleU->SetVisibility(0);
2812
2813
2814 // First Lower part of the undulation
2815 TGeoTorus* vacPlieTorusL = new TGeoTorus(rMin + rPlie, 0., rPlie);
2816 snprintf(nameA, 64, "%svacTorusL", ext);
2817 vacPlieTorusL->SetName(nameA);
2818
2819 TGeoTube* vacPlieTubeL = new TGeoTube (0., rMin + rPlie, rPlie/2.);
2820 snprintf(nameB, 64, "%svacTubeL", ext);
2821 vacPlieTubeL->SetName(nameB);
2822 snprintf(name, 64, "%svacLowerPlie", ext);
2823 snprintf(bools, 64, "%s:t1-%s", nameB, nameA);
2824 TGeoCompositeShape* vacLowerPlie1 = new TGeoCompositeShape(name, bools);
2825
2826 TGeoVolume* voVacWiggleL1 = new TGeoVolume(name, vacLowerPlie1, kMedVac);
2827 voVacWiggleL1->SetVisibility(0);
2828
2829
2830 // Second Lower part of the undulation
2831
2832 snprintf(bools, 64, "%s:t2-%s", nameB, nameA);
2833 TGeoCompositeShape* vacLowerPlie2 = new TGeoCompositeShape(name, bools);
2834
2835 TGeoVolume* voVacWiggleL2 = new TGeoVolume(name, vacLowerPlie2, kMedVac);
2836 voVacWiggleL2->SetVisibility(0);
2837
04a81bb0
AM
2838
2839 // One wiggle
2840 Float_t dz = rPlie - dPlie / 2.;
2841 Float_t z0 = 2.*rPlie;
2842 snprintf(name, 64, "%sWiggle", ext);
2843 TGeoVolumeAssembly* asWiggle = new TGeoVolumeAssembly(name);
2844
2845 asWiggle->AddNode(voWiggleL1 , 1 , new TGeoTranslation(0., 0., z0));
236bcb7a 2846 asWiggle->AddNode(voVacWiggleL1 , 1 , new TGeoTranslation(0., 0., z0));
04a81bb0
AM
2847 z0 -= dz;
2848 asWiggle->AddNode(voWiggleC1, 1 , new TGeoTranslation(0., 0., z0));
2849 z0 -= dz;
2850 asWiggle->AddNode(voWiggleU, 1 , new TGeoTranslation(0., 0., z0));
236bcb7a 2851 asWiggle->AddNode(voVacWiggleU, 1 , new TGeoTranslation(0., 0., z0));
04a81bb0
AM
2852 z0 -= dz;
2853 asWiggle->AddNode(voWiggleC1, 2 , new TGeoTranslation(0., 0., z0));
2854 z0 -= dz;
2855 asWiggle->AddNode(voWiggleL2 , 1 , new TGeoTranslation(0., 0., z0));
236bcb7a 2856 asWiggle->AddNode(voVacWiggleL2 , 1 , new TGeoTranslation(0., 0., z0));
2857
2858
2859
2860
04a81bb0
AM
2861 // Positioning of the volumes
2862 z0 = + dU / 2.;
2863 Float_t zsh = 4. * dz;
236bcb7a 2864 //for (Int_t iw = 0; iw < 1; iw++) {
2865 for (Int_t iw = 0; iw < nc; iw++) {
04a81bb0
AM
2866 Float_t zpos = z0 - iw * zsh;
2867 voBellow->AddNode(asWiggle, iw + 1, new TGeoTranslation(0., 0., zpos));
2868
2869 }
2870 return voBellow;
2871}
236bcb7a 2872//TGeoVolume* AliPIPEv4::MakeBellowCside(const char* ext, Int_t nc, Float_t rMin, Float_t rMax, Float_t dU, Float_t rPlie, Float_t dPlie)
2873//{
2874// // nc Number of convolution
2875// // rMin Inner radius of the bellow
2876// // rMax Outer radius of the bellow
2877// // dU Undulation length
2878// // rPlie Plie radius
2879// // dPlie Plie thickness
2880// const TGeoMedium* kMedVac = gGeoManager->GetMedium("PIPE_VACUUM");
2881// //const TGeoMedium* kMedSteel = gGeoManager->GetMedium("PIPE_INOX");
2882// const TGeoMedium* kMedAlu5083 = gGeoManager->GetMedium("PIPE_AA5083"); //fm
2883//
2884// char name[64], nameA[64], nameB[64], bools[64];
2885// snprintf(name, 64, "%sBellowUS", ext);
2886// TGeoVolume* voBellow = new TGeoVolume(name, new TGeoTube(rMin, rMax, dU/2.), kMedVac);
2887// //
2888// // Upper part of the undulation
2889// //
2890//
2891// TGeoTorus* shPlieTorusU = new TGeoTorus(rMax - rPlie, rPlie - dPlie, rPlie);
2892// snprintf(nameA, 64, "%sTorusU", ext);
2893// shPlieTorusU->SetName(nameA);
2894// TGeoTube* shPlieTubeU = new TGeoTube (rMax - rPlie, rMax, rPlie);
2895// snprintf(nameB, 64, "%sTubeU", ext);
2896// shPlieTubeU->SetName(nameB);
2897// snprintf(name, 64, "%sUpperPlie", ext);
2898// snprintf(bools, 64, "%s*%s", nameA, nameB);
2899// TGeoCompositeShape* shUpperPlie = new TGeoCompositeShape(name, bools);
2900//
2901// TGeoVolume* voWiggleU = new TGeoVolume(name, shUpperPlie, kMedAlu5083);
2902// voWiggleU->SetLineColor(kOrange); // fm
2903//
2904// // First Lower part of the undulation
2905// TGeoTorus* shPlieTorusL = new TGeoTorus(rMin + rPlie, rPlie - dPlie, rPlie);
2906// snprintf(nameA, 64, "%sTorusL", ext);
2907// shPlieTorusL->SetName(nameA);
2908// TGeoTranslation *t1 = new TGeoTranslation("t1",0,0,-rPlie/2.);
2909// t1->RegisterYourself();
2910//
2911// TGeoTube* shPlieTubeL = new TGeoTube (rMin, rMin + rPlie, rPlie/2.);
2912// snprintf(nameB, 64, "%sTubeL", ext);
2913// shPlieTubeL->SetName(nameB);
2914// snprintf(name, 64, "%sLowerPlie", ext);
2915// snprintf(bools, 64, "%s*%s:t1", nameA, nameB);
2916// TGeoCompositeShape* shLowerPlie1 = new TGeoCompositeShape(name, bools);
2917//
2918// TGeoVolume* voWiggleL1 = new TGeoVolume(name, shLowerPlie1, kMedAlu5083);
2919// voWiggleL1->SetLineColor(kOrange); // fm
2920//
2921// // Second Lower part of the undulation
2922// TGeoTranslation *t2 = new TGeoTranslation("t2",0,0,rPlie/2.);
2923// t2->RegisterYourself();
2924//
2925// snprintf(bools, 64, "%s*%s:t2", nameA, nameB);
2926// TGeoCompositeShape* shLowerPlie2 = new TGeoCompositeShape(name, bools);
2927//
2928// TGeoVolume* voWiggleL2 = new TGeoVolume(name, shLowerPlie2, kMedAlu5083);
2929// voWiggleL2->SetLineColor(kOrange); // fm
2930//
2931// // Connection between upper and lower part of undulation
2932// snprintf(name, 64, "%sPlieConn1", ext);
2933// TGeoVolume* voWiggleC1 = new TGeoVolume(name, new TGeoTube(rMin + rPlie, rMax - rPlie, dPlie/2.), kMedAlu5083);
2934// voWiggleC1->SetLineColor(kOrange); // fm
2935//
2936// // One wiggle
2937// Float_t dz = rPlie - dPlie / 2.;
2938// Float_t z0 = 2.*rPlie;
2939// snprintf(name, 64, "%sWiggle", ext);
2940// TGeoVolumeAssembly* asWiggle = new TGeoVolumeAssembly(name);
2941//
2942// asWiggle->AddNode(voWiggleL1 , 1 , new TGeoTranslation(0., 0., z0));
2943// // z0 -= dz;
2944// // asWiggle->AddNode(voWiggleC1, 1 , new TGeoTranslation(0., 0., z0));
2945// // z0 -= dz;
2946// // asWiggle->AddNode(voWiggleU, 1 , new TGeoTranslation(0., 0., z0));
2947// // z0 -= dz;
2948// // asWiggle->AddNode(voWiggleC1, 2 , new TGeoTranslation(0., 0., z0));
2949// // z0 -= dz;
2950// // asWiggle->AddNode(voWiggleL2 , 1 , new TGeoTranslation(0., 0., z0));
2951// // Positioning of the volumes
2952// z0 = + dU / 2.;
2953// Float_t zsh = 4. * dz;
2954// for (Int_t iw = 0; iw < 1; iw++) {
2955// // for (Int_t iw = 0; iw < nc; iw++) {
2956// Float_t zpos = z0 - iw * zsh;
2957// voBellow->AddNode(asWiggle, iw + 1, new TGeoTranslation(0., 0., zpos));
2958//
2959// }
2960// return voBellow;
2961//}
04a81bb0 2962
820b4d9e 2963