#include "TGeoManager.h"
void cheongwadae()
{
// Drawing the Cheongwadae building which is the Presidential Residence of the Republic of Korea, using ROOT geometry class.
//
// Author: Hee Jun Shin (s-heejun@hanmail.net), Dept. of Physics, Univ. of Seoul
// Reviewed by Sunman Kim (sunman98@hanmail.net)
// Supervisor: Prof. Inkyu Park (icpark@physics.uos.ac.kr)
//
// 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
//
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(0);
// 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(0);
top->Draw("ogl");
}