Logo ROOT  

Reference Guide
 
Loading...
Searching...
No Matches
cheongwadae.C File Reference
Tutorials » Geometry tutorials

Detailed Description

Drawing the Cheongwadae building which is the Presidential Residence of the Republic of Korea, 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 cheongwadae.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 cheongwadae()
{
TGeoManager *geom = new TGeoManager("geom","My first 3D geometry");
//material
TGeoMaterial *vacuum = new TGeoMaterial("vacuum",0,0,0);
TGeoMaterial *Fe = new TGeoMaterial("Fe",55.845,26,7.87);
//creat media
TGeoMedium *Air = new TGeoMedium("Vacuum",0,vacuum);
TGeoMedium *Iron = new TGeoMedium("Iron",1,Fe);
//creat volume
TGeoVolume *top = geom->MakeBox("top",Air,300,300,300);
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);
char nBlocks[100];
int N = 0;
int f=0;
int di[2]; di[0] = 0; di[1] = 30;
TGeoVolume *mBlock;
for(int k=0;k<7;k++){
for(int i=0;i<20;i++){
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeBox(nBlocks, Iron, 0.6,1.8,63);
mBlock->SetLineColor(20);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(-10.6-(2.6*i),-17.8+(6*k),0));
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeBox(nBlocks, Iron, 0.7,1.8,58);
mBlock->SetLineColor(12);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(-11.9-(2.6*i),-17.8+(6*k),0));
}
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeBox(nBlocks, Iron, 26,1.2,63);
mBlock->SetLineColor(20);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(-36,-14.8+(6*k),0));
}
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeBox(nBlocks, Iron, 26,2,63);
mBlock->SetLineColor(20);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(-36,-21.6,0));
for(int k=0;k<7;k++){
for(int i=0;i<20;i++){
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeBox(nBlocks, Iron, 0.6,1.8,63);
mBlock->SetLineColor(20);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(-10.6-(2.6*i),-17.8+(6*k),0));
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeBox(nBlocks, Iron, 0.7,1.8,58);
mBlock->SetLineColor(12);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(-11.9-(2.6*i),-17.8+(6*k),0));
}
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeBox(nBlocks, Iron, 26,1.2,63);
mBlock->SetLineColor(20);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(-36,-14.8+(6*k),0));
}
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeBox(nBlocks, Iron, 10,22,58);
mBlock->SetLineColor(2);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(0,0,0));
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeBox(nBlocks, Iron, 3.5,8,0.1);
mBlock->SetLineColor(13);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(4,-14,60));
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeBox(nBlocks, Iron, 3.5,8,0.1);
mBlock->SetLineColor(13);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(-4,-14,60));
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeBox(nBlocks, Iron, 10,0.2,0.1);
mBlock->SetLineColor(1);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(0,20,60));
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeBox(nBlocks, Iron, 10,0.2,0.1);
mBlock->SetLineColor(1);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(0,17,60));
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeBox(nBlocks, Iron, 10,0.2,0.1);
mBlock->SetLineColor(1);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(0,14,60));
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeBox(nBlocks, Iron, 10,0.2,0.1);
mBlock->SetLineColor(1);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(0,11,60));
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeBox(nBlocks, Iron, 10,0.2,0.1);
mBlock->SetLineColor(1);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(0,8,60));
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeBox(nBlocks, Iron, 10,0.2,0.1);
mBlock->SetLineColor(1);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(0,5,60));
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeBox(nBlocks, Iron, 10,0.2,0.1);
mBlock->SetLineColor(1);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(0,2,60));
for(int k=0;k<7;k++){
for(int i=0;i<20;i++){
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeBox(nBlocks, Iron, 0.6,1.8,63);
mBlock->SetLineColor(20);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(10.6+(2.6*i),-17.8+(6*k),0));
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeBox(nBlocks, Iron, 0.7,1.8,58);
mBlock->SetLineColor(12);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(11.9+(2.6*i),-17.8+(6*k),0));
}
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeBox(nBlocks, Iron, 26,1.2,63);
mBlock->SetLineColor(20);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(36,-14.8+(6*k),0));
}
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeBox(nBlocks, Iron, 26,2,63);
mBlock->SetLineColor(20);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(36,-21.6,0));
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeBox(nBlocks, Iron, 82,2,82);
mBlock->SetLineColor(18);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(0,24,0));
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeBox(nBlocks, Iron, 85,0.5,85);
mBlock->SetLineColor(18);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(0,26,0));
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeBox(nBlocks, Iron, 88,2,88);
mBlock->SetLineColor(18);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(0,-24,0));
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeSphere(nBlocks, Iron, 0, 30, 0, 180, 0, 180);
mBlock->SetLineColor(32);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(0,24,0));
sprintf(nBlocks,"ab%d",N++);
mBlock = geom->MakeBox(nBlocks,Iron, 0.1,30,0.1);
mBlock->SetLineColor(10);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(0,40,0));
sprintf(nBlocks,"ab%d",N++);
mBlock = geom->MakeTubs(nBlocks,Iron, 0,30,4,360,360);
mBlock->SetLineColor(10);
top->AddNodeOverlap(mBlock,1,new TGeoCombiTrans(0,27,0, new TGeoRotation("r1",0,90,0)));
for(int i=0;i<8;i++){
sprintf(nBlocks,"ab%d",N++);
mBlock = geom->MakeBox(nBlocks,Iron, 2,22,2);
mBlock->SetLineColor(18);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(-70+(20*i),0,80));
}
for(int i=0;i<8;i++){
sprintf(nBlocks,"ab%d",N++);
mBlock = geom->MakeBox(nBlocks,Iron, 2,22,2);
mBlock->SetLineColor(18);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(-70+(20*i),0,-80));
}
for(int i=0;i<7;i++){
sprintf(nBlocks,"ab%d",N++);
mBlock = geom->MakeBox(nBlocks,Iron, 2,22,2);
mBlock->SetLineColor(18);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(-70,0,-80+(23*i)));
}
for(int i=0;i<7;i++){
sprintf(nBlocks,"ab%d",N++);
mBlock = geom->MakeBox(nBlocks,Iron, 2,22,2);
mBlock->SetLineColor(18);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(70,0,-80+(23*i)));
}
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeBox(nBlocks, Iron, 100,0.5,160);
mBlock->SetLineColor(41);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(0,-26,40));
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeBox(nBlocks, Iron, 10,0.01,160);
mBlock->SetLineColor(19);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(0,-25,40));
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeSphere(nBlocks, Iron, 0, 5, 0, 180, 0, 180);
mBlock->SetLineColor(8);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(15,-22,170));
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeSphere(nBlocks, Iron, 0, 5, 0, 180, 0, 180);
mBlock->SetLineColor(8);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(15,-25,170));
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeSphere(nBlocks, Iron, 0, 5, 0, 180, 0, 180);
mBlock->SetLineColor(8);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(15,-22,150));
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeSphere(nBlocks, Iron, 0, 5, 0, 180, 0, 180);
mBlock->SetLineColor(8);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(15,-25,150));
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeSphere(nBlocks, Iron, 0, 5, 0, 180, 0, 180);
mBlock->SetLineColor(8);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(15,-22,130));
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeSphere(nBlocks, Iron, 0, 5, 0, 180, 0, 180);
mBlock->SetLineColor(8);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(15,-25,130));
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeSphere(nBlocks, Iron, 0, 5, 0, 180, 0, 180);
mBlock->SetLineColor(8);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(15,-22,110));
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeSphere(nBlocks, Iron, 0, 5, 0, 180, 0, 180);
mBlock->SetLineColor(8);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(15,-25,110));
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeSphere(nBlocks, Iron, 0, 5, 0, 180, 0, 180);
mBlock->SetLineColor(8);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(-15,-22,170));
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeSphere(nBlocks, Iron, 0, 5, 0, 180, 0, 180);
mBlock->SetLineColor(8);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(-15,-25,170));
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeSphere(nBlocks, Iron, 0, 5, 0, 180, 0, 180);
mBlock->SetLineColor(8);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(-15,-22,150));
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeSphere(nBlocks, Iron, 0, 5, 0, 180, 0, 180);
mBlock->SetLineColor(8);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(-15,-25,150));
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeSphere(nBlocks, Iron, 0, 5, 0, 180, 0, 180);
mBlock->SetLineColor(8);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(-15,-22,130));
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeSphere(nBlocks, Iron, 0, 5, 0, 180, 0, 180);
mBlock->SetLineColor(8);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(-15,-25,130));
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeSphere(nBlocks, Iron, 0, 5, 0, 180, 0, 180);
mBlock->SetLineColor(8);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(-15,-22,110));
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeSphere(nBlocks, Iron, 0, 5, 0, 180, 0, 180);
mBlock->SetLineColor(8);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(-15,-25,110));
sprintf(nBlocks,"ab%d",N++);
mBlock = geom->MakeBox(nBlocks,Iron, 0.1,10,0.1);
mBlock->SetLineColor(12);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(20,-15,110));
sprintf(nBlocks,"ab%d",N++);
mBlock = geom->MakeBox(nBlocks,Iron, 5,3,0.1);
mBlock->SetLineColor(10);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(25,-8,110));
sprintf(nBlocks,"ab%d",N++);
mBlock = geom->MakeBox(nBlocks,Iron, 0.1,10,0.1);
mBlock->SetLineColor(12);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(-20,-15,110));
sprintf(nBlocks,"ab%d",N++);
mBlock = geom->MakeBox(nBlocks,Iron, 5,3,0.1);
mBlock->SetLineColor(10);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(-15,-8,110));
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeBox(nBlocks, Iron, 7,1.5,5);
mBlock->SetLineColor(18);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(0,-24,88));
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeBox(nBlocks, Iron, 7,1,5);
mBlock->SetLineColor(18);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(0,-24,92));
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeBox(nBlocks, Iron, 7,0.5,5);
mBlock->SetLineColor(18);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(0,-24,96));
sprintf(nBlocks,"f%d_bg%d",f,N++);
mBlock = geom->MakeBox(nBlocks, Iron, 7,0.1,5);
mBlock->SetLineColor(18);
top->AddNodeOverlap(mBlock,1,new TGeoTranslation(0,-24,100));
geom->CloseGeometry();
top->SetVisibility(false);
top->Draw("ogl");
}
f
#define f(i)
Definition RSha256.hxx:104
TRangeDynCast
ROOT::Detail::TRangeCast< T, true > TRangeDynCast
TRangeDynCast is an adapter class that allows the typed iteration through a TCollection.
Definition TCollection.h:358
N
#define N
TGeoCombiTrans
Class describing rotation + translation.
Definition TGeoMatrix.h:317
TGeoManager
The manager class for any TGeo geometry.
Definition TGeoManager.h:44
TGeoMaterial
Base class describing materials.
Definition TGeoMaterial.h:34
TGeoMedium
Media are used to store properties related to tracking and which are useful only when using geometry ...
Definition TGeoMedium.h:23
TGeoRotation
Class describing rotations.
Definition TGeoMatrix.h:168
TGeoTranslation
Class describing translations.
Definition TGeoMatrix.h:116
TGeoVolume
TGeoVolume, TGeoVolumeMulti, TGeoVolumeAssembly are the volume classes.
Definition TGeoVolume.h:43
TGeoVolume::SetVisibility
void SetVisibility(Bool_t vis=kTRUE) override
set visibility of this volume
Definition TGeoVolume.cxx:2347
TGeoVolume::Draw
void Draw(Option_t *option="") override
draw top volume according to option
Definition TGeoVolume.cxx:1206
TGeoVolume::AddNodeOverlap
virtual void AddNodeOverlap(TGeoVolume *vol, Int_t copy_no, TGeoMatrix *mat=nullptr, Option_t *option="")
Add a TGeoNode to the list of nodes.
Definition TGeoVolume.cxx:1044
Author
Hee Jun Shin (s-hee.nosp@m.jun@.nosp@m.hanma.nosp@m.il.n.nosp@m.et), Dept. of Physics, Univ. of Seoul

Definition in file cheongwadae.C.