FairRoot/PandaRoot
hypGe_TripleGeo.C
Go to the documentation of this file.
1 #include "TGeoManager.h"
2 double Pi = TMath::Pi();
3 
4 void ClusterPlacer(double x, double y, double z, double GlobalZOffset , int ObjectNo, TGeoVolume *top,TGeoVolume *Clusterr,double phi, double theta, double psi)
5 {
6 TGeoRotation *Clusterrot = new TGeoRotation("Clusterrot",phi, theta,psi-phi);
7  Clusterrot->RegisterYourself();
8 top->AddNode(Clusterr,ObjectNo,new TGeoCombiTrans(x,y,GlobalZOffset+z,Clusterrot));
9 }
10 
11 int hypGe_TripleGeo()
12 {
13 
14 gROOT->Macro("$VMCWORKDIR/gconfig/rootlogon.C");
15 
16  // Load this libraries
17  gSystem->Load("libGeoBase");
18  gSystem->Load("libParBase");
19  gSystem->Load("libBase");
20  gSystem->Load("libPndData");
21  gSystem->Load("libPassive");
22  TString outfile= "../../../../geometry/hypGe_GeoMarcell.root";
23 
24  TFile* fi = new TFile(outfile,"RECREATE");
25 
26 
27 
28  FairGeoLoader* geoLoad = new FairGeoLoader("TGeo","FairGeoLoader");
29  FairGeoInterface *geoFace = geoLoad->getGeoInterface();
30  geoFace->setMediaFile("../../../../geometry/media_pnd.geo");
31  geoFace->readMedia();
32  geoFace->print();
33 
34  FairGeoMedia *Media = geoFace->getMedia();
35  FairGeoBuilder *geobuild=geoLoad->getGeoBuilder();
36 
37  FairGeoMedium *medGe = Media->getMedium("germanium");
38  Int_t nmedGe=geobuild->createMedium(medGe);
39 
40  // aluminum capsule
41  FairGeoMedium *medcap = Media->getMedium("HYPaluminium");
42  Int_t nmedcap=geobuild->createMedium(medcap);
43 
44  TGeoMedium *Ge = gGeoManager->GetMedium("germanium");
45 
46  //capsule
47  TGeoMedium *Al_sol = gGeoManager->GetMedium("HYPaluminium");
48 
49 
50 /*TGeoElementTable *table = gGeoManager->GetElementTable();
51 TGeoElement *Alum = table->GetElement(13);
52 TGeoElement *Germ = table->GetElement(32);
53  //Materials
54  TGeoMaterial *mat = new TGeoMaterial("Vacuum",0,0,0);
55  TGeoMaterial *Ger = new TGeoMaterial("Ger",Germ,5.3);
56  TGeoMaterial *Alu = new TGeoMaterial("Alu",Alum,2.7);
57  TGeoMaterial *Alusol = new TGeoMaterial("Alu",Alum,2.7);
58  //Media
59  TGeoMedium *vac = new TGeoMedium("Vacuum",0,mat);
60  TGeoMedium *Ge = new TGeoMedium("Ge",1,Ger);
61  TGeoMedium *Al = new TGeoMedium("Al",2,Alu);
62  TGeoMedium *Al_sol = new TGeoMedium("Alsol",2,Alusol);
63 
64 //Germ->Print();*/
65 
66  TGeoManager *geom = (TGeoManager*)gROOT->FindObject("FAIRGeom");
67  TGeoVolume *top = new TGeoVolumeAssembly("hpGe");
68  TGeoVolume *strange = new TGeoVolumeAssembly("strange");
69 
70 cout<<" geom "<<geom<<endl;
71  geom->SetTopVolume(top);
72 
73  Double_t GlobalZOffset = -55;
74  int ObjektNr = 0;
75  //Strahlrohr
76  TGeoVolume *Strahlrohr =geom->MakeTube("Strahlrohr",Al_sol,0,7.5,150);
77  Strahlrohr->SetLineColor(kBlue);
78  //top->AddNode(Strahlrohr,ObjektNr);
79  ObjektNr++;
80  //Außenhülle
81  TGeoVolume *Außenhülle =geom->MakeTube("Außenhülle",Al_sol,45,45.5,100);
82  Außenhülle->SetLineColor(kBlue);
83  //top->AddNode(Außenhülle,ObjektNr,new TGeoTranslation(0,0,GlobalZOffset));
84  ObjektNr++;
85 
86 //Halteplatte
87  TGeoVolume *HoldingStruc =geom->MakeBox("HoldingStruc",Al_sol,1,45,50);
88  HoldingStruc->SetLineColor(kGreen);
89  //top->AddNode(HoldingStruc,ObjektNr,new TGeoTranslation(0,0,GlobalZOffset));
90  ObjektNr++;
91 //Hilfskugel
92  TGeoSphere *Kugel = new TGeoSphere("Kugel",1,30,0,180,0,360);
93  TGeoVolume *Kugell = new TGeoVolume("Kugell",Kugel,Al_sol);
94  Kugell->SetLineColor(kBlue);
95  Kugell->SetTransparency(50);
96  //top->AddNode(Kugell,ObjektNr,new TGeoTranslation(0,0,GlobalZOffset));
97  ObjektNr++;
98 
99 //Nullpunktmarker
100  TGeoSphere *Origin = new TGeoSphere("Origin",0,1,0,180,0,360);
101  TGeoVolume *Originn = new TGeoVolume("Originn",Origin,Al_sol);
102  Originn->SetLineColor(kGreen);
103  //top->AddNode(Originn,ObjektNr);
104  ObjektNr++;
105 //Kappe
106  //Außen
107  double orangexy[16]={5.1211,-3.9906 , 2.5436,-7.0666 , -1.40290,-6.3389 , -2.5693,-3.0401 , 6.5119,-5.1633 , 3.2569,-9.0480 , -1.7273,-8.1289 , -3.2443,-3.8387};
108  double gelbxy[16]={-2.5693,-3.0401 , -6.0164,-2.4697 , -7.3916,1.3304 , -4.7881,4.3843 , -3.2443,-3.8387 , -7.7274,-3.0578 , -9.4641,1.7035 , -6.1761,5.5903};
109  double braunxy[16]={5.1211,-3.9906 , -2.5693,-3.0401 , -4.7881,4.3843 , -1.3481,3.7451 , 6.5119,-5.1633 , -3.2443,-3.8387 , -6.1761,5.5603 , -1.7021,4.7290};
110  double grauxy[16]={5.1211,-3.9906 , -1.3481,3.7451 , 0.8954,6.4301 , 3.9175,-0.7051 , 6.5119,-5.1633 , -1.7021,4.7290 , 1.2156,8.2210 , 4.9466,-0.8905};
111  double gruenxy[16]={0.8954,6.4301 , 4.8480,5.7361 , 6.1610,1.9544 , 3.9175,-0.7051 , 1.2156,8.2210 , 6.2074,7.3444 , 7.9036,2.5685 , 4.9466,-0.8905};
112 
113  TGeoArb8 *orange = new TGeoArb8 ("orange",12.95/2, orangexy );
114  TGeoArb8 *gelb = new TGeoArb8 ("gelb",12.95/2, gelbxy );
115  TGeoArb8 *braun = new TGeoArb8 ("braun",12.95/2, braunxy );
116  TGeoArb8 *grau = new TGeoArb8 ("grau",12.95/2, grauxy );
117  TGeoArb8 *gruen = new TGeoArb8 ("gruen",12.95/2, gruenxy );
118 
119  //Innen
120  double orangeInxy[16]={-2.4572,-2.9075 , -1.2905,-6.2072 , 2.4849,-6.9033 , 4.9505,-3.9607 , -3.1322,-3.7060 , -1.6148,-7.9972 , 3.1981,-8.8846 , 6.3412,-5.1334};
121  double gelbInxy[16]={-2.4572,-2.9075 , -4.7303,4.2211 , -7.2208,1.2997 , -5.9053,-2.3069 , -3.1322,-3.7060 , -6.1183,5.3971 , -9.2933,1.6727 , -7.6162,-2.9250};
122  double braunInxy[16]={-2.4572,-2.9075 , 4.9505,-3.9607 , -1.2893,3.5817 , -4.7303,4.2211 , -3.1322,-3.7060 , 6.3412,-5.1334 , -1.6433,4.5655 , -6.1183,5.3971};
123  double grauInxy[16]={4.9505,-3.9607 , 3.7466,-0.4744 , 0.9549,6.2675 , -1.2893,3.5817 , 6.3412,-5.1334 , 4.7756,-0.8597 , 1.2751,8.0583 , -1.6433,4.5655};
124  double gruenInxy[16]={3.7466,-0.4744 , 6.0208,1.9859 , 4.7360,5.6035 , 0.9549,6.2675 , 4.7756,-0.8597 , 7.7332,2.6000 , 6.0953,7.2118 , 1.2751,8.0583};
125 
126  TGeoTranslation *In_Z = new TGeoTranslation("In_Z",0,0,0.15);
127  In_Z->RegisterYourself();
128  TGeoArb8 *orangeIn = new TGeoArb8 ("orangeIn",(12.95-0.15)/2, orangeInxy );
129  TGeoArb8 *gelbIn = new TGeoArb8 ("gelbIn",(12.95-0.15)/2, gelbInxy );
130  TGeoArb8 *braunIn = new TGeoArb8 ("braunIn",(12.95-0.15)/2, braunInxy );
131  TGeoArb8 *grauIn = new TGeoArb8 ("grauIn",(12.95-0.15)/2, grauInxy );
132  TGeoArb8 *gruenIn = new TGeoArb8 ("gruenIn",(12.95-0.15)/2, gruenInxy );
133 
134 
135 //UpperPart
136  double up1_2xy[16]={6.5119,-5.1633 , 3.2569,-9.0480 , 3.1981,-8.8846 , 6.3412,-5.1334 ,
137  6.9097,-5.2329 , 3.3939,-9.4288 , 3.1981,-8.8846 , 6.3412,-5.1334};
138  double up2_3xy[16]={3.2569,-9.0480 , -1.7273,-8.1289 , -1.6148,-7.9972 , 3.1981,-8.8846 ,
139  3.3939,-9.4288 , -1.9896,-8.4361 , -1.6148,-7.9972 , 3.1981,-8.8846};
140  double up3_4xy[16]={-1.7273,-8.1289 , -3.2443,-3.8387 , -3.1322,-3.7060 , -1.6148,-7.9972 ,
141 -1.9896,-8.4361 , -3.5050,-4.1481 , -3.1322,-3.7060 , -1.6148,-7.9972};
142  double up4_5xy[16]={-3.2443,-3.8387 , -7.7274,-3.0578 , -7.6162,-2.9250 ,
143 -3.1322,-3.7060 ,
144 -3.5050,-4.1481 , -7.9866,-3.3866 , -7.6162,-2.9250 , -3.1322,-3.7060};
145  double up5_6xy[16]={-7.7274,-3.0578 , -9.4641,1.7035 , -9.2933,1.6727 , -7.6162,-2.9250 , -7.9866,-3.3866 , -9.8624,1.7752 , -9.2933,1.6727 ,-7.6162,-2.9250 };
146  double up6_7xy[16]={-9.4641,1.7035 , -6.1761,5.5603 , -6.1183,5.3971 , -9.2933,1.6727 , -9.8624,1.7752 , -6.3111,5.9410 , -6.1183,5.3971 , -9.2933,1.6727};
147  double up7_8xy[16]={-6.1761,5.5603 , -1.7021,4.7290 , -1.6433,4.5655 , -6.1183,5.3971 , -6.3111,5.9410 , -1.8393,5.1101 , -1.6433,4.5655 , -6.1183,5.3971};
148  double up8_9xy[16]={-1.7021,4.7290 , 1.2156,8.2210 , 1.2751,8.0583 , -1.6433,4.5655 , -1.8393,5.1101 , 1.0770,8.6003 , 1.2751,8.0583 , -1.6433,4.5655};
149  double up9_10xy[16]={1.2156,8.2210 , 6.2074,7.3444 , 6.0953,7.2118 , 1.2751,8.0583 , 1.0770,8.6003 , 6.4686,7.6535 , 6.0953,7.2118 , 1.2751,8.0583};
150  double up10_11xy[16]={6.2074,7.3444 , 7.9036,2.5685 , 7.7332,2.6000 , 6.0953,7.2118 , 6.4686,7.6535 , 8.3007,2.4950 , 7.7332,2.6000 , 6.0953,7.2118};
151  double up11_12xy[16]={7.9036,2.5685 , 4.9466,-0.8905 , 4.7756,-0.8597 , 7.7332,2.6000 , 8.3007,2.4950 , 5.3453,-0.9622 , 4.7756,-0.8597 , 7.7332,2.6000};
152 
153  double up12_1xy[16]={4.9466,-0.8905 , 6.5119,-5.1633 , 6.3412,-5.1334 , 4.7756,-0.8597 , 5.3453,-0.9622 , 6.9097,-5.2329 ,6.3412,-5.1334 , 4.7756,-0.8597};
154  TGeoArb8 *up1_2 = new TGeoArb8 ("up1_2",2.5, up1_2xy );
155  TGeoArb8 *up2_3 = new TGeoArb8 ("up2_3",2.5, up2_3xy );
156  TGeoArb8 *up3_4 = new TGeoArb8 ("up3_4",2.5, up3_4xy );
157  TGeoArb8 *up4_5 = new TGeoArb8 ("up4_5",2.5, up4_5xy );
158  TGeoArb8 *up5_6 = new TGeoArb8 ("up5_6",2.5, up5_6xy );
159  TGeoArb8 *up6_7 = new TGeoArb8 ("up6_7",2.5, up6_7xy );
160  TGeoArb8 *up7_8 = new TGeoArb8 ("up7_8",2.5, up7_8xy );
161  TGeoArb8 *up8_9 = new TGeoArb8 ("up8_9",2.5, up8_9xy );
162  TGeoArb8 *up9_10 = new TGeoArb8 ("up9_10",2.5, up9_10xy );
163  TGeoArb8 *up10_11 = new TGeoArb8 ("up10_11",2.5, up10_11xy );
164  TGeoArb8 *up11_12 = new TGeoArb8 ("up11_12",2.5, up11_12xy );
165  TGeoArb8 *up12_1 = new TGeoArb8 ("up12_1",2.5, up12_1xy );
166  TGeoTranslation *up_Z = new TGeoTranslation("up_Z",0,0,12.95/2+5./2);
167  up_Z->RegisterYourself();
168 
169 
170 
171  //Deckel
172 double xDeckel[12] = {6.9097,3.3939,-1.9896,-3.5058,-7.9866,-9.8624,-6.3111,-1.8393,1.0770,6.4686,8.3007,5.3453};
173 
174  double yDeckel[12] ={-5.2329,-9.4288,-8.4361,-4.1481,-3.3866,1.7752,5.9410,5.1101,8.6003,7.6535,2.4950,-0.9622};
175 
176 TGeoXtru *Deckel = new TGeoXtru (2);
177  Deckel->DefinePolygon(12,xDeckel,yDeckel);
178  Deckel->DefineSection(0,12.95/2+5,0,0,1);
179  Deckel->DefineSection(1,12.95/2+5+0.5,0,0,1);
180 TGeoVolume *DeckelVol = new TGeoVolume("Deckel",Deckel,Al_sol);
181  TGeoCompositeShape *KappenForm = new TGeoCompositeShape("KappenForm","(orange+gelb+braun+grau+gruen)-orangeIn:In_Z-gelbIn:In_Z-braunIn:In_Z-grauIn:In_Z-gruenIn:In_Z+(up1_2:up_Z+up2_3:up_Z+up3_4:up_Z+up4_5:up_Z+up5_6:up_Z+up6_7:up_Z+up7_8:up_Z+up8_9:up_Z+up9_10:up_Z+up10_11:up_Z+up11_12:up_Z+up12_1:up_Z)");
182  TGeoVolume *Kappe = new TGeoVolume("Kappe",KappenForm,Al_sol);
183  //Kappe->SetTransparency(50);
184  //DeckelVol->SetTransparency(50);
185  TGeoRotation *rightPosRot = new TGeoRotation("rightPosRot", 0,0,180+130.2);
186  rightPosRot->RegisterYourself();
187  TGeoCombiTrans *rightPos = new TGeoCombiTrans(0,0,6.475,rightPosRot);
188  rightPos->RegisterYourself();
189  //top->AddNode(Kappe,0,rightPos);
190  //top->AddNode(DeckelVol,1,rightPos);
191 
192  //TGeoVolume *Kappe2 = new TGeoVolume("Kappe2",KappenForm,Al_sol);
193  Kappe->SetLineColor(kBlue);
194  DeckelVol->SetLineColor(kBlue);
195  //top->AddNode(Kappe2,2);
196 
197  //Kapsel
198  TGeoTube *CapTube = new TGeoTube("CapTube",0,3.75,5);
199  TGeoPgon *CapPgon = new TGeoPgon("CapPgon",0.,360.,6,2);
200  CapPgon->DefineSection(0,-5,3.0571,3.8);
201  CapPgon->DefineSection(1,4.9094,3.75,3.8);
202 double CapThickness = 0.08;
203 double Scale = (3.0571-CapThickness)/3.0571;
204  TGeoTube *CapInTube = new TGeoTube("CapInTube",0,3.75*Scale,5);
205  TGeoPgon *CapInPgon = new TGeoPgon("CapInPgon",0.,360.,6,2);
206  CapInPgon->DefineSection(0,-5,3.0571-CapThickness,3.8);
207  CapInPgon->DefineSection(1,4.9094,3.75-CapThickness,3.8);
208  TGeoTube *CapTop = new TGeoTube("CapTop",0,3.75,CapThickness/2);
209  TGeoTranslation *CapTop_Z = new TGeoTranslation("CapTop_Z",0,0,5-CapThickness/2);
210  CapTop_Z->RegisterYourself();
211  TGeoCompositeShape *CapOutForm= new TGeoCompositeShape("CapOutForm","(CapTube-CapPgon)");
212  TGeoCompositeShape *CapInForm= new TGeoCompositeShape("CapInForm","(CapInTube-CapInPgon)");
213  TGeoTranslation *Cap_Z = new TGeoTranslation("Cap_Z",0,0,CapThickness);
214  Cap_Z->RegisterYourself();
215  TGeoCompositeShape *CapForm = new TGeoCompositeShape("CapForm","CapOutForm-CapInForm:Cap_Z+CapTop:CapTop_Z");
216 
217  TGeoVolume *Capsule = new TGeoVolume("Capsule",CapForm,Al_sol);
218  Capsule->SetTransparency(10);
219  Capsule->SetLineColor(10);
220  //top->AddNode(Capsule,5);
221 
222 //Kristall
223 double VacThickness = 0.08;
224 double Scale2 =(3.0571-CapThickness-VacThickness)/3.0571;
225 cout << Scale2 << endl;
226  TGeoTube *CryTube = new TGeoTube("CryTube",0,3.75,3.5);
227  TGeoPgon *CryPgon = new TGeoPgon("CryPgon",0.,360.,6,2);
228  CryPgon->DefineSection(0,-3.5,3.0571-CapThickness-VacThickness,3.8);
229  CryPgon->DefineSection(1,1.5+4.9094,3.75-CapThickness-VacThickness,3.8);
230  TGeoTube *CryHole = new TGeoTube("CryHole",0,0.5,3);
231  TGeoTranslation *CryHoleTrans = new TGeoTranslation("CryHoleTrans",0,0,0.5);
232 CryHoleTrans->RegisterYourself();
233 TGeoCompositeShape *CryForm= new TGeoCompositeShape("CryForm","((CryTube-CryPgon)-CryHole:CryHoleTrans)");
234 
235 TGeoVolume *Crystal = new TGeoVolume("Crystal",CryForm,Ge);
236  Crystal->SetLineColor(kGreen);
237  Crystal->SetFillColor(kGreen);
238 TGeoTranslation *CryTrans = new TGeoTranslation(0,0,5-1.5+CapThickness+VacThickness);
239 CryTrans->RegisterYourself();
240 //top->AddNode(Crystal,6,CryTrans); //top->AddNode(Kappe,0,rightPos);
241  //top->AddNode(DeckelVol,1,rightPos);
242 
243 
244 // Capsule+Crystal
245 TGeoVolume *CryCap[48];
246 for(int i = 0; i < 48;i++)
247 {
248  CryCap[i]= new TGeoVolumeAssembly("CryCap");
249  CryCap[i]->AddNode(Crystal,i+1,CryTrans);
250  CryCap[i]->AddNode(Capsule,100+i+1,new TGeoTranslation(0,0,5));
251 }
252 double TripleXOffset = -3.53/2-0.2;
253 double TripleYOffset = -6.1386/2;
254 TGeoRotation *CryCap1rot = new TGeoRotation("CryCap1rot", -0,+4,0);
255  CryCap1rot->RegisterYourself();
256  TGeoCombiTrans *CryCap1Combitrans = new TGeoCombiTrans(0+TripleXOffset,0+TripleYOffset,+0.45+0.2843,CryCap1rot);
257  CryCap1Combitrans->RegisterYourself();
258 
259 TGeoRotation *CryCap2rot = new TGeoRotation("CryCap2rot", -90,-6.93,90);
260  CryCap2rot->RegisterYourself();
261  TGeoCombiTrans *CryCap2Combitrans = new TGeoCombiTrans(5.3162+TripleXOffset,6.1386/2+TripleYOffset,0.45+0.2843,CryCap2rot);
262  CryCap2Combitrans->RegisterYourself();
263 
264 TGeoRotation *CryCap3rot = new TGeoRotation("CryCap3rot", 0,-4,0);
265  CryCap3rot->RegisterYourself();
266  TGeoCombiTrans *CryCap3Combitrans = new TGeoCombiTrans(0+TripleXOffset,6.1386+TripleYOffset,0.45+0.2843,CryCap3rot);
267  CryCap3Combitrans->RegisterYourself();
268 
269 cout << "test0" << endl;
270 
271 //Building 1 ClusterTemp
272 TGeoVolume *Pointingtube = geom->MakeTube("Pointer",Al_sol,0,0.2,150);
273 
274 TGeoVolume *Cluster[16];
275 TGeoVolume *ClusterTemp[16];
276 for(int i = 0;i < 16; i ++)
277 {
278  char buffer[20];
279  sprintf(buffer, "Cluster_%d",i);
280  Cluster[i]= new TGeoVolumeAssembly(buffer);
281  ClusterTemp[i]= new TGeoVolumeAssembly("ClusterTemp");
282 // ClusterTemp->AddNode(Pointingtube,5);
283  ClusterTemp[i]->AddNode(Kappe,200+i+1,rightPos);
284  ClusterTemp[i]->AddNode(DeckelVol,300+i+1,rightPos);
285  ClusterTemp[i]->AddNode(CryCap[3*i],1000+i+1,CryCap1Combitrans);
286 
287  ClusterTemp[i]->AddNode(CryCap[3*i+1],2000+i+1,CryCap2Combitrans);
288  ClusterTemp[i]->AddNode(CryCap[3*i+2],3000+i+1,CryCap3Combitrans);
289  cout <<"winkel "<< 180/Pi*TMath::ATan(12.5751/-26.3009) << endl;
290  cout <<"winkel2 "<< 180/Pi*TMath::ATan(7.0801/-26.3009) << endl;
291  TGeoRotation *ClusterTemprot = new TGeoRotation;
292  ClusterTemprot-> RotateZ(0);
293  ClusterTemprot->RotateY(180);
294  ClusterTemprot->RotateX(0);
295  ClusterTemprot->RegisterYourself();
296 
297 Cluster[i]->AddNode(ClusterTemp[i],400+i,new TGeoCombiTrans(0,0,0,ClusterTemprot));
298 }
299 //top->AddNode(ClusterTemp,100,new TGeoCombiTrans(-7.5713,12.0363,GlobalZOffset-26.416,ClusterTemprot2));
300 
301 cout << "test1" << endl;
302 
303 /*TGeoRotation *Clusterrot = new TGeoRotation;
304  Clusterrot-> RotateZ(0);
305  Clusterrot->RotateX(-180/Pi*TMath::ATan(12.5751/-26.3009));
306  Clusterrot->RotateY(180/Pi*TMath::ATan(7.0801/-26.3009));
307  Clusterrot->RegisterYourself();
308 top->AddNode(Cluster,1000,new TGeoCombiTrans(7.0801,12.5751,GlobalZOffset-26.3009,Clusterrot));*/
309 
310 //Inner Ring
311 ClusterPlacer(7.0801,12.5751, -26.3009, GlobalZOffset , 1000, strange,Cluster[0],180+150.62, 28.73,180-4.95);
312 cout << "test2" << endl;
313 ClusterPlacer(12.8422,-0, -27.1123, GlobalZOffset , 2000, strange,Cluster[1],180+90,25.34,180);
314  ClusterPlacer(7.0801,-12.5751, -26.3009, GlobalZOffset , 3000, strange,Cluster[2],180+29.38,28.73,180+4.95);
315 
316  ClusterPlacer(-7.0801,-12.5751, -26.3009, GlobalZOffset , 4000, strange,Cluster[3],150.62,28.73,4.95);
317  ClusterPlacer(-12.8422,0, -27.1123, GlobalZOffset , 5000, strange,Cluster[4],90,25.34,0);
318  ClusterPlacer(-7.0801,12.5751, -26.3009, GlobalZOffset , 6000, strange,Cluster[5],29.38,28.73,4.95);
319 
320 //Outer Ring
321 ClusterPlacer(8.421,20.8241, -19.8857, GlobalZOffset , 7000, strange,Cluster[6],180+157.98,48.5,0); //Patient
322 ClusterPlacer(17.9414,10.7963, -21.4836, GlobalZOffset , 8000, strange,Cluster[7],180+121.04,44.29,180-4.95);
323 ClusterPlacer(24.3173,0, -17.5689, GlobalZOffset , 9000, strange,Cluster[8],180+90,54.51,0);
324 ClusterPlacer(17.9414,-10.7963, -21.4836, GlobalZOffset , 10000, strange,Cluster[9],180+180-121.04,44.29,180+4.95);
325 ClusterPlacer(8.421,-20.8241, -19.8857, GlobalZOffset , 11000, strange,Cluster[10],180+180-157.98,48.5,0); //Patient
326 
327 ClusterPlacer(-8.421,-20.8241, -19.8857, GlobalZOffset , 12000, strange,Cluster[11],+157.98,48.5,180); //Patient
328 ClusterPlacer(-17.9414,-10.7963, -21.4836, GlobalZOffset , 13000, strange,Cluster[12],+121.04,44.29,-4.95);
329 ClusterPlacer(-24.3173,0, -17.5689, GlobalZOffset , 14000, strange,Cluster[13],+90,54.51,180);
330 ClusterPlacer(-17.9414,10.7963, -21.4836, GlobalZOffset , 15000, strange,Cluster[14],+180-121.04,44.29,4.95);
331 ClusterPlacer(-8.421,20.8241, -19.8857, GlobalZOffset , 16000, strange,Cluster[15],+180-157.98,48.5,180); //Patient
332 
333  top->AddNode(strange,0);
334  geom->CloseGeometry();
335  top->Write();
336  fi->Close();
337  //strange->Draw(""); strange->Raytrace();
338  strange->Draw("ogl");
339  geom->CheckOverlaps(0.0000000001,option="d"); //1 Overlap von Strahlrohr + Kugel ist normal!!!!!!!!!!!!!!!!!
340  geom->PrintOverlaps();
341  cout <<"Anzahl der Objekte: "<< ObjektNr << endl;
342  //Paket[1]->Print();
343  return 0;
344 }
345 
346 
theta
double theta
Definition: macro/detectors/emc/dedicated/reco_analys.C:24
geoLoad
FairGeoLoader * geoLoad
Definition: createRootGeoFile.C:16
i
Int_t i
Definition: run_full.C:25
Media
FairGeoMedia * Media
Definition: createRootFscGeometryFile.C:68
gGeoManager
TGeoManager * gGeoManager
Definition: createRootGeoFile.C:97
top
TGeoVolume * top
Definition: createRootFscGeometryFile.C:89
RotateZ
rotCave RotateZ(90)
geom
void geom()
Definition: createRootGeoFileBarrel_2018v1.C:79
geobuild
FairGeoBuilder * geobuild
Definition: createRootFscGeometryFile.C:69
fi
TFile * fi
Definition: createRootGeoFile.C:238
Double_t
Double_t
Definition: SimCompleteLinkDef.h:6
z
Double_t z
Definition: createRootFscGeometryFile.C:146
phi
double phi
Definition: macro/detectors/emc/dedicated/reco_analys.C:24
Scale
f1 Scale(1/numero)
hypGe_TripleGeo
int hypGe_TripleGeo()
Definition: hypGe_TripleGeo.C:11
x
Double_t x
Definition: createRootFscGeometryFile.C:159
ClusterPlacer
void ClusterPlacer(double x, double y, double z, double GlobalZOffset, int ObjectNo, TGeoVolume *top, TGeoVolume *Clusterr, double phi, double theta, double psi)
Definition: hypGe_DoubleGeo.C:4
geoFace
FairGeoInterface * geoFace
Definition: createRootGeoFile.C:17
y
Double_t y
Definition: createRootFscGeometryFile.C:159
Pi
Double_t Pi
Definition: PndGeoHypGeCrystal.cxx:10
outfile
TString outfile
Definition: createRootFscGeometryFile.C:57
Generated on Wed Apr 3 2019 10:02:47 for FairRoot/PandaRoot by   doxygen 1.8.5