Logo ROOT   6.18/05
Reference Guide
south_gate.C
Go to the documentation of this file.
1/// \file
2/// \ingroup tutorial_geom
3/// Drawing a famous Korean gate, the South gate, called Namdeamoon in Korean, 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 south_gate.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 Lan Hee Yang(yangd5d5@hotmail.com), Dept. of Physics, Univ. of Seoul
16
17#include "TGeoManager.h"
18
19void south_gate()
20{
21
22
23 TGeoManager *geom=new TGeoManager("geom","My first 3D geometry");
24
25
26 TGeoMaterial *vacuum=new TGeoMaterial("vacuum",0,0,0);//a,z,rho
27 TGeoMaterial *Fe=new TGeoMaterial("Fe",55.845,26,7.87);
28
29 //Create media
30
31 TGeoMedium *Air = new TGeoMedium("Vacuum",0,vacuum);
32 TGeoMedium *Iron = new TGeoMedium("Iron",1,Fe);
33
34 //Create volume
35
36 TGeoVolume *top = geom->MakeBox("top",Air,1000,1000,1000);
37 geom->SetTopVolume(top);
38 geom->SetTopVisible(0);
39 // If you want to see the boundary, please input the number, 1 instead of 0.
40 // Like this, geom->SetTopVisible(1);
41
42
43//base
44
45char nBlocks[100];
46int i=1;
47int N = 0;
48int f=0;
49int di[2]; di[0] = 0; di[1] = 30;
50TGeoVolume *mBlock;
51
52while (f<11){
53while (i<14){
54 if (i==6 && f<8){
55 i = i+3;
56 }
57
58 sprintf(nBlocks,"f%d_bg%d",f,N++);
59 mBlock = geom->MakeBox(nBlocks, Iron, 29,149,9);
60 mBlock->SetLineColor(20);
61 if (f<8){
62 if (i<=5 && f<8){
63 top->AddNodeOverlap(mBlock,1,new TGeoTranslation(-120-((i-1)*60)-di[f%2],5,5+(20*f)));
64 } else if (i>5 && f<8){
65 top->AddNodeOverlap(mBlock,1,new TGeoTranslation(120+((i-9)*60) +di[f%2],5,5+(20*f)));
66 }
67 } else {
68 top->AddNodeOverlap(mBlock,1,new TGeoTranslation(-420+(i*60)-di[f%2],5,5+(20*f)));
69 }
70 i++;
71 if (i>=14 && f>=8 && f%2 == 1){
72 sprintf(nBlocks,"f%d_bg%d",f,N++);
73 mBlock = geom->MakeBox(nBlocks, Iron, 29,149,9);
74 mBlock->SetLineColor(20);
75 top->AddNodeOverlap(mBlock,1,new TGeoTranslation(-420+(i*60)-di[f%2],5,5+(20*f)));
76 i++;
77 }
78 if (f%2 ==0){
79 sprintf(nBlocks,"f%d_bg%d",f,N++);
80 mBlock = geom->MakeBox(nBlocks, Iron, 14.5,149,9);
81 mBlock->SetLineColor(20);
82 top->AddNodeOverlap(mBlock,1,new TGeoTranslation(-405,5,5+(20*f)));
83 sprintf(nBlocks,"f%d_bg%d",f,N++);
84 mBlock = geom->MakeBox(nBlocks, Iron, 14.5,149,9);
85 mBlock->SetLineColor(20);
86 top->AddNodeOverlap(mBlock,1,new TGeoTranslation(405,5,5+(20*f)));
87 } else if (f<5){
88 sprintf(nBlocks,"f%d_bg%d",f,N++);
89 mBlock = geom->MakeBox(nBlocks, Iron, 14.5,149,9);
90 mBlock->SetLineColor(20);
91 top->AddNodeOverlap(mBlock,1,new TGeoTranslation(-105,5,5+(20*f)));
92 sprintf(nBlocks,"f%d_bg%d",f,N++);
93 mBlock = geom->MakeBox(nBlocks, Iron, 14.5,149,9);
94 mBlock->SetLineColor(20);
95 top->AddNodeOverlap(mBlock,1,new TGeoTranslation(105,5,5+(20*f)));
96
97 }
98}
99 sprintf(nBlocks,"f%d_bg%d",8,N++);
100 mBlock = geom->MakeBox(nBlocks, Iron, 40,149,9);
101 mBlock->SetLineColor(20);
102 top->AddNodeOverlap(mBlock,1,new TGeoTranslation(-80,5,145));
103 sprintf(nBlocks,"f%d_bg%d",8,N++);
104 mBlock = geom->MakeBox(nBlocks, Iron, 40,149,9);
105 mBlock->SetLineColor(20);
106 top->AddNodeOverlap(mBlock,1,new TGeoTranslation(80,5,145));
107
108 sprintf(nBlocks,"f%d_bg%d",7,N++);
109 mBlock = geom->MakeBox(nBlocks, Iron, 15,149,9);
110 mBlock->SetLineColor(20);
111 top->AddNodeOverlap(mBlock,1,new TGeoTranslation(-75,5,125));
112 sprintf(nBlocks,"f%d_bg%d",7,N++);
113 mBlock = geom->MakeBox(nBlocks, Iron, 15,149,9);
114 mBlock->SetLineColor(20);
115 top->AddNodeOverlap(mBlock,1,new TGeoTranslation(75,5,125));
116
117 sprintf(nBlocks,"f%d_bg%d",6,N++);
118 mBlock = geom->MakeBox(nBlocks, Iron, 24,149,9);
119 mBlock->SetLineColor(20);
120 top->AddNodeOverlap(mBlock,1,new TGeoTranslation(-95,5,105));
121 sprintf(nBlocks,"f%d_bg%d",6,N++);
122 mBlock = geom->MakeBox(nBlocks, Iron, 24,149,9);
123 mBlock->SetLineColor(20);
124 top->AddNodeOverlap(mBlock,1,new TGeoTranslation(95,5,105));
125
126
127
128i=1;f++;
129
130}
131
132
133
134
135//wall
136
137f=0;
138while (f<5){
139i=0;
140while (i<65){
141 sprintf(nBlocks,"f%d_bg%d",f,N++);
142 mBlock = geom->MakeBox(nBlocks, Iron, 5.8,3,3.8);
143 mBlock->SetLineColor(25);
144 top->AddNodeOverlap(mBlock,1,new TGeoTranslation(-384+(i*12),137,218+(f*8)));
145i++;
146
147}
148f++;
149}
150
151
152
153f=0;
154while (f<5){
155i=0;
156while (i<65){
157 sprintf(nBlocks,"f%d_bg%d",f,N++);
158 mBlock = geom->MakeBox(nBlocks, Iron, 5.8,3,3.8);
159 mBlock->SetLineColor(25);
160 top->AddNodeOverlap(mBlock,1,new TGeoTranslation(-384+(i*12),-137,218+(f*8)));
161i++;
162
163}
164f++;
165}
166
167
168
169
170
171f=0;
172while (f<7){
173i=0;
174while (i<22){
175 sprintf(nBlocks,"f%d_bg%d",f,N++);
176 mBlock = geom->MakeBox(nBlocks, Iron, 3,5.8,3.8);
177 mBlock->SetLineColor(25);
178 top->AddNodeOverlap(mBlock,1,new TGeoTranslation(-384,-126+(i*12),218+(f*8)));
179i++;
180
181}
182f++;
183}
184
185
186
187f=0;
188while (f<7){
189i=0;
190while (i<22){
191 sprintf(nBlocks,"f%d_bg%d",f,N++);
192 mBlock = geom->MakeBox(nBlocks, Iron, 3,5.8,3.8);
193 mBlock->SetLineColor(25);
194 top->AddNodeOverlap(mBlock,1,new TGeoTranslation(384,-126+(i*12),218+(f*8)));
195i++;
196
197}
198f++;
199}
200
201
202// arch
203
204
205int k;
206k=0; i=0;
207
208while (i<5){
209while(k<10){
210 sprintf(nBlocks,"ab%d",N++);
211 mBlock = geom->MakeTubs(nBlocks,Iron, 70,89,14, (i*36)+0.5, (i+1)*36-0.5);
212 mBlock->SetLineColor(20);
213 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(0,-130+(k*30),70, new TGeoRotation("r1",0,90,0)));
214 k++;
215}
216 i++; k=0;
217}
218
219 sprintf(nBlocks,"ab%d",N++);
220 mBlock = geom->MakeBox(nBlocks, Iron, 9,149,17);
221 mBlock->SetLineColor(20);
222 top->AddNodeOverlap(mBlock,1,new TGeoTranslation(80,5,14));
223 sprintf(nBlocks,"ab%d",N++);
224 mBlock = geom->MakeBox(nBlocks, Iron, 9,149,18);
225 mBlock->SetLineColor(20);
226 top->AddNodeOverlap(mBlock,1,new TGeoTranslation(80,5,51));
227
228 sprintf(nBlocks,"ab%d",N++);
229 mBlock = geom->MakeBox(nBlocks, Iron, 9,149,17);
230 mBlock->SetLineColor(20);
231 top->AddNodeOverlap(mBlock,1,new TGeoTranslation(-80,5,14));
232 sprintf(nBlocks,"ab%d",N++);
233 mBlock = geom->MakeBox(nBlocks, Iron, 9,149,18);
234 mBlock->SetLineColor(20);
235 top->AddNodeOverlap(mBlock,1,new TGeoTranslation(-80,5,51));
236
237
238
239
240
241//wall's kiwa
242
243k=0; i=0;
244
245while (i<5){
246while(k<52){
247 sprintf(nBlocks,"ab%d",N++);
248 mBlock = geom->MakeTubs(nBlocks,Iron, 1,3,7, 0, 180);
249 mBlock->SetLineColor(12);
250 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-382+(k*15),137,255, new TGeoRotation("r1",90,90,0)));
251 k++;
252}
253 i++; k=0;
254}
255
256
257
258
259
260k=0; i=0;
261
262while (i<5){
263while(k<52){
264 sprintf(nBlocks,"ab%d",N++);
265 mBlock = geom->MakeTubs(nBlocks,Iron, 2.5,3,7, 0, 180);
266 mBlock->SetLineColor(12);
267 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-382+(k*15),-137,255, new TGeoRotation("r1",90,90,0)));
268 k++;
269}
270 i++; k=0;
271}
272
273
274
275k=0; i=0;
276
277while (i<5){
278while(k<20){
279 sprintf(nBlocks,"ab%d",N++);
280 mBlock = geom->MakeTubs(nBlocks,Iron, 2.5,3,6, 0, 180);
281 mBlock->SetLineColor(12);
282 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-382,-123+(k*13),271, new TGeoRotation("r1",0,90,0)));
283 k++;
284}
285 i++; k=0;
286}
287
288
289
290
291k=0; i=0;
292
293while (i<5){
294while(k<20){
295 sprintf(nBlocks,"ab%d",N++);
296 mBlock = geom->MakeTubs(nBlocks,Iron, 2.5,3,7, 0, 180);
297 mBlock->SetLineColor(12);
298 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(382,-123+(k*13),271, new TGeoRotation("r1",0,90,0)));
299 k++;
300}
301 i++; k=0;
302}
303
304
305////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
306
307
308// 1 floor
309
310
311k=0; i=0;
312
313while (i<5){
314while(k<7){
315 sprintf(nBlocks,"ab%d",N++);
316 mBlock = geom->MakeTubs(nBlocks,Iron, 0,5,56, 0, 360);
317 mBlock->SetLineColor(50);
318
319 if (k<=2){
320
321 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-300+(k*100),80,260, new TGeoRotation("r1",0,0,0)));
322 }else if (k>=4){
323 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-300+(k*100),80,260, new TGeoRotation("r1",0,0,0)));
324 }
325
326
327k++;
328}
329 i++; k=0;
330}
331
332
333
334k=0; i=0;
335
336while (i<5){
337while(k<7){
338 sprintf(nBlocks,"ab%d",N++);
339 mBlock = geom->MakeTubs(nBlocks,Iron, 0,5,56, 0, 360);
340 mBlock->SetLineColor(50);
341
342 if (k<=2){
343
344 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-300+(k*100),-80,260, new TGeoRotation("r1",0,0,0)));
345 }else if (k>=4){
346 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-300+(k*100),-80,260, new TGeoRotation("r1",0,0,0)));
347 }
348
349
350k++;
351}
352 i++; k=0;
353}
354
355// ||=====||======||=====||=====||=====||=====||
356
357
358 sprintf(nBlocks,"ab%d",N++);
359 mBlock = geom->MakeBox(nBlocks, Iron, 298,78,8);
360 mBlock->SetLineColor(42);
361 top->AddNodeOverlap(mBlock,1,new TGeoTranslation(0,0,300));
362
363
364
365 sprintf(nBlocks,"ab%d",N++);
366 mBlock = geom->MakeBox(nBlocks, Iron, 298,78,5);
367 mBlock->SetLineColor(42);
368 top->AddNodeOverlap(mBlock,1,new TGeoTranslation(0,0,320));
369
370
371
372//1
373k=0; i=0;
374
375while (i<5){
376while(k<6){
377 sprintf(nBlocks,"ab%d",N++);
378 mBlock = geom->MakeBox(nBlocks,Iron,18,10,8);
379 mBlock->SetLineColor(8);
380 {
381 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-250+(k*100),70,300, new TGeoRotation("r1",0,0,0)));
382
383 }
384 k++;
385}
386 i++; k=0;
387}
388
389
390
391
392
393k=0; i=0;
394
395while (i<5){
396while(k<6){
397 sprintf(nBlocks,"ab%d",N++);
398 mBlock = geom->MakeBox(nBlocks,Iron,18,10,8);
399 mBlock->SetLineColor(8);
400 {
401 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-250+(k*100),-70,300, new TGeoRotation("r1",0,0,0)));
402
403 }
404 k++;
405}
406 i++; k=0;
407}
408
409
410
411
412 sprintf(nBlocks,"ab%d",N++);
413 mBlock = geom->MakeBox(nBlocks, Iron, 18,10,8);
414 mBlock->SetLineColor(8);
415 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-290,0,300, new TGeoRotation("r1",90,0,0)));
416
417
418
419 sprintf(nBlocks,"ab%d",N++);
420 mBlock = geom->MakeBox(nBlocks, Iron, 18,10,8);
421 mBlock->SetLineColor(8);
422 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(290,0,300, new TGeoRotation("r1",90,0,0)));
423
424
425
426
427
428
429//2
430k=0; i=0;
431
432while (i<5){
433while(k<6){
434 sprintf(nBlocks,"ab%d",N++);
435 mBlock = geom->MakeBox(nBlocks,Iron,18,10,5);
436 mBlock->SetLineColor(8);
437 {
438 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-250+(k*100),70,320, new TGeoRotation("r1",0,0,0)));
439
440 }
441 k++;
442}
443 i++; k=0;
444}
445
446
447
448
449
450k=0; i=0;
451
452while (i<5){
453 while(k<6){
454 sprintf(nBlocks,"ab%d",N++);
455 mBlock = geom->MakeBox(nBlocks,Iron,18,10,5);
456 mBlock->SetLineColor(8);
457 {
458 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-250+(k*100),-70,320, new TGeoRotation("r1",0,0,0)));
459
460 }
461 k++;
462 }
463 i++; k=0;
464}
465
466
467
468 sprintf(nBlocks,"ab%d",N++);
469 mBlock = geom->MakeBox(nBlocks, Iron, 18,10,5);
470 mBlock->SetLineColor(8);
471 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-290,0,320, new TGeoRotation("r1",90,0,0)));
472
473
474
475 sprintf(nBlocks,"ab%d",N++);
476 mBlock = geom->MakeBox(nBlocks, Iron, 18,10,5);
477 mBlock->SetLineColor(8);
478 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(290,0,320, new TGeoRotation("r1",90,0,0)));
479
480
481
482
483
484
485
486
487
488
489
490//___||____||_____||____||____||____||____||
491
492
493k=0; i=0;
494
495while (i<5){
496 while(k<19){
497 sprintf(nBlocks,"ab%d",N++);
498 mBlock = geom->MakeBox(nBlocks,Iron, 5,5,20);
499 mBlock->SetLineColor(50);
500 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-300+(k*33.3),78,345, new TGeoRotation("r1",0,0,0)));
501 k++;
502 }
503 i++; k=0;
504}
505
506
507
508
509k=0; i=0;
510
511while (i<5){
512 while(k<19){
513 sprintf(nBlocks,"ab%d",N++);
514 mBlock = geom->MakeBox(nBlocks,Iron, 5,5,20);
515 mBlock->SetLineColor(50);
516 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-300+(k*33.3),-78,345, new TGeoRotation("r1",0,0,0)));
517 k++;
518 }
519 i++; k=0;
520}
521
522
523
524k=0; i=0;
525
526while (i<5){
527 while(k<5){
528 sprintf(nBlocks,"ab%d",N++);
529 mBlock = geom->MakeBox(nBlocks,Iron, 5,5,20);
530 mBlock->SetLineColor(50);
531 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-300,-78+(k*33),345, new TGeoRotation("r1",0,0,0)));
532 k++;
533 }
534 i++; k=0;
535}
536
537
538
539
540
541k=0; i=0;
542
543while (i<5){
544 while(k<5){
545 sprintf(nBlocks,"ab%d",N++);
546 mBlock = geom->MakeBox(nBlocks,Iron, 5,5,20);
547 mBlock->SetLineColor(50);
548 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(300,-78+(k*33),345, new TGeoRotation("r1",0,0,0)));
549 k++;
550 }
551 i++; k=0;
552}
553
554// ||// ||// ||// ||//
555
556
557
558k=0; i=0;
559
560while (i<5){
561 while(k<19){
562 sprintf(nBlocks,"ab%d",N++);
563 mBlock = geom->MakeBox(nBlocks,Iron, 5,5,20);
564 mBlock->SetLineColor(50);
565 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-300+(k*33.3),90,342, new TGeoRotation("r1",0,-45,0)));
566 k++;
567 }
568 i++; k=0;
569}
570
571
572
573
574k=0; i=0;
575
576while (i<5){
577 while(k<19){
578 sprintf(nBlocks,"ab%d",N++);
579 mBlock = geom->MakeBox(nBlocks,Iron, 5,5,20);
580 mBlock->SetLineColor(50);
581 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-300+(k*33.3),-90,342, new TGeoRotation("r1",0,45,0)));
582 k++;
583 }
584 i++; k=0;
585}
586
587
588
589k=0; i=0;
590
591while (i<5){
592 while(k<5){
593 sprintf(nBlocks,"ab%d",N++);
594 mBlock = geom->MakeBox(nBlocks,Iron, 5,5,20);
595 mBlock->SetLineColor(50);
596 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-318,-78+(k*33),345, new TGeoRotation("r1",-90,45,0)));
597 k++;
598 }
599 i++; k=0;
600}
601
602
603
604
605
606k=0; i=0;
607
608while (i<5){
609 while(k<5){
610 sprintf(nBlocks,"ab%d",N++);
611 mBlock = geom->MakeBox(nBlocks,Iron, 5,5,20);
612 mBlock->SetLineColor(50);
613 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(318,-78+(k*33),345, new TGeoRotation("r1",90,45,0)));
614 k++;
615 }
616 i++; k=0;
617}
618
619
620// /// || / / / / / / / || / / / / / / / / || / / / / / / / / / / /
621
622
623
624 sprintf(nBlocks,"ab%d",N++);
625 mBlock = geom->MakeBox(nBlocks, Iron, 330,10,2);
626 mBlock->SetLineColor(42);
627 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(0,-107,362, new TGeoRotation("r1",0,-45,0)));
628
629
630
631
632 sprintf(nBlocks,"ab%d",N++);
633 mBlock = geom->MakeBox(nBlocks, Iron, 330,10,2);
634 mBlock->SetLineColor(42);
635 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(0,107,362, new TGeoRotation("r1",0,45,0)));
636
637
638
639 sprintf(nBlocks,"ab%d",N++);
640 mBlock = geom->MakeBox(nBlocks, Iron, 110,10,2);
641 mBlock->SetLineColor(42);
642 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(330,0,362, new TGeoRotation("r1",90,-45,0)));
643
644
645
646 sprintf(nBlocks,"ab%d",N++);
647 mBlock = geom->MakeBox(nBlocks, Iron, 110,10,2);
648 mBlock->SetLineColor(42);
649 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-330,0,362, new TGeoRotation("r1",90,45,0)));
650
651
652
653
654/////////////////////// add box
655
656
657
658
659k=0; i=0;
660
661while (i<5){
662 while(k<6){
663 sprintf(nBlocks,"ab%d",N++);
664 mBlock = geom->MakeBox(nBlocks,Iron,18,10,2);
665 mBlock->SetLineColor(8);
666 {
667 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-270+(k*100),-108,362, new TGeoRotation("r1",0,-45,0)));
668
669 }
670 k++;
671 }
672 i++; k=0;
673}
674
675
676
677
678k=0; i=0;
679
680while (i<5){
681 while(k<6){
682 sprintf(nBlocks,"ab%d",N++);
683 mBlock = geom->MakeBox(nBlocks,Iron,18,10,2);
684 mBlock->SetLineColor(8);
685 {
686 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-270+(k*100),108,362, new TGeoRotation("r1",0,45,0)));
687
688 }
689 k++;
690 }
691 i++; k=0;
692}
693
694
695 sprintf(nBlocks,"ab%d",N++);
696 mBlock = geom->MakeBox(nBlocks, Iron, 18,10,2);
697 mBlock->SetLineColor(8);
698 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(331,0,362, new TGeoRotation("r1",90,-45,0)));
699
700
701
702 sprintf(nBlocks,"ab%d",N++);
703 mBlock = geom->MakeBox(nBlocks, Iron, 18,10,2);
704 mBlock->SetLineColor(8);
705 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-331,0,362, new TGeoRotation("r1",90,45,0)));
706
707
708
709
710
711/////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
712
713// 2nd floor
714
715
716k=0; i=0;
717
718while (i<5){
719 while(k<7){
720 sprintf(nBlocks,"ab%d",N++);
721 mBlock = geom->MakeTubs(nBlocks,Iron, 0,5,30, 0, 360);
722 mBlock->SetLineColor(50);
723
724 if (k<=2){
725
726 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-300+(k*100),80,465, new TGeoRotation("r1",0,0,0)));
727 }else if (k>=4){
728 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-300+(k*100),80,465, new TGeoRotation("r1",0,0,0)));
729 }
730
731
732 k++;
733 }
734 i++; k=0;
735}
736
737
738
739k=0; i=0;
740
741while (i<5){
742 while(k<7){
743 sprintf(nBlocks,"ab%d",N++);
744 mBlock = geom->MakeTubs(nBlocks,Iron, 0,5,30, 0, 360);
745 mBlock->SetLineColor(50);
746
747 if (k<=2){
748
749 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-300+(k*100),-80,465, new TGeoRotation("r1",0,0,0)));
750 }else if (k>=4){
751 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-300+(k*100),-80,465, new TGeoRotation("r1",0,0,0)));
752 }
753
754
755 k++;
756 }
757 i++; k=0;
758}
759
760
761
762
763// ||=====||======||=====||=====||=====||=====||
764
765
766 sprintf(nBlocks,"ab%d",N++);
767 mBlock = geom->MakeBox(nBlocks, Iron, 302,80,8);
768 mBlock->SetLineColor(42);
769 top->AddNodeOverlap(mBlock,1,new TGeoTranslation(0,0,480));
770
771
772
773 sprintf(nBlocks,"ab%d",N++);
774 mBlock = geom->MakeBox(nBlocks, Iron, 302,80,5);
775 mBlock->SetLineColor(42);
776 top->AddNodeOverlap(mBlock,1,new TGeoTranslation(0,0,500));
777
778
779 sprintf(nBlocks,"ab%d",N++);
780 mBlock = geom->MakeBox(nBlocks, Iron, 305,80,2.5);
781 mBlock->SetLineColor(50);
782 top->AddNodeOverlap(mBlock,1,new TGeoTranslation(0,0,465));
783
784
785///////////////////////add box
786
787
788
789
790
791
792//1
793k=0; i=0;
794
795while (i<5){
796 while(k<6){
797 sprintf(nBlocks,"ab%d",N++);
798 mBlock = geom->MakeBox(nBlocks,Iron,18,10,8);
799 mBlock->SetLineColor(8);
800 {
801 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-250+(k*100),71,480, new TGeoRotation("r1",0,0,0)));
802
803 }
804 k++;
805 }
806 i++; k=0;
807}
808
809
810
811
812
813k=0; i=0;
814
815while (i<5){
816 while(k<6){
817 sprintf(nBlocks,"ab%d",N++);
818 mBlock = geom->MakeBox(nBlocks,Iron,18,10,8);
819 mBlock->SetLineColor(8);
820 {
821 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-250+(k*100),-71,480, new TGeoRotation("r1",0,0,0)));
822
823 }
824 k++;
825 }
826 i++; k=0;
827}
828
829
830
831
832 sprintf(nBlocks,"ab%d",N++);
833 mBlock = geom->MakeBox(nBlocks, Iron, 18,10,8);
834 mBlock->SetLineColor(8);
835 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-293,0,480, new TGeoRotation("r1",90,0,0)));
836
837
838
839 sprintf(nBlocks,"ab%d",N++);
840 mBlock = geom->MakeBox(nBlocks, Iron, 18,10,8);
841 mBlock->SetLineColor(8);
842 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(293,0,480, new TGeoRotation("r1",90,0,0)));
843
844
845
846
847
848
849//2
850k=0; i=0;
851
852while (i<5){
853 while(k<6){
854 sprintf(nBlocks,"ab%d",N++);
855 mBlock = geom->MakeBox(nBlocks,Iron,18,10,5);
856 mBlock->SetLineColor(8);
857 {
858 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-250+(k*100),71,500, new TGeoRotation("r1",0,0,0)));
859
860 }
861 k++;
862 }
863 i++; k=0;
864}
865
866
867
868
869
870k=0; i=0;
871
872while (i<5){
873 while(k<6){
874 sprintf(nBlocks,"ab%d",N++);
875 mBlock = geom->MakeBox(nBlocks,Iron,18,10,5);
876 mBlock->SetLineColor(8);
877 {
878 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-250+(k*100),-71,500, new TGeoRotation("r1",0,0,0)));
879
880 }
881 k++;
882 }
883 i++; k=0;
884}
885
886
887
888 sprintf(nBlocks,"ab%d",N++);
889 mBlock = geom->MakeBox(nBlocks, Iron, 18,10,5);
890 mBlock->SetLineColor(8);
891 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-293,0,500, new TGeoRotation("r1",90,0,0)));
892
893
894
895 sprintf(nBlocks,"ab%d",N++);
896 mBlock = geom->MakeBox(nBlocks, Iron, 18,10,5);
897 mBlock->SetLineColor(8);
898 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(293,0,500, new TGeoRotation("r1",90,0,0)));
899
900
901
902
903
904
905
906
907
908
909
910// 1 ___||____||_____||____||____||____||____||
911
912
913k=0; i=0;
914
915while (i<5){
916 while(k<25){
917 sprintf(nBlocks,"ab%d",N++);
918 mBlock = geom->MakeBox(nBlocks,Iron, 1.5,5,15);
919 mBlock->SetLineColor(50);
920 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-300+(k*25),78,450, new TGeoRotation("r1",0,0,0)));
921 k++;
922 }
923 i++; k=0;
924}
925
926
927
928
929k=0; i=0;
930
931while (i<5){
932 while(k<25){
933 sprintf(nBlocks,"ab%d",N++);
934 mBlock = geom->MakeBox(nBlocks,Iron, 1.5,5,15);
935 mBlock->SetLineColor(50);
936 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-300+(k*25),-78,450, new TGeoRotation("r1",0,0,0)));
937 k++;
938 }
939 i++; k=0;
940}
941
942
943
944k=0; i=0;
945
946while (i<5){
947 while(k<7){
948 sprintf(nBlocks,"ab%d",N++);
949 mBlock = geom->MakeBox(nBlocks,Iron, 5,1.5,15);
950 mBlock->SetLineColor(50);
951 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-300,-78+(k*25),450, new TGeoRotation("r1",0,0,0)));
952 k++;
953 }
954 i++; k=0;
955}
956
957
958
959
960
961k=0; i=0;
962
963while (i<5){
964 while (k<7){
965 sprintf(nBlocks,"ab%d",N++);
966 mBlock = geom->MakeBox(nBlocks,Iron, 5,1.5,15);
967 mBlock->SetLineColor(50);
968 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(300,-78+(k*25),450, new TGeoRotation("r1",0,0,0)));
969 k++;
970 }
971 i++; k=0;
972}
973
974
975
976
977// 2 ___||____||_____||____||____||____||____||
978
979
980k=0; i=0;
981
982while (i<5){
983while(k<19){
984 sprintf(nBlocks,"ab%d",N++);
985 mBlock = geom->MakeBox(nBlocks,Iron, 5,5,20);
986 mBlock->SetLineColor(50);
987 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-300+(k*33.3),78,525, new TGeoRotation("r1",0,0,0)));
988 k++;
989}
990 i++; k=0;
991}
992
993
994
995
996k=0; i=0;
997
998while (i<5){
999 while(k<19){
1000 sprintf(nBlocks,"ab%d",N++);
1001 mBlock = geom->MakeBox(nBlocks,Iron, 5,5,20);
1002 mBlock->SetLineColor(50);
1003 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-300+(k*33.3),-78,525, new TGeoRotation("r1",0,0,0)));
1004 k++;
1005 }
1006 i++; k=0;
1007}
1008
1009
1010
1011k=0; i=0;
1012
1013while (i<5){
1014 while(k<5){
1015 sprintf(nBlocks,"ab%d",N++);
1016 mBlock = geom->MakeBox(nBlocks,Iron, 5,5,20);
1017 mBlock->SetLineColor(50);
1018 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-300,-78+(k*33),525, new TGeoRotation("r1",0,0,0)));
1019 k++;
1020 }
1021 i++; k=0;
1022}
1023
1024
1025
1026
1027
1028k=0; i=0;
1029
1030while (i<5){
1031 while(k<5){
1032 sprintf(nBlocks,"ab%d",N++);
1033 mBlock = geom->MakeBox(nBlocks,Iron, 5,5,20);
1034 mBlock->SetLineColor(50);
1035 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(300,-78+(k*33),525, new TGeoRotation("r1",0,0,0)));
1036 k++;
1037 }
1038 i++; k=0;
1039}
1040
1041
1042
1043
1044// ||// ||// ||// ||//
1045
1046//down
1047
1048k=0; i=0;
1049
1050while (i<5){
1051 while(k<19){
1052 sprintf(nBlocks,"ab%d",N++);
1053 mBlock = geom->MakeBox(nBlocks,Iron, 5,5,20);
1054 mBlock->SetLineColor(50);
1055 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-300+(k*33.3),90,522, new TGeoRotation("r1",0,-45,0)));
1056 k++;
1057 }
1058 i++; k=0;
1059}
1060
1061
1062
1063
1064k=0; i=0;
1065
1066while (i<5){
1067 while(k<19){
1068 sprintf(nBlocks,"ab%d",N++);
1069 mBlock = geom->MakeBox(nBlocks,Iron, 5,5,20);
1070 mBlock->SetLineColor(50);
1071 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-300+(k*33.3),-90,522, new TGeoRotation("r1",0,45,0)));
1072 k++;
1073 }
1074 i++; k=0;
1075}
1076
1077
1078k=0; i=0;
1079
1080while (i<5){
1081 while(k<5){
1082 sprintf(nBlocks,"ab%d",N++);
1083 mBlock = geom->MakeBox(nBlocks,Iron, 5,5,20);
1084 mBlock->SetLineColor(50);
1085 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-318,-78+(k*33.3),525, new TGeoRotation("r1",-90,45,0)));
1086 k++;
1087 }
1088 i++; k=0;
1089}
1090
1091
1092
1093
1094
1095k=0; i=0;
1096
1097while (i<5){
1098 while(k<5){
1099 sprintf(nBlocks,"ab%d",N++);
1100 mBlock = geom->MakeBox(nBlocks,Iron, 5,5,20);
1101 mBlock->SetLineColor(50);
1102 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(318,-78+(k*33.3),525, new TGeoRotation("r1",90,45,0)));
1103 k++;
1104 }
1105 i++; k=0;
1106}
1107
1108
1109// up
1110
1111
1112k=0; i=0;
1113
1114while (i<5){
1115 while(k<50){
1116 sprintf(nBlocks,"ab%d",N++);
1117 mBlock = geom->MakeBox(nBlocks,Iron, 2.5,2.5,20);
1118 mBlock->SetLineColor(50);
1119 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-320+(k*13),115,562, new TGeoRotation("r1",0,-115,0)));
1120 k++;
1121 }
1122 i++; k=0;
1123}
1124
1125
1126
1127
1128k=0; i=0;
1129
1130while (i<5){
1131 while(k<50){
1132 sprintf(nBlocks,"ab%d",N++);
1133 mBlock = geom->MakeBox(nBlocks,Iron, 2.5,2.5,20);
1134 mBlock->SetLineColor(50);
1135 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-320+(k*13),-115,562, new TGeoRotation("r1",0,115,0)));
1136 k++;
1137 }
1138 i++; k=0;
1139}
1140
1141
1142
1143k=0; i=0;
1144
1145while (i<5){
1146 while(k<17){
1147 sprintf(nBlocks,"ab%d",N++);
1148 mBlock = geom->MakeBox(nBlocks,Iron, 2.5,2.5,20);
1149 mBlock->SetLineColor(50);
1150 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-340,-98+(k*13),565, new TGeoRotation("r1",-90,115,0)));
1151 k++;
1152 }
1153 i++; k=0;
1154}
1155
1156
1157
1158
1159
1160k=0; i=0;
1161
1162while (i<5){
1163 while(k<17){
1164 sprintf(nBlocks,"ab%d",N++);
1165 mBlock = geom->MakeBox(nBlocks,Iron, 2.5,2.5,20);
1166 mBlock->SetLineColor(50);
1167 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(340,-98+(k*13),565, new TGeoRotation("r1",90,115,0)));
1168 k++;
1169 }
1170 i++; k=0;
1171}
1172
1173
1174//up2
1175
1176
1177
1178k=0; i=0;
1179
1180while (i<5){
1181 while(k<50){
1182 sprintf(nBlocks,"ab%d",N++);
1183 mBlock = geom->MakeBox(nBlocks,Iron, 2.5,2.5,20);
1184 mBlock->SetLineColor(50);
1185 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-320+(k*13),115,375, new TGeoRotation("r1",0,-115,0)));
1186 k++;
1187 }
1188 i++; k=0;
1189}
1190
1191
1192
1193
1194k=0; i=0;
1195
1196while (i<5){
1197 while(k<50){
1198 sprintf(nBlocks,"ab%d",N++);
1199 mBlock = geom->MakeBox(nBlocks,Iron, 2.5,2.5,20);
1200 mBlock->SetLineColor(50);
1201 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-320+(k*13),-115,375, new TGeoRotation("r1",0,115,0)));
1202 k++;
1203 }
1204 i++; k=0;
1205}
1206
1207
1208
1209k=0; i=0;
1210
1211while (i<5){
1212 while(k<17){
1213 sprintf(nBlocks,"ab%d",N++);
1214 mBlock = geom->MakeBox(nBlocks,Iron, 2.5,2.5,20);
1215 mBlock->SetLineColor(50);
1216 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-340,-98+(k*13),375, new TGeoRotation("r1",-90,115,0)));
1217 k++;
1218 }
1219 i++; k=0;
1220}
1221
1222
1223
1224
1225
1226k=0; i=0;
1227
1228while (i<5){
1229 while(k<17){
1230 sprintf(nBlocks,"ab%d",N++);
1231 mBlock = geom->MakeBox(nBlocks,Iron, 2.5,2.5,20);
1232 mBlock->SetLineColor(50);
1233 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(340,-98+(k*13),375, new TGeoRotation("r1",90,115,0)));
1234 k++;
1235 }
1236 i++; k=0;
1237}
1238
1239
1240//up 3
1241
1242k=0; i=0;
1243
1244while (i<5){
1245 while(k<50){
1246 sprintf(nBlocks,"ab%d",N++);
1247 mBlock = geom->MakeBox(nBlocks,Iron, 2.5,2.5,20);
1248 mBlock->SetLineColor(44);
1249 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-320+(k*13),115,568, new TGeoRotation("r1",0,-115,0)));
1250 k++;
1251 }
1252 i++; k=0;
1253}
1254
1255
1256
1257
1258k=0; i=0;
1259
1260while (i<5){
1261 while(k<50){
1262 sprintf(nBlocks,"ab%d",N++);
1263 mBlock = geom->MakeBox(nBlocks,Iron, 2.5,2.5,20);
1264 mBlock->SetLineColor(44);
1265 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-320+(k*13),-115,568, new TGeoRotation("r1",0,115,0)));
1266 k++;
1267 }
1268 i++; k=0;
1269}
1270
1271
1272k=0; i=0;
1273
1274while (i<5){
1275 while(k<17){
1276 sprintf(nBlocks,"ab%d",N++);
1277 mBlock = geom->MakeBox(nBlocks,Iron, 2.5,2.5,20);
1278 mBlock->SetLineColor(44);
1279 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-340,-98+(k*13),568, new TGeoRotation("r1",-90,115,0)));
1280 k++;
1281 }
1282 i++; k=0;
1283}
1284
1285
1286
1287
1288
1289k=0; i=0;
1290
1291while (i<5){
1292 while(k<17){
1293 sprintf(nBlocks,"ab%d",N++);
1294 mBlock = geom->MakeBox(nBlocks,Iron, 2.5,2.5,20);
1295 mBlock->SetLineColor(44);
1296 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(340,-98+(k*13),568, new TGeoRotation("r1",90,115,0)));
1297 k++;
1298 }
1299 i++; k=0;
1300}
1301
1302
1303
1304
1305
1306//up4
1307
1308
1309k=0; i=0;
1310
1311while (i<5){
1312 while(k<50){
1313 sprintf(nBlocks,"ab%d",N++);
1314 mBlock = geom->MakeBox(nBlocks,Iron, 2.5,2.5,20);
1315 mBlock->SetLineColor(44);
1316 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-320+(k*13),115,385, new TGeoRotation("r1",0,-115,0)));
1317 k++;
1318 }
1319 i++; k=0;
1320}
1321
1322
1323
1324
1325k=0; i=0;
1326
1327while (i<5){
1328 while(k<50){
1329 sprintf(nBlocks,"ab%d",N++);
1330 mBlock = geom->MakeBox(nBlocks,Iron, 2.5,2.5,20);
1331 mBlock->SetLineColor(44);
1332 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-320+(k*13),-115,385, new TGeoRotation("r1",0,115,0)));
1333 k++;
1334 }
1335 i++; k=0;
1336}
1337
1338
1339k=0; i=0;
1340
1341while (i<5){
1342 while(k<17){
1343 sprintf(nBlocks,"ab%d",N++);
1344 mBlock = geom->MakeBox(nBlocks,Iron, 2.5,2.5,20);
1345 mBlock->SetLineColor(44);
1346 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-340,-98+(k*13),385, new TGeoRotation("r1",-90,115,0)));
1347 k++;
1348 }
1349 i++; k=0;
1350}
1351
1352
1353
1354
1355
1356k=0; i=0;
1357
1358while (i<5){
1359 while(k<17){
1360 sprintf(nBlocks,"ab%d",N++);
1361 mBlock = geom->MakeBox(nBlocks,Iron, 2.5,2.5,20);
1362 mBlock->SetLineColor(44);
1363 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(340,-98+(k*13),385, new TGeoRotation("r1",90,115,0)));
1364 k++;
1365 }
1366 i++; k=0;
1367}
1368
1369
1370// up kiwa
1371 //=========
1372 sprintf(nBlocks,"ab%d",N++);
1373 mBlock = geom->MakeBox(nBlocks,Iron, 270,15,20);
1374 mBlock->SetLineColor(10);
1375 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(0,0,620, new TGeoRotation("r1",0,0,0)));
1376 //===============//2
1377 sprintf(nBlocks,"ab%d",N++);
1378 mBlock = geom->MakeBox(nBlocks,Iron, 75,15,20);
1379 mBlock->SetLineColor(10);
1380 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(300,-50,600, new TGeoRotation("r1",0,20,-40)));
1381
1382 sprintf(nBlocks,"ab%d",N++);
1383 mBlock = geom->MakeBox(nBlocks,Iron, 75,15,20);
1384 mBlock->SetLineColor(10);
1385 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(300,50,600, new TGeoRotation("r1",0,-20,40)));
1386
1387 sprintf(nBlocks,"ab%d",N++);
1388 mBlock = geom->MakeBox(nBlocks,Iron, 75,15,20);
1389 mBlock->SetLineColor(10);
1390 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-300,50,600, new TGeoRotation("r1",0,-20,-40)));
1391
1392 sprintf(nBlocks,"ab%d",N++);
1393 mBlock = geom->MakeBox(nBlocks,Iron, 75,15,20);
1394 mBlock->SetLineColor(10);
1395 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-300,-50,600, new TGeoRotation("r1",0,20,40)));
1396
1397
1398
1399
1400 //===============//1
1401 sprintf(nBlocks,"ab%d",N++);
1402 mBlock = geom->MakeBox(nBlocks,Iron, 50,15,20);
1403 mBlock->SetLineColor(10);
1404 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(300,-80,413, new TGeoRotation("r1",0,20,-40)));
1405
1406 sprintf(nBlocks,"ab%d",N++);
1407 mBlock = geom->MakeBox(nBlocks,Iron, 50,15,20);
1408 mBlock->SetLineColor(10);
1409 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(300,80,413, new TGeoRotation("r1",0,-20,40)));
1410
1411 sprintf(nBlocks,"ab%d",N++);
1412 mBlock = geom->MakeBox(nBlocks,Iron, 50,15,20);
1413 mBlock->SetLineColor(10);
1414 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-300,80,413, new TGeoRotation("r1",0,-20,-40)));
1415
1416 sprintf(nBlocks,"ab%d",N++);
1417 mBlock = geom->MakeBox(nBlocks,Iron, 50,15,20);
1418 mBlock->SetLineColor(10);
1419 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-300,-80,413, new TGeoRotation("r1",0,20,40)));
1420
1421
1422
1423
1424// _1_
1425
1426//front
1427
1428k=0; i=0;
1429while (i<7){
1430 while(k<44){
1431 sprintf(nBlocks,"ab%d",N++);
1432 mBlock = geom->MakeTubs(nBlocks,Iron, 3,6,6,10,170);
1433 mBlock->SetLineColor(13);
1434 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-280+(k*13),70+(i*12.5),425-(i*5), new TGeoRotation("r1",0,60,0)));
1435 k++;
1436 }
1437 i++; k=0;
1438}
1439
1440
1441
1442k=0; i=0;
1443
1444while (i<7){
1445 while(k<44){
1446 sprintf(nBlocks,"ab%d",N++);
1447 mBlock = geom->MakeTubs(nBlocks,Iron, 3,6,6,10,170);
1448 mBlock->SetLineColor(13);
1449 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-280+(k*13),-70-(i*12.5),425-(i*5), new TGeoRotation("r1",0,120,0)));
1450 k++;
1451 }
1452 i++; k=0;
1453}
1454
1455//_2_
1456
1457
1458
1459
1460k=0; i=0;
1461while (i<11){
1462 while(k<43){
1463 sprintf(nBlocks,"ab%d",N++);
1464 mBlock = geom->MakeTubs(nBlocks,Iron, 3,6,6,10,170);
1465 mBlock->SetLineColor(13);
1466
1467 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-270+(k*13),15+(i*12.5),620-(i*5), new TGeoRotation("r1",0,60,0)));
1468 k++;
1469 }
1470 i++; k=0;
1471}
1472
1473
1474
1475k=0; i=0;
1476
1477while (i<11){
1478 while(k<43){
1479 sprintf(nBlocks,"ab%d",N++);
1480 mBlock = geom->MakeTubs(nBlocks,Iron, 3,6,6,10,170);
1481 mBlock->SetLineColor(13);
1482 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-270+(k*13),-15-(i*12.5),620-(i*5), new TGeoRotation("r1",0,120,0)));
1483 k++;
1484 }
1485 i++; k=0;
1486}
1487
1488
1489
1490
1491//////left
1492k=0; i=0;
1493
1494while (i<6){
1495 while(k<11){
1496 sprintf(nBlocks,"ab%d",N++);
1497 mBlock = geom->MakeTubs(nBlocks,Iron, 3,6,6,10,170);
1498 mBlock->SetLineColor(13);
1499 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-335,81.25+(i*12.5),592.5-(i*2), new TGeoRotation("r1",0,60,0)));
1500 k++;
1501 }
1502 i++; k=0;
1503}
1504
1505k=0; i=0;
1506
1507while (i<7){
1508 while(k<11){
1509 sprintf(nBlocks,"ab%d",N++);
1510 mBlock = geom->MakeTubs(nBlocks,Iron, 3,6,6,10,170);
1511 mBlock->SetLineColor(13);
1512 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-322,69.75+(i*12.5),595-(i*2), new TGeoRotation("r1",0,60,0)));
1513 k++;
1514 }
1515 i++; k=0;
1516}
1517
1518
1519k=0; i=0;
1520
1521while (i<8){
1522 while(k<11){
1523 sprintf(nBlocks,"ab%d",N++);
1524 mBlock = geom->MakeTubs(nBlocks,Iron, 3,6,6,10,170);
1525 mBlock->SetLineColor(13);
1526 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-309,56.25+(i*12.5),605-(i*4), new TGeoRotation("r1",0,60,0)));
1527 k++;
1528 }
1529 i++; k=0;
1530}
1531
1532k=0; i=0;
1533
1534while (i<9){
1535 while(k<11){
1536 sprintf(nBlocks,"ab%d",N++);
1537 mBlock = geom->MakeTubs(nBlocks,Iron, 3,6,6,10,170);
1538 mBlock->SetLineColor(13);
1539 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-296,50+(i*12.5),610-(i*4), new TGeoRotation("r1",0,60,0)));
1540 k++;
1541 }
1542 i++; k=0;
1543}
1544
1545
1546k=0; i=0;
1547
1548while (i<10){
1549 while(k<11){
1550 sprintf(nBlocks,"ab%d",N++);
1551 mBlock = geom->MakeTubs(nBlocks,Iron, 3,6,6,10,170);
1552 mBlock->SetLineColor(13);
1553 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-283,37.5+(i*12.5),615-(i*4), new TGeoRotation("r1",0,60,0)));
1554 k++;
1555 }
1556 i++; k=0;
1557}
1558
1559
1560
1561k=0; i=0;
1562
1563while (i<6){
1564 while(k<11){
1565 sprintf(nBlocks,"ab%d",N++);
1566 mBlock = geom->MakeTubs(nBlocks,Iron, 3,6,6,10,170);
1567 mBlock->SetLineColor(13);
1568 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-335,-81.25-(i*12.5),592.5-(i*2), new TGeoRotation("r1",0,120,0)));
1569 k++;
1570 }
1571 i++; k=0;
1572}
1573
1574k=0; i=0;
1575
1576while (i<7){
1577 while(k<11){
1578 sprintf(nBlocks,"ab%d",N++);
1579 mBlock = geom->MakeTubs(nBlocks,Iron, 3,6,6,10,170);
1580 mBlock->SetLineColor(13);
1581 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-322,-69.75-(i*12.5),595-(i*2), new TGeoRotation("r1",0,120,0)));
1582 k++;
1583 }
1584 i++; k=0;
1585}
1586
1587
1588k=0; i=0;
1589
1590while (i<8){
1591 while(k<11){
1592 sprintf(nBlocks,"ab%d",N++);
1593 mBlock = geom->MakeTubs(nBlocks,Iron, 3,6,6,10,170);
1594 mBlock->SetLineColor(13);
1595 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-309,-56.25-(i*12.5),605-(i*4), new TGeoRotation("r1",0,120,0)));
1596 k++;
1597 }
1598 i++; k=0;
1599}
1600
1601k=0; i=0;
1602
1603while (i<9){
1604 while(k<11){
1605 sprintf(nBlocks,"ab%d",N++);
1606 mBlock = geom->MakeTubs(nBlocks,Iron, 3,6,6,10,170);
1607 mBlock->SetLineColor(13);
1608 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-296,-50-(i*12.5),610-(i*4), new TGeoRotation("r1",0,120,0)));
1609 k++;
1610 }
1611 i++; k=0;
1612}
1613
1614
1615k=0; i=0;
1616
1617while (i<10){
1618 while(k<11){
1619 sprintf(nBlocks,"ab%d",N++);
1620 mBlock = geom->MakeTubs(nBlocks,Iron, 3,6,6,10,170);
1621 mBlock->SetLineColor(13);
1622 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-283,-37.5-(i*12.5),615-(i*4), new TGeoRotation("r1",0,120,0)));
1623 k++;
1624 }
1625 i++; k=0;
1626}
1627
1628
1629
1630//////////right
1631
1632
1633
1634k=0; i=0;
1635
1636while (i<6){
1637 while(k<11){
1638 sprintf(nBlocks,"ab%d",N++);
1639 mBlock = geom->MakeTubs(nBlocks,Iron, 3,6,6,10,170);
1640 mBlock->SetLineColor(13);
1641 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(335,81.25+(i*12.5),592.5-(i*2), new TGeoRotation("r1",0,60,0)));
1642 k++;
1643 }
1644 i++; k=0;
1645}
1646
1647k=0; i=0;
1648
1649while (i<7){
1650 while(k<11){
1651 sprintf(nBlocks,"ab%d",N++);
1652 mBlock = geom->MakeTubs(nBlocks,Iron, 3,6,6,10,170);
1653 mBlock->SetLineColor(13);
1654 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(322,69.75+(i*12.5),595-(i*2), new TGeoRotation("r1",0,60,0)));
1655 k++;
1656 }
1657 i++; k=0;
1658}
1659
1660
1661k=0; i=0;
1662
1663while (i<8){
1664 while(k<11){
1665 sprintf(nBlocks,"ab%d",N++);
1666 mBlock = geom->MakeTubs(nBlocks,Iron, 3,6,6,10,170);
1667 mBlock->SetLineColor(13);
1668 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(309,56.25+(i*12.5),605-(i*4), new TGeoRotation("r1",0,60,0)));
1669 k++;
1670 }
1671 i++; k=0;
1672}
1673
1674k=0; i=0;
1675
1676while (i<9){
1677 while(k<11){
1678 sprintf(nBlocks,"ab%d",N++);
1679 mBlock = geom->MakeTubs(nBlocks,Iron, 3,6,6,10,170);
1680 mBlock->SetLineColor(13);
1681 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(296,50+(i*12.5),610-(i*4), new TGeoRotation("r1",0,60,0)));
1682 k++;
1683 }
1684 i++; k=0;
1685}
1686
1687
1688k=0; i=0;
1689
1690while (i<10){
1691 while(k<11){
1692 sprintf(nBlocks,"ab%d",N++);
1693 mBlock = geom->MakeTubs(nBlocks,Iron, 3,6,6,10,170);
1694 mBlock->SetLineColor(13);
1695 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(283,37.5+(i*12.5),615-(i*4), new TGeoRotation("r1",0,60,0)));
1696 k++;
1697 }
1698 i++; k=0;
1699}
1700
1701
1702//
1703
1704
1705
1706
1707
1708k=0; i=0;
1709
1710while (i<6){
1711 while(k<11){
1712 sprintf(nBlocks,"ab%d",N++);
1713 mBlock = geom->MakeTubs(nBlocks,Iron, 3,6,6,10,170);
1714 mBlock->SetLineColor(13);
1715 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(335,-81.25-(i*12.5),592.5-(i*2), new TGeoRotation("r1",0,120,0)));
1716 k++;
1717 }
1718 i++; k=0;
1719}
1720
1721k=0; i=0;
1722
1723while (i<7){
1724 while(k<11){
1725 sprintf(nBlocks,"ab%d",N++);
1726 mBlock = geom->MakeTubs(nBlocks,Iron, 3,6,6,10,170);
1727 mBlock->SetLineColor(13);
1728 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(322,-69.75-(i*12.5),595-(i*2), new TGeoRotation("r1",0,120,0)));
1729 k++;
1730 }
1731 i++; k=0;
1732}
1733
1734
1735k=0; i=0;
1736
1737while (i<8){
1738 while(k<11){
1739 sprintf(nBlocks,"ab%d",N++);
1740 mBlock = geom->MakeTubs(nBlocks,Iron, 3,6,6,10,170);
1741 mBlock->SetLineColor(13);
1742 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(309,-56.25-(i*12.5),605-(i*4), new TGeoRotation("r1",0,120,0)));
1743 k++;
1744 }
1745 i++; k=0;
1746}
1747
1748k=0; i=0;
1749
1750while (i<9){
1751 while(k<11){
1752 sprintf(nBlocks,"ab%d",N++);
1753 mBlock = geom->MakeTubs(nBlocks,Iron, 3,6,6,10,170);
1754 mBlock->SetLineColor(13);
1755 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(296,-50-(i*12.5),610-(i*4), new TGeoRotation("r1",0,120,0)));
1756 k++;
1757 }
1758 i++; k=0;
1759}
1760
1761
1762k=0; i=0;
1763
1764while (i<10){
1765 while(k<11){
1766 sprintf(nBlocks,"ab%d",N++);
1767 mBlock = geom->MakeTubs(nBlocks,Iron, 3,6,6,10,170);
1768 mBlock->SetLineColor(13);
1769 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(283,-37.5-(i*12.5),615-(i*4), new TGeoRotation("r1",0,120,0)));
1770 k++;
1771 }
1772 i++; k=0;
1773}
1774
1775
1776// /// || / / / / / / / || / / / / / / / / || / / / / / / / / / / /
1777
1778
1779 sprintf(nBlocks,"ab%d",N++);
1780 mBlock = geom->MakeBox(nBlocks, Iron, 330,10,2);
1781 mBlock->SetLineColor(42);
1782 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(0,-110,550, new TGeoRotation("r1",0,-45,0)));
1783
1784
1785
1786
1787 sprintf(nBlocks,"ab%d",N++);
1788 mBlock = geom->MakeBox(nBlocks, Iron, 330,10,2);
1789 mBlock->SetLineColor(42);
1790 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(0,110,550, new TGeoRotation("r1",0,45,0)));
1791
1792
1793
1794 sprintf(nBlocks,"ab%d",N++);
1795 mBlock = geom->MakeBox(nBlocks, Iron, 110,10,2);
1796 mBlock->SetLineColor(42);
1797 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(335,0,550, new TGeoRotation("r1",90,-45,0)));
1798
1799
1800
1801 sprintf(nBlocks,"ab%d",N++);
1802 mBlock = geom->MakeBox(nBlocks, Iron, 110,10,2);
1803 mBlock->SetLineColor(42);
1804 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-335,0,550, new TGeoRotation("r1",90,45,0)));
1805
1806
1807
1808////////////////////////////////add box
1809
1810
1811
1812
1813
1814k=0; i=0;
1815
1816while (i<5){
1817while(k<6){
1818 sprintf(nBlocks,"ab%d",N++);
1819 mBlock = geom->MakeBox(nBlocks,Iron,18,10,2);
1820 mBlock->SetLineColor(8);
1821 {
1822 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-270+(k*100),-111,550, new TGeoRotation("r1",0,-45,0)));
1823
1824 }
1825 k++;
1826}
1827 i++; k=0;
1828}
1829
1830
1831
1832
1833k=0; i=0;
1834
1835while (i<5){
1836while(k<6){
1837 sprintf(nBlocks,"ab%d",N++);
1838 mBlock = geom->MakeBox(nBlocks,Iron,18,10,2);
1839 mBlock->SetLineColor(8);
1840 {
1841 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-270+(k*100),111,550, new TGeoRotation("r1",0,45,0)));
1842
1843 }
1844 k++;
1845}
1846 i++; k=0;
1847}
1848
1849
1850 sprintf(nBlocks,"ab%d",N++);
1851 mBlock = geom->MakeBox(nBlocks, Iron, 18,10,2);
1852 mBlock->SetLineColor(8);
1853 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(336,0,550, new TGeoRotation("r1",90,-45,0)));
1854
1855
1856
1857 sprintf(nBlocks,"ab%d",N++);
1858 mBlock = geom->MakeBox(nBlocks, Iron, 18,10,2);
1859 mBlock->SetLineColor(8);
1860 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(-336,0,550, new TGeoRotation("r1",90,45,0)));
1861
1862
1863
1864
1865// | | | | |
1866
1867
1868
1869 sprintf(nBlocks,"ab%d",N++);
1870 mBlock = geom->MakeBox(nBlocks, Iron, 300,75,40);
1871 mBlock->SetLineColor(45);
1872 top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(0,0,450, new TGeoRotation("r1",0,0,0)));
1873
1874
1875
1876//kiwa
1877 sprintf(nBlocks,"ab%d",N++);
1878 mBlock = geom->MakeBox(nBlocks, Iron, 305,80,2.5);
1879 mBlock->SetLineColor(10);
1880 top->AddNodeOverlap(mBlock,1,new TGeoTranslation(0,0,430));
1881
1882
1883
1884
1885 top->SetVisibility(0);
1886 geom->CloseGeometry();
1887
1888 top->Draw("ogl");
1889}
f
#define f(i)
Definition: RSha256.hxx:104
N
#define N
TGeoCombiTrans
Class describing rotation + translation.
Definition: TGeoMatrix.h:292
TGeoManager
The manager class for any TGeo geometry.
Definition: TGeoManager.h:43
TGeoManager::CloseGeometry
void CloseGeometry(Option_t *option="d")
Closing geometry implies checking the geometry validity, fixing shapes with negative parameters (run-...
Definition: TGeoManager.cxx:1535
TGeoManager::MakeBox
TGeoVolume * MakeBox(const char *name, TGeoMedium *medium, Double_t dx, Double_t dy, Double_t dz)
Make in one step a volume pointing to a box shape with given medium.
Definition: TGeoManager.cxx:3089
TGeoManager::SetTopVolume
void SetTopVolume(TGeoVolume *vol)
Set the top volume and corresponding node as starting point of the geometry.
Definition: TGeoManager.cxx:3568
TGeoManager::SetTopVisible
void SetTopVisible(Bool_t vis=kTRUE)
make top volume visible on screen
Definition: TGeoManager.cxx:2369
TGeoManager::MakeTubs
TGeoVolume * MakeTubs(const char *name, TGeoMedium *medium, Double_t rmin, Double_t rmax, Double_t dz, Double_t phi1, Double_t phi2)
Make in one step a volume pointing to a tube segment shape with given medium.
Definition: TGeoManager.cxx:3137
TGeoMaterial
Base class describing materials.
Definition: TGeoMaterial.h:31
TGeoMedium
Media are used to store properties related to tracking and which are useful only when using geometry ...
Definition: TGeoMedium.h:24
TGeoRotation
Class describing rotations.
Definition: TGeoMatrix.h:175
TGeoTranslation
Class describing translations.
Definition: TGeoMatrix.h:122
TGeoVolume
TGeoVolume, TGeoVolumeMulti, TGeoVolumeAssembly are the volume classes.
Definition: TGeoVolume.h:53
TGeoVolume::SetVisibility
virtual void SetVisibility(Bool_t vis=kTRUE)
set visibility of this volume
Definition: TGeoVolume.cxx:2249
TGeoVolume::AddNodeOverlap
virtual void AddNodeOverlap(TGeoVolume *vol, Int_t copy_no, TGeoMatrix *mat=0, Option_t *option="")
Add a TGeoNode to the list of nodes.
Definition: TGeoVolume.cxx:1048
TGeoVolume::Draw
virtual void Draw(Option_t *option="")
draw top volume according to option
Definition: TGeoVolume.cxx:1198
TGeoVolume::SetLineColor
virtual void SetLineColor(Color_t lcolor)
Set the line color.
Definition: TGeoVolume.cxx:2087