Minors
[u/mrichter/AliRoot.git] / ITS / Displayv11.C
1 //----------------------------------------------------------------------
2 void Displayv11(const char* filename=""){
3     // Display AliITSv11 Geometry
4     // Inputs:
5     //    const char* filename output file with the display in it
6     // Outputs:
7     //    none.
8     // Retrurn:
9     //    none.
10
11     gSystem->Load("libGeom");
12     //
13     if(gGeoManager) delete gGeoManager;
14     mgr2 = gGeoManager = new TGeoManager("ITSGeometry",
15                                          " ITS Simulation Geometry Manager");
16     //
17     TGeoMaterial *vacmat = new TGeoMaterial("Vacume",0,0,0);
18     TGeoMedium   *vacmed = new TGeoMedium("Vacume_med",1,vacmat);
19     TGeoVolume *ALIC = mgr2->MakeBox("ALIC",vacmed,100.,100.,200.);
20     mgr2->SetTopVolume(ALIC);
21     TGeoVolume *ITS = mgr2->MakeBox("ITSV",vacmed,99.,99.,199.);
22     ALIC->AddNode(ITS,1);
23     //
24     /*
25     AliITSv11 *its = new AliITSv11(0,3);
26     its->SetDebug(ISetits(0,-1));
27     its->GetSPDGeometry()->SetDebug(ISetits(0,-1));
28     its->GetSupGeometry()->SetDebug(ISetits(0,-1));
29     its->CreateMaterials();
30     its->CreateGeometry();
31     */
32     AliITSv11GeometrySPD *gspd = new AliITSv11GeometrySPD();
33     //AliITSv11GeometrySDD *gsdd = new AliITSgeometrySDD();
34     //AliITSv11GeometrySupport *gsupp = new AliITSGeometrySupport();
35     //AliITSv11GeometrySSD *gssd = new AliITSGeometrySSD();
36     //
37     Int_t imat=1,imed=1;
38     gspd->CreateSPDCenteralMaterials(imed,imat);
39     gspd->SPDSector(ITS,mgr2);
40     //
41     mgr2->CloseGeometry();
42     //
43     TControlBar *bar=new TControlBar("vertical","ITS Geometry Display",10,10);
44     bar->AddButton("Set Clipping on","ISetits(2,1)","Clipping on");
45     bar->AddButton("Set Cllipping off","ISetits(2,0)","Clipping off");
46     bar->AddButton("Set axis on","ISetits(3,1)","Show Axis on");
47     bar->AddButton("Set axis off","ISetits(3,0)","Show Axis off");
48     bar->AddButton("Set perspective on","ISetits(4,1)","Perspective on");
49     bar->AddButton("Set perspective off","ISetits(4,0)","Perspective off");
50     bar->AddButton("Set RayTrace on","ISetits(5,1)","Perspective on");
51     bar->AddButton("Set RayTrace off","ISetits(5,0)","Perspective off");
52     bar->AddButton("Set circle/80","ISetits(1,80)","circles ~ by 80 lines");
53     bar->AddButton("Display Geometry","Displayit()","Run Displayit");
54     bar->AddButton("Display SPD Ceneral Volume","EngineeringSPDLayer()",
55                    "Run EngineeringSPDLayer");
56     bar->AddButton("Display SPD Thermal Sheald","EngineeringSPDThS()",
57                    "Run EngineeringSPDThS");
58     bar->AddButton("Display SDD Layer 3","EngineeringSDDLayer3()",
59                    "Run EngineeringSDDLayer3");
60     bar->AddButton("Display SDD Layer 4","EngineeringSDDLayer4()",
61                    "Run EngineeringSDDLayer4");
62     bar->AddButton("Display SDD Cone","EngineeringSDDCone()",
63                    "Run EngineeringSDDCone");
64     bar->AddButton("Display SDD Centeral Cylinder","EngineeringSDDCylinder()",
65                    "Run EngineeringSDDCylinder");
66     bar->AddButton("Display SSD Cone","EngineeringSSDCone()",
67                    "Run EngineeringSSDCone");
68     bar->AddButton("Display SSD Centeral Cylinder","EngineeringSSDCylinder()",
69                    "Run EngineeringSSDCylinder");
70     bar->AddButton("Display SUP RB24 side","EngineeringSupRB24()",
71                    "Run EngineeringSupRB24");
72     bar->AddButton("Display Cable Trays RB24 side","EngineeringSupTrayRB24()",
73                    "Run EngineeringSupTrayRB24");
74     bar->AddButton("Display SUP RB26 side","EngineeringSupRB26()",
75                    "Run EngineeringSupRB26");
76     bar->AddButton("Save Geometry to File","ExportToFile()",
77                    "Run ExportToFile");
78     bar->AddButton("Quit/Exit",".q","Exit");
79     bar->Show();
80     gROOT->SaveContext();
81          //Displayit();
82 }
83 //----------------------------------------------------------------------
84 void ExportToFile(){
85     // Quirry file name and write geometry to a root file.
86     // Inputs:
87     //    const char* filename output file with the display in it
88     // Outputs:
89     //    none.
90     // Retrurn:
91     //    none.
92     Char_t filename[100];
93
94     printf("Eneter File name:");
95     scanf("%s",filename);
96
97     gGeoManager->Export(filename);
98 }
99 //----------------------------------------------------------------------
100 Int_t ISetits(Int_t t,Int_t v){
101     static Int_t itsdebug=1,nsegments=80,cut=0,axis=1,perspective=0,ray=0;
102
103     switch (t) {
104     case 0:
105         if(v<0) return itsdebug;
106         itsdebug = v;
107         break;
108     case 1:
109         if(v<0) return nsegments;
110         nsegments= v;
111         break;
112     case 2:
113         if(v<0) return cut;
114         cut = v;
115         break;
116     case 3:
117         if(v<0) return axis;
118         axis = v;
119         break;
120     case 4:
121         if(v<0) return perspective;
122         perspective = v;
123         break;
124     case 5:
125         if(v<0) return ray;
126         ray = v;
127         break;
128     }// end switch
129     return 0;
130 }
131 //----------------------------------------------------------------------
132 Double_t DSetits(Int_t t,Double_t v){
133     static Double_t phimincut=0.0,phimaxcut=180.0;
134     static Double_t longitude=90.0,latitude=0.0,psi=0.0;
135
136     switch (t) {
137     case 0:
138         if(v<0.) return phimincut;
139         phimincut = v;
140         break;
141     case 1:
142         if(v<0.) return phimaxcut;
143         phimaxcut = v;
144         break;
145     case 2:
146         if(v<0.) return longitude;
147         longitude = v;
148         break;
149     case 3:
150         if(v<0.) return latitude;
151         latitude = v;
152         break;
153     case 4:
154         if(v<0.) return latitude;
155         latitude = v;
156         break;
157     }// end switch
158     return 0;
159 }
160 //----------------------------------------------------------------------
161 void Displayit(){
162     // Display AliITSv11 Geometry
163     // Inputs:
164     //    const char* filename output file with the display in it
165     // Outputs:
166     //    none.
167     // Retrurn:
168     //    none.
169     Int_t irr;
170     //
171     TGeoManager *mgr2 = gGeoManager;
172     TGeoVolume *ALIC = mgr2->GetTopVolume();
173     TCanvas *c1;
174     if(!(c1 = (TCanvas*)gROOT->FindObject("C1")))
175         c1 = new TCanvas("C1","ITS Simulation Geometry",900,900);
176     c1->Divide(2,2);
177     //
178     mgr2->SetNsegments(ISetits(1,-1));
179     //
180     mgr2->SetVisLevel(6);
181     mgr2->SetVisOption(0);
182     //mgr2->CheckOverlaps(0.01);
183     //mgr2->PrintOverlaps();
184     if(ISetits(2,-1)==1){
185         TGeoShape *clip = new TGeoTubeSeg(0, 1000, 2000, 45, 90);
186         mgr2->SetClippingShape(clip);
187     } // end if
188     if(ISetits(2,-1)!=0) mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
189     //
190     c1->cd(1);
191     ALIC->Draw();
192     TPad *p1 = c1->GetPad(1);
193     TView *view1 = p1->GetView();
194     if(view1){
195         view1->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
196         if(irr) cout <<"error="<<irr<<endl;
197         if(ISetits(4,-1)==0) view1->SetParallel();
198         else  view1->SetPerspective();
199         view1->Front();
200         if(ISetits(3,-1)!=0) view1->ShowAxis();
201     } // end if view1
202     if(ISetits(5,-1)==1) ALIC->Raytrace();
203     c1->cd(2);
204     ALIC->Draw();
205     TPad *p2 = c1->GetPad(2);
206     TView *view2 = p2->GetView();
207     if(view2){
208         view2->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
209         if(irr) cout <<"error="<<irr<<endl;
210         if(ISetits(4,-1)==0) view2->SetParallel();
211         else  view2->SetPerspective();
212         view2->RotateView(60.,30.);
213         if(ISetits(3,-1)!=0) view2->ShowAxis();
214     } // end if view2
215     if(ISetits(5,-1)==1) ALIC->Raytrace();
216     c1->cd(3);
217     ALIC->Draw();
218     c1->SetPhi(90.0); c1->SetTheta(90.0);
219     TPad *p3 = c1->GetPad(3);
220     TView *view3 = p3->GetView();
221     if(view3){
222         view3->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
223         if(irr) cout <<"error="<<irr<<endl;
224         if(ISetits(4,-1)==0) view3->SetParallel();
225         else  view3->SetPerspective();
226         view3->Top();
227         if(ISetits(3,-1)!=0) view3->ShowAxis();
228     } // end if view3
229     if(ISetits(5,-1)==1) ALIC->Raytrace();
230     c1->cd(4);
231     ALIC->Draw();
232     TPad *p4 = c1->GetPad(4);
233     TView *view4 = p4->GetView();
234     if(view4){
235         view4->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
236         if(irr) cout <<"error="<<irr<<endl;
237         if(ISetits(4,-1)==0) view4->SetParallel();
238         else  view4->SetPerspective();
239         view4->Side();
240         if(ISetits(3,-1)!=0) view4->ShowAxis();
241     } // end if view4
242     if(ISetits(5,-1)==1) ALIC->Raytrace();
243     //
244 }
245 //----------------------------------------------------------------------
246 void EngineeringSPDLayer(){
247     // Display SPD Layer Geometry
248     // Inputs:
249     //    none.
250     // Outputs:
251     //    none.
252     // Retrurn:
253     //    none.
254     Int_t irr;
255     //
256     TGeoManager *mgr2 = gGeoManager;
257     TGeoVolume *ALIC = mgr2->GetTopVolume();
258     TCanvas *c4;
259     if(!(c4 = (TCanvas*)gROOT->FindObject("C4")))
260         c4 = new TCanvas("C4","ITS SPD Layer Geometry Side View",500,500);
261     TGeoVolume *ITS,*SPDLay=0;
262     TGeoNode *node;
263     TArrow *arrow=new TArrow();
264     //
265     node = ALIC->FindNode("ITSV_1");
266     ITS = node->GetVolume();
267     node = ITS->FindNode("ITSSPDCarbonFiberSectorV_1");
268     //node = ITS->FindNode("ITSSPDTempSPDMotherVolume_1");
269     SPDLay = node->GetVolume();
270     //
271     mgr2->SetNsegments(ISetits(1,-1));
272     //
273     mgr2->SetVisLevel(6);
274     mgr2->SetVisOption(0);
275     //mgr2->CheckOverlaps(0.01);
276     //mgr2->PrintOverlaps();
277     if(ISetits(2,-1)==1){
278         TGeoShape *clip = new TGeoTubeSeg(0, 1000, 2000, 45, 90);
279         mgr2->SetClippingShape(clip);
280     } // end if
281     mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
282     if(ISetits(2,-1)!=0) mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
283     //
284     SPDLay->Draw();
285     TView *view1 = c4->GetView();
286     if(view1){
287         view1->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
288         if(irr) cout <<"error="<<irr<<endl;
289         if(ISetits(4,-1)==0) view1->SetParralel();
290         else  view1->SetPerspective();
291         view1->Front();
292         if(ISetits(3,-1)!=0) view1->ShowAxis();
293     } // end if view1
294     if(ISetits(5,-1)==1) SPDLay->Raytrace();
295     //
296     if(!(c5 = (TCanvas*)gROOT->FindObject("C5")))
297         c5 = new TCanvas("C5","ITS SPD Layer Geometry End View",500,500);
298     SPDLay->Draw();
299     TView *view2 = c5->GetView();
300     if(view2){
301         view2->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
302         if(irr) cout <<"error="<<irr<<endl;
303         if(ISetits(4,-1)==0) view2->SetParralel();
304         else  view2->SetPerspective();
305         view2->Top();
306         if(ISetits(3,-1)!=0) view2->ShowAxis();
307     } // end if view2
308     if(ISetits(5,-1)==1) SPDLay->Raytrace();
309     //
310 }
311 //----------------------------------------------------------------------
312 void EngineeringSPDThS(){
313     // Display SPD Thermal Sheald Geometry
314     // Inputs:
315     //    none.
316     // Outputs:
317     //    none.
318     // Retrurn:
319     //    none.
320     Int_t irr;
321     //
322     TGeoManager *mgr2 = gGeoManager;
323     TGeoVolume *ALIC = mgr2->GetTopVolume();
324     TCanvas *c4;
325     if(!(c4 = (TCanvas*)gROOT->FindObject("C4")))
326         c4 = new TCanvas("C4","ITS SDD Cylinder Geometry",900,450);
327     c4->Divide(2,1);
328     TGeoVolume *ITS,*SPDThS=0;
329     TGeoNode *node;
330     TArrow *arrow=new TArrow();
331     //
332     node = ALIC->FindNode("ITSV_1");
333     ITS = node->GetVolume();
334     node = ITS->FindNode("ITSspdThermalShealdMotherM_1");
335     SPDThS = node->GetVolume();
336     //
337     mgr2->SetNsegments(ISetits(1,-1));
338     //
339     mgr2->SetVisLevel(6);
340     mgr2->SetVisOption(0);
341     //mgr2->CheckOverlaps(0.01);
342     //mgr2->PrintOverlaps();
343     if(ISetits(2,-1)==1){
344         TGeoShape *clip = new TGeoTubeSeg(0, 1000, 2000, 45, 90);
345         mgr2->SetClippingShape(clip);
346     } // end if
347     mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
348     if(ISetits(2,-1)!=0) mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
349     //
350     c4->cd(1);
351     SPDThS->Draw();
352     TPad *p1 = c4->GetPad(1);
353     TView *view1 = p1->GetView();
354     if(view1){
355         view1->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
356         if(irr) cout <<"error="<<irr<<endl;
357         if(ISetits(4,-1)==0) view1->SetParralel();
358         else  view1->SetPerspective();
359         view1->Front();
360         if(ISetits(3,-1)!=0) view1->ShowAxis();
361     } // end if view1
362     if(ISetits(5,-1)==1) SPDThS->Raytrace();
363     //
364     c4->cd(2);
365     SPDThS->Draw();
366     TPad *p2 = c4->GetPad(2);
367     TView *view2 = p2->GetView();
368     if(view2){
369         view2->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
370         if(irr) cout <<"error="<<irr<<endl;
371         if(ISetits(4,-1)==0) view2->SetParralel();
372         else  view2->SetPerspective();
373         view2->Top();
374         if(ISetits(3,-1)!=0) view2->ShowAxis();
375     } // end if view2
376     if(ISetits(5,-1)==1) SPDThS->Raytrace();
377     //
378 }
379 //----------------------------------------------------------------------
380 void EngineeringSDDLayer3(){
381     // Display SDD Layer 3 Geometry
382     // Inputs:
383     //    none.
384     // Outputs:
385     //    none.
386     // Retrurn:
387     //    none.
388     Int_t irr;
389     //
390     TGeoManager *mgr2 = gGeoManager;
391     TGeoVolume *ALIC = mgr2->GetTopVolume();
392     TCanvas *c4;
393     if(!(c4 = (TCanvas*)gROOT->FindObject("C4")))
394         c4 = new TCanvas("C4","ITS SDD Layer 3 Geometry Side View",500,500);
395     TGeoVolume *ITS,*SDDLay3=0;
396     TGeoNode *node;
397     TArrow *arrow=new TArrow();
398     //
399     node = ALIC->FindNode("ITSV_1");
400     ITS = node->GetVolume();
401     node = ITS->FindNode("ITSsddLayer3_1");
402     SDDLay3 = node->GetVolume();
403     //
404     mgr2->SetNsegments(ISetits(1,-1));
405     //
406     mgr2->SetVisLevel(6);
407     mgr2->SetVisOption(0);
408     //mgr2->CheckOverlaps(0.01);
409     //mgr2->PrintOverlaps();
410     if(ISetits(2,-1)==1){
411         TGeoShape *clip = new TGeoTubeSeg(0, 1000, 2000, 45, 90);
412         mgr2->SetClippingShape(clip);
413     } // end if
414     mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
415     if(ISetits(2,-1)!=0) mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
416     //
417     SDDLay3->Draw();
418     TView *view1 = c4->GetView();
419     if(view1){
420         view1->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
421         if(irr) cout <<"error="<<irr<<endl;
422         if(ISetits(4,-1)==0) view1->SetParralel();
423         else  view1->SetPerspective();
424         view1->Front();
425         if(ISetits(3,-1)!=0) view1->ShowAxis();
426     } // end if view1
427     if(ISetits(5,-1)==1) SDDLay3->Raytrace();
428     //
429     if(!(c5 = (TCanvas*)gROOT->FindObject("C5")))
430         c5 = new TCanvas("C5","ITS SDD Layer 3 Geometry End View",500,500);
431     SDDLay3->Draw();
432     TView *view2 = c5->GetView();
433     if(view2){
434         view2->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
435         if(irr) cout <<"error="<<irr<<endl;
436         if(ISetits(4,-1)==0) view2->SetParralel();
437         else  view2->SetPerspective();
438         view2->Top();
439         if(ISetits(3,-1)!=0) view2->ShowAxis();
440     } // end if view2
441     if(ISetits(5,-1)==1) SDDLay3->Raytrace();
442     //
443 }
444 //----------------------------------------------------------------------
445 void EngineeringSDDLayer4(){
446     // Display SDD Layer 4 Geometry
447     // Inputs:
448     //    none.
449     // Outputs:
450     //    none.
451     // Retrurn:
452     //    none.
453     Int_t irr;
454     //
455     TGeoManager *mgr2 = gGeoManager;
456     TGeoVolume *ALIC = mgr2->GetTopVolume();
457     TCanvas *c4;
458     if(!(c4 = (TCanvas*)gROOT->FindObject("C4")))
459         c4 = new TCanvas("C4","ITS SDD Layer 4 Geometry Side View",500,500);
460     TGeoVolume *ITS,*SDDLay4=0;
461     TGeoNode *node;
462     TArrow *arrow=new TArrow();
463     //
464     node = ALIC->FindNode("ITSV_1");
465     ITS = node->GetVolume();
466     node = ITS->FindNode("ITSsddLayer4_1");
467     SDDLay4 = node->GetVolume();
468     //
469     mgr2->SetNsegments(ISetits(1,-1));
470     //
471     mgr2->SetVisLevel(6);
472     mgr2->SetVisOption(0);
473     //mgr2->CheckOverlaps(0.01);
474     //mgr2->PrintOverlaps();
475     if(ISetits(2,-1)==1){
476         TGeoShape *clip = new TGeoTubeSeg(0, 1000, 2000, 45, 90);
477         mgr2->SetClippingShape(clip);
478     } // end if
479     mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
480     if(ISetits(2,-1)!=0) mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
481     //
482     SDDLay4->Draw();
483     TView *view1 = c4->GetView();
484     if(view1){
485         view1->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
486         if(irr) cout <<"error="<<irr<<endl;
487         if(ISetits(4,-1)==0) view1->SetParralel();
488         else  view1->SetPerspective();
489         view1->Front();
490         if(ISetits(3,-1)!=0) view1->ShowAxis();
491     } // end if view1
492     if(ISetits(5,-1)==1) SDDLay4->Raytrace();
493     //
494     if(!(c5 = (TCanvas*)gROOT->FindObject("C5")))
495         c5 = new TCanvas("C5","ITS SDD Layer 4 Geometry End View",500,500);
496     SDDLay4->Draw();
497     TView *view2 = c5->GetView();
498     if(view2){
499         view2->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
500         if(irr) cout <<"error="<<irr<<endl;
501         if(ISetits(4,-1)==0) view2->SetParralel();
502         else  view2->SetPerspective();
503         view2->Top();
504         if(ISetits(3,-1)!=0) view2->ShowAxis();
505     } // end if view2
506     if(ISetits(5,-1)==1) SDDLay4->Raytrace();
507     //
508 }
509 //----------------------------------------------------------------------
510 void EngineeringSDDCone(){
511     // Display SDD Cone Geometry
512     // Inputs:
513     //    none.
514     // Outputs:
515     //    none.
516     // Retrurn:
517     //    none.
518     Int_t irr;
519     //
520     TGeoManager *mgr2 = gGeoManager;
521     TGeoVolume *ALIC = mgr2->GetTopVolume();
522     TCanvas *c2;
523     if(!(c2 = (TCanvas*)gROOT->FindObject("C2")))
524         c2 = new TCanvas("C2","ITS SDD Cone Geometry",900,450);
525     c2->Divide(2,1);
526     TGeoVolume *ITS,*SDD=0;
527     TGeoNode *node;
528     //
529     node = ALIC->FindNode("ITSV_1");
530     ITS = node->GetVolume();
531     node = ITS->FindNode("ITSsddConeL_1");
532     SDD = node->GetVolume();
533     //
534     mgr2->SetNsegments(ISetits(1,-1));
535     //
536     mgr2->SetVisLevel(6);
537     mgr2->SetVisOption(0);
538     //mgr2->CheckOverlaps(0.01);
539     //mgr2->PrintOverlaps();
540     if(ISetits(2,-1)==1){
541         TGeoShape *clip = new TGeoTubeSeg(0, 1000, 2000, 45, 90);
542         mgr2->SetClippingShape(clip);
543     } // end if
544     mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
545     if(ISetits(2,-1)!=0) mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
546     //
547     c2->cd(1);
548     SDD->Draw();
549     TPad *p1 = c2->GetPad(1);
550     TView *view1 = p1->GetView();
551     if(view1){
552         view1->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
553         if(irr) cout <<"error="<<irr<<endl;
554         if(ISetits(4,-1)==0) view1->SetParralel();
555         else  view1->SetPerspective();
556         view1->Front();
557         if(ISetits(3,-1)!=0) view1->ShowAxis();
558     } // end if view1
559     if(ISetits(5,-1)==1) SDD->Raytrace();
560     //
561     c2->cd(2);
562     SDD->Draw();
563     TPad *p2 = c2->GetPad(2);
564     TView *view2 = p2->GetView();
565     if(view2){
566         view2->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
567         if(irr) cout <<"error="<<irr<<endl;
568         if(ISetits(4,-1)==0) view2->SetParralel();
569         else  view2->SetPerspective();
570         view2->Top();
571         if(ISetits(3,-1)!=0) view2->ShowAxis();
572     } // end if view1
573     if(ISetits(5,-1)==1) SDD->Raytrace();
574     //
575 }
576 //----------------------------------------------------------------------
577 void EngineeringSDDCylinder(){
578     // Display SDD Cylinder Geometry
579     // Inputs:
580     //    none.
581     // Outputs:
582     //    none.
583     // Retrurn:
584     //    none.
585     Int_t irr;
586     //
587     TGeoManager *mgr2 = gGeoManager;
588     TGeoVolume *ALIC = mgr2->GetTopVolume();
589     TCanvas *c3;
590     if(!(c3 = (TCanvas*)gROOT->FindObject("C3")))
591         c3 = new TCanvas("C3","ITS SDD Cylinder Geometry",900,450);
592     c3->Divide(2,1);
593     TGeoVolume *ITS,*SDD=0;
594     TGeoNode *node;
595     TArrow *arrow=new TArrow();
596     //
597     node = ALIC->FindNode("ITSV_1");
598     ITS = node->GetVolume();
599     node = ITS->FindNode("ITSsddCentCylCF_1");
600     SDD = node->GetVolume();
601     Double_t Rmin = ((TGeoTube*)(SDD->GetShape()))->GetRmin();
602     Double_t Rmax = ((TGeoTube*)(SDD->GetShape()))->GetRmax();
603     Double_t Dz   = ((TGeoTube*)(SDD->GetShape()))->GetDz();
604     //
605     mgr2->SetNsegments(ISetits(1,-1));
606     //
607     mgr2->SetVisLevel(6);
608     mgr2->SetVisOption(0);
609     //mgr2->CheckOverlaps(0.01);
610     //mgr2->PrintOverlaps();
611     if(ISetits(2,-1)==1){
612         TGeoShape *clip = new TGeoTubeSeg(0, 1000, 2000, 45, 90);
613         mgr2->SetClippingShape(clip);
614     } // end if
615     mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
616     if(ISetits(2,-1)!=0) mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
617     //
618     c3->cd(1);
619     SDD->Draw();
620     TPad *p1 = c3->GetPad(1);
621     TView *view1 = p1->GetView();
622     if(view1){
623         view1->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
624         if(irr) cout <<"error="<<irr<<endl;
625         if(ISetits(4,-1)==0) view1->SetParralel();
626         else  view1->SetPerspective();
627         view1->Front();
628         if(ISetits(3,-1)!=0) view1->ShowAxis();
629     } // end if view1
630     if(ISetits(5,-1)==1) SDD->Raytrace();
631     //arrow->DrawArrow(1.01*Rmax,-Dz,1.01*Rmax,+Dz);
632     //
633     c3->cd(2);
634     SDD->Draw();
635     TPad *p2 = c3->GetPad(2);
636     TView *view2 = p2->GetView();
637     if(view2){
638         view2->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
639         if(irr) cout <<"error="<<irr<<endl;
640         if(ISetits(4,-1)==0) view2->SetParralel();
641         else  view2->SetPerspective();
642         view2->Top();
643         if(ISetits(3,-1)!=0) view2->ShowAxis();
644     } // end if view2
645     if(ISetits(5,-1)==1) SDD->Raytrace();
646     //arrow->DrawArrow(Rmax,0.0,Rmax,0.0);
647     //Double_t s = TMath::Sin(0.7),c = TMath::Cos(0.7);
648     //arrow->DrawArrow(-Rmin*c,-Rmin*s,Rmin*c,Rmin*s);
649     //
650 }
651 //----------------------------------------------------------------------
652 void EngineeringSSDCone(){
653     // Display SSD Cone Geometry
654     // Inputs:
655     //    none.
656     // Outputs:
657     //    none.
658     // Retrurn:
659     //    none.
660     Int_t irr;
661     //
662     TGeoManager *mgr2 = gGeoManager;
663     TGeoVolume *ALIC = mgr2->GetTopVolume();
664     TCanvas *c2;
665     if(!(c2 = (TCanvas*)gROOT->FindObject("C2")))
666         c2 = new TCanvas("C2","ITS SSD Cone Geometry",900,450);
667     c2->Divide(2,1);
668     TGeoVolume *ITS,*SSD=0;
669     TGeoNode *node;
670     //
671     node = ALIC->FindNode("ITSV_1");
672     ITS = node->GetVolume();
673     node = ITS->FindNode("ITSssdConeA_1");
674     SSD = node->GetVolume();
675     //
676     mgr2->SetNsegments(ISetits(1,-1));
677     //
678     mgr2->SetVisLevel(6);
679     mgr2->SetVisOption(0);
680     //mgr2->CheckOverlaps(0.01);
681     //mgr2->PrintOverlaps();
682     if(ISetits(2,-1)==1){
683         TGeoShape *clip = new TGeoTubeSeg(0, 1000, 2000, 45, 90);
684         mgr2->SetClippingShape(clip);
685     } // end if
686     mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
687     if(ISetits(2,-1)!=0) mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
688     //
689     c2->cd(1);
690     SSD->Draw();
691     TPad *p1 = c2->GetPad(1);
692     TView *view1 = p1->GetView();
693     if(view1){
694         view1->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
695         if(irr) cout <<"error="<<irr<<endl;
696         if(ISetits(4,-1)==0) view1->SetParralel();
697         else  view1->SetPerspective();
698         view1->Top();
699         if(ISetits(3,-1)!=0) view1->ShowAxis();
700     } // end if view1
701     if(ISetits(5,-1)==1) SSD->Raytrace();
702     //
703     c2->cd(2);
704     SSD->Draw();
705     TPad *p2 = c2->GetPad(2);
706     TView *view2 = p2->GetView();
707     if(view2){
708         view2->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
709         if(irr) cout <<"error="<<irr<<endl;
710         if(ISetits(4,-1)==0) view2->SetParralel();
711         else  view2->SetPerspective();
712         view2->Front();
713         if(ISetits(3,-1)!=0) view2->ShowAxis();
714     } // end if view1
715     if(ISetits(5,-1)==1) SSD->Raytrace();
716     //
717 }
718 //----------------------------------------------------------------------
719 void EngineeringSSDCylinder(){
720     // Display SSD Cylinder Geometry
721     // Inputs:
722     //    none.
723     // Outputs:
724     //    none.
725     // Retrurn:
726     //    none.
727     Int_t irr;
728     //
729     TGeoManager *mgr2 = gGeoManager;
730     TGeoVolume *ALIC = mgr2->GetTopVolume();
731     TCanvas *c3;
732     if(!(c3 = (TCanvas*)gROOT->FindObject("C3")))
733         c3 = new TCanvas("C3","ITS SSD Cylinder Geometry",900,450);
734     c3->Divide(2,1);
735     TGeoVolume *ITS,*SSD=0;
736     TGeoNode *node;
737     TArrow *arrow=new TArrow();
738     //
739     node = ALIC->FindNode("ITSV_1");
740     ITS = node->GetVolume();
741     node = ITS->FindNode("ITSssdCentCylCA_1");
742     SSD = node->GetVolume();
743     //Double_t Rmin = ((TGeoPcon*)(SSD->GetShape()))->GetRmin();
744     //Double_t Rmax = ((TGeoPcon*)(SSD->GetShape()))->GetRmax();
745     //Double_t Dz   = ((TGeoPcon*)(SSD->GetShape()))->GetDz();
746     //
747     mgr2->SetNsegments(ISetits(1,-1));
748     //
749     mgr2->SetVisLevel(6);
750     mgr2->SetVisOption(0);
751     //mgr2->CheckOverlaps(0.01);
752     //mgr2->PrintOverlaps();
753     if(ISetits(2,-1)==1){
754         TGeoShape *clip = new TGeoTubeSeg(0, 1000, 2000, 45, 90);
755         mgr2->SetClippingShape(clip);
756     } // end if
757     mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
758     if(ISetits(2,-1)!=0) mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
759     //
760     c3->cd(1);
761     SSD->Draw();
762     TPad *p1 = c3->GetPad(1);
763     TView *view1 = p1->GetView();
764     if(view1){
765         view1->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
766         if(irr) cout <<"error="<<irr<<endl;
767         if(ISetits(4,-1)==0) view1->SetParralel();
768         else  view1->SetPerspective();
769         view1->Front();
770         if(ISetits(3,-1)!=0) view1->ShowAxis();
771     } // end if view1
772     if(ISetits(5,-1)==1) SSD->Raytrace();
773     //arrow->DrawArrow(1.01*Rmax,-Dz,1.01*Rmax,+Dz);
774     //
775     c3->cd(2);
776     SSD->Draw();
777     TPad *p2 = c3->GetPad(2);
778     TView *view2 = p2->GetView();
779     if(view2){
780         view2->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
781         if(irr) cout <<"error="<<irr<<endl;
782         if(ISetits(4,-1)==0) view2->SetParralel();
783         else  view2->SetPerspective();
784         view2->Top();
785         if(ISetits(3,-1)!=0) view2->ShowAxis();
786     } // end if view2
787     if(ISetits(5,-1)==1) SSD->Raytrace();
788     //arrow->DrawArrow(Rmax,0.0,Rmax,0.0);
789     //Double_t s = TMath::Sin(0.7),c = TMath::Cos(0.7);
790     //arrow->DrawArrow(-Rmin*c,-Rmin*s,Rmin*c,Rmin*s);
791     //
792 }
793 //----------------------------------------------------------------------
794 void EngineeringSupRB24(){
795     // Display  RB 24 side cable tray support structure Geometry
796     // Inputs:
797     //    none.
798     // Outputs:
799     //    none.
800     // Retrurn:
801     //    none.
802     Int_t irr;
803     //
804     TGeoManager *mgr2 = gGeoManager;
805     TGeoVolume *ALIC = mgr2->GetTopVolume();
806     TCanvas *c4;
807     if(!(c4 = (TCanvas*)gROOT->FindObject("C4")))
808         c4 = new TCanvas("C4","ITS SDD Cylinder Geometry",900,450);
809     c4->Divide(2,1);
810     TGeoVolume *ITS,*SUPRB24=0;
811     TGeoNode *node;
812     TArrow *arrow=new TArrow();
813     //
814     node = ALIC->FindNode("ITSV_1");
815     ITS = node->GetVolume();
816     node = ITS->FindNode("ITSsupFrameM24_1");
817     SUPRB24 = node->GetVolume();
818     //
819     mgr2->SetNsegments(ISetits(1,-1));
820     //
821     mgr2->SetVisLevel(6);
822     mgr2->SetVisOption(0);
823     //mgr2->CheckOverlaps(0.01);
824     //mgr2->PrintOverlaps();
825     if(ISetits(2,-1)==1){
826         TGeoShape *clip = new TGeoTubeSeg(0, 1000, 2000, 45, 90);
827         mgr2->SetClippingShape(clip);
828     } // end if
829     mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
830     if(ISetits(2,-1)!=0) mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
831     //
832     c4->cd(1);
833     SUPRB24->Draw();
834     TPad *p1 = c4->GetPad(1);
835     TView *view1 = p1->GetView();
836     if(view1){
837         view1->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
838         if(irr) cout <<"error="<<irr<<endl;
839         if(ISetits(4,-1)==0) view1->SetParralel();
840         else  view1->SetPerspective();
841         view1->Front();
842         if(ISetits(3,-1)!=0) view1->ShowAxis();
843     } // end if view1
844     if(ISetits(5,-1)==1) SUPRB24->Raytrace();
845     //
846     c4->cd(2);
847     SUPRB24->Draw();
848     TPad *p2 = c4->GetPad(2);
849     TView *view2 = p2->GetView();
850     if(view2){
851         view2->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
852         if(irr) cout <<"error="<<irr<<endl;
853         if(ISetits(4,-1)==0) view2->SetParralel();
854         else  view2->SetPerspective();
855         view2->Top();
856         if(ISetits(3,-1)!=0) view2->ShowAxis();
857     } // end if view2
858     if(ISetits(5,-1)==1) SUPRB24->Raytrace();
859     //
860 }
861 //----------------------------------------------------------------------
862 void EngineeringSupTrayRB24(){
863     // Display  RB 24 side cable tray support structure Geometry
864     // Inputs:
865     //    none.
866     // Outputs:
867     //    none.
868     // Retrurn:
869     //    none.
870     Int_t irr;
871     //
872     TGeoManager *mgr2 = gGeoManager;
873     TGeoVolume *ALIC = mgr2->GetTopVolume();
874     TCanvas *c4,*c5;
875     //if(!(c4 = (TCanvas*)gROOT->FindObject("C4")))
876         c4 = new TCanvas("C4","ITS RB24 Cable Trays and Patch Pannels",500,500);
877     //c4->Divide(2,1);
878     TGeoVolume *ITS,*SUPRB24=0;
879     TGeoNode *node;
880     TArrow *arrow=new TArrow();
881     //
882     node = ALIC->FindNode("ITSV_1");
883     ITS = node->GetVolume();
884     node = ITS->FindNode("ITSsupCableTrayMotherMT24_1");
885     SUPRB24 = node->GetVolume();
886     //
887     mgr2->SetNsegments(ISetits(1,-1));
888     //
889     mgr2->SetVisLevel(6);
890     mgr2->SetVisOption(0);
891     //mgr2->CheckOverlaps(0.01);
892     //mgr2->PrintOverlaps();
893     if(ISetits(2,-1)==1){
894         TGeoShape *clip = new TGeoTubeSeg(0, 1000, 2000, 45, 90);
895         mgr2->SetClippingShape(clip);
896     } // end if
897     mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
898     if(ISetits(2,-1)!=0) mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
899     //
900     c4->cd(1);
901     SUPRB24->Draw();
902     //TPad *p1 = c4->GetPad(1);
903     //TView *view1 = p1->GetView();
904     TView *view1 = c4->GetView();
905     if(view1){
906         view1->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
907         if(irr) cout <<"error="<<irr<<endl;
908         if(ISetits(4,-1)==0) view1->SetParralel();
909         else  view1->SetPerspective();
910         view1->Front();
911         if(ISetits(3,-1)!=0) view1->ShowAxis();
912     } // end if view1
913     if(ISetits(5,-1)==1) SUPRB24->Raytrace();
914     //
915     //c4->cd(2);
916         c5 = new TCanvas("C5","ITS RB24 Cable Trays and Patch Pannels",500,500);
917     c5->cd(1);
918     SUPRB24->Draw();
919     //TPad *p2 = c5->GetPad(1);
920     //TView *view2 = p2->GetView();
921     TView *view2 = c5->GetView();
922     if(view2){
923         view2->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
924         if(irr) cout <<"error="<<irr<<endl;
925         if(ISetits(4,-1)==0) view2->SetParralel();
926         else  view2->SetPerspective();
927         view2->Top();
928         if(ISetits(3,-1)!=0) view2->ShowAxis();
929     } // end if view2
930     if(ISetits(5,-1)==1) SUPRB24->Raytrace();
931     //
932 }
933 //----------------------------------------------------------------------
934 void EngineeringSupRB26(){
935     // Display RB 26 side cable tray support structure
936     // Inputs:
937     //    none.
938     // Outputs:
939     //    none.
940     // Retrurn:
941     //    none.
942     Int_t irr;
943     //
944     TGeoManager *mgr2 = gGeoManager;
945     TGeoVolume *ALIC = mgr2->GetTopVolume();
946     TCanvas *c5;
947     if(!(c5 = (TCanvas*)gROOT->FindObject("C5")))
948         c5 = new TCanvas("C5","ITS SDD Cylinder Geometry",900,450);
949     c5->Divide(2,1);
950     TGeoVolume *ITS,*SUPRB26=0;
951     TGeoNode *node;
952     TArrow *arrow=new TArrow();
953     //
954     node = ALIC->FindNode("ITSV_1");
955     ITS = node->GetVolume();
956     node = ITS->FindNode("ITSsupFrameM26_1");
957     SUPRB26 = node->GetVolume();
958     //
959     mgr2->SetNsegments(ISetits(1,-1));
960     //
961     mgr2->SetVisLevel(6);
962     mgr2->SetVisOption(0);
963     //mgr2->CheckOverlaps(0.01);
964     //mgr2->PrintOverlaps();
965     if(ISetits(2,-1)==1){
966         TGeoShape *clip = new TGeoTubeSeg(0, 1000, 2000, 45, 90);
967         mgr2->SetClippingShape(clip);
968     } // end if
969     mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
970     if(ISetits(2,-1)!=0) mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
971     //
972     c5->cd(1);
973     SUPRB26->Draw();
974     TPad *p1 = c5->GetPad(1);
975     TView *view1 = p1->GetView();
976     if(view1){
977         view1->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
978         if(irr) cout <<"error="<<irr<<endl;
979         if(ISetits(4,-1)==0) view1->SetParralel();
980         else  view1->SetPerspective();
981         view1->Front();
982         if(ISetits(3,-1)!=0) view1->ShowAxis();
983     } // end if view1
984     if(ISetits(5,-1)==1) SUPRB26->Raytrace();
985     //
986     c5->cd(2);
987     SUPRB26->Draw();
988     TPad *p2 = c5->GetPad(2);
989     TView *view2 = p2->GetView();
990     if(view2){
991         view2->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
992         if(irr) cout <<"error="<<irr<<endl;
993         if(ISetits(4,-1)==0) view2->SetParralel();
994         else  view2->SetPerspective();
995         view2->Top();
996         if(ISetits(3,-1)!=0) view2->ShowAxis();
997     } // end if view2
998     if(ISetits(5,-1)==1) SUPRB26->Raytrace();
999     //
1000 }