Updated error definitions.

[u/mrichter/AliRoot.git] / ITS / AliITSBaseGeometry.cxx

diff --git a/ITS/AliITSBaseGeometry.cxx b/ITS/AliITSBaseGeometry.cxx
index 03044848f1f47bb7eaa1808bad46af2bfe2ac665..bb4267feebfe995e7e35acaea3104cc112088c84 100644 (file)
--- a/ITS/AliITSBaseGeometry.cxx
+++ b/ITS/AliITSBaseGeometry.cxx
@@ -19,36 +19,10 @@
 #include <TObject.h>
 #include <TObjArray.h>
 
-//#include <TGeoArb8.h>
-//#include <TGeoEltu.h>
-//#include <TGeoOverlap.h>
-//#include <TGeoTrack.h>
-//#include <TGeoAtt.h>
-//#include <TGeoManager.h>
-//#include <TGeoPainter.h>
-//#include <TGeoTrd1.h>
-//#include <TGeoBBox.h>
 #include <TGeoMaterial.h>
-//#include <TGeoPara.h>
-//#include <TGeoTrd2.h>
-//#include <TGeoBoolNode.h>
 #include <TGeoMatrix.h>
-//#include <TGeoPatternFinder.h>
-//#include <TGeoTube.h>
-//#include <TGeoCache.h>
-//#include <TGeoMCGeometry.h>
-//#include <TGeoPcon.h>
-//#include <TGeoVolume.h>
-//#include <TGeoChecker.h>
 #include <TGeoMedium.h>
 #include <TGeoPgon.h>
-//#include <TGeoVoxelFinder.h>
-//#include <TGeoCompositeShape.h>
-//#include <TGeometry.h>
-//#include <TGeoShape.h>
-//#include <TGeoCone.h>
-//#include <TGeoNode.h>
-//#include <TGeoSphere.h>
 #include "AliITSBaseGeometry.h"
 
 
@@ -150,21 +124,103 @@ AliITSMixture::AliITSMixture(const char *name,Int_t N,Double_t *w,TObjArray *m,
     } // end for i
     delete[] nw;
 }
+//======================================================================
+ClassImp(AliITSGeoCable)
+;
+AliITSGeoCable::AliITSGeoCable(){
+    //
 
-//ClassImp(AliITSArb8)
-
-//ClassImp(AliITSBBox)
-
-//ClassImp(AliITSCone)
-
-//ClassImp(AliITSConeSeg)
-
-//ClassImp(AliITSCtub)
-
-//ClassImp(AliITSEltu)
-
-//ClassImp(AliITSGtra)
+    fRmin = fRmax = 0.0;
+    fNs.SetXYZ(0.0,0.0,0.0);
+    fNe.SetXYZ(0.0,0.0,0.0);
+    fTubes    = 0;
+    fTranRot  = 0;
+}
+//----------------------------------------------------------------------
+AliITSGeoCable::AliITSGeoCable(const char *name,const TObjArray *vect,
+                               const Double_t Rmin,const Double_t Rmax,
+                               const TVector3 ns,const TVector3 ne){
+    //
+    // Inputs:
+    //    char      *name   Name of this compound object
+    //    TObjArray *vect   Array of TVector3's of points representing the
+    //                      path of the cable
+    //    TVector3  ns=0    Normal vector representing the angle of the 
+    //                      starting surface, default perpendicular
+    //    TVector3  ne=0    Normal vector representing the angle of the 
+    //                      ending surface, default perpendicular
+    // Outputs:
+    //    none.
+    // Return:
+    //    A fully initilized and created AliITSGeoCable class.
+    Char_t nam[500];
+    Int_t i,n;
+    Double_t s,th,ph;
+    TVector3 x0,x1,x2,d,t,n0,n1;
 
-//ClassImp(AliITSpCon)
+    fRmin = Rmin;
+    fRmax = Rmax;
+    fRmin = fRmax = 0.0;
+    fNs.SetXYZ(0.0,0.0,-1.0);
+    fNe.SetXYZ(0.0,0.0,1.0);
+    n = vect->GetEntries();
+    fTubes = new TObjArray(n-1);
+    fTranRot  = new TObjArray(n-1);
+    fNs = ns*(1./ns.Mag());
+    fNe = ne*(1./ne.Mag());
+    //
+    x0 = *((TVector3 *)(vect->At(0)));
+    n0 = ns;
+    for(i=1;i<n;i++){
+        x1 = *((TVector3 *)(vect->At(i)));
+        d  = x1 - x0;
+        if(i<n-1) {
+            x2 = *((TVector3 *)(vect->At(i+1)));
+            n1 = d + (x2-x1);
+            n1 *= 1./n1.Mag();
+        }else{
+            n1 = fNe;
+        } // end if
+        t  = 0.5*(x1 + x0);
+        th =  TMath::ATan2(TMath::Sqrt(2.*d.Mag2()-d.z()-2.*d.Mag()*d.y()-
+                                        2.*d.Mag()*d.x()),
+                           TMath::Sqrt(d.z()*d.z()-2.*d.Mag()*d.z()+d.Mag2()));
+        th *= TMath::RadToDeg();
+        ph =  TMath::ATan2(d.y()-d.Mag(),d.x()-d.Mag());
+        ph *= TMath::RadToDeg();
+        sprintf(nam,"%sCombiTrans%dCable",name,i-1);
+        fTranRot->AddAt(new TGeoCombiTrans(nam,t.x(),t.y(),t.z(),
+                                        new TGeoRotation("",ph,th,0.0)),i-1);
+        s  = d.Mag();
+        sprintf(nam,"%sPart%dCable",name,i-1);
+        fTubes->AddAt( new TGeoCtub(nam,fRmin,fRmax,0.5*s,0.0,360.0,n0.x(),
+                                   n0.y(),n0.z(),n1.x(),n1.y(),n1.z()),i-1);
+        n0 = -n1;
+        x0 = x1;
+    } // end for i
+}
+//----------------------------------------------------------------------
+AliITSGeoCable::~AliITSGeoCable(){
+    //
+    // Inputs:
+    //    none.
+    // Outputs:
+    //    none.
+    // Return:
+    //    none.
+    Int_t i;
 
-//ClassImp(AliITSTube)
+    if(fTubes){
+        for(i=0;i<fTubes->GetEntries();i++) 
+            delete (TGeoCtub*)(fTubes->At(i));
+        delete fTubes;
+    } // end if
+    fTubes = 0;
+    if(fTranRot){
+        for(i=0;i<fTranRot->GetEntries();i++) 
+            delete (TGeoCombiTrans*)(fTranRot->At(i));
+        delete fTranRot;
+    } // end if
+    fTranRot  = 0;
+}
+//----------------------------------------------------------------------