Start of a TTask to replace Displayv11.C

diff --git a/ITS/DisplayITSv11.cxx b/ITS/DisplayITSv11.cxx
new file mode 100644 (file)
index 0000000..ab2cb7d
--- /dev/null
+++ b/ITS/DisplayITSv11.cxx
@@ -0,0 +1,452 @@
+/**************************************************************************
+ * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ *                                                                        *
+ * Author: The ALICE Off-line Project.                                    *
+ * Contributors are mentioned in the code where appropriate.              *
+ *                                                                        *
+ * Permission to use, copy, modify and distribute this software and its   *
+ * documentation strictly for non-commercial purposes is hereby granted   *
+ * without fee, provided that the above copyright notice appears in all   *
+ * copies and that both the copyright notice and this permission notice   *
+ * appear in the supporting documentation. The authors make no claims     *
+ * about the suitability of this software for any purpose. It is          *
+ * provided "as is" without express or implied warranty.                  *
+ **************************************************************************/
+
+/*
+$Id$
+*/
+#include <TSystem.h>
+#include <TROOT.h>
+#include <TGeoManager.h>
+#include <TGeoNode.h>
+#include <TGeoVolume.h>
+#include <TGeoShape.h>
+#include <TGeoTube.h>
+#include <TGeoMaterial.h>
+#include <TGeoMedium.h>
+#include <TCanvas.h>
+#include <TControlBar.h>
+#include <TArrow.h>
+#include <TPad.h>
+#include <TView.h>
+
+#include "AliITS.h"
+#include "AliITSv11.h"
+#include "DisplayITSv11.h"
+
+ClassImp(DisplayITSv11)
+//______________________________________________________________________
+DisplayITSv11::DisplayITSv11(){
+    // Default Constructor.
+    // Inputs:
+    //    none.
+    // Outputs:
+    //    none.
+    // Return:
+    //    A Default constructed DispalyITSv11 class/task.
+
+    fmgr         = 0;
+    fits         = 0;
+    fALICE       = 0;
+    fITS         = 0;
+    fClip        = 0;
+    fITSdebug    = 0;
+    fNsegments   = 80;
+    fCut         = 0;
+    fAxis        = 1;
+    fPerspective = 0;
+    fSolid       = 0;
+    fPhimincut   = 0.0;
+    fPhimaxcut   = 180.0;
+    fRmin[0] = fRmin[1] = fRmin[2] = fRmax[0] = fRmax[1] = fRmax[2] = 0.0;
+    fLongitude   = 90.0;
+    fLatitude    = 0.0;
+    fPsi         = 0.0;
+}
+//______________________________________________________________________
+DisplayITSv11::~DisplayITSv11(){
+    // Destructor.
+    // Inputs:
+    //    none.
+    // Outputs:
+    //    none.
+    // Return:
+    //    none..
+
+    if(fits) delete fits;
+    if(fClip) delete fClip;
+    fITS   = 0;
+    fALICE = 0;
+    fmgr   = 0;
+}
+//______________________________________________________________________
+void DisplayITSv11::Exec(Option_t* opt){
+    // Main display routine
+    // Inputs:
+    //      Option_t * opt   Not presently used.
+    // Outputs:
+    //      none.
+    // Return:
+    //      none.
+    Int_t i;
+
+    gSystem->Load("libGeom");
+    //
+    if(gGeoManager) delete gGeoManager;
+    fmgr = gGeoManager = new TGeoManager("ITSGeometry",
+                                         " ITS Simulation Geometry Manager");
+    //
+    TGeoMaterial *vacmat = new TGeoMaterial("Vacume",0,0,0);
+    TGeoMedium   *vacmed = new TGeoMedium("Vacume_med",1,vacmat);
+    fRmax[0] = 100.;
+    fRmax[1] = 100.;
+    fRmax[2] = 200.;
+    for(i=0;i<3;i++) fRmin[i] = -fRmax[i];
+    fALICE = fmgr->MakeBox("ALIC",vacmed,fRmax[0],fRmax[1],fRmax[2]);
+    fmgr->SetTopVolume(fALICE);
+    //
+    AliITSv11 *fits = new AliITSv11();
+    fits->SetDebug(GetDebugITS());
+    fits->CreateMaterials();
+    fits->CreateGeometry();
+    //
+    fmgr->CloseGeometry();
+    fITS = fALICE->FindNode("ITSV_1")->GetVolume();
+    //
+    TControlBar *bar=new TControlBar("vertical","ITS Geometry Display",10,10);
+    bar->AddButton("Set ITS Debug level 1","SetITSdebugOn()","Debug on");
+    bar->AddButton("Set ITS Debug level 0","SetITSdebugOff()","Debug off");
+    bar->AddButton("Set ITS Theta,Phi cut on","SetCylindericalClipVolume()",
+                   "Cut on");
+    bar->AddButton("Set ITS Theta,Phi cut off","SetCylindericalCutOff()",
+                   "Cut off");
+    bar->AddButton("Set axis on","SetAxisOn()","Show Axis on");
+    bar->AddButton("Set axis off","SetAxisOff()","Show Axis off");
+    bar->AddButton("Set perspective on","SetPerspectiveOn()",
+                   "Perspective on");
+    bar->AddButton("Set perspective off","SetPerspectiveOff()",
+                   "Perspective off");
+    bar->AddButton("Draw volumes as Solid","SetSolid()","Solid Volumes");
+    bar->AddButton("Draw volumes as wire","SetWire()","Wire Volumes");
+    bar->AddButton("Set circle/80","SetCircleSegments(80)",
+                   "circles ~ by 80 lines");
+    bar->AddButton("Check Overlaps","fmgr->CheckOverlaps()",
+                   "Check for overlaps");
+    bar->AddButton("Display Geometry","Displayit()","Run Displayit");
+    bar->AddButton("Display SPD Thermal Sheald","EngineeringSPDThS()",
+                   "Run EngineeringSPDThS");
+    bar->AddButton("Display SDD Cone","EngineeringSDDCone()",
+                   "Run EngineeringSDDCone");
+    bar->AddButton("Display SDD Centeral Cylinder","EngineeringSDDCylinder()",
+                   "Run EngineeringSDDCylinder");
+    bar->AddButton("Display SUP RB24 side","EngineeringSupRB24()",
+                   "Run EngineeringSDDCylinder");
+    bar->AddButton("Display SUP RB26 side","EngineeringSupRB26()",
+                   "Run EngineeringSupRB26");
+    bar->AddButton("Quit/Exit",".q","Exit");
+    bar->Show();
+    gROOT->SaveContext();
+         //Displayit();
+}
+//______________________________________________________________________
+void DisplayITSv11::Displaying(TGeoVolume *v,TCanvas *c,Int_t ipad){
+    // Display Volume according to existing values
+    // Inputs:
+    //    TGeoVolume *v volume to be drawn.
+    //    TCanvas    *p Pad where drawing is to be done
+    //    Int_t      ipad subpad to draw on
+    // Outputs:
+    //    none.
+    // Return:
+    //    none.
+    Int_t irr;
+    TVirtualPad *p = 0;
+    TView *view=0;
+
+    c->cd(ipad);
+    p = c->GetPad(ipad);
+    if(!p) return;
+    view = p->GetView();
+    if(!view){view = new TView(fRmin,fRmax,1);p->SetView(view);}
+    view->SetRange(fRmin,fRmax);
+    view->SetView(fLongitude,fLatitude,fPsi,irr);
+    if(fPerspective) view->SetPerspective();
+    else view->SetParralel();
+    if(fAxis) view->ShowAxis();
+    //
+    if(fClip) fmgr->SetClippingShape(fClip);
+    v->Draw();
+    if(fSolid) v->Raytrace();
+    p->Modified(1);
+    p->Update();
+}
+//----------------------------------------------------------------------
+void DisplayITSv11::DisplayITS(){
+    // Display AliITSv11 Geometry
+    // Inputs:
+    //    const char* filename output file with the display in it
+    // Outputs:
+    //    none.
+    // Retrurn:
+    //    none.
+    Double_t lon,lat,psi;
+    //
+    TCanvas *c1;
+    if(!(c1 = (TCanvas*)gROOT->FindObject("C1")))
+        c1 = new TCanvas("C1","ITS Simulation Geometry",900,900);
+    c1->Divide(2,2);
+    //
+    fmgr->SetNsegments(fNsegments);
+    //
+    fmgr->SetVisLevel(6);
+    fmgr->SetVisOption(0);
+    // Perspective
+    Displaying(fALICE,c1,2);
+    lon = fLongitude;
+    lat = fLatitude;
+    psi = fPsi;
+    fLongitude = 270.0; // Front
+    fLatitude  =  90.0;
+    fPsi       =   0.0;
+    Displaying(fALICE,c1,1);
+    fLongitude = 270.0; // Top
+    fLatitude  =   0.0;
+    fPsi       =   0.0;
+    Displaying(fALICE,c1,3);
+    fLongitude =   0.0; // Side
+    fLatitude  =  90.0;
+    fPsi       =   0.0;
+    Displaying(fALICE,c1,4);
+    fLongitude = lon;
+    fLatitude  = lat;
+    fPsi       = psi;
+    //
+}
+//----------------------------------------------------------------------
+void DisplayITSv11::EngineeringSPDThS(){
+    // Display SPD Thermal Sheald Geometry
+    // Inputs:
+    //    none.
+    // Outputs:
+    //    none.
+    // Retrurn:
+    //    none.
+    Double_t lon,lat,psi;
+    //
+    TCanvas *c4;
+    TGeoNode *node;
+    if(!(c4 = (TCanvas*)gROOT->FindObject("C4")))
+        c4 = new TCanvas("C4","ITS SDD Cylinder Geometry",900,450);
+    c4->Divide(2,1);
+    TGeoVolume *SPDThS=0;
+         //TArrow *arrow=new TArrow();
+    //
+    node = fITS->FindNode("ITSspdThermalSheald_1");
+    SPDThS = node->GetVolume();
+    //
+    fmgr->SetNsegments(fNsegments);
+    //
+    fmgr->SetVisLevel(6);
+    fmgr->SetVisOption(0);
+    //
+    lon = fLongitude;
+    lat = fLatitude;
+    psi = fPsi;
+    fLongitude = 270.0; // Front
+    fLatitude  =  90.0;
+    fPsi       =   0.0;
+    Displaying(SPDThS,c4,1);
+    //
+    fLongitude = 270.0; // Top
+    fLatitude  =   0.0;
+    fPsi       =   0.0;
+    Displaying(SPDThS,c4,2);
+    //
+    fLongitude = lon;
+    fLatitude  = lat;
+    fPsi       = psi;
+}
+//----------------------------------------------------------------------
+void DisplayITSv11::EngineeringSDDCone(){
+    // Display SDD Cone Geometry
+    // Inputs:
+    //    none.
+    // Outputs:
+    //    none.
+    // Retrurn:
+    //    none.
+    Double_t lon,lat,psi;
+    //
+    TCanvas *c2;
+    if(!(c2 = (TCanvas*)gROOT->FindObject("C2")))
+        c2 = new TCanvas("C2","ITS SDD Cone Geometry",900,450);
+    c2->Divide(2,1);
+    TGeoVolume *SDD=0;
+    TGeoNode *node;
+    //
+    node = fITS->FindNode("ITSsddConeL_1");
+    SDD = node->GetVolume();
+    //
+    fmgr->SetNsegments(fNsegments);
+    //
+    fmgr->SetVisLevel(6);
+    fmgr->SetVisOption(0);
+    //
+    lon = fLongitude;
+    lat = fLatitude;
+    psi = fPsi;
+    fLongitude = 270.0; // Front
+    fLatitude  =  90.0;
+    fPsi       =   0.0;
+    Displaying(SDD,c2,1);
+    //
+    fLongitude = 270.0; // Top
+    fLatitude  =   0.0;
+    fPsi       =   0.0;
+    Displaying(SDD,c2,2);
+    //
+    fLongitude = lon;
+    fLatitude  = lat;
+    fPsi       = psi;
+}
+//----------------------------------------------------------------------
+void DisplayITSv11::EngineeringSDDCylinder(){
+    // Display SDD Cylinder Geometry
+    // Inputs:
+    //    none.
+    // Outputs:
+    //    none.
+    // Retrurn:
+    //    none.
+    Double_t lon,lat,psi;
+    //
+    TCanvas *c3;
+    if(!(c3 = (TCanvas*)gROOT->FindObject("C3")))
+        c3 = new TCanvas("C3","ITS SDD Cylinder Geometry",900,450);
+    c3->Divide(2,1);
+    TGeoVolume *SDD=0;
+    TGeoNode *node;
+    TArrow *arrow=new TArrow();
+    //
+    node = fITS->FindNode("ITSsddCentCylCF_1");
+    SDD = node->GetVolume();
+    Double_t Rmin = ((TGeoTube*)(SDD->GetShape()))->GetRmin();
+    Double_t Rmax = ((TGeoTube*)(SDD->GetShape()))->GetRmax();
+    Double_t Dz   = ((TGeoTube*)(SDD->GetShape()))->GetDz();
+    //
+    fmgr->SetNsegments(fNsegments);
+    //
+    fmgr->SetVisLevel(6);
+    fmgr->SetVisOption(0);
+    //
+    lon = fLongitude;
+    lat = fLatitude;
+    psi = fPsi;
+    fLongitude = 270.0; // Front
+    fLatitude  =  90.0;
+    fPsi       =   0.0;
+    Displaying(SDD,c3,1);
+    arrow->DrawArrow(1.01*Rmax,-Dz,1.01*Rmax,+Dz);
+    //
+    fLongitude = 270.0; // Top
+    fLatitude  =   0.0;
+    fPsi       =   0.0;
+    Displaying(SDD,c3,2);
+    arrow->DrawArrow(Rmax,0.0,Rmax,0.0);
+    Double_t s = TMath::Sin(0.7),c = TMath::Cos(0.7);
+    arrow->DrawArrow(-Rmin*c,-Rmin*s,Rmin*c,Rmin*s);
+    //
+    fLongitude = lon;
+    fLatitude  = lat;
+    fPsi       = psi;
+}
+//----------------------------------------------------------------------
+void DisplayITSv11::EngineeringSupRB24(){
+    // Display SDD Cylinder Geometry
+    // Inputs:
+    //    none.
+    // Outputs:
+    //    none.
+    // Retrurn:
+    //    none.
+    Double_t lon,lat,psi;
+    //
+    TCanvas *c4;
+    if(!(c4 = (TCanvas*)gROOT->FindObject("C4")))
+        c4 = new TCanvas("C4","ITS SDD Cylinder Geometry",900,450);
+    c4->Divide(2,1);
+    TGeoVolume *SUPRB24=0;
+    TGeoNode *node;
+         //TArrow *arrow=new TArrow();
+    //
+    node = fITS->FindNode("ITSsupFrameM24_1");
+    SUPRB24 = node->GetVolume();
+    //
+    fmgr->SetNsegments(fNsegments);
+    //
+    fmgr->SetVisLevel(6);
+    fmgr->SetVisOption(0);
+    //
+    lon = fLongitude;
+    lat = fLatitude;
+    psi = fPsi;
+    fLongitude = 270.0; // Front
+    fLatitude  =  90.0;
+    fPsi       =   0.0;
+    Displaying(SUPRB24,c4,1);
+    //
+    fLongitude = 270.0; // Top
+    fLatitude  =   0.0;
+    fPsi       =   0.0;
+    Displaying(SUPRB24,c4,2);
+    //
+    fLongitude = lon;
+    fLatitude  = lat;
+    fPsi       = psi;
+}
+
+//----------------------------------------------------------------------
+void DisplayITSv11::EngineeringSupRB26(){
+    // Display SDD Cylinder Geometry
+    // Inputs:
+    //    none.
+    // Outputs:
+    //    none.
+    // Retrurn:
+    //    none.
+    Double_t lon,lat,psi;
+    //
+    TCanvas *c5;
+    if(!(c5 = (TCanvas*)gROOT->FindObject("C5")))
+        c5 = new TCanvas("C5","ITS SDD Cylinder Geometry",900,450);
+    c5->Divide(2,1);
+    TGeoVolume *SUPRB26=0;
+    TGeoNode *node;
+         //TArrow *arrow=new TArrow();
+    //
+    node = fITS->FindNode("ITSsupFrameM26_1");
+    SUPRB26 = node->GetVolume();
+    //
+    fmgr->SetNsegments(fNsegments);
+    //
+    fmgr->SetVisLevel(6);
+    fmgr->SetVisOption(0);
+    //
+    lon = fLongitude;
+    lat = fLatitude;
+    psi = fPsi;
+    fLongitude = 270.0; // Front
+    fLatitude  =  90.0;
+    fPsi       =   0.0;
+    Displaying(SUPRB26,c5,1);
+    //
+    fLongitude = 270.0; // Top
+    fLatitude  =   0.0;
+    fPsi       =   0.0;
+    Displaying(SUPRB26,c5,2);
+    //
+    fLongitude = lon;
+    fLatitude  = lat;
+    fPsi       = psi;
+}

diff --git a/ITS/DisplayITSv11.h b/ITS/DisplayITSv11.h
new file mode 100644 (file)
index 0000000..a04418a
--- /dev/null
+++ b/ITS/DisplayITSv11.h
@@ -0,0 +1,72 @@
+#ifndef DISPLAYITSV11_H
+#define DISPLAYITSV11_H
+/* Copyright (c) 1998-2001, ALICE Experiment at CERN, All rights reserved *
+ * See cxx source for full Copyright notice                               */
+
+/*
+  $Id$
+ */
+#include <TTask.h>
+#include <TGeoTube.h>
+
+class TGeoManager;
+class AliITS;
+class TGeoVolume;
+class TCanvas;
+
+class DisplayITSv11 : public TTask {
+  public:
+    DisplayITSv11();
+    virtual ~DisplayITSv11();
+    //
+    virtual void Exec(Option_t* opt="");
+    virtual void DisplayITS();
+    virtual void EngineeringSPDThS();
+    virtual void EngineeringSDDCone();
+    virtual void EngineeringSDDCylinder();
+    virtual void EngineeringSupRB24();
+    virtual void EngineeringSupRB26();
+    //
+    virtual Int_t GetDebugITS(){return fITSdebug;}
+    //
+    void SetITSdebugOn(){fITSdebug=1;}
+    void SetITSdebugOff(){fITSdebug=0;}
+    void SetCircleSegments(Int_t i=80){fNsegments=i;}
+    void SetCylindericalCutOn(Double_t phimin=0.0,Double_t phimax=180.){
+        fPhimincut = phimin;fPhimaxcut=phimax;fCut=1;}
+    void SetCylindericalCutOff(){fCut=0;if(fClip) delete fClip;fClip=0;}
+    void SetCylindericalClipVolume(){if(fClip)delete fClip;
+         fClip = new TGeoTubeSeg(0.0,fRmax[0],fRmax[2],fPhimincut,fPhimaxcut);}
+    void SetPerspectiveOn(){fPerspective=1;}
+    void SetPerspectiveOff(){fPerspective=0;}
+    void SetSolid(){fSolid=1;}
+    void SetWire(){fSolid=0;}
+    void SetAxisOn(){fAxis=1;}
+    void SetAxisOff(){fAxis=0;}
+    void SetDisplayAngles(Double_t lon,Double_t lat,Double_t psi){
+        fLongitude=lon;fLatitude=lat;fPsi=psi;}
+  private:
+    void Displaying(TGeoVolume *v,TCanvas *c,Int_t ipad);
+    TGeoManager *fmgr;  // pointer to the geometry manager.
+    AliITS      *fits;  // pointer to AliITS
+    TGeoVolume *fALICE; // Pointer to the true mother volume
+    TGeoVolume *fITS;   // Pointer to the ITS mother volume
+    TGeoShape  *fClip;  // Clipping Volume.
+    //
+    Int_t fITSdebug;
+    Int_t fNsegments;
+    Int_t fCut;
+    Int_t fAxis;
+    Int_t fPerspective;
+    Int_t fSolid;
+    Double_t fRmin[3],fRmax[3];
+    Double_t fPhimincut;
+    Double_t fPhimaxcut;
+    Double_t fLongitude;
+    Double_t fLatitude;
+    Double_t fPsi;
+
+    ClassDef(DisplayITSv11,1) // Task to display ITS v11 Geometry
+};
+
+#endif