Logo ROOT   6.30.04
Reference Guide
 All Namespaces Files Pages
station1.C
Go to the documentation of this file.
1 /// \file
2 /// \ingroup tutorial_geom
3 /// Drawing a space station, using ROOT geometry class.
4 ///
5 /// Reviewed by Sunman Kim (sunman98@hanmail.net)
6 /// Supervisor: Prof. Inkyu Park (icpark@physics.uos.ac.kr)
7 ///
8 /// How to run: `.x station1.C` in ROOT terminal, then use OpenGL
9 ///
10 /// This macro was created for the evaluation of Computational Physics course in 2006.
11 /// We thank to Prof. Inkyu Park for his special lecture on ROOT and to all of ROOT team
12 ///
13 /// \macro_code
14 ///
15 /// \author Chang Yeol Lee, Dept. of Physics, Univ. of Seoul
16 
17 #include "TGeoManager.h"
18 
19 void station1()
20 {
21 
22 
23  TGeoManager *geom=new TGeoManager("geom","My first 3D geometry");
24 
25  TGeoMaterial *vacuum=new TGeoMaterial("vacuum",0,0,0);
26  TGeoMaterial *Fe=new TGeoMaterial("Fe",55.845,26,7.87);
27  TGeoMaterial *Cu=new TGeoMaterial("Cu",63.549,29,8.92);
28 
29  TGeoMedium *Air=new TGeoMedium("Vacuum",0,vacuum);
30  TGeoMedium *Iron=new TGeoMedium("Iron",1,Fe);
31  TGeoMedium *Copper=new TGeoMedium("Copper",2,Cu);
32 
33  TGeoVolume *top=geom->MakeBox("top",Air,1000,1000,1000);
34  geom->SetTopVolume(top);
35  geom->SetTopVisible(0);
36  // If you want to see the boundary, please input the number, 1 instead of 0.
37  // Like this, geom->SetTopVisible(1);
38 
39 
40 
41  TGeoVolume *Cone1=geom->MakeCone("Cone1",Copper,650,0,20,0,20);
42  Cone1->SetFillColor(35);
43  Cone1->SetLineColor(35);
44  top->AddNodeOverlap(Cone1,1,new TGeoTranslation(0,0,0));
45 
46  TGeoVolume *Cone2=geom->MakeCone("Cone2",Copper,25,0,30,0,30);
47  Cone2->SetFillColor(7);
48  Cone2->SetLineColor(7);
49  top->AddNodeOverlap(Cone2,1,new TGeoTranslation(0,0,630));
50 
51  TGeoVolume *Cone21=geom->MakeCone("Cone21",Copper,30,0,30,0,30);
52  Cone21->SetFillColor(29);
53  Cone21->SetLineColor(29);
54  top->AddNodeOverlap(Cone21,1,new TGeoTranslation(0,0,550));
55 
56  TGeoVolume *Cone22=geom->MakeCone("Cone22",Copper,5,0,50,0,50);
57  Cone22->SetFillColor(2);
58  Cone22->SetLineColor(2);
59  top->AddNodeOverlap(Cone22,1,new TGeoTranslation(0,0,500));
60 
61 for(int i=0;i<28;i++){
62  TGeoVolume *Cone00=geom->MakeCone("Cone00",Copper,3,0,25,0,25);
63  Cone00->SetFillColor(1);
64  Cone00->SetLineColor(1);
65  top->AddNodeOverlap(Cone00,1,new TGeoTranslation(0,0,-100+20*i));
66 }
67 
68  TGeoVolume *Cone3=geom->MakeCone("Cone3",Copper,60,0,70,0,0);
69  Cone3->SetFillColor(13);
70  Cone3->SetLineColor(13);
71  top->AddNodeOverlap(Cone3,1,new TGeoTranslation(-60,0,-110));
72 
73  TGeoVolume *Cone31=geom->MakeCone("Cone31",Copper,230,0,70,0,70);
74  Cone31->SetFillColor(13);
75  Cone31->SetLineColor(13);
76  top->AddNodeOverlap(Cone31,1,new TGeoTranslation(-60,0,-400));
77 
78 for(int i=0;i<5;i++){
79  Cone31=geom->MakeCone("Cone31",Copper,7,0,73,0,73);
80  Cone31->SetFillColor(21);
81  Cone31->SetLineColor(21);
82  top->AddNodeOverlap(Cone31,1,new TGeoTranslation(-60,0,-170-(500/6*(i+1))));
83  top->AddNodeOverlap(Cone31,1,new TGeoTranslation(60,0,-170-(500/6*(i+1))));
84  top->AddNodeOverlap(Cone31,1,new TGeoTranslation(0,-60,-170-(500/6*(i+1))));
85  top->AddNodeOverlap(Cone31,1,new TGeoTranslation(0,60,-170-(500/6*(i+1))));
86 }
87 
88  TGeoVolume *Cone32=geom->MakeCone("Cone32",Copper,30,60,50,0,70);
89  Cone32->SetFillColor(35);
90  Cone32->SetLineColor(35);
91  top->AddNodeOverlap(Cone32,1,new TGeoTranslation(-60,0,-650));
92 
93  TGeoVolume *Cone321=geom->MakeCone("Cone321",Copper,5,60,50,0,50);
94  Cone321->SetFillColor(2);
95  Cone321->SetLineColor(2);
96  top->AddNodeOverlap(Cone321,1,new TGeoTranslation(-60,0,-680));
97 
98  TGeoVolume *Cone4=geom->MakeCone("Cone4",Copper,60,0,70,0,0);
99  Cone4->SetFillColor(13);
100  Cone4->SetLineColor(13);
101  top->AddNodeOverlap(Cone3,1,new TGeoTranslation(60,0,-110));
102 
103 for(int i=1;i<=8;i++){
104  TGeoVolume *Torus2=geom->MakeTorus("Torus2",Iron,120,20,40,45*i-4,8);
105  Torus2->SetFillColor(18);
106  Torus2->SetLineColor(18);
107  top->AddNodeOverlap(Torus2,1,new TGeoTranslation(0,0,610));
108 
109  TGeoVolume *Tubs=geom->MakeTubs("Line",Iron,0,190,5,45*i-1,45*i+1);
110  Tubs->SetFillColor(18);
111  Tubs->SetLineColor(18);
112  top->AddNodeOverlap(Tubs,1,new TGeoTranslation(0,0,610));
113 }
114 
115  TGeoVolume *Cone41=geom->MakeCone("Cone41",Copper,230,0,70,0,70);
116  Cone41->SetFillColor(13);
117  Cone41->SetLineColor(13);
118  top->AddNodeOverlap(Cone41,1,new TGeoTranslation(60,0,-400));
119 
120  TGeoVolume *Cone42=geom->MakeCone("Cone42",Copper,30,60,50,0,70);
121  Cone42->SetFillColor(35);
122  Cone42->SetLineColor(35);
123  top->AddNodeOverlap(Cone42,1,new TGeoTranslation(60,0,-650));
124 
125  TGeoVolume *Cone421=geom->MakeCone("Cone421",Copper,5,60,50,0,50);
126  Cone421->SetFillColor(2);
127  Cone421->SetLineColor(2);
128  top->AddNodeOverlap(Cone421,1,new TGeoTranslation(60,0,-680));
129 
130  TGeoVolume *Cone5=geom->MakeCone("Cone5",Copper,60,0,70,0,0);
131  Cone5->SetFillColor(13);
132  Cone5->SetLineColor(13);
133  top->AddNodeOverlap(Cone3,1,new TGeoTranslation(0,-60,-110));
134 
135  TGeoVolume *Cone51=geom->MakeCone("Cone51",Copper,230,0,70,0,70);
136  Cone51->SetFillColor(13);
137  Cone51->SetLineColor(13);
138  top->AddNodeOverlap(Cone51,1,new TGeoTranslation(0,-60,-400));
139 
140  TGeoVolume *Cone52=geom->MakeCone("Cone52",Copper,30,60,50,0,70);
141  Cone52->SetFillColor(35);
142  Cone52->SetLineColor(35);
143  top->AddNodeOverlap(Cone52,1,new TGeoTranslation(0,-60,-650));
144 
145  TGeoVolume *Cone521=geom->MakeCone("Cone521",Copper,5,60,50,0,50);
146  Cone521->SetFillColor(2);
147  Cone521->SetLineColor(2);
148  top->AddNodeOverlap(Cone521,1,new TGeoTranslation(0,-60,-680));
149 
150  TGeoVolume *Cone6=geom->MakeCone("Cone6",Copper,60,0,70,0,0);
151  Cone6->SetFillColor(13);
152  Cone6->SetLineColor(13);
153  top->AddNodeOverlap(Cone3,1,new TGeoTranslation(0,60,-110));
154 
155  TGeoVolume *Cone61=geom->MakeCone("Cone61",Copper,230,0,70,0,70);
156  Cone61->SetFillColor(13);
157  Cone61->SetLineColor(13);
158  top->AddNodeOverlap(Cone61,1,new TGeoTranslation(0,60,-400));
159 
160  TGeoVolume *Cone62=geom->MakeCone("Cone62",Copper,30,60,50,0,70);
161  Cone62->SetFillColor(35);
162  Cone62->SetLineColor(35);
163  top->AddNodeOverlap(Cone62,1,new TGeoTranslation(0,60,-650));
164 
165  TGeoVolume *Cone621=geom->MakeCone("Cone621",Copper,5,60,50,0,50);
166  Cone621->SetFillColor(2);
167  Cone621->SetLineColor(2);
168  top->AddNodeOverlap(Cone621,1,new TGeoTranslation(0,60,-680));
169 
170  TGeoVolume *Cone7=geom->MakeCone("Cone7",Copper,50,0,40,0,5);
171  Cone7->SetFillColor(13);
172  Cone7->SetLineColor(13);
173  top->AddNodeOverlap(Cone7,1,new TGeoCombiTrans(-90,-60,10,new TGeoRotation("Cone7",90,-90,-90)));
174 
175  TGeoVolume *Cone71=geom->MakeCone("Cone71",Copper,50,0,60,0,40);
176  Cone71->SetFillColor(16);
177  Cone71->SetLineColor(16);
178  top->AddNodeOverlap(Cone71,1,new TGeoCombiTrans(10,-60,10,new TGeoRotation("Cone7",90,-90,-90)));
179 
180  TGeoVolume *Cone711=geom->MakeCone("Cone711",Copper,10,0,10,0,60);
181  Cone711->SetFillColor(13);
182  Cone711->SetLineColor(13);
183  top->AddNodeOverlap(Cone711,1,new TGeoCombiTrans(70,-60,10,new TGeoRotation("Cone7",90,-90,-90)));
184 
185  TGeoVolume *Torus1=geom->MakeTorus("Torus1",Iron,120,30,20);
186  Torus1->SetFillColor(33);
187  Torus1->SetLineColor(33);
188  top->AddNodeOverlap(Torus1,1,new TGeoTranslation(0,0,610));
189 
190  TGeoVolume *Cone8=geom->MakeCone("Cone8",Copper,50,0,40,0,5);
191  Cone8->SetFillColor(13);
192  Cone8->SetLineColor(13);
193  top->AddNodeOverlap(Cone8,1,new TGeoCombiTrans(100,60,10,new TGeoRotation("Cone8",90,90,0)));
194 
195  TGeoVolume *Cone81=geom->MakeCone("Cone81",Copper,50,0,60,0,40);
196  Cone81->SetFillColor(16);
197  Cone81->SetLineColor(16);
198  top->AddNodeOverlap(Cone81,1,new TGeoCombiTrans(0,60,10,new TGeoRotation("Cone8",90,90,0)));
199 
200  TGeoVolume *Cone811=geom->MakeCone("Cone811",Copper,10,0,10,0,60);
201  Cone811->SetFillColor(13);
202  Cone811->SetLineColor(13);
203  top->AddNodeOverlap(Cone811,1,new TGeoCombiTrans(-60,60,10,new TGeoRotation("Cone8",90,90,0)));
204 
205  TGeoVolume *Box1=geom->MakeBox("Box1",Copper,10,10,3);
206  Box1->SetFillColor(3);
207  Box1->SetLineColor(3);
208  top->AddNodeOverlap(Box1,1,new TGeoCombiTrans(-110,-50,645,new TGeoRotation("Box1",0,0,30)));
209 
210  TGeoVolume *Box2=geom->MakeBox("Box2",Copper,10,10,3);
211  Box2->SetFillColor(3);
212  Box2->SetLineColor(3);
213  top->AddNodeOverlap(Box2,1,new TGeoCombiTrans(110,45,645,new TGeoRotation("Box2",0,0,30)));
214 
215  TGeoVolume *Box3=geom->MakeBox("Box3",Copper,10,10,3);
216  Box3->SetFillColor(3);
217  Box3->SetLineColor(3);
218  top->AddNodeOverlap(Box3,1,new TGeoCombiTrans(-45,-110,645,new TGeoRotation("Box3",0,0,70)));
219 
220  TGeoVolume *Box4=geom->MakeBox("Box4",Copper,10,10,3);
221  Box4->SetFillColor(3);
222  Box4->SetLineColor(3);
223  top->AddNodeOverlap(Box4,1,new TGeoCombiTrans(45,110,645,new TGeoRotation("Box4",0,0,70)));
224 
225  TGeoVolume *Box5=geom->MakeBox("Box5",Copper,10,10,3);
226  Box5->SetFillColor(3);
227  Box5->SetLineColor(3);
228  top->AddNodeOverlap(Box5,1,new TGeoCombiTrans(45,-110,645,new TGeoRotation("Box5",0,0,30)));
229 
230  TGeoVolume *Box6=geom->MakeBox("Box6",Copper,10,10,3);
231  Box6->SetFillColor(3);
232  Box6->SetLineColor(3);
233  top->AddNodeOverlap(Box6,1,new TGeoCombiTrans(-45,110,645,new TGeoRotation("Box6",0,0,25)));
234 
235  TGeoVolume *Box7=geom->MakeBox("Box7",Copper,10,10,3);
236  Box7->SetFillColor(3);
237  Box7->SetLineColor(3);
238  top->AddNodeOverlap(Box7,1,new TGeoCombiTrans(110,-50,645,new TGeoRotation("Box7",0,0,60)));
239 
240  TGeoVolume *Box8=geom->MakeBox("Box8",Copper,10,10,3);
241  Box8->SetFillColor(3);
242  Box8->SetLineColor(3);
243  top->AddNodeOverlap(Box8,1,new TGeoCombiTrans(-110,45,645,new TGeoRotation("Box8",0,0,60)));
244 
245  Torus1=geom->MakeTorus("Torus1",Iron,120,30,20);
246  Torus1->SetFillColor(33);
247  Torus1->SetLineColor(33);
248  top->AddNodeOverlap(Torus1,1,new TGeoTranslation(0,0,610));
249 
250 for(int i=1;i<=8;i++){
251  TGeoVolume *Torus2=geom->MakeTorus("Torus2",Iron,120,20,40,45*i-4,8);
252  Torus2->SetFillColor(18);
253  Torus2->SetLineColor(18);
254  top->AddNodeOverlap(Torus2,1,new TGeoTranslation(0,0,610));
255 
256  TGeoVolume *Tubs=geom->MakeTubs("Line",Iron,0,190,5,45*i-1,45*i+1);
257  Tubs->SetFillColor(18);
258  Tubs->SetLineColor(18);
259  top->AddNodeOverlap(Tubs,1,new TGeoTranslation(0,0,610));
260 }
261 
262  TGeoVolume *Sphere00=geom->MakeSphere("Sphere00",Iron,0,15,0,45,0);
263  Sphere00->SetFillColor(2);
264  Sphere00->SetLineColor(2);
265  top->AddNodeOverlap(Sphere00,1,new TGeoTranslation(-145,-145,600));
266 
267  TGeoVolume *Sphere01=geom->MakeSphere("Sphere01",Iron,0,15,0,45,0);
268  Sphere01->SetFillColor(2);
269  Sphere01->SetLineColor(2);
270  top->AddNodeOverlap(Sphere01,1,new TGeoTranslation(0,-210,600));
271 
272  TGeoVolume *Sphere02=geom->MakeSphere("Sphere02",Iron,0,15,0,45,0);
273  Sphere02->SetFillColor(2);
274  Sphere02->SetLineColor(2);
275  top->AddNodeOverlap(Sphere02,1,new TGeoTranslation(145,145,600));
276 
277  TGeoVolume *Sphere03=geom->MakeSphere("Sphere03",Iron,0,15,0,45,0);
278  Sphere03->SetFillColor(2);
279  Sphere03->SetLineColor(2);
280  top->AddNodeOverlap(Sphere03,1,new TGeoTranslation(0,210,600));
281 
282  TGeoVolume *Sphere04=geom->MakeSphere("Sphere04",Iron,0,15,0,45,0);
283  Sphere04->SetFillColor(2);
284  Sphere04->SetLineColor(2);
285  top->AddNodeOverlap(Sphere04,1,new TGeoTranslation(145,-145,600));
286 
287  TGeoVolume *Sphere05=geom->MakeSphere("Sphere05",Iron,0,15,0,45,0);
288  Sphere05->SetFillColor(2);
289  Sphere05->SetLineColor(2);
290  top->AddNodeOverlap(Sphere05,1,new TGeoTranslation(-210,0,600));
291 
292  TGeoVolume *Sphere06=geom->MakeSphere("Sphere06",Iron,0,15,0,45,0);
293  Sphere06->SetFillColor(2);
294  Sphere06->SetLineColor(2);
295  top->AddNodeOverlap(Sphere06,1,new TGeoTranslation(210,0,600));
296 
297  TGeoVolume *Sphere07=geom->MakeSphere("Sphere07",Iron,0,15,0,45,0);
298  Sphere07->SetFillColor(2);
299  Sphere07->SetLineColor(2);
300  top->AddNodeOverlap(Sphere07,1,new TGeoTranslation(-145,145,600));
301 
302  TGeoVolume *Torus3=geom->MakeTorus("Torus3",Iron,190,0,10);
303  Torus3->SetFillColor(18);
304  Torus3->SetLineColor(18);
305  top->AddNodeOverlap(Torus3,1,new TGeoTranslation(0,0,610));
306 
307  TGeoVolume *Sphere1=geom->MakeSphere("Sphere1",Iron,0,20,0,180,0,360);
308  Sphere1->SetFillColor(2);
309  Sphere1->SetLineColor(2);
310  top->AddNodeOverlap(Sphere1,1,new TGeoTranslation(0,0,650));
311 
312  TGeoVolume *Tubs=geom->MakeTubs("Tubs",Iron,0,40,50,0,360);
313  Tubs->SetFillColor(29);
314  Tubs->SetLineColor(29);
315  top->AddNodeOverlap(Tubs,1,new TGeoTranslation(0,0,500));
316 
317  TGeoVolume *Tubs1=geom->MakeTubs("Tubs1",Iron,50,60,230,40,150);
318  Tubs1->SetFillColor(18);
319  Tubs1->SetLineColor(18);
320  top->AddNodeOverlap(Tubs1,1,new TGeoTranslation(-170,-30,-400));
321 
322  TGeoVolume *Tubs11=geom->MakeTubs("Tubs11",Iron,50,60,230,220,330);
323  Tubs11->SetFillColor(18);
324  Tubs11->SetLineColor(18);
325  top->AddNodeOverlap(Tubs11,1,new TGeoTranslation(-260,35,-400));
326 
327  TGeoVolume *Sphere111=geom->MakeSphere("Sphere111",Iron,0,10,0,180,0,360);
328  Sphere111->SetFillColor(2);
329  Sphere111->SetLineColor(2);
330  top->AddNodeOverlap(Sphere111,1,new TGeoTranslation(-310,0,-165));
331 
332  TGeoVolume *Sphere112=geom->MakeSphere("Sphere112",Iron,0,10,0,180,0,360);
333  Sphere112->SetFillColor(2);
334  Sphere112->SetLineColor(2);
335  top->AddNodeOverlap(Sphere112,1,new TGeoTranslation(-310,0,-400));
336 
337  TGeoVolume *Sphere113=geom->MakeSphere("Sphere113",Iron,0,10,0,180,0,360);
338  Sphere113->SetFillColor(2);
339  Sphere113->SetLineColor(2);
340  top->AddNodeOverlap(Sphere113,1,new TGeoTranslation(-310,0,-635));
341 
342  TGeoVolume *Tubs2=geom->MakeTubs("Tubs2",Iron,50,60,230,220,330);
343  Tubs2->SetFillColor(18);
344  Tubs2->SetLineColor(18);
345  top->AddNodeOverlap(Tubs2,1,new TGeoTranslation(170,30,-400));
346 
347  TGeoVolume *Tubs21=geom->MakeTubs("Tubs21",Iron,50,60,230,400,510);
348  Tubs21->SetFillColor(18);
349  Tubs21->SetLineColor(18);
350  top->AddNodeOverlap(Tubs21,1,new TGeoTranslation(265,-25,-400));
351 
352  TGeoVolume *Sphere211=geom->MakeSphere("Sphere211",Iron,0,10,0,180,0,360);
353  Sphere211->SetFillColor(2);
354  Sphere211->SetLineColor(2);
355  top->AddNodeOverlap(Sphere211,1,new TGeoTranslation(310,0,-165));
356 
357  TGeoVolume *Sphere212=geom->MakeSphere("Sphere212",Iron,0,10,0,180,0,360);
358  Sphere212->SetFillColor(2);
359  Sphere212->SetLineColor(2);
360  top->AddNodeOverlap(Sphere212,1,new TGeoTranslation(310,0,-400));
361 
362  TGeoVolume *Sphere213=geom->MakeSphere("Sphere213",Iron,0,10,0,180,0,360);
363  Sphere213->SetFillColor(2);
364  Sphere213->SetLineColor(2);
365  top->AddNodeOverlap(Sphere213,1,new TGeoTranslation(310,0,-635));
366 
367  TGeoVolume *Tubs3=geom->MakeTubs("Tubs3",Iron,50,60,230,130,260);
368  Tubs3->SetFillColor(18);
369  Tubs3->SetLineColor(18);
370  top->AddNodeOverlap(Tubs3,1,new TGeoTranslation(30,-170,-400));
371 
372  TGeoVolume *Tubs31=geom->MakeTubs("Tubs31",Iron,50,60,230,310,440);
373  Tubs31->SetFillColor(18);
374  Tubs31->SetLineColor(18);
375  top->AddNodeOverlap(Tubs31,1,new TGeoTranslation(0,-275,-400));
376 
377  TGeoVolume *Sphere311=geom->MakeSphere("Sphere311",Iron,0,10,0,180,0,360);
378  Sphere311->SetFillColor(2);
379  Sphere311->SetLineColor(2);
380  top->AddNodeOverlap(Sphere311,1,new TGeoTranslation(-35,320,-165));
381 
382  TGeoVolume *Sphere312=geom->MakeSphere("Sphere312",Iron,0,10,0,180,0,360);
383  Sphere312->SetFillColor(2);
384  Sphere312->SetLineColor(2);
385  top->AddNodeOverlap(Sphere312,1,new TGeoTranslation(-35,320,-400));
386 
387  TGeoVolume *Sphere313=geom->MakeSphere("Sphere313",Iron,0,10,0,180,0,360);
388  Sphere313->SetFillColor(2);
389  Sphere313->SetLineColor(2);
390  top->AddNodeOverlap(Sphere313,1,new TGeoTranslation(-35,320,-635));
391 
392  TGeoVolume *Tubs4=geom->MakeTubs("Tubs4",Iron,50,60,230,310,440);
393  Tubs4->SetFillColor(18);
394  Tubs4->SetLineColor(18);
395  top->AddNodeOverlap(Tubs4,1,new TGeoTranslation(-30,170,-400));
396 
397  TGeoVolume *Tubs41=geom->MakeTubs("Tubs41",Iron,50,60,230,490,620);
398  Tubs41->SetFillColor(18);
399  Tubs41->SetLineColor(18);
400  top->AddNodeOverlap(Tubs41,1,new TGeoTranslation(0,275,-400));
401 
402  TGeoVolume *Sphere411=geom->MakeSphere("Sphere411",Iron,0,10,0,180,0,360);
403  Sphere411->SetFillColor(2);
404  Sphere411->SetLineColor(2);
405  top->AddNodeOverlap(Sphere411,1,new TGeoTranslation(30,-320,-165));
406 
407  TGeoVolume *Sphere412=geom->MakeSphere("Sphere412",Iron,0,10,0,180,0,360);
408  Sphere412->SetFillColor(2);
409  Sphere412->SetLineColor(2);
410  top->AddNodeOverlap(Sphere412,1,new TGeoTranslation(30,-320,-400));
411 
412  TGeoVolume *Sphere413=geom->MakeSphere("Sphere413",Iron,0,10,0,180,0,360);
413  Sphere413->SetFillColor(2);
414  Sphere413->SetLineColor(2);
415  top->AddNodeOverlap(Sphere413,1,new TGeoTranslation(30,-320,-635));
416 
417  TGeoVolume *Cone010=geom->MakeCone("Cone010",Iron,30,0,30,0,30);
418  Cone010->SetFillColor(2);
419  Cone010->SetLineColor(2);
420  top->AddNodeOverlap(Cone010,1,new TGeoTranslation(0,0,250));
421 
422  TGeoVolume *Torus010=geom->MakeTorus("Torus010",Iron,300,50,40);
423  Torus010->SetFillColor(33);
424  Torus010->SetLineColor(33);
425  top->AddNodeOverlap(Torus010,1,new TGeoTranslation(0,0,250));
426 
427  TGeoVolume *Torus011=geom->MakeTorus("Torus011",Iron,400,10,10);
428  Torus011->SetFillColor(33);
429  Torus011->SetLineColor(33);
430  top->AddNodeOverlap(Torus011,1,new TGeoTranslation(0,0,250));
431 
432  TGeoVolume *Torus012=geom->MakeTorus("Torus012",Iron,200,10,10);
433  Torus012->SetFillColor(33);
434  Torus012->SetLineColor(33);
435  top->AddNodeOverlap(Torus012,1,new TGeoTranslation(0,0,250));
436 
437  TGeoVolume *Sphere010=geom->MakeSphere("Sphere010",Iron,0,10,0,180,0,360);
438  Sphere010->SetFillColor(2);
439  Sphere010->SetLineColor(2);
440  top->AddNodeOverlap(Sphere010,1,new TGeoTranslation(-290,-290,250));
441 
442  TGeoVolume *Sphere011=geom->MakeSphere("Sphere011",Iron,0,10,0,180,0,360);
443  Sphere011->SetFillColor(2);
444  Sphere011->SetLineColor(2);
445  top->AddNodeOverlap(Sphere011,1,new TGeoTranslation(290,290,250));
446 
447  TGeoVolume *Sphere012=geom->MakeSphere("Sphere012",Iron,0,10,0,180,0,360);
448  Sphere012->SetFillColor(2);
449  Sphere012->SetLineColor(2);
450  top->AddNodeOverlap(Sphere012,1,new TGeoTranslation(0,-410,250));
451 
452  TGeoVolume *Sphere013=geom->MakeSphere("Sphere013",Iron,0,10,0,180,0,360);
453  Sphere013->SetFillColor(2);
454  Sphere013->SetLineColor(2);
455  top->AddNodeOverlap(Sphere013,1,new TGeoTranslation(0,410,250));
456 
457  TGeoVolume *Sphere014=geom->MakeSphere("Sphere014",Iron,0,10,0,180,0,360);
458  Sphere014->SetFillColor(2);
459  Sphere014->SetLineColor(2);
460  top->AddNodeOverlap(Sphere014,1,new TGeoTranslation(290,-290,250));
461 
462  TGeoVolume *Sphere015=geom->MakeSphere("Sphere015",Iron,0,10,0,180,0,360);
463  Sphere015->SetFillColor(2);
464  Sphere015->SetLineColor(2);
465  top->AddNodeOverlap(Sphere015,1,new TGeoTranslation(-290,290,250));
466 
467  TGeoVolume *Sphere016=geom->MakeSphere("Sphere016",Iron,0,10,0,180,0,360);
468  Sphere016->SetFillColor(2);
469  Sphere016->SetLineColor(2);
470  top->AddNodeOverlap(Sphere016,1,new TGeoTranslation(410,0,250));
471 
472  TGeoVolume *Sphere017=geom->MakeSphere("Sphere017",Iron,0,10,0,180,0,360);
473  Sphere017->SetFillColor(2);
474  Sphere017->SetLineColor(2);
475  top->AddNodeOverlap(Sphere017,1,new TGeoTranslation(-410,0,250));
476 
477  TGeoVolume *Box010=geom->MakeBox("Box010",Copper,10,10,3);
478  Box010->SetFillColor(3);
479  Box010->SetLineColor(3);
480  top->AddNodeOverlap(Box1,1,new TGeoCombiTrans(-120,-280,300,new TGeoRotation("Box010",0,0,70)));
481 
482  TGeoVolume *Box011=geom->MakeBox("Box011",Copper,10,10,3);
483  Box011->SetFillColor(3);
484  Box011->SetLineColor(3);
485  top->AddNodeOverlap(Box011,1,new TGeoCombiTrans(120,280,300,new TGeoRotation("Box011",0,0,70)));
486 
487  TGeoVolume *Box012=geom->MakeBox("Box012",Copper,10,10,3);
488  Box012->SetFillColor(3);
489  Box012->SetLineColor(3);
490  top->AddNodeOverlap(Box012,1,new TGeoCombiTrans(120,-280,300,new TGeoRotation("Box012",0,0,30)));
491 
492  TGeoVolume *Box013=geom->MakeBox("Box013",Copper,10,10,3);
493  Box013->SetFillColor(3);
494  Box013->SetLineColor(3);
495  top->AddNodeOverlap(Box013,1,new TGeoCombiTrans(-120,280,300,new TGeoRotation("Box013",0,0,30)));
496 
497  TGeoVolume *Box014=geom->MakeBox("Box010",Copper,10,10,3);
498  Box014->SetFillColor(3);
499  Box014->SetLineColor(3);
500  top->AddNodeOverlap(Box014,1,new TGeoCombiTrans(270,-120,300,new TGeoRotation("Box014",0,0,70)));
501 
502  TGeoVolume *Box015=geom->MakeBox("Box015",Copper,10,10,3);
503  Box015->SetFillColor(3);
504  Box015->SetLineColor(3);
505  top->AddNodeOverlap(Box015,1,new TGeoCombiTrans(-270,120,300,new TGeoRotation("Box015",0,0,70)));
506 
507  TGeoVolume *Box016=geom->MakeBox("Box016",Copper,10,10,3);
508  Box016->SetFillColor(3);
509  Box016->SetLineColor(3);
510  top->AddNodeOverlap(Box016,1,new TGeoCombiTrans(270,100,300,new TGeoRotation("Box016",0,0,30)));
511 
512  TGeoVolume *Box017=geom->MakeBox("Box017",Copper,10,10,3);
513  Box017->SetFillColor(3);
514  Box017->SetLineColor(3);
515  top->AddNodeOverlap(Box017,1,new TGeoCombiTrans(-270,-120,300,new TGeoRotation("Box017",0,0,30)));
516 
517 for(int i=1;i<=8;i++){
518  TGeoVolume *Torus0101=geom->MakeTorus("Torus0101",Iron,300,70,40,45*i-4,8);
519  Torus0101->SetFillColor(18);
520  Torus0101->SetLineColor(18);
521  top->AddNodeOverlap(Torus0101,1,new TGeoTranslation(0,0,250));
522 
523  TGeoVolume *Tubs0101=geom->MakeTubs("Line",Iron,0,400,5,45*i-1,45*i+1);
524  Tubs0101->SetFillColor(18);
525  Tubs0101->SetLineColor(18);
526  top->AddNodeOverlap(Tubs0101,1,new TGeoTranslation(0,0,250));
527 }
528 
529  Cone31->SetFillColor(38);
530  top->SetVisibility(0);
531  geom->CloseGeometry();
532 
533  top->Draw("ogl");
534 
535 }