station1.C File Reference

Detailed Description

Drawing a space station, using ROOT geometry class.

Reviewed by Sunman Kim (sunma.nosp@m.n98@.nosp@m.hanma.nosp@m.il.n.nosp@m.et) Supervisor: Prof. Inkyu Park (icpar.nosp@m.k@ph.nosp@m.ysics.nosp@m..uos.nosp@m..ac.k.nosp@m.r)

How to run: .x station1.C in ROOT terminal, then use OpenGL

This macro was created for the evaluation of Computational Physics course in 2006. We thank to Prof. Inkyu Park for his special lecture on ROOT and to all of ROOT team

 
#include "TGeoManager.h"
 
void station1()
{
 
   TGeoManager *geom = new TGeoManager("geom", "My first 3D geometry");
 
   TGeoMaterial *vacuum = new TGeoMaterial("vacuum", 0, 0, 0);
   TGeoMaterial *Fe = new TGeoMaterial("Fe", 55.845, 26, 7.87);
   TGeoMaterial *Cu = new TGeoMaterial("Cu", 63.549, 29, 8.92);
 
   TGeoMedium *Air = new TGeoMedium("Vacuum", 0, vacuum);
   TGeoMedium *Iron = new TGeoMedium("Iron", 1, Fe);
   TGeoMedium *Copper = new TGeoMedium("Copper", 2, Cu);
 
   TGeoVolume *top = geom->MakeBox("top", Air, 1000, 1000, 1000);
   geom->SetTopVolume(top);
   geom->SetTopVisible(false);
   // If you want to see the boundary, please input the number, 1 instead of 0.
   // Like this, geom->SetTopVisible(1);
 
   TGeoVolume *Cone1 = geom->MakeCone("Cone1", Copper, 650, 0, 20, 0, 20);
   Cone1->SetFillColor(35);
   Cone1->SetLineColor(35);
   top->AddNodeOverlap(Cone1, 1, new TGeoTranslation(0, 0, 0));
 
   TGeoVolume *Cone2 = geom->MakeCone("Cone2", Copper, 25, 0, 30, 0, 30);
   Cone2->SetFillColor(7);
   Cone2->SetLineColor(7);
   top->AddNodeOverlap(Cone2, 1, new TGeoTranslation(0, 0, 630));
 
   TGeoVolume *Cone21 = geom->MakeCone("Cone21", Copper, 30, 0, 30, 0, 30);
   Cone21->SetFillColor(29);
   Cone21->SetLineColor(29);
   top->AddNodeOverlap(Cone21, 1, new TGeoTranslation(0, 0, 550));
 
   TGeoVolume *Cone22 = geom->MakeCone("Cone22", Copper, 5, 0, 50, 0, 50);
   Cone22->SetFillColor(2);
   Cone22->SetLineColor(2);
   top->AddNodeOverlap(Cone22, 1, new TGeoTranslation(0, 0, 500));
 
   for (int i = 0; i < 28; i++) {
      TGeoVolume *Cone00 = geom->MakeCone("Cone00", Copper, 3, 0, 25, 0, 25);
      Cone00->SetFillColor(1);
      Cone00->SetLineColor(1);
      top->AddNodeOverlap(Cone00, 1, new TGeoTranslation(0, 0, -100 + 20 * i));
   }
 
   TGeoVolume *Cone3 = geom->MakeCone("Cone3", Copper, 60, 0, 70, 0, 0);
   Cone3->SetFillColor(13);
   Cone3->SetLineColor(13);
   top->AddNodeOverlap(Cone3, 1, new TGeoTranslation(-60, 0, -110));
 
   TGeoVolume *Cone31 = geom->MakeCone("Cone31", Copper, 230, 0, 70, 0, 70);
   Cone31->SetFillColor(13);
   Cone31->SetLineColor(13);
   top->AddNodeOverlap(Cone31, 1, new TGeoTranslation(-60, 0, -400));
 
   for (int i = 0; i < 5; i++) {
      Cone31 = geom->MakeCone("Cone31", Copper, 7, 0, 73, 0, 73);
      Cone31->SetFillColor(21);
      Cone31->SetLineColor(21);
      top->AddNodeOverlap(Cone31, 1, new TGeoTranslation(-60, 0, -170 - (500 / 6 * (i + 1))));
      top->AddNodeOverlap(Cone31, 1, new TGeoTranslation(60, 0, -170 - (500 / 6 * (i + 1))));
      top->AddNodeOverlap(Cone31, 1, new TGeoTranslation(0, -60, -170 - (500 / 6 * (i + 1))));
      top->AddNodeOverlap(Cone31, 1, new TGeoTranslation(0, 60, -170 - (500 / 6 * (i + 1))));
   }
 
   TGeoVolume *Cone32 = geom->MakeCone("Cone32", Copper, 30, 60, 50, 0, 70);
   Cone32->SetFillColor(35);
   Cone32->SetLineColor(35);
   top->AddNodeOverlap(Cone32, 1, new TGeoTranslation(-60, 0, -650));
 
   TGeoVolume *Cone321 = geom->MakeCone("Cone321", Copper, 5, 60, 50, 0, 50);
   Cone321->SetFillColor(2);
   Cone321->SetLineColor(2);
   top->AddNodeOverlap(Cone321, 1, new TGeoTranslation(-60, 0, -680));
 
   TGeoVolume *Cone4 = geom->MakeCone("Cone4", Copper, 60, 0, 70, 0, 0);
   Cone4->SetFillColor(13);
   Cone4->SetLineColor(13);
   top->AddNodeOverlap(Cone3, 1, new TGeoTranslation(60, 0, -110));
 
   for (int i = 1; i <= 8; i++) {
      TGeoVolume *Torus2 = geom->MakeTorus("Torus2", Iron, 120, 20, 40, 45 * i - 4, 8);
      Torus2->SetFillColor(18);
      Torus2->SetLineColor(18);
      top->AddNodeOverlap(Torus2, 1, new TGeoTranslation(0, 0, 610));
 
      TGeoVolume *Tubs = geom->MakeTubs("Line", Iron, 0, 190, 5, 45 * i - 1, 45 * i + 1);
      Tubs->SetFillColor(18);
      Tubs->SetLineColor(18);
      top->AddNodeOverlap(Tubs, 1, new TGeoTranslation(0, 0, 610));
   }
 
   TGeoVolume *Cone41 = geom->MakeCone("Cone41", Copper, 230, 0, 70, 0, 70);
   Cone41->SetFillColor(13);
   Cone41->SetLineColor(13);
   top->AddNodeOverlap(Cone41, 1, new TGeoTranslation(60, 0, -400));
 
   TGeoVolume *Cone42 = geom->MakeCone("Cone42", Copper, 30, 60, 50, 0, 70);
   Cone42->SetFillColor(35);
   Cone42->SetLineColor(35);
   top->AddNodeOverlap(Cone42, 1, new TGeoTranslation(60, 0, -650));
 
   TGeoVolume *Cone421 = geom->MakeCone("Cone421", Copper, 5, 60, 50, 0, 50);
   Cone421->SetFillColor(2);
   Cone421->SetLineColor(2);
   top->AddNodeOverlap(Cone421, 1, new TGeoTranslation(60, 0, -680));
 
   TGeoVolume *Cone5 = geom->MakeCone("Cone5", Copper, 60, 0, 70, 0, 0);
   Cone5->SetFillColor(13);
   Cone5->SetLineColor(13);
   top->AddNodeOverlap(Cone3, 1, new TGeoTranslation(0, -60, -110));
 
   TGeoVolume *Cone51 = geom->MakeCone("Cone51", Copper, 230, 0, 70, 0, 70);
   Cone51->SetFillColor(13);
   Cone51->SetLineColor(13);
   top->AddNodeOverlap(Cone51, 1, new TGeoTranslation(0, -60, -400));
 
   TGeoVolume *Cone52 = geom->MakeCone("Cone52", Copper, 30, 60, 50, 0, 70);
   Cone52->SetFillColor(35);
   Cone52->SetLineColor(35);
   top->AddNodeOverlap(Cone52, 1, new TGeoTranslation(0, -60, -650));
 
   TGeoVolume *Cone521 = geom->MakeCone("Cone521", Copper, 5, 60, 50, 0, 50);
   Cone521->SetFillColor(2);
   Cone521->SetLineColor(2);
   top->AddNodeOverlap(Cone521, 1, new TGeoTranslation(0, -60, -680));
 
   TGeoVolume *Cone6 = geom->MakeCone("Cone6", Copper, 60, 0, 70, 0, 0);
   Cone6->SetFillColor(13);
   Cone6->SetLineColor(13);
   top->AddNodeOverlap(Cone3, 1, new TGeoTranslation(0, 60, -110));
 
   TGeoVolume *Cone61 = geom->MakeCone("Cone61", Copper, 230, 0, 70, 0, 70);
   Cone61->SetFillColor(13);
   Cone61->SetLineColor(13);
   top->AddNodeOverlap(Cone61, 1, new TGeoTranslation(0, 60, -400));
 
   TGeoVolume *Cone62 = geom->MakeCone("Cone62", Copper, 30, 60, 50, 0, 70);
   Cone62->SetFillColor(35);
   Cone62->SetLineColor(35);
   top->AddNodeOverlap(Cone62, 1, new TGeoTranslation(0, 60, -650));
 
   TGeoVolume *Cone621 = geom->MakeCone("Cone621", Copper, 5, 60, 50, 0, 50);
   Cone621->SetFillColor(2);
   Cone621->SetLineColor(2);
   top->AddNodeOverlap(Cone621, 1, new TGeoTranslation(0, 60, -680));
 
   TGeoVolume *Cone7 = geom->MakeCone("Cone7", Copper, 50, 0, 40, 0, 5);
   Cone7->SetFillColor(13);
   Cone7->SetLineColor(13);
   top->AddNodeOverlap(Cone7, 1, new TGeoCombiTrans(-90, -60, 10, new TGeoRotation("Cone7", 90, -90, -90)));
 
   TGeoVolume *Cone71 = geom->MakeCone("Cone71", Copper, 50, 0, 60, 0, 40);
   Cone71->SetFillColor(16);
   Cone71->SetLineColor(16);
   top->AddNodeOverlap(Cone71, 1, new TGeoCombiTrans(10, -60, 10, new TGeoRotation("Cone7", 90, -90, -90)));
 
   TGeoVolume *Cone711 = geom->MakeCone("Cone711", Copper, 10, 0, 10, 0, 60);
   Cone711->SetFillColor(13);
   Cone711->SetLineColor(13);
   top->AddNodeOverlap(Cone711, 1, new TGeoCombiTrans(70, -60, 10, new TGeoRotation("Cone7", 90, -90, -90)));
 
   TGeoVolume *Torus1 = geom->MakeTorus("Torus1", Iron, 120, 30, 20);
   Torus1->SetFillColor(33);
   Torus1->SetLineColor(33);
   top->AddNodeOverlap(Torus1, 1, new TGeoTranslation(0, 0, 610));
 
   TGeoVolume *Cone8 = geom->MakeCone("Cone8", Copper, 50, 0, 40, 0, 5);
   Cone8->SetFillColor(13);
   Cone8->SetLineColor(13);
   top->AddNodeOverlap(Cone8, 1, new TGeoCombiTrans(100, 60, 10, new TGeoRotation("Cone8", 90, 90, 0)));
 
   TGeoVolume *Cone81 = geom->MakeCone("Cone81", Copper, 50, 0, 60, 0, 40);
   Cone81->SetFillColor(16);
   Cone81->SetLineColor(16);
   top->AddNodeOverlap(Cone81, 1, new TGeoCombiTrans(0, 60, 10, new TGeoRotation("Cone8", 90, 90, 0)));
 
   TGeoVolume *Cone811 = geom->MakeCone("Cone811", Copper, 10, 0, 10, 0, 60);
   Cone811->SetFillColor(13);
   Cone811->SetLineColor(13);
   top->AddNodeOverlap(Cone811, 1, new TGeoCombiTrans(-60, 60, 10, new TGeoRotation("Cone8", 90, 90, 0)));
 
   TGeoVolume *Box1 = geom->MakeBox("Box1", Copper, 10, 10, 3);
   Box1->SetFillColor(3);
   Box1->SetLineColor(3);
   top->AddNodeOverlap(Box1, 1, new TGeoCombiTrans(-110, -50, 645, new TGeoRotation("Box1", 0, 0, 30)));
 
   TGeoVolume *Box2 = geom->MakeBox("Box2", Copper, 10, 10, 3);
   Box2->SetFillColor(3);
   Box2->SetLineColor(3);
   top->AddNodeOverlap(Box2, 1, new TGeoCombiTrans(110, 45, 645, new TGeoRotation("Box2", 0, 0, 30)));
 
   TGeoVolume *Box3 = geom->MakeBox("Box3", Copper, 10, 10, 3);
   Box3->SetFillColor(3);
   Box3->SetLineColor(3);
   top->AddNodeOverlap(Box3, 1, new TGeoCombiTrans(-45, -110, 645, new TGeoRotation("Box3", 0, 0, 70)));
 
   TGeoVolume *Box4 = geom->MakeBox("Box4", Copper, 10, 10, 3);
   Box4->SetFillColor(3);
   Box4->SetLineColor(3);
   top->AddNodeOverlap(Box4, 1, new TGeoCombiTrans(45, 110, 645, new TGeoRotation("Box4", 0, 0, 70)));
 
   TGeoVolume *Box5 = geom->MakeBox("Box5", Copper, 10, 10, 3);
   Box5->SetFillColor(3);
   Box5->SetLineColor(3);
   top->AddNodeOverlap(Box5, 1, new TGeoCombiTrans(45, -110, 645, new TGeoRotation("Box5", 0, 0, 30)));
 
   TGeoVolume *Box6 = geom->MakeBox("Box6", Copper, 10, 10, 3);
   Box6->SetFillColor(3);
   Box6->SetLineColor(3);
   top->AddNodeOverlap(Box6, 1, new TGeoCombiTrans(-45, 110, 645, new TGeoRotation("Box6", 0, 0, 25)));
 
   TGeoVolume *Box7 = geom->MakeBox("Box7", Copper, 10, 10, 3);
   Box7->SetFillColor(3);
   Box7->SetLineColor(3);
   top->AddNodeOverlap(Box7, 1, new TGeoCombiTrans(110, -50, 645, new TGeoRotation("Box7", 0, 0, 60)));
 
   TGeoVolume *Box8 = geom->MakeBox("Box8", Copper, 10, 10, 3);
   Box8->SetFillColor(3);
   Box8->SetLineColor(3);
   top->AddNodeOverlap(Box8, 1, new TGeoCombiTrans(-110, 45, 645, new TGeoRotation("Box8", 0, 0, 60)));
 
   Torus1 = geom->MakeTorus("Torus1", Iron, 120, 30, 20);
   Torus1->SetFillColor(33);
   Torus1->SetLineColor(33);
   top->AddNodeOverlap(Torus1, 1, new TGeoTranslation(0, 0, 610));
 
   for (int i = 1; i <= 8; i++) {
      TGeoVolume *Torus2 = geom->MakeTorus("Torus2", Iron, 120, 20, 40, 45 * i - 4, 8);
      Torus2->SetFillColor(18);
      Torus2->SetLineColor(18);
      top->AddNodeOverlap(Torus2, 1, new TGeoTranslation(0, 0, 610));
 
      TGeoVolume *Tubs = geom->MakeTubs("Line", Iron, 0, 190, 5, 45 * i - 1, 45 * i + 1);
      Tubs->SetFillColor(18);
      Tubs->SetLineColor(18);
      top->AddNodeOverlap(Tubs, 1, new TGeoTranslation(0, 0, 610));
   }
 
   TGeoVolume *Sphere00 = geom->MakeSphere("Sphere00", Iron, 0, 15, 0, 45, 0);
   Sphere00->SetFillColor(2);
   Sphere00->SetLineColor(2);
   top->AddNodeOverlap(Sphere00, 1, new TGeoTranslation(-145, -145, 600));
 
   TGeoVolume *Sphere01 = geom->MakeSphere("Sphere01", Iron, 0, 15, 0, 45, 0);
   Sphere01->SetFillColor(2);
   Sphere01->SetLineColor(2);
   top->AddNodeOverlap(Sphere01, 1, new TGeoTranslation(0, -210, 600));
 
   TGeoVolume *Sphere02 = geom->MakeSphere("Sphere02", Iron, 0, 15, 0, 45, 0);
   Sphere02->SetFillColor(2);
   Sphere02->SetLineColor(2);
   top->AddNodeOverlap(Sphere02, 1, new TGeoTranslation(145, 145, 600));
 
   TGeoVolume *Sphere03 = geom->MakeSphere("Sphere03", Iron, 0, 15, 0, 45, 0);
   Sphere03->SetFillColor(2);
   Sphere03->SetLineColor(2);
   top->AddNodeOverlap(Sphere03, 1, new TGeoTranslation(0, 210, 600));
 
   TGeoVolume *Sphere04 = geom->MakeSphere("Sphere04", Iron, 0, 15, 0, 45, 0);
   Sphere04->SetFillColor(2);
   Sphere04->SetLineColor(2);
   top->AddNodeOverlap(Sphere04, 1, new TGeoTranslation(145, -145, 600));
 
   TGeoVolume *Sphere05 = geom->MakeSphere("Sphere05", Iron, 0, 15, 0, 45, 0);
   Sphere05->SetFillColor(2);
   Sphere05->SetLineColor(2);
   top->AddNodeOverlap(Sphere05, 1, new TGeoTranslation(-210, 0, 600));
 
   TGeoVolume *Sphere06 = geom->MakeSphere("Sphere06", Iron, 0, 15, 0, 45, 0);
   Sphere06->SetFillColor(2);
   Sphere06->SetLineColor(2);
   top->AddNodeOverlap(Sphere06, 1, new TGeoTranslation(210, 0, 600));
 
   TGeoVolume *Sphere07 = geom->MakeSphere("Sphere07", Iron, 0, 15, 0, 45, 0);
   Sphere07->SetFillColor(2);
   Sphere07->SetLineColor(2);
   top->AddNodeOverlap(Sphere07, 1, new TGeoTranslation(-145, 145, 600));
 
   TGeoVolume *Torus3 = geom->MakeTorus("Torus3", Iron, 190, 0, 10);
   Torus3->SetFillColor(18);
   Torus3->SetLineColor(18);
   top->AddNodeOverlap(Torus3, 1, new TGeoTranslation(0, 0, 610));
 
   TGeoVolume *Sphere1 = geom->MakeSphere("Sphere1", Iron, 0, 20, 0, 180, 0, 360);
   Sphere1->SetFillColor(2);
   Sphere1->SetLineColor(2);
   top->AddNodeOverlap(Sphere1, 1, new TGeoTranslation(0, 0, 650));
 
   TGeoVolume *Tubs = geom->MakeTubs("Tubs", Iron, 0, 40, 50, 0, 360);
   Tubs->SetFillColor(29);
   Tubs->SetLineColor(29);
   top->AddNodeOverlap(Tubs, 1, new TGeoTranslation(0, 0, 500));
 
   TGeoVolume *Tubs1 = geom->MakeTubs("Tubs1", Iron, 50, 60, 230, 40, 150);
   Tubs1->SetFillColor(18);
   Tubs1->SetLineColor(18);
   top->AddNodeOverlap(Tubs1, 1, new TGeoTranslation(-170, -30, -400));
 
   TGeoVolume *Tubs11 = geom->MakeTubs("Tubs11", Iron, 50, 60, 230, 220, 330);
   Tubs11->SetFillColor(18);
   Tubs11->SetLineColor(18);
   top->AddNodeOverlap(Tubs11, 1, new TGeoTranslation(-260, 35, -400));
 
   TGeoVolume *Sphere111 = geom->MakeSphere("Sphere111", Iron, 0, 10, 0, 180, 0, 360);
   Sphere111->SetFillColor(2);
   Sphere111->SetLineColor(2);
   top->AddNodeOverlap(Sphere111, 1, new TGeoTranslation(-310, 0, -165));
 
   TGeoVolume *Sphere112 = geom->MakeSphere("Sphere112", Iron, 0, 10, 0, 180, 0, 360);
   Sphere112->SetFillColor(2);
   Sphere112->SetLineColor(2);
   top->AddNodeOverlap(Sphere112, 1, new TGeoTranslation(-310, 0, -400));
 
   TGeoVolume *Sphere113 = geom->MakeSphere("Sphere113", Iron, 0, 10, 0, 180, 0, 360);
   Sphere113->SetFillColor(2);
   Sphere113->SetLineColor(2);
   top->AddNodeOverlap(Sphere113, 1, new TGeoTranslation(-310, 0, -635));
 
   TGeoVolume *Tubs2 = geom->MakeTubs("Tubs2", Iron, 50, 60, 230, 220, 330);
   Tubs2->SetFillColor(18);
   Tubs2->SetLineColor(18);
   top->AddNodeOverlap(Tubs2, 1, new TGeoTranslation(170, 30, -400));
 
   TGeoVolume *Tubs21 = geom->MakeTubs("Tubs21", Iron, 50, 60, 230, 400, 510);
   Tubs21->SetFillColor(18);
   Tubs21->SetLineColor(18);
   top->AddNodeOverlap(Tubs21, 1, new TGeoTranslation(265, -25, -400));
 
   TGeoVolume *Sphere211 = geom->MakeSphere("Sphere211", Iron, 0, 10, 0, 180, 0, 360);
   Sphere211->SetFillColor(2);
   Sphere211->SetLineColor(2);
   top->AddNodeOverlap(Sphere211, 1, new TGeoTranslation(310, 0, -165));
 
   TGeoVolume *Sphere212 = geom->MakeSphere("Sphere212", Iron, 0, 10, 0, 180, 0, 360);
   Sphere212->SetFillColor(2);
   Sphere212->SetLineColor(2);
   top->AddNodeOverlap(Sphere212, 1, new TGeoTranslation(310, 0, -400));
 
   TGeoVolume *Sphere213 = geom->MakeSphere("Sphere213", Iron, 0, 10, 0, 180, 0, 360);
   Sphere213->SetFillColor(2);
   Sphere213->SetLineColor(2);
   top->AddNodeOverlap(Sphere213, 1, new TGeoTranslation(310, 0, -635));
 
   TGeoVolume *Tubs3 = geom->MakeTubs("Tubs3", Iron, 50, 60, 230, 130, 260);
   Tubs3->SetFillColor(18);
   Tubs3->SetLineColor(18);
   top->AddNodeOverlap(Tubs3, 1, new TGeoTranslation(30, -170, -400));
 
   TGeoVolume *Tubs31 = geom->MakeTubs("Tubs31", Iron, 50, 60, 230, 310, 440);
   Tubs31->SetFillColor(18);
   Tubs31->SetLineColor(18);
   top->AddNodeOverlap(Tubs31, 1, new TGeoTranslation(0, -275, -400));
 
   TGeoVolume *Sphere311 = geom->MakeSphere("Sphere311", Iron, 0, 10, 0, 180, 0, 360);
   Sphere311->SetFillColor(2);
   Sphere311->SetLineColor(2);
   top->AddNodeOverlap(Sphere311, 1, new TGeoTranslation(-35, 320, -165));
 
   TGeoVolume *Sphere312 = geom->MakeSphere("Sphere312", Iron, 0, 10, 0, 180, 0, 360);
   Sphere312->SetFillColor(2);
   Sphere312->SetLineColor(2);
   top->AddNodeOverlap(Sphere312, 1, new TGeoTranslation(-35, 320, -400));
 
   TGeoVolume *Sphere313 = geom->MakeSphere("Sphere313", Iron, 0, 10, 0, 180, 0, 360);
   Sphere313->SetFillColor(2);
   Sphere313->SetLineColor(2);
   top->AddNodeOverlap(Sphere313, 1, new TGeoTranslation(-35, 320, -635));
 
   TGeoVolume *Tubs4 = geom->MakeTubs("Tubs4", Iron, 50, 60, 230, 310, 440);
   Tubs4->SetFillColor(18);
   Tubs4->SetLineColor(18);
   top->AddNodeOverlap(Tubs4, 1, new TGeoTranslation(-30, 170, -400));
 
   TGeoVolume *Tubs41 = geom->MakeTubs("Tubs41", Iron, 50, 60, 230, 490, 620);
   Tubs41->SetFillColor(18);
   Tubs41->SetLineColor(18);
   top->AddNodeOverlap(Tubs41, 1, new TGeoTranslation(0, 275, -400));
 
   TGeoVolume *Sphere411 = geom->MakeSphere("Sphere411", Iron, 0, 10, 0, 180, 0, 360);
   Sphere411->SetFillColor(2);
   Sphere411->SetLineColor(2);
   top->AddNodeOverlap(Sphere411, 1, new TGeoTranslation(30, -320, -165));
 
   TGeoVolume *Sphere412 = geom->MakeSphere("Sphere412", Iron, 0, 10, 0, 180, 0, 360);
   Sphere412->SetFillColor(2);
   Sphere412->SetLineColor(2);
   top->AddNodeOverlap(Sphere412, 1, new TGeoTranslation(30, -320, -400));
 
   TGeoVolume *Sphere413 = geom->MakeSphere("Sphere413", Iron, 0, 10, 0, 180, 0, 360);
   Sphere413->SetFillColor(2);
   Sphere413->SetLineColor(2);
   top->AddNodeOverlap(Sphere413, 1, new TGeoTranslation(30, -320, -635));
 
   TGeoVolume *Cone010 = geom->MakeCone("Cone010", Iron, 30, 0, 30, 0, 30);
   Cone010->SetFillColor(2);
   Cone010->SetLineColor(2);
   top->AddNodeOverlap(Cone010, 1, new TGeoTranslation(0, 0, 250));
 
   TGeoVolume *Torus010 = geom->MakeTorus("Torus010", Iron, 300, 50, 40);
   Torus010->SetFillColor(33);
   Torus010->SetLineColor(33);
   top->AddNodeOverlap(Torus010, 1, new TGeoTranslation(0, 0, 250));
 
   TGeoVolume *Torus011 = geom->MakeTorus("Torus011", Iron, 400, 10, 10);
   Torus011->SetFillColor(33);
   Torus011->SetLineColor(33);
   top->AddNodeOverlap(Torus011, 1, new TGeoTranslation(0, 0, 250));
 
   TGeoVolume *Torus012 = geom->MakeTorus("Torus012", Iron, 200, 10, 10);
   Torus012->SetFillColor(33);
   Torus012->SetLineColor(33);
   top->AddNodeOverlap(Torus012, 1, new TGeoTranslation(0, 0, 250));
 
   TGeoVolume *Sphere010 = geom->MakeSphere("Sphere010", Iron, 0, 10, 0, 180, 0, 360);
   Sphere010->SetFillColor(2);
   Sphere010->SetLineColor(2);
   top->AddNodeOverlap(Sphere010, 1, new TGeoTranslation(-290, -290, 250));
 
   TGeoVolume *Sphere011 = geom->MakeSphere("Sphere011", Iron, 0, 10, 0, 180, 0, 360);
   Sphere011->SetFillColor(2);
   Sphere011->SetLineColor(2);
   top->AddNodeOverlap(Sphere011, 1, new TGeoTranslation(290, 290, 250));
 
   TGeoVolume *Sphere012 = geom->MakeSphere("Sphere012", Iron, 0, 10, 0, 180, 0, 360);
   Sphere012->SetFillColor(2);
   Sphere012->SetLineColor(2);
   top->AddNodeOverlap(Sphere012, 1, new TGeoTranslation(0, -410, 250));
 
   TGeoVolume *Sphere013 = geom->MakeSphere("Sphere013", Iron, 0, 10, 0, 180, 0, 360);
   Sphere013->SetFillColor(2);
   Sphere013->SetLineColor(2);
   top->AddNodeOverlap(Sphere013, 1, new TGeoTranslation(0, 410, 250));
 
   TGeoVolume *Sphere014 = geom->MakeSphere("Sphere014", Iron, 0, 10, 0, 180, 0, 360);
   Sphere014->SetFillColor(2);
   Sphere014->SetLineColor(2);
   top->AddNodeOverlap(Sphere014, 1, new TGeoTranslation(290, -290, 250));
 
   TGeoVolume *Sphere015 = geom->MakeSphere("Sphere015", Iron, 0, 10, 0, 180, 0, 360);
   Sphere015->SetFillColor(2);
   Sphere015->SetLineColor(2);
   top->AddNodeOverlap(Sphere015, 1, new TGeoTranslation(-290, 290, 250));
 
   TGeoVolume *Sphere016 = geom->MakeSphere("Sphere016", Iron, 0, 10, 0, 180, 0, 360);
   Sphere016->SetFillColor(2);
   Sphere016->SetLineColor(2);
   top->AddNodeOverlap(Sphere016, 1, new TGeoTranslation(410, 0, 250));
 
   TGeoVolume *Sphere017 = geom->MakeSphere("Sphere017", Iron, 0, 10, 0, 180, 0, 360);
   Sphere017->SetFillColor(2);
   Sphere017->SetLineColor(2);
   top->AddNodeOverlap(Sphere017, 1, new TGeoTranslation(-410, 0, 250));
 
   TGeoVolume *Box010 = geom->MakeBox("Box010", Copper, 10, 10, 3);
   Box010->SetFillColor(3);
   Box010->SetLineColor(3);
   top->AddNodeOverlap(Box1, 1, new TGeoCombiTrans(-120, -280, 300, new TGeoRotation("Box010", 0, 0, 70)));
 
   TGeoVolume *Box011 = geom->MakeBox("Box011", Copper, 10, 10, 3);
   Box011->SetFillColor(3);
   Box011->SetLineColor(3);
   top->AddNodeOverlap(Box011, 1, new TGeoCombiTrans(120, 280, 300, new TGeoRotation("Box011", 0, 0, 70)));
 
   TGeoVolume *Box012 = geom->MakeBox("Box012", Copper, 10, 10, 3);
   Box012->SetFillColor(3);
   Box012->SetLineColor(3);
   top->AddNodeOverlap(Box012, 1, new TGeoCombiTrans(120, -280, 300, new TGeoRotation("Box012", 0, 0, 30)));
 
   TGeoVolume *Box013 = geom->MakeBox("Box013", Copper, 10, 10, 3);
   Box013->SetFillColor(3);
   Box013->SetLineColor(3);
   top->AddNodeOverlap(Box013, 1, new TGeoCombiTrans(-120, 280, 300, new TGeoRotation("Box013", 0, 0, 30)));
 
   TGeoVolume *Box014 = geom->MakeBox("Box010", Copper, 10, 10, 3);
   Box014->SetFillColor(3);
   Box014->SetLineColor(3);
   top->AddNodeOverlap(Box014, 1, new TGeoCombiTrans(270, -120, 300, new TGeoRotation("Box014", 0, 0, 70)));
 
   TGeoVolume *Box015 = geom->MakeBox("Box015", Copper, 10, 10, 3);
   Box015->SetFillColor(3);
   Box015->SetLineColor(3);
   top->AddNodeOverlap(Box015, 1, new TGeoCombiTrans(-270, 120, 300, new TGeoRotation("Box015", 0, 0, 70)));
 
   TGeoVolume *Box016 = geom->MakeBox("Box016", Copper, 10, 10, 3);
   Box016->SetFillColor(3);
   Box016->SetLineColor(3);
   top->AddNodeOverlap(Box016, 1, new TGeoCombiTrans(270, 100, 300, new TGeoRotation("Box016", 0, 0, 30)));
 
   TGeoVolume *Box017 = geom->MakeBox("Box017", Copper, 10, 10, 3);
   Box017->SetFillColor(3);
   Box017->SetLineColor(3);
   top->AddNodeOverlap(Box017, 1, new TGeoCombiTrans(-270, -120, 300, new TGeoRotation("Box017", 0, 0, 30)));
 
   for (int i = 1; i <= 8; i++) {
      TGeoVolume *Torus0101 = geom->MakeTorus("Torus0101", Iron, 300, 70, 40, 45 * i - 4, 8);
      Torus0101->SetFillColor(18);
      Torus0101->SetLineColor(18);
      top->AddNodeOverlap(Torus0101, 1, new TGeoTranslation(0, 0, 250));
 
      TGeoVolume *Tubs0101 = geom->MakeTubs("Line", Iron, 0, 400, 5, 45 * i - 1, 45 * i + 1);
      Tubs0101->SetFillColor(18);
      Tubs0101->SetLineColor(18);
      top->AddNodeOverlap(Tubs0101, 1, new TGeoTranslation(0, 0, 250));
   }
 
   Cone31->SetFillColor(38);
   top->SetVisibility(false);
   geom->CloseGeometry();
 
   top->Draw("ogl");
}

Author

Chang Yeol Lee, Dept. of Physics, Univ. of Seoul

Definition in file station1.C.