RC12, RC17 violation: suppression
[u/mrichter/AliRoot.git] / ITS / Displayv11.C
1 //----------------------------------------------------------------------
2
3 static AliITSv11GeometrySPD *gspd;
4 static AliITSv11GeometrySDD *gsdd;
5 static AliITSv11GeometrySSD *gssd;
6 static AliITSv11GeometrySupport *gsupp;
7 static AliITSgeom *geom;
8 //
9 //----------------------------------------------------------------------
10 void Displayv11(const char* filename=""){
11     // Display AliITSv11 Geometry
12     // Inputs:
13     //    const char* filename output file with the display in it
14     // Outputs:
15     //    none.
16     // Retrurn:
17     //    none.
18
19     gSystem->Load("libGeom");
20     //
21     if(gGeoManager) delete gGeoManager;
22     mgr2 = gGeoManager = new TGeoManager("ITSGeometry",
23                                          " ITS Simulation Geometry Manager");
24     //
25     const AliITSVersion_t kv11=110;
26     const Char_t *cvsDate="$Date$";
27     const Char_t *cvsRevision="$Revision$";
28     const Int_t kLength=100;
29     Char_t vstrng[kLength];
30     AliITSInitGeometry initgeom(kv11,1);
31     //
32     TGeoMaterial *vacmat = new TGeoMaterial("Vacume",0,0,0);
33     TGeoMedium   *vacmed = new TGeoMedium("Vacume_med",1,vacmat);
34     TGeoVolume *ALIC = mgr2->MakeBox("ALIC",vacmed,1000.,1000.,2000.);
35     mgr2->SetTopVolume(ALIC);
36     TGeoVolume *ITS = mgr2->MakeBox("ITSV",vacmed,990.,990.,1990.);
37     if(initgeom.WriteVersionString(vstrng,kLength,kv11,1,cvsDate,cvsRevision))
38         ITS->SetTitle(vstrng);
39     //TGeoVolumeAssembly *ITSSPD = new TGeoVolumeAssembly("ITSSPD");
40     //ITS->AddNode(ITSSPD,1);
41     ALIC->AddNode(ITS,1);
42     //
43     /*
44     AliITSv11 *its = new AliITSv11(0,3);
45     its->SetDebug(ISetits(0,-1));
46     its->GetSPDGeometry()->SetDebug(ISetits(0,-1));
47     its->GetSupGeometry()->SetDebug(ISetits(0,-1));
48     its->CreateMaterials();
49     its->CreateGeometry();
50     */
51     gspd  = new AliITSv11GeometrySPD(0);
52     gsdd  = new AliITSv11GeometrySDD();
53     gsupp = new AliITSv11GeometrySupport(0);
54     gssd  = new AliITSv11GeometrySSD();
55     //
56     Int_t imat=1,imed=1,ireturn=0;
57     ireturn = gspd->CreateSPDCentralMaterials(imed,imat);
58     gspd->SPDSector(ITS,mgr2);
59     gsupp->SPDCone(ITS,mgr2);
60     gsupp->SetDebug(0);
61     gsupp->SDDCone(ITS,mgr2);
62     gsdd->Layer3(ITS);
63     gsdd->Layer4(ITS);
64     gsdd->ForwardLayer3(ITS);// in Hybrid its in IS02
65     gsdd->ForwardLayer4(ITS);// in Hybrid its in IS02
66     gssd->Layer5(ITS);
67     gssd->Layer6(ITS);
68     gssd->LadderSupportLayer5(ITS);
69     gssd->LadderSupportLayer6(ITS);
70     gssd->EndCapSupportSystemLayer6(ITS);
71     gssd->EndCapSupportSystemLayer5(ITS);
72     gsupp->SSDCone(ITS,mgr2);
73     gsupp->ServicesCableSupport(ITS);
74     //
75     mgr2->CloseGeometry();
76     //
77     geom = new AliITSgeom();
78     initgeom.InitAliITSgeom(geom);
79     geom->WriteNewFile("ITSgeomV11.det");
80     //
81     TControlBar *bar=new TControlBar("vertical","ITS Geometry Display",10,10);
82     bar->AddButton("Set Clipping on","ISetits(2,1)","Clipping on");
83     bar->AddButton("Set Cllipping off","ISetits(2,0)","Clipping off");
84     bar->AddButton("Set axis on","ISetits(3,1)","Show Axis on");
85     bar->AddButton("Set axis off","ISetits(3,0)","Show Axis off");
86     bar->AddButton("Set perspective on","ISetits(4,1)","Perspective on");
87     bar->AddButton("Set perspective off","ISetits(4,0)","Perspective off");
88     bar->AddButton("Set RayTrace on","ISetits(5,1)","Perspective on");
89     bar->AddButton("Set RayTrace off","ISetits(5,0)","Perspective off");
90     bar->AddButton("Set circle/80","ISetits(1,80)","circles ~ by 80 lines");
91     bar->AddButton("Display Geometry","Displayit()","Run Displayit");
92     bar->AddButton("Display SPD Sector Volume","EngineeringSPDSector()",
93                    "Run EngineeringSPDSector");
94     bar->AddButton("Print SPD Sector Volume data xfig","PrintSPDSectorData()",
95                    "Run PrintSPDSectorData");
96     bar->AddButton("Display SPD Ceneral Volume","EngineeringSPDCenter()",
97                    "Run EngineeringSPDCenter");
98     bar->AddButton("Display SPD Thermal Sheald","EngineeringSPDThS()",
99                    "Run EngineeringSPDThS");
100     bar->AddButton("Display SPD Sector Cross Sections","EngineeringSPDSCS()",
101                    "Run EngineeringSPDSCS");
102     bar->AddButton("Display SDD Layer 3","EngineeringSDDLayer3()",
103                    "Run EngineeringSDDLayer3");
104     bar->AddButton("Display SDD Layer 4","EngineeringSDDLayer4()",
105                    "Run EngineeringSDDLayer4");
106     bar->AddButton("Display SDD Cone","EngineeringSDDCone()",
107                    "Run EngineeringSDDCone");
108     bar->AddButton("Display SDD Centeral Cylinder","EngineeringSDDCylinder()",
109                    "Run EngineeringSDDCylinder");
110     bar->AddButton("Display SSD Cone","EngineeringSSDCone()",
111                    "Run EngineeringSSDCone");
112     bar->AddButton("Display SSD Centeral Cylinder","EngineeringSSDCylinder()",
113                    "Run EngineeringSSDCylinder");
114     bar->AddButton("Display SUP RB24 side","EngineeringSupRB24()",
115                    "Run EngineeringSupRB24");
116     bar->AddButton("Display Cable Trays RB24 side","EngineeringSupTrayRB24()",
117                    "Run EngineeringSupTrayRB24");
118     bar->AddButton("Display SUP RB26 side","EngineeringSupRB26()",
119                    "Run EngineeringSupRB26");
120     bar->AddButton("Save Geometry to File","ExportToFile()",
121                    "Run ExportToFile");
122     bar->AddButton("Quit/Exit",".q","Exit");
123     bar->Show();
124     gROOT->SaveContext();
125          //Displayit();
126 }
127 //----------------------------------------------------------------------
128 void ExportToFile(){
129     // Quirry file name and write geometry to a root file.
130     // Inputs:
131     //    const char* filename output file with the display in it
132     // Outputs:
133     //    none.
134     // Retrurn:
135     //    none.
136     Char_t filename[100];
137
138     printf("Eneter File name:");
139     scanf("%s",filename);
140
141     gGeoManager->Export(filename);
142 }
143 //----------------------------------------------------------------------
144 Int_t ISetits(Int_t t,Int_t v){
145     static Int_t itsdebug=1,nsegments=80,cut=0,axis=1,perspective=0,ray=0;
146
147     switch (t) {
148     case 0:
149         if(v<0) return itsdebug;
150         itsdebug = v;
151         break;
152     case 1:
153         if(v<0) return nsegments;
154         nsegments= v;
155         break;
156     case 2:
157         if(v<0) return cut;
158         cut = v;
159         break;
160     case 3:
161         if(v<0) return axis;
162         axis = v;
163         break;
164     case 4:
165         if(v<0) return perspective;
166         perspective = v;
167         break;
168     case 5:
169         if(v<0) return ray;
170         ray = v;
171         break;
172     }// end switch
173     return 0;
174 }
175 //----------------------------------------------------------------------
176 Double_t DSetits(Int_t t,Double_t v){
177     static Double_t phimincut=0.0,phimaxcut=180.0;
178     static Double_t longitude=90.0,latitude=0.0,psi=0.0;
179
180     switch (t) {
181     case 0:
182         if(v<0.) return phimincut;
183         phimincut = v;
184         break;
185     case 1:
186         if(v<0.) return phimaxcut;
187         phimaxcut = v;
188         break;
189     case 2:
190         if(v<0.) return longitude;
191         longitude = v;
192         break;
193     case 3:
194         if(v<0.) return latitude;
195         latitude = v;
196         break;
197     case 4:
198         if(v<0.) return latitude;
199         latitude = v;
200         break;
201     }// end switch
202     return 0;
203 }
204 //----------------------------------------------------------------------
205 void Displayit(){
206     // Display AliITSv11 Geometry
207     // Inputs:
208     //    const char* filename output file with the display in it
209     // Outputs:
210     //    none.
211     // Retrurn:
212     //    none.
213     Int_t irr;
214     //
215     TGeoManager *mgr2 = gGeoManager;
216     TGeoVolume *ALIC = mgr2->GetTopVolume();
217     TCanvas *c1;
218     if(!(c1 = (TCanvas*)gROOT->FindObject("C1")))
219         c1 = new TCanvas("C1","ITS Simulation Geometry",900,900);
220     //c1->Divide(2,2);
221     //
222     mgr2->SetNsegments(ISetits(1,-1));
223     //
224     mgr2->SetVisLevel(6);
225     mgr2->SetVisOption(0);
226     //mgr2->CheckOverlaps(0.01);
227     //mgr2->PrintOverlaps();
228     if(ISetits(2,-1)==1){
229         TGeoShape *clip = new TGeoTubeSeg(0, 1000, 2000, 45, 90);
230         mgr2->SetClippingShape(clip);
231     } // end if
232     if(ISetits(2,-1)!=0) mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
233     //
234     c1->cd(1);
235     ALIC->Draw();
236     TPad *p1 = c1->GetPad(1);
237     TView *view1 = p1->GetView();
238     if(view1){
239         view1->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
240         if(irr) cout <<"error="<<irr<<endl;
241         if(ISetits(4,-1)==0) view1->SetParallel();
242         else  view1->SetPerspective();
243         view1->Front();
244         if(ISetits(3,-1)!=0) view1->ShowAxis();
245     } // end if view1
246     if(ISetits(5,-1)==1) ALIC->Raytrace();
247     /*c1->cd(2);
248     ALIC->Draw();
249     TPad *p2 = c1->GetPad(2);
250     TView *view2 = p2->GetView();
251     if(view2){
252         view2->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
253         if(irr) cout <<"error="<<irr<<endl;
254         if(ISetits(4,-1)==0) view2->SetParallel();
255         else  view2->SetPerspective();
256         view2->RotateView(60.,30.);
257         if(ISetits(3,-1)!=0) view2->ShowAxis();
258     } // end if view2
259     if(ISetits(5,-1)==1) ALIC->Raytrace();
260     c1->cd(3);
261     ALIC->Draw();
262     c1->SetPhi(90.0); c1->SetTheta(90.0);
263     TPad *p3 = c1->GetPad(3);
264     TView *view3 = p3->GetView();
265     if(view3){
266         view3->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
267         if(irr) cout <<"error="<<irr<<endl;
268         if(ISetits(4,-1)==0) view3->SetParallel();
269         else  view3->SetPerspective();
270         view3->Top();
271         if(ISetits(3,-1)!=0) view3->ShowAxis();
272     } // end if view3
273     if(ISetits(5,-1)==1) ALIC->Raytrace();
274     c1->cd(4);
275     ALIC->Draw();
276     TPad *p4 = c1->GetPad(4);
277     TView *view4 = p4->GetView();
278     if(view4){
279         view4->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
280         if(irr) cout <<"error="<<irr<<endl;
281         if(ISetits(4,-1)==0) view4->SetParallel();
282         else  view4->SetPerspective();
283         view4->Side();
284         if(ISetits(3,-1)!=0) view4->ShowAxis();
285     } // end if view4
286     if(ISetits(5,-1)==1) ALIC->Raytrace();
287     *///
288 }
289 //----------------------------------------------------------------------
290 void EngineeringSPDSCS(){
291     // Display SPD Carbon Fiber Sector Cross sections A and B
292     // Inputs:
293     //    none.
294     // Outputs:
295     //    none.
296     // Return:
297     //    none.
298     TPolyLine *a0,*a1,*b0,*b1;
299     TPolyMarker *p;
300     TCanvas *cSPDSCS=0;
301     Int_t i;
302     Double_t max=0.0;
303
304     a0 = new TPolyLine();
305     a1 = new TPolyLine();
306     b0 = new TPolyLine();
307     b1 = new TPolyLine();
308     p = new TPolyMarker();
309     a0->SetLineColor(1);
310     a1->SetLineColor(4);
311     b0->SetLineColor(3);
312     b0->SetLineStyle(2); // dashed
313     b1->SetLineColor(6);
314     b1->SetLineStyle(2); // dashed
315     p->SetMarkerColor(2);
316     p->SetMarkerStyle(5);
317     if(gspd==0||gspd->Make2DcrossSections(*a0,*a1,*b0,*b1,*p)==kFALSE)return;
318     for(i=0;i<a0->GetN();i++) {
319       if(TMath::Abs(a0->GetX()[i])>max) max = TMath::Abs(a0->GetX()[i]);
320       if(TMath::Abs(a0->GetY()[i])>max) max = TMath::Abs(a0->GetY()[i]);
321     } // end for i
322     for(i=0;i<a1->GetN();i++) {
323       if(TMath::Abs(a1->GetX()[i])>max) max = TMath::Abs(a0->GetX()[i]);
324       if(TMath::Abs(a1->GetY()[i])>max) max = TMath::Abs(a1->GetY()[i]);
325     } // end for i
326     for(i=0;i<b0->GetN();i++) {
327       if(TMath::Abs(b0->GetX()[i])>max) max = TMath::Abs(b0->GetX()[i]);
328       if(TMath::Abs(b0->GetY()[i])>max) max = TMath::Abs(b0->GetY()[i]);
329     } // end for i
330     for(i=0;i<b1->GetN();i++) {
331       if(TMath::Abs(b1->GetX()[i])>max) max = TMath::Abs(b1->GetX()[i]);
332       if(TMath::Abs(b1->GetY()[i])>max) max = TMath::Abs(b1->GetY()[i]);
333     } // end for i
334     max *= 1.05;
335     cSPDSCS = gROOT->FindObject("cSPDSCS");
336     if(cSPDSCS==0) delete cSPDSCS;
337     cSPDSCS = new TCanvas("cSPDSCS","SPD Carbon Fiber Sector Cross sections",2);
338     cSPDSCS->Range(-max,-max,max,max);
339     //cSPDSCS->RangeAxis();
340     cSPDSCS->SetFixedAspectRatio(kTRUE);
341     //
342     a0->Draw("");
343     a1->Draw("same");
344     p->Draw("same");
345     b0->Draw("same");
346     b1->Draw("same");
347     cSPDSCS->Update();
348     return;
349 }
350 //----------------------------------------------------------------------
351 void EngineeringSPDSector(){
352     // Display SPD Sector Geometry
353     // Inputs:
354     //    none.
355     // Outputs:
356     //    none.
357     // Retrurn:
358     //    none.
359     Int_t irr;
360     //
361     TGeoManager *mgr2 = gGeoManager;
362     TGeoVolume *ALIC = mgr2->GetTopVolume();
363     TCanvas *c4=0;
364     if(!(c4 = (TCanvas*)gROOT->FindObject("C4")))
365         c4 = new TCanvas("C4","ITS SPD Layer Geometry Side View",500,500);
366     TGeoVolume *ITS=0,*ITSSPD=0,*SPDLay=0;
367     TGeoNode *node=0;
368     TArrow *arrow=new TArrow();
369     //
370     node = ALIC->FindNode("ITSV_1");
371     ITS = node->GetVolume();
372     node = ITS->FindNode("ITSSPD_1");
373     ITSSPD = node->GetVolume();
374     node = ITSSPD->FindNode("ITSSPDCarbonFiberSectorV_1");
375     if(node==0)Error("EngineeringSPDSector","could not find node %s",
376                      "ITSSPDCarbonFiberSectorV_1");
377     SPDLay = node->GetVolume();
378     if(SPDLay==0)Error("EngineeringSPDSector","could not find volume SPDLay");
379     //
380     mgr2->SetNsegments(ISetits(1,-1));
381     //
382     mgr2->SetVisLevel(6);
383     mgr2->SetVisOption(0);
384     //mgr2->CheckOverlaps(0.01);
385     //mgr2->PrintOverlaps();
386     if(ISetits(2,-1)==1){
387         TGeoShape *clip = new TGeoTubeSeg(0, 1000, 2000, 45, 90);
388         mgr2->SetClippingShape(clip);
389     } // end if
390     mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
391     if(ISetits(2,-1)!=0) mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
392     //
393     SPDLay->Draw();
394     TView *view1 = c4->GetView();
395     if(view1){
396         view1->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
397         if(irr) cout <<"error="<<irr<<endl;
398         if(ISetits(4,-1)==0) view1->SetParallel();
399         else  view1->SetPerspective();
400         view1->Front();
401         if(ISetits(3,-1)!=0) view1->ShowAxis();
402     } // end if view1
403     if(ISetits(5,-1)==1) SPDLay->Raytrace();
404     //
405     if(!(c5 = (TCanvas*)gROOT->FindObject("C5")))
406         c5 = new TCanvas("C5","ITS SPD Sector Geometry End View",500,500);
407     SPDLay->Draw();
408     TView *view2 = c5->GetView();
409     if(view2){
410         view2->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
411         if(irr) cout <<"error="<<irr<<endl;
412         if(ISetits(4,-1)==0) view2->SetParallel();
413         else  view2->SetPerspective();
414         view2->Top();
415         if(ISetits(3,-1)!=0) view2->ShowAxis();
416     } // end if view2
417     if(ISetits(5,-1)==1) SPDLay->Raytrace();
418     //
419 }
420 //----------------------------------------------------------------------
421 void PrintSPDSectorData(){
422     // Print SPD Sector Data
423     // Inputs:
424     //    none.
425     // Outputs:
426     //    none.
427     // Retrurn:
428     //    none.
429     Int_t irr,i;
430     //
431     TGeoManager *mgr2 = gGeoManager;
432     TGeoXtru * sA0;
433     TGeoVolume *ITS,*ITSSPD,*vA0=0;
434     TGeoNode *node;
435
436     //mgr2->PushPath();
437     //mgr2->cd("ITSSPDCarbonFiberSupportSectorA0_1");
438     vA0 = mgr2->FindVolumeFast("ITSSPDCarbonFiberSupportSectorA0");
439     sA0 = (TGeoXtru*) vA0->GetShape();
440     irr = sA0->GetNvert();
441     Double_t x,y;
442     cout <<endl;
443     cout <<"2 3 0 1 0 7 50 -1 -1 0.000 0 0 -1 0 0 "<<irr;
444     for(i=0;i<irr;i++){
445         x = sA0->GetX(i)+2.5;
446         y = sA0->GetY(i)+2.5;
447         if(!(i%6)) { cout << endl; cout <<"        ";}
448         cout<<" "<<TMath::Nint(x*450.)<<" "<<TMath::Nint(y*450);
449         //cout<<" "<<x<<" "<<y;
450     } // end for i
451     x = sA0->GetX(0)+2.5;
452     y = sA0->GetY(0)+2.5;
453     if(!(i%6)) { cout << endl; cout <<"        ";}
454     cout<<" "<<TMath::Nint(x*450.)<<" "<<TMath::Nint(y*450);
455     //cout<<" "<<x<<" "<<y;
456     cout << endl;
457     //
458 }
459 //----------------------------------------------------------------------
460 void EngineeringSPDCenter(){
461     // Display SPD Centeral Geometry
462     // Inputs:
463     //    none.
464     // Outputs:
465     //    none.
466     // Retrurn:
467     //    none.
468     Int_t irr;
469     //
470     TGeoManager *mgr2 = gGeoManager;
471     TGeoVolume *ALIC = mgr2->GetTopVolume();
472     TCanvas *c4;
473     if(!(c4 = (TCanvas*)gROOT->FindObject("C4")))
474         c4 = new TCanvas("C4","ITS SPD Layer Geometry Side View",500,500);
475     TGeoVolume *ITS,*ITSSPD,*SPDLay=0;
476     TGeoNode *node;
477     TArrow *arrow=new TArrow();
478     //
479     node = ALIC->FindNode("ITSV_1");
480     ITS = node->GetVolume();
481     node = ITS->FindNode("ITSSPD_1");
482     ITSSPD = node->GetVolume();
483     //
484     mgr2->SetNsegments(ISetits(1,-1));
485     //
486     mgr2->SetVisLevel(6);
487     mgr2->SetVisOption(0);
488     //mgr2->CheckOverlaps(0.01);
489     //mgr2->PrintOverlaps();
490     if(ISetits(2,-1)==1){
491         TGeoShape *clip = new TGeoTubeSeg(0, 1000, 2000, 45, 90);
492         mgr2->SetClippingShape(clip);
493     } // end if
494     mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
495     if(ISetits(2,-1)!=0) mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
496     //
497     ITSSPD->Draw();
498     TView *view1 = c4->GetView();
499     if(view1){
500         view1->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
501         if(irr) cout <<"error="<<irr<<endl;
502         if(ISetits(4,-1)==0) view1->SetParallel();
503         else  view1->SetPerspective();
504         view1->Front();
505         if(ISetits(3,-1)!=0) view1->ShowAxis();
506     } // end if view1
507     if(ISetits(5,-1)==1) SPDLay->Raytrace();
508     //
509     if(!(c5 = (TCanvas*)gROOT->FindObject("C5")))
510         c5 = new TCanvas("C5","ITS SPD Centeral Geometry End View",500,500);
511     ITSSPD->Draw();
512     TView *view2 = c5->GetView();
513     if(view2){
514         view2->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
515         if(irr) cout <<"error="<<irr<<endl;
516         if(ISetits(4,-1)==0) view2->SetParallel();
517         else  view2->SetPerspective();
518         view2->Top();
519         if(ISetits(3,-1)!=0) view2->ShowAxis();
520     } // end if view2
521     if(ISetits(5,-1)==1) ITSSPD->Raytrace();
522     //
523 }
524 //----------------------------------------------------------------------
525 void EngineeringSPDThS(){
526     // Display SPD Thermal Sheald Geometry
527     // Inputs:
528     //    none.
529     // Outputs:
530     //    none.
531     // Retrurn:
532     //    none.
533     Int_t irr;
534     //
535     TGeoManager *mgr2 = gGeoManager;
536     TGeoVolume *ALIC = mgr2->GetTopVolume();
537     TCanvas *c4;
538     if(!(c4 = (TCanvas*)gROOT->FindObject("C4")))
539         c4 = new TCanvas("C4","ITS SDD Cylinder Geometry",900,450);
540     c4->Divide(2,1);
541     TGeoVolume *ITS,*SPDThS=0;
542     TGeoNode *node;
543     TArrow *arrow=new TArrow();
544     //
545     node = ALIC->FindNode("ITSV_1");
546     ITS = node->GetVolume();
547     node = ITS->FindNode("ITSspdThermalShealdMother_1");
548     SPDThS = node->GetVolume();
549     //
550     mgr2->SetNsegments(ISetits(1,-1));
551     //
552     mgr2->SetVisLevel(6);
553     mgr2->SetVisOption(0);
554     //mgr2->CheckOverlaps(0.01);
555     //mgr2->PrintOverlaps();
556     if(ISetits(2,-1)==1){
557         TGeoShape *clip = new TGeoTubeSeg(0, 1000, 2000, 45, 90);
558         mgr2->SetClippingShape(clip);
559     } // end if
560     mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
561     if(ISetits(2,-1)!=0) mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
562     //
563     c4->cd(1);
564     SPDThS->Draw();
565     TPad *p1 = c4->GetPad(1);
566     TView *view1 = p1->GetView();
567     if(view1){
568         view1->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
569         if(irr) cout <<"error="<<irr<<endl;
570         if(ISetits(4,-1)==0) view1->SetParallel();
571         else  view1->SetPerspective();
572         view1->Front();
573         if(ISetits(3,-1)!=0) view1->ShowAxis();
574     } // end if view1
575     if(ISetits(5,-1)==1) SPDThS->Raytrace();
576     //
577     c4->cd(2);
578     SPDThS->Draw();
579     TPad *p2 = c4->GetPad(2);
580     TView *view2 = p2->GetView();
581     if(view2){
582         view2->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
583         if(irr) cout <<"error="<<irr<<endl;
584         if(ISetits(4,-1)==0) view2->SetParallel();
585         else  view2->SetPerspective();
586         view2->Top();
587         if(ISetits(3,-1)!=0) view2->ShowAxis();
588     } // end if view2
589     if(ISetits(5,-1)==1) SPDThS->Raytrace();
590     //
591 }
592 //----------------------------------------------------------------------
593 void EngineeringSDDLayer3(){
594     // Display SDD Layer 3 Geometry
595     // Inputs:
596     //    none.
597     // Outputs:
598     //    none.
599     // Retrurn:
600     //    none.
601     Int_t irr;
602     //
603     TGeoManager *mgr2 = gGeoManager;
604     TGeoVolume *ALIC = mgr2->GetTopVolume();
605     TCanvas *c4;
606     if(!(c4 = (TCanvas*)gROOT->FindObject("C4")))
607         c4 = new TCanvas("C4","ITS SDD Layer 3 Geometry Side View",500,500);
608     TGeoVolume *ITS,*SDDLay3=0;
609     TGeoNode *node;
610     TArrow *arrow=new TArrow();
611     //
612     node = ALIC->FindNode("ITSV_1");
613     ITS = node->GetVolume();
614     node = ITS->FindNode("ITSsddLayer3_1");
615     SDDLay3 = node->GetVolume();
616     //
617     mgr2->SetNsegments(ISetits(1,-1));
618     //
619     mgr2->SetVisLevel(6);
620     mgr2->SetVisOption(0);
621     //mgr2->CheckOverlaps(0.01);
622     //mgr2->PrintOverlaps();
623     if(ISetits(2,-1)==1){
624         TGeoShape *clip = new TGeoTubeSeg(0, 1000, 2000, 45, 90);
625         mgr2->SetClippingShape(clip);
626     } // end if
627     mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
628     if(ISetits(2,-1)!=0) mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
629     //
630     SDDLay3->Draw();
631     TView *view1 = c4->GetView();
632     if(view1){
633         view1->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
634         if(irr) cout <<"error="<<irr<<endl;
635         if(ISetits(4,-1)==0) view1->SetParallel();
636         else  view1->SetPerspective();
637         view1->Front();
638         if(ISetits(3,-1)!=0) view1->ShowAxis();
639     } // end if view1
640     if(ISetits(5,-1)==1) SDDLay3->Raytrace();
641     //
642     if(!(c5 = (TCanvas*)gROOT->FindObject("C5")))
643         c5 = new TCanvas("C5","ITS SDD Layer 3 Geometry End View",500,500);
644     SDDLay3->Draw();
645     TView *view2 = c5->GetView();
646     if(view2){
647         view2->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
648         if(irr) cout <<"error="<<irr<<endl;
649         if(ISetits(4,-1)==0) view2->SetParallel();
650         else  view2->SetPerspective();
651         view2->Top();
652         if(ISetits(3,-1)!=0) view2->ShowAxis();
653     } // end if view2
654     if(ISetits(5,-1)==1) SDDLay3->Raytrace();
655     //
656 }
657 //----------------------------------------------------------------------
658 void EngineeringSDDLayer4(){
659     // Display SDD Layer 4 Geometry
660     // Inputs:
661     //    none.
662     // Outputs:
663     //    none.
664     // Retrurn:
665     //    none.
666     Int_t irr;
667     //
668     TGeoManager *mgr2 = gGeoManager;
669     TGeoVolume *ALIC = mgr2->GetTopVolume();
670     TCanvas *c4;
671     if(!(c4 = (TCanvas*)gROOT->FindObject("C4")))
672         c4 = new TCanvas("C4","ITS SDD Layer 4 Geometry Side View",500,500);
673     TGeoVolume *ITS,*SDDLay4=0;
674     TGeoNode *node;
675     TArrow *arrow=new TArrow();
676     //
677     node = ALIC->FindNode("ITSV_1");
678     ITS = node->GetVolume();
679     node = ITS->FindNode("ITSsddLayer4_1");
680     SDDLay4 = node->GetVolume();
681     //
682     mgr2->SetNsegments(ISetits(1,-1));
683     //
684     mgr2->SetVisLevel(6);
685     mgr2->SetVisOption(0);
686     //mgr2->CheckOverlaps(0.01);
687     //mgr2->PrintOverlaps();
688     if(ISetits(2,-1)==1){
689         TGeoShape *clip = new TGeoTubeSeg(0, 1000, 2000, 45, 90);
690         mgr2->SetClippingShape(clip);
691     } // end if
692     mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
693     if(ISetits(2,-1)!=0) mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
694     //
695     SDDLay4->Draw();
696     TView *view1 = c4->GetView();
697     if(view1){
698         view1->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
699         if(irr) cout <<"error="<<irr<<endl;
700         if(ISetits(4,-1)==0) view1->SetParallel();
701         else  view1->SetPerspective();
702         view1->Front();
703         if(ISetits(3,-1)!=0) view1->ShowAxis();
704     } // end if view1
705     if(ISetits(5,-1)==1) SDDLay4->Raytrace();
706     //
707     if(!(c5 = (TCanvas*)gROOT->FindObject("C5")))
708         c5 = new TCanvas("C5","ITS SDD Layer 4 Geometry End View",500,500);
709     SDDLay4->Draw();
710     TView *view2 = c5->GetView();
711     if(view2){
712         view2->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
713         if(irr) cout <<"error="<<irr<<endl;
714         if(ISetits(4,-1)==0) view2->SetParallel();
715         else  view2->SetPerspective();
716         view2->Top();
717         if(ISetits(3,-1)!=0) view2->ShowAxis();
718     } // end if view2
719     if(ISetits(5,-1)==1) SDDLay4->Raytrace();
720     //
721 }
722 //----------------------------------------------------------------------
723 void EngineeringSDDCone(){
724     // Display SDD Cone Geometry
725     // Inputs:
726     //    none.
727     // Outputs:
728     //    none.
729     // Retrurn:
730     //    none.
731     Int_t irr;
732     //
733     TGeoManager *mgr2 = gGeoManager;
734     TGeoVolume *ALIC = mgr2->GetTopVolume();
735     TCanvas *c2;
736     if(!(c2 = (TCanvas*)gROOT->FindObject("C2")))
737         c2 = new TCanvas("C2","ITS SDD Cone Geometry",900,450);
738     c2->Divide(2,1);
739     TGeoVolume *ITS,*SDD=0;
740     TGeoNode *node;
741     //
742     node = ALIC->FindNode("ITSV_1");
743     ITS = node->GetVolume();
744     node = ITS->FindNode("ITSsddConeL_1");
745     SDD = node->GetVolume();
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     c2->cd(1);
761     SDD->Draw();
762     TPad *p1 = c2->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->SetParallel();
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) SDD->Raytrace();
773     //
774     c2->cd(2);
775     SDD->Draw();
776     TPad *p2 = c2->GetPad(2);
777     TView *view2 = p2->GetView();
778     if(view2){
779         view2->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
780         if(irr) cout <<"error="<<irr<<endl;
781         if(ISetits(4,-1)==0) view2->SetParallel();
782         else  view2->SetPerspective();
783         view2->Top();
784         if(ISetits(3,-1)!=0) view2->ShowAxis();
785     } // end if view1
786     if(ISetits(5,-1)==1) SDD->Raytrace();
787     //
788 }
789 //----------------------------------------------------------------------
790 void EngineeringSDDCylinder(){
791     // Display SDD Cylinder Geometry
792     // Inputs:
793     //    none.
794     // Outputs:
795     //    none.
796     // Retrurn:
797     //    none.
798     Int_t irr;
799     //
800     TGeoManager *mgr2 = gGeoManager;
801     TGeoVolume *ALIC = mgr2->GetTopVolume();
802     TCanvas *c3;
803     if(!(c3 = (TCanvas*)gROOT->FindObject("C3")))
804         c3 = new TCanvas("C3","ITS SDD Cylinder Geometry",900,450);
805     c3->Divide(2,1);
806     TGeoVolume *ITS,*SDD=0;
807     TGeoNode *node;
808     TArrow *arrow=new TArrow();
809     //
810     node = ALIC->FindNode("ITSV_1");
811     ITS = node->GetVolume();
812     node = ITS->FindNode("ITSsddCentCylCF_1");
813     SDD = node->GetVolume();
814     Double_t Rmin = ((TGeoTube*)(SDD->GetShape()))->GetRmin();
815     Double_t Rmax = ((TGeoTube*)(SDD->GetShape()))->GetRmax();
816     Double_t Dz   = ((TGeoTube*)(SDD->GetShape()))->GetDz();
817     //
818     mgr2->SetNsegments(ISetits(1,-1));
819     //
820     mgr2->SetVisLevel(6);
821     mgr2->SetVisOption(0);
822     //mgr2->CheckOverlaps(0.01);
823     //mgr2->PrintOverlaps();
824     if(ISetits(2,-1)==1){
825         TGeoShape *clip = new TGeoTubeSeg(0, 1000, 2000, 45, 90);
826         mgr2->SetClippingShape(clip);
827     } // end if
828     mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
829     if(ISetits(2,-1)!=0) mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
830     //
831     c3->cd(1);
832     SDD->Draw();
833     TPad *p1 = c3->GetPad(1);
834     TView *view1 = p1->GetView();
835     if(view1){
836         view1->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
837         if(irr) cout <<"error="<<irr<<endl;
838         if(ISetits(4,-1)==0) view1->SetParallel();
839         else  view1->SetPerspective();
840         view1->Front();
841         if(ISetits(3,-1)!=0) view1->ShowAxis();
842     } // end if view1
843     if(ISetits(5,-1)==1) SDD->Raytrace();
844     //arrow->DrawArrow(1.01*Rmax,-Dz,1.01*Rmax,+Dz);
845     //
846     c3->cd(2);
847     SDD->Draw();
848     TPad *p2 = c3->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->SetParallel();
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) SDD->Raytrace();
859     //arrow->DrawArrow(Rmax,0.0,Rmax,0.0);
860     //Double_t s = TMath::Sin(0.7),c = TMath::Cos(0.7);
861     //arrow->DrawArrow(-Rmin*c,-Rmin*s,Rmin*c,Rmin*s);
862     //
863 }
864 //----------------------------------------------------------------------
865 void EngineeringSSDCone(){
866     // Display SSD Cone Geometry
867     // Inputs:
868     //    none.
869     // Outputs:
870     //    none.
871     // Retrurn:
872     //    none.
873     Int_t irr;
874     //
875     TGeoManager *mgr2 = gGeoManager;
876     TGeoVolume *ALIC = mgr2->GetTopVolume();
877     TCanvas *c2;
878     if(!(c2 = (TCanvas*)gROOT->FindObject("C2")))
879         c2 = new TCanvas("C2","ITS SSD Cone Geometry",900,450);
880     c2->Divide(2,1);
881     TGeoVolume *ITS,*SSD=0;
882     TGeoNode *node;
883     //
884     node = ALIC->FindNode("ITSV_1");
885     ITS = node->GetVolume();
886     node = ITS->FindNode("ITSssdConeA_1");
887     SSD = node->GetVolume();
888     //
889     mgr2->SetNsegments(ISetits(1,-1));
890     //
891     mgr2->SetVisLevel(6);
892     mgr2->SetVisOption(0);
893     //mgr2->CheckOverlaps(0.01);
894     //mgr2->PrintOverlaps();
895     if(ISetits(2,-1)==1){
896         TGeoShape *clip = new TGeoTubeSeg(0, 1000, 2000, 45, 90);
897         mgr2->SetClippingShape(clip);
898     } // end if
899     mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
900     if(ISetits(2,-1)!=0) mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
901     //
902     c2->cd(1);
903     SSD->Draw();
904     TPad *p1 = c2->GetPad(1);
905     TView *view1 = p1->GetView();
906     if(view1){
907         view1->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
908         if(irr) cout <<"error="<<irr<<endl;
909         if(ISetits(4,-1)==0) view1->SetParallel();
910         else  view1->SetPerspective();
911         view1->Top();
912         if(ISetits(3,-1)!=0) view1->ShowAxis();
913     } // end if view1
914     if(ISetits(5,-1)==1) SSD->Raytrace();
915     //
916     c2->cd(2);
917     SSD->Draw();
918     TPad *p2 = c2->GetPad(2);
919     TView *view2 = p2->GetView();
920     if(view2){
921         view2->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
922         if(irr) cout <<"error="<<irr<<endl;
923         if(ISetits(4,-1)==0) view2->SetParallel();
924         else  view2->SetPerspective();
925         view2->Front();
926         if(ISetits(3,-1)!=0) view2->ShowAxis();
927     } // end if view1
928     if(ISetits(5,-1)==1) SSD->Raytrace();
929     //
930 }
931 //----------------------------------------------------------------------
932 void EngineeringSSDCylinder(){
933     // Display SSD Cylinder Geometry
934     // Inputs:
935     //    none.
936     // Outputs:
937     //    none.
938     // Retrurn:
939     //    none.
940     Int_t irr;
941     //
942     TGeoManager *mgr2 = gGeoManager;
943     TGeoVolume *ALIC = mgr2->GetTopVolume();
944     TCanvas *c3;
945     if(!(c3 = (TCanvas*)gROOT->FindObject("C3")))
946         c3 = new TCanvas("C3","ITS SSD Cylinder Geometry",900,450);
947     c3->Divide(2,1);
948     TGeoVolume *ITS,*SSD=0;
949     TGeoNode *node;
950     TArrow *arrow=new TArrow();
951     //
952     node = ALIC->FindNode("ITSV_1");
953     ITS = node->GetVolume();
954     node = ITS->FindNode("ITSssdCentCylCA_1");
955     SSD = node->GetVolume();
956     //Double_t Rmin = ((TGeoPcon*)(SSD->GetShape()))->GetRmin();
957     //Double_t Rmax = ((TGeoPcon*)(SSD->GetShape()))->GetRmax();
958     //Double_t Dz   = ((TGeoPcon*)(SSD->GetShape()))->GetDz();
959     //
960     mgr2->SetNsegments(ISetits(1,-1));
961     //
962     mgr2->SetVisLevel(6);
963     mgr2->SetVisOption(0);
964     //mgr2->CheckOverlaps(0.01);
965     //mgr2->PrintOverlaps();
966     if(ISetits(2,-1)==1){
967         TGeoShape *clip = new TGeoTubeSeg(0, 1000, 2000, 45, 90);
968         mgr2->SetClippingShape(clip);
969     } // end if
970     mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
971     if(ISetits(2,-1)!=0) mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
972     //
973     c3->cd(1);
974     SSD->Draw();
975     TPad *p1 = c3->GetPad(1);
976     TView *view1 = p1->GetView();
977     if(view1){
978         view1->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
979         if(irr) cout <<"error="<<irr<<endl;
980         if(ISetits(4,-1)==0) view1->SetParallel();
981         else  view1->SetPerspective();
982         view1->Front();
983         if(ISetits(3,-1)!=0) view1->ShowAxis();
984     } // end if view1
985     if(ISetits(5,-1)==1) SSD->Raytrace();
986     //arrow->DrawArrow(1.01*Rmax,-Dz,1.01*Rmax,+Dz);
987     //
988     c3->cd(2);
989     SSD->Draw();
990     TPad *p2 = c3->GetPad(2);
991     TView *view2 = p2->GetView();
992     if(view2){
993         view2->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
994         if(irr) cout <<"error="<<irr<<endl;
995         if(ISetits(4,-1)==0) view2->SetParallel();
996         else  view2->SetPerspective();
997         view2->Top();
998         if(ISetits(3,-1)!=0) view2->ShowAxis();
999     } // end if view2
1000     if(ISetits(5,-1)==1) SSD->Raytrace();
1001     //arrow->DrawArrow(Rmax,0.0,Rmax,0.0);
1002     //Double_t s = TMath::Sin(0.7),c = TMath::Cos(0.7);
1003     //arrow->DrawArrow(-Rmin*c,-Rmin*s,Rmin*c,Rmin*s);
1004     //
1005 }
1006 //----------------------------------------------------------------------
1007 void EngineeringSupRB24(){
1008     // Display  RB 24 side cable tray support structure Geometry
1009     // Inputs:
1010     //    none.
1011     // Outputs:
1012     //    none.
1013     // Retrurn:
1014     //    none.
1015     Int_t irr;
1016     //
1017     TGeoManager *mgr2 = gGeoManager;
1018     TGeoVolume *ALIC = mgr2->GetTopVolume();
1019     TCanvas *c4;
1020     if(!(c4 = (TCanvas*)gROOT->FindObject("C4")))
1021         c4 = new TCanvas("C4","ITS SDD Cylinder Geometry",900,450);
1022     c4->Divide(2,1);
1023     TGeoVolume *ITS,*SUPRB24=0;
1024     TGeoNode *node;
1025     TArrow *arrow=new TArrow();
1026     //
1027     node = ALIC->FindNode("ITSV_1");
1028     ITS = node->GetVolume();
1029     node = ITS->FindNode("ITSsupFrameM24_1");
1030     SUPRB24 = node->GetVolume();
1031     //
1032     mgr2->SetNsegments(ISetits(1,-1));
1033     //
1034     mgr2->SetVisLevel(6);
1035     mgr2->SetVisOption(0);
1036     //mgr2->CheckOverlaps(0.01);
1037     //mgr2->PrintOverlaps();
1038     if(ISetits(2,-1)==1){
1039         TGeoShape *clip = new TGeoTubeSeg(0, 1000, 2000, 45, 90);
1040         mgr2->SetClippingShape(clip);
1041     } // end if
1042     mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
1043     if(ISetits(2,-1)!=0) mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
1044     //
1045     c4->cd(1);
1046     SUPRB24->Draw();
1047     TPad *p1 = c4->GetPad(1);
1048     TView *view1 = p1->GetView();
1049     if(view1){
1050         view1->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
1051         if(irr) cout <<"error="<<irr<<endl;
1052         if(ISetits(4,-1)==0) view1->SetParallel();
1053         else  view1->SetPerspective();
1054         view1->Front();
1055         if(ISetits(3,-1)!=0) view1->ShowAxis();
1056     } // end if view1
1057     if(ISetits(5,-1)==1) SUPRB24->Raytrace();
1058     //
1059     c4->cd(2);
1060     SUPRB24->Draw();
1061     TPad *p2 = c4->GetPad(2);
1062     TView *view2 = p2->GetView();
1063     if(view2){
1064         view2->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
1065         if(irr) cout <<"error="<<irr<<endl;
1066         if(ISetits(4,-1)==0) view2->SetParallel();
1067         else  view2->SetPerspective();
1068         view2->Top();
1069         if(ISetits(3,-1)!=0) view2->ShowAxis();
1070     } // end if view2
1071     if(ISetits(5,-1)==1) SUPRB24->Raytrace();
1072     //
1073 }
1074 //----------------------------------------------------------------------
1075 void EngineeringSupTrayRB24(){
1076     // Display  RB 24 side cable tray support structure Geometry
1077     // Inputs:
1078     //    none.
1079     // Outputs:
1080     //    none.
1081     // Retrurn:
1082     //    none.
1083     Int_t irr;
1084     //
1085     TGeoManager *mgr2 = gGeoManager;
1086     TGeoVolume *ALIC = mgr2->GetTopVolume();
1087     TCanvas *c4,*c5;
1088     //if(!(c4 = (TCanvas*)gROOT->FindObject("C4")))
1089         c4 = new TCanvas("C4","ITS RB24 Cable Trays and Patch Pannels",500,500);
1090     //c4->Divide(2,1);
1091     TGeoVolume *ITS,*SUPRB24=0;
1092     TGeoNode *node;
1093     TArrow *arrow=new TArrow();
1094     //
1095     node = ALIC->FindNode("ITSV_1");
1096     ITS = node->GetVolume();
1097     node = ITS->FindNode("ITSsupCableTrayMotherMT24_1");
1098     SUPRB24 = node->GetVolume();
1099     //
1100     mgr2->SetNsegments(ISetits(1,-1));
1101     //
1102     mgr2->SetVisLevel(6);
1103     mgr2->SetVisOption(0);
1104     //mgr2->CheckOverlaps(0.01);
1105     //mgr2->PrintOverlaps();
1106     if(ISetits(2,-1)==1){
1107         TGeoShape *clip = new TGeoTubeSeg(0, 1000, 2000, 45, 90);
1108         mgr2->SetClippingShape(clip);
1109     } // end if
1110     mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
1111     if(ISetits(2,-1)!=0) mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
1112     //
1113     c4->cd(1);
1114     SUPRB24->Draw();
1115     //TPad *p1 = c4->GetPad(1);
1116     //TView *view1 = p1->GetView();
1117     TView *view1 = c4->GetView();
1118     if(view1){
1119         view1->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
1120         if(irr) cout <<"error="<<irr<<endl;
1121         if(ISetits(4,-1)==0) view1->SetParallel();
1122         else  view1->SetPerspective();
1123         view1->Front();
1124         if(ISetits(3,-1)!=0) view1->ShowAxis();
1125     } // end if view1
1126     if(ISetits(5,-1)==1) SUPRB24->Raytrace();
1127     //
1128     //c4->cd(2);
1129         c5 = new TCanvas("C5","ITS RB24 Cable Trays and Patch Pannels",500,500);
1130     c5->cd(1);
1131     SUPRB24->Draw();
1132     //TPad *p2 = c5->GetPad(1);
1133     //TView *view2 = p2->GetView();
1134     TView *view2 = c5->GetView();
1135     if(view2){
1136         view2->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
1137         if(irr) cout <<"error="<<irr<<endl;
1138         if(ISetits(4,-1)==0) view2->SetParallel();
1139         else  view2->SetPerspective();
1140         view2->Top();
1141         if(ISetits(3,-1)!=0) view2->ShowAxis();
1142     } // end if view2
1143     if(ISetits(5,-1)==1) SUPRB24->Raytrace();
1144     //
1145 }
1146 //----------------------------------------------------------------------
1147 void EngineeringSupRB26(){
1148     // Display RB 26 side cable tray support structure
1149     // Inputs:
1150     //    none.
1151     // Outputs:
1152     //    none.
1153     // Retrurn:
1154     //    none.
1155     Int_t irr;
1156     //
1157     TGeoManager *mgr2 = gGeoManager;
1158     TGeoVolume *ALIC = mgr2->GetTopVolume();
1159     TCanvas *c5;
1160     if(!(c5 = (TCanvas*)gROOT->FindObject("C5")))
1161         c5 = new TCanvas("C5","ITS SDD Cylinder Geometry",900,450);
1162     c5->Divide(2,1);
1163     TGeoVolume *ITS,*SUPRB26=0;
1164     TGeoNode *node;
1165     TArrow *arrow=new TArrow();
1166     //
1167     node = ALIC->FindNode("ITSV_1");
1168     ITS = node->GetVolume();
1169     node = ITS->FindNode("ITSsupFrameM26_1");
1170     SUPRB26 = node->GetVolume();
1171     //
1172     mgr2->SetNsegments(ISetits(1,-1));
1173     //
1174     mgr2->SetVisLevel(6);
1175     mgr2->SetVisOption(0);
1176     //mgr2->CheckOverlaps(0.01);
1177     //mgr2->PrintOverlaps();
1178     if(ISetits(2,-1)==1){
1179         TGeoShape *clip = new TGeoTubeSeg(0, 1000, 2000, 45, 90);
1180         mgr2->SetClippingShape(clip);
1181     } // end if
1182     mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
1183     if(ISetits(2,-1)!=0) mgr2->SetPhiRange(DSetits(1,-1.),DSetits(0,-1.));
1184     //
1185     c5->cd(1);
1186     SUPRB26->Draw();
1187     TPad *p1 = c5->GetPad(1);
1188     TView *view1 = p1->GetView();
1189     if(view1){
1190         view1->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
1191         if(irr) cout <<"error="<<irr<<endl;
1192         if(ISetits(4,-1)==0) view1->SetParallel();
1193         else  view1->SetPerspective();
1194         view1->Front();
1195         if(ISetits(3,-1)!=0) view1->ShowAxis();
1196     } // end if view1
1197     if(ISetits(5,-1)==1) SUPRB26->Raytrace();
1198     //
1199     c5->cd(2);
1200     SUPRB26->Draw();
1201     TPad *p2 = c5->GetPad(2);
1202     TView *view2 = p2->GetView();
1203     if(view2){
1204         view2->SetView(DSetits(2,-1.),DSetits(3,-1.),DSetits(4,-1.),irr);
1205         if(irr) cout <<"error="<<irr<<endl;
1206         if(ISetits(4,-1)==0) view2->SetParallel();
1207         else  view2->SetPerspective();
1208         view2->Top();
1209         if(ISetits(3,-1)!=0) view2->ShowAxis();
1210     } // end if view2
1211     if(ISetits(5,-1)==1) SUPRB26->Raytrace();
1212     //
1213 }