Coding convetion violations fixed.
[u/mrichter/AliRoot.git] / STRUCT / AliDIPOv3.cxx
CommitLineData
7778505e 1/**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3 * *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
6 * *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
15
16/* $Id$ */
17
18#include <TVirtualMC.h>
19#include <TGeoManager.h>
20#include <TGeoVolume.h>
21#include <TGeoMedium.h>
22#include <TGeoMatrix.h>
23#include <TGeoBBox.h>
24#include <TGeoTube.h>
25#include <TGeoCone.h>
26#include <TGeoPcon.h>
d4916cbf 27#include <TGeoCompositeShape.h>
7778505e 28
29#include "AliConst.h"
30#include "AliDIPOv3.h"
31#include "AliMagF.h"
32#include "AliRun.h"
33
34ClassImp(AliDIPOv3)
35
36//_____________________________________________________________________________
37AliDIPOv3::AliDIPOv3()
38{
39 //
d4916cbf 40 // Last design of magnetic dipole version 3
7778505e 41 //
42}
43
44//_____________________________________________________________________________
45AliDIPOv3::AliDIPOv3(const char *name, const char *title)
cf0c8fa1 46 : AliDIPOv2(name,title)
7778505e 47{
48 //
d4916cbf 49 // Standard constructor for the magnetic dipole version 3
7778505e 50}
51
52
53//_____________________________________________________________________________
54void AliDIPOv3::CreateSpectrometerDipole()
55{
56// Detailed dipole geometry as built
57//
d4916cbf 58// Drawing: ALIP2A__0026
59// Geometer measurements: EDMS 596079
60// EDMS 584963
61
7778505e 62//
63// The top volume
64//
65 TGeoVolume* top = gGeoManager->GetVolume("ALIC");
66//
67// Media
68//
d4916cbf 69 TGeoMedium* kMedSteel = gGeoManager->GetMedium("DIPO_ST_C3");
70 TGeoMedium* kMedCoil = gGeoManager->GetMedium("DIPO_Coil_C1");
71 TGeoMedium* kMedCoilSh = gGeoManager->GetMedium("DIPO_Coil_C3");
72 TGeoMedium* kMedCable = gGeoManager->GetMedium("DIPO_ALU_C2");
73 TGeoMedium* kMedAlu = gGeoManager->GetMedium("DIPO_ALU_C2");
74 TGeoMedium* kMedAir = gGeoManager->GetMedium("DIPO_AIR_MUON");
7778505e 75//
d4916cbf 76// Rotations
7778505e 77//
d4916cbf 78 Float_t alhc = 0.794;
79
cf0c8fa1 80 TGeoRotation* rotxz = new TGeoRotation("rotxz", 270., 0., 90., 90., 180., 0.);
81 TGeoRotation* rotiz = new TGeoRotation("rotiz", 90., 0., 90., 90., 180., 0.);
d4916cbf 82 TGeoRotation* rotxzlhc = new TGeoRotation("rotxzlhc", 180., 180. + alhc, 0.);
d4916cbf 83
7778505e 84 TGeoRotation* rotxz108 = new TGeoRotation("rotxz108", 90., 108., 90., 198., 180., 0.);
85 TGeoRotation* rotxz180 = new TGeoRotation("rotxz180", 90., 180., 90., 270., 180., 0.);
86 TGeoRotation* rotxz288 = new TGeoRotation("rotxz288", 90., 288., 90., 18., 180., 0.);
87
88 TGeoRotation* rotxy180 = new TGeoRotation("rotxy180", 90., 180., 90., 270., 0., 0.);
89 TGeoRotation* rotxy108 = new TGeoRotation("rotxy108", 90., 108., 90., 198., 0., 0.);
90 TGeoRotation* rotxy288 = new TGeoRotation("rotxy288", 90., 288., 90., 18., 0., 0.);
91
d4916cbf 92 TGeoRotation* rot00 = new TGeoRotation("rot00", 180., 0., 90., 151., 90., 61.);
93 TGeoRotation* rot01 = new TGeoRotation("rot01", 180., 0., 90., 29.,- 90., -61.);
94 TGeoRotation* rot02 = new TGeoRotation("rot02", 0., 0., 90., 151., 90., 61.);
95 TGeoRotation* rot03 = new TGeoRotation("rot03", 0., 0., 90., 29.,- 90., -61.);
96 TGeoRotation* rot04 = new TGeoRotation("rot04", 90., 61., 90., 151., 0., 0.);
97 TGeoRotation* rot05 = new TGeoRotation("rot05", 90., -61., 90.,-151., 0., 0.);
98 TGeoRotation* rot06 = new TGeoRotation("rot06", 90., 119., 90., 209., 0., 0.);
99 TGeoRotation* rot07 = new TGeoRotation("rot07", 90.,-119., 90.,-209., 0., 0.);
100
101 const Float_t dipoleL = 498.;
102 const Float_t kZDipoleR = 1244.;
103 const Float_t kZDipole = kZDipoleR - dipoleL/2.;
104 const Float_t kZDipoleF = kZDipoleR - dipoleL;
105 const Float_t yokeLength = 309.4;
106 const Float_t blockLength = yokeLength / 7.;
107 const Float_t gapWidthFront = 297.6;
108 const Float_t gapWidthRear = 395.4;
109 const Float_t dGap = (gapWidthRear - gapWidthFront) / 12.;
110 const Float_t gapHeight = 609.1;
111 const Float_t blockHeight = 145.45;
112 const Float_t dzCoil = 4.45;
7778505e 113
7778505e 114 Float_t dx, dy, dz;
d4916cbf 115
116
7778505e 117 //
cf0c8fa1 118 // Mother volume for muon spectrometer tracking station 3
119
120 Float_t z30 = 825.;
121 Float_t zst = 1052.;
122
123 Float_t rcD0 = (kZDipoleF - 5.) * TMath::Tan(9. * kDegrad);
d4916cbf 124 Float_t rcD1 = kZDipole * TMath::Tan(9. * kDegrad);
d4916cbf 125 Float_t rcD2 = rcD1 + dipoleL/2. * TMath::Tan(10.1 * kDegrad);
cf0c8fa1 126 Float_t rc30 = z30 * TMath::Tan(9. * kDegrad);
127 Float_t rcst = rcD1 + (zst - kZDipole) * TMath::Tan(10.1 * kDegrad);
128
129
130 Float_t riD0 = (kZDipoleF - 5.) * TMath::Tan(2. * kDegrad) + 0.2;
131 Float_t riD1 = 30.;
132 Float_t riD2 = 35.8;
133 Float_t riD3 = riD2 + (kZDipoleR - zst) * TMath::Tan(2. * kDegrad);
134 Float_t riD4 = riD2 + (kZDipoleR - zst + 5.) * TMath::Tan(2. * kDegrad);
135
136
137 TGeoPcon* shDDIP1 = new TGeoPcon("shDDIP1", 0., 360., 7);
138
139 shDDIP1->DefineSection(0, (kZDipoleF - 5.), riD0, rcD0);
140 shDDIP1->DefineSection(1, z30 , riD1, rc30);
141 shDDIP1->DefineSection(2, kZDipole , riD1, rcD1);
142 shDDIP1->DefineSection(3, zst , riD1, rcst);
143 shDDIP1->DefineSection(4, zst , riD2, rcst);
144 shDDIP1->DefineSection(5, kZDipoleR , riD3, rcD2);
145 shDDIP1->DefineSection(6, (kZDipoleR + 5.), riD4, rcD2);
d4916cbf 146
cf0c8fa1 147 TGeoBBox* shDDIP2 = new TGeoBBox(164., 182., 36.);
5c8f914e 148 shDDIP2->SetName("shDDIP2");
c52bffc9 149 TGeoTranslation* trDDIP2 = new TGeoTranslation("trDDIP2", 0., 0., kZDipole - 12.);
cf0c8fa1 150 trDDIP2->RegisterYourself();
151
c52bffc9 152 TGeoTube* shDDIP3 = new TGeoTube(0., 30., 40.);
153 shDDIP3->SetName("shDDIP3");
d4916cbf 154
c52bffc9 155
156 TGeoCompositeShape* shDDIP = new TGeoCompositeShape("shDDIP", "shDDIP1+(shDDIP2:trDDIP2-shDDIP3:trDDIP2)");
d4916cbf 157 TGeoVolume* voDDIP = new TGeoVolume("DDIP", shDDIP, kMedAir);
7778505e 158//
159// Yoke
160//
7778505e 161
162 TGeoVolumeAssembly* asYoke = new TGeoVolumeAssembly("DYoke");
163// Base
164 char name[32];
165 Float_t lx0 = gapWidthFront + 2. * blockHeight;
166 Float_t lx = lx0;
167
168 TGeoVolumeAssembly* asYokeBase = new TGeoVolumeAssembly("DYokeBase");
169 for (Int_t i = 0; i < 7; i++) {
170 sprintf(name, "DYokeBaseBlock%1d", i);
171 TGeoVolume* voBaseBlock = new TGeoVolume(name,
172 new TGeoBBox(lx/2., blockHeight/2., blockLength/2.),
173 kMedSteel);
174 asYokeBase->AddNode(voBaseBlock, 1, new TGeoTranslation(0., 0., Float_t(i - 3) * blockLength));
175 lx += 2. * dGap;
176 }
177
178 asYoke->AddNode(asYokeBase, 1, new TGeoTranslation(0., -(gapHeight + blockHeight)/2. , 0.));
179 asYoke->AddNode(asYokeBase, 2, new TGeoTranslation(0., +(gapHeight + blockHeight)/2. , 0.));
180
181
182// Side Wall
183 TGeoVolumeAssembly* asYokeSide = new TGeoVolumeAssembly("DYokeSide");
184 TGeoVolume* voSideBlock = new TGeoVolume("DSideBlock",
185 new TGeoBBox(blockHeight/2., gapHeight/2., blockLength/2.),
186 kMedSteel);
187
188 for (Int_t i = 0; i < 7; i++) {
189 asYokeSide->AddNode(voSideBlock, i, new TGeoTranslation(Float_t(i - 3) * dGap, 0., Float_t(i - 3) * blockLength));
190 }
191
192
193 asYoke->AddNode(asYokeSide, 1, new TGeoTranslation(+lx0/2. + 3. * dGap - blockHeight/2., 0., 0.));
cf0c8fa1 194 asYoke->AddNode(asYokeSide, 2, new TGeoCombiTrans( -lx0/2. - 3. * dGap + blockHeight/2., 0., 0., rotiz));
7778505e 195
196//
197// Coils
198//
d4916cbf 199 Float_t coilRi = 206.;
7778505e 200 Float_t coilD = 70.;
201 Float_t coilRo = coilRi + coilD;
202 Float_t coilH = 77.;
203 Float_t phiMin = -61.;
204 Float_t phiMax = 61.;
d4916cbf 205 Float_t lengthSt = 240. + 33.9;
7778505e 206 Float_t phiKnee = phiMax * kDegrad;
d4916cbf 207 Float_t rKnee = 31.5;
7778505e 208
209// Circular sections
210 TGeoVolumeAssembly* asCoil = new TGeoVolumeAssembly("DCoil");
211
212 TGeoVolume* voDC1 = new TGeoVolume("DC1",
213 new TGeoTubeSeg(coilRi, coilRo, coilH / 2., phiMin, phiMax),
214 kMedCoil);
215 TGeoVolume* voDC2 = new TGeoVolume("DC2",
216 new TGeoTubeSeg(coilRi + 5., coilRo - 5., coilH / 2., phiMin, phiMax),
217 kMedCoilSh);
218
219 voDC1->AddNode(voDC2, 1, gGeoIdentity);
220 voDC2->SetVisibility(0);
221
222 dz = lengthSt / 2. + coilH / 2. + rKnee;
223 dx = 0.;
224
225 asCoil->AddNode(voDC1, 1, new TGeoTranslation(-dx, 0., -dz));
226 asCoil->AddNode(voDC1, 2, new TGeoCombiTrans( dx, 0., -dz, rotxy180));
227 asCoil->AddNode(voDC1, 3, new TGeoTranslation(-dx, 0., dz));
228 asCoil->AddNode(voDC1, 4, new TGeoCombiTrans( dx, 0., dz, rotxz180));
229
230
231// 90deg Knees
232
233
234 TGeoVolume* voDC11 = new TGeoVolume("DC11",
235 new TGeoTubeSeg(rKnee, rKnee + coilH, coilD/2., 270., 360.),
236 kMedCoil);
237
238
239 dx = - TMath::Cos(phiKnee) * (coilRi + coilD/2.);
240 dy = - TMath::Sin(phiKnee) * (coilRi + coilD/2.);
241 dz = lengthSt / 2.;
242
243 asCoil->AddNode(voDC11, 1, new TGeoCombiTrans( dx, dy, -dz, rot00));
244 asCoil->AddNode(voDC11, 2, new TGeoCombiTrans( dx, dy, dz, rot02));
245 asCoil->AddNode(voDC11, 3, new TGeoCombiTrans(-dx, dy, -dz, rot01));
246 asCoil->AddNode(voDC11, 4, new TGeoCombiTrans(-dx, dy, dz, rot03));
247
248 TGeoVolume* voDC12 = new TGeoVolume("DC12",
249 new TGeoTubeSeg(rKnee, rKnee + coilH, coilD/2., 0., 90.),
250 kMedCoil);
251
252
253 asCoil->AddNode(voDC12, 1, new TGeoCombiTrans( dx, -dy, -dz, rot01));
254 asCoil->AddNode(voDC12, 2, new TGeoCombiTrans( dx, -dy, dz, rot03));
255 asCoil->AddNode(voDC12, 3, new TGeoCombiTrans(-dx, -dy, -dz, rot00));
256 asCoil->AddNode(voDC12, 4, new TGeoCombiTrans(-dx, -dy, dz, rot02));
257
258// Straight sections
259
260
261 TGeoVolume* voDL0 = new TGeoVolume("DL0",
262 new TGeoBBox(coilD / 2. + 2., coilH / 2. + 2., lengthSt / 2.),
263 kMedCoil);
264
265 TGeoVolume* voDL1 = new TGeoVolume("DL1",
266 new TGeoBBox(coilD / 2., coilH / 2., lengthSt / 2.),
267 kMedCoil);
268
269
270 TGeoVolume* voDL2 = new TGeoVolume("DL2",
271 new TGeoBBox(coilD / 2. - 5., coilH / 2. - 5., lengthSt / 2. - 5.),
272 kMedCoilSh);
273 // Sleeves
274 TGeoVolume* voDL3 = new TGeoVolume("DL3",
275 new TGeoBBox(1., coilH / 2., 120.),
276 kMedAlu);
277
278 TGeoVolume* voDL4 = new TGeoVolume("DL4",
279 new TGeoBBox(coilD/2., 1., 120.),
280 kMedAlu);
281
282 voDL0->SetVisibility(0);
283 voDL1->AddNode(voDL2, 1, gGeoIdentity);
284 voDL0->AddNode(voDL1, 1, gGeoIdentity);
285 voDL0->AddNode(voDL3, 1, new TGeoTranslation(-coilD/2. - 1., 0., 0.));
286 voDL0->AddNode(voDL3, 2, new TGeoTranslation(+coilD/2. + 1., 0., 0.));
287 voDL0->AddNode(voDL4, 1, new TGeoTranslation(0., -coilH/2. - 1., 0.));
288 voDL0->AddNode(voDL4, 2, new TGeoTranslation(0., +coilH/2. + 1., 0.));
289
290
291 dx += (rKnee + coilH/2.) * TMath::Sin(phiKnee);
292 dy -= (rKnee + coilH/2.) * TMath::Cos(phiKnee);
293 dz = 0.;
294
295 asCoil->AddNode(voDL0, 1, new TGeoCombiTrans( dx, dy, dz, rot04));
296 asCoil->AddNode(voDL0, 2, new TGeoCombiTrans( dx, -dy, dz, rot05));
297 asCoil->AddNode(voDL0, 3, new TGeoCombiTrans(-dx, dy, dz, rot06));
298 asCoil->AddNode(voDL0, 4, new TGeoCombiTrans(-dx, -dy, dz, rot07));
299
300// Contactor
301// Outer face planes
302
303 TGeoVolumeAssembly* asContactor = new TGeoVolumeAssembly("DContactor");
304 dx = -5.;
305 TGeoVolume* voDC10 = new TGeoVolume("DC10",
d4916cbf 306 new TGeoTubeSeg(coilRo + 5.1, coilRo + 73.5, 1., -20., 20.),
7778505e 307 kMedCable);
308 asContactor->AddNode(voDC10, 1, new TGeoTranslation(dx, 0, -32.325));
309 asContactor->AddNode(voDC10, 2, new TGeoTranslation(dx, 0, +32.325));
310
311
312// Coil Support
313//
d4916cbf 314 Float_t sW = 83.;
7778505e 315
316 TGeoVolumeAssembly* asDCoilSupport = new TGeoVolumeAssembly("DCoilSupport");
317
318 // Steel fixed to the yoke
319 TGeoVolume* voDCS01 = new TGeoVolume("DCS01",
320 new TGeoTubeSeg(coilRo, 325., 1., 21., 51.),
321 kMedAlu);
322
323 // Steel on the coil
324 TGeoVolume* voDCS02 = new TGeoVolume("DCS02",
d4916cbf 325 new TGeoTubeSeg(coilRo, coilRo + 3.125, sW/2., 21., 51.),
7778505e 326 kMedAlu);
327 TGeoVolume* voDCS021 = new TGeoVolume("DCS021",
968388a6 328 new TGeoConeSeg(sW/2., coilRo + 3.124, 320., coilRo + 3.125, coilRo + 5.125, 21., 21.4),
7778505e 329 kMedAlu);
330
331
332 // Sleeves
333 TGeoVolume* voDCS03 = new TGeoVolume("DCS03",
d4916cbf 334 new TGeoTubeSeg(coilRi - 3.125, coilRo + 3.125, 3.125/2., 21., 51.),
7778505e 335 kMedAlu);
336
337 TGeoVolume* voDCS04 = new TGeoVolume("DCS04",
d4916cbf 338 new TGeoTubeSeg(coilRi - 3.125, coilRi, coilH/2., 21., 51.),
7778505e 339 kMedAlu);
340
341
342 TGeoVolume* voDCS05 = new TGeoVolume("DCS05",
d4916cbf 343 new TGeoTubeSeg(coilRi - 3.125, coilRo, 3.125/2., 21., 51.),
7778505e 344 kMedAlu);
345 //
7778505e 346 asDCoilSupport->AddNode(voDCS02, 1, new TGeoTranslation(0., 0., -(sW - coilH)/2.));
347 asDCoilSupport->AddNode(voDCS04, 1, gGeoIdentity);
348 for (Int_t i = 0; i < 9; i++)
349 {
3b7df2c4 350 char nameR[16];
351 sprintf(nameR, "rotdcs%1d", i);
7778505e 352 Float_t phi = Float_t(i) * 3.75;
3b7df2c4 353 TGeoRotation* rot = new TGeoRotation(nameR, 90., phi, 90., 90. + phi, 0., 0.);
dd6154b2 354 asDCoilSupport->AddNode(voDCS021, i, new TGeoCombiTrans(0., 0.004, -(sW - coilH)/2., rot));
7778505e 355 }
356
357
358
d4916cbf 359 asDCoilSupport->AddNode(voDCS01, 1, new TGeoTranslation(0., 0., -sW/2. - (sW - coilH)/2. - 3.125/2.));
360 asDCoilSupport->AddNode(voDCS03, 1, new TGeoTranslation(0., 0., +coilH/2. + 3.125/2.));
361 asDCoilSupport->AddNode(voDCS05, 1, new TGeoTranslation(0., 0., -coilH/2. - 3.125/2.));
7778505e 362
d4916cbf 363
364 //
365 // SAA1 Support: Hanger 1
366 //
367 TGeoTranslation* trHanger = new TGeoTranslation("trHanger", 0., 250., 0.);
368 trHanger->RegisterYourself();
7778505e 369
d4916cbf 370 Float_t rmin1, rmin2, rmax1, rmax2;
371
372 Float_t zHanger1 = 811.9;
5c8f914e 373 TGeoBBox* shHanger11 = new TGeoBBox(2.5/2., 250., 25./2.);
374 shHanger11->SetName("shHanger11");
375
d4916cbf 376 rmin1 = (zHanger1 - 13.) * TMath::Tan(2. * kDegrad);
377 rmin2 = rmin1 + 26. * TMath::Tan( 2.0 * kDegrad);
cf0c8fa1 378
d4916cbf 379 rmax1 = (zHanger1 - 13.) * TMath::Tan(9. * kDegrad);
380 rmax2 = rmax1 + 26. * TMath::Tan(9. * kDegrad);
381
5c8f914e 382 TGeoCone* shHanger12 = new TGeoCone(13., rmin1, rmax1, rmin2, rmax2);
383 shHanger12->SetName("shHanger12");
d4916cbf 384 TGeoCompositeShape* shHanger1 = new TGeoCompositeShape("shHanger1", "shHanger12*shHanger11:trHanger");
385 TGeoVolume* voHanger1 = new TGeoVolume("DHanger1", shHanger1, kMedSteel);
386 //
387 // SAA1 Support: Hanger 2
388 //
389 Float_t zHanger2 = 1171.9;
5c8f914e 390 TGeoBBox* shHanger21 = new TGeoBBox(3.5/2., 250., 25./2.);
391 shHanger21->SetName("shHanger21");
d4916cbf 392
cf0c8fa1 393 rmin1 = 35.8 + (zHanger2 - 13. - zst) * TMath::Tan(2. * kDegrad);
d4916cbf 394 rmin2 = rmin1 + 26. * TMath::Tan( 2.0 * kDegrad);
cf0c8fa1 395
d4916cbf 396 rmax1 = rcD1 + (zHanger2 - 13. - kZDipole) * TMath::Tan(10.1 * kDegrad);
397 rmax2 = rmax1 + 26. * TMath::Tan(10.1 * kDegrad);
5c8f914e 398 TGeoCone* shHanger22 = new TGeoCone(13., rmin1, rmax1, rmin2, rmax2);
399 shHanger22->SetName("shHanger22");
400
d4916cbf 401 TGeoCompositeShape* shHanger2 = new TGeoCompositeShape("shHanger2", "shHanger22*shHanger21:trHanger");
402
403 TGeoVolume* voHanger2 = new TGeoVolume("DHanger2", shHanger2, kMedSteel);
404 //
405 // Hanger support
406 Float_t hsLength = yokeLength + (zHanger2 - kZDipole - yokeLength/2.) + 25./2.;
407
408 TGeoVolume* voHS1 = new TGeoVolume("DHS1", new TGeoBBox( 1.5, 12.5, hsLength/2.), kMedSteel);
409 TGeoVolume* voHS2 = new TGeoVolume("DHS2", new TGeoBBox(12.5, 1.5, hsLength/2.), kMedSteel);
c52bffc9 410 Float_t hsH = gapHeight/2. + blockHeight - (rmax1+rmax2)/2. - 2.;
d4916cbf 411
412 TGeoVolume* voHS3 = new TGeoVolume("DHS3", new TGeoBBox(3.5/2., hsH/2., 25./2.), kMedSteel);
413
414 TGeoVolumeAssembly* asHS = new TGeoVolumeAssembly("asHS");
415 asHS->AddNode(voHS1, 1, gGeoIdentity);
416 asHS->AddNode(voHS2, 1, new TGeoTranslation(0., +14., 0.));
417 asHS->AddNode(voHS2, 2, new TGeoTranslation(0., -14., 0.));
ca3e9926 418 asHS->AddNode(voHS3, 1, new TGeoTranslation(0., -hsH/2. - 14. -1.5, hsLength/2. - 25./2.));
d4916cbf 419
420
421
cf0c8fa1 422 dz = zHanger1;
423 voDDIP->AddNode(voHanger1, 1, new TGeoTranslation(0., 0., dz));
d4916cbf 424
cf0c8fa1 425 dz = zHanger2;
426 voDDIP->AddNode(voHanger2, 1, new TGeoTranslation(0., 0., dz));
d4916cbf 427
428
429
430
431 // Assembly everything
432
433 TGeoVolumeAssembly* asDipole = new TGeoVolumeAssembly("Dipole");
434 // Yoke
435 asDipole->AddNode(asYoke, 1, new TGeoTranslation(0., 0., -dzCoil));
7778505e 436 asDipole->AddNode(asCoil, 1, gGeoIdentity);
d4916cbf 437 // Contactor
7778505e 438 dz = lengthSt / 2. + coilH / 2. + rKnee;
d4916cbf 439 asDipole->AddNode(asContactor, 1, new TGeoTranslation(0., 0., dz + dzCoil));
440 asDipole->AddNode(asContactor, 2, new TGeoCombiTrans( 0., 0., dz - dzCoil, rotxy180));
441 // Coil support
7778505e 442 asDipole->AddNode(asDCoilSupport, 1, new TGeoTranslation(0., 0., dz));
443 asDipole->AddNode(asDCoilSupport, 2, new TGeoCombiTrans( 0., 0., dz, rotxy180));
444 asDipole->AddNode(asDCoilSupport, 3, new TGeoCombiTrans( 0., 0., dz, rotxy108));
445 asDipole->AddNode(asDCoilSupport, 4, new TGeoCombiTrans( 0., 0., dz, rotxy288));
d4916cbf 446
cf0c8fa1 447 asDipole->AddNode(asDCoilSupport, 5, new TGeoCombiTrans( 0., 0., -dz, rotiz));
7778505e 448 asDipole->AddNode(asDCoilSupport, 6, new TGeoCombiTrans( 0., 0., -dz, rotxz108));
449 asDipole->AddNode(asDCoilSupport, 7, new TGeoCombiTrans( 0., 0., -dz, rotxz180));
450 asDipole->AddNode(asDCoilSupport, 8, new TGeoCombiTrans( 0., 0., -dz, rotxz288));
d4916cbf 451
452 // Hanger (Support)
453 dy = gapHeight/2. + blockHeight + 14.;
454
ca3e9926 455 asDipole->AddNode(asHS, 1, new TGeoTranslation(0., dy + 1.5, ((zHanger2 - kZDipole - yokeLength/2.) + 25./2.)/2.));
d4916cbf 456
457
458 asDipole->SetVisContainers(1);
7778505e 459 voDDIP->SetVisibility(0);
7778505e 460
cf0c8fa1 461 top->AddNode(asDipole, 1, new TGeoCombiTrans(0., dipoleL / 2. * TMath::Tan(alhc * kDegrad), -kZDipole, rotxzlhc));
462 top->AddNode(voDDIP, 1, new TGeoCombiTrans(0., 0., 0., rotxz));
7778505e 463}
464
465
466
467
468
469
470