Finished geometry objects defninistions, and fixed bugs in SSD cone geometry.
authornilsen <nilsen@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 22 Apr 2003 20:41:08 +0000 (20:41 +0000)
committernilsen <nilsen@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 22 Apr 2003 20:41:08 +0000 (20:41 +0000)
ITS/AliITSBaseGeometry.cxx
ITS/AliITSBaseGeometry.h
ITS/AliITSGeometrySSDCone.cxx
ITS/ITSLinkDef.h

index e01e52a..59ed015 100644 (file)
@@ -15,6 +15,9 @@
 
 /*
 $Log$
+Revision 1.4  2003/04/17 22:29:23  nilsen
+Geometry bug fixes and the like. Work still progressing.
+
 Revision 1.3  2003/03/27 08:49:39  hristov
 Initialization of static data members in the implementation file
 
@@ -47,6 +50,7 @@ Created February 2003.
 #include <TTUBE.h>
 #include <TTUBS.h>
 #include <TPCON.h>
+#include <TVector3.h>
 #include <TFile.h>    // only required for Tracking function?
 #include <TCanvas.h>
 #include <TObjArray.h>
@@ -211,7 +215,7 @@ Int_t AliITSBaseGeometry::ITSIndexToITSG3name(const Int_t i){
     return k;
 }
 //______________________________________________________________________
-Int_t AliITSBaseGeometry::ITSG3VnameToIndex(const char *name)const{
+Int_t AliITSBaseGeometry::ITSG3VnameToIndex(const char *name){
     // Given the last three characters of the ITS Geant3 volume name,
     // this returns the index. The valid characters must be in the range
     // '0' through 'Z'. This will include all upper case letter and the
@@ -277,7 +281,7 @@ Int_t AliITSBaseGeometry::GetVolumeIndex(const TString &a){
 //______________________________________________________________________
 void AliITSBaseGeometry::Box(const char *gnam,const TString &dis,
                             Double_t dx,Double_t dy,Double_t dz,Int_t med){
-    // Interface to TMC->Gsvolu() for ITS bos geometries. Box with faces
+    // Interface to TMC->Gsvolu() for ITS box geometries. Box with faces
     // perpendicular to the axes. It has 3 paramters. See SetScale() for
     // units. Default units are geant 3 [cm].
     // Inputs:
@@ -303,6 +307,32 @@ void AliITSBaseGeometry::Box(const char *gnam,const TString &dis,
     gMC->Gsvolu(name,"BOX ",GetMed(med),param,3);
 }
 //______________________________________________________________________
+void AliITSBaseGeometry::Box(AliITSBoxData &d,Int_t med){
+    // Interface to TMC->Gsvolu() for ITS box geometries. Box with faces
+    // perpendicular to the axes. It has 3 paramters. See SetScale() for
+    // units. Default units are geant 3 [cm].
+    // Inputs:
+    //    AliITSBoxData &d   Structure with the Box parameters defined.
+    //    Int_t         med  media index number.
+    // Output:
+    //    none.
+    // Return.
+    //    none.
+    char name[5];
+    Float_t param[3];
+    Int_t i,k;
+    char *j = (char *) &k;
+
+    param[0] = fScale*d.DxAt();
+    param[1] = fScale*d.DyAt();
+    param[2] = fScale*d.DzAt();
+    d.SetVid(AddVolName((d.GetName())->Data()));
+    k = ITSIndexToITSG3name(d.GetVid());
+    for(i=0;i<4;i++) name[i] = j[i];
+    name[4] = '\0';
+    gMC->Gsvolu(name,"BOX ",GetMed(med),param,3);
+}
+//______________________________________________________________________
 void AliITSBaseGeometry::Trapezoid1(const char *gnam,const TString &dis,
                                    Double_t dxn,Double_t dxp,Double_t dy,
                                    Double_t dz,Int_t med){
@@ -336,6 +366,33 @@ void AliITSBaseGeometry::Trapezoid1(const char *gnam,const TString &dis,
     gMC->Gsvolu(name,"TRD1",GetMed(med),param,4);
 }
 //______________________________________________________________________
+void AliITSBaseGeometry::Trapezoid1(AliITSTrapezoid1Data &d,Int_t med){
+    // Interface to TMC->Gsvolu() for ITS TRD1 geometries. Trapezoid with the 
+    // x dimension varing along z. It has 4 parameters. See SetScale() for
+    // units. Default units are geant 3 [cm].
+    // Inputs:
+    //    AliITSTrapezoid1Data &d   Structure with the Trapazoid data in it.
+    //    Int_t                med  media index number.
+    // Output:
+    //    none.
+    // Return.
+    //    none.
+    char name[5];
+    Float_t param[4];
+    Int_t i,k;
+    char *j = (char *) &k;
+
+    param[0] = fScale*d.DxAt(0);
+    param[1] = fScale*d.DxAt(1);
+    param[2] = fScale*d.DyAt();
+    param[3] = fScale*d.DzAt();
+    d.SetVid(AddVolName((d.GetName())->Data()));
+    k = ITSIndexToITSG3name(d.GetVid());
+    for(i=0;i<4;i++) name[i] = j[i];
+    name[4] = '\0';
+    gMC->Gsvolu(name,"TRD1",GetMed(med),param,4);
+}
+//______________________________________________________________________
 void AliITSBaseGeometry::Trapezoid2(const char *gnam,const TString &dis,
                                    Double_t dxn,Double_t dxp,Double_t dyn,
                                    Double_t dyp,Double_t dz,Int_t med){
@@ -373,6 +430,34 @@ void AliITSBaseGeometry::Trapezoid2(const char *gnam,const TString &dis,
     gMC->Gsvolu(name,"TRD2",GetMed(med),param,5);
 }
 //______________________________________________________________________
+void AliITSBaseGeometry::Trapezoid2(AliITSTrapezoid2Data &d,Int_t med){
+    // Interface to TMC->Gsvolu() for ITS TRD2 geometries. Trapezoid with the 
+    // x and y dimension varing along z. It has 5 parameters. See SetScale() 
+    // for units. Default units are geant 3 [cm].
+    // Inputs:
+    //    AliITSTrapezoid2Data &d   Structure with the Trapazoid data in it.
+    //    Int_t                med  media index number.
+    // Output:
+    //    none.
+    // Return.
+    //    none.
+    char name[5];
+    Float_t param[5];
+    Int_t i,k;
+    char *j = (char *) &k;
+
+    param[0] = fScale*d.DxAt(0);
+    param[1] = fScale*d.DxAt(1);
+    param[2] = fScale*d.DyAt(0);
+    param[3] = fScale*d.DyAt(1);
+    param[4] = fScale*d.DzAt();
+    d.SetVid(AddVolName((d.GetName())->Data()));
+    k = ITSIndexToITSG3name(d.GetVid());
+    for(i=0;i<4;i++) name[i] = j[i];
+    name[4] = '\0';
+    gMC->Gsvolu(name,"TRD2",GetMed(med),param,5);
+}
+//______________________________________________________________________
 void AliITSBaseGeometry::Trapezoid(const char *gnam,const TString &dis,
                                   Double_t dz,Double_t thet,Double_t phi,
                                   Double_t h1,Double_t bl1,Double_t tl1,
@@ -439,6 +524,252 @@ void AliITSBaseGeometry::Trapezoid(const char *gnam,const TString &dis,
     gMC->Gsvolu(name,"TRAP",GetMed(med),param,11);
 }
 //______________________________________________________________________
+void AliITSBaseGeometry::Trapezoid(AliITSTrapezoidData &d,Int_t med){
+    // Interface to TMC->Gsvolu() for ITS TRAP geometries. General Trapezoid, 
+    // The faces perpendicular to z are trapezia and their centers are not 
+    // necessarily on a line parallel to the z axis. This shape has 11 
+    // parameters, but only cosidering that the faces should be planar, only
+    // 9 are really independent. A check is performed on the user parameters 
+    // and a message is printed in case of non-planar faces. Ignoring this
+    // warning may cause unpredictable effects at tracking time. See 
+    // SetScale() for units. Default units are geant 3 [cm].
+    // Inputs:
+    //    AliITSTrapezoidData &d   Structure with the Trapazoid data in it.
+    //    Int_t    med        media index number.
+    // Output:
+    //    none.
+    // Return.
+    //    none.
+    char name[5];
+    Float_t param[11];
+    Int_t i,k;
+    char *j = (char *) &k;
+
+    param[0] = fScale*d.DzAt();
+    param[1] = d.Theta();
+    param[2] = d.Phi();
+    param[3] = fScale*d.HAt(0);
+    param[4] = fScale*d.Bl(0);
+    param[5] = fScale*d.Tl(0);
+    param[6] = d.Alpha(0);
+    param[7] = fScale*d.HAt(1);
+    param[8] = fScale*d.Bl(1);
+    param[9] = fScale*d.Tl(1);
+    param[10] = d.Alpha(1);
+    d.SetVid(AddVolName((d.GetName())->Data()));
+    k = ITSIndexToITSG3name(d.GetVid());
+    for(i=0;i<4;i++) name[i] = j[i];
+    name[4] = '\0';
+    gMC->Gsvolu(name,"TRAP",GetMed(med),param,11);
+}
+//______________________________________________________________________
+void AliITSBaseGeometry::TwistedTrapezoid(const char *gnam,
+                                         const TString &dis,
+                                Double_t dz,Double_t thet,Double_t phi,
+                                Double_t twist,Double_t h1,Double_t bl1,
+                                Double_t tl1,Double_t apl1,Double_t h2,
+                                Double_t bl2,Double_t tl2,Double_t apl2,
+                                Int_t med){
+    // Interface to TMC->Gsvolu() for ITS GTRA geometries. General twisted 
+    // trapazoid. The faces perpendicular to z are trapazia and their centers 
+    // are not necessarily on a line parallel to the z axis as the TRAP. 
+    // Additionally, the faces may be twisted so that none of their edges are 
+    // parallel. It is a TRAP shape, exept that it is twisted in the x-y 
+    // plane as a function of z. The parallel sides perpendicular to the x 
+    // axis are rotated with respect to the x axis by an angle TWIST, which 
+    // is one of the parameters. The shape is defined by the eight corners 
+    // and is assumed to be constructed of straight lines joingin points on 
+    // the boundry of the trapezoidal face at Z=-dz to the coresponding 
+    // points on the face at z=+dz. Divisions are not allowed. It has 12 
+    // parameters. See SetScale() for units. Default units are geant 3 [cm].
+    // Note: This shape suffers from the same limitations than the TRAP. The
+    // tracking routines assume that the faces are planar, but htis
+    // constraint is not easily expressed in terms of the 12 parameters. 
+    // Additionally, no check on th efaces is performed in this case. Users 
+    // should avoid to use this shape as much as possible, and if they have
+    // to do so, they should make sure that the faces are really planes. 
+    // If this is not the case, the result of the trasport is unpredictable. 
+    // To accelerat ethe computations necessary for trasport, 18 additioanl 
+    // parameters are calculated for this shape are 1 DXODZ dx/dz of the 
+    // line joing the centers of the faces at z=+_dz. 2 DYODZ dy/dz of the 
+    // line joing the centers of the faces at z=+_dz.
+    // 3 XO1    x at z=0 for line joing the + on parallel side, perpendicular 
+    //          corners at z=+_dz.
+    // 4 YO1    y at z=0 for line joing the + on parallel side, + on 
+    //          perpendicular corners at z=+-dz.
+    // 5 DXDZ1  dx/dz for line joing the + on parallel side, + on 
+    //          perpendicular corners at z=+-dz.
+    // 6 DYDZ1  dy/dz for line joing the + on parallel side, + on 
+    //          perpendicular corners at z=+-dz.
+    // 7 X02    x at z=0 for line joing the - on parallel side, + on 
+    //          perpendicular corners at z=+-dz.
+    // 8 YO2    y at z=0 for line joing the - on parallel side, + on 
+    //          perpendicular corners at z=+-dz.
+    // 9 DXDZ2  dx/dz for line joing the - on parallel side, + on 
+    //          perpendicular corners at z=+-dz.
+    // 10 DYDZ2dy/dz for line joing the - on parallel side, + on 
+    //          perpendicular corners at z=+-dz.
+    // 11 XO3   x at z=0 for line joing the - on parallel side, - on 
+    //          perpendicular corners at z=+-dz.
+    // 12 YO3   y at z=0 for line joing the - on parallel side, - on 
+    //          perpendicular corners at z=+-dz.
+    // 13 DXDZ3 dx/dzfor line joing the - on parallel side, - on 
+    //          perpendicular corners at z=+-dz.
+    // 14 DYDZ3 dydz for line joing the - on parallel side, - on 
+    //          perpendicular corners at z=+-dz.
+    // 15 XO4   x at z=0 for line joing the + on parallel side, - on 
+    //          perpendicular corners at z=+-dz.
+    // 16 YO4   y at z=0 for line joing the + on parallel side, - on 
+    //          perpendicular corners at z=+-dz.
+    // 17 DXDZ4 dx/dz for line joing the + on parallel side, - on 
+    //          perpendicular corners at z=+-dz.
+    // 18 DYDZ4 dydz for line joing the + on parallel side, - on 
+    //          perpendicular corners at z=+-dz.
+    // Inputs:
+    //    const char *gnam  3 character geant volume name. The letter "I"
+    //                        is appended to the front to indecate that this
+    //                        is an ITS volume.
+    //    TString &dis        String containging part discription.
+    //    Double_t dz         half-length along the z axis.
+    //    Double_t thet       polar angle of the line joing the center of the 
+    //                        face at -dz to the center of the one at +dz 
+    //                        [degrees].
+    //    Double_t phi        Azymuthal angle of teh line joing the centre of 
+    //                        the face at -dz to the center of the one at +dz 
+    //                        [degrees].
+    //    Double_t twist      Twist angle of the faces parallel to the x-y 
+    //                        plane at z=+-dz around an axis parallel to z 
+    //                        passing through their centre [degrees].
+    //    Double_t h1         Half-length along y of the face at -dz.
+    //    Double_t bl1        half-length along x of the side -h1 in y of the 
+    //                        face at -dz in z.
+    //    Double_t tl1        half-length along x of the side at +h1 in y of 
+    //                        the face at -dz in z.
+    //    Double_t apl1       Angle with respect to the y ais from the center 
+    //                        of the side at -h1 in y to the centere of the 
+    //                        side at +h1 in y of the face at -dz in z 
+    //                        [degrees].
+    //    Double_t h2         half-length along the face at +dz.
+    //    Double_t bl2        half-length along x of the side at -h2 in y of 
+    //                        the face at -dz in z.
+    //    Double_t tl2        half-length along x of the side at +h2 in y of 
+    //                        the face at +dz in z.
+    //    Double_t apl2       angle with respect to the y axis from the 
+    //                        center of the side at -h2 in y to the center 
+    //                        of the side at +h2 in y of the face at +dz in 
+    //                        z [degrees].
+    //    Int_t    med        media index number.
+    // Output:
+    //    none.
+    // Return.
+    //    none.
+    char name[5];
+    Float_t param[12];
+
+    param[0] = fScale*dz;
+    param[1] = thet;
+    param[2] = phi;
+    param[3] = twist;
+    param[4] = fScale*h1;
+    param[5] = fScale*bl1;
+    param[6] = fScale*tl1;
+    param[7] = apl1;
+    param[8] = fScale*h2;
+    param[9] = fScale*bl2;
+    param[10] = fScale*tl2;
+    param[11] = apl2;
+    G3name(gnam,name);
+    gMC->Gsvolu(name,"GTRA",GetMed(med),param,12);
+}
+//______________________________________________________________________
+void AliITSBaseGeometry::TwistedTrapezoid(AliITSTrapezoidTwistedData &d,
+                                         Int_t med){
+    // Interface to TMC->Gsvolu() for ITS GTRA geometries. General twisted 
+    // trapazoid. The faces perpendicular to z are trapazia and their centers 
+    // are not necessarily on a line parallel to the z axis as the TRAP. 
+    // Additionally, the faces may be twisted so that none of their edges are 
+    // parallel. It is a TRAP shape, exept that it is twisted in the x-y 
+    // plane as a function of z. The parallel sides perpendicular to the x 
+    // axis are rotated with respect to the x axis by an angle TWIST, which 
+    // is one of the parameters. The shape is defined by the eight corners 
+    // and is assumed to be constructed of straight lines joingin points on 
+    // the boundry of the trapezoidal face at Z=-dz to the coresponding 
+    // points on the face at z=+dz. Divisions are not allowed. It has 12 
+    // parameters. See SetScale() for units. Default units are geant 3 [cm].
+    // Note: This shape suffers from the same limitations than the TRAP. The
+    // tracking routines assume that the faces are planar, but htis
+    // constraint is not easily expressed in terms of the 12 parameters. 
+    // Additionally, no check on th efaces is performed in this case. Users 
+    // should avoid to use this shape as much as possible, and if they have
+    // to do so, they should make sure that the faces are really planes. 
+    // If this is not the case, the result of the trasport is unpredictable. 
+    // To accelerat ethe computations necessary for trasport, 18 additioanl 
+    // parameters are calculated for this shape are 1 DXODZ dx/dz of the 
+    // line joing the centers of the faces at z=+_dz. 2 DYODZ dy/dz of the 
+    // line joing the centers of the faces at z=+_dz.
+    // 3 XO1    x at z=0 for line joing the + on parallel side, perpendicular 
+    //          corners at z=+_dz.
+    // 4 YO1    y at z=0 for line joing the + on parallel side, + on 
+    //          perpendicular corners at z=+-dz.
+    // 5 DXDZ1  dx/dz for line joing the + on parallel side, + on 
+    //          perpendicular corners at z=+-dz.
+    // 6 DYDZ1  dy/dz for line joing the + on parallel side, + on 
+    //          perpendicular corners at z=+-dz.
+    // 7 X02    x at z=0 for line joing the - on parallel side, + on 
+    //          perpendicular corners at z=+-dz.
+    // 8 YO2    y at z=0 for line joing the - on parallel side, + on 
+    //          perpendicular corners at z=+-dz.
+    // 9 DXDZ2  dx/dz for line joing the - on parallel side, + on 
+    //          perpendicular corners at z=+-dz.
+    // 10 DYDZ2dy/dz for line joing the - on parallel side, + on 
+    //          perpendicular corners at z=+-dz.
+    // 11 XO3   x at z=0 for line joing the - on parallel side, - on 
+    //          perpendicular corners at z=+-dz.
+    // 12 YO3   y at z=0 for line joing the - on parallel side, - on 
+    //          perpendicular corners at z=+-dz.
+    // 13 DXDZ3 dx/dzfor line joing the - on parallel side, - on 
+    //          perpendicular corners at z=+-dz.
+    // 14 DYDZ3 dydz for line joing the - on parallel side, - on 
+    //          perpendicular corners at z=+-dz.
+    // 15 XO4   x at z=0 for line joing the + on parallel side, - on 
+    //          perpendicular corners at z=+-dz.
+    // 16 YO4   y at z=0 for line joing the + on parallel side, - on 
+    //          perpendicular corners at z=+-dz.
+    // 17 DXDZ4 dx/dz for line joing the + on parallel side, - on 
+    //          perpendicular corners at z=+-dz.
+    // 18 DYDZ4 dydz for line joing the + on parallel side, - on 
+    //          perpendicular corners at z=+-dz.
+    // Inputs:
+    //    AliITSTrapezoidTwistedData &d   Structure with the tube parameters
+    //    Int_t                      med  media index number.
+    // Output:
+    //    none.
+    // Return.
+    //    none.
+    char name[5];
+    Float_t param[12];
+    Int_t i,k;
+    char *j = (char *) &k;
+
+    param[0] = fScale*d.DzAt();
+    param[1] = d.Theta();
+    param[2] = d.Phi();
+    param[3] = d.Twist();
+    param[4] = fScale*d.HAt(0);
+    param[5] = fScale*d.Bl(0);
+    param[6] = fScale*d.Tl(0);
+    param[7] = d.Alpha(0);
+    param[8] = fScale*d.HAt(1);
+    param[9] = fScale*d.Bl(1);
+    param[10] = fScale*d.Tl(1);
+    param[11] = d.Alpha(1);
+    d.SetVid(AddVolName((d.GetName())->Data()));
+    k = ITSIndexToITSG3name(d.GetVid());
+    for(i=0;i<4;i++) name[i] = j[i];
+    name[4] = '\0';
+    gMC->Gsvolu(name,"GTRA",GetMed(med),param,12);
+}
+//______________________________________________________________________
 void AliITSBaseGeometry::Tube(const char *gnam,const TString &dis,
                              Double_t rmin,Double_t rmax,Double_t dz,
                              Int_t med){
@@ -473,8 +804,8 @@ void AliITSBaseGeometry::Tube(AliITSTubeData &d,Int_t med){
     // 3 parameters. See SetScale() 
     // for units. Default units are geant 3 [cm].
     // Inputs:
-    //    AliITSTubeData     Structure with the tube parameters
-    //    Int_t    med        media index number.
+    //    AliITSTubeData &d    Structure with the tube parameters
+    //    Int_t          med   media index number.
     // Output:
     //    none.
     // Return.
@@ -528,6 +859,249 @@ void AliITSBaseGeometry::TubeSegment(const char *gnam,const TString &dis,
     gMC->Gsvolu(name,"TUBS",GetMed(med),param,5);
 }
 //______________________________________________________________________
+void AliITSBaseGeometry::TubeSegment(AliITSTubeSegData &d,Int_t med){
+    // Interface to TMC->Gsvolu() for ITS TUBE geometries. Phi segment of a 
+    // tube. It has 5  parameters. Phi1 should be smaller than phi2. If this
+    // is not the case, the system adds 360 degrees to phi2. See SetScale() 
+    // for units. Default units are geant 3 [cm].
+    // Inputs:
+    //    AliITSTubeSegData &d   Structure with the tube parameters
+    //    Int_t             med  media index number.
+    // Output:
+    //    none.
+    // Return.
+    //    none.
+    char name[5];
+    Float_t param[5];
+    Int_t i,k;
+    char *j = (char *) &k;
+
+    param[0] = fScale*d.Rmin();
+    param[1] = fScale*d.Rmax();
+    param[2] = fScale*d.DzAt();
+    param[3] = d.Phi0();
+    param[4] = d.Phi1();
+    d.SetVid(AddVolName((d.GetName())->Data()));
+    k = ITSIndexToITSG3name(d.GetVid());
+    for(i=0;i<4;i++) name[i] = j[i];
+    name[4] = '\0';
+    gMC->Gsvolu(name,"TUBS",GetMed(med),param,5);
+}
+//______________________________________________________________________
+void AliITSBaseGeometry::CutTube(const char *gnam,const TString &dis,
+                                Double_t rmin,Double_t rmax,Double_t dz,
+                                Double_t phi1,Double_t phi2,Double_t lx,
+                                Double_t ly,Double_t lz,Double_t hx,
+                                Double_t hy,Double_t hz,Int_t med){
+    // Interface to TMC->Gsvolu() for ITS CTUB geometries. Cut tube. A tube 
+    // cut at the extremities with planes not necessarily perpendicular to 
+    // the z axis. It has 11 parameters. See SetScale() for units. Default 
+    // units are geant 3 [cm]. phi1 should be smaller than phi2. If this is 
+    // not the case, the system adds 360 degrees to phi2.
+    // Inputs:
+    //    const char *gnam  3 character geant volume name. The letter "I"
+    //                        is appended to the front to indecate that this
+    //                        is an ITS volume.
+    //    TString &dis        String containging part discription.
+    //    Double_t rmin       Inner radius at z=0 where tube is narrowest.
+    //    Double_t rmax       Outer radius at z=0 where tube is narrowest.
+    //    Double_t dz         half-length along the z-axis
+    //    Double_t phi1       Starting angle of the segment [degree].
+    //    Double_t phi2       Ending angle of the segment [degree].
+    //    Double_t lx         x component of a unit vector perpendicular to 
+    //                        the face at -dz.
+    //    Double_t ly         y component of a unit vector perpendicular to 
+    //                        the face at -dz.
+    //    Double_t lz         z component of a unit vector perpendicular to 
+    //                        the face at -dz.
+    //    Double_t hx         x component of a unit vector perpendicular to 
+    //                        the face at +dz.
+    //    Double_t hy         y component of a unit vector perpendicular to 
+    //                        the face at +dz.
+    //    Double_t hz         z component of a unit vector perpendicular to 
+    //                        the face at +dz.
+    //    Int_t    med        media index number.
+    // Output:
+    //    none.
+    // Return.
+    //    none.
+    char name[5];
+    Float_t param[11];
+
+    param[0] = fScale*rmin;
+    param[1] = fScale*rmax;
+    param[2] = fScale*dz;
+    param[3] = phi1;
+    param[4] = phi2;
+    param[5] = lx;
+    param[6] = ly;
+    param[7] = lz;
+    param[8] = hx;
+    param[9] = hy;
+    param[10] = hz;
+    G3name(gnam,name);
+    gMC->Gsvolu(name,"CTUB",GetMed(med),param,11);
+}
+//______________________________________________________________________
+void AliITSBaseGeometry::CutTube(AliITSTubeCutData &d,Int_t med){
+    // Interface to TMC->Gsvolu() for ITS CTUB geometries. Cut tube. A tube 
+    // cut at the extremities with planes not necessarily perpendicular to 
+    // the z axis. It has 11 parameters. See SetScale() for units. Default 
+    // units are geant 3 [cm]. phi1 should be smaller than phi2. If this is 
+    // not the case, the system adds 360 degrees to phi2.
+    // Inputs:
+    //    AliITSTubeCutData &d    Structure with the tube parameters
+    //    Int_t             med    media index number.
+    // Output:
+    //    none.
+    // Return.
+    //    none.
+    char name[5];
+    Float_t param[11];
+    Int_t i,k;
+    char *j = (char *) &k;
+
+    param[0] = fScale*d.Rmin();
+    param[1] = fScale*d.Rmax();
+    param[2] = fScale*d.DzAt();
+    param[3] = d.Phi0();
+    param[4] = d.Phi1();
+    param[5] = d.Normal(0,0);
+    param[6] = d.Normal(0,1);
+    param[7] = d.Normal(0,2);
+    param[8] = d.Normal(1,0);
+    param[9] = d.Normal(1,1);
+    param[10] = d.Normal(1,2);
+    d.SetVid(AddVolName((d.GetName())->Data()));
+    k = ITSIndexToITSG3name(d.GetVid());
+    for(i=0;i<4;i++) name[i] = j[i];
+    name[4] = '\0';
+    gMC->Gsvolu(name,"CTUB",GetMed(med),param,11);
+}
+//______________________________________________________________________
+void AliITSBaseGeometry::TubeElliptical(const char *gnam,const TString &dis,
+                              Double_t p1,Double_t p2,Double_t dz,Int_t med){
+    // Interface to TMC->Gsvolu() for ITS ELTU geometries. Elliptical 
+    // cross-section Tube. It has 3 parameters. See SetScale() 
+    // for units. Default units are geant 3 [cm]. The equation of the surface 
+    // is x^2 * p1^-2 + y^2 * p2^-2 = 1.
+    // Inputs:
+    //    const char *gnam  3 character geant volume name. The letter "I"
+    //                        is appended to the front to indecate that this
+    //                        is an ITS volume.
+    //    TString &dis        String containging part discription.
+    //    Double_t p1         semi-axis of the elipse along x.
+    //    Double_t p2         semi-axis of the elipse along y.
+    //    Double_t dz         half-length along the z-axis
+    //    Int_t    med        media index number.
+    // Output:
+    //    none.
+    // Return.
+    //    none.
+    char name[5];
+    Float_t param[3];
+
+    param[0] = fScale*p1;
+    param[1] = fScale*p2;
+    param[2] = fScale*dz;
+    G3name(gnam,name);
+    gMC->Gsvolu(name,"ELTU",GetMed(med),param,3);
+}
+//______________________________________________________________________
+void AliITSBaseGeometry::TubeElliptical(AliITSTubeEllipticalData &d,
+                                       Int_t med){
+    // Interface to TMC->Gsvolu() for ITS ELTU geometries. Elliptical 
+    // cross-section Tube. It has 3 parameters. See SetScale() 
+    // for units. Default units are geant 3 [cm]. The equation of the surface 
+    // is x^2 * p1^-2 + y^2 * p2^-2 = 1.
+    // Inputs:
+    //    AliITSTubeElipticData &d  Structure with the tube parameters
+    //    Int_t                med  media index number.
+    // Output:
+    //    none.
+    // Return.
+    //    none.
+    char name[5];
+    Float_t param[3];
+    Int_t i,k;
+    char *j = (char *) &k;
+
+    param[0] = fScale*d.P0();
+    param[1] = fScale*d.P1();
+    param[2] = fScale*d.DzAt();
+    d.SetVid(AddVolName((d.GetName())->Data()));
+    k = ITSIndexToITSG3name(d.GetVid());
+    for(i=0;i<4;i++) name[i] = j[i];
+    name[4] = '\0';
+    gMC->Gsvolu(name,"ELTU",GetMed(med),param,3);
+}
+//______________________________________________________________________
+void AliITSBaseGeometry::HyperbolicTube(const char *gnam,const TString &dis,
+                              Double_t rmin,Double_t rmax,Double_t dz,
+                              Double_t thet,Int_t med){
+    // Interface to TMC->Gsvolu() for ITS HYPE geometries. Hyperbolic tube. 
+    // Fore example the inner and outer surfaces are hyperboloids, as would 
+    // be foumed by a system of cylinderical wires which were then rotated 
+    // tangentially about their centers. It has 4 parameters. See SetScale() 
+    // for units. Default units are geant 3 [cm]. The hyperbolic surfaces are 
+    // given by r^2 = (ztan(thet)^2 + r(z=0)^2.
+    // Inputs:
+    //    const char *gnam  3 character geant volume name. The letter "I"
+    //                        is appended to the front to indecate that this
+    //                        is an ITS volume.
+    //    TString &dis        String containging part discription.
+    //    Double_t rmin       Inner radius at z=0 where tube is narrowest.
+    //    Double_t rmax       Outer radius at z=0 where tube is narrowest.
+    //    Double_t dz         half-length along the z-axis
+    //    Double_t thet       stero angel of rotation of the two faces 
+    //                       [degrees].
+    //    Int_t    med        media index number.
+    // Output:
+    //    none.
+    // Return.
+    //    none.
+    char name[5];
+    Float_t param[4];
+
+    param[0] = fScale*rmin;
+    param[1] = fScale*rmax;
+    param[2] = fScale*dz;
+    param[3] = thet;
+    G3name(gnam,name);
+    gMC->Gsvolu(name,"HYPE",GetMed(med),param,4);
+}
+//______________________________________________________________________
+void AliITSBaseGeometry::HyperbolicTube(AliITSTubeHyperbolicData &d,
+                                       Int_t med){
+    // Interface to TMC->Gsvolu() for ITS HYPE geometries. Hyperbolic tube. 
+    // Fore example the inner and outer surfaces are hyperboloids, as would 
+    // be foumed by a system of cylinderical wires which were then rotated 
+    // tangentially about their centers. It has 4 parameters. See SetScale() 
+    // for units. Default units are geant 3 [cm]. The hyperbolic surfaces are 
+    // given by r^2 = (ztan(thet)^2 + r(z=0)^2.
+    // Inputs:
+    //    AliITSTubeHyperbolicData &d  Structure with the tube parameters
+    //    Int_t                    med  media index number.
+    // Output:
+    //    none.
+    // Return.
+    //    none.
+    char name[5];
+    Float_t param[4];
+    Int_t i,k;
+    char *j = (char *) &k;
+
+    param[0] = fScale*d.Rmin();
+    param[1] = fScale*d.Rmax();
+    param[2] = fScale*d.DzAt();
+    param[3] = d.Theta();
+    d.SetVid(AddVolName((d.GetName())->Data()));
+    k = ITSIndexToITSG3name(d.GetVid());
+    for(i=0;i<4;i++) name[i] = j[i];
+    name[4] = '\0';
+    gMC->Gsvolu(name,"HYPE",GetMed(med),param,4);
+}
+//______________________________________________________________________
 void AliITSBaseGeometry::Cone(const char *gnam,const TString &dis,
                              Double_t dz,Double_t rmin1,Double_t rmax1,
                              Double_t rmin2,Double_t rmax2,Int_t med){
@@ -561,6 +1135,34 @@ void AliITSBaseGeometry::Cone(const char *gnam,const TString &dis,
     gMC->Gsvolu(name,"CONS",GetMed(med),param,5);
 }
 //______________________________________________________________________
+void AliITSBaseGeometry::Cone(AliITSConeData &d,Int_t med){
+    // Interface to TMC->Gsvolu() for ITS Cone geometries. Conical tube. It 
+    // has 5 parameters. See SetScale() 
+    // for units. Default units are geant 3 [cm].
+    // Inputs:
+    //    AliITSConeData &d  Structure with the tube parameters
+    //    Int_t         med  media index number.
+    // Output:
+    //    none.
+    // Return.
+    //    none.
+    char name[5];
+    Float_t param[5];
+    Int_t i,k;
+    char *j = (char *) &k;
+
+    param[0] = fScale*d.DzAt();
+    param[1] = fScale*d.Rmin0();
+    param[2] = fScale*d.Rmax0();
+    param[3] = fScale*d.Rmin1();
+    param[4] = fScale*d.Rmax1();
+    d.SetVid(AddVolName((d.GetName())->Data()));
+    k = ITSIndexToITSG3name(d.GetVid());
+    for(i=0;i<4;i++) name[i] = j[i];
+    name[4] = '\0';
+    gMC->Gsvolu(name,"CONS",GetMed(med),param,5);
+}
+//______________________________________________________________________
 void AliITSBaseGeometry::ConeSegment(const char *gnam,const TString &dis,
                                     Double_t dz,Double_t rmin1,
                                     Double_t rmax1,Double_t rmin2,
@@ -601,6 +1203,122 @@ void AliITSBaseGeometry::ConeSegment(const char *gnam,const TString &dis,
     gMC->Gsvolu(name,"CONS",GetMed(med),param,7);
 }
 //______________________________________________________________________
+void AliITSBaseGeometry::ConeSegment(AliITSConeSegData &d,Int_t med){
+    // Interface to TMC->Gsvolu() for ITS ConS geometries. One segment of a 
+    // conical tube. It has 7 parameters. Phi1 should be smaller than phi2. 
+    // If this is not the case, the system adds 360 degrees to phi2. See 
+    // SetScale() for units. Default units are geant 3 [cm].
+    // Inputs:
+    //    AliITSConeSegData &d   Structure with the tube parameters
+    //    Int_t             med  media index number.
+    // Output:
+    //    none.
+    // Return.
+    //    none.
+    char name[5];
+    Float_t param[7];
+    Int_t i,k;
+    char *j = (char *) &k;
+
+    param[0] = fScale*d.DzAt();
+    param[1] = fScale*d.Rmin0();
+    param[2] = fScale*d.Rmax0();
+    param[3] = fScale*d.Rmin1();
+    param[4] = fScale*d.Rmax1();
+    param[5] = d.Phi0();
+    param[6] = d.Phi1();
+    d.SetVid(AddVolName((d.GetName())->Data()));
+    k = ITSIndexToITSG3name(d.GetVid());
+    for(i=0;i<4;i++) name[i] = j[i];
+    name[4] = '\0';
+    gMC->Gsvolu(name,"CONS",GetMed(med),param,7);
+}
+//______________________________________________________________________
+void AliITSBaseGeometry::PolyCone(const char *gnam,const TString &dis,
+                                 Double_t phi1,Double_t dphi,Int_t nz,
+                                 Double_t *z,Double_t *rmin,Double_t *rmax,
+                                 Int_t med){
+    // Interface to TMC->Gsvolu() for ITS PCON geometry. Poly-cone It has 9 
+    // parameters or more. See SetScale() for units. Default units are geant
+    // 3 [cm].
+    // Inputs:
+    //    const char *gnam  3 character geant volume name. The letter "I"
+    //                        is appended to the front to indecate that this
+    //                        is an ITS volume.
+    //    TString &dis        String containging part discription.
+    //    Double_t phi1       the azimuthal angle at which the volume begins 
+    //                        (angles are counted clouterclockwise) [degrees].
+    //    Double_t dphi       opening angle of the volume, which extends from 
+    //                        phi1 to phi1+dphi [degree].
+    //    Int_t nz            number of planes perpendicular to the z axis 
+    //                        where the dimension of the section is given - 
+    //                        this number should be at least 2 and NP triples 
+    //                        of number must follow.
+    //    Double_t *z         Array [nz] of z coordinate of the section.
+    //    Double_t *rmin      Array [nz] of radius of teh inner circle in the 
+    //                        cross-section.
+    //    Double_t *rmax      Array [nz] of radius of the outer circle in the 
+    //                        cross-section.
+    //    Int_t    med        media index number.
+    // Output:
+    //    none.
+    // Return.
+    //    none.
+    char name[5];
+    Float_t *param;
+    Int_t n,i;
+
+    n = 3+3*nz;
+    param = new Float_t[n];
+    param[0] = phi1;
+    param[1] = dphi;
+    param[2] = (Float_t) nz;
+    for(i=0;i<nz;i++){
+       param[3+3*i] = fScale*z[i];
+       param[4+3*i] = fScale*rmin[i];
+       param[5+3*i] = fScale*rmax[i];
+    } // end for i
+    G3name(gnam,name);
+    gMC->Gsvolu(name,"PCON",GetMed(med),param,n);
+
+    delete[] param;
+}
+//______________________________________________________________________
+void AliITSBaseGeometry::PolyCone(AliITSPConeData &d,Int_t med){
+    // Interface to TMC->Gsvolu() for ITS PCON geometry. Poly-cone It has 9 
+    // parameters or more. See SetScale() for units. Default units are geant
+    // 3 [cm].
+    // Inputs:
+    //    AliITSPConeData &d  Object with poly cone data stored in it.
+    //    Int_t    med        media index number.
+    // Output:
+    //    none.
+    // Return.
+    //    none.
+    char name[5];
+    Float_t *param;
+    Int_t n,i,k;
+    char *j = (char *) &k;
+
+    n = 3+3*d.Nz();
+    param = new Float_t[n];
+    param[0] = d.Phi0();
+    param[1] = d.DPhi();
+    param[2] = (Float_t) d.Nz();
+    for(i=0;i<d.Nz();i++){
+       param[3+3*i] = fScale*d.ZAt(i);
+       param[4+3*i] = fScale*d.Rmin(i);
+       param[5+3*i] = fScale*d.Rmax(i);
+    } // end for if
+    d.SetVid(AddVolName((d.GetName())->Data()));
+    k = ITSIndexToITSG3name(d.GetVid());
+    for(i=0;i<4;i++) name[i] = j[i];
+    name[4] = '\0';
+    gMC->Gsvolu(name,"PCON",GetMed(med),param,n);
+
+    delete[] param;
+}
+//______________________________________________________________________
 void AliITSBaseGeometry::Sphere(const char *gnam,const TString &dis,
                                Double_t rmin,Double_t rmax,Double_t the1,
                                Double_t the2,Double_t phi1,Double_t phi2,
@@ -637,6 +1355,35 @@ void AliITSBaseGeometry::Sphere(const char *gnam,const TString &dis,
     gMC->Gsvolu(name,"SPHE",GetMed(med),param,6);
 }
 //______________________________________________________________________
+void AliITSBaseGeometry::Sphere(AliITSSphereData &d,Int_t med){
+    // Interface to TMC->Gsvolu() for ITS SPHE geometries. Segment of a 
+    // sphereical shell. It has 6 parameters. See SetScale() 
+    // for units. Default units are geant 3 [cm].
+    // Inputs:
+    //    AliITSSphereData &d   Structure with the tube parameters
+    //    Int_t            med  media index number.
+    // Output:
+    //    none.
+    // Return.
+    //    none.
+    char name[5];
+    Float_t param[6];
+    Int_t i,k;
+    char *j = (char *) &k;
+
+    param[0] = fScale*d.Rmin();
+    param[1] = fScale*d.Rmax();
+    param[2] = d.Theta0();
+    param[3] = d.Theta1();
+    param[4] = d.Phi0();
+    param[5] = d.Phi1();
+    d.SetVid(AddVolName((d.GetName())->Data()));
+    k = ITSIndexToITSG3name(d.GetVid());
+    for(i=0;i<4;i++) name[i] = j[i];
+    name[4] = '\0';
+    gMC->Gsvolu(name,"SPHE",GetMed(med),param,6);
+}
+//______________________________________________________________________
 void AliITSBaseGeometry::Parallelepiped(const char *gnam,const TString &dis,
                                        Double_t dx,Double_t dy,Double_t dz,
                                        Double_t alpha,Double_t thet,
@@ -677,6 +1424,35 @@ void AliITSBaseGeometry::Parallelepiped(const char *gnam,const TString &dis,
     gMC->Gsvolu(name,"PARA",GetMed(med),param,6);
 }
 //______________________________________________________________________
+void AliITSBaseGeometry::Parallelepiped(AliITSParallelpipedData &d,Int_t med){
+    // Interface to TMC->Gsvolu() for ITS PARA geometries. Parallelepiped. It 
+    // has 6 parameters. See SetScale() for units. Default units are geant 3 
+    // [cm].
+    // Inputs:
+    //    AliITSParrellepipedData &d  Structre witht the volume data in it.
+    //    Int_t                   med  media index number.
+    // Output:
+    //    none.
+    // Return.
+    //    none.
+    char name[5];
+    Float_t param[6];
+    Int_t i,k;
+    char *j = (char *) &k;
+
+    param[0] = fScale*d.DxAt();
+    param[1] = fScale*d.DyAt();
+    param[2] = fScale*d.DzAt();
+    param[3] = d.Alpha();
+    param[4] = d.Theta();
+    param[5] = d.Phi();
+    d.SetVid(AddVolName((d.GetName())->Data()));
+    k = ITSIndexToITSG3name(d.GetVid());
+    for(i=0;i<4;i++) name[i] = j[i];
+    name[4] = '\0';
+    gMC->Gsvolu(name,"PARA",GetMed(med),param,6);
+}
+//______________________________________________________________________
 void AliITSBaseGeometry::PolyGon(const char *gnam,const TString &dis,
                                 Double_t phi1,Double_t dphi,Int_t npdv,
                                 Int_t nz,Double_t *z,Double_t *rmin,
@@ -768,330 +1544,6 @@ void AliITSBaseGeometry::PolyGon(AliITSPGonData &d,Int_t med){
     delete[] param;
 }
 //______________________________________________________________________
-void AliITSBaseGeometry::PolyCone(const char *gnam,const TString &dis,
-                                 Double_t phi1,Double_t dphi,Int_t nz,
-                                 Double_t *z,Double_t *rmin,Double_t *rmax,
-                                 Int_t med){
-    // Interface to TMC->Gsvolu() for ITS PCON geometry. Poly-cone It has 9 
-    // parameters or more. See SetScale() for units. Default units are geant
-    // 3 [cm].
-    // Inputs:
-    //    const char *gnam  3 character geant volume name. The letter "I"
-    //                        is appended to the front to indecate that this
-    //                        is an ITS volume.
-    //    TString &dis        String containging part discription.
-    //    Double_t phi1       the azimuthal angle at which the volume begins 
-    //                        (angles are counted clouterclockwise) [degrees].
-    //    Double_t dphi       opening angle of the volume, which extends from 
-    //                        phi1 to phi1+dphi [degree].
-    //    Int_t nz            number of planes perpendicular to the z axis 
-    //                        where the dimension of the section is given - 
-    //                        this number should be at least 2 and NP triples 
-    //                        of number must follow.
-    //    Double_t *z         Array [nz] of z coordinate of the section.
-    //    Double_t *rmin      Array [nz] of radius of teh inner circle in the 
-    //                        cross-section.
-    //    Double_t *rmax      Array [nz] of radius of the outer circle in the 
-    //                        cross-section.
-    //    Int_t    med        media index number.
-    // Output:
-    //    none.
-    // Return.
-    //    none.
-    char name[5];
-    Float_t *param;
-    Int_t n,i;
-
-    n = 3+3*nz;
-    param = new Float_t[n];
-    param[0] = phi1;
-    param[1] = dphi;
-    param[2] = (Float_t) nz;
-    for(i=0;i<nz;i++){
-       param[3+3*i] = fScale*z[i];
-       param[4+3*i] = fScale*rmin[i];
-       param[5+3*i] = fScale*rmax[i];
-    } // end for i
-    G3name(gnam,name);
-    gMC->Gsvolu(name,"PCON",GetMed(med),param,n);
-
-    delete[] param;
-}
-//______________________________________________________________________
-void AliITSBaseGeometry::PolyCone(AliITSPConeData &d,Int_t med){
-    // Interface to TMC->Gsvolu() for ITS PCON geometry. Poly-cone It has 9 
-    // parameters or more. See SetScale() for units. Default units are geant
-    // 3 [cm].
-    // Inputs:
-    //    AliITSPConeData &d  Object with poly cone data stored in it.
-    //    Int_t    med        media index number.
-    // Output:
-    //    none.
-    // Return.
-    //    none.
-    char name[5];
-    Float_t *param;
-    Int_t n,i,k;
-    char *j = (char *) &k;
-
-    n = 3+3*d.Nz();
-    param = new Float_t[n];
-    param[0] = d.Phi0();
-    param[1] = d.DPhi();
-    param[2] = (Float_t) d.Nz();
-    for(i=0;i<d.Nz();i++){
-       param[3+3*i] = fScale*d.ZAt(i);
-       param[4+3*i] = fScale*d.Rmin(i);
-       param[5+3*i] = fScale*d.Rmax(i);
-    } // end for if
-    d.SetVid(AddVolName((d.GetName())->Data()));
-    k = ITSIndexToITSG3name(d.GetVid());
-    for(i=0;i<4;i++) name[i] = j[i];
-    name[4] = '\0';
-    gMC->Gsvolu(name,"PCON",GetMed(med),param,n);
-
-    delete[] param;
-}
-//______________________________________________________________________
-void AliITSBaseGeometry::TubeElliptical(const char *gnam,const TString &dis,
-                              Double_t p1,Double_t p2,Double_t dz,Int_t med){
-    // Interface to TMC->Gsvolu() for ITS ELTU geometries. Elliptical 
-    // cross-section Tube. It has 3 parameters. See SetScale() 
-    // for units. Default units are geant 3 [cm]. The equation of the surface 
-    // is x^2 * p1^-2 + y^2 * p2^-2 = 1.
-    // Inputs:
-    //    const char *gnam  3 character geant volume name. The letter "I"
-    //                        is appended to the front to indecate that this
-    //                        is an ITS volume.
-    //    TString &dis        String containging part discription.
-    //    Double_t p1         semi-axis of the elipse along x.
-    //    Double_t p2         semi-axis of the elipse along y.
-    //    Double_t dz         half-length along the z-axis
-    //    Int_t    med        media index number.
-    // Output:
-    //    none.
-    // Return.
-    //    none.
-    char name[5];
-    Float_t param[3];
-
-    param[0] = fScale*p1;
-    param[1] = fScale*p2;
-    param[2] = fScale*dz;
-    G3name(gnam,name);
-    gMC->Gsvolu(name,"ELTU",GetMed(med),param,3);
-}
-//______________________________________________________________________
-void AliITSBaseGeometry::HyperbolicTube(const char *gnam,const TString &dis,
-                              Double_t rmin,Double_t rmax,Double_t dz,
-                              Double_t thet,Int_t med){
-    // Interface to TMC->Gsvolu() for ITS HYPE geometries. Hyperbolic tube. 
-    // Fore example the inner and outer surfaces are hyperboloids, as would 
-    // be foumed by a system of cylinderical wires which were then rotated 
-    // tangentially about their centers. It has 4 parameters. See SetScale() 
-    // for units. Default units are geant 3 [cm]. The hyperbolic surfaces are 
-    // given by r^2 = (ztan(thet)^2 + r(z=0)^2.
-    // Inputs:
-    //    const char *gnam  3 character geant volume name. The letter "I"
-    //                        is appended to the front to indecate that this
-    //                        is an ITS volume.
-    //    TString &dis        String containging part discription.
-    //    Double_t rmin       Inner radius at z=0 where tube is narrowest.
-    //    Double_t rmax       Outer radius at z=0 where tube is narrowest.
-    //    Double_t dz         half-length along the z-axis
-    //    Double_t thet       stero angel of rotation of the two faces 
-    //                       [degrees].
-    //    Int_t    med        media index number.
-    // Output:
-    //    none.
-    // Return.
-    //    none.
-    char name[5];
-    Float_t param[4];
-
-    param[0] = fScale*rmin;
-    param[1] = fScale*rmax;
-    param[2] = fScale*dz;
-    param[3] = thet;
-    G3name(gnam,name);
-    gMC->Gsvolu(name,"HYPE",GetMed(med),param,4);
-}
-//______________________________________________________________________
-void AliITSBaseGeometry::TwistedTrapezoid(const char *gnam,
-                                         const TString &dis,
-                                Double_t dz,Double_t thet,Double_t phi,
-                                Double_t twist,Double_t h1,Double_t bl1,
-                                Double_t tl1,Double_t apl1,Double_t h2,
-                                Double_t bl2,Double_t tl2,Double_t apl2,
-                                Int_t med){
-    // Interface to TMC->Gsvolu() for ITS GTRA geometries. General twisted 
-    // trapazoid. The faces perpendicular to z are trapazia and their centers 
-    // are not necessarily on a line parallel to the z axis as the TRAP. 
-    // Additionally, the faces may be twisted so that none of their edges are 
-    // parallel. It is a TRAP shape, exept that it is twisted in the x-y 
-    // plane as a function of z. The parallel sides perpendicular to the x 
-    // axis are rotated with respect to the x axis by an angle TWIST, which 
-    // is one of the parameters. The shape is defined by the eight corners 
-    // and is assumed to be constructed of straight lines joingin points on 
-    // the boundry of the trapezoidal face at Z=-dz to the coresponding 
-    // points on the face at z=+dz. Divisions are not allowed. It has 12 
-    // parameters. See SetScale() for units. Default units are geant 3 [cm].
-    // Note: This shape suffers from the same limitations than the TRAP. The
-    // tracking routines assume that the faces are planar, but htis
-    // constraint is not easily expressed in terms of the 12 parameters. 
-    // Additionally, no check on th efaces is performed in this case. Users 
-    // should avoid to use this shape as much as possible, and if they have
-    // to do so, they should make sure that the faces are really planes. 
-    // If this is not the case, the result of the trasport is unpredictable. 
-    // To accelerat ethe computations necessary for trasport, 18 additioanl 
-    // parameters are calculated for this shape are 1 DXODZ dx/dz of the 
-    // line joing the centers of the faces at z=+_dz. 2 DYODZ dy/dz of the 
-    // line joing the centers of the faces at z=+_dz.
-    // 3 XO1    x at z=0 for line joing the + on parallel side, perpendicular 
-    //          corners at z=+_dz.
-    // 4 YO1    y at z=0 for line joing the + on parallel side, + on 
-    //          perpendicular corners at z=+-dz.
-    // 5 DXDZ1  dx/dz for line joing the + on parallel side, + on 
-    //          perpendicular corners at z=+-dz.
-    // 6 DYDZ1  dy/dz for line joing the + on parallel side, + on 
-    //          perpendicular corners at z=+-dz.
-    // 7 X02    x at z=0 for line joing the - on parallel side, + on 
-    //          perpendicular corners at z=+-dz.
-    // 8 YO2    y at z=0 for line joing the - on parallel side, + on 
-    //          perpendicular corners at z=+-dz.
-    // 9 DXDZ2  dx/dz for line joing the - on parallel side, + on 
-    //          perpendicular corners at z=+-dz.
-    // 10 DYDZ2dy/dz for line joing the - on parallel side, + on 
-    //          perpendicular corners at z=+-dz.
-    // 11 XO3   x at z=0 for line joing the - on parallel side, - on 
-    //          perpendicular corners at z=+-dz.
-    // 12 YO3   y at z=0 for line joing the - on parallel side, - on 
-    //          perpendicular corners at z=+-dz.
-    // 13 DXDZ3 dx/dzfor line joing the - on parallel side, - on 
-    //          perpendicular corners at z=+-dz.
-    // 14 DYDZ3 dydz for line joing the - on parallel side, - on 
-    //          perpendicular corners at z=+-dz.
-    // 15 XO4   x at z=0 for line joing the + on parallel side, - on 
-    //          perpendicular corners at z=+-dz.
-    // 16 YO4   y at z=0 for line joing the + on parallel side, - on 
-    //          perpendicular corners at z=+-dz.
-    // 17 DXDZ4 dx/dz for line joing the + on parallel side, - on 
-    //          perpendicular corners at z=+-dz.
-    // 18 DYDZ4 dydz for line joing the + on parallel side, - on 
-    //          perpendicular corners at z=+-dz.
-    // Inputs:
-    //    const char *gnam  3 character geant volume name. The letter "I"
-    //                        is appended to the front to indecate that this
-    //                        is an ITS volume.
-    //    TString &dis        String containging part discription.
-    //    Double_t dz         half-length along the z axis.
-    //    Double_t thet       polar angle of the line joing the center of the 
-    //                        face at -dz to the center of the one at +dz 
-    //                        [degrees].
-    //    Double_t phi        Azymuthal angle of teh line joing the centre of 
-    //                        the face at -dz to the center of the one at +dz 
-    //                        [degrees].
-    //    Double_t twist      Twist angle of the faces parallel to the x-y 
-    //                        plane at z=+-dz around an axis parallel to z 
-    //                        passing through their centre [degrees].
-    //    Double_t h1         Half-length along y of the face at -dz.
-    //    Double_t bl1        half-length along x of the side -h1 in y of the 
-    //                        face at -dz in z.
-    //    Double_t tl1        half-length along x of the side at +h1 in y of 
-    //                        the face at -dz in z.
-    //    Double_t apl1       Angle with respect to the y ais from the center 
-    //                        of the side at -h1 in y to the centere of the 
-    //                        side at +h1 in y of the face at -dz in z 
-    //                        [degrees].
-    //    Double_t h2         half-length along the face at +dz.
-    //    Double_t bl2        half-length along x of the side at -h2 in y of 
-    //                        the face at -dz in z.
-    //    Double_t tl2        half-length along x of the side at +h2 in y of 
-    //                        the face at +dz in z.
-    //    Double_t apl2       angle with respect to the y axis from the 
-    //                        center of the side at -h2 in y to the center 
-    //                        of the side at +h2 in y of the face at +dz in 
-    //                        z [degrees].
-    //    Int_t    med        media index number.
-    // Output:
-    //    none.
-    // Return.
-    //    none.
-    char name[5];
-    Float_t param[12];
-
-    param[0] = fScale*dz;
-    param[1] = thet;
-    param[2] = phi;
-    param[3] = twist;
-    param[4] = fScale*h1;
-    param[5] = fScale*bl1;
-    param[6] = fScale*tl1;
-    param[7] = apl1;
-    param[8] = fScale*h2;
-    param[9] = fScale*bl2;
-    param[10] = fScale*tl2;
-    param[11] = apl2;
-    G3name(gnam,name);
-    gMC->Gsvolu(name,"GTRA",GetMed(med),param,12);
-}
-//______________________________________________________________________
-void AliITSBaseGeometry::CutTube(const char *gnam,const TString &dis,
-                                Double_t rmin,Double_t rmax,Double_t dz,
-                                Double_t phi1,Double_t phi2,Double_t lx,
-                                Double_t ly,Double_t lz,Double_t hx,
-                                Double_t hy,Double_t hz,Int_t med){
-    // Interface to TMC->Gsvolu() for ITS CTUB geometries. Cut tube. A tube 
-    // cut at the extremities with planes not necessarily perpendicular to 
-    // the z axis. It has 11 parameters. See SetScale() for units. Default 
-    // units are geant 3 [cm]. phi1 should be smaller than phi2. If this is 
-    // not the case, the system adds 360 degrees to phi2.
-    // Inputs:
-    //    const char *gnam  3 character geant volume name. The letter "I"
-    //                        is appended to the front to indecate that this
-    //                        is an ITS volume.
-    //    TString &dis        String containging part discription.
-    //    Double_t rmin       Inner radius at z=0 where tube is narrowest.
-    //    Double_t rmax       Outer radius at z=0 where tube is narrowest.
-    //    Double_t dz         half-length along the z-axis
-    //    Double_t dz         half-length along the z-axis
-    //    Double_t phi1       Starting angle of the segment [degree].
-    //    Double_t phi2       Ending angle of the segment [degree].
-    //    Double_t lx         x component of a unit vector perpendicular to 
-    //                        the face at -dz.
-    //    Double_t ly         y component of a unit vector perpendicular to 
-    //                        the face at -dz.
-    //    Double_t lz         z component of a unit vector perpendicular to 
-    //                        the face at -dz.
-    //    Double_t hx         x component of a unit vector perpendicular to 
-    //                        the face at +dz.
-    //    Double_t hy         y component of a unit vector perpendicular to 
-    //                        the face at +dz.
-    //    Double_t hz         z component of a unit vector perpendicular to 
-    //                        the face at +dz.
-    //    Int_t    med        media index number.
-    // Output:
-    //    none.
-    // Return.
-    //    none.
-    char name[5];
-    Float_t param[11];
-
-    param[0] = fScale*rmin;
-    param[1] = fScale*rmax;
-    param[2] = fScale*dz;
-    param[3] = phi1;
-    param[4] = phi2;
-    param[5] = lx;
-    param[6] = ly;
-    param[7] = lz;
-    param[8] = hx;
-    param[9] = hy;
-    param[10] = hz;
-    G3name(gnam,name);
-    gMC->Gsvolu(name,"CTUB",GetMed(med),param,11);
-}
-//______________________________________________________________________
 void AliITSBaseGeometry::Pos(AliITSBaseVolParams &v,Int_t cn,
                             AliITSBaseVolParams &m,
                             TVector3 &t,Int_t irot){
@@ -1190,7 +1642,7 @@ void AliITSBaseGeometry::Matrix(Int_t irot,Double_t thet1,Double_t phi1,
     //    none.
     // Return:
     //    none.
-    Float_t t1,p1,t2,p2,t3,p3;
+    Float_t t1=0.0,p1=0.0,t2=0.0,p2=0.0,t3=0.0,p3=0.0;
 
     if(thet1==90.0&&phi1== 0.0&&
        thet2==90.0&&phi2==90.0&&
@@ -1205,6 +1657,8 @@ void AliITSBaseGeometry::Matrix(Int_t irot,Double_t thet1,Double_t phi1,
        p3 = phi3;
        fits->AliMatrix(fidrot[irot],t1,p1,t2,p2,t3,p3);
     } // end if
+    cout << "Matrix:: fidrot["<<irot<<"]="<<fidrot[irot];
+    cout <<" angles="<<t1<<" "<<p1<<" "<<t2<<" "<<p2<<" "<<t3<< " "<<p3<<endl;
 }
 //______________________________________________________________________
 void AliITSBaseGeometry::Matrix(Int_t irot,Int_t axis,Double_t thet){
@@ -1228,18 +1682,39 @@ void AliITSBaseGeometry::Matrix(Int_t irot,Int_t axis,Double_t thet){
        switch (axis) {
        case 0: //Rotate about x-axis, x-axis does not change.
            fits->AliMatrix(fidrot[irot],90.0,0.0,90.0+thet,90.0,thet,90.0);
+           /*
+           cout << "Matrix:: axis="<<axis<<" fidrot["<<irot<<"]=";
+           cout <<fidrot[irot];
+           cout <<" angles="<<90.0<<" "<<0.0<<" "<<90.0+thet<<" "<<90.0;
+           cout <<" "<<thet<< " "<<90.0<<endl;
+           */
            break;
        case 1: //Rotate about y-axis, y-axis does not change.
-           fits->AliMatrix(fidrot[irot],360.-90.0-thet,0.0,90.0,90.0,
-                           thet,90.0);
+           fits->AliMatrix(fidrot[irot],90.0-thet,0.0,90.0,90.0,-thet,0.0);
+           /*
+           cout << "Matrix:: axis="<<axis<<" fidrot["<<irot<<"]=";
+           cout << fidrot[irot];
+           cout <<" angles="<<90.-thet<<" "<<0.0<<" "<<90.0<<" "<<90.0;
+           cout <<" "<<-thet<< " "<<0.0<<endl;
+           */
            break;
        case 2: //Rotate about z-axis, z-axis does not change.
-           fits->AliMatrix(fidrot[irot],90.0,thet,90.0,360.-thet-90.0,
-                           0.0,0.0);
+           fits->AliMatrix(fidrot[irot],90.0,thet,90.0,90.+thet,0.0,0.0);
+           /*
+           cout << "Matrix:: axis="<<axis<<" fidrot["<<irot<<"]=";
+           cout <<fidrot[irot];
+           cout <<" angles="<<90.0<<" "<<thet<<" "<<90.0<<" "<<90.0+thet;
+           cout <<" "<<0.0<< " "<<0.0<<endl;
+           */
            break;
        default:
            Error("Matrix","axis must be either 0, 1, or 2. for matrix=%d",
                  irot);
+           /*
+           cout << "Matrix:: axis="<<axis<<" fidrot["<<irot<<"]=";
+           cout <<fidrot[irot];
+           cout <<" thet=" << thet<< endl;
+           */
            break;
        } // end switch
     } // end if
@@ -1258,15 +1733,14 @@ void AliITSBaseGeometry::Matrix(Int_t irot,Double_t rot[3][3]){
     //    none.
     // Return:
     //    none.
+       Double_t si,c=180./TMath::Pi();
+       Double_t ang[6]={90.0,0.0,90.0,90.0,0.0,0.0};
 
     if(rot[0][0]==1.0&&rot[1][1]==1.0&&rot[2][2]==1.0&&
        rot[0][1]==0.0&&rot[0][2]==0.0&&rot[1][0]==0.0&&
        rot[1][2]==0.0&&rot[2][0]==0.0&&rot[2][1]==0.0){
        fidrot[irot] = 0; // Unit matrix
     }else{
-       Double_t si,c=180./TMath::Pi();
-       Double_t ang[6];
-
        ang[1] = TMath::ATan2(rot[0][1],rot[0][0]);
        if(TMath::Cos(ang[1])!=0.0) si = rot[0][0]/TMath::Cos(ang[1]);
        else si = rot[0][1]/TMath::Sin(ang[1]);
@@ -1286,6 +1760,9 @@ void AliITSBaseGeometry::Matrix(Int_t irot,Double_t rot[3][3]){
        fits->AliMatrix(fidrot[irot],ang[0],ang[1],ang[2],ang[3],
                        ang[4],ang[5]);
     } // end if
+    cout << "Matrix rot[3][3]:: fidrot["<<irot<<"]="<<fidrot[irot];
+    cout <<" angles="<<ang[0]<<" "<<ang[1]<<" "<<ang[2]<<" "<<
+       ang[3]<<" "<<ang[4]<< " "<<ang[5]<<endl;
 }
 //______________________________________________________________________
 Float_t AliITSBaseGeometry::GetA(Int_t z){
@@ -1570,6 +2047,1109 @@ Double_t AliITSBaseGeometry::RadLength(Int_t iz,Double_t a){
     return xz; // [gm/cm^2]
 }
 //======================================================================
+ClassImp(AliITSBaseVolParams)
+//______________________________________________________________________
+void AliITSBaseVolParams::Print(ostream *os){
+    // Prints out the data kept in this class
+    // Inputs:
+    //    ostream *os The output stream pointer
+    // Outputs:
+    //    none.
+    // Return:
+    //    none.
+
+    *os<<"Volume Id="<<fVol<<" Copy="<<fCpn<<" Name: "<<fName<<endl;
+}
+//______________________________________________________________________
+void AliITSBaseVolParams::Read(istream *is){
+    // Read in data kept in this class
+    // Inputs:
+    //   istream *is  the input stream
+    // Outputs:
+    //   none.
+    // Return:
+    //   none.
+    char s[50];
+
+    is->get(s,10);
+    *is >> fVol;
+    is->get(s,6);
+    *is >> fCpn;
+    is->get(s,7);
+    *is >> fName;
+}
+//______________________________________________________________________
+ostream &operator<<(ostream &os,AliITSBaseVolParams &p){
+    // Operator << for C++ like output
+    // Inputs:
+    //    ostream &os            The output stream
+    //    AliITSBaseVolParams &p The class to be outputed
+    // Output:
+    //    none.
+    // Return:
+    //    ostream &os        The output stream
+
+    p.Print(&os);
+    return os;
+}
+//______________________________________________________________________
+istream &operator>>(istream &is,AliITSBaseVolParams &r){
+    // Operator << for C++ like output
+    // Inputs:
+    //    istream &is            The input stream
+    //    AliITSBaseVolParams &r The class to be read in
+    // Output:
+    //    none.
+    // Return:
+    //    istream &is        The input stream
+
+    r.Read(&is);
+    return is;
+}
+//======================================================================
+ClassImp(AliITSBoxData)
+//______________________________________________________________________
+void AliITSBoxData::Print(ostream *os){
+    // Prints out the data kept in this class
+    // Inputs:
+    //    ostream *os The output stream pointer
+    // Outputs:
+    //    none.
+    // Return:
+    //    none.
+
+#if defined __GNUC__
+#if __GNUC__ > 2
+    ios::fmtflags fmt;
+#else
+    Int_t fmt;
+#endif
+#else
+#if defined __ICC || defined __ECC
+    ios::fmtflags fmt;
+#else
+    Int_t fmt;
+#endif
+#endif
+
+    AliITSBaseVolParams::Print(os);
+    fmt = os->setf(ios::scientific);  // set scientific floating point output
+    *os << "fDx=" << fDx << " fDy=" << fDy << " fDz=" << fDz << endl;
+    os->flags(fmt); // reset back to old formating.
+    return;
+}
+//______________________________________________________________________
+void AliITSBoxData::Read(istream *is){
+    // Read in data kept in this class
+    // Inputs:
+    //   istream *is  the input stream
+    // Outputs:
+    //   none.
+    // Return:
+    //   none.
+    char s[50];
+
+    AliITSBaseVolParams::Read(is);
+    is->get(s,4);
+    *is >> fDx;
+    is->get(s,5);
+    *is >> fDy;
+    is->get(s,5);
+    *is >> fDz;
+}
+//______________________________________________________________________
+ostream &operator<<(ostream &os,AliITSBoxData &p){
+    // Operator << for C++ like output
+    // Inputs:
+    //    ostream &os      The output stream
+    //    AliITSBoxData &p The class to be outputed
+    // Output:
+    //    none.
+    // Return:
+    //    ostream &os        The output stream
+
+    p.Print(&os);
+    return os;
+}
+//______________________________________________________________________
+istream &operator>>(istream &is,AliITSBoxData &r){
+    // Operator << for C++ like output
+    // Inputs:
+    //    istream &is      The input stream
+    //    AliITSBoxData &r The class to be read in
+    // Output:
+    //    none.
+    // Return:
+    //    istream &is        The input stream
+
+    r.Read(&is);
+    return is;
+}
+//======================================================================
+ClassImp(AliITSTrapezoid1Data)
+//______________________________________________________________________
+void AliITSTrapezoid1Data::Print(ostream *os){
+    // Prints out the data kept in this class
+    // Inputs:
+    //    ostream *os The output stream pointer
+    // Outputs:
+    //    none.
+    // Return:
+    //    none.
+
+#if defined __GNUC__
+#if __GNUC__ > 2
+    ios::fmtflags fmt;
+#else
+    Int_t fmt;
+#endif
+#else
+#if defined __ICC || defined __ECC
+    ios::fmtflags fmt;
+#else
+    Int_t fmt;
+#endif
+#endif
+
+    AliITSBaseVolParams::Print(os);
+    fmt = os->setf(ios::scientific);  // set scientific floating point output
+    *os << "fDx[0]=" << fDx[0]<< " fDx[1]=" << fDx[1]  << " fDy=" << fDy;
+    *os << " fDz=" << fDz << endl;
+    os->flags(fmt); // reset back to old formating.
+    return;
+}
+//______________________________________________________________________
+void AliITSTrapezoid1Data::Read(istream *is){
+    // Read in data kept in this class
+    // Inputs:
+    //   istream *is  the input stream
+    // Outputs:
+    //   none.
+    // Return:
+    //   none.
+    char s[50];
+
+    AliITSBaseVolParams::Read(is);
+    is->get(s,7);
+    *is >> fDx[0];
+    is->get(s,8);
+    *is >> fDx[1];
+    is->get(s,5);
+    *is >> fDy;
+    is->get(s,5);
+    *is >> fDz;
+}
+//______________________________________________________________________
+ostream &operator<<(ostream &os,AliITSTrapezoid1Data &p){
+    // Operator << for C++ like output
+    // Inputs:
+    //    ostream &os      The output stream
+    //    AliITSBoxData &p The class to be outputed
+    // Output:
+    //    none.
+    // Return:
+    //    ostream &os        The output stream
+
+    p.Print(&os);
+    return os;
+}
+//______________________________________________________________________
+istream &operator>>(istream &is,AliITSTrapezoid1Data &r){
+    // Operator << for C++ like output
+    // Inputs:
+    //    istream &is       The input stream
+    //    AliITSPGonData &r The class to be read in
+    // Output:
+    //    none.
+    // Return:
+    //    istream &is        The input stream
+
+    r.Read(&is);
+    return is;
+}
+//======================================================================
+ClassImp(AliITSTrapezoid2Data)
+//______________________________________________________________________
+void AliITSTrapezoid2Data::Print(ostream *os){
+    // Prints out the data kept in this class
+    // Inputs:
+    //    ostream *os The output stream pointer
+    // Outputs:
+    //    none.
+    // Return:
+    //    none.
+
+#if defined __GNUC__
+#if __GNUC__ > 2
+    ios::fmtflags fmt;
+#else
+    Int_t fmt;
+#endif
+#else
+#if defined __ICC || defined __ECC
+    ios::fmtflags fmt;
+#else
+    Int_t fmt;
+#endif
+#endif
+
+    AliITSBaseVolParams::Print(os);
+    fmt = os->setf(ios::scientific);  // set scientific floating point output
+    *os <<  "fDx[0]=" << fDx[0]<< " fDx[1]=" << fDx[1];
+    *os << " fDy[0]=" << fDy[0] << " fDy[1]=" << fDy[1];
+    *os << " fDz=" << fDz << endl;
+    os->flags(fmt); // reset back to old formating.
+    return;
+}
+//______________________________________________________________________
+void AliITSTrapezoid2Data::Read(istream *is){
+    // Read in data kept in this class
+    // Inputs:
+    //   istream *is  the input stream
+    // Outputs:
+    //   none.
+    // Return:
+    //   none.
+    char s[50];
+
+    AliITSBaseVolParams::Read(is);
+    is->get(s,7);
+    *is >> fDx[0];
+    is->get(s,8);
+    *is >> fDx[1];
+    is->get(s,8);
+    *is >> fDy[0];
+    is->get(s,8);
+    *is >> fDy[1];
+    is->get(s,5);
+    *is >> fDz;
+}
+//______________________________________________________________________
+ostream &operator<<(ostream &os,AliITSTrapezoid2Data &p){
+    // Operator << for C++ like output
+    // Inputs:
+    //    ostream &os      The output stream
+    //    AliITSBoxData &p The class to be outputed
+    // Output:
+    //    none.
+    // Return:
+    //    ostream &os        The output stream
+
+    p.Print(&os);
+    return os;
+}
+//______________________________________________________________________
+istream &operator>>(istream &is,AliITSTrapezoid2Data &r){
+    // Operator << for C++ like output
+    // Inputs:
+    //    istream &is       The input stream
+    //    AliITSPGonData &r The class to be read in
+    // Output:
+    //    none.
+    // Return:
+    //    istream &is        The input stream
+
+    r.Read(&is);
+    return is;
+}
+//======================================================================
+ClassImp(AliITSTrapezoidData)
+//______________________________________________________________________
+void AliITSTrapezoidData::Print(ostream *os){
+    // Prints out the data kept in this class
+    // Inputs:
+    //    ostream *os The output stream pointer
+    // Outputs:
+    //    none.
+    // Return:
+    //    none.
+
+#if defined __GNUC__
+#if __GNUC__ > 2
+    ios::fmtflags fmt;
+#else
+    Int_t fmt;
+#endif
+#else
+#if defined __ICC || defined __ECC
+    ios::fmtflags fmt;
+#else
+    Int_t fmt;
+#endif
+#endif
+
+    AliITSBaseVolParams::Print(os);
+    fmt = os->setf(ios::scientific);  // set scientific floating point output
+    *os << "fTheta=" << fTheta << " fPhi=" << fPhi << " fDz=" << fDz;
+    *os << " fH[0]=" << fH[0]<< " fH[1]=" << fH[1];
+    *os << " fBl[0]=" << fBl[0] << " fBl[1]=" << fBl[1];
+    *os << " fTl[0]=" << fTl[0] << " fTl[1]=" << fTl[1];
+    *os << " fAlp[0]=" << fAlp[0] << " fAlp[1]=" << fAlp[1];
+    *os << endl;
+    os->flags(fmt); // reset back to old formating.
+    return;
+}
+//______________________________________________________________________
+void AliITSTrapezoidData::Read(istream *is){
+    // Read in data kept in this class
+    // Inputs:
+    //   istream *is  the input stream
+    // Outputs:
+    //   none.
+    // Return:
+    //   none.
+    char s[50];
+
+    AliITSBaseVolParams::Read(is);
+    is->get(s,6);
+    *is >> fTheta;
+    is->get(s,6);
+    *is >> fPhi;
+    is->get(s,5);
+    *is >> fDz;
+    is->get(s,7);
+    *is >> fH[0];
+    is->get(s,7);
+    *is >> fH[1];
+    is->get(s,8);
+    *is >> fBl[0];
+    is->get(s,8);
+    *is >> fBl[1];
+    is->get(s,8);
+    *is >> fTl[0];
+    is->get(s,8);
+    *is >> fTl[1];
+    is->get(s,9);
+    *is >> fAlp[0];
+    is->get(s,9);
+    *is >> fAlp[1];
+}
+//______________________________________________________________________
+ostream &operator<<(ostream &os,AliITSTrapezoidData &p){
+    // Operator << for C++ like output
+    // Inputs:
+    //    ostream &os      The output stream
+    //    AliITSBoxData &p The class to be outputed
+    // Output:
+    //    none.
+    // Return:
+    //    ostream &os        The output stream
+
+    p.Print(&os);
+    return os;
+}
+//______________________________________________________________________
+istream &operator>>(istream &is,AliITSTrapezoidData &r){
+    // Operator << for C++ like output
+    // Inputs:
+    //    istream &is       The input stream
+    //    AliITSPGonData &r The class to be read in
+    // Output:
+    //    none.
+    // Return:
+    //    istream &is        The input stream
+
+    r.Read(&is);
+    return is;
+}
+//======================================================================
+ClassImp(AliITSTrapezoidTwistedData)
+//______________________________________________________________________
+void AliITSTrapezoidTwistedData::Print(ostream *os){
+    // Prints out the data kept in this class
+    // Inputs:
+    //    ostream *os The output stream pointer
+    // Outputs:
+    //    none.
+    // Return:
+    //    none.
+
+#if defined __GNUC__
+#if __GNUC__ > 2
+    ios::fmtflags fmt;
+#else
+    Int_t fmt;
+#endif
+#else
+#if defined __ICC || defined __ECC
+    ios::fmtflags fmt;
+#else
+    Int_t fmt;
+#endif
+#endif
+
+    AliITSBaseVolParams::Print(os);
+    fmt = os->setf(ios::scientific);  // set scientific floating point output
+    *os << "fTheta=" << fTheta << " fPhi=" << fPhi << " fDz=" << fDz;
+    *os << " fTwist=" << fTwist;
+    *os << " fH[0]=" << fH[0]<< " fH[1]=" << fH[1];
+    *os << " fBl[0]=" << fBl[0] << " fBl[1]=" << fBl[1];
+    *os << " fTl[0]=" << fTl[0] << " fTl[1]=" << fTl[1];
+    *os << " fAlp[0]=" << fAlp[0] << " fAlp[1]=" << fAlp[1];
+    *os << endl;
+    os->flags(fmt); // reset back to old formating.
+    return;
+}
+//______________________________________________________________________
+void AliITSTrapezoidTwistedData::Read(istream *is){
+    // Read in data kept in this class
+    // Inputs:
+    //   istream *is  the input stream
+    // Outputs:
+    //   none.
+    // Return:
+    //   none.
+    char s[50];
+
+    AliITSBaseVolParams::Read(is);
+    is->get(s,6);
+    *is >> fTheta;
+    is->get(s,6);
+    *is >> fPhi;
+    is->get(s,5);
+    *is >> fDz;
+    is->get(s,8);
+    *is >> fTwist;
+    is->get(s,7);
+    *is >> fH[0];
+    is->get(s,7);
+    *is >> fH[1];
+    is->get(s,8);
+    *is >> fBl[0];
+    is->get(s,8);
+    *is >> fBl[1];
+    is->get(s,8);
+    *is >> fTl[0];
+    is->get(s,8);
+    *is >> fTl[1];
+    is->get(s,9);
+    *is >> fAlp[0];
+    is->get(s,9);
+    *is >> fAlp[1];
+}
+//______________________________________________________________________
+ostream &operator<<(ostream &os,AliITSTrapezoidTwistedData &p){
+    // Operator << for C++ like output
+    // Inputs:
+    //    ostream &os      The output stream
+    //    AliITSBoxData &p The class to be outputed
+    // Output:
+    //    none.
+    // Return:
+    //    ostream &os        The output stream
+
+    p.Print(&os);
+    return os;
+}
+//______________________________________________________________________
+istream &operator>>(istream &is,AliITSTrapezoidTwistedData &r){
+    // Operator << for C++ like output
+    // Inputs:
+    //    istream &is       The input stream
+    //    AliITSPGonData &r The class to be read in
+    // Output:
+    //    none.
+    // Return:
+    //    istream &is        The input stream
+
+    r.Read(&is);
+    return is;
+}
+//======================================================================
+ClassImp(AliITSTubeData)
+//______________________________________________________________________
+void AliITSTubeData::Print(ostream *os){
+    // Prints out the data kept in this class
+    // Inputs:
+    //    ostream *os The output stream pointer
+    // Outputs:
+    //    none.
+    // Return:
+    //    none.
+
+#if defined __GNUC__
+#if __GNUC__ > 2
+    ios::fmtflags fmt;
+#else
+    Int_t fmt;
+#endif
+#else
+#if defined __ICC || defined __ECC
+    ios::fmtflags fmt;
+#else
+    Int_t fmt;
+#endif
+#endif
+
+    AliITSBaseVolParams::Print(os);
+    fmt = os->setf(ios::scientific);  // set scientific floating point output
+    *os <<"       Z        ,      Rmin      ,      Rmax      " << endl;
+    fmt = os->setf(ios::scientific);  // set scientific floating point output
+    *os << setprecision(16) << fDz <<"\t";
+    *os << setprecision(16) << fRmin << "\t";
+    *os << setprecision(16) << fRmax << endl;
+    os->flags(fmt); // reset back to old formating.
+    return;
+}
+
+//______________________________________________________________________
+void AliITSTubeData::Read(istream *is){
+    // Read in data kept in this class
+    // Inputs:
+    //   istream *is  the input stream
+    // Outputs:
+    //   none.
+    // Return:
+    //   none.
+    char s[50];
+
+    AliITSBaseVolParams::Read(is);
+
+    is->getline(s,49);
+        *is >> fDz >> fRmin >> fRmax;
+}
+//______________________________________________________________________
+ostream &operator<<(ostream &os,AliITSTubeData &p){
+    // Operator << for C++ like output
+    // Inputs:
+    //    ostream &os       The output stream
+    //    AliITSTubeData &p The class to be outputed
+    // Output:
+    //    none.
+    // Return:
+    //    ostream &os        The output stream
+
+    p.Print(&os);
+    return os;
+}
+//______________________________________________________________________
+istream &operator>>(istream &is,AliITSTubeData &r){
+    // Operator << for C++ like output
+    // Inputs:
+    //    istream &is       The input stream
+    //    AliITSTubeData &r The class to be read in
+    // Output:
+    //    none.
+    // Return:
+    //    istream &is        The input stream
+
+    r.Read(&is);
+    return is;
+}
+//======================================================================
+ClassImp(AliITSTubeSegData)
+//______________________________________________________________________
+void AliITSTubeSegData::Print(ostream *os){
+    // Prints out the data kept in this class
+    // Inputs:
+    //    ostream *os The output stream pointer
+    // Outputs:
+    //    none.
+    // Return:
+    //    none.
+
+#if defined __GNUC__
+#if __GNUC__ > 2
+    ios::fmtflags fmt;
+#else
+    Int_t fmt;
+#endif
+#else
+#if defined __ICC || defined __ECC
+    ios::fmtflags fmt;
+#else
+    Int_t fmt;
+#endif
+#endif
+
+    AliITSBaseVolParams::Print(os);
+    fmt = os->setf(ios::scientific);  // set scientific floating point output
+    *os << "fPhi0=" << fPhi0 << " fPhi1=" << fPhi1 << endl;
+    *os <<"       Z        ,      Rmin      ,      Rmax      " << endl;
+    fmt = os->setf(ios::scientific);  // set scientific floating point output
+    *os << setprecision(16) << fDz <<"\t";
+    *os << setprecision(16) << fRmin << "\t";
+    *os << setprecision(16) << fRmax << endl;
+    os->flags(fmt); // reset back to old formating.
+    return;
+}
+//______________________________________________________________________
+void AliITSTubeSegData::Read(istream *is){
+    // Read in data kept in this class
+    // Inputs:
+    //   istream *is  the input stream
+    // Outputs:
+    //   none.
+    // Return:
+    //   none.
+    char s[50];
+
+    AliITSBaseVolParams::Read(is);
+
+    is->get(s,6);
+    *is >> fPhi0;
+    is->get(s,7);
+    *is >> fPhi1;
+    is->getline(s,49);
+       *is >> fDz >> fRmin >> fRmax;
+}
+//______________________________________________________________________
+ostream &operator<<(ostream &os,AliITSTubeSegData &p){
+    // Operator << for C++ like output
+    // Inputs:
+    //    ostream &os       The output stream
+    //    AliITSTubeData &p The class to be outputed
+    // Output:
+    //    none.
+    // Return:
+    //    ostream &os        The output stream
+
+    p.Print(&os);
+    return os;
+}
+//______________________________________________________________________
+istream &operator>>(istream &is,AliITSTubeSegData &r){
+    // Operator << for C++ like output
+    // Inputs:
+    //    istream &is       The input stream
+    //    AliITSTubeData &r The class to be read in
+    // Output:
+    //    none.
+    // Return:
+    //    istream &is        The input stream
+
+    r.Read(&is);
+    return is;
+}
+//======================================================================
+ClassImp(AliITSTubeCutData)
+//______________________________________________________________________
+void AliITSTubeCutData::Print(ostream *os){
+    // Prints out the data kept in this class
+    // Inputs:
+    //    ostream *os The output stream pointer
+    // Outputs:
+    //    none.
+    // Return:
+    //    none.
+
+#if defined __GNUC__
+#if __GNUC__ > 2
+    ios::fmtflags fmt;
+#else
+    Int_t fmt;
+#endif
+#else
+#if defined __ICC || defined __ECC
+    ios::fmtflags fmt;
+#else
+    Int_t fmt;
+#endif
+#endif
+
+    AliITSBaseVolParams::Print(os);
+    fmt = os->setf(ios::scientific);  // set scientific floating point output
+    *os << "fPhi0=" << fPhi0 << " fPhi1=" << fPhi1;
+    *os << " Norm0=("<<(fNorm[0])[0]<<" "<<(fNorm[0])[1]<<" "<<(fNorm[0])[2];
+    *os << ") Norm1=("<<(fNorm[1])[0]<<" "<<(fNorm[1])[1]<<" "<<(fNorm[1])[2];
+    *os << ")"<< endl;
+    *os <<"       Z        ,      Rmin      ,      Rmax      " << endl;
+    *os << setprecision(16) << fDz <<"\t";
+    *os << setprecision(16) << fRmin << "\t";
+    *os << setprecision(16) << fRmax << endl;
+    os->flags(fmt); // reset back to old formating.
+    return;
+}
+//______________________________________________________________________
+void AliITSTubeCutData::Read(istream *is){
+    // Read in data kept in this class
+    // Inputs:
+    //   istream *is  the input stream
+    // Outputs:
+    //   none.
+    // Return:
+    //   none.
+    char s[50];
+
+    AliITSBaseVolParams::Read(is);
+
+    is->get(s,6);
+    *is >> fPhi0;
+    is->get(s,7);
+    *is >> fPhi1;
+    is->get(s,8);
+    *is >> (fNorm[0])[0]>>(fNorm[0])[1]>>(fNorm[0])[2];
+    is->get(s,9);
+    *is >> (fNorm[1])[0]>>(fNorm[1])[1]>>(fNorm[1])[2];
+    is->get(s,1);
+    is->getline(s,49);
+    *is >> fDz >> fRmin >> fRmax;
+}
+//______________________________________________________________________
+ostream &operator<<(ostream &os,AliITSTubeCutData &p){
+    // Operator << for C++ like output
+    // Inputs:
+    //    ostream &os       The output stream
+    //    AliITSTubeData &p The class to be outputed
+    // Output:
+    //    none.
+    // Return:
+    //    ostream &os        The output stream
+
+    p.Print(&os);
+    return os;
+}
+//______________________________________________________________________
+istream &operator>>(istream &is,AliITSTubeCutData &r){
+    // Operator << for C++ like output
+    // Inputs:
+    //    istream &is       The input stream
+    //    AliITSTubeData &r The class to be read in
+    // Output:
+    //    none.
+    // Return:
+    //    istream &is        The input stream
+
+    r.Read(&is);
+    return is;
+}
+
+//======================================================================
+ClassImp(AliITSTubeEllipticalData)
+//______________________________________________________________________
+void AliITSTubeEllipticalData::Print(ostream *os){
+    // Prints out the data kept in this class
+    // Inputs:
+    //    ostream *os The output stream pointer
+    // Outputs:
+    //    none.
+    // Return:
+    //    none.
+
+#if defined __GNUC__
+#if __GNUC__ > 2
+    ios::fmtflags fmt;
+#else
+    Int_t fmt;
+#endif
+#else
+#if defined __ICC || defined __ECC
+    ios::fmtflags fmt;
+#else
+    Int_t fmt;
+#endif
+#endif
+
+    AliITSBaseVolParams::Print(os);
+    fmt = os->setf(ios::scientific);  // set scientific floating point output
+    *os <<"       Z        ,  Semi-axis-x      ,  Semi-axis-y      " << endl;
+    fmt = os->setf(ios::scientific);  // set scientific floating point output
+    *os << setprecision(16) << fDz <<"\t";
+    *os << setprecision(16) << fP0 << "\t";
+    *os << setprecision(16) << fP1 << endl;
+    os->flags(fmt); // reset back to old formating.
+    return;
+}
+
+//______________________________________________________________________
+void AliITSTubeEllipticalData::Read(istream *is){
+    // Read in data kept in this class
+    // Inputs:
+    //   istream *is  the input stream
+    // Outputs:
+    //   none.
+    // Return:
+    //   none.
+    char s[50];
+
+    AliITSBaseVolParams::Read(is);
+
+    is->getline(s,49);
+    *is >> fDz >> fP0 >> fP1;
+}
+//______________________________________________________________________
+ostream &operator<<(ostream &os,AliITSTubeEllipticalData &p){
+    // Operator << for C++ like output
+    // Inputs:
+    //    ostream &os       The output stream
+    //    AliITSTubeData &p The class to be outputed
+    // Output:
+    //    none.
+    // Return:
+    //    ostream &os        The output stream
+
+    p.Print(&os);
+    return os;
+}
+//______________________________________________________________________
+istream &operator>>(istream &is,AliITSTubeEllipticalData &r){
+    // Operator << for C++ like output
+    // Inputs:
+    //    istream &is       The input stream
+    //    AliITSTubeData &r The class to be read in
+    // Output:
+    //    none.
+    // Return:
+    //    istream &is        The input stream
+
+    r.Read(&is);
+    return is;
+}
+//======================================================================
+ClassImp(AliITSTubeHyperbolicData)
+//______________________________________________________________________
+void AliITSTubeHyperbolicData::Print(ostream *os){
+    // Prints out the data kept in this class
+    // Inputs:
+    //    ostream *os The output stream pointer
+    // Outputs:
+    //    none.
+    // Return:
+    //    none.
+
+#if defined __GNUC__
+#if __GNUC__ > 2
+    ios::fmtflags fmt;
+#else
+    Int_t fmt;
+#endif
+#else
+#if defined __ICC || defined __ECC
+    ios::fmtflags fmt;
+#else
+    Int_t fmt;
+#endif
+#endif
+
+    AliITSBaseVolParams::Print(os);
+    fmt = os->setf(ios::scientific);  // set scientific floating point output
+    *os <<"       Z               Rmin             Rmax         Theta"<<endl;
+    fmt = os->setf(ios::scientific);  // set scientific floating point output
+    *os << setprecision(16) << fDz <<"\t";
+    *os << setprecision(16) << fRmin << "\t";
+    *os << setprecision(16) << fRmax << "\t";
+    *os << setprecision(16) << fTheta << endl;
+    os->flags(fmt); // reset back to old formating.
+    return;
+}
+
+//______________________________________________________________________
+void AliITSTubeHyperbolicData::Read(istream *is){
+    // Read in data kept in this class
+    // Inputs:
+    //   istream *is  the input stream
+    // Outputs:
+    //   none.
+    // Return:
+    //   none.
+    char s[50];
+
+    AliITSBaseVolParams::Read(is);
+
+    is->getline(s,49);
+    *is >> fDz >> fRmin >> fRmax >> fTheta;
+}
+//______________________________________________________________________
+ostream &operator<<(ostream &os,AliITSTubeHyperbolicData &p){
+    // Operator << for C++ like output
+    // Inputs:
+    //    ostream &os       The output stream
+    //    AliITSTubeData &p The class to be outputed
+    // Output:
+    //    none.
+    // Return:
+    //    ostream &os        The output stream
+
+    p.Print(&os);
+    return os;
+}
+//______________________________________________________________________
+istream &operator>>(istream &is,AliITSTubeHyperbolicData &r){
+    // Operator << for C++ like output
+    // Inputs:
+    //    istream &is       The input stream
+    //    AliITSTubeData &r The class to be read in
+    // Output:
+    //    none.
+    // Return:
+    //    istream &is        The input stream
+
+    r.Read(&is);
+    return is;
+}
+//======================================================================
+ClassImp(AliITSConeData)
+//______________________________________________________________________
+void AliITSConeData::Print(ostream *os){
+    // Prints out the data kept in this class
+    // Inputs:
+    //    ostream *os The output stream pointer
+    // Outputs:
+    //    none.
+    // Return:
+    //    none.
+
+#if defined __GNUC__
+#if __GNUC__ > 2
+    ios::fmtflags fmt;
+#else
+    Int_t fmt;
+#endif
+#else
+#if defined __ICC || defined __ECC
+    ios::fmtflags fmt;
+#else
+    Int_t fmt;
+#endif
+#endif
+
+    AliITSBaseVolParams::Print(os);
+    fmt = os->setf(ios::scientific);  // set scientific floating point output
+    *os <<"       Z               Rmin             Rmax" << endl;
+    fmt = os->setf(ios::scientific);  // set scientific floating point output
+    *os << setprecision(16) << fDz <<"\t";
+    *os << setprecision(16) << fRmin0 << "\t";
+    *os << setprecision(16) << fRmax0 << endl;
+    *os << setprecision(16) << fDz <<"\t";
+    *os << setprecision(16) << fRmin1 << "\t";
+    *os << setprecision(16) << fRmax1 << endl;
+    os->flags(fmt); // reset back to old formating.
+    return;
+}
+//______________________________________________________________________
+void AliITSConeData::Read(istream *is){
+    // Read in data kept in this class
+    // Inputs:
+    //   istream *is  the input stream
+    // Outputs:
+    //   none.
+    // Return:
+    //   none.
+    char s[50];
+
+    AliITSBaseVolParams::Read(is);
+
+    is->getline(s,49);
+    *is >> fDz >> fRmin0 >> fRmax0;
+    *is >> fDz >> fRmin1 >> fRmax1;
+}
+//______________________________________________________________________
+ostream &operator<<(ostream &os,AliITSConeData &p){
+    // Operator << for C++ like output
+    // Inputs:
+    //    ostream &os       The output stream
+    //    AliITSTubeData &p The class to be outputed
+    // Output:
+    //    none.
+    // Return:
+    //    ostream &os        The output stream
+
+    p.Print(&os);
+    return os;
+}
+//______________________________________________________________________
+istream &operator>>(istream &is,AliITSConeData &r){
+    // Operator << for C++ like output
+    // Inputs:
+    //    istream &is       The input stream
+    //    AliITSTubeData &r The class to be read in
+    // Output:
+    //    none.
+    // Return:
+    //    istream &is        The input stream
+
+    r.Read(&is);
+    return is;
+}
+//======================================================================
+ClassImp(AliITSConeSegData)
+//______________________________________________________________________
+void AliITSConeSegData::Print(ostream *os){
+    // Prints out the data kept in this class
+    // Inputs:
+    //    ostream *os The output stream pointer
+    // Outputs:
+    //    none.
+    // Return:
+    //    none.
+
+#if defined __GNUC__
+#if __GNUC__ > 2
+    ios::fmtflags fmt;
+#else
+    Int_t fmt;
+#endif
+#else
+#if defined __ICC || defined __ECC
+    ios::fmtflags fmt;
+#else
+    Int_t fmt;
+#endif
+#endif
+
+    AliITSBaseVolParams::Print(os);
+    fmt = os->setf(ios::scientific);  // set scientific floating point output
+    *os << "fPhi0=" << fPhi0 << " fPhi1=" << fPhi1 << endl;
+    *os <<"       Z        ,      Rmin      ,      Rmax      " << endl;
+    fmt = os->setf(ios::scientific);  // set scientific floating point output
+    *os << setprecision(16) << fDz <<"\t";
+    *os << setprecision(16) << fRmin0 << "\t";
+    *os << setprecision(16) << fRmax0 << endl;
+    *os << setprecision(16) << fDz <<"\t";
+    *os << setprecision(16) << fRmin1 << "\t";
+    *os << setprecision(16) << fRmax1 << endl;
+    os->flags(fmt); // reset back to old formating.
+    return;
+}
+//______________________________________________________________________
+void AliITSConeSegData::Read(istream *is){
+    // Read in data kept in this class
+    // Inputs:
+    //   istream *is  the input stream
+    // Outputs:
+    //   none.
+    // Return:
+    //   none.
+    char s[50];
+
+    AliITSBaseVolParams::Read(is);
+    is->get(s,6);
+    *is >> fPhi0;
+    is->get(s,7);
+    *is >> fPhi1;
+    is->getline(s,49);
+    *is >> fDz >> fRmin0 >> fRmax0;
+    *is >> fDz >> fRmin1 >> fRmax1;
+}
+//______________________________________________________________________
+ostream &operator<<(ostream &os,AliITSConeSegData &p){
+    // Operator << for C++ like output
+    // Inputs:
+    //    ostream &os        The output stream
+    //    AliITSConeSegData &p The class to be outputed
+    // Output:
+    //    none.
+    // Return:
+    //    ostream &os        The output stream
+
+    p.Print(&os);
+    return os;
+}
+//______________________________________________________________________
+istream &operator>>(istream &is,AliITSConeSegData &r){
+    // Operator << for C++ like output
+    // Inputs:
+    //    istream &is        The input stream
+    //    AliITSConeSegData &r The class to be read in
+    // Output:
+    //    none.
+    // Return:
+    //    istream &is        The input stream
+
+    r.Read(&is);
+    return is;
+}
+//======================================================================
 ClassImp(AliITSPConeData)
 //______________________________________________________________________
 void AliITSPConeData::Print(ostream *os){
@@ -1596,7 +3176,8 @@ void AliITSPConeData::Print(ostream *os){
 #endif
 #endif
 
-    *os << "Volume "<< GetVid() << " Name: " << *GetName() << endl;
+    AliITSBaseVolParams::Print(os);
+    fmt = os->setf(ios::scientific);  // set scientific floating point output
     *os << "fNz=" << fNz << " fPhi0=" << fPhi0 << " fdPhi=" << fDphi << endl;
     *os <<"       Z        ,      Rmin      ,      Rmax      " << endl;
     fmt = os->setf(ios::scientific);  // set scientific floating point output
@@ -1619,12 +3200,8 @@ void AliITSPConeData::Read(istream *is){
     //   none.
     Int_t i;
     char s[50];
-    TString t;
 
-    is->get(s,7);
-    *is >> i; SetVid(i);
-    is->get(s,7);
-    *is >> t; SetName(t.Data());
+    AliITSBaseVolParams::Read(is);
     is->get(s,4);
     *is >> fNz;
     is->get(s,6);
@@ -1666,9 +3243,9 @@ istream &operator>>(istream &is,AliITSPConeData &r){
     return is;
 }
 //======================================================================
-ClassImp(AliITSPGonData)
+ClassImp(AliITSSphereData)
 //______________________________________________________________________
-void AliITSPGonData::Print(ostream *os){
+void AliITSSphereData::Print(ostream *os){
     // Prints out the data kept in this class
     // Inputs:
     //    ostream *os The output stream pointer
@@ -1676,7 +3253,6 @@ void AliITSPGonData::Print(ostream *os){
     //    none.
     // Return:
     //    none.
-    Int_t i;
 
 #if defined __GNUC__
 #if __GNUC__ > 2
@@ -1692,21 +3268,19 @@ void AliITSPGonData::Print(ostream *os){
 #endif
 #endif
 
-    *os << "Volume "<< GetVid() << " Name: " << *GetName() << endl;
-    *os << "fNz=" << fNz << " fNphi=" << fNphi << " fPhi0=" << fPhi0;
-    *os << " fdPhi=" << fDphi << endl;
-    *os <<"       Z        ,      Rmin      ,      Rmax      " << endl;
+    AliITSBaseVolParams::Print(os);
     fmt = os->setf(ios::scientific);  // set scientific floating point output
-    for(i=0;i<fNz;i++){
-       *os << setprecision(16) << fZ[i] <<"\t";
-       *os << setprecision(16) << fRmin[i] << "\t";
-       *os << setprecision(16) << fRmax[i] << endl;
-    } // end for i
+    *os << "fTheta[0]=" << fTheta[0] << " fTheta[1]=" << fTheta[1] << endl;
+    *os << "fPhi[0]=" << fPhi[0] << " fPhi[1]=" << fPhi[1] << endl;
+    *os <<"      Rmin      ,      Rmax      " << endl;
+    fmt = os->setf(ios::scientific);  // set scientific floating point output
+    *os << setprecision(16) << fRmin << "\t";
+    *os << setprecision(16) << fRmax << endl;
     os->flags(fmt); // reset back to old formating.
     return;
 }
 //______________________________________________________________________
-void AliITSPGonData::Read(istream *is){
+void AliITSSphereData::Read(istream *is){
     // Read in data kept in this class
     // Inputs:
     //   istream *is  the input stream
@@ -1714,36 +3288,112 @@ void AliITSPGonData::Read(istream *is){
     //   none.
     // Return:
     //   none.
-    Int_t i;
     char s[50];
-    TString t;
 
-    is->get(s,7);
-    *is >> i;SetVid(i);
-    is->get(s,7);
-    *is >> t; SetName(t.Data());
-    
-    is->get(s,4);
-    *is >> fNz;
-    is->get(s,6);
-    *is >> fNphi;
-    is->get(s,6);
-    *is >> fPhi0;
-    is->get(s,6);
-    *is >> fDphi;
+    AliITSBaseVolParams::Read(is);
+    is->get(s,10);
+    *is >> fTheta[0];
+    is->get(s,11);
+    *is >> fTheta[1];
+    is->get(s,8);
+    *is >> fPhi[0];
+    is->get(s,9);
+    *is >> fPhi[1];
     is->getline(s,49);
+    *is >>fRmin >> fRmax;
+}
+//______________________________________________________________________
+ostream &operator<<(ostream &os,AliITSSphereData &p){
+    // Operator << for C++ like output
+    // Inputs:
+    //    ostream &os        The output stream
+    //    AliITSPConeData &p The class to be outputed
+    // Output:
+    //    none.
+    // Return:
+    //    ostream &os        The output stream
 
-    Size(fNz);
-    for(i=0;i<fNz;i++){
-       *is >> fZ[i] >> fRmin[i] >> fRmax[i];
-    } // end for i
+    p.Print(&os);
+    return os;
 }
 //______________________________________________________________________
-ostream &operator<<(ostream &os,AliITSPGonData &p){
+istream &operator>>(istream &is,AliITSSphereData &r){
     // Operator << for C++ like output
     // Inputs:
-    //    ostream &os       The output stream
-    //    AliITSPGonData &p The class to be outputed
+    //    istream &is        The input stream
+    //    AliITSPConeData &r The class to be read in
+    // Output:
+    //    none.
+    // Return:
+    //    istream &is        The input stream
+
+    r.Read(&is);
+    return is;
+}
+//======================================================================
+ClassImp(AliITSParallelpipedData)
+//______________________________________________________________________
+void AliITSParallelpipedData::Print(ostream *os){
+    // Prints out the data kept in this class
+    // Inputs:
+    //    ostream *os The output stream pointer
+    // Outputs:
+    //    none.
+    // Return:
+    //    none.
+
+#if defined __GNUC__
+#if __GNUC__ > 2
+    ios::fmtflags fmt;
+#else
+    Int_t fmt;
+#endif
+#else
+#if defined __ICC || defined __ECC
+    ios::fmtflags fmt;
+#else
+    Int_t fmt;
+#endif
+#endif
+
+    AliITSBaseVolParams::Print(os);
+    fmt = os->setf(ios::scientific);  // set scientific floating point output
+    *os << "fDx=" << fDx << " fDy=" << fDy << " fDz=" << fDz << endl;
+    *os << "fAlpha=" << fAlpha << " fTheta=" << fTheta <<" fPhi="<<fPhi<<endl;
+    os->flags(fmt); // reset back to old formating.
+    return;
+}
+//______________________________________________________________________
+void AliITSParallelpipedData::Read(istream *is){
+    // Read in data kept in this class
+    // Inputs:
+    //   istream *is  the input stream
+    // Outputs:
+    //   none.
+    // Return:
+    //   none.
+    char s[50];
+
+    AliITSBaseVolParams::Read(is);
+    is->get(s,4);
+    *is >> fDx;
+    is->get(s,5);
+    *is >> fDy;
+    is->get(s,5);
+    *is >> fDz;
+    is->get(s,7);
+    *is >> fAlpha;
+    is->get(s,8);
+    *is >> fTheta;
+    is->get(s,6);
+    *is >> fPhi;
+}
+//______________________________________________________________________
+ostream &operator<<(ostream &os,AliITSParallelpipedData &p){
+    // Operator << for C++ like output
+    // Inputs:
+    //    ostream &os      The output stream
+    //    AliITSBoxData &p The class to be outputed
     // Output:
     //    none.
     // Return:
@@ -1753,11 +3403,11 @@ ostream &operator<<(ostream &os,AliITSPGonData &p){
     return os;
 }
 //______________________________________________________________________
-istream &operator>>(istream &is,AliITSPGonData &r){
+istream &operator>>(istream &is,AliITSParallelpipedData &r){
     // Operator << for C++ like output
     // Inputs:
-    //    istream &is       The input stream
-    //    AliITSPGonData &r The class to be read in
+    //    istream &is      The input stream
+    //    AliITSBoxData &r The class to be read in
     // Output:
     //    none.
     // Return:
@@ -1767,9 +3417,9 @@ istream &operator>>(istream &is,AliITSPGonData &r){
     return is;
 }
 //======================================================================
-ClassImp(AliITSTubeData)
+ClassImp(AliITSPGonData)
 //______________________________________________________________________
-void AliITSTubeData::Print(ostream *os){
+void AliITSPGonData::Print(ostream *os){
     // Prints out the data kept in this class
     // Inputs:
     //    ostream *os The output stream pointer
@@ -1777,6 +3427,7 @@ void AliITSTubeData::Print(ostream *os){
     //    none.
     // Return:
     //    none.
+    Int_t i;
 
 #if defined __GNUC__
 #if __GNUC__ > 2
@@ -1792,17 +3443,22 @@ void AliITSTubeData::Print(ostream *os){
 #endif
 #endif
 
-    *os << "Volume "<< GetVid() << " Name: " << *GetName() << endl;
-    *os <<"       Z        ,      Rmin      ,      Rmax      " << endl;
+    AliITSBaseVolParams::Print(os);
     fmt = os->setf(ios::scientific);  // set scientific floating point output
-    *os << setprecision(16) << fDz <<"\t";
-    *os << setprecision(16) << fRmin << "\t";
-    *os << setprecision(16) << fRmax << endl;
+    fmt = os->setf(ios::scientific);  // set scientific floating point output
+    *os << "fNz=" << fNz << " fNphi=" << fNphi << " fPhi0=" << fPhi0;
+    *os << " fdPhi=" << fDphi << endl;
+    *os <<"       Z        ,      Rmin      ,      Rmax      " << endl;
+    for(i=0;i<fNz;i++){
+       *os << setprecision(16) << fZ[i] <<"\t";
+       *os << setprecision(16) << fRmin[i] << "\t";
+       *os << setprecision(16) << fRmax[i] << endl;
+    } // end for i
     os->flags(fmt); // reset back to old formating.
     return;
 }
 //______________________________________________________________________
-void AliITSTubeData::Read(istream *is){
+void AliITSPGonData::Read(istream *is){
     // Read in data kept in this class
     // Inputs:
     //   istream *is  the input stream
@@ -1812,22 +3468,30 @@ void AliITSTubeData::Read(istream *is){
     //   none.
     Int_t i;
     char s[50];
-    TString t;
 
-    is->get(s,7);
-    *is >> i;SetVid(i);
-    is->get(s,7);
-    *is >> t; SetName(t.Data());
+    AliITSBaseVolParams::Read(is);
     
+    is->get(s,4);
+    *is >> fNz;
+    is->get(s,6);
+    *is >> fNphi;
+    is->get(s,6);
+    *is >> fPhi0;
+    is->get(s,6);
+    *is >> fDphi;
     is->getline(s,49);
-       *is >> fDz >> fRmin >> fRmax;
+
+    Size(fNz);
+    for(i=0;i<fNz;i++){
+       *is >> fZ[i] >> fRmin[i] >> fRmax[i];
+    } // end for i
 }
 //______________________________________________________________________
-ostream &operator<<(ostream &os,AliITSTubeData &p){
+ostream &operator<<(ostream &os,AliITSPGonData &p){
     // Operator << for C++ like output
     // Inputs:
     //    ostream &os       The output stream
-    //    AliITSTubeData &p The class to be outputed
+    //    AliITSPGonData &p The class to be outputed
     // Output:
     //    none.
     // Return:
@@ -1837,11 +3501,11 @@ ostream &operator<<(ostream &os,AliITSTubeData &p){
     return os;
 }
 //______________________________________________________________________
-istream &operator>>(istream &is,AliITSTubeData &r){
+istream &operator>>(istream &is,AliITSPGonData &r){
     // Operator << for C++ like output
     // Inputs:
     //    istream &is       The input stream
-    //    AliITSTubeData &r The class to be read in
+    //    AliITSPGonData &r The class to be read in
     // Output:
     //    none.
     // Return:
index 4423ee4..7bbd93c 100644 (file)
@@ -15,24 +15,279 @@ class TString;
 
 class AliITSBaseVolParams : public TObject {
  public:
-    AliITSBaseVolParams(){fVol=-1;} // Default constructor
+    AliITSBaseVolParams(){fVol=-1;fCpn=-1;} // Default constructor
     virtual ~AliITSBaseVolParams(){} // Destructor
     virtual void SetVid(Int_t i){fVol = i;}// Sets the volume id from the next
                                           // available one
     virtual Int_t GetVid(){return fVol;} // Returns volume id
     virtual void SetName(const char *c){fName=c;}//Sets name of this volume
     virtual TString* GetName(){return &fName;} // Returns volume name
-    virtual void Print(ostream *os){} // Prints output content of this class
-    virtual void Read(istream *is){} // Reads output created by Print above.
+    virtual void Print(ostream *os); // Prints output content of this class
+    virtual void Read(istream *is); // Reads output created by Print above.
+    virtual Int_t GetCopyNumber(){return fCpn;}//Returns existing max copy no.
+    // incoment and return copy number.
+    virtual Int_t GetNextcpn(){fCpn++;return fCpn;}
+    // Same as above, but add's 1 for FORTRAN/Geant3 compliance.
+    virtual Int_t GetG3cpn(){return GetNextcpn()+1;}
+    /*
+    Int_t ITSIndexToITSG3name(const Int_t i)
+       {return AliITSBaseGeometry::ITSIndexToITSG3name(i);};
+    Int_t ITSG3VnameToIndex(const char *name)const
+       {return AliITSBaseGeometry::ITSG3VnameToIndex(name);};
+    */
  private:
     Int_t   fVol;  // Volume index number
+    Int_t   fCpn;  // max Copy number
     TString fName; // Volume name
 
     ClassDef(AliITSBaseVolParams,1) // Basic ITS volume parameters class
 };
+// Input and output function for standard C++ input/output.
+ostream &operator<<(ostream &os,AliITSBaseVolParams &source);
+istream &operator>>(istream &os,AliITSBaseVolParams &source);
+#endif
+
+#ifndef ALIITSBOXDATA_H
+#define ALIITSBOXDATA_H
+/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ * See cxx source for full Copyright notice                               */
+
+/*
+  $Id$
+ */
+
+/////////////////////////////////////////////////////////////////////////
+//  Geant 3 Box data structure.
+/////////////////////////////////////////////////////////////////////////
+class TString;
+
+class AliITSBoxData : public AliITSBaseVolParams {
+ public:
+    AliITSBoxData() : AliITSBaseVolParams() // Default constructor
+       {fDx=0.0;fDy=0.0;fDz=0;}
+    virtual ~AliITSBoxData(){} // Standard destructor
+    // Getters
+    Double_t DxAt() {return fDx;} // Returm arrau of rmin values
+    Double_t DyAt() {return fDy;} // Returm arrau of rmax values
+    Double_t DzAt() {return fDz;} // Return fDz coordiante
+    Double_t& Dx(){return fDx;}// Returns address of fRmax
+    Double_t& Dy(){return fDy;}// Returns address of fRmin
+    Double_t& Dz(){return fDz;}// Returns address of fDz
+    void Print(ostream *os); // Prints output content of this class
+    void Read(istream *is); // Reads output created by Print above.
+ private:
+    Double_t fDx; // X half length
+    Double_t fDy; // Y half length
+    Double_t fDz; // Z half length.
+
+    ClassDef(AliITSBoxData,1) // Box data class
+};
+// Input and output function for standard C++ input/output.
+ostream &operator<<(ostream &os,AliITSBoxData &source);
+istream &operator>>(istream &os,AliITSBoxData &source);
+#endif
+
+#ifndef ALIITSTRAPEZOID1DATA_H
+#define ALIITSTRAPEZOID1DATA_H
+/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ * See cxx source for full Copyright notice                               */
+
+/*
+  $Id$
+ */
+
+/////////////////////////////////////////////////////////////////////////
+//  Geant 3 Trapezoid 1 data structure.
+/////////////////////////////////////////////////////////////////////////
+class TString;
+
+class AliITSTrapezoid1Data : public AliITSBaseVolParams {
+ public:
+    AliITSTrapezoid1Data() : AliITSBaseVolParams() // Default constructor
+       {fDx[0]=0.0;fDx[1]=0.0;fDy=0.0;fDz=0;}
+    virtual ~AliITSTrapezoid1Data(){} // Standard destructor
+    // Getters
+    Double_t DxAt(Int_t i) {return fDx[i];} // Returm arrau of rmin values
+    Double_t DyAt() {return fDy;} // Returm arrau of rmax values
+    Double_t DzAt() {return fDz;} // Return fDz coordiante
+    Double_t& Dx(Int_t i){return fDx[i];}// Returns address of fRmax
+    Double_t& Dy(){return fDy;}// Returns address of fRmin
+    Double_t& Dz(){return fDz;}// Returns address of fDz
+    void Print(ostream *os); // Prints output content of this class
+    void Read(istream *is); // Reads output created by Print above.
+ private:
+    Double_t fDx[2]; // X half length
+    Double_t fDy; // Y half length
+    Double_t fDz; // Z half length.
+
+    ClassDef(AliITSTrapezoid1Data,1) // Trapezoid 1 data class
+};
+// Input and output function for standard C++ input/output.
+ostream &operator<<(ostream &os,AliITSTrapezoid1Data &source);
+istream &operator>>(istream &os,AliITSTrapezoid1Data &source);
+#endif
+
+#ifndef ALIITSTRAPEZOID2DATA_H
+#define ALIITSTRAPEZOID2DATA_H
+/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ * See cxx source for full Copyright notice                               */
+
+/*
+  $Id$
+ */
+
+/////////////////////////////////////////////////////////////////////////
+//  Geant 3 Trapezoid 2 data structure.
+/////////////////////////////////////////////////////////////////////////
+class TString;
 
+class AliITSTrapezoid2Data : public AliITSBaseVolParams {
+ public:
+    AliITSTrapezoid2Data() : AliITSBaseVolParams() // Default constructor
+       {fDx[0]=0.0;fDx[1]=0.0;fDy[0]=0.0;fDy[1]=0.0;fDz=0.0;}
+    virtual ~AliITSTrapezoid2Data(){} // Standard destructor
+    // Getters
+    Double_t DxAt(Int_t i) {return fDx[i];} // Returm arrau of rmin values
+    Double_t DyAt(Int_t i) {return fDy[i];} // Returm arrau of rmax values
+    Double_t DzAt() {return fDz;} // Return fDz coordiante
+    Double_t& Dx(Int_t i){return fDx[i];}// Returns address of fRmax
+    Double_t& Dy(Int_t i){return fDy[i];}// Returns address of fRmin
+    Double_t& Dz(){return fDz;}// Returns address of fDz
+    void Print(ostream *os); // Prints output content of this class
+    void Read(istream *is); // Reads output created by Print above.
+ private:
+    Double_t fDx[2]; // X half length
+    Double_t fDy[2]; // Y half length
+    Double_t fDz; // Z half length.
+
+    ClassDef(AliITSTrapezoid2Data,1) // Trapezoid 2 data class
+};
+// Input and output function for standard C++ input/output.
+ostream &operator<<(ostream &os,AliITSTrapezoid2Data &source);
+istream &operator>>(istream &os,AliITSTrapezoid2Data &source);
 #endif
 
+#ifndef ALIITSTRAPEZOIDDATA_H
+#define ALIITSTRAPEZOIDDATA_H
+/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ * See cxx source for full Copyright notice                               */
+
+/*
+  $Id$
+ */
+
+/////////////////////////////////////////////////////////////////////////
+//  Geant 3 Trapezoid General data structure.
+/////////////////////////////////////////////////////////////////////////
+class TString;
+
+class AliITSTrapezoidData : public AliITSBaseVolParams {
+ public:
+    AliITSTrapezoidData() : AliITSBaseVolParams() // Default constructor
+       {fDz=0.0;fTheta=0.0;fPhi=0.0;fH[0]=0.0;fH[1]=0.0;fBl[0]=0.0;
+       fBl[1]=0.0;fTl[0]=0.0;fTl[1]=0.0;fAlp[0]=0.0;fAlp[1]=0.0;}
+    virtual ~AliITSTrapezoidData(){} // Standard destructor
+    // Getters
+    Double_t DzAt() {return fDz;} // Returm arrau of dZ values
+    Double_t Theta() {return fTheta;} // Returm arrau of Theta values
+    Double_t Phi() {return fPhi;} // Return Phi coordiante
+    Double_t HAt(Int_t i){return fH[i];} // Return fH[]
+    Double_t Bl(Int_t i){return fBl[i];} // Return fBl[]
+    Double_t Tl(Int_t i){return fTl[i];} // Return fTl[]
+    Double_t Alpha(Int_t i){return fAlp[i];} // Return fAlpha[]
+    Double_t& Dz(){return fDz;}// Returns address of fDz
+    Double_t& Th(){return fTheta;}// Returns address of fTheta
+    Double_t& Ph(){return fPhi;}// Returns address of fPhi
+    Double_t& H(Int_t i){return fH[i];}// Returns address of fH
+    Double_t& B(Int_t i){return fBl[i];}// Returns address of fBl
+    Double_t& T(Int_t i){return fTl[i];}// Returns address of fTl
+    Double_t& A(Int_t i){return fAlp[i];}// Returns address of fAlp
+    void Print(ostream *os); // Prints output content of this class
+    void Read(istream *is); // Reads output created by Print above.
+ private:
+    Double_t fTheta; // Polar angle of the line joing the center of the 
+                     // face at -dz to the center of the one at dz [degree].
+    Double_t fPhi;   // aximuthal angle of the line joing the center of the
+                     //  face at -dz to the center of the one at +dz [degree].
+    Double_t fDz;    // Half-length along the z-asix
+    Double_t fH[2];  // half-length along y of the face at -dz & _dz.
+    Double_t fBl[2]; // half-length along x of the side at -h1 in y of 
+                     // the face at -dz & +dz in z.
+    Double_t fTl[2]; // half-length along x of teh side at +h1 in y of 
+                     // the face at -dz & +dz in z.
+    Double_t fAlp[2];// angle with respect to the y axis from the center of
+                     // the side at -fH[0] & +fH[1] in y to the cetner of the 
+                     // side at +h1 in y of the face at -dz in z [degree].
+
+    ClassDef(AliITSTrapezoidData,1) // General Trapezoid data class
+};
+// Input and output function for standard C++ input/output.
+ostream &operator<<(ostream &os,AliITSTrapezoidData &source);
+istream &operator>>(istream &os,AliITSTrapezoidData &source);
+#endif
+
+#ifndef ALIITSTRAPEZOIDTWISTEDDATA_H
+#define ALIITSTRAPEZOIDTWISTEDDATA_H
+/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ * See cxx source for full Copyright notice                               */
+
+/*
+  $Id$
+ */
+
+/////////////////////////////////////////////////////////////////////////
+//  Geant 3 Trapezoid General data structure.
+/////////////////////////////////////////////////////////////////////////
+class TString;
+
+class AliITSTrapezoidTwistedData : public AliITSBaseVolParams {
+ public:
+    AliITSTrapezoidTwistedData() : AliITSBaseVolParams() //Default constructor
+       {fTwist=0.0;fDz=0.0;fTheta=0.0;fPhi=0.0;fH[0]=0.0;fH[1]=0.0;
+       fBl[0]=0.0;fBl[1]=0.0;fTl[0]=0.0;fTl[1]=0.0;fAlp[0]=0.0;fAlp[1]=0.0;}
+    virtual ~AliITSTrapezoidTwistedData(){} // Standard destructor
+    // Getters
+    Double_t DzAt() {return fDz;} // Returm arrau of dZ values
+    Double_t Theta() {return fTheta;} // Returm arrau of Theta values
+    Double_t Phi() {return fPhi;} // Return Phi angle
+    Double_t Twist(){return fTwist;} // Returns Twist angle
+    Double_t HAt(Int_t i){return fH[i];} // Return fH[]
+    Double_t Bl(Int_t i){return fBl[i];} // Return fBl[]
+    Double_t Tl(Int_t i){return fTl[i];} // Return fTl[]
+    Double_t Alpha(Int_t i){return fAlp[i];} // Return fAlpha[]
+    Double_t& Dz(){return fDz;}// Returns address of fDz
+    Double_t& Th(){return fTheta;}// Returns address of fTheta
+    Double_t& Ph(){return fPhi;}// Returns address of fPhi
+    Double_t& Tw(){return fTwist;}// Returns address of fTwist
+    Double_t& H(Int_t i){return fH[i];}// Returns address of fH
+    Double_t& B(Int_t i){return fBl[i];}// Returns address of fBl
+    Double_t& T(Int_t i){return fTl[i];}// Returns address of fTl
+    Double_t& A(Int_t i){return fAlp[i];}// Returns address of fAlp
+    void Print(ostream *os); // Prints output content of this class
+    void Read(istream *is); // Reads output created by Print above.
+ private:
+    Double_t fTwist; // Twist angle of the faces parallel to the x-y 
+                     // plane at z=+-dz around an axis parallel to z
+    Double_t fTheta; // Polar angle of the line joing the center of the 
+                     // face at -dz to the center of the one at dz [degree].
+    Double_t fPhi;   // aximuthal angle of the line joing the center of the
+                     //  face at -dz to the center of the one at +dz [degree].
+    Double_t fDz;    // Half-length along the z-asix
+    Double_t fH[2];  // half-length along y of the face at -dz & _dz.
+    Double_t fBl[2]; // half-length along x of the side at -h1 in y of 
+                     // the face at -dz & +dz in z.
+    Double_t fTl[2]; // half-length along x of teh side at +h1 in y of 
+                     // the face at -dz & +dz in z.
+    Double_t fAlp[2];// angle with respect to the y axis from the center of
+                     // the side at -fH[0] & +fH[1] in y to the cetner of the 
+                     // side at +h1 in y of the face at -dz in z [degree].
+
+    ClassDef(AliITSTrapezoidTwistedData,1) // Twisted Trapezoid data class
+};
+// Input and output function for standard C++ input/output.
+ostream &operator<<(ostream &os,AliITSTrapezoidTwistedData &source);
+istream &operator>>(istream &os,AliITSTrapezoidTwistedData &source);
+#endif
 
 #ifndef ALIITSTUBEDATA_H
 #define ALIITSTUBEDATA_H
@@ -44,9 +299,8 @@ class AliITSBaseVolParams : public TObject {
  */
 
 /////////////////////////////////////////////////////////////////////////
-//  Geant 3 Poly-Cone data structure.
+//  Geant 3 Tube data structure.
 /////////////////////////////////////////////////////////////////////////
-#include <math.h> // for the definision of NAN.
 class TString;
 
 class AliITSTubeData : public AliITSBaseVolParams {
@@ -55,16 +309,12 @@ class AliITSTubeData : public AliITSBaseVolParams {
        {fDz=0;fRmin=0;fRmax=0;}
     virtual ~AliITSTubeData(){} // Standard destructor
     // Getters
-    Double_t DzAt() // Return fDz coordiante
-       {return fDz;}
+    Double_t DzAt(){return fDz;} // Return fDz coordiante
     Double_t Rmin() {return fRmin;} // Returm arrau of rmin values
     Double_t Rmax() {return fRmax;} // Returm arrau of rmax values
-    Double_t& Z()// Returns address of fDz
-       {return fDz;}
-    Double_t& Rn()// Returns address of fRmin
-       {return fRmin;}
-    Double_t& Rx()// Returns address of fRmax
-       {return fRmax;}
+    Double_t& Z(){return fDz;}// Returns address of fDz
+    Double_t& Rn(){return fRmin;}// Returns address of fRmin
+    Double_t& Rx(){return fRmax;}// Returns address of fRmax
     void Print(ostream *os); // Prints output content of this class
     void Read(istream *is); // Reads output created by Print above.
  private:
@@ -72,15 +322,295 @@ class AliITSTubeData : public AliITSBaseVolParams {
     Double_t fRmax; // Outer Radius
     Double_t fDz;   // Z half length.
 
-    ClassDef(AliITSTubeData,1) // Poly Cone data class
+    ClassDef(AliITSTubeData,1) // Simple Tube data class
 };
 // Input and output function for standard C++ input/output.
 ostream &operator<<(ostream &os,AliITSTubeData &source);
 istream &operator>>(istream &os,AliITSTubeData &source);
+#endif
 
+#ifndef ALIITSTUBESEGDATA_H
+#define ALIITSTUBESEGDATA_H
+/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ * See cxx source for full Copyright notice                               */
+
+/*
+  $Id$
+ */
+
+/////////////////////////////////////////////////////////////////////////
+//  Geant 3 Tube Segment data structure.
+/////////////////////////////////////////////////////////////////////////
+class TString;
+
+class AliITSTubeSegData : public AliITSBaseVolParams {
+ public:
+    AliITSTubeSegData() : AliITSBaseVolParams() // Default constructor
+       {fDz=0;fRmin=0;fRmax=0;}
+    virtual ~AliITSTubeSegData(){} // Standard destructor
+    // Getters
+    Double_t Phi0() {return fPhi0;} // Returns starting Phi value
+    Double_t Phi1() {return fPhi1;} // Returns endinging Phi value
+    Double_t DzAt(){return fDz;} // Return fDz coordiante
+    Double_t Rmin() {return fRmin;} // Returm arrau of rmin values
+    Double_t Rmax() {return fRmax;} // Returm arrau of rmax values
+    Double_t& P0(){return fPhi0;}// Returns addres of fPhi0
+    Double_t& P1(){return fPhi1;}// Returns addres of fPhi1
+    Double_t& Z(){return fDz;}// Returns address of fDz
+    Double_t& Rn(){return fRmin;}// Returns address of fRmin
+    Double_t& Rx(){return fRmax;}// Returns address of fRmax
+    void Print(ostream *os); // Prints output content of this class
+    void Read(istream *is); // Reads output created by Print above.
+ private:
+    Double_t fPhi0,fPhi1;   // Starting and ending phi angles [degrees]
+    Double_t fRmin; // Inner Radius
+    Double_t fRmax; // Outer Radius
+    Double_t fDz;   // Z half length.
 
+    ClassDef(AliITSTubeSegData,1) // Segment of a Tube data class
+};
+// Input and output function for standard C++ input/output.
+ostream &operator<<(ostream &os,AliITSTubeSegData &source);
+istream &operator>>(istream &os,AliITSTubeSegData &source);
+#endif
+
+#ifndef ALIITSTUBECUTDATA_H
+#define ALIITSTUBECUTDATA_H
+/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ * See cxx source for full Copyright notice                               */
+
+/*
+  $Id$
+ */
+
+#include <TVector3.h>
+
+/////////////////////////////////////////////////////////////////////////
+//  Geant 3 Tube Cut data structure.
+/////////////////////////////////////////////////////////////////////////
+class TString;
+
+class AliITSTubeCutData : public AliITSBaseVolParams {
+ public:
+    AliITSTubeCutData() : AliITSBaseVolParams() // Default constructor
+       {fDz=0;fRmin=0;fRmax=0;fPhi0=0.0,fPhi1=0.0;}
+    virtual ~AliITSTubeCutData(){} // Standard destructor
+    // Getters
+    Double_t Phi0() {return fPhi0;} // Returns starting Phi value
+    Double_t Phi1() {return fPhi1;} // Returns endinging Phi value
+    Double_t DzAt(){return fDz;} // Return fDz coordiante
+    Double_t Rmin() {return fRmin;} // Returm arrau of rmin values
+    Double_t Rmax() {return fRmax;} // Returm arrau of rmax values
+    TVector3 Normal(Int_t i){return (fNorm[i]).Unit();}// Returns unit normal
+                                                     //at z serface -dz & +dz.
+    // Returns jth unit normal at z serface -dz & +dz.
+    Double_t Normal(Int_t i,Int_t j){return ((fNorm[i]).Unit())[j];}
+    Double_t& P0(){return fPhi0;}// Returns addres of fPhi0
+    Double_t& P1(){return fPhi1;}// Returns addres of fPhi1
+    Double_t& Z(){return fDz;}// Returns address of fDz
+    Double_t& Rn(){return fRmin;}// Returns address of fRmin
+    Double_t& Rx(){return fRmax;}// Returns address of fRmax
+    // Returns address of normal at z serface -dz & +dz.
+    TVector3& Nv(Int_t i){return fNorm[i];}
+    Double_t& Nx(Int_t i){return (fNorm[i])[0];}// return address of x
+                                                //component of the nomal vect.
+    Double_t& Ny(Int_t i){return (fNorm[i])[1];}// return address of y
+                                                //component of the nomal vect.
+    Double_t& Nz(Int_t i){return (fNorm[i])[2];}// return address of z
+                                                //component of the nomal vect.
+    void Print(ostream *os); // Prints output content of this class
+    void Read(istream *is); // Reads output created by Print above.
+ private:
+    Double_t fPhi0,fPhi1;   // Starting and ending phi angles [degrees]
+    Double_t fRmin; // Inner Radius
+    Double_t fRmax; // Outer Radius
+    Double_t fDz;   // Z half length.
+    TVector3 fNorm[2]; // unit vector normal to -dz and +dz surfaces.
+
+    ClassDef(AliITSTubeCutData,1) // A tube segment with cut ends data class
+};
+// Input and output function for standard C++ input/output.
+ostream &operator<<(ostream &os,AliITSTubeCutData &source);
+istream &operator>>(istream &os,AliITSTubeCutData &source);
 #endif
-   
+
+#ifndef ALIITSTUBEELLIPTICALDATA_H
+#define ALIITSTUBEELLIPTICALDATA_H
+/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ * See cxx source for full Copyright notice                               */
+
+/*
+  $Id$
+ */
+
+/////////////////////////////////////////////////////////////////////////
+//  Geant 3 Tube Elliptical data structure.
+/////////////////////////////////////////////////////////////////////////
+class TString;
+
+class AliITSTubeEllipticalData : public AliITSBaseVolParams {
+ public:
+    AliITSTubeEllipticalData() : AliITSBaseVolParams() // Default constructor
+       {fDz=0;fP0=0;fP1=0;}
+    virtual ~AliITSTubeEllipticalData(){} // Standard destructor
+    // Getters
+    Double_t DzAt(){return fDz;} // Return fDz coordiante
+    Double_t P0At() {return fP0;} // Returm fP0 values
+    Double_t P1At() {return fP1;} // Returm fP1 values
+    Double_t& Z(){return fDz;}// Returns address of fDz
+    Double_t& P0(){return fP0;}// Returns address of fP0
+    Double_t& P1(){return fP1;}// Returns address of fP1
+    void Print(ostream *os); // Prints output content of this class
+    void Read(istream *is); // Reads output created by Print above.
+ private:
+    Double_t fP0; // semi-axis of the elipse along x
+    Double_t fP1; // semi-asis of the elipse along y
+    Double_t fDz;   // Z half length.
+
+    ClassDef(AliITSTubeEllipticalData,1) // Tube with an elliptical cross
+                                        // section data class
+};
+// Input and output function for standard C++ input/output.
+ostream &operator<<(ostream &os,AliITSTubeEllipticalData &source);
+istream &operator>>(istream &os,AliITSTubeEllipticalData &source);
+#endif
+
+#ifndef ALIITSTUBEHYPERBOLICDATA_H
+#define ALIITSTUBEHYPERBOLICDATA_H
+/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ * See cxx source for full Copyright notice                               */
+
+/*
+  $Id$
+ */
+
+/////////////////////////////////////////////////////////////////////////
+//  Geant 3 Tube Hyperbolic data structure.
+/////////////////////////////////////////////////////////////////////////
+class TString;
+
+class AliITSTubeHyperbolicData : public AliITSBaseVolParams {
+ public:
+    AliITSTubeHyperbolicData() : AliITSBaseVolParams() // Default constructor
+       {fTheta=0.0;fDz=0;fRmin=0;fRmax=0;}
+    virtual ~AliITSTubeHyperbolicData(){} // Standard destructor
+    // Getters
+    Double_t ThetaAt() {return fTheta;} // Returns starting Theta value
+    Double_t DzAt(){return fDz;} // Return fDz coordiante
+    Double_t Rmin() {return fRmin;} // Returm arrau of rmin values
+    Double_t Rmax() {return fRmax;} // Returm arrau of rmax values
+    Double_t& Theta(){return fTheta;}// Returns addres of fTheta
+    Double_t& Z(){return fDz;}// Returns address of fDz
+    Double_t& Rn(){return fRmin;}// Returns address of fRmin
+    Double_t& Rx(){return fRmax;}// Returns address of fRmax
+    void Print(ostream *os); // Prints output content of this class
+    void Read(istream *is); // Reads output created by Print above.
+ private:
+    Double_t fTheta;   // stero angle of rotation of the tetwo faces [degrees]
+    Double_t fRmin; // Inner Radius
+    Double_t fRmax; // Outer Radius
+    Double_t fDz;   // Z half length.
+
+    ClassDef(AliITSTubeHyperbolicData,1) // Hyperbolic Tube data class
+};
+// Input and output function for standard C++ input/output.
+ostream &operator<<(ostream &os,AliITSTubeHyperbolicData &source);
+istream &operator>>(istream &os,AliITSTubeHyperbolicData &source);
+#endif
+
+#ifndef ALIITSCONEDATA_H
+#define ALIITSCONEDATA_H
+/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ * See cxx source for full Copyright notice                               */
+
+/*
+  $Id$
+ */
+
+/////////////////////////////////////////////////////////////////////////
+//  Geant 3 Cone data structure.
+/////////////////////////////////////////////////////////////////////////
+class TString;
+
+class AliITSConeData : public AliITSBaseVolParams {
+ public:
+    AliITSConeData() : AliITSBaseVolParams() // Default constructor
+       {fDz=0;fRmin0=0;fRmax0=0;fRmin1=0;fRmax1=0;}
+    virtual ~AliITSConeData(){} // Standard destructor
+    // Getters
+    Double_t DzAt(){return fDz;} // Return fDz coordiante
+    Double_t Rmin0() {return fRmin0;} // Returm arrau of rmin0 values
+    Double_t Rmax0() {return fRmax0;} // Returm arrau of rmax0 values
+    Double_t Rmin1() {return fRmin1;} // Returm arrau of rmin1 values
+    Double_t Rmax1() {return fRmax1;} // Returm arrau of rmax1 values
+    Double_t& Z(){return fDz;}// Returns address of fDz
+    Double_t& Rn0(){return fRmin0;}// Returns address of fRmin0
+    Double_t& Rx0(){return fRmax0;}// Returns address of fRmax0
+    Double_t& Rn1(){return fRmin1;}// Returns address of fRmin1
+    Double_t& Rx1(){return fRmax1;}// Returns address of fRmax1
+    void Print(ostream *os); // Prints output content of this class
+    void Read(istream *is); // Reads output created by Print above.
+ private:
+    Double_t fRmin0,fRmin1; // Inner Radius
+    Double_t fRmax0,fRmax1; // Outer Radius
+    Double_t fDz;   // Z half length.
+
+    ClassDef(AliITSConeData,1) // Simple Cone data class
+};
+// Input and output function for standard C++ input/output.
+ostream &operator<<(ostream &os,AliITSConeData &source);
+istream &operator>>(istream &os,AliITSConeData &source);
+#endif
+  
+#ifndef ALIITSCONESEGDATA_H
+#define ALIITSCONESEGDATA_H
+/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ * See cxx source for full Copyright notice                               */
+
+/*
+  $Id$
+ */
+
+/////////////////////////////////////////////////////////////////////////
+//  Geant 3 Cone Segment data structure.
+/////////////////////////////////////////////////////////////////////////
+class TString;
+
+class AliITSConeSegData : public AliITSBaseVolParams {
+ public:
+    AliITSConeSegData() : AliITSBaseVolParams() // Default constructor
+       {fPhi0=0.0;fPhi1=360.0;fDz=0;fRmin0=0;fRmax0=0;fRmin1=0;fRmax1=0;}
+    virtual ~AliITSConeSegData(){} // Standard destructor
+    // Getters
+    Double_t Phi0() {return fPhi0;} // Returns starting Phi value
+    Double_t Phi1() {return fPhi1;} // Returns endinging Phi value
+    Double_t DzAt() {return fDz;} // Return fDz coordiante
+    Double_t Rmin0() {return fRmin0;} // Returm arrau of rmin0 values
+    Double_t Rmax0() {return fRmax0;} // Returm arrau of rmax0 values
+    Double_t Rmin1() {return fRmin1;} // Returm arrau of rmin1 values
+    Double_t Rmax1() {return fRmax1;} // Returm arrau of rmax1 values
+    Double_t& P0(){return fPhi0;}// Returns addres of fPhi0
+    Double_t& P1(){return fPhi1;}// Returns addres of fPhi1
+    Double_t& Z(){return fDz;}// Returns address of fDz
+    Double_t& Rn0(){return fRmin0;}// Returns address of fRmin0
+    Double_t& Rx0(){return fRmax0;}// Returns address of fRmax0
+    Double_t& Rn1(){return fRmin1;}// Returns address of fRmin1
+    Double_t& Rx1(){return fRmax1;}// Returns address of fRmax1
+    void Print(ostream *os); // Prints output content of this class
+    void Read(istream *is); // Reads output created by Print above.
+ private:
+    Double_t fPhi0,fPhi1;   // Starting and ending phi angles [degrees]
+    Double_t fRmin0,fRmin1; // Inner Radius
+    Double_t fRmax0,fRmax1; // Outer Radius
+    Double_t fDz;   // Z half length.
+
+    ClassDef(AliITSConeSegData,1) // Cone segment data class
+};
+// Input and output function for standard C++ input/output.
+ostream &operator<<(ostream &os,AliITSConeSegData &source);
+istream &operator>>(istream &os,AliITSConeSegData &source);
+#endif
+
 #ifndef ALIITSPCONEDATA_H
 #define ALIITSPCONEDATA_H
 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
@@ -166,9 +696,107 @@ class AliITSPConeData : public AliITSBaseVolParams {
 ostream &operator<<(ostream &os,AliITSPConeData &source);
 istream &operator>>(istream &os,AliITSPConeData &source);
 
+#endif
+
+#ifndef ALIITSSPHEREDATA_H
+#define ALIITSSPHEREDATA_H
+/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ * See cxx source for full Copyright notice                               */
+
+/*
+  $Id$
+ */
+
+/////////////////////////////////////////////////////////////////////////
+//  Geant 3 Tube Segment data structure.
+/////////////////////////////////////////////////////////////////////////
+class TString;
+
+class AliITSSphereData : public AliITSBaseVolParams {
+ public:
+    AliITSSphereData() : AliITSBaseVolParams() // Default constructor
+       {fTheta[0]=0.0,fTheta[1]=0.0;fPhi[0]=0.0,fPhi[1]=0.0;fRmin=0;fRmax=0;}
+    virtual ~AliITSSphereData(){} // Standard destructor
+    // Getters
+    Double_t Theta0() {return fTheta[0];} // Returns starting Phi value
+    Double_t Theta1() {return fTheta[1];} // Returns endinging Phi value
+    Double_t Phi0() {return fPhi[0];} // Returns starting Phi value
+    Double_t Phi1() {return fPhi[1];} // Returns endinging Phi value
+    Double_t Rmin() {return fRmin;} // Returm arrau of rmin values
+    Double_t Rmax() {return fRmax;} // Returm arrau of rmax values
+    Double_t& T0(){return fTheta[0];}// Returns addres of fTheta[0]
+    Double_t& T1(){return fTheta[1];}// Returns addres of fTheta[1]
+    Double_t& P0(){return fPhi[0];}// Returns addres of fPhi[0]
+    Double_t& P1(){return fPhi[1];}// Returns addres of fPhi[1]
+    Double_t& Rn(){return fRmin;}// Returns address of fRmin
+    Double_t& Rx(){return fRmax;}// Returns address of fRmax
+    void Print(ostream *os); // Prints output content of this class
+    void Read(istream *is); // Reads output created by Print above.
+ private:
+    Double_t fTheta[2]; // Starting and ending theta angles [degree]
+    Double_t fPhi[2];   // Starting and ending phi angles [degrees]
+    Double_t fRmin; // Inner Radius
+    Double_t fRmax; // Outer Radius
 
+    ClassDef(AliITSSphereData,1) // sphere data class
+};
+// Input and output function for standard C++ input/output.
+ostream &operator<<(ostream &os,AliITSSphereData &source);
+istream &operator>>(istream &os,AliITSSphereData &source);
 #endif
-   
+
+#ifndef ALIITSPARALLELEPIPEDDATA_H
+#define ALIITSPARALLELEPIPEDDATA_H
+/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ * See cxx source for full Copyright notice                               */
+
+/*
+  $Id$
+ */
+
+/////////////////////////////////////////////////////////////////////////
+//  Geant 3 Box data structure.
+/////////////////////////////////////////////////////////////////////////
+class TString;
+
+class AliITSParallelpipedData : public AliITSBaseVolParams {
+ public:
+    AliITSParallelpipedData() : AliITSBaseVolParams() // Default constructor
+       {fDx=0.0;fDy=0.0;fDz=0;}
+    virtual ~AliITSParallelpipedData(){} // Standard destructor
+    // Getters
+    Double_t DxAt() {return fDx;} // Returm arrau of rmin values
+    Double_t DyAt() {return fDy;} // Returm arrau of rmax values
+    Double_t DzAt() {return fDz;} // Return fDz coordiante
+    Double_t Theta() {return fTheta;} // Returm arrau of Theta values
+    Double_t Phi() {return fPhi;} // Return Phi coordiante
+    Double_t Alpha(){return fAlpha;} // Return fAlph
+    Double_t& Dx(){return fDx;}// Returns address of fRmax
+    Double_t& Dy(){return fDy;}// Returns address of fRmin
+    Double_t& Dz(){return fDz;}// Returns address of fDz
+    Double_t& Th(){return fTheta;}// Returns address of fTheta
+    Double_t& Ph(){return fPhi;}// Returns address of fPhi
+    Double_t& A(){return fAlpha;}// Returns address of fAlpha
+    void Print(ostream *os); // Prints output content of this class
+    void Read(istream *is); // Reads output created by Print above.
+ private:
+    Double_t fDx;    // X half length
+    Double_t fDy;    // Y half length
+    Double_t fDz;    // Z half length.
+    Double_t fAlpha; // angle formed by the y axis and by the plane 
+                     // joining the center of teh faces parallel to the 
+                     // z-x plane at -dY and +dy [degree].
+    Double_t fTheta; //polar angle of the line joining the centers of 
+                     // the faces at -dz and +dz in z [degree].
+    Double_t fPhi;   // azimuthal angle of teh line joing the centers 
+                     // of the faaces at -dz and +dz in z [degree].
+
+    ClassDef(AliITSParallelpipedData,1) // Parallel piped data class
+};
+// Input and output function for standard C++ input/output.
+ostream &operator<<(ostream &os,AliITSParallelpipedData &source);
+istream &operator>>(istream &os,AliITSParallelpipedData &source);
+#endif   
 
 #ifndef ALIITSPGONDATA_H
 #define ALIITSPGONDATA_H
@@ -288,8 +916,8 @@ class AliITSBaseGeometry : public TObject {
     virtual void CreateG3Materials(){};//Calls Geant3 interface for materials
     virtual Int_t IsVersion() const{return 11;}// return version of geometry.
     // Get Index for Geant3 v name
-    Int_t ITSG3VnameToIndex(const char name[3])const; 
-    Int_t ITSIndexToITSG3name(const Int_t i); // Get Geant3 volume name
+    static Int_t ITSG3VnameToIndex(const char name[3]);
+    static Int_t ITSIndexToITSG3name(const Int_t i); // Get Geant3 volume name
     Int_t AddVolName(const TString name); // Add volumen name to list
     TString GetVolName(const Int_t i)const; // Return volume name at index
     Int_t GetVolumeIndex(const TString &a);
@@ -303,86 +931,95 @@ class AliITSBaseGeometry : public TObject {
     // Create a Box
     virtual void Box(const char *gnam,const TString &dis,
              Double_t dx,Double_t dy,Double_t dz,Int_t med);
+    virtual void Box(AliITSBoxData &d,Int_t med);
     // Greate A Trapizoid with the x dimension varing along z.
     virtual void Trapezoid1(const char *gnam,const TString &dis,Double_t dxn,
                     Double_t dxp,Double_t dy,Double_t dz,Int_t med);
+    virtual void Trapezoid1(AliITSTrapezoid1Data &d,Int_t med);
     // Greate A Trapizoid with the x and y dimension varing along z.
     virtual void Trapezoid2(const char *gnam,const TString &dis,Double_t dxn,
                     Double_t dxp,Double_t dyn,Double_t dyp,Double_t dz,
                     Int_t med);
+    virtual void Trapezoid2(AliITSTrapezoid2Data &d,Int_t med);
     // General trapazoid.
     virtual void Trapezoid(const char *gnam,const TString &dis,Double_t dz,
                    Double_t thet,Double_t phi,Double_t h1,Double_t bl1,
                    Double_t tl1,Double_t alp1,Double_t h2,Double_t bl2,
                    Double_t tl2,Double_t alp2,Int_t med);
+    virtual void Trapezoid(AliITSTrapezoidData &d,Int_t med);
+    // Twisted genral trapezoid.
+    virtual void TwistedTrapezoid(const char *gnam,const TString &dis,
+                                 Double_t dz,
+                                 Double_t thet,Double_t phi,Double_t twist,
+                                 Double_t h1,Double_t bl1,Double_t tl1,
+                                 Double_t apl1,Double_t h2,Double_t bl2,
+                                 Double_t tl2,Double_t apl2,Int_t med);
+    virtual void TwistedTrapezoid(AliITSTrapezoidTwistedData &d,Int_t med);
     // Simple Tube.
     virtual void Tube(const char *gnam,const TString &dis,Double_t rmin,
               Double_t rmax,Double_t dz,Int_t med);
     virtual void Tube(AliITSTubeData &d,Int_t med);
     // Tube segment.
-    virtual void TubeSegment(const char *gnam,const TString &dis,Double_t rmin,
-                     Double_t rmax,Double_t dz,Double_t phi1,Double_t phi2,
-                     Int_t med);
+    virtual void TubeSegment(const char *gnam,const TString &dis,
+                            Double_t rmin,Double_t rmax,Double_t dz,
+                            Double_t phi1,Double_t phi2,Int_t med);
+    virtual void TubeSegment(AliITSTubeSegData &v,Int_t med);
+    // Cut tube.
+    virtual void CutTube(const char *gnam,const TString &dis,Double_t rmin,
+                 Double_t rmax,Double_t dz,Double_t phi1,Double_t phi2,
+                 Double_t lx,Double_t ly,Double_t lz,Double_t hx,Double_t hy,
+                 Double_t hz,Int_t med);
+    virtual void CutTube(AliITSTubeCutData &d,Int_t med);
+    // Ellliptical cross-sectino tube
+    virtual void TubeElliptical(const char *gnam,const TString &dis,
+                               Double_t p1,Double_t p2,Double_t dz,
+                               Int_t med);
+    virtual void TubeElliptical(AliITSTubeEllipticalData &v,Int_t med);
+    // Hyperbolic tube
+    virtual void HyperbolicTube(const char *gnam,const TString &dis,
+                               Double_t rmin,Double_t rmax,Double_t dz,
+                               Double_t thet,Int_t med);
+    virtual void HyperbolicTube(AliITSTubeHyperbolicData &d,Int_t med);
     // Simple Cone.
     virtual void Cone(const char *gnam,const TString &dis,Double_t dz,
              Double_t rmin1,Double_t rmax1,Double_t rmin2,Double_t rmax2,
              Int_t med);
+    virtual void Cone(AliITSConeData &d,Int_t med);
     // Segment of a Cone.
     virtual void ConeSegment(const char *gnam,const TString &dis,Double_t dz,
                      Double_t rmin1,Double_t rmax1,Double_t rmin2,
                      Double_t rmax2,Double_t phi1,Double_t phi2,Int_t med);
+    virtual void ConeSegment(AliITSConeSegData &d,Int_t med);
+    //Poly-Cone
+    virtual void PolyCone(const char *gnam,const TString &dis,Double_t phi1,
+                  Double_t dphi,Int_t nz,Double_t *z,Double_t *rmin,
+                  Double_t *rmax,Int_t med);
+    virtual void PolyCone(AliITSPConeData &d,Int_t med);
     // Spherical shell segment.
     virtual void Sphere(const char *gnam,const TString &dis,Double_t rmin,
                 Double_t rmax,Double_t the1,Double_t the2,Double_t phi1,
                 Double_t phi2,Int_t med);
+    virtual void Sphere(AliITSSphereData &d,Int_t med);
     // Parallelepiped.
     virtual void Parallelepiped(const char *gnam,const TString &dis,
                                Double_t dx,Double_t dy,Double_t dz,
                                Double_t alph,Double_t thet,
                                Double_t phi,Int_t med);
+    virtual void Parallelepiped(AliITSParallelpipedData &d,Int_t med);
     // Polygon.
     virtual void PolyGon(const char *gnam,const TString &dis,Double_t phi1,
                  Double_t dphi,Int_t npdv,Int_t nz,Double_t *z,Double_t *rmin,
                  Double_t *rmax,Int_t med);
-    virtual void PolyGon(const char *gnam,const TString &dis,AliITSPGonData &d,
-                 Int_t med){PolyGon(gnam,dis,d.Phi0(),d.DPhi(),d.NPhi(),
-                                     d.Nz(),d.Z(),d.Rmin(),d.Rmax(),med);}
     virtual void PolyGon(AliITSPGonData &d,Int_t med);
-    //Poly-Cone
-    virtual void PolyCone(const char *gnam,const TString &dis,Double_t phi1,
-                  Double_t dphi,Int_t nz,Double_t *z,Double_t *rmin,
-                  Double_t *rmax,Int_t med);
-    virtual void PolyCone(const char *gnam,const TString &dis,
-                         AliITSPConeData &d,
-                         Int_t med){PolyCone(gnam,dis,d.Phi0(),
-                                             d.DPhi(),d.Nz(),d.Z(),
-                                             d.Rmin(),d.Rmax(),med);}
-    virtual void PolyCone(AliITSPConeData &d,Int_t med);
-    // Ellliptical cross-sectino tube
-    virtual void TubeElliptical(const char *gnam,const TString &dis,
-                               Double_t p1,Double_t p2,Double_t dz,Int_t med);
-    // Hyperbolic tube
-    virtual void HyperbolicTube(const char *gnam,const TString &dis,
-                               Double_t rmin,Double_t rmax,Double_t dz,
-                               Double_t thet,Int_t med);
-    // Twisted genral trapezoid.
-    virtual void TwistedTrapezoid(const char *gnam,const TString &dis,
-                                 Double_t dz,
-                                 Double_t thet,Double_t phi,Double_t twist,
-                                 Double_t h1,Double_t bl1,Double_t tl1,
-                                 Double_t apl1,Double_t h2,Double_t bl2,
-                                 Double_t tl2,Double_t apl2,Int_t med);
-    // Cut tube.
-    virtual void CutTube(const char *gnam,const TString &dis,Double_t rmin,
-                 Double_t rmax,Double_t dz,Double_t phi1,Double_t phi2,
-                 Double_t lx,Double_t ly,Double_t lz,Double_t hx,Double_t hy,
-                 Double_t hz,Int_t med);
     // Position one volume inside another
     virtual void Pos(const char *vol,Int_t cn,const char *moth,Double_t x,
              Double_t y,Double_t z,Int_t irot);
     // Position one volume inside another
     virtual void Pos(AliITSBaseVolParams &v,Int_t cn,
             AliITSBaseVolParams &m,TVector3 &t,Int_t irot);
+    // Position one volume inside another
+    virtual void Pos(AliITSBaseVolParams &v,AliITSBaseVolParams &m,
+                    TVector3 &t,Int_t irot){Pos(v,m.GetG3cpn(),m,t,irot);};
     Int_t GetMed(Int_t med){return (fits->GetIdtmed())->At(med);}
     // Define rotation matrix
     void Matrix(Int_t irot,Double_t thet1,Double_t phi1,Double_t thet2,
@@ -426,6 +1063,9 @@ class AliITSBaseGeometry : public TObject {
     Double_t ATand(Double_t t){return 180.0*TMath::ATan(t)/TMath::Pi();}
     Double_t ATand2(Double_t y,Double_t x){return 180.0*TMath::ATan2(y,x)/
                                               TMath::Pi();}
+    // gives angles in degree between 0.<=t<360.
+    Double_t Mod360(Double_t t){if(t>=360.) return Mod360(t-360.);
+        if(t<0.0) return Mod360(t+360.);return t;}
     Double_t RadLength(Int_t iz,Double_t a); // Computes radiation length
                                              // for an element
  private:
index 969fa0a..8bca77c 100644 (file)
@@ -15,6 +15,9 @@
 
 /*
 $Log$
+Revision 1.3  2003/04/17 22:29:23  nilsen
+Geometry bug fixes and the like. Work still progressing.
+
 Revision 1.2  2003/03/25 23:27:19  nilsen
 ITS new Geometry files. Not yet ready for uses, committed to allow additional
 development.
@@ -440,7 +443,7 @@ AliITSGeometrySSDCone::AliITSGeometrySSDCone(AliITS *its,TVector3 &tran,
     fS.Rx() = fR.Rmax();
     // SDD support plate, SSD side.
     fT.dP()  = 180.0*fWsddSupportPlate/(fRsddSupportPlate*TMath::Pi());
-    fT.P0()  = fPhi0SDDsupports=90.0;
+    fT.P0()  = fPhi0SDDsupports=-0.5*fT.DPhi();
     fT.Z(0)  = fK.ZAt(2);
     fT.Rn(0) = fI.Rmin(4);
     fT.Rx(0) = fRsddSupportPlate;
@@ -622,6 +625,48 @@ void AliITSGeometrySSDCone::PositionG3Geometry(AliITSBaseVolParams &moth,
        Pos(fN,i+1,fI,zero,irotPost+i);
        Pos(fM,i+1,fJ,zero,irotPost+i);
     } // end for i
+    // Add SDD thermal sheald screws, and SDD suppport bracket.
+    t = fPhi0Screws;
+    v.SetX(fRcylinderScrews*Sind(t));
+    v.SetY(fRcylinderScrews*Cosd(t));
+    v.SetZ(0.0);
+    Pos(fQ,1,fK,v,0);
+    v.SetZ(fK.ZAt(2)-fR.DzAt());
+    Pos(fR,1,fK,v,0);
+    v.SetZ(fI.ZAt(4)-fS.DzAt());
+    Pos(fS,1,fI,v,0);
+    Pos(fT,1,fK,zero,0);
+    Pos(fU,1,fI,zero,0);
+    Int_t ir=2,it=2;
+    Double_t d = 360./((Double_t)fNssdSupports);
+    cout << fPhi0SDDsupports<<" " << fNssdSupports<< " " << fPhi0Screws<< endl;
+    for(i=1;i<fNinScrews;i++){
+       t = fPhi0Screws + ((Double_t)i)*360./((Double_t)fNinScrews);
+       v.SetX(fRcylinderScrews*Sind(t));
+       v.SetY(fRcylinderScrews*Cosd(t));
+       v.SetZ(0.0);
+       Pos(fQ,i+1,fK,v,0);
+       cout << " t="<<t<<endl;
+       if((t>=-0.5*fT.DPhi()&&t<=+0.5*fT.DPhi())||
+          (t>=d-0.5*fT.DPhi()&&t<=d+0.5*fT.DPhi())||
+          (t>=2.*d-0.5*fT.DPhi()&&t<=2.*d+0.5*fT.DPhi())){
+           // SDD support here
+           if((t>=-0.5*fT.DPhi()&&t<=+0.5*fT.DPhi())) continue;
+           t = ((Double_t)(it-1))*d;
+           cout << " Zmatrix t=" << t << endl;
+           ZMatrix(++irot,t);
+           Pos(fT,it,fK,zero,irot);
+           Pos(fU,it,fI,zero,irot);
+           it++;
+       }else{ // Air here
+           v.SetZ(fK.ZAt(2)-fR.DzAt());
+           Pos(fR,ir,fK,v,0);
+           v.SetZ(fI.ZAt(4)-fS.DzAt());
+           Pos(fS,ir,fI,v,0);
+           ir++;
+       } // end if
+    } // end for i
+    // Air holes and bolts in SSD mounting studs, volumes I,M,N,O.
     return;
 }
 //______________________________________________________________________
index 9246a53..c7f54aa 100644 (file)
 //
 #pragma link C++ class  AliITSBaseGeometry+;
 #pragma link C++ class  AliITSBaseVolParams+;
+#pragma link C++ class  AliITSBoxData+;
+#pragma link C++ class  AliITSTrapezoid1Data+;
+#pragma link C++ class  AliITSTrapezoid2Data+;
+#pragma link C++ class  AliITSTrapezoidData+;
+#pragma link C++ class  AliITSTrapezoidTwistedData+;
 #pragma link C++ class  AliITSTubeData+;
+#pragma link C++ class  AliITSTubeSegData+;
+#pragma link C++ class  AliITSTubeCutData+;
+#pragma link C++ class  AliITSTubeEllipticalData+;
+#pragma link C++ class  AliITSTubeHyperbolicData+;
+#pragma link C++ class  AliITSConeData+;
+#pragma link C++ class  AliITSConeSegData+;
 #pragma link C++ class  AliITSPConeData+;
+#pragma link C++ class  AliITSSphereData+;
+#pragma link C++ class  AliITSParallelpipedData+;
 #pragma link C++ class  AliITSPGonData+;
 #pragma link C++ class  AliITSGeometryITSV+;
 #pragma link C++ class  AliITSGeometrySSDCone+;