- fixing warnings/coverity
[u/mrichter/AliRoot.git] / ITS / AliITSv11GeometrySPD.cxx
index 79e6cd9..1f37c6f 100644 (file)
@@ -14,7 +14,7 @@
  **************************************************************************/
 //
 // This class Defines the Geometry for the ITS services and support cones
-// outside of the central volume (except for the Central support 
+// outside of the central volume (except for the Central support
 // cylinders). Other classes define the rest of the ITS, specifically the
 // SSD support cone, the SSD Support central cylinder, the SDD support cone,
 // the SDD support central cylinder, the SPD Thermal Shield, The supports
 #include <TPolyMarker.h>
 
 // Root Geometry includes
-#include <TGeoVolume.h>
-#include <TGeoTube.h> // contains TGeoTubeSeg
+#include <TGeoCompositeShape.h>
 #include <TGeoEltu.h>
-#include <TGeoXtru.h>
-#include <TGeoMatrix.h>
+#include <TGeoGlobalMagField.h>
 #include <TGeoMaterial.h>
+#include <TGeoMatrix.h>
 #include <TGeoMedium.h>
-#include <TGeoCompositeShape.h>
+#include <TGeoTube.h> // contains TGeoTubeSeg
+#include <TGeoVolume.h>
+#include <TGeoXtru.h>
+#include <TGeoPcon.h>
 
 // AliRoot includes
 #include "AliLog.h"
@@ -87,9 +89,9 @@
 #include "AliITSv11GeometrySPD.h"
 
 // Constant definistions
-const Double_t AliITSv11GeometrySPD::fgkGapLadder    = 
+const Double_t AliITSv11GeometrySPD::fgkGapLadder    =
                       AliITSv11Geometry::fgkmicron*75.; //  75 microns
-const Double_t AliITSv11GeometrySPD::fgkGapHalfStave = 
+const Double_t AliITSv11GeometrySPD::fgkGapHalfStave =
                      AliITSv11Geometry::fgkmicron*120.; // 120 microns
 
 ClassImp(AliITSv11GeometrySPD)
@@ -102,11 +104,11 @@ fSPDsectorX0(0),    // X of first edge of sector plane for stave
 fSPDsectorY0(0),    // Y of first edge of sector plane for stave
 fSPDsectorX1(0),    // X of second edge of sector plane for stave
 fSPDsectorY1(0),    // Y of second edge of sector plane for stave
-fTubeEndSector()    // coordinate of cooling tube ends 
+fTubeEndSector()    // coordinate of cooling tube ends
 {
     //
     // Default constructor.
-    // This does not initialize anything and is provided just for 
+    // This does not initialize anything and is provided just for
     // completeness. It is recommended to use the other one.
     // The alignment gap is specified as argument (default = 0.0075 cm).
     // Inputs:
@@ -133,7 +135,7 @@ fSPDsectorX0(0),    // X of first edge of sector plane for stave
 fSPDsectorY0(0),    // Y of first edge of sector plane for stave
 fSPDsectorX1(0),    // X of second edge of sector plane for stave
 fSPDsectorY1(0),    // Y of second edge of sector plane for stave
-fTubeEndSector()    // coordinate of cooling tube ends 
+fTubeEndSector()    // coordinate of cooling tube ends
 {
     //
     // Constructor with debug setting argument
@@ -183,7 +185,7 @@ fSPDsectorY1(s.fSPDsectorY1)     // Y of second edge of sector plane for stave
     } // end for i,j
 }
 //______________________________________________________________________
-AliITSv11GeometrySPD& AliITSv11GeometrySPD::operator=(const 
+AliITSv11GeometrySPD& AliITSv11GeometrySPD::operator=(const
                                                AliITSv11GeometrySPD &s)
 {
     //
@@ -214,265 +216,30 @@ TGeoMedium* AliITSv11GeometrySPD::GetMedium(const char* mediumName,
                                             TGeoManager *mgr) const
 {
     //
-    // This function is used to recovery any medium 
-    // used to build the geometry volumes. 
-    // If the required medium does not exists, 
+    // This function is used to recovery any medium
+    // used to build the geometry volumes.
+    // If the required medium does not exists,
     // a NULL pointer is returned, and an error message is written.
     //
      Char_t itsMediumName[30];
 
-     sprintf(itsMediumName, "ITS_%s", mediumName);
+     snprintf(itsMediumName, 30, "ITS_%s", mediumName);
      TGeoMedium* medium = mgr->GetMedium(itsMediumName);
      if (!medium) AliError(Form("Medium <%s> not found", mediumName));
 
      return medium;
 }
-//______________________________________________________________________
-Int_t AliITSv11GeometrySPD::CreateSPDCentralMaterials(Int_t &medOffset,
-                                                      Int_t &matOffset) const
-{
-    //
-    // Define the specific materials used for the ITS SPD central detectors.
-    // ---
-    // NOTE: These are the same old names. 
-    //       By the ALICE naming conventions, they start with "ITS SPD ...."
-    //       Data taken from ** AliITSvPPRasymmFMD::CreateMaterials() **.
-    // ---
-    // Arguments [the ones passed by reference contain output values]:
-    // - medOffset --> (by ref) starting number of the list of media
-    // - matOffset --> (by ref) starting number of the list of Materials
-    // ---
-    // Inputs:
-    //   Int_t &medOffset  Starting number of the list of media
-    //   Int_t &matOffset  Starting number of the list of materials
-    // Outputs:
-    //   Int_t &medOffset  Ending number of the list of media
-    //   Int_t &matOffset  Ending number of the list of materials
-    // Return:
-    //   The last material indexused +1. (= next avaiable material index)
-    //
-    const Double_t ktmaxfd    = 0.1 * fgkDegree; // Degree
-    const Double_t kstemax    = 1.0 * fgkcm; // cm
-    const Double_t kdeemax    = 0.1;//Fraction of particle's energy 0<deemax<=1
-    const Double_t kepsil     = 1.0E-4; //
-    const Double_t kstmin     = 0.0 * fgkcm; // cm "Default value used"
-    const Double_t ktmaxfdAir = 0.1 * fgkDegree; // Degree
-    const Double_t kstemaxAir = 1.0000E+00 * fgkcm; // cm
-    const Double_t kdeemaxAir = 0.1;//Fraction of particle's energy 0<deemax<=1
-    const Double_t kepsilAir  = 1.0E-4;//
-    const Double_t kstminAir  = 0.0 * fgkcm; // cm "Default value used"
-    const Double_t ktmaxfdSi  = 0.1 * fgkDegree; // .10000E+01; // Degree
-    const Double_t kstemaxSi  = 0.0075 * fgkcm; //  .10000E+01; // cm
-    const Double_t kdeemaxSi  = 0.1;//Fraction of particle's energy 0<deemax<=1
-    const Double_t kepsilSi   = 1.0E-4;//
-    const Double_t kstminSi   = 0.0 * fgkcm; // cm "Default value used"
-    //
-    Int_t matindex = matOffset;
-    Int_t medindex = medOffset;
-    TGeoMaterial *mat;
-    TGeoMixture  *mix;
-    TGeoMedium   *med;
-    //
-    Int_t    ifield = (gAlice->Field()->Integ());
-    Double_t fieldm = (gAlice->Field()->Max());
-    Double_t params[8] = {8 * 0.0};
-
-    params[1] = (Double_t) ifield;
-    params[2] = fieldm;
-    params[3] = ktmaxfdSi;
-    params[4] = kstemaxSi;
-    params[5] = kdeemaxSi;
-    params[6] = kepsilSi;
-    params[7] = kstminSi;
-
-    // Definition of materials and mediums.
-    // Last argument in material definition is its pressure,
-    // which is initialized to ZERO.
-    // For better readability, it is simply set to zero.
-    // Then the writing "0.0 * fgkPascal" is replaced by "0."
-    // (Alberto)
-    
-    // silicon definition for ITS (overall)
-    mat = new TGeoMaterial("ITS_SI", 28.086, 14.0, 2.33 * fgkgcm3,
-                           TGeoMaterial::kMatStateSolid, 25.0*fgkCelsius, 0.);
-    mat->SetIndex(matindex);
-    med = new TGeoMedium("SI", medindex++, mat, params);
-    
-    // silicon for ladder chips
-    mat = new TGeoMaterial("SPD SI CHIP", 28.086, 14.0, 2.33 * fgkgcm3,
-                           TGeoMaterial::kMatStateSolid, 25.0*fgkCelsius, 0.);
-    mat->SetIndex(matindex);
-    med = new TGeoMedium("SPD SI CHIP", medindex++, mat, params);
-    
-    // silicon for pixel bus
-    mat = new TGeoMaterial("SPD SI BUS", 28.086, 14.0, 2.33 * fgkgcm3,
-                           TGeoMaterial::kMatStateSolid, 25.0*fgkCelsius, 0.);
-    mat->SetIndex(matindex);
-    med = new TGeoMedium("SPD SI BUS", medindex++, mat, params);
-    
-    // carbon fiber material is defined as a mix of C-O-N-H
-    // defined in terms of fractional weights according to 'C (M55J)'
-    // it is used for the support and clips
-    mix = new TGeoMixture("C (M55J)", 4, 1.9866 * fgkgcm3);
-    mix->SetIndex(matindex);
-    mix->DefineElement(0, 12.01070, 6.0, 0.908508078);// C by fractional weight
-    mix->DefineElement(1, 14.00670, 7.0, 0.010387573);// N by fractional weight
-    mix->DefineElement(2, 15.99940, 8.0, 0.055957585);// O by fractional weight
-    mix->DefineElement(3,  1.00794, 1.0, 0.025146765);// H by fractional weight
-    mix->SetPressure(0.0 * fgkPascal);
-    mix->SetTemperature(25.0 * fgkCelsius);
-    mix->SetState(TGeoMaterial::kMatStateSolid);
-    params[3] = ktmaxfd;
-    params[4] = kstemax;
-    params[5] = kdeemax;
-    params[6] = kepsil;
-    params[7] = kstmin;
-    med = new TGeoMedium("ITSspdCarbonFiber", medindex++, mix, params);
-
-    // air defined as a mixture of C-N-O-Ar: 
-    // it is used to fill all containers
-    mix = new TGeoMixture("Air", 4, 1.20479E-3 * fgkgcm3);
-    mix->SetIndex(matindex);
-    mix->DefineElement(0, 12.0107,  6.0, 0.000124); // C by fractional weight
-    mix->DefineElement(1, 14.0067,  7.0, 0.755267); // N by fractional weight
-    mix->DefineElement(2, 15.9994,  8.0, 0.231781); // O by fractional weight
-    mix->DefineElement(3, 39.9480, 18.0, 0.012827); // Ar by fractional weight
-    mix->SetPressure(101325.0 * fgkPascal); // = 1 atmosphere
-    mix->SetTemperature(25.0 * fgkCelsius);
-    mix->SetState(TGeoMaterial::kMatStateGas);
-    params[3] = ktmaxfdAir;
-    params[4] = kstemaxAir;
-    params[5] = kdeemaxAir;
-    params[6] = kepsilAir;
-    params[7] = kstminAir;
-    med = new TGeoMedium("ITSspdAir", medindex++, mix, params);
-
-    // inox stainless steel, defined as a mixture
-    // used for all metallic parts
-    mix = new TGeoMixture("INOX", 9, 8.03 * fgkgcm3);
-    mix->SetIndex(matindex);
-    mix->DefineElement(0, 12.0107,  6., .0003);  // C  by fractional weight
-    mix->DefineElement(1, 54.9380, 25., .02);    // Fe by fractional weight
-    mix->DefineElement(2, 28.0855, 14., .01);    // Na by fractional weight
-    mix->DefineElement(3, 30.9738, 15., .00045); // P  by fractional weight
-    mix->DefineElement(4, 32.066 , 16., .0003);  // S  by fractional weight
-    mix->DefineElement(5, 58.6928, 28., .12);    // Ni by fractional weight
-    mix->DefineElement(6, 55.9961, 24., .17);    //    by fractional weight
-    mix->DefineElement(7, 95.84  , 42., .025);   //    by fractional weight
-    mix->DefineElement(8, 55.845 , 26., .654);   //    by fractional weight
-    mix->SetPressure(0.0 * fgkPascal);
-    mix->SetTemperature(25.0 * fgkCelsius);
-    mix->SetState(TGeoMaterial::kMatStateSolid);
-    params[3] = ktmaxfdAir;
-    params[4] = kstemaxAir;
-    params[5] = kdeemaxAir;
-    params[6] = kepsilAir;
-    params[7] = kstminAir;
-    med = new TGeoMedium("ITSspdStainlessSteel", medindex++, mix, params);
-
-    // freon gas which fills the cooling system (C+F)
-    mix = new TGeoMixture("Freon", 2, 1.63 * fgkgcm3);
-    mix->SetIndex(matindex);
-    mix->DefineElement(0, 12.0107   , 6.0,  4);  // C by fractional weight
-    mix->DefineElement(1, 18.9984032, 9.0, 10); // F by fractional weight
-    mix->SetPressure(101325.0 * fgkPascal); // = 1 atmosphere
-    mix->SetTemperature(25.0 * fgkCelsius);
-    mix->SetState(TGeoMaterial::kMatStateLiquid);
-    params[3] = ktmaxfdAir;
-    params[4] = kstemaxAir;
-    params[5] = kdeemaxAir;
-    params[6] = kepsilAir;
-    params[7] = kstminAir;
-    med = new TGeoMedium("ITSspdCoolingFluid", medindex++, mix, params);
-
-    // return the next index to be used in case of adding new materials
-    medOffset = medindex;
-    matOffset = matindex;
-    return matOffset;
-}
-//______________________________________________________________________
-void AliITSv11GeometrySPD::InitSPDCentral(Int_t offset, TVirtualMC *vmc) const
-{
-     //
-     // Do all SPD Central detector initializations (e.g.: transport cuts).
-     // ---
-     // Here follow some GEANT3 physics switches, which are interesting 
-     // for these settings to be defined:
-     // - "MULTS" (MULtiple Scattering):
-     //   the variable IMULS controls this process. See [PHYS320/325/328]
-     //   0 - No multiple scattering.
-     //   1 - (DEFAULT) Multiple scattering according to Moliere theory.
-     //   2 - Same as 1. Kept for backward compatibility.
-     //   3 - Pure Gaussian scattering according to the Rossi formula.
-     // - "DRAY" (Delta RAY production)
-     //   The variable IDRAY controls this process. See [PHYS430]
-     //   0 - No delta rays production.
-     //   1 - (DEFAULT) Delta rays production with generation of.
-     //   2 - Delta rays production without generation of.
-     // - "LOSS" (continuous energy loss)
-     //   The variable ILOSS controls this process.
-     //   0 - No continuous energy loss, IDRAY is set to 0.
-     //   1 - Continuous energy loss with generation of delta rays above 
-     //       DCUTE (common/GCUTS/) and restricted Landau fluctuations 
-     //        below DCUTE.
-     //   2 - (DEFAULT) Continuous energy loss without generation of 
-     //       delta rays 
-     //       and full Landau-Vavilov-Gauss fluctuations.
-     //       In this case the variable IDRAY is forced to 0 to avoid
-     //       double counting of fluctuations.
-     //   3 - Same as 1, kept for backward compatibility.
-     //   4 - Energy loss without fluctuation.
-     //       The value obtained from the tables is used directly.
-     // ---
-     // Arguments:
-     //    Int_t offset    --> the material/medium index offset
-     //    TVirtualMC *vmc --> pointer to the virtual Monte Carlo default gMC
-     //
-
-     Int_t i, n = 4;
-     
-     for(i=0;i<n;i++) {
-          vmc->Gstpar(i+offset, "CUTGAM", 30.0 * fgkKeV);
-          vmc->Gstpar(i+offset, "CUTELE", 30.0 * fgkKeV);
-          vmc->Gstpar(i+offset, "CUTNEU", 30.0 * fgkKeV);
-          vmc->Gstpar(i+offset, "CUTHAD", 30.0 * fgkKeV);
-          vmc->Gstpar(i+offset, "CUTMUO", 30.0 * fgkKeV);
-          vmc->Gstpar(i+offset, "BCUTE",  30.0 * fgkKeV);
-          vmc->Gstpar(i+offset, "BCUTM",  30.0 * fgkKeV);
-          vmc->Gstpar(i+offset, "DCUTE",  30.0 * fgkKeV);
-          vmc->Gstpar(i+offset, "DCUTM",  30.0 * fgkKeV);
-          //vmc->Gstpar(i+offset, "PPCUTM", );
-          //vmc->Gstpar(i+offset, "PAIR", );
-          //vmc->Gstpar(i+offset, "COMPT", );
-          //vmc->Gstpar(i+offset, "PHOT", );
-          //vmc->Gstpar(i+offset, "PFIS", );
-          vmc->Gstpar(i+offset, "DRAY", 1);
-          //vmc->Gstpar(i+offset, "ANNI", );
-          //vmc->Gstpar(i+offset, "BREM", );
-          //vmc->Gstpar(i+offset, "HADR", );
-          //vmc->Gstpar(i+offset, "MUNU", );
-          //vmc->Gstpar(i+offset, "DCAY", );
-          vmc->Gstpar(i+offset, "LOSS", 1);
-          //vmc->Gstpar(i+offset, "MULS", );
-          //vmc->Gstpar(i+offset, "GHCOR1", );
-          //vmc->Gstpar(i+offset, "BIRK1", );
-          //vmc->Gstpar(i+offset, "BRIK2", );
-          //vmc->Gstpar(i+offset, "BRIK3", );
-          //vmc->Gstpar(i+offset, "LABS", );
-          //vmc->Gstpar(i+offset, "SYNC", );
-          //vmc->Gstpar(i+offset, "STRA", );
-     }
-}
+
 //______________________________________________________________________
 void AliITSv11GeometrySPD::SPDSector(TGeoVolume *moth, TGeoManager *mgr)
 {
     //
-    // Creates a single SPD carbon fiber sector and places it 
+    // Creates a single SPD carbon fiber sector and places it
     // in a container volume passed as first argument ('moth').
     // Second argument points to the TGeoManager which coordinates
     // the overall volume creation.
-    // The position of the sector is based on distance of 
-    // closest point of SPD stave to beam pipe 
+    // The position of the sector is based on distance of
+    // closest point of SPD stave to beam pipe
     // (figures all-sections-modules.ps) of 7.22mm at section A-A.
     //
 
@@ -480,7 +247,7 @@ void AliITSv11GeometrySPD::SPDSector(TGeoVolume *moth, TGeoManager *mgr)
     /*
      <img src="http://alice.pd.infn.it/latestdr/Geometric-Revision/assembly.ps"
      title="SPD     Sector    drawing   with all  cross     sections  defined">
-     <p>The    SPD  Sector    definition.    In   
+     <p>The    SPD  Sector    definition.    In
      <a   href="http://alice.pd.infn.it/latestdr/Geometric-Revision/assembly.hpgl">HPGL</a>    format.
      <img src="http://alice.pd.infn.it/latestdr/Geometric-Revision/assembly-10-modules.ps"
      titile="SPD    All  Sectors   end  view with thermal   sheald">
@@ -517,7 +284,7 @@ void AliITSv11GeometrySPD::SPDSector(TGeoVolume *moth, TGeoManager *mgr)
     // Return:
     //    none.
     // Updated values for kSPDclossesStaveAA, kBeamPipeRadius, and
-    // staveThicknessAA are taken from 
+    // staveThicknessAA are taken from
     // http://physics.mps.ohio-state.edu/~nilsen/ITSfigures/Sezione_layerAA.pdf
     //
     const Double_t kSPDclossesStaveAA   =   7.25* fgkmm;
@@ -527,19 +294,22 @@ void AliITSv11GeometrySPD::SPDSector(TGeoVolume *moth, TGeoManager *mgr)
     const Double_t kBeamPipeRadius      =   0.5 * 59.6 * fgkmm; // diam. = 59.6 mm
   //const Double_t staveThicknessAA     =   0.9 *fgkmm;         // nominal thickness
     const Double_t staveThicknessAA     =   1.02 * fgkmm;       // get from stave geometry.
-    
+
     Int_t i, j, k;
     Double_t angle, radiusSector, xAAtubeCenter0, yAAtubeCenter0;
     TGeoCombiTrans *secRot = new TGeoCombiTrans(), *comrot;
-    TGeoVolume *vCarbonFiberSector;
+    TGeoVolume *vCarbonFiberSector[10];
     TGeoMedium *medSPDcf;
 
-    // Define an assembly and fill it with the support of 
+    // Define an assembly and fill it with the support of
     // a single carbon fiber sector and staves in it
     medSPDcf = GetMedium("SPD C (M55J)$", mgr);
-    vCarbonFiberSector = new TGeoVolumeAssembly("ITSSPDCarbonFiberSectorV");
-    vCarbonFiberSector->SetMedium(medSPDcf);
-    CarbonFiberSector(vCarbonFiberSector, xAAtubeCenter0, yAAtubeCenter0, mgr);
+    for(Int_t is=0; is<10; is++)
+    {
+           vCarbonFiberSector[is] = new TGeoVolumeAssembly("ITSSPDCarbonFiberSectorV");
+           vCarbonFiberSector[is]->SetMedium(medSPDcf);
+           CarbonFiberSector(vCarbonFiberSector[is], is, xAAtubeCenter0, yAAtubeCenter0, mgr);
+    }
 
     // Compute the radial shift out of the sectors
     radiusSector = kBeamPipeRadius + kSPDclossesStaveAA + staveThicknessAA;
@@ -548,13 +318,13 @@ void AliITSv11GeometrySPD::SPDSector(TGeoVolume *moth, TGeoManager *mgr)
   //radiusSector *= radiusSector; // squaring;
   //radiusSector -= xAAtubeCenter0 * xAAtubeCenter0;
   //radiusSector  = -yAAtubeCenter0 + TMath::Sqrt(radiusSector);
-    
+
     AliDebug(1, Form("SPDSector : radiusSector=%f\n",radiusSector));
     i = 1;
     AliDebug(1, Form("i= %d x0=%f y0=%f x1=%f y1=%f\n", i,
                      fSPDsectorX0.At(i), fSPDsectorY0.At(i),
                      fSPDsectorX1.At(i),fSPDsectorY1.At(i)));
-    
+
     // add 10 single sectors, by replicating the virtual sector defined above
     // and placing at different angles
     Double_t shiftX, shiftY, tub[2][6][3];
@@ -566,7 +336,7 @@ void AliITSv11GeometrySPD::SPDSector(TGeoVolume *moth, TGeoManager *mgr)
     for(i = 0; i < kNSectorsTotal; i++) {
         shiftX = -radiusSector * TMath::Sin(angle/fgkRadian);
         shiftY =  radiusSector * TMath::Cos(angle/fgkRadian);
-        //cout << "ANGLE = " << angle << endl; 
+        //cout << "ANGLE = " << angle << endl;
         shiftX += 0.1094 * TMath::Cos((angle + 196.)/fgkRadian);
         shiftY += 0.1094 * TMath::Sin((angle + 196.)/fgkRadian);
         //shiftX -= 0.105;
@@ -576,7 +346,7 @@ void AliITSv11GeometrySPD::SPDSector(TGeoVolume *moth, TGeoManager *mgr)
         secRot->SetDx(shiftX);
         secRot->SetDy(shiftY);
         comrot  = new TGeoCombiTrans(*secRot);
-        vcenteral->AddNode(vCarbonFiberSector,i+1,comrot);
+        vcenteral->AddNode(vCarbonFiberSector[i],i+1,comrot);
         for(j=0;j<2;j++)for(k=0;k<6;k++) // Transform Tube ends for each sector
             comrot->LocalToMaster(tub[j][k],fTubeEndSector[i][j][k]);
         if(GetDebug(5)) {
@@ -589,13 +359,17 @@ void AliITSv11GeometrySPD::SPDSector(TGeoVolume *moth, TGeoManager *mgr)
     } // end for i
     if(GetDebug(3)) moth->PrintNodes();
     delete secRot;
-    
+
     CreateCones(moth);
 }
 //______________________________________________________________________
-void AliITSv11GeometrySPD::CarbonFiberSector(TGeoVolume *moth,
+void AliITSv11GeometrySPD::CarbonFiberSector(TGeoVolume *moth, Int_t sect,
      Double_t &xAAtubeCenter0, Double_t &yAAtubeCenter0, TGeoManager *mgr)
 {
+    // The method has been modified in order to build a support sector
+    // whose shape is dependent on the sector number; the aim is to get
+    // as close as possible to the shape inferred from alignment
+    // and avoid as much as possible overlaps generated by alignment.
     //
     // Define the detail SPD Carbon fiber support Sector geometry.
     // Based on the drawings:
@@ -626,7 +400,7 @@ void AliITSv11GeometrySPD::CarbonFiberSector(TGeoVolume *moth,
     // center of the #0 detector and the beam pipe.
     // Measurements are taken at cross section A-A.
     //
-     
+
     //TGeoMedium *medSPDfs      = 0;//SPD support cone inserto stesalite 4411w
     //TGeoMedium *medSPDfo      = 0;//SPD support cone foam, Rohacell 50A.
     //TGeoMedium *medSPDal      = 0;//SPD support cone SDD mounting bracket Al
@@ -649,67 +423,89 @@ void AliITSv11GeometrySPD::CarbonFiberSector(TGeoVolume *moth,
     // 1) the (ksecX, ksecY) points are plotted
     //    and circles of the specified radii are drawn around them.
     // 2) each pair of consecutive circles is connected by a line
-    //    tangent to them, in accordance with the radii being "internal" 
-    //    or "external" with respect to the closed shape which describes 
+    //    tangent to them, in accordance with the radii being "internal"
+    //    or "external" with respect to the closed shape which describes
     //    the sector itself.
-    // The resulting connected shape is the section 
+    // The resulting connected shape is the section
     // of the SPD sector surface in the transverse plane (XY).
     //
     const Double_t ksecX0   = -10.725 * fgkmm;
     const Double_t ksecY0   = -14.853 * fgkmm;
     const Double_t ksecR0   =  -0.8   * fgkmm; // external
-    const Double_t ksecX1   = -13.187 * fgkmm;
-    const Double_t ksecY1   = -19.964 * fgkmm;
-    const Double_t ksecR1   =  +0.6   * fgkmm; // internal  // (modif. by Alberto)
-    //const Double_t ksecR1   =  +0.8   * fgkmm; // internal  // (modif. by Alberto)
 
-    // const Double_t ksecDip0 = 5.9 * fgkmm;
-    //
-    //const Double_t ksecX2   =  -3.883 * fgkmm;
-    const Double_t ksecX2   =  -3.833 * fgkmm; // (corr. by Alberto)
-    const Double_t ksecY2   = -17.805 * fgkmm;
-    const Double_t ksecR2   =  +0.6  * fgkmm; // internal (guess)
-    const Double_t ksecX3   =  -3.123 * fgkmm;
-    const Double_t ksecY3   = -14.618 * fgkmm;
-    const Double_t ksecR3   =  -0.6   * fgkmm; // external
-    //const Double_t ksecDip1 = 8.035 * fgkmm;
-    //
-    const Double_t ksecX4   = +11.280 * fgkmm;
-    const Double_t ksecY4   = -14.473 * fgkmm;
-    const Double_t ksecR4   =  +0.8   * fgkmm; // internal
-    const Double_t ksecX5   = +19.544 * fgkmm;
-    const Double_t ksecY5   = +10.961 * fgkmm;
-    const Double_t ksecR5   =  +0.8   * fgkmm; // internal
-    //const Double_t ksecDip2 = 4.553 * fgkmm;
-    // 
-    const Double_t ksecX6   = +10.830 * fgkmm;
-    const Double_t ksecY6   = +16.858 * fgkmm;
-    const Double_t ksecR6   =  +0.6   * fgkmm; // internal
-    const Double_t ksecX7   = +11.581 * fgkmm;
-    const Double_t ksecY7   = +13.317 * fgkmm;
-    const Double_t ksecR7   =  -0.6   * fgkmm; // external
-    //const Double_t ksecDip3 = 6.978 * fgkmm;
-    //
-    const Double_t ksecX8   =  -0.733 * fgkmm;
-    const Double_t ksecY8   = +17.486 * fgkmm;
-    const Double_t ksecR8   =  +0.6   * fgkmm; // internal
-    const Double_t ksecX9   =  +0.562 * fgkmm;
-    //const Double_t ksecY9 = +14.486 * fgkmm; // correction by
-    const Double_t ksecY9   = +14.107 * fgkmm; // Alberto
-    const Double_t ksecR9   =  -0.6   * fgkmm; // external
-    //const Double_t ksecDip4 = 6.978 * fgkmm;
-    //
-    const Double_t ksecX10  = -12.252 * fgkmm;
-    const Double_t ksecY10  = +16.298 * fgkmm;
-    const Double_t ksecR10  =  +0.6   * fgkmm; // internal
-    const Double_t ksecX11  = -10.445 * fgkmm;
-    const Double_t ksecY11  = +13.162 * fgkmm;
-    const Double_t ksecR11  =  -0.6   * fgkmm; // external
-    //const Double_t ksecDip5 = 6.978 * fgkmm;
-    //
-    const Double_t ksecX12  = -22.276 * fgkmm;
-    const Double_t ksecY12  = +12.948 * fgkmm;
-    const Double_t ksecR12  =  +0.85  * fgkmm; // internal
+    const Double_t ksecR1   =  +0.6   * fgkmm;
+    const Double_t ksecR2   =  +0.6   * fgkmm;
+    const Double_t ksecR3   =  -0.6   * fgkmm;
+    const Double_t ksecR4   =  +0.8   * fgkmm;
+    const Double_t ksecR5   =  +0.8   * fgkmm;
+    const Double_t ksecR6   =  +0.6   * fgkmm;
+    const Double_t ksecR7   =  -0.6   * fgkmm;
+    const Double_t ksecR8   =  +0.6   * fgkmm;
+    const Double_t ksecR9   =  -0.6   * fgkmm;
+    const Double_t ksecR10   =  +0.6   * fgkmm;
+    const Double_t ksecR11   =  -0.6   * fgkmm;
+    const Double_t ksecR12   =  +0.85   * fgkmm;
+
+//    // IDEAL GEOMETRY
+//     const Double_t ksecX1[10] ={-1.3187,-1.3187,-1.3187,-1.3187,-1.3187,-1.3187,-1.3187,-1.3187,-1.3187,-1.3187};
+//     const Double_t ksecY1[10] ={-1.9964,-1.9964,-1.9964,-1.9964,-1.9964,-1.9964,-1.9964,-1.9964,-1.9964,-1.9964};
+//     const Double_t ksecX2[10] ={-0.3833,-0.3833,-0.3833,-0.3833,-0.3833,-0.3833,-0.3833,-0.3833,-0.3833,-0.3833};
+//     const Double_t ksecY2[10] ={-1.7805,-1.7805,-1.7805,-1.7805,-1.7805,-1.7805,-1.7805,-1.7805,-1.7805,-1.7805};
+//     const Double_t ksecX3[10] ={-0.3123,-0.3123,-0.3123,-0.3123,-0.3123,-0.3123,-0.3123,-0.3123,-0.3123,-0.3123};
+//     const Double_t ksecY3[10] ={-1.4618,-1.4618,-1.4618,-1.4618,-1.4618,-1.4618,-1.4618,-1.4618,-1.4618,-1.4618};
+//     const Double_t ksecX4[10] ={+1.1280,+1.1280,+1.1280,+1.1280,+1.1280,+1.1280,+1.1280,+1.1280,+1.1280,+1.1280};
+//     const Double_t ksecY4[10] ={-1.4473,-1.4473,-1.4473,-1.4473,-1.4473,-1.4473,-1.4473,-1.4473,-1.4473,-1.4473};
+//     const Double_t ksecX5[10] ={+1.9544,+1.9544,+1.9544,+1.9544,+1.9544,+1.9544,+1.9544,+1.9544,+1.9544,+1.9544};
+//     const Double_t ksecY5[10] ={+1.0961,+1.0961,+1.0961,+1.0961,+1.0961,+1.0961,+1.0961,+1.0961,+1.0961,+1.0961};
+//     const Double_t ksecX6[10] ={+1.0830,+1.0830,+1.0830,+1.0830,+1.0830,+1.0830,+1.0830,+1.0830,+1.0830,+1.0830};
+//     const Double_t ksecY6[10] ={+1.6868,+1.6868,+1.6868,+1.6868,+1.6868,+1.6868,+1.6868,+1.6868,+1.6868,+1.6868};
+//     const Double_t ksecX7[10] ={+1.1581,+1.1581,+1.1581,+1.1581,+1.1581,+1.1581,+1.1581,+1.1581,+1.1581,+1.1581};
+//     const Double_t ksecY7[10] ={+1.3317,+1.3317,+1.3317,+1.3317,+1.3317,+1.3317,+1.3317,+1.3317,+1.3317,+1.3317};
+//     const Double_t ksecX8[10] ={-0.0733,-0.0733,-0.0733,-0.0733,-0.0733,-0.0733,-0.0733,-0.0733,-0.0733,-0.0733};
+//     const Double_t ksecY8[10] ={+1.7486,+1.7486,+1.7486,+1.7486,+1.7486,+1.7486,+1.7486,+1.7486,+1.7486,+1.7486};
+//     const Double_t ksecX9[10] ={+0.0562,+0.0562,+0.0562,+0.0562,+0.0562,+0.0562,+0.0562,+0.0562,+0.0562,+0.0562};
+//     const Double_t ksecY9[10] ={+1.4107,+1.4107,+1.4107,+1.4107,+1.4107,+1.4107,+1.4107,+1.4107,+1.4107,+1.4107};
+//     const Double_t ksecX10[10]={-1.2252,-1.2252,-1.2252,-1.2252,-1.2252,-1.2252,-1.2252,-1.2252,-1.2252,-1.2252};
+//     const Double_t ksecY10[10]={+1.6298,+1.6298,+1.6298,+1.6298,+1.6298,+1.6298,+1.6298,+1.6298,+1.6298,+1.6298};
+//     const Double_t ksecX11[10]={-1.0445,-1.0445,-1.0445,-1.0445,-1.0445,-1.0445,-1.0445,-1.0445,-1.0445,-1.0445};
+//     const Double_t ksecY11[10]={+1.3162,+1.3162,+1.3162,+1.3162,+1.3162,+1.3162,+1.3162,+1.3162,+1.3162,+1.3162};
+//     const Double_t ksecX12[10]={-2.2276,-2.2276,-2.2276,-2.2276,-2.2276,-2.2276,-2.2276,-2.2276,-2.2276,-2.2276};
+//     const Double_t ksecY12[10]={+1.2948,+1.2948,+1.2948,+1.2948,+1.2948,+1.2948,+1.2948,+1.2948,+1.2948,+1.2948};
+  
+
+//    MODIFIED GEOMETRY according with partial alignment of Staves relative to Sectors
+//    last numbers: 2010/06/11 (ML)
+
+    const Double_t ksecX1[10]={-1.305917, -1.322242, -1.300649, -1.298700, -1.290830, -1.274307, -1.276433, -1.286468, -1.274381, -1.314864};
+    const Double_t ksecY1[10]={-1.997857, -2.018611, -2.005854, -2.004897, -1.995517, -2.002552, -1.995860, -2.021062, -2.012931, -2.043967};
+    const Double_t ksecX2[10]={-0.366115, -0.385562, -0.372689, -0.365682, -0.348432, -0.348442, -0.342468, -0.354071, -0.346900, -0.381275};
+    const Double_t ksecY2[10]={-1.801679, -1.808306, -1.759315, -1.778851, -1.811655, -1.747888, -1.773811, -1.792427, -1.764514, -1.820324};
+//     const Double_t ksecX1[10]={-1.305917, -1.322242, -1.300649, -1.298700, -1.290830, -1.274307, -1.276433, -1.286468, -1.274381, -1.325864};
+//     const Double_t ksecY1[10]={-1.997857, -2.018611, -2.005854, -2.004897, -1.995517, -2.002552, -1.995860, -2.021062, -2.012931, -2.032967};
+//     const Double_t ksecX2[10]={-0.366115, -0.385562, -0.372689, -0.365682, -0.348432, -0.348442, -0.342468, -0.354071, -0.346900, -0.392275};
+//     const Double_t ksecY2[10]={-1.801679, -1.808306, -1.759315, -1.778851, -1.811655, -1.747888, -1.773811, -1.792427, -1.764514, -1.809324};
+    const Double_t ksecX3[10]={-0.314030, -0.315531, -0.347521, -0.337675, -0.300420, -0.378487, -0.330729, -0.330850, -0.362360, -0.321097};
+    const Double_t ksecY3[10]={-1.452488, -1.460418, -1.447060, -1.443146, -1.472410, -1.430019, -1.469073, -1.472048, -1.462010, -1.444355};
+    const Double_t ksecX4[10]={1.124299, 1.124162, 1.089523, 1.095520, 1.136171, 1.058616, 1.105626, 1.106433, 1.077455, 1.117946};
+    const Double_t ksecY4[10]={-1.458714, -1.452649, -1.465297, -1.492717, -1.494665, -1.447732, -1.493369, -1.488126, -1.452925, -1.443447};
+    const Double_t ksecX5[10]={1.951621, 1.939284, 1.931830, 1.935235, 1.952206, 1.939082, 1.924822, 1.940114, 1.918160, 1.960017};
+    const Double_t ksecY5[10]={1.092731, 1.118870, 1.129765, 1.129422, 1.081511, 1.127387, 1.103960, 1.101784, 1.121428, 1.150110};
+    const Double_t ksecX6[10]={1.070070, 1.048297, 1.035920, 1.049049, 1.083621, 1.045882, 1.050399, 1.067823, 1.037967, 1.070850};
+    const Double_t ksecY6[10]={1.667590, 1.678571, 1.681383, 1.696892, 1.676520, 1.683470, 1.689988, 1.691111, 1.698432, 1.712770};
+    const Double_t ksecX7[10]={1.139398, 1.150471, 1.150074, 1.132807, 1.150192, 1.124064, 1.124335, 1.137723, 1.143056, 1.130568};
+    const Double_t ksecY7[10]={1.345588, 1.356062, 1.342468, 1.320467, 1.335807, 1.334477, 1.328622, 1.347184, 1.319861, 1.308420};
+    const Double_t ksecX8[10]={-0.096963, -0.098603, -0.095286, -0.099990, -0.075132, -0.121593, -0.108673, -0.104237, -0.092082, -0.104044};
+    const Double_t ksecY8[10]={1.751207, 1.731467, 1.726908, 1.734219, 1.766159, 1.718203, 1.741891, 1.739743, 1.728288, 1.718046};
+    const Double_t ksecX9[10]={0.047615, 0.087875, 0.034917, 0.071603, 0.026468, 0.091619, 0.051994, 0.059947, 0.079785, 0.043443};
+    const Double_t ksecY9[10]={1.414699, 1.403187, 1.399061, 1.403430, 1.435056, 1.384557, 1.397692, 1.420269, 1.391372, 1.398954};
+    const Double_t ksecX10[10]={-1.233255, -1.186874, -1.246702, -1.213368, -1.259425, -1.190067, -1.225655, -1.224171, -1.197833, -1.237182};
+    const Double_t ksecY10[10]={1.635767, 1.646249, 1.617336, 1.608928, 1.636944, 1.602583, 1.630504, 1.629065, 1.624295, 1.620934};
+    const Double_t ksecX11[10]={-1.018270, -1.031317, -0.960524, -1.001155, -1.045437, -0.986867, -1.002685, -1.017369, -1.005614, -0.985385};
+    const Double_t ksecY11[10]={1.318108, 1.330683, 1.301572, 1.314410, 1.326680, 1.295226, 1.306372, 1.309414, 1.306542, 1.307086};
+    const Double_t ksecX12[10]={-2.199004, -2.214964, -2.139247, -2.180547, -2.224505, -2.165324, -2.175883, -2.193485, -2.183227, -2.161570};
+    const Double_t ksecY12[10]={1.317677, 1.303982, 1.317057, 1.324766, 1.339537, 1.312715, 1.359642, 1.343638, 1.330234, 1.340836};
+
+
     const Double_t ksecR13  =  -0.8   * fgkmm; // external
     const Double_t ksecAngleSide13 = 36.0 * fgkDegree;
     //
@@ -749,24 +545,24 @@ void AliITSv11GeometrySPD::CarbonFiberSector(TGeoVolume *moth,
     // in the format of arrays (???)
     const Int_t ksecNPoints = (ksecNPointsPerRadii + 1) * ksecNRadii + 8;
     Double_t secX[ksecNRadii] = {
-        ksecX0,  ksecX1,  -1000.0,
-        ksecX2,  ksecX3,  -1000.0,
-        ksecX4,  ksecX5,  -1000.0,
-        ksecX6,  ksecX7,  -1000.0,
-        ksecX8,  ksecX9,  -1000.0,
-        ksecX10, ksecX11, -1000.0,
-        ksecX12, -1000.0
+        ksecX0,  ksecX1[sect],  -1000.0,
+        ksecX2[sect],  ksecX3[sect],  -1000.0,
+        ksecX4[sect],  ksecX5[sect],  -1000.0,
+        ksecX6[sect],  ksecX7[sect],  -1000.0,
+        ksecX8[sect],  ksecX9[sect],  -1000.0,
+        ksecX10[sect], ksecX11[sect], -1000.0,
+        ksecX12[sect], -1000.0
     };
     Double_t secY[ksecNRadii] = {
-        ksecY0,  ksecY1,  -1000.0,
-        ksecY2,  ksecY3,  -1000.0,
-        ksecY4,  ksecY5,  -1000.0,
-        ksecY6,  ksecY7,  -1000.0,
-        ksecY8,  ksecY9,  -1000.0,
-        ksecY10, ksecY11, -1000.0,
-        ksecY12, -1000.0
+        ksecY0,  ksecY1[sect],  -1000.0,
+        ksecY2[sect],  ksecY3[sect],  -1000.0,
+        ksecY4[sect],  ksecY5[sect],  -1000.0,
+        ksecY6[sect],  ksecY7[sect],  -1000.0,
+        ksecY8[sect],  ksecY9[sect],  -1000.0,
+        ksecY10[sect], ksecY11[sect], -1000.0,
+        ksecY12[sect], -1000.0
     };
-    Double_t secR[ksecNRadii] = { 
+    Double_t secR[ksecNRadii] = {
         ksecR0,  ksecR1,  -.5 * ksecDipLength - ksecDipRadii,
         ksecR2,  ksecR3,  -.5 * ksecDipLength - ksecDipRadii,
         ksecR4,  ksecR5,  -.5 * ksecDipLength - ksecDipRadii,
@@ -775,14 +571,7 @@ void AliITSv11GeometrySPD::CarbonFiberSector(TGeoVolume *moth,
         ksecR10, ksecR11, -.5 * ksecDipLength - ksecDipRadii,
         ksecR12, ksecR13
     };
-    /*
-      Double_t secDip[ksecNRadii] = {
-      0., 0., ksecDip0, 0., 0., ksecDip1,
-      0., 0., ksecDip2, 0., 0., ksecDip3,
-      0., 0., ksecDip4, 0., 0., ksecDip5,
-      0., 0.
-      };
-    */
+
     Double_t secX2[ksecNRadii];
     Double_t secY2[ksecNRadii];
     Double_t secR2[ksecNRadii] = {
@@ -794,15 +583,16 @@ void AliITSv11GeometrySPD::CarbonFiberSector(TGeoVolume *moth,
         ksecR10, ksecR11, ksecRCoolOut,
         ksecR12, ksecR13
     };
-    Double_t secDip2[ksecNCoolingTubeDips] = { 
-        ksecDl1, ksecDl2, ksecDl3, 
-        ksecDl4, ksecDl5, ksecDl6 
+    Double_t secDip2[ksecNCoolingTubeDips] = {
+        ksecDl1, ksecDl2, ksecDl3,
+        ksecDl4, ksecDl5, ksecDl6
     };
     Double_t secX3[ksecNRadii];
     Double_t secY3[ksecNRadii];
     const Int_t ksecDipIndex[ksecNCoolingTubeDips] = {2, 5, 8, 11, 14, 17};
     Double_t secAngleStart[ksecNRadii];
     Double_t secAngleEnd[ksecNRadii];
+    for(Int_t i = 0; i < ksecNRadii; i++)secAngleEnd[i] = 0.;
     Double_t secAngleStart2[ksecNRadii];
     Double_t secAngleEnd2[ksecNRadii];
     Double_t secAngleTurbo[ksecNCoolingTubeDips] = {0., 0., 0., 0., 0., 0.0};
@@ -886,7 +676,7 @@ void AliITSv11GeometrySPD::CarbonFiberSector(TGeoVolume *moth,
     secY2[ksecNRadii-1] = secY[ksecNRadii-1];
     secX3[ksecNRadii-1] = secX[ksecNRadii-1];
     secY3[ksecNRadii-1] = secY[ksecNRadii-1];
-     
+
     // find location of cooling tube centers
     for(i = 0; i < ksecNCoolingTubeDips; i++) {
         j = ksecDipIndex[i];
@@ -898,7 +688,7 @@ void AliITSv11GeometrySPD::CarbonFiberSector(TGeoVolume *moth,
         t  = secDip2[i] / t0;
         a  = x0+(x1-x0) * t;
         b  = y0+(y1-y0) * t;
-        if(i == 0) { 
+        if(i == 0) {
             // get location of tube center->Surface for locating
             // this sector around the beam pipe.
             // This needs to be double checked, but I need my notes for that.
@@ -911,7 +701,7 @@ void AliITSv11GeometrySPD::CarbonFiberSector(TGeoVolume *moth,
             secY[j]  = b - TMath::Sign(2.0*ksecDipRadii,y1-y0) * (x1-x0)/t0;
             secX2[j] = a + TMath::Abs(y1-y0) * ksecTl/t0;
             secY2[j] = b - TMath::Sign(ksecTl,y1-y0) * (x1-x0) / t0;
-            secX3[j] = a + TMath::Abs(y1-y0) * 
+            secX3[j] = a + TMath::Abs(y1-y0) *
                        (2.0*ksecDipRadii-0.5*ksecCoolTubeFlatY)/t0;
             secY3[j] = b - TMath::Sign(2.0*ksecDipRadii-0.5*ksecCoolTubeFlatY,
                                        y1-y0)*(x1-x0)/t0;
@@ -925,7 +715,7 @@ void AliITSv11GeometrySPD::CarbonFiberSector(TGeoVolume *moth,
             secY3[j] = b + TMath::Sign(2.0*ksecDipRadii-0.5*ksecCoolTubeFlatY,
                                        y1-y0)*(x1-x0)/t0;
         } // end if(a+b*(a-x0)/(b-y0)>0.0)
-          
+
           // Set up Start and End angles to correspond to start/end of dips.
         t1 = (secDip2[i]-TMath::Abs(secR[j])) / t0;
         secAngleStart[j] =TMath::RadToDeg()*TMath::ATan2(y0+(y1-y0)*t1-secY[j],
@@ -940,7 +730,7 @@ void AliITSv11GeometrySPD::CarbonFiberSector(TGeoVolume *moth,
         if (secAngleEnd[j]>secAngleStart[j]) secAngleEnd[j] -= 360.0;
         secR[j] = TMath::Sqrt(secR[j]*secR[j]+4.0*ksecDipRadii*ksecDipRadii);
     } // end for i
-     
+
     // Special cases
     secAngleStart2[8] -= 360.;
     secAngleStart2[11] -= 360.;
@@ -967,7 +757,7 @@ void AliITSv11GeometrySPD::CarbonFiberSector(TGeoVolume *moth,
             yp[j][k] = y0+(y1-y0) * t;
         } // end for k
         secAngleTurbo[i] = -TMath::RadToDeg() * TMath::ATan2(y1-y0, x1-x0);
-        if(GetDebug(3)) { 
+        if(GetDebug(3)) {
             AliInfo(
                 Form("i=%d -- angle=%f -- x0,y0=(%f, %f) -- x1,y1=(%f, %f)",
                      i, secAngleTurbo[i], x0, y0, x1, y1));
@@ -978,7 +768,7 @@ void AliITSv11GeometrySPD::CarbonFiberSector(TGeoVolume *moth,
     sA0->DefinePolygon(m, xpp, ypp);
     sA0->DefineSection(0, -ksecDz);
     sA0->DefineSection(1,  ksecDz);
-     
+
     // store the edges of each XY segment which defines
     // one of the plane zones where staves will have to be placed
     fSPDsectorX0.Set(ksecNCoolingTubeDips);
@@ -990,7 +780,7 @@ void AliITSv11GeometrySPD::CarbonFiberSector(TGeoVolume *moth,
         // Find index in xpp[] and ypp[] corresponding to where the
         // SPD ladders are to be attached. Order them according to
         // the ALICE numbering schema. Using array of indexes (+-1 for
-        // cooling tubes. For any "bend/dip/edge, there are 
+        // cooling tubes. For any "bend/dip/edge, there are
         // ksecNPointsPerRadii+1 points involved.
         if(i == 0) j = 1;
         else if (i == 1) j = 0;
@@ -1003,7 +793,7 @@ void AliITSv11GeometrySPD::CarbonFiberSector(TGeoVolume *moth,
         fSPDsectorX1[i] = sA0->GetX(ixy1);
         fSPDsectorY1[i] = sA0->GetY(ixy1);
     } // end for i
-     
+
     //printf("SectorA#%d ",0);
     InsidePoint(xpp[m-1],ypp[m-1],xpp[0],ypp[0],xpp[1],ypp[1],ksecCthick,
                 xpp2[0],ypp2[0]);
@@ -1042,7 +832,7 @@ void AliITSv11GeometrySPD::CarbonFiberSector(TGeoVolume *moth,
     // Error in TGeoEltu. Semi-axis X must be < Semi-axis Y (?).
     sTA0 = new TGeoEltu("ITS SPD Cooling Tube TA0", 0.5 * ksecCoolTubeFlatY,
                         0.5 * ksecCoolTubeFlatX, ksecDz);
-    sTA1 = new TGeoEltu("ITS SPD Cooling Tube coolant TA1", 
+    sTA1 = new TGeoEltu("ITS SPD Cooling Tube coolant TA1",
                         sTA0->GetA() - ksecCoolTubeThick,
                         sTA0->GetB()-ksecCoolTubeThick,ksecDz);
     SPDsectorShape(ksecNRadii,secX2,secY2,secR2,secAngleStart2,secAngleEnd2,
@@ -1054,8 +844,12 @@ void AliITSv11GeometrySPD::CarbonFiberSector(TGeoVolume *moth,
     sB0->DefineSection(1, ksecDz + ksecZEndLen);
 
     //printf("SectorB#%d ",0);
+  // Points around the most sharpened tips have to be avoided - M.S. 24 feb 09
+    const Int_t nSpecialPoints = 5;
+    const Int_t kSpecialPoints[nSpecialPoints] = {7, 17, 47, 62, 77};
+    Int_t i2 = 0;
     InsidePoint(xpp[m-1],ypp[m-1],xpp[0],ypp[0],xpp[1],ypp[1],
-                ksecCthick2,xpp2[0],ypp2[0]);
+                ksecCthick2,xpp2[i2],ypp2[i2]);
     for(i = 1; i < m - 1; i++) {
         t = ksecCthick2;
         for(k = 0; k < ksecNCoolingTubeDips; k++)
@@ -1065,15 +859,23 @@ void AliITSv11GeometrySPD::CarbonFiberSector(TGeoVolume *moth,
                      ksecNPointsPerRadii == i))
                     t = ksecRCoolOut-ksecRCoolIn;
         //printf("SectorB#%d ",i);
-        InsidePoint(xpp[i-1],ypp[i-1],xpp[i],ypp[i],xpp[i+1],ypp[i+1],t,
-                    xpp2[i],ypp2[i]);
+       Bool_t useThisPoint = kTRUE;
+       for(Int_t ii = 0; ii < nSpecialPoints; ii++)
+         if ( (i == kSpecialPoints[ii] - 1) ||
+              (i == kSpecialPoints[ii] + 1)   ) useThisPoint = kFALSE;
+       if (useThisPoint) {
+         i2++;
+         InsidePoint(xpp[i-1],ypp[i-1],xpp[i],ypp[i],xpp[i+1],ypp[i+1],t,
+                     xpp2[i2],ypp2[i2]);
+       }
     }// end for i
     //printf("SectorB#%d ",m);
+    i2++;
     InsidePoint(xpp[m-2],ypp[m-2],xpp[m-1],ypp[m-1],xpp[0],ypp[0],
-                ksecCthick2,xpp2[m-1],ypp2[m-1]);
+                ksecCthick2,xpp2[i2],ypp2[i2]);
     sB1 = new TGeoXtru(2);
     sB1->SetName("ITS SPD Carbon fiber support Sector Air End B1");
-    sB1->DefinePolygon(m, xpp2, ypp2);
+    sB1->DefinePolygon(i2+1, xpp2, ypp2);
     sB1->DefineSection(0,sB0->GetZ(0));
     sB1->DefineSection(1,sB0->GetZ(1)-ksecCthick2);
     const Double_t kspdEndHoleRadius1=5.698*fgkmm;
@@ -1123,7 +925,7 @@ void AliITSv11GeometrySPD::CarbonFiberSector(TGeoVolume *moth,
         sB1->InspectShape();
         sB2->InspectShape();
     } // end if(GetDebug(3))
-     
+
     // create the assembly of the support and place staves on it
     TGeoVolumeAssembly *vM0 = new TGeoVolumeAssembly(
                                          "ITSSPDSensitiveVirtualvolumeM0");
@@ -1197,7 +999,7 @@ void AliITSv11GeometrySPD::CarbonFiberSector(TGeoVolume *moth,
     vTB1->SetLineWidth(1);
     vTB1->SetFillColor(vTB1->GetLineColor());
     vTB1->SetFillStyle(4050); // 0% transparent
-     
+
     // add volumes to mother container passed as argument of this method
     moth->AddNode(vM0,1,0); // Add virtual volume to mother
     vA0->AddNode(vA1,1,0); // Put air inside carbon fiber.
@@ -1240,26 +1042,9 @@ void AliITSv11GeometrySPD::CarbonFiberSector(TGeoVolume *moth,
     vM0->AddNode(vB3,1,rotrans); // Put Mounting bracket on sector
     rotrans = new TGeoCombiTrans("",x0,y0,-z0,rot);
     vM0->AddNode(vB3,2,rotrans); // Put Mounting bracket on sector
-    /*
-    j = 0; // right side, find point with largest x value
-    x1 = sB0->GetX(0);
-    for(i=1;i<sB0->GetNvert();i++)if(sB0->GetX(i)>x1) {j=i;x1=sB0->GetX(i);}
-    j--; // Too big by 1
-    //t = -TMath::RadToDeg()*TMath::ATan2(
-    //                               sB0->GetX(j)-sB0->GetX(j-1),
-    //                               sB0->GetY(j)-sB0->GetY(j-1));
-    */
     t *= -1.0;
     rot = new TGeoRotation("",t,0.0,0.0); // z axis rotation
-    /*  // this way gets correct orientation but wrong "height"
-    x0 = 0.5*(sB0->GetX(j)+sB0->GetX(j-1))+
-        sB3->GetDX()*TMath::Cos(t*TMath::DegToRad());
-    y0 = 0.5*(sB0->GetY(j)+sB0->GetY(j-1))+
-        sB3->GetDX()*TMath::Sin(t*TMath::DegToRad());
-    z0 = sB0->GetZ(0)+sB3->GetDZ();
-    */ // I don't understand the need for this factor 3.5.
-    // posibly the SPD sector as coded isn't symetric which the
-    // plans would suggest.
+  
     x0 = -0.5*(sB0->GetX(0)+sB0->GetX(sB0->GetNvert()-1))-3.5*
         sB3->GetDX()*TMath::Cos(t*TMath::DegToRad());
     y0 = 0.5*(sB0->GetY(0)+sB0->GetY(sB0->GetNvert()-1))-3.5*
@@ -1310,7 +1095,7 @@ Bool_t AliITSv11GeometrySPD::CFHolePoints(Double_t s,Double_t r1,
     bc = 2.*l*r2;
     if(bc==0.0) {printf("bc=0 l=%e r2=%e\n",l,r2);return kFALSE;}
     betac = TMath::ACos(ac/bc);
-    alphac = TMath::Sqrt(bc*bc-ac*ac)/(2.*l*r1);
+    alphac = TMath::Sqrt((bc-ac)*(bc+ac))/(2.*l*r1);
     scb = r2*betac;
     sca = r1*alphac;
     t = r1*0.5*TMath::Pi() - sca + scb;
@@ -1359,13 +1144,13 @@ Bool_t AliITSv11GeometrySPD::GetSectorMountingPoints(Int_t index,Double_t &x0,
     // ---
     // Returns kTRUE if no problems encountered.
     // Returns kFALSE if a problem was encountered (e.g.: shape not found).
-    // 
+    //
     Int_t isize = fSPDsectorX0.GetSize();
 
     x0 = x1 = y0 = y1 = 0.0;
     if(index < 0 || index > isize) {
-        AliError(Form("index = %d: allowed 0 --> %", index, isize));
-        return kFALSE;
+      AliError(Form("index = %d: allowed 0 --> %d", index, isize));
+      return kFALSE;
     } // end if(index<0||index>isize)
     x0 = fSPDsectorX0[index];
     x1 = fSPDsectorX1[index];
@@ -1374,9 +1159,9 @@ Bool_t AliITSv11GeometrySPD::GetSectorMountingPoints(Int_t index,Double_t &x0,
     return kTRUE;
 }
 //______________________________________________________________________
-void AliITSv11GeometrySPD::SPDsectorShape(Int_t n,const Double_t *xc, 
+void AliITSv11GeometrySPD::SPDsectorShape(Int_t n,const Double_t *xc,
                               const Double_t *yc,  const Double_t *r,
-                              const Double_t *ths, const Double_t *the, 
+                              const Double_t *ths, const Double_t *the,
                       Int_t npr, Int_t &m, Double_t **xp, Double_t **yp) const
 {
     //
@@ -1391,7 +1176,7 @@ void AliITSv11GeometrySPD::SPDsectorShape(Int_t n,const Double_t *xc,
     //   Double_t *the  array of ending angles [degrees].
     //   Int_t     npr  the number of lines segments to aproximate the arc.
     // Outputs (arguments passed by reference):
-    //   Int_t       m    the number of enetries in the arrays *xp[npr+1] 
+    //   Int_t       m    the number of enetries in the arrays *xp[npr+1]
     //                    and *yp[npr+1].
     //   Double_t **xp    array of x coordinate values of the line segments
     //                    which make up the SPD support sector shape.
@@ -1457,11 +1242,11 @@ TGeoVolume* AliITSv11GeometrySPD::CreateLadder(Int_t layer,TArrayD &sizes,
     //      dimensions related to this object:
     //      size[0] = 'thickness' (the smallest dimension)
     //      size[1] = 'length' (the direction along the ALICE Z axis)
-    //      size[2] = 'width' (extension in the direction perp. to the 
+    //      size[2] = 'width' (extension in the direction perp. to the
     //                         above ones)
     //  3 - the used TGeoManager
 
-    // ** CRITICAL CHECK **    
+    // ** CRITICAL CHECK **
     // layer number can be ONLY 1 or 2
     if (layer != 1 && layer != 2) AliFatal("Layer number MUST be 1 or 2");
 
@@ -1470,8 +1255,8 @@ TGeoVolume* AliITSv11GeometrySPD::CreateLadder(Int_t layer,TArrayD &sizes,
     TGeoMedium *medSPDSiChip = GetMedium("SPD SI CHIP$",mgr); // SPD SI CHIP
     TGeoMedium *medSi        = GetMedium("SI$",mgr);
     TGeoMedium *medBumpBond  = GetMedium("COPPER$",mgr);  // ??? BumpBond
-    
-    // ** SIZES **    
+
+    // ** SIZES **
     Double_t chipThickness  = fgkmm *  0.150;
     Double_t chipWidth      = fgkmm * 15.950;
     Double_t chipLength     = fgkmm * 13.600;
@@ -1479,7 +1264,7 @@ TGeoVolume* AliITSv11GeometrySPD::CreateLadder(Int_t layer,TArrayD &sizes,
     Double_t sensThickness  = fgkmm *  0.200;
     Double_t sensLength     = fgkmm * 69.600;
     Double_t sensWidth      = fgkmm * 12.800;
-    Double_t guardRingWidth = fgkmm *  0.560; // a border of this thickness 
+    Double_t guardRingWidth = fgkmm *  0.560; // a border of this thickness
                                               // all around the sensor
     Double_t bbLength       = fgkmm * 0.042;
     Double_t bbWidth        = sensWidth;
@@ -1501,9 +1286,9 @@ TGeoVolume* AliITSv11GeometrySPD::CreateLadder(Int_t layer,TArrayD &sizes,
     // While creating this volume, since it is a sensitive volume,
     // we must respect some standard criteria for its local reference frame.
     // Local X must correspond to x coordinate of the sensitive volume:
-    // this means that we are going to create the container with a local 
+    // this means that we are going to create the container with a local
     // reference system that is **not** in the middle of the box.
-    // This is accomplished by calling the shape constructor with an 
+    // This is accomplished by calling the shape constructor with an
     // additional option ('originShift'):
     Double_t xSens = 0.5 * (width - sensWidth - 2.0*guardRingWidth);
     Double_t originShift[3] = {-xSens, 0., 0.};
@@ -1518,7 +1303,7 @@ TGeoVolume* AliITSv11GeometrySPD::CreateLadder(Int_t layer,TArrayD &sizes,
     // the sensor as well
     TGeoVolume *volSens = mgr->MakeBox(GetSenstiveVolumeName(layer),medSi,
                              0.5*sensWidth,0.5*sensThickness,0.5*sensLength);
-    // the guard ring shape is the subtraction of two boxes with the 
+    // the guard ring shape is the subtraction of two boxes with the
     // same center.
     TGeoBBox  *shIn = new TGeoBBox(0.5*sensWidth,sensThickness,0.5*sensLength);
     TGeoBBox  *shOut = new TGeoBBox(0.5*sensWidth+guardRingWidth,
@@ -1541,7 +1326,7 @@ TGeoVolume* AliITSv11GeometrySPD::CreateLadder(Int_t layer,TArrayD &sizes,
     // sensor is translated along thickness (X) and width (Y)
     Double_t ySens = 0.5 * (thickness - sensThickness);
     Double_t zSens = 0.0;
-    // we want that the x of the ladder is the same as the one of 
+    // we want that the x of the ladder is the same as the one of
     // its sensitive volume
     TGeoTranslation *trSens = new TGeoTranslation(0.0, ySens, zSens);
     // bump bonds are translated along all axes:
@@ -1567,11 +1352,11 @@ TGeoVolume* AliITSv11GeometrySPD::CreateLadder(Int_t layer,TArrayD &sizes,
     y = 0.5 * (chipThickness - thickness);
     z = 0.0;
     for (i = 0; i < 5; i++) {
-        z = -0.5*length + guardRingWidth 
+        z = -0.5*length + guardRingWidth
             + (Double_t)i*chipSpacing + ((Double_t)(i) + 0.5)*chipLength;
         trChip[i] = new TGeoTranslation(x, y, z);
     } // end ofr i
-    
+
     // add nodes to container
     container->AddNode(volSens, 1, trSens);
     container->AddNode(volBorder, 1, trSens);
@@ -1581,158 +1366,6 @@ TGeoVolume* AliITSv11GeometrySPD::CreateLadder(Int_t layer,TArrayD &sizes,
     return container;
 }
 
-/*
-//______________________________________________________________________
-TGeoVolume* AliITSv11GeometrySPD::CreateLadder
-        (Int_t layer, TArrayD &sizes, TGeoManager *mgr) const
-{
-    //
-    // Creates the "ladder" = silicon sensor + 5 chips.
-    // Returns a TGeoVolume containing the following components:
-    //  - the sensor (TGeoBBox), whose name depends on the layer
-    //  - 5 identical chips (TGeoBBox)
-    //  - a guard ring around the sensor (subtraction of TGeoBBoxes),
-    //    which is separated from the rest of sensor because it is not
-    //    a sensitive part
-    //  - bump bondings (TGeoBBox stripes for the whole width of the
-    //    sensor, one per column).
-    // ---
-    // Arguments:
-    //  1 - the owner layer (MUST be 1 or 2 or a fatal error is raised)
-    //  2 - a TArrayD passed by reference, which will contain relevant
-    //      dimensions related to this object:
-    //      size[0] = 'thickness' (the smallest dimension)
-    //      size[1] = 'length' (the direction along the ALICE Z axis)
-    //      size[2] = 'width' (extension in the direction perp. to the 
-    //                         above ones)
-    //  3 - the used TGeoManager
-
-    // ** CRITICAL CHECK ******************************************************
-    // layer number can be ONLY 1 or 2
-    if (layer != 1 && layer != 2) AliFatal("Layer number MUST be 1 or 2");
-
-    // ** MEDIA ***************************************************************
-    
-    TGeoMedium *medAir       = GetMedium("AIR$",mgr);
-    TGeoMedium *medSPDSiChip = GetMedium("SPD SI CHIP$",mgr); // SPD SI CHIP
-    TGeoMedium *medSi        = GetMedium("SI$",mgr);
-    TGeoMedium *medBumpBond  = GetMedium("COPPER$",mgr);  // ??? BumpBond
-
-    // ** SIZES ***************************************************************
-        
-    Double_t chipThickness  = fgkmm *  0.150;
-    Double_t chipWidth      = fgkmm * 15.950;
-    Double_t chipLength     = fgkmm * 13.600;
-    Double_t chipSpacing    = fgkmm *  0.400; // separation of chips along Z
-    Double_t sensThickness  = fgkmm *  0.200;
-    Double_t sensLength     = fgkmm * 69.600;
-    Double_t sensWidth      = fgkmm * 12.800;
-    Double_t guardRingWidth = fgkmm *  0.560; // guard ring around sensor
-    Double_t bbLength       = fgkmm * 0.042;
-    Double_t bbWidth        = sensWidth;
-    Double_t bbThickness    = fgkmm * 0.012;
-    Double_t bbPos          = 0.080;          // Z position w.r. to left pixel edge
-    
-    // the three dimensions of the box which contains the ladder
-    // are returned in the 'sizes' argument, and are used for volumes positionement
-    // for readability purpose, they are linked by reference to a more meaningful name
-    sizes.Set(3);
-    Double_t &thickness = sizes[0];
-    Double_t &length = sizes[1];
-    Double_t &width = sizes[2];
-    // the container is a box which exactly enclose all the stuff
-    width = chipWidth;
-    length = sensLength + 2.0*guardRingWidth;
-    thickness = sensThickness + chipThickness + bbThickness;
-
-    // ** VOLUMES *************************************************************
-    
-    // This is a sensitive volume.
-    // Local X must correspond to x coordinate of the sensitive volume:
-    // to respect this, the origin of the local reference system 
-    // must be shifted from the middle of the box, using 
-    // an additional option ('originShift') when creating the container shape:
-    Double_t xSens = 0.5 * (width - sensWidth - 2.0*guardRingWidth);
-    Double_t originShift[3] = {-xSens, 0., 0.};
-    
-    // now the container is a TGeoBBox with this shift,
-    // and the volume is made of air (it does not exist in reality)
-    TGeoBBox *shLadder = new TGeoBBox(0.5*width, 0.5*thickness, 0.5*length, originShift);
-    TGeoVolume *vLadder = new TGeoVolume(Form("ITSSPDlay%d-Ladder", layer), shLadder, medAir);
-    
-    // the chip is a common box
-    TGeoVolume *vChip = mgr->MakeBox("ITSSPDchip", medSPDSiChip, 
-                                     0.5*chipWidth, 0.5*chipThickness, 0.5*chipLength);
-    
-    // to build the sensor with its guard ring, we create a TGeoBBox with the size
-    // of the sensor + guard ring, and we insert the true sensor into it as an 
-    // internal node: this simplifies the implementation with the same result
-    TGeoVolume *vSensGuard = mgr->MakeBox(Form("%s-guardRing", GetSenstiveVolumeName(layer)),
-                                          medSi, 
-                                          0.5*sensWidth + guardRingWidth,
-                                          0.5*sensThickness,
-                                          0.5*sensLength + guardRingWidth);
-    TGeoVolume *vSens = mgr->MakeBox(GetSenstiveVolumeName(layer), medSi,
-                                     0.5*sensWidth,0.5*sensThickness,0.5*sensLength);
-    vSensGuard->AddNode(vSens, 0);
-    vSensGuard->SetTransparency(50);
-    
-    // bump bond is a common box for one whole column
-    TGeoVolume *vBB = mgr->MakeBox("ITSSPDbb", medBumpBond,
-                                   0.5*bbWidth, 0.5*bbThickness, 0.5*bbLength);
-    
-    // set colors of all objects for visualization
-    vLadder->SetLineColor(kRed);
-    vSens->SetLineColor(kYellow + 1);
-    vChip->SetLineColor(kGreen);
-    vSensGuard->SetLineColor(kYellow + 3);
-    vBB->SetLineColor(kGray);
-
-    // ** MOVEMENTS **
-    // sensor is translated along thickness (Y) and width (X)
-    Double_t ySens = 0.5 * (thickness - sensThickness);
-    Double_t zSens = 0.0;
-    // we want that the x of the ladder is the same as the one of 
-    // its sensitive volume
-    TGeoTranslation *trSens = new TGeoTranslation(0.0, ySens, zSens);
-    // bump bonds are translated along all axes:
-    // keep same Y used for sensors, but change the Z
-    TGeoTranslation *trBB[160];
-    Double_t x =  0.0;
-    Double_t y =  0.5 * (thickness - bbThickness) - sensThickness;
-    Double_t z = -0.5 * sensLength + guardRingWidth + fgkmm*0.425 - bbPos;
-    Int_t i;
-    for (i = 0; i < 160; i++) {
-        trBB[i] = new TGeoTranslation(x, y, z);
-        switch(i) {
-            case  31:case  63:case  95:case 127:
-                z += fgkmm * 0.625 + fgkmm * 0.2;
-                break;
-            default:
-                z += fgkmm * 0.425;
-        } // end switch
-    } // end for i
-    // the chips are translated along the length (Z) and thickness (X)
-    TGeoTranslation *trChip[5] = {0, 0, 0, 0, 0};
-    x = -xSens;
-    y = 0.5 * (chipThickness - thickness);
-    z = 0.0;
-    for (i = 0; i < 5; i++) {
-        z = -0.5*length + guardRingWidth 
-                + (Double_t)i*chipSpacing + ((Double_t)(i) + 0.5)*chipLength;
-        trChip[i] = new TGeoTranslation(x, y, z);
-    } // end ofr i
-    
-    // add nodes to container
-    vLadder->AddNode(vSensGuard, 1, trSens);
-    //vLadderAddNode(volBorder, 1, trSens);
-    for (i = 0; i < 160; i++) vLadder->AddNode(vBB,i+1,trBB[i]);
-    for (i = 0; i < 5; i++) vLadder->AddNode(vChip,i+3,trChip[i]);
-    // return the container
-    return vLadder;
-}
-*/
-
 //______________________________________________________________________
 TGeoVolume* AliITSv11GeometrySPD::CreateClip(TArrayD &sizes,Bool_t isDummy,
                                              TGeoManager *mgr) const
@@ -1742,16 +1375,20 @@ TGeoVolume* AliITSv11GeometrySPD::CreateClip(TArrayD &sizes,Bool_t isDummy,
     // They have a complicated shape which is approximated by a TGeoXtru
     // Implementation of a single clip over an half-stave.
     // It has a complicated shape which is approximated to a section like this:
-    //   
+    //
     //     6
     //     /\   .
     //  7 //\\  5
     //    / 1\\___________________4
     //   0    \___________________
     //        2                   3
-    // with a finite thickness for all the shape 
+    // with a finite thickness for all the shape
     // Its local reference frame is such that point A corresponds to origin.
-    // 
+    //
+
+  // MODIFIED geometry
+    Double_t sposty = fgkmm * -0.5; // lower internal side to avoid overlaps with modified geometry
+
     Double_t fullLength      = fgkmm * 12.6;    // = x4 - x0
     Double_t flatLength      = fgkmm *  5.4;    // = x4 - x3
     Double_t inclLongLength  = fgkmm *  5.0;    // = 5-6
@@ -1759,22 +1396,22 @@ TGeoVolume* AliITSv11GeometrySPD::CreateClip(TArrayD &sizes,Bool_t isDummy,
     Double_t fullHeight      = fgkmm *  2.8;    // = y6 - y3
     Double_t thickness       = fgkmm *  0.18;    // thickness
     Double_t totalLength     = fgkmm * 52.0;    // total length in Z
-    Double_t holeSize        = fgkmm *  5.0;    // dimension of cubic 
+    Double_t holeSize        = fgkmm *  5.0;    // dimension of cubic
                                                 // hole inserted for pt1000
     Double_t angle1          = 27.0;            // supplementary of angle DCB
     Double_t angle2;                            // angle DCB
     Double_t angle3;                            // angle of GH with vertical
+
     angle2 = 0.5 * (180.0 - angle1);
-    angle3 = 90.0 - TMath::ACos(fullLength - flatLength - 
-                                inclLongLength*TMath::Cos(angle1)) * 
+    angle3 = 90.0 - TMath::ACos(fullLength - flatLength -
+                                inclLongLength*TMath::Cos(angle1)) *
                                 TMath::RadToDeg();
     angle1 *= TMath::DegToRad();
     angle2 *= TMath::DegToRad();
     angle3 *= TMath::DegToRad();
 
     Double_t x[8], y[8];
+
     x[0] =  0.0;
     x[1] = x[0] + fullLength - flatLength - inclLongLength*TMath::Cos(angle1);
     x[2] = x[0] + fullLength - flatLength;
@@ -1783,7 +1420,7 @@ TGeoVolume* AliITSv11GeometrySPD::CreateClip(TArrayD &sizes,Bool_t isDummy,
     x[5] = x[4] - flatLength + thickness * TMath::Cos(angle2);
     x[6] = x[1];
     x[7] = x[0];
+
     y[0] = 0.0;
     y[1] = y[0] + inclShortLength * TMath::Cos(angle3);
     y[2] = y[1] - inclLongLength * TMath::Sin(angle1);
@@ -1792,7 +1429,10 @@ TGeoVolume* AliITSv11GeometrySPD::CreateClip(TArrayD &sizes,Bool_t isDummy,
     y[5] = y[4];
     y[6] = y[1] + thickness;
     y[7] = y[0] + thickness;
+
+    y[0] += sposty;
+    y[7] += sposty;
+
     sizes.Set(7);
     sizes[0] = totalLength;
     sizes[1] = fullHeight;
@@ -1812,7 +1452,7 @@ TGeoVolume* AliITSv11GeometrySPD::CreateClip(TArrayD &sizes,Bool_t isDummy,
     shClip->DefinePolygon(8, x, y);
     shClip->DefineSection(0, -0.5*totalLength, 0., 0., 1.0);
     shClip->DefineSection(1,  0.5*totalLength, 0., 0., 1.0);
+
     TGeoBBox *shHole = new TGeoBBox("ITSSPDSHClipHole",0.5*holeSize,
                                     0.5*holeSize,0.5*holeSize);
     TGeoTranslation *tr1 = new TGeoTranslation("ITSSPDTRClipHole1",x[2],0.0,
@@ -1838,14 +1478,83 @@ TGeoVolume* AliITSv11GeometrySPD::CreateClip(TArrayD &sizes,Bool_t isDummy,
     TGeoVolume *vClip = new TGeoVolume("ITSSPDclip", shClipHole, mat);
     vClip->SetLineColor(kGray + 2);
     return vClip;
-}//______________________________________________________________________
+}
+
+//______________________________________________________________________
+TGeoVolume* AliITSv11GeometrySPD::CreatePatchPanel(TArrayD &sizes,
+                                                  TGeoManager *mgr) const
+{
+    //
+    // Creates the patch panel approximated with a "L"-shaped TGeoXtru
+    // with a finite thickness for all the shape
+    // Its local reference frame is such that point A corresponds to origin.
+    //
+    Double_t hLength         = fgkmm *  50.0;    // horizontal length
+    Double_t vLength         = fgkmm *  50.0;    // vertical length
+    Double_t angle           = 88.3;             // angle between hor and vert
+    Double_t thickness       = fgkmm *   4.0;    // thickness
+    Double_t width           = fgkmm * 100.0;    // width looking from cone
+
+    Double_t x[7], y[7];
+
+    y[0] =  0.0;
+    y[1] = y[0] + hLength;
+    y[2] = y[1];
+    y[3] = y[0] + thickness;
+    y[4] = y[3] + vLength * TMath::Cos(angle*TMath::DegToRad());
+    y[5] = y[4] - thickness / TMath::Sin(angle*TMath::DegToRad());
+    y[6] = y[0];
+
+    x[0] = 0.0;
+    x[1] = x[0];
+    x[2] = x[1] + thickness;
+    x[3] = x[2];
+    x[4] = x[3] + vLength * TMath::Sin(angle*TMath::DegToRad());
+    x[5] = x[4];
+    x[6] = x[0] + thickness;
+
+    sizes.Set(3);
+    sizes[0] = hLength;
+    sizes[1] = vLength;
+    sizes[2] = thickness;
+
+    TGeoXtru *shPatch = new TGeoXtru(2);
+    shPatch->SetName("ITSSPDpatchShape1");
+    shPatch->DefinePolygon(7, x, y);
+    shPatch->DefineSection(0, -0.5*width, 0., 0., 1.0);
+    shPatch->DefineSection(1,  0.5*width, 0., 0., 1.0);
+    
+    /*
+    Double_t subThickness = 10.0 * fgkmm;
+    Double_t subWidth     = 55.0 * fgkmm;
+    new TGeoBBox("ITSSPDpatchShape2", 0.5*subThickness, 60.0 * fgkmm, 0.5*subWidth);
+    TGeoRotation *rotSub = new TGeoRotation(*gGeoIdentity);
+    rotSub->SetName("shPatchSubRot");
+    rotSub->RotateZ(50.0);
+    rotSub->RegisterYourself();
+    TGeoCombiTrans *trSub = new TGeoCombiTrans(0.26*hLength, 0.26*vLength, 0.0, rotSub);
+    trSub->SetName("shPatchSubTr");
+    trSub->RegisterYourself();
+    
+    TGeoCompositeShape *shPatchFinal = new TGeoCompositeShape("ITSSPDpatchShape1-(ITSSPDpatchShape2:shPatchSubTr)");
+    */
+
+    TGeoMedium *mat = GetMedium("AL$", mgr);
+    //TGeoVolume *vPatch = new TGeoVolume("ITSSPDpatchPanel", shPatchFinal, mat);
+    TGeoVolume *vPatch = new TGeoVolume("ITSSPDpatchPanel", shPatch, mat);
+    vPatch->SetLineColor(kAzure);
+    
+    return vPatch;
+}
+
+//___________________________________________________________________
 TGeoCompositeShape* AliITSv11GeometrySPD::CreateGroundingFoilShape
                        (Int_t itype,Double_t &length,Double_t &width,
                         Double_t thickness,TArrayD &sizes)
 {
     //
-    // Creates the typical composite shape of the grounding foil: 
-    // 
+    // Creates the typical composite shape of the grounding foil:
+    //
     //  +---------------------------------------------------------+
     //  |                         5           6      9            |
     //  |                         +-----------+      +------------+ 10
@@ -1857,38 +1566,38 @@ TGeoCompositeShape* AliITSv11GeometrySPD::CreateGroundingFoilShape
     //       0
     //       Z                                                    + 11
     //
-    // This shape is used 4 times: two layers of glue, one in kapton 
-    // and one in aluminum, taking into account that the aliminum 
+    // This shape is used 4 times: two layers of glue, one in kapton
+    // and one in aluminum, taking into account that the aliminum
     // layer has small differences in the size of some parts.
     // ---
-    // In order to overcome problems apparently due to a large number 
-    // of points, the shape creation is done according the following 
+    // In order to overcome problems apparently due to a large number
+    // of points, the shape creation is done according the following
     // steps:
-    //    1) a TGeoBBox is created with a size right enough to contain 
+    //    1) a TGeoBBox is created with a size right enough to contain
     //       the whole shape (0-1-X-13)
-    //    2) holes are defined as other TGeoBBox which are subtracted 
+    //    2) holes are defined as other TGeoBBox which are subtracted
     //       from the main shape
-    //    3) a TGeoXtru is defined connecting the points (0-->11-->0) 
+    //    3) a TGeoXtru is defined connecting the points (0-->11-->0)
     //       and is also subtracted from the main shape
     // ---
-    // The argument ("type") is used to choose between all these 
+    // The argument ("type") is used to choose between all these
     // possibilities:
     //   - type = 0 --> kapton layer
     //   - type = 1 --> aluminum layer
     //   - type = 2 --> glue layer between support and GF
     //   - type = 3 --> glue layer between GF and ladders
-    // Returns: a TGeoCompositeShape which will then be used to shape 
-    // several volumes. Since TGeoXtru is used, the local reference 
+    // Returns: a TGeoCompositeShape which will then be used to shape
+    // several volumes. Since TGeoXtru is used, the local reference
     // frame of this object has X horizontal and Y vertical w.r to
     // the shape drawn above, and Z axis going perpendicularly to the screen.
-    // This is not the correct reference for the half stave, for which 
-    // the "long" dimension is Z and the "short" is X, while Y goes in 
-    // the direction of thickness. This will imply some rotations when 
+    // This is not the correct reference for the half stave, for which
+    // the "long" dimension is Z and the "short" is X, while Y goes in
+    // the direction of thickness. This will imply some rotations when
     // using the volumes created with this shape.
-    
+
     // suffix to differentiate names
     Char_t type[10];
-    
+
     // size of the virtual box containing exactly this volume
     length = fgkmm * 243.18;
     width  = fgkmm *  15.95;
@@ -1898,20 +1607,20 @@ TGeoCompositeShape* AliITSv11GeometrySPD::CreateGroundingFoilShape
     } // end if itype==1
     switch (itype) {
     case 0:
-        sprintf(type,"Kap");
+        snprintf(type,10,"Kap");
         break;
     case 1:
-        sprintf(type,"Alu");
+        snprintf(type,10, "Alu");
         break;
     case 2:
-        sprintf(type,"Glue1");
+        snprintf(type,10,"Glue1");
         break;
     case 3:
-        sprintf(type,"Glue2");
+        snprintf(type,10,"Glue2");
         break;
     }
-    // we divide the shape in several slices along the horizontal 
-    // direction (local X) here we define define the length of all 
+    // we divide the shape in several slices along the horizontal
+    // direction (local X) here we define define the length of all
     // sectors (from leftmost to rightmost)
     Int_t i;
     Double_t sliceLength[] = { 140.71,  2.48,  26.78,   4.00,
@@ -1923,8 +1632,8 @@ TGeoCompositeShape* AliITSv11GeometrySPD::CreateGroundingFoilShape
         sliceLength[5] += fgkmm * 0.4;
         sliceLength[6] -= fgkmm * 0.4;
     } // end if itype ==1
-    
-    // as shown in the drawing, we have four different widths 
+
+    // as shown in the drawing, we have four different widths
     // (along local Y) in this shape:
     Double_t widthMax  = fgkmm * 15.95;
     Double_t widthMed1 = fgkmm * 15.00;
@@ -1936,20 +1645,22 @@ TGeoCompositeShape* AliITSv11GeometrySPD::CreateGroundingFoilShape
         widthMed2 -= fgkmm * 0.4;
         widthMin  -= fgkmm * 0.4;
     } // end if itype==1
-    
+
     // create the main shape
     TGeoBBox *shGroundFull = 0;
     shGroundFull = new TGeoBBox(Form("ITSSPDSHgFoil%sFull", type),
                                 0.5*length,0.5*width, 0.5*thickness);
-    
-    // create the polygonal shape to be subtracted to give the correct 
-    // shape to the borders its vertices are defined in sugh a way that 
-    // this polygonal will be placed in the correct place considered 
-    // that the origin of the local reference frame is in the center 
-    // of the main box: we fix the starting point at the lower-left 
-    // edge of the shape (point 12), and add all points in order, 
+
+    if(GetDebug(5)) shGroundFull->Print(); // Avoid Coverity warning
+
+    // create the polygonal shape to be subtracted to give the correct
+    // shape to the borders its vertices are defined in sugh a way that
+    // this polygonal will be placed in the correct place considered
+    // that the origin of the local reference frame is in the center
+    // of the main box: we fix the starting point at the lower-left
+    // edge of the shape (point 12), and add all points in order,
     // following a clockwise rotation
-    
+
     Double_t x[13], y[13];
     x[ 0] = -0.5 * length + sliceLength[0];
     y[ 0] = -0.5 * widthMax;
@@ -1996,22 +1707,22 @@ TGeoCompositeShape* AliITSv11GeometrySPD::CreateGroundingFoilShape
     shGroundXtru->DefinePolygon(13, x, y);
     shGroundXtru->DefineSection(0, -thickness, 0., 0., 1.0);
     shGroundXtru->DefineSection(1,  thickness, 0., 0., 1.0);
-    
+
     // define a string which will express the algebric operations among volumes
     // and add the subtraction of this shape from the main one
     TString strComposite(Form("ITSSPDSHgFoil%sFull-(%s+", type,
                               shGroundXtru->GetName()));
-    
+
     // define the holes according to size information coming from drawings:
     Double_t holeLength = fgkmm * 10.00;
     Double_t holeWidth  = fgkmm *  7.50;
-    Double_t holeSepX0  = fgkmm *  7.05;  // separation between center 
+    Double_t holeSepX0  = fgkmm *  7.05;  // separation between center
                                           // of first hole and left border
-    Double_t holeSepXC  = fgkmm * 14.00;  // separation between the centers 
+    Double_t holeSepXC  = fgkmm * 14.00;  // separation between the centers
                                           // of two consecutive holes
-    Double_t holeSepX1  = fgkmm * 15.42;  // separation between centers of 
+    Double_t holeSepX1  = fgkmm * 15.42;  // separation between centers of
                                           // 5th and 6th hole
-    Double_t holeSepX2  = fgkmm * 22.00;  // separation between centers of 
+    Double_t holeSepX2  = fgkmm * 22.00;  // separation between centers of
                                           // 10th and 11th hole
     if (itype == 1) {
         holeSepX0  -= fgkmm * 0.2;
@@ -2026,19 +1737,18 @@ TGeoCompositeShape* AliITSv11GeometrySPD::CreateGroundingFoilShape
     sizes[4] = holeSepX1;
     sizes[5] = holeSepX2;
     sizes[6] = fgkmm * 4.40;
-    
+
     // X position of hole center (will change for each hole)
     Double_t holeX = -0.5*length;
     // Y position of center of all holes (= 4.4 mm from upper border)
     Double_t holeY = 0.5*(width - holeWidth) - widthMin;
-    
+
     // create a shape for the holes (common)
-    TGeoBBox *shHole = 0;
-    shHole = new TGeoBBox(Form("ITSSPD%sGfoilHole", type),0.5*holeLength,
-                          0.5*holeWidth, thickness);
-    
+    new TGeoBBox(Form("ITSSPD%sGfoilHole", type),0.5*holeLength,
+                       0.5*holeWidth, thickness);
+
     // insert the holes in the XTRU shape:
-    // starting from the first value of X, they are simply 
+    // starting from the first value of X, they are simply
     // shifted along this axis
     char name[200];
     TGeoTranslation *transHole[11];
@@ -2056,13 +1766,13 @@ TGeoCompositeShape* AliITSv11GeometrySPD::CreateGroundingFoilShape
             holeX += holeSepX2;
         } // end if else if's
         //cout << i << " --> X = " << holeX << endl;
-        sprintf(name,"ITSSPDTRgFoil%sHole%d", type, i);
+        snprintf(name,200,"ITSSPDTRgFoil%sHole%d", type, i);
         transHole[i] = new TGeoTranslation(name, holeX, holeY, 0.0);
         transHole[i]->RegisterYourself();
         strComposite.Append(Form("ITSSPD%sGfoilHole:%s", type, name));
         if (i < 10) strComposite.Append("+"); else strComposite.Append(")");
     } // end for i
-    
+
     // create composite shape
     TGeoCompositeShape *shGround = new TGeoCompositeShape(
         Form("ITSSPDSHgFoil%s", type), strComposite.Data());
@@ -2070,12 +1780,12 @@ TGeoCompositeShape* AliITSv11GeometrySPD::CreateGroundingFoilShape
     return shGround;
 }
 //______________________________________________________________________
-TGeoVolume* AliITSv11GeometrySPD::CreateGroundingFoil(Bool_t isRight,
+TGeoVolumeAssembly* AliITSv11GeometrySPD::CreateGroundingFoil(Bool_t isRight,
                                    TArrayD &sizes, TGeoManager *mgr)
 {
     //
-    // Create a volume containing all parts of the grounding foil a 
-    // for a half-stave. 
+    // Create a volume containing all parts of the grounding foil a
+    // for a half-stave.
     // It consists of 4 layers with the same shape but different thickness:
     // 1) a layer of glue
     // 2) the aluminum layer
@@ -2083,34 +1793,34 @@ TGeoVolume* AliITSv11GeometrySPD::CreateGroundingFoil(Bool_t isRight,
     // 4) another layer of glue
     // ---
     // Arguments:
-    //  1: a boolean value to know if it is the grounding foir for 
+    //  1: a boolean value to know if it is the grounding foir for
     //     the right or left side
     //  2: a TArrayD which will contain the dimension of the container box:
     //       - size[0] = length along Z (the beam line direction)
-    //       - size[1] = the 'width' of the stave, which defines, together 
+    //       - size[1] = the 'width' of the stave, which defines, together
     //                   with Z, the plane of the carbon fiber support
-    //       - size[2] = 'thickness' (= the direction along which all 
+    //       - size[2] = 'thickness' (= the direction along which all
     //                    stave components are superimposed)
     //  3: the TGeoManager
     // ---
-    // The return value is a TGeoBBox volume containing all grounding 
+    // The return value is a TGeoBBox volume containing all grounding
     // foil components.
     // to avoid strange behaviour of the geometry manager,
     // create a suffix to be used in the names of all shapes
     //
     char suf[5];
-    if (isRight) strcpy(suf, "R"); else strcpy(suf, "L");
-    // this volume will be created in order to ease its placement in 
-    // the half-stave; then, it is added here the small distance of 
-    // the "central" edge of each volume from the Z=0 plane in the stave 
+    if (isRight) strncpy(suf, "R", 5); else strncpy(suf, "L", 5);
+    // this volume will be created in order to ease its placement in
+    // the half-stave; then, it is added here the small distance of
+    // the "central" edge of each volume from the Z=0 plane in the stave
     // reference (which coincides with ALICE one)
     Double_t dist = fgkmm * 0.71;
-    
+
     // define materials
     TGeoMedium *medKap  = GetMedium("SPD KAPTON(POLYCH2)$", mgr);
     TGeoMedium *medAlu  = GetMedium("AL$", mgr);
     TGeoMedium *medGlue = GetMedium("EPOXY$", mgr); //??? GLUE_GF_SUPPORT
-    
+
     // compute the volume shapes (thicknesses change from one to the other)
     Double_t kpLength, kpWidth, alLength, alWidth;
     TArrayD  kpSize, alSize, glSize;
@@ -2129,8 +1839,8 @@ TGeoVolume* AliITSv11GeometrySPD::CreateGroundingFoil(Bool_t isRight,
                                                           g0Thickness, glSize);
     TGeoCompositeShape *g1Shape = CreateGroundingFoilShape(3,kpLength,kpWidth,
                                                           g1Thickness, glSize);
-    // create the component volumes and register their sizes in the 
-    // passed arrays for readability reasons, some reference variables 
+    // create the component volumes and register their sizes in the
+    // passed arrays for readability reasons, some reference variables
     // explicit the meaning of the array slots
     TGeoVolume *kpVol = new TGeoVolume(Form("ITSSPDgFoilKap%s",suf),
                                        kpShape, medKap);
@@ -2150,22 +1860,23 @@ TGeoVolume* AliITSv11GeometrySPD::CreateGroundingFoil(Bool_t isRight,
     Double_t &fullThickness = sizes[0];
     Double_t &fullLength = sizes[1];
     Double_t &fullWidth = sizes[2];
-    // kapton leads the larger dimensions of the foil 
+    // kapton leads the larger dimensions of the foil
     // (including the cited small distance from Z=0 stave reference plane)
     // the thickness is the sum of the ones of all components
     fullLength    = kpLength + dist;
     fullWidth     = kpWidth;
     fullThickness = kpThickness + alThickness + g0Thickness + g1Thickness;
     // create the container
-    TGeoMedium *air = GetMedium("AIR$", mgr);
-    TGeoVolume *container = mgr->MakeBox(Form("ITSSPDgFOIL-%s",suf),
-                 air, 0.5*fullThickness, 0.5*fullWidth, 0.5*fullLength);
-    // create the common correction rotation (which depends of what side 
+//    TGeoMedium *air = GetMedium("AIR$", mgr);
+    TGeoVolumeAssembly *container = new TGeoVolumeAssembly(Form("ITSSPDgFOIL-%s",suf));
+//    TGeoVolume *container = mgr->MakeBox(Form("ITSSPDgFOIL-%s",suf),
+//                 air, 0.5*fullThickness, 0.5*fullWidth, 0.5*fullLength);
+    // create the common correction rotation (which depends of what side
     // we are building)
     TGeoRotation *rotCorr = new TGeoRotation(*gGeoIdentity);
     if (isRight) rotCorr->RotateY(90.0);
-    else rotCorr->RotateY(-90.0);        
-    // compute the translations, which are in the length and 
+    else rotCorr->RotateY(-90.0);
+    // compute the translations, which are in the length and
     // thickness directions
     Double_t x, y, z, shift = 0.0;
     if (isRight) shift = dist;
@@ -2184,7 +1895,7 @@ TGeoVolume* AliITSv11GeometrySPD::CreateGroundingFoil(Bool_t isRight,
     x += 0.5*(alThickness + g1Thickness);
     z  = 0.5*(fullLength - kpLength) - shift;
     TGeoCombiTrans *glTrans1 = new TGeoCombiTrans(x, 0.0, z, rotCorr);
-    
+
     //cout << fgkGapHalfStave << endl;
     //cout << g0Thickness << endl;
     //cout << kpThickness << endl;
@@ -2196,8 +1907,8 @@ TGeoVolume* AliITSv11GeometrySPD::CreateGroundingFoil(Bool_t isRight,
     container->AddNode(kpVol, 1, kpTrans);
     container->AddNode(alVol, 1, alTrans);
     container->AddNode(g0Vol, 1, glTrans0);
-    container->AddNode(g1Vol, 2, glTrans1);    
-    // to add the grease we remember the sizes of the holes, stored as 
+    container->AddNode(g1Vol, 2, glTrans1);
+    // to add the grease we remember the sizes of the holes, stored as
     // additional parameters in the kapton layer size:
     //   - sizes[3] = hole length
     //   - sizes[4] = hole width
@@ -2205,7 +1916,7 @@ TGeoVolume* AliITSv11GeometrySPD::CreateGroundingFoil(Bool_t isRight,
     //   - sizes[6] = standard separation between holes
     //   - sizes[7] = separation between 5th and 6th hole
     //   - sizes[8] = separation between 10th and 11th hole
-    //   - sizes[9] = separation between the upper hole border and 
+    //   - sizes[9] = separation between the upper hole border and
     //                the foil border
     Double_t holeLength      = kpSize[0];
     Double_t holeWidth       = kpSize[1];
@@ -2253,37 +1964,37 @@ TGeoVolumeAssembly* AliITSv11GeometrySPD::CreateMCM(Bool_t isRight,
     //  - the MCM chips (specifications from EDMS)
     //  - the cap which covers the zone where chips are bound to MCM
     // ---
-    // The local reference frame of this assembly is defined in such a way 
-    // that all volumes are contained in a virtual box whose center 
-    // is placed exactly in the middle of the occupied space w.r to all 
-    // directions. This will ease the positioning of this object in the 
-    // half-stave. The sizes of this virtual box are stored in 
+    // The local reference frame of this assembly is defined in such a way
+    // that all volumes are contained in a virtual box whose center
+    // is placed exactly in the middle of the occupied space w.r to all
+    // directions. This will ease the positioning of this object in the
+    // half-stave. The sizes of this virtual box are stored in
     // the array passed by reference.
     // ---
     // Arguments:
-    //  - a boolean flag to know if this is the "left" or "right" MCM, when 
-    //    looking at the stave from above (i.e. the direction from which 
-    //    one sees bus over ladders over grounding foil) and keeping the 
-    //    continuous border in the upper part, one sees the thicker part 
+    //  - a boolean flag to know if this is the "left" or "right" MCM, when
+    //    looking at the stave from above (i.e. the direction from which
+    //    one sees bus over ladders over grounding foil) and keeping the
+    //    continuous border in the upper part, one sees the thicker part
     //    on the left or right.
-    //  - an array passed by reference which will contain the size of 
+    //  - an array passed by reference which will contain the size of
     //    the virtual container.
     //  - a pointer to the used TGeoManager.
     //
 
     // to distinguish the "left" and "right" objects, a suffix is created
     char suf[5];
-    if (isRight) strcpy(suf, "R"); else strcpy(suf, "L");
+    if (isRight) strncpy(suf, "R", 5); else strncpy(suf, "L", 5);
 
     // ** MEDIA **
     TGeoMedium *medBase = GetMedium("SPD KAPTON(POLYCH2)$",mgr);// ??? MCM BASE
     TGeoMedium *medChip = GetMedium("SPD SI CHIP$",mgr);
     TGeoMedium *medCap  = GetMedium("AL$",mgr);
 
-    // The shape of the MCM is divided into 3 sectors with different 
+    // The shape of the MCM is divided into 3 sectors with different
     // widths (Y) and lengths (X), like in this sketch:
     //
-    //   0                      1                                   2 
+    //   0                      1                                   2
     //    +---------------------+-----------------------------------+
     //    |                                    4       sect 2       |
     //    |                    6      sect 1    /-------------------+
@@ -2294,7 +2005,7 @@ TGeoVolumeAssembly* AliITSv11GeometrySPD::CreateMCM(Bool_t isRight,
     // the inclination of all oblique borders (6-7, 4-5) is always 45 degrees.
     // From drawings we can parametrize the dimensions of all these sectors,
     // then the shape of this part of the MCM is implemented as a
-    // TGeoXtru centerd in the virtual XY space. 
+    // TGeoXtru centerd in the virtual XY space.
     // The first step is definig the relevant sizes of this shape:
     Int_t i, j;
     Double_t mcmThickness  = fgkmm * 0.35;
@@ -2331,9 +2042,9 @@ TGeoVolumeAssembly* AliITSv11GeometrySPD::CreateMCM(Bool_t isRight,
     width = sizeYsector[0];
     thickness = mcmThickness + capHeight;
 
-    // define all the relevant vertices of the polygon 
+    // define all the relevant vertices of the polygon
     // which defines the transverse shape of the MCM.
-    // These values are used to several purposes, and 
+    // These values are used to several purposes, and
     // for each one, some points must be excluded
     Double_t xRef[9], yRef[9];
     xRef[0] = -0.5*sizeXtot;
@@ -2355,12 +2066,12 @@ TGeoVolumeAssembly* AliITSv11GeometrySPD::CreateMCM(Bool_t isRight,
     xRef[8] =  xRef[0];
     yRef[8] = -yRef[0];
 
-    // the above points are defined for the "right" MCM (if ve view the 
-    // stave from above) in order to change to the "left" one, we must 
+    // the above points are defined for the "right" MCM (if ve view the
+    // stave from above) in order to change to the "left" one, we must
     // change the sign to all X values:
     if (isRight) for (i = 0; i < 9; i++) xRef[i] = -xRef[i];
-    
-    // the shape of the MCM and glue layer are done excluding point 1, 
+
+    // the shape of the MCM and glue layer are done excluding point 1,
     // which is not necessary and cause the geometry builder to get confused
     j = 0;
     Double_t xBase[8], yBase[8];
@@ -2380,7 +2091,7 @@ TGeoVolumeAssembly* AliITSv11GeometrySPD::CreateMCM(Bool_t isRight,
         j++;
     } // end for i
 
-    // define positions of chips, 
+    // define positions of chips,
     // which must be added to the bottom-left corner of MCM
     // and divided by 1E4;
     Double_t chipX[5], chipY[5];
@@ -2422,7 +2133,7 @@ TGeoVolumeAssembly* AliITSv11GeometrySPD::CreateMCM(Bool_t isRight,
         chipThickness[i] *= fgkmm;
     } // end for i
 
-    // create shapes for MCM 
+    // create shapes for MCM
     Double_t z1, z2;
     TGeoXtru *shBase = new TGeoXtru(2);
     z1 = -0.5*thickness;
@@ -2435,8 +2146,8 @@ TGeoVolumeAssembly* AliITSv11GeometrySPD::CreateMCM(Bool_t isRight,
     TGeoVolume *volBase = new TGeoVolume("ITSSPDbase", shBase, medBase);
     volBase->SetLineColor(kRed);
 
-    // to create the border of the MCM cover, it is required the 
-    // subtraction of two shapes the outer is created using the 
+    // to create the border of the MCM cover, it is required the
+    // subtraction of two shapes the outer is created using the
     // reference points defined here
     TGeoXtru *shCapOut = new TGeoXtru(2);
     shCapOut->SetName(Form("ITSSPDshCAPOUT%s", suf));
@@ -2486,15 +2197,15 @@ TGeoVolumeAssembly* AliITSv11GeometrySPD::CreateMCM(Bool_t isRight,
     shCapIn->DefineSection(1, z2 + 0.01, 0., 0., 1.0);
     // compose shapes
     TGeoCompositeShape *shCapBorder = new TGeoCompositeShape(
-                            Form("ITSSPDshBORDER%s", suf), 
+                            Form("ITSSPDshBORDER%s", suf),
                             Form("%s-%s", shCapOut->GetName(),
                                  shCapIn->GetName()));
     // create volume
     TGeoVolume *volCapBorder = new TGeoVolume("ITSSPDcapBoarder",
                                               shCapBorder,medCap);
     volCapBorder->SetLineColor(kGreen);
-    // finally, we create the top of the cover, which has the same 
-    // shape of outer border and a thickness equal of the one othe 
+    // finally, we create the top of the cover, which has the same
+    // shape of outer border and a thickness equal of the one othe
     // cover border one
     TGeoXtru *shCapTop = new TGeoXtru(2);
     z1 = z2;
@@ -2523,236 +2234,26 @@ TGeoVolumeAssembly* AliITSv11GeometrySPD::CreateMCM(Bool_t isRight,
     // add cap border
     mcmAssembly->AddNode(volCapBorder, 1, gGeoIdentity);
     // add cap top
-    mcmAssembly->AddNode(volCapTop, 1, gGeoIdentity);    
+    mcmAssembly->AddNode(volCapTop, 1, gGeoIdentity);
 
     return mcmAssembly;
 }
 
-/*
-//__________________________________________________________________________________________
-TGeoVolumeAssembly* AliITSv11GeometrySPD::CreatePixelBus
-(Bool_t isRight, TArrayD &sizes, TGeoManager *mgr) const
-{
-    //
-    // The pixel bus is implemented as a TGeoBBox with some objects on it, 
-    // which could affect the particle energy loss.
-    // ---
-    // In order to avoid confusion, the bus is directly displaced 
-    // according to the axis orientations which are used in the final stave:
-    // X --> thickness direction
-    // Y --> width direction
-    // Z --> length direction
-    //
-  
-    
-    // ** MEDIA **
-    
-    //PIXEL BUS
-    TGeoMedium *medBus     = GetMedium("SPDBUS(AL+KPT+EPOX)$",mgr);
-    TGeoMedium *medPt1000  = GetMedium("CERAMICS$",mgr); // ??? PT1000
-    // Capacity
-    TGeoMedium *medCap     = GetMedium("SDD X7R capacitors$",mgr);
-    // ??? Resistance
-    // TGeoMedium *medRes     = GetMedium("SDD X7R capacitors$",mgr);
-    TGeoMedium *medRes     = GetMedium("ALUMINUM$",mgr);
-    TGeoMedium *medExt     = GetMedium("SDDKAPTON (POLYCH2)$", mgr);
-    // ** SIZES & POSITIONS **
-    Double_t busLength          = 170.501 * fgkmm; // length of plane part
-    Double_t busWidth           =  13.800 * fgkmm; // width
-    Double_t busThickness       =   0.280 * fgkmm; // thickness
-    Double_t pt1000Length       = fgkmm * 1.50;
-    Double_t pt1000Width        = fgkmm * 3.10;
-    Double_t pt1000Thickness    = fgkmm * 0.60;
-    Double_t pt1000Y, pt1000Z[10];// position of the pt1000's along the bus
-    Double_t capLength          = fgkmm * 2.55;
-    Double_t capWidth           = fgkmm * 1.50;
-    Double_t capThickness       = fgkmm * 1.35;
-    Double_t capY[2], capZ[2];
-    
-    Double_t resLength          = fgkmm * 2.20;
-    Double_t resWidth           = fgkmm * 0.80;
-    Double_t resThickness       = fgkmm * 0.35;
-    Double_t resY[2], resZ[2];
-    
-    Double_t extThickness       = fgkmm * 0.25;
-    Double_t ext1Length         = fgkmm * (26.7 - 10.0);
-    Double_t ext2Length         = fgkmm * (285.0 - ext1Length + extThickness);
-    Double_t extWidth           = fgkmm * 11.0;
-    Double_t extHeight          = fgkmm * 2.5;
-    
-            
-    // position of pt1000, resistors and capacitors depends on the 
-    // bus if it's left or right one
-    if (!isRight) {
-        pt1000Y    =   64400.;
-        pt1000Z[0] =   66160.;
-        pt1000Z[1] =  206200.;
-        pt1000Z[2] =  346200.;
-        pt1000Z[3] =  486200.;
-        pt1000Z[4] =  626200.;
-        pt1000Z[5] =  776200.;
-        pt1000Z[6] =  916200.;
-        pt1000Z[7] = 1056200.;
-        pt1000Z[8] = 1196200.;
-        pt1000Z[9] = 1336200.;    
-        resZ[0]    = 1397500.;
-        resY[0]    =   26900.;
-        resZ[1]    =  682500.;
-        resY[1]    =   27800.;
-        capZ[0]    = 1395700.;
-        capY[0]    =   45700.;
-        capZ[1]    =  692600.;
-        capY[1]    =   45400.;
-    } else {
-        pt1000Y    =   66100.;
-        pt1000Z[0] =  319700.;
-        pt1000Z[1] =  459700.;
-        pt1000Z[2] =  599700.;
-        pt1000Z[3] =  739700.;
-        pt1000Z[4] =  879700.;
-        pt1000Z[5] = 1029700.;
-        pt1000Z[6] = 1169700.;
-        pt1000Z[7] = 1309700.;
-        pt1000Z[8] = 1449700.;
-        pt1000Z[9] = 1589700.;    
-        capY[0]    =   44500.;
-        capZ[0]    =  266700.;
-        capY[1]    =   44300.;
-        capZ[1]    =  974700.;
-        resZ[0]    =  266500.;
-        resY[0]    =   29200.;
-        resZ[1]    =  974600.;
-        resY[1]    =   29900.;
-    } // end if isRight
-    Int_t i;
-    pt1000Y *= 1E-4 * fgkmm;
-    for (i = 0; i < 10; i++) {
-        pt1000Z[i] *= 1E-4 * fgkmm;
-        if (i < 2) {
-            capZ[i] *= 1E-4 * fgkmm;
-            capY[i] *= 1E-4 * fgkmm;
-            resZ[i] *= 1E-4 * fgkmm;
-            resY[i] *= 1E-4 * fgkmm;
-        }  // end if iM2
-    } // end for i
-    
-    Double_t &fullLength = sizes[1];
-    Double_t &fullWidth = sizes[2];
-    Double_t &fullThickness = sizes[0];
-    fullLength = busLength;
-    fullWidth = busWidth;
-    // add the thickness of the thickest component on bus (capacity)
-    fullThickness = busThickness + capThickness; 
-    // ** VOLUMES **
-    TGeoVolumeAssembly *container = new TGeoVolumeAssembly("PixelBus");
-    TGeoVolume *bus = mgr->MakeBox("Bus", medBus, 0.5*busThickness, 0.5*busWidth, 0.5*busLength);
-    TGeoVolume *pt1000 = mgr->MakeBox("PT1000", medPt1000, 0.5*pt1000Thickness, 0.5*pt1000Width, 0.5*pt1000Length);
-    TGeoVolume *res = mgr->MakeBox("Resistor", medRes, 0.5*resThickness, 0.5*resWidth, 0.5*resLength);
-    TGeoVolume *cap = mgr->MakeBox("Capacitor", medCap, 0.5*capThickness, 0.5*capWidth, 0.5*capLength);
-    TGeoVolume *ext1 = mgr->MakeBox("Extender1", medExt, 0.5*extThickness, 0.5*extWidth, 0.5*ext1Length);
-    TGeoVolume *ext2 = mgr->MakeBox("Extender2", medExt, 0.5*extHeight - extThickness, 0.5*extWidth, 0.5*extThickness);
-    TGeoVolume *ext3 = mgr->MakeBox("Extender3", medExt, extThickness, 0.5*extWidth, 0.5*ext2Length);
-    bus->SetLineColor(kYellow + 2);
-    pt1000->SetLineColor(kGreen + 3);
-    res->SetLineColor(kRed + 1);
-    cap->SetLineColor(kBlue - 7);
-    ext1->SetLineColor(kGray);
-    ext2->SetLineColor(kGray);
-    ext3->SetLineColor(kGray);
-    
-    // ** MOVEMENTS AND POSITIONEMENT **
-    // bus
-    TGeoTranslation *trBus = new TGeoTranslation(0.5 * (busThickness - 
-                                                        fullThickness), 0.0, 0.0);
-    container->AddNode(bus, 0, trBus);
-    Double_t zRef, yRef, x, y, z;
-    if (isRight) {
-        zRef = -0.5*fullLength;
-        yRef = -0.5*fullWidth;
-    } else {
-        zRef = -0.5*fullLength;
-        yRef = -0.5*fullWidth;
-    } // end if isRight
-    // pt1000
-    x = 0.5*(pt1000Thickness - fullThickness) + busThickness;
-    for (i = 0; i < 10; i++) {
-        y = yRef + pt1000Y;
-        z = zRef + pt1000Z[i];
-        TGeoTranslation *tr = new TGeoTranslation(x, y, z);
-        container->AddNode(pt1000, i, tr);
-    } // end for i
-    // capacitors
-    x = 0.5*(capThickness - fullThickness) + busThickness;
-    for (i = 0; i < 2; i++) {
-        y = yRef + capY[i];
-        z = zRef + capZ[i];
-        TGeoTranslation *tr = new TGeoTranslation(x, y, z);
-        container->AddNode(cap, i, tr);
-    } // end for i
-    // resistors
-    x = 0.5*(resThickness - fullThickness) + busThickness;
-    for (i = 0; i < 2; i++) {
-        y = yRef + resY[i];
-        z = zRef + resZ[i];
-        TGeoTranslation *tr = new TGeoTranslation(x, y, z);
-        container->AddNode(res, i, tr);
-    } // end for i
-    // extender
-    if (isRight) {
-        y = 0.5 * (-fullWidth + extWidth);
-        z = 0.5 * (-fullLength + fgkmm * 10.0);
-    }
-    else {
-        y = 0.5 * (fullWidth - extWidth);
-        z = 0.5 * ( fullLength - fgkmm * 10.0);
-    }
-    x = 0.5 * (extThickness - fullThickness) + busThickness;
-    //y = 0.5 * (fullWidth - extWidth);
-    TGeoTranslation *trExt1 = new TGeoTranslation(x, y, z);
-    if (isRight) {
-        z -= 0.5 * (ext1Length - extThickness);
-    }
-    else {
-        z += 0.5 * (ext1Length - extThickness);
-    }
-    x += 0.5*(extHeight - extThickness);
-    TGeoTranslation *trExt2 = new TGeoTranslation(x, y, z);
-    if (isRight) {
-        z -= 0.5 * (ext2Length - extThickness);
-    }
-    else {
-        z += 0.5 * (ext2Length - extThickness);
-    }
-    x += 0.5*(extHeight - extThickness) + extThickness;
-    TGeoTranslation *trExt3 = new TGeoTranslation(x, y, z);
-    container->AddNode(ext1, 0, trExt1);
-    container->AddNode(ext2, 0, trExt2);
-    container->AddNode(ext3, 0, trExt3);
-    
-    
-    sizes[3] = yRef + pt1000Y;
-    sizes[4] = zRef + pt1000Z[2];
-    sizes[5] = zRef + pt1000Z[7];
-    
-    return container;
-}
-*/
-
 //______________________________________________________________________
 TGeoVolumeAssembly* AliITSv11GeometrySPD::CreatePixelBus
 (Bool_t isRight, Int_t ilayer, TArrayD &sizes, TGeoManager *mgr) const
 {
     //
-    // The pixel bus is implemented as a TGeoBBox with some objects on it, 
+    // The pixel bus is implemented as a TGeoBBox with some objects on it,
     // which could affect the particle energy loss.
     // ---
-    // In order to avoid confusion, the bus is directly displaced 
+    // In order to avoid confusion, the bus is directly displaced
     // according to the axis orientations which are used in the final stave:
     // X --> thickness direction
     // Y --> width direction
     // Z --> length direction
     //
-    
+
     // ** CRITICAL CHECK ******************************************************
     // layer number can be ONLY 1 or 2
     if (ilayer != 1 && ilayer != 2) AliFatal("Layer number MUST be 1 or 2");
@@ -2764,9 +2265,10 @@ TGeoVolumeAssembly* AliITSv11GeometrySPD::CreatePixelBus
     // Capacity
     TGeoMedium *medCap     = GetMedium("SDD X7R capacitors$",mgr);
     // ??? Resistance
-    //TGeoMedium *medRes     = GetMedium("SDD X7R capacitors$",mgr); 
+    //TGeoMedium *medRes     = GetMedium("SDD X7R capacitors$",mgr);
     TGeoMedium *medRes     = GetMedium("ALUMINUM$",mgr);
-    TGeoMedium *medExt     = GetMedium("SDDKAPTON (POLYCH2)$", mgr);
+    //TGeoMedium *medExt     = GetMedium("SDDKAPTON (POLYCH2)$", mgr);
+    TGeoMedium *medExt     = GetMedium("SPD-MIX CU KAPTON$", mgr);
     // ** SIZES & POSITIONS **
     Double_t busLength          = 170.501 * fgkmm; // length of plane part
     Double_t busWidth           =  13.800 * fgkmm; // width
@@ -2779,19 +2281,19 @@ TGeoVolumeAssembly* AliITSv11GeometrySPD::CreatePixelBus
     Double_t capWidth           = fgkmm * 1.50;
     Double_t capThickness       = fgkmm * 1.35;
     Double_t capY[2], capZ[2];
-     
+
     Double_t resLength          = fgkmm * 2.20;
     Double_t resWidth           = fgkmm * 0.80;
     Double_t resThickness       = fgkmm * 0.35;
     Double_t resY[2], resZ[2];
-    
+
     Double_t extThickness       = fgkmm * 0.25;
     Double_t ext1Length         = fgkmm * (26.7 - 10.0);
-    Double_t ext2Length         = fgkmm * (284.0 - ext1Length + extThickness);
+    Double_t ext2Length         = fgkmm * 284.0 - ext1Length + extThickness;
     Double_t extWidth           = fgkmm * 11.0;
     Double_t extHeight          = fgkmm * 2.5;
-               
-    // position of pt1000, resistors and capacitors depends on the 
+
+    // position of pt1000, resistors and capacitors depends on the
     // bus if it's left or right one
     if (!isRight) {
         pt1000Y    =   64400.;
@@ -2804,7 +2306,7 @@ TGeoVolumeAssembly* AliITSv11GeometrySPD::CreatePixelBus
         pt1000Z[6] =  916200.;
         pt1000Z[7] = 1056200.;
         pt1000Z[8] = 1196200.;
-        pt1000Z[9] = 1336200.;   
+        pt1000Z[9] = 1336200.;
         resZ[0]    = 1397500.;
         resY[0]    =   26900.;
         resZ[1]    =  682500.;
@@ -2824,7 +2326,7 @@ TGeoVolumeAssembly* AliITSv11GeometrySPD::CreatePixelBus
         pt1000Z[6] = 1169700.;
         pt1000Z[7] = 1309700.;
         pt1000Z[8] = 1449700.;
-        pt1000Z[9] = 1589700.;   
+        pt1000Z[9] = 1589700.;
         capY[0]    =   44500.;
         capZ[0]    =  266700.;
         capY[1]    =   44300.;
@@ -2845,18 +2347,18 @@ TGeoVolumeAssembly* AliITSv11GeometrySPD::CreatePixelBus
             resY[i] *= 1E-4 * fgkmm;
         }  // end if iM2
     } // end for i
-     
+
     Double_t &fullLength = sizes[1];
     Double_t &fullWidth = sizes[2];
     Double_t &fullThickness = sizes[0];
     fullLength = busLength;
     fullWidth = busWidth;
     // add the thickness of the thickest component on bus (capacity)
-    fullThickness = busThickness + capThickness; 
+    fullThickness = busThickness + capThickness;
 
     // ** VOLUMES **
     TGeoVolumeAssembly *container = new TGeoVolumeAssembly("ITSSPDpixelBus");
-    TGeoVolume *bus = mgr->MakeBox("ITSSPDbus", medBus, 0.5*busThickness, 
+    TGeoVolume *bus = mgr->MakeBox("ITSSPDbus", medBus, 0.5*busThickness,
                                    0.5*busWidth, 0.5*busLength);
     TGeoVolume *pt1000 = mgr->MakeBox("ITSSPDpt1000",medPt1000,
                         0.5*pt1000Thickness,0.5*pt1000Width, 0.5*pt1000Length);
@@ -2864,10 +2366,37 @@ TGeoVolumeAssembly* AliITSv11GeometrySPD::CreatePixelBus
                                    0.5*resWidth, 0.5*resLength);
     TGeoVolume *cap = mgr->MakeBox("ITSSPDcapacitor", medCap, 0.5*capThickness,
                                    0.5*capWidth, 0.5*capLength);
-                                   
-    TGeoVolume *ext1 = mgr->MakeBox("Extender1", medExt, 0.5*extThickness, 0.5*extWidth, 0.5*ext1Length);
-    TGeoVolume *ext2 = mgr->MakeBox("Extender2", medExt, 0.5*extHeight - 2.*extThickness, 0.5*extWidth, 0.5*extThickness);
-    TGeoVolume *ext3 = mgr->MakeBox("Extender3", medExt, 0.5*extThickness, 0.5*(extWidth-0.8*fgkmm), 0.5*ext2Length + extThickness); // Hardcode fix of a small overlap
+
+    char extname[12];
+    snprintf(extname,12,"Extender1l%d",ilayer);
+    TGeoVolume *ext1 = mgr->MakeBox(extname, medExt, 0.5*extThickness, 0.5*extWidth, 0.5*ext1Length);
+    snprintf(extname,12,"Extender2l%d",ilayer);
+    TGeoVolume *ext2 = mgr->MakeBox(extname, medExt, 0.5*extHeight - 2.*extThickness, 0.5*extWidth, 0.5*extThickness);
+    TGeoVolume *ext3=0;
+    snprintf(extname,12,"Extender3l%d",ilayer);
+    if (ilayer==1) {
+      Double_t halflen=(0.5*ext2Length + extThickness);
+      Double_t xprof[6],yprof[6];
+      Double_t alpha=24;
+      xprof[0] = -halflen;
+      yprof[0] = -0.5*extThickness;
+      xprof[1] = halflen/2;
+      yprof[1] = yprof[0];
+      xprof[2] = xprof[1] + 0.5*halflen*CosD(alpha);
+      yprof[2] = yprof[1] + 0.5*halflen*SinD(alpha);
+      xprof[3] = xprof[2] - extThickness*SinD(alpha);
+      yprof[3] = yprof[2] + extThickness*CosD(alpha);
+      InsidePoint(xprof[0], yprof[0], xprof[1], yprof[1], xprof[2], yprof[2],
+                 extThickness, xprof[4], yprof[4]);
+      xprof[5] = xprof[0];
+      yprof[5] = 0.5*extThickness;
+      TGeoXtru *ext3sh = new TGeoXtru(2);
+      ext3sh->DefinePolygon(6, xprof, yprof);
+      ext3sh->DefineSection(0, -0.5*(extWidth-0.8*fgkmm));
+      ext3sh->DefineSection(1,  0.5*(extWidth-0.8*fgkmm));
+      ext3 = new TGeoVolume(extname, ext3sh, medExt);
+    } else
+      ext3 = mgr->MakeBox(extname, medExt, 0.5*extThickness, 0.5*(extWidth-0.8*fgkmm), 0.5*ext2Length + extThickness); // Hardcode fix of a small overlap
     bus->SetLineColor(kYellow + 2);
     pt1000->SetLineColor(kGreen + 3);
     res->SetLineColor(kRed + 1);
@@ -2878,7 +2407,7 @@ TGeoVolumeAssembly* AliITSv11GeometrySPD::CreatePixelBus
 
     // ** MOVEMENTS AND POSITIONEMENT **
     // bus
-    TGeoTranslation *trBus = new TGeoTranslation(0.5 * (busThickness - 
+    TGeoTranslation *trBus = new TGeoTranslation(0.5 * (busThickness -
                                                    fullThickness), 0.0, 0.0);
     container->AddNode(bus, 1, trBus);
     Double_t zRef, yRef, x, y, z;
@@ -2913,7 +2442,7 @@ TGeoVolumeAssembly* AliITSv11GeometrySPD::CreatePixelBus
         TGeoTranslation *tr = new TGeoTranslation(x, y, z);
         container->AddNode(res, i+1, tr);
     } // end for i
-    
+
     // extender
         if (ilayer == 2) {
        if (isRight) {
@@ -2953,139 +2482,536 @@ TGeoVolumeAssembly* AliITSv11GeometrySPD::CreatePixelBus
         z += 0.5 * (ext2Length - extThickness) + 2.5*extThickness;
     }
     x += 0.5*(extHeight - extThickness) - 2.*extThickness;
-    TGeoTranslation *trExt3 = new TGeoTranslation(x, y, z);
+    TGeoCombiTrans *trExt3=0;
+    if (ilayer==1) {
+      if (isRight)
+       trExt3 = new TGeoCombiTrans(x, y, z, new TGeoRotation("",0.,-90.,90.));
+      else
+       trExt3 = new TGeoCombiTrans(x, y, z, new TGeoRotation("",0., 90.,90.));
+    } else
+      trExt3 = new TGeoCombiTrans(x, y, z, 0);
     container->AddNode(ext1, 0, trExt1);
     container->AddNode(ext2, 0, trExt2);
     container->AddNode(ext3, 0, trExt3);
-    
+
     sizes[3] = yRef + pt1000Y;
     sizes[4] = zRef + pt1000Z[2];
     sizes[5] = zRef + pt1000Z[7];
-    
+
     return container;
 }
 
 //______________________________________________________________________
-TGeoVolumeAssembly* AliITSv11GeometrySPD::CreateConeModule(TGeoManager *mgr) const
+TList* AliITSv11GeometrySPD::CreateConeModule(Bool_t sideC, const Double_t angrot,
+                                             TGeoManager *mgr) const
 {
+    //
+    // Creates all services modules and places them in a TList
+    // angrot is the rotation angle (passed as an argument to avoid
+    // defining the same quantity in two different places)
+    //
+    // Created:      ?? ??? 2008  A. Pulvirenti
+    // Updated:      03 May 2010  M. Sitta
+    // Updated:      20 Jun 2010  A. Pulvirenti  Optical patch panels
+    // Updated:      22 Jun 2010  M. Sitta  Fiber cables
+    // Updated:      04 Jul 2010  M. Sitta  Water cooling
+    // Updated:      08 Jul 2010  A. Pulvirenti  Air cooling on Side C
+    //
+
     TGeoMedium *medInox  = GetMedium("INOX$",mgr);
-    TGeoMedium *medExt   = GetMedium("SDDKAPTON (POLYCH2)$", mgr);
+    //TGeoMedium *medExt   = GetMedium("SDDKAPTON (POLYCH2)$", mgr);
+    TGeoMedium *medExtB  = GetMedium("SPD-BUS CU KAPTON$", mgr);
+    TGeoMedium *medExtM  = GetMedium("SPD-MCM CU KAPTON$", mgr);
     TGeoMedium *medPlate = GetMedium("SPD C (M55J)$", mgr);
-    
+    TGeoMedium *medFreon = GetMedium("Freon$", mgr);
+    TGeoMedium *medGas   = GetMedium("GASEOUS FREON$", mgr);
+    TGeoMedium *medFibs  = GetMedium("SDD OPTICFIB$",mgr);
+    TGeoMedium *medCopper= GetMedium("COPPER$",mgr);
+    TGeoMedium *medPVC   = GetMedium("PVC$",mgr);
+
     Double_t extThickness = fgkmm * 0.25;
     Double_t ext1Length   = fgkmm * (26.7 - 10.0);
-    Double_t ext2Length   = fgkmm * (285.0 - ext1Length + extThickness);
-    
-    Double_t cableThickness = 1.5 * fgkmm;
-    Double_t cableL1 = 350.0 * fgkmm - extThickness - ext1Length - ext2Length;
-    Double_t cableL2 = 426.0 * fgkmm;
+//    Double_t ext2Length   = fgkmm * (285.0 - ext1Length + extThickness);
+    Double_t ext2Length   = fgkmm * 285.0 - ext1Length + extThickness;
+
+    const Double_t kCableThickness  =   1.5  *fgkmm;
+    Double_t cableL0 =  10.0 * fgkmm;
+    Double_t cableL1 = 340.0 * fgkmm - extThickness - ext1Length - ext2Length;
+    Double_t cableL2 = 300.0 * fgkmm;
     //Double_t cableL3 = 570.0 * fgkmm;
     Double_t cableL3 = 57.0 * fgkmm;
     Double_t cableW1 =  11.0 * fgkmm;
     Double_t cableW2 =  30.0 * fgkmm;
     Double_t cableW3 =  50.0 * fgkmm;
-    
-    Double_t mcmThickness = 1.2 *fgkmm;
-    Double_t mcmLength = cableL1 + cableL2 + cableL3;
-    Double_t mcmWidth = cableW1;
-    
-    Double_t plateLength    = 200.0 * fgkmm;
-    Double_t plateWidth     =  50.0 * fgkmm;
-    Double_t plateThickness =   5.0 * fgkmm;
-    
+
+    const Double_t kMCMLength       =   cableL0 + cableL1 + cableL2 + cableL3;
+    const Double_t kMCMWidth        =   cableW1;
+    const Double_t kMCMThickness    =   1.2  *fgkmm;
+
+    const Double_t kPlateLength     = 200.0  *fgkmm;
+    const Double_t kPlateWidth      =  50.0  *fgkmm;
+    const Double_t kPlateThickness  =   5.0  *fgkmm;
+
+    const Double_t kConeTubeRmin    =   2.0  *fgkmm;
+    const Double_t kConeTubeRmax    =   3.0  *fgkmm;
+
+    const Double_t kHorizTubeLen    = 150.0  *fgkmm;
+    const Double_t kYtoHalfStave    =   9.5  *fgkmm;
+
+    const Double_t kWaterCoolRMax   =   2.6  *fgkmm;
+    const Double_t kWaterCoolThick  =   0.04 *fgkmm;
+    const Double_t kWaterCoolLen    = 250.0  *fgkmm;
+    const Double_t kWCPlateThick    =   0.5  *fgkmm;
+    const Double_t kWCPlateWide     =  33.0  *fgkmm;
+    const Double_t kWCPlateLen      = 230.0  *fgkmm;
+    const Double_t kWCFittingRext1  =   2.4  *fgkmm;
+    const Double_t kWCFittingRext2  =   3.7  *fgkmm;
+    const Double_t kWCFittingRint1  =   1.9  *fgkmm;
+    const Double_t kWCFittingRint2  = kWaterCoolRMax;
+    const Double_t kWCFittingLen1   =   7.0  *fgkmm;
+    const Double_t kWCFittingLen2   =   8.0  *fgkmm;
+    
+    const Double_t kCollWidth       =  40.0  *fgkmm;
+    const Double_t kCollLength      =  60.0  *fgkmm;
+    const Double_t kCollThickness   =  10.0  *fgkmm;
+    const Double_t kCollTubeThick   =   1.0  *fgkmm;
+    const Double_t kCollTubeRadius  =   7.0  *fgkmm;
+    const Double_t kCollTubeLength  = 190.0  *fgkmm;
+
+    const Double_t kOptFibDiamet    =   4.5  *fgkmm;
+
     Double_t x[12], y[12];
-    
-    x[0] = 7.5;
+    Double_t xloc, yloc, zloc;
+
+    Int_t kPurple = 6; // Purple (Root does not define it)
+
+    TGeoVolumeAssembly* container[5];
+    container[0] = new TGeoVolumeAssembly("ITSSPDConeModule");
+    container[1] = new TGeoVolumeAssembly("ITSSPDCoolingModuleSideA");
+    container[2] = new TGeoVolumeAssembly("ITSSPDCoolingModuleSideC");
+    container[3] = new TGeoVolumeAssembly("ITSSPDPatchPanelModule");
+    container[4] = new TGeoVolumeAssembly("ITSSPDWaterCooling");
+
+    // The extender on the cone as a Xtru
+    x[0] = -cableL0;
     y[0] = 0.0 + 0.5 * cableW1;
-    
-    x[1] = x[0] + cableL1 - 0.5*(cableW2 - cableW1);
+
+    x[1] = x[0] + cableL0 + cableL1 - 0.5*(cableW2 - cableW1);
     y[1] = y[0];
-    
-    x[2] = x[0] + cableL1; 
+
+    x[2] = x[0] + cableL0 + cableL1;
     y[2] = y[1] + 0.5*(cableW2 - cableW1);
-    
+
     x[3] = x[2] + cableL2;
     y[3] = y[2];
-    
+
     x[4] = x[3] + 0.5*(cableW3 - cableW2);
     y[4] = y[3] + 0.5*(cableW3 - cableW2);
-    
+
     x[5] = x[4] + cableL3 - 0.5*(cableW3 - cableW2);
     y[5] = y[4];
-    
+
     for (Int_t i = 6; i < 12; i++) {
         x[i] =  x[11 - i];
         y[i] = -y[11 - i];
     }
-    
-    TGeoVolumeAssembly *container = new TGeoVolumeAssembly("ITSSPDConeModule");
-    
+
     TGeoXtru *shCable = new TGeoXtru(2);
     shCable->DefinePolygon(12, x, y);
-    shCable->DefineSection(0, 0., 0., 0., 1.0);
-    shCable->DefineSection(1, cableThickness, 0., 0., 1.0);
-    
-    TGeoVolume *volCable = new TGeoVolume("ITSSPDExtender", shCable, medExt);
+    shCable->DefineSection(0, 0.0);
+    shCable->DefineSection(1, kCableThickness);
+
+    TGeoVolume *volCable = new TGeoVolume("ITSSPDExtender", shCable, medExtB);
     volCable->SetLineColor(kGreen);
-    
-    TGeoVolume *volTube = gGeoManager->MakeTube("ITSSPDCoolingTubeCone", medInox, 5.*fgkmm, 6.*fgkmm, 0.5*(x[5] - x[0]));
-    volTube->SetLineColor(kGray);
-    
-    Double_t thickness = cableThickness + mcmThickness;
-    TGeoBBox *shOut = new TGeoBBox("ITSSPD_shape_plateout", 0.5*plateThickness, 0.5*plateLength, 0.5*plateWidth);
-    TGeoBBox *shIn = new TGeoBBox("ITSSPD_shape_platein", 0.5*thickness, 0.52*plateLength, 0.5*cableW2);
+
+    // The MCM extender on the cone as a Xtru
+    TGeoBBox *shMCMExt = new TGeoBBox(0.5*kMCMLength,
+                                     0.5*kMCMWidth,
+                                     0.5*kMCMThickness);
+
+    TGeoVolume *volMCMExt = new TGeoVolume("ITSSPDExtenderMCM",
+                                          shMCMExt, medExtM);
+    volMCMExt->SetLineColor(kGreen+3);
+
+    // The support plate on the cone as a composite shape
+    Double_t thickness = kCableThickness + kMCMThickness;
+    TGeoBBox *shOut = new TGeoBBox("ITSSPD_shape_plateout",
+                                  0.5*kPlateLength,
+                                  0.5*kPlateWidth,
+                                  0.5*kPlateThickness);
+    TGeoBBox *shIn  = new TGeoBBox("ITSSPD_shape_platein" ,
+                                  0.5*kPlateLength,
+                                  0.5*cableW2,
+                                  0.5*thickness);
     Char_t string[255];
-    sprintf(string, "%s-%s", shOut->GetName(), shIn->GetName());
-    TGeoCompositeShape *shPlate = new TGeoCompositeShape("ITSSPDPlate_shape", string);
-    TGeoVolume *volPlate = new TGeoVolume("ITSSPDPlate", shPlate, medPlate);
+    snprintf(string, 255, "%s-%s", shOut->GetName(), shIn->GetName());
+    TGeoCompositeShape *shPlate = new TGeoCompositeShape("ITSSPDPlate_shape",
+                                string);
+
+    TGeoVolume *volPlate = new TGeoVolume("ITSSPDPlate",
+                                         shPlate, medPlate);
     volPlate->SetLineColor(kRed);
     
-    TGeoVolume *volMCMExt = gGeoManager->MakeBox("ITSSPDextenderMCM", medExt, 0.5*mcmThickness, 0.5*mcmLength, 0.5*mcmWidth);
-    volMCMExt->SetLineColor(kGreen+3);
-    
-    TGeoRotation *rot = new TGeoRotation(*gGeoIdentity);
-    rot->RotateX(90.0);
-    rot->RotateZ(90.0);
-    container->AddNode(volCable, 0, rot);
-    
-    TGeoTranslation *combi = new TGeoTranslation(cableThickness + 0.5*mcmThickness, x[0] + 0.5*mcmLength, 0.0);
-    container->AddNode(volMCMExt, 0, combi);
-    
-    TGeoRotation *rot1 = new TGeoRotation(*gGeoIdentity);
-    rot1->RotateX(88.5);
-    TGeoCombiTrans *tr = new TGeoCombiTrans(1.0, x[0] + 0.5*(x[5] - x[0]), -2.95, rot1);
-    container->AddNode(volTube, 0, tr);
-    
-    TGeoTranslation *tr1 = new TGeoTranslation(0.5*plateThickness - 0.5*(plateThickness-thickness), x[3] - x[0] - 0.52*plateLength, 0.0);
-    container->AddNode(volPlate, 0, tr1);
-    
-    return container;
+    // The air cooling tubes
+    TGeoBBox   *shCollBox   = new TGeoBBox("ITSSPD_shape_collector_box", 0.5*kCollLength, 0.5*kCollWidth, 0.5*kCollThickness);
+    TGeoTube   *shCollTube  = new TGeoTube("ITSSPD_shape_collector_tube",kCollTubeRadius - kCollTubeThick, kCollTubeRadius, 0.5*kCollTubeLength);
+    TGeoVolume *volCollBox  = new TGeoVolume("ITSSPDCollectorBox", shCollBox, medPVC);
+    TGeoVolume *volCollTube = new TGeoVolume("ITSSPDCollectorTube", shCollTube, medPVC);
+    volCollBox->SetLineColor(kAzure);
+    volCollTube->SetLineColor(kAzure);
+
+    // The cooling tube on the cone as a Ctub
+    Double_t tubeLength = shCable->GetX(5) - shCable->GetX(0) + kYtoHalfStave;
+    TGeoCtub *shTube = new TGeoCtub(0, kConeTubeRmax, 0.5*tubeLength, 0, 360,
+                                   0, SinD(angrot/2), -CosD(angrot/2),
+                                   0,              0,              1);
+
+    TGeoVolume *volTubeA = new TGeoVolume("ITSSPDCoolingTubeOnConeA",
+                                         shTube, medInox);
+    volTubeA->SetLineColor(kGray);
+
+    TGeoVolume *volTubeC = new TGeoVolume("ITSSPDCoolingTubeOnConeC",
+                                         shTube, medInox);
+    volTubeC->SetLineColor(kGray);
+
+    // The freon in the cooling tubes on the cone as a Ctub
+    TGeoCtub *shFreon = new TGeoCtub(0, kConeTubeRmin, 0.5*tubeLength, 0, 360,
+                                    0, SinD(angrot/2), -CosD(angrot/2),
+                                    0,              0,              1);
+
+    TGeoVolume *volFreon = new TGeoVolume("ITSSPDCoolingFreonOnCone",
+                                         shFreon, medFreon);
+    volFreon->SetLineColor(kPurple);
+
+    TGeoVolume *volGasFr = new TGeoVolume("ITSSPDCoolingFreonGasOnCone",
+                                         shFreon, medGas);
+    volGasFr->SetLineColor(kPurple);
+
+    // The cooling tube inside the cylinder as a Ctub
+    TGeoCtub *shCylTub = new TGeoCtub(0, kConeTubeRmax,
+                                     0.5*kHorizTubeLen, 0, 360,
+                                     0,            0,           -1,
+                                     0, SinD(angrot/2), CosD(angrot/2));
+
+    TGeoVolume *volCylTubA = new TGeoVolume("ITSSPDCoolingTubeOnCylA",
+                                           shCylTub, medInox);
+    volCylTubA->SetLineColor(kGray);
+
+    TGeoVolume *volCylTubC = new TGeoVolume("ITSSPDCoolingTubeOnCylC",
+                                           shCylTub, medInox);
+    volCylTubC->SetLineColor(kGray);
+
+    // The freon in the cooling tubes in the cylinder as a Ctub
+    TGeoCtub *shCylFr = new TGeoCtub(0, kConeTubeRmin,
+                                    0.5*kHorizTubeLen, 0, 360,
+                                    0,            0,           -1,
+                                    0, SinD(angrot/2), CosD(angrot/2));
+
+    TGeoVolume *volCylFr = new TGeoVolume("ITSSPDCoolingFreonOnCyl",
+                                         shCylFr, medFreon);
+    volCylFr->SetLineColor(kPurple);
+
+    TGeoVolume *volCylGasFr = new TGeoVolume("ITSSPDCoolingFreonGasOnCyl",
+                                            shCylFr, medGas);
+    volCylGasFr->SetLineColor(kPurple);
+
+    // The optical fibers bundle on the cone as a Tube
+    Double_t optLength = shCable->GetX(5) - shCable->GetX(0) + kYtoHalfStave;
+    TGeoTube *shOptFibs = new TGeoTube(0., 0.5*kOptFibDiamet, 0.5*optLength);
+
+    TGeoVolume *volOptFibs = new TGeoVolume("ITSSPDOpticalFibersOnCone",
+                                           shOptFibs, medFibs);
+    volOptFibs->SetLineColor(kOrange);
+
+    // The optical patch panels
+    TArrayD psizes;
+    TGeoVolume *volPatch = CreatePatchPanel(psizes, mgr);
+
+    // The water cooling tube as a Tube
+    TGeoTube *shWatCool = new TGeoTube(kWaterCoolRMax-kWaterCoolThick,
+                                      kWaterCoolRMax, kWaterCoolLen/2);
+
+    TGeoVolume *volWatCool = new TGeoVolume("ITSSPDWaterCoolingOnCone",
+                                           shWatCool, medInox);
+    volWatCool->SetLineColor(kGray);
+
+    // The support plate for the water tubes: a Tubs and a BBox
+    TGeoTubeSeg *shWCPltT = new TGeoTubeSeg(kWaterCoolRMax,
+                                           kWaterCoolRMax+kWCPlateThick,
+                                           kWCPlateLen/2, 180., 360.);
+
+    Double_t plateBoxWide = (kWCPlateWide - 2*kWaterCoolRMax)/2;
+    TGeoBBox *shWCPltB = new TGeoBBox(plateBoxWide/2,
+                                     kWCPlateThick/2,
+                                     kWCPlateLen/2);
+
+    TGeoVolume *volWCPltT = new TGeoVolume("ITSSPDWaterCoolingTubsPlate",
+                                         shWCPltT, medPlate);
+    volWCPltT->SetLineColor(kRed);
+
+    TGeoVolume *volWCPltB = new TGeoVolume("ITSSPDWaterCoolingBoxPlate",
+                                         shWCPltB, medPlate);
+    volWCPltB->SetLineColor(kRed);
+
+    // The fitting for the water cooling tube: a Pcon
+    TGeoPcon *shFitt = new TGeoPcon(0., 360., 4);
+    shFitt->Z(0)    = -kWCFittingLen1;
+    shFitt->Rmin(0) =  kWCFittingRint1;
+    shFitt->Rmax(0) =  kWCFittingRext1;
+
+    shFitt->Z(1)    =  0;
+    shFitt->Rmin(1) =  kWCFittingRint1;
+    shFitt->Rmax(1) =  kWCFittingRext1;
+
+    shFitt->Z(2)    =  0;
+    shFitt->Rmin(2) =  kWCFittingRint2;
+    shFitt->Rmax(2) =  kWCFittingRext2;
+
+    shFitt->Z(3)    =  kWCFittingLen2;
+    shFitt->Rmin(3) =  kWCFittingRint2;
+    shFitt->Rmax(3) =  kWCFittingRext2;
+
+    TGeoVolume *volFitt = new TGeoVolume("ITSSPDWaterCoolingFitting",
+                                        shFitt, medCopper);
+    volFitt->SetLineColor(kOrange);
+
+    // Now place everything in the containers
+    volTubeA->AddNode(volGasFr, 1, 0);
+    volTubeC->AddNode(volFreon, 1, 0);
+
+    volCylTubA->AddNode(volCylGasFr, 1, 0);
+    volCylTubC->AddNode(volCylFr   , 1, 0);
+
+    container[0]->AddNode(volCable, 1, 0);
+
+    xloc = shMCMExt->GetDX() - cableL0;
+    zloc = shMCMExt->GetDZ();
+    container[0]->AddNode(volMCMExt, 1,
+                         new TGeoTranslation( xloc, 0.,-zloc));
+
+    xloc = shMCMExt->GetDX();
+    zloc = shCable->GetZ(1)/2 - shMCMExt->GetDZ();
+    container[0]->AddNode(volPlate, 1,
+                         new TGeoTranslation( xloc, 0., zloc));
+
+    TGeoRotation *rot2 = new TGeoRotation(*gGeoIdentity);
+    rot2->SetName("rotPatch");
+    rot2->RotateX(90.0);
+    rot2->RotateY(163.0);
+    //rot2->RotateZ(132.5);
+    
+    // add collectors only on side C
+    if (sideC)
+    {
+      TGeoTranslation *trCollBox   = new TGeoTranslation(xloc - 0.5*kPlateLength + 0.5*kCollLength, 0.0, +0.5*(kPlateThickness+1.1*kCollThickness));
+      TGeoRotation    *rotCollTube = new TGeoRotation(*gGeoIdentity);
+      rotCollTube->RotateY(90.0);
+      TGeoCombiTrans  *trCollTube  = new TGeoCombiTrans(xloc + 0.5*kCollTubeLength - (0.5*kPlateLength - kCollLength), 0.0, +0.5*(kPlateThickness+2.0*kCollTubeRadius+kCollTubeThick), rotCollTube);
+      container[0]->AddNode(volCollBox, 1, trCollBox);
+      container[0]->AddNode(volCollTube, 1, trCollTube);
+    }
+        
+    Double_t dxPatch = 2.9;
+    Double_t dzPatch = 2.8;
+    TGeoCombiTrans *tr2 = new TGeoCombiTrans(1.7*ext2Length - dxPatch, 0.0, dzPatch, rot2);
+    container[3]->AddNode(volPatch, 0, tr2);
+
+    xloc = shTube->GetRmax();
+    yloc = shTube->GetRmax();
+    zloc = shTube->GetDz() - shTube->GetRmax() - kYtoHalfStave;
+    container[1]->AddNode(volTubeA, 1,
+                         new TGeoTranslation(-xloc, -yloc, zloc));
+    container[2]->AddNode(volTubeC, 1,
+                         new TGeoTranslation(-xloc, -yloc, zloc));
+
+    xloc = shTube->GetRmax();
+    yloc = (shCylTub->GetDz())*SinD(angrot) - shTube->GetRmax();
+    zloc = (shCylTub->GetDz())*CosD(angrot) + shTube->GetRmax() +kYtoHalfStave;
+    container[1]->AddNode(volCylTubA, 1,
+                         new TGeoCombiTrans(-xloc, yloc,-zloc,
+                                    new TGeoRotation("",0.,angrot,0.)));
+    container[2]->AddNode(volCylTubC, 1,
+                         new TGeoCombiTrans(-xloc, yloc,-zloc,
+                                    new TGeoRotation("",0.,angrot,0.)));
+
+    xloc = shOptFibs->GetRmax() + 2*shTube->GetRmax();
+    yloc = 1.6*shOptFibs->GetRmax();
+    zloc = shOptFibs->GetDZ() - shTube->GetRmax() - kYtoHalfStave;
+    container[1]->AddNode(volOptFibs, 1,
+                         new TGeoTranslation(-xloc, -yloc, zloc));
+    container[2]->AddNode(volOptFibs, 1,
+                         new TGeoTranslation(-xloc, -yloc, zloc));
+
+    yloc = shWatCool->GetRmax();
+    zloc = (2*shTube->GetDz() - shTube->GetRmax() - kYtoHalfStave)/2;
+    container[4]->AddNode(volWatCool, 1,
+                         new TGeoTranslation(0, -yloc, zloc));
+
+    container[4]->AddNode(volWCPltT, 1,
+                         new TGeoTranslation(0, -yloc, zloc));
+
+    yloc -= shWCPltB->GetDY();
+    xloc = shWatCool->GetRmax() + shWCPltB->GetDX();
+    container[4]->AddNode(volWCPltB, 1,
+                         new TGeoTranslation( xloc, -yloc, zloc));
+    container[4]->AddNode(volWCPltB, 2,
+                         new TGeoTranslation(-xloc, -yloc, zloc));
+
+    yloc = shWatCool->GetRmax();
+    zloc -= shWatCool->GetDz();
+    container[4]->AddNode(volFitt, 1,
+                         new TGeoTranslation(0, -yloc, zloc));
+
+    // Finally create the list of assemblies and return it to the caller
+    TList* conemodulelist = new TList();
+    conemodulelist->Add(container[0]);
+    conemodulelist->Add(container[1]);
+    conemodulelist->Add(container[2]);
+    conemodulelist->Add(container[3]);
+    conemodulelist->Add(container[4]);
+
+    return conemodulelist;
 }
 
 //______________________________________________________________________
 void AliITSv11GeometrySPD::CreateCones(TGeoVolume *moth) const
 {
-    
-    TGeoVolumeAssembly *module = CreateConeModule(gGeoManager);
-         
+    //
+    // Places all services modules in the mother reference system
+    //
+    // Created:      ?? ??? 2008  Alberto Pulvirenti
+    // Updated:      03 May 2010  Mario Sitta
+    // Updated:      04 Jul 2010  Mario Sitta  Water cooling
+    //
+
+    const Int_t kNumberOfModules    =  10;
+
+    const Double_t kInnerRadius     =  80.775*fgkmm;
+    const Double_t kZTrans          = 452.000*fgkmm;
+    const Double_t kAlphaRot        =  46.500*fgkDegree;
+    const Double_t kAlphaSpaceCool  =   9.200*fgkDegree;
+
+    TList*  modulelistA = CreateConeModule(kFALSE, 90-kAlphaRot);
+    TList*  modulelistC = CreateConeModule(kTRUE , 90-kAlphaRot);
+    TList* &modulelist  = modulelistC;
+    TGeoVolumeAssembly* module, *moduleA, *moduleC;
+
+    Double_t xloc, yloc, zloc;
+
     //Double_t angle[10] = {18., 54., 90., 126., 162., -18., -54., -90., -126., -162.};
-    Double_t angle[10] = {18., 54., 90., 126., 162., 198.0, 234.0, 270.0, 306.0, 342.0};
-    for (Int_t i = 0; i < 10; i++) {
+    // anglem for cone modules (cables and cooling tubes)
+    // anglep for pathc panels
+    Double_t anglem[10] = {18., 54., 90., 126., 162., 198., 234., 270., 306., 342.};
+    Double_t anglep[10] = {18., 62., 90., 115., 162., 198., 242., 270., 295., 342.};
+//    Double_t angle1m[10] = {23., 53., 90., 127., 157., 203.0, 233.0, 270.0, 307.0, 337.0};
+//    Double_t angle2m[10] = {18., 53., 90., 126., 162., 198.0, 233.0, 270.0, 309.0, 342.0};
+//    Double_t angle1c[10] = {23., 53., 90., 124., 157., 203.0, 233.0, 270.0, 304.0, 337.0};
+//    Double_t angle2c[10] = {18., 44., 90., 126., 162., 198.0, 223.0, 270.0, 309.0, 342.0};
+
+    // First add the cables
+    moduleA = (TGeoVolumeAssembly*)modulelistA->At(0);
+    moduleC = (TGeoVolumeAssembly*)modulelistC->At(0);
+    for (Int_t i = 0; i < kNumberOfModules; i++) {
+        TGeoRotation *rot1 = new TGeoRotation(*gGeoIdentity);
+       rot1->RotateY(-kAlphaRot);
+       rot1->RotateZ(anglem[i]);
+        xloc = kInnerRadius*CosD(anglem[i]);
+        yloc = kInnerRadius*SinD(anglem[i]);
+       zloc = kZTrans;
+        moth->AddNode(moduleA, 2*i+2,
+                     new TGeoCombiTrans( xloc, yloc, zloc, rot1));
+
+        TGeoRotation *rot2 = new TGeoRotation(*gGeoIdentity);
+       rot2->RotateY(180.-kAlphaRot);
+       rot2->RotateZ(anglem[i]);
+        xloc = kInnerRadius*CosD(anglem[i]);
+        yloc = kInnerRadius*SinD(anglem[i]);
+       zloc = kZTrans;
+        moth->AddNode(moduleC, 2*i+1,
+                     new TGeoCombiTrans(-xloc,-yloc,-zloc, rot2));
+    }
+
+    // Then the cooling tubes on Side A
+    module = (TGeoVolumeAssembly*)modulelist->At(1);
+    Double_t anglec;
+    for (Int_t i = 0; i < kNumberOfModules; i++) {
+        anglec = anglem[i] + kAlphaSpaceCool;
+        TGeoRotation *rot1 = new TGeoRotation(*gGeoIdentity);
+        rot1->RotateX(-90.0+kAlphaRot-0.04); // 0.04 fixes small overlap
+       rot1->RotateZ(-90.0+anglec);
+        xloc = kInnerRadius*CosD(anglec);
+        yloc = kInnerRadius*SinD(anglec);
+       zloc = kZTrans+0.162; // 0.162 fixes small overlap
+        moth->AddNode(module, 2*i+2, 
+                     new TGeoCombiTrans( xloc, yloc, zloc, rot1));
+    }
+
+    // And the cooling tubes on Side C
+    module = (TGeoVolumeAssembly*)modulelist->At(2);
+    for (Int_t i = 0; i < kNumberOfModules; i++) {
+        anglec = anglem[i] - kAlphaSpaceCool;
+        TGeoRotation *rot2 = new TGeoRotation(*gGeoIdentity);
+        rot2->RotateX(-90.0+kAlphaRot-0.04); // 0.04 fixes small overlap
+       rot2->RotateY(180.);
+       rot2->RotateZ(90.0+anglec);
+        xloc = kInnerRadius*CosD(anglec);
+        yloc = kInnerRadius*SinD(anglec);
+       zloc = kZTrans+0.162; // 0.162 fixes small overlap
+        moth->AddNode(module, 2*i+1,
+                     new TGeoCombiTrans(-xloc,-yloc,-zloc, rot2));
+    }
+
+    // Then the water cooling tubes
+    module = (TGeoVolumeAssembly*)modulelist->At(4);
+    for (Int_t i = 1; i < kNumberOfModules; i++) { // i = 1,2,...,9
+        if (i != 5) { // There is no tube in this position
+         anglec = (anglem[i-1]+anglem[i])/2;
+           TGeoRotation *rot1 = new TGeoRotation(*gGeoIdentity);
+           rot1->RotateX(-90.0+kAlphaRot);
+           rot1->RotateZ(-90.0+anglec);
+           xloc = kInnerRadius*CosD(anglec);
+           yloc = kInnerRadius*SinD(anglec);
+           zloc = kZTrans;
+           moth->AddNode(module, 2*i+2,
+                         new TGeoCombiTrans( xloc, yloc, zloc, rot1));
+
+           TGeoRotation *rot2 = new TGeoRotation(*gGeoIdentity);
+           rot2->RotateX(-90.0+kAlphaRot);
+           rot2->RotateY(180.);
+           rot2->RotateZ(90.0+anglec);
+           xloc = kInnerRadius*CosD(anglec);
+           yloc = kInnerRadius*SinD(anglec);
+           zloc = kZTrans;
+           moth->AddNode(module, 2*i+1,
+                         new TGeoCombiTrans(-xloc,-yloc,-zloc, rot2));
+       }
+    }
+
+    // Finally the optical patch panels
+    module = (TGeoVolumeAssembly*)modulelist->At(3);
+    for (Int_t i = 0; i < kNumberOfModules; i++) {
         TGeoRotation *rot1 = new TGeoRotation(*gGeoIdentity);
-        rot1->RotateY(-90.0);
-        rot1->RotateX(44.7);
-        rot1->RotateZ(90.0 - angle[i]);
-        TGeoCombiTrans *tr1 = new TGeoCombiTrans(0.0, 0.0, 40.5, rot1);
-        moth->AddNode(module, 2*i, tr1);
+       rot1->RotateY(-kAlphaRot);
+       rot1->RotateZ(anglep[i]);
+        xloc = kInnerRadius*CosD(anglep[i]);
+        yloc = kInnerRadius*SinD(anglep[i]);
+       zloc = kZTrans;
+        moth->AddNode(module, 2*i+2,
+                     new TGeoCombiTrans( xloc, yloc, zloc, rot1));
+
         TGeoRotation *rot2 = new TGeoRotation(*gGeoIdentity);
-        rot2->RotateY(90.0);
-        rot2->RotateX(-44.7);
-        rot2->RotateZ(90.0 - angle[i]);
-        TGeoCombiTrans *tr2 = new TGeoCombiTrans(0.0, 0.0, -40.5, rot2);
-        moth->AddNode(module, 2*i+1, tr2);
+       rot2->RotateY(180.-kAlphaRot);
+       rot2->RotateZ(anglep[i]);
+        xloc = kInnerRadius*CosD(anglep[i]);
+        yloc = kInnerRadius*SinD(anglep[i]);
+       zloc = kZTrans;
+        moth->AddNode(module, 2*i+1,
+                     new TGeoCombiTrans(-xloc,-yloc,-zloc, rot2));
     }
+
 }
 
+
 //______________________________________________________________________
 TGeoVolume* AliITSv11GeometrySPD::CreateExtender(
     const Double_t *extenderParams, const TGeoMedium *extenderMedium,
@@ -3096,8 +3022,8 @@ TGeoVolume* AliITSv11GeometrySPD::CreateExtender(
     //
     // This function creates the following picture (in plane xOy)
     // Should be useful for the definition of the pixel bus and MCM extenders
-    // The origin corresponds to point 0 on the picture, at half-width 
-    // in Z direction 
+    // The origin corresponds to point 0 on the picture, at half-width
+    // in Z direction
     //
     //   Y                         7     6                      5
     //   ^                           +---+---------------------+
@@ -3124,21 +3050,21 @@ TGeoVolume* AliITSv11GeometrySPD::CreateExtender(
     //   |--> par 4 : outer length [3-4] / [6-5]
     //   |--> par 5 : width in local Z direction
     //
-    Double_t slopeDeltaX = (extenderParams[3] - extenderParams[1] 
-                            * TMath::Cos(extenderParams[2])) / 
+    Double_t slopeDeltaX = (extenderParams[3] - extenderParams[1]
+                            * TMath::Cos(extenderParams[2])) /
                             TMath::Tan(extenderParams[2]);
     Double_t extenderXtruX[10] = {
         0 ,
         extenderParams[0] ,
-        extenderParams[0]+extenderParams[1]*TMath::Sin(extenderParams[2]) , 
+        extenderParams[0]+extenderParams[1]*TMath::Sin(extenderParams[2]) ,
         extenderParams[0]+extenderParams[1]*TMath::Sin(extenderParams[2])+
                                                               slopeDeltaX ,
         extenderParams[0]+extenderParams[1]*TMath::Sin(extenderParams[2])+
-                                           slopeDeltaX + extenderParams[4], 
+                                           slopeDeltaX + extenderParams[4],
         extenderParams[0]+extenderParams[1]*TMath::Sin(extenderParams[2])+
-                                           slopeDeltaX + extenderParams[4], 
+                                           slopeDeltaX + extenderParams[4],
         extenderParams[0]+extenderParams[1]*TMath::Sin(extenderParams[2])+
-                                                              slopeDeltaX , 
+                                                              slopeDeltaX ,
         extenderParams[0]+extenderParams[1]*TMath::Sin(extenderParams[2])+
           slopeDeltaX - extenderParams[1] * TMath::Sin(extenderParams[2]) ,
         extenderParams[0] ,
@@ -3176,203 +3102,16 @@ TGeoVolume* AliITSv11GeometrySPD::CreateExtender(
     extenderXtru->DefineSection(1, 0.5*extenderParams[4]);
     return extenderXtruVol;
 }
-//______________________________________________________________________
-TGeoVolumeAssembly* AliITSv11GeometrySPD::CreatePixelBusAndExtensions
-(Bool_t /*zpos*/, TGeoManager *mgr) const
-{
-    //
-    // Creates an assembly which contains the pixel bus and its extension
-    // and the extension of the MCM.
-    // By: Renaud Vernet
-    // NOTE: to be defined its material and its extension in the outside 
-    // direction
-    //
-    // ====   constants   =====
-    //get the media
-    // PIXEL BUS
-    //TGeoMedium   *medPixelBus    = GetMedium("SPDBUS(AL+KPT+EPOX)$",mgr);
-    // IXEL BUS EXTENDER
-    TGeoMedium *medPBExtender  = GetMedium("SDDKAPTON (POLYCH2)$",mgr);
-    //MCM EXTENDER
-    TGeoMedium *medMCMExtender = GetMedium("SDDKAPTON (POLYCH2)$",mgr);
-    //   //geometrical constants
-    const Double_t kPbextenderThickness     =   0.07 * fgkmm;
-    //design=?? 70 deg. seems OK
-    const Double_t kPbExtenderSlopeAngle    =  70.0  * TMath::Pi()/180.;
-    // = 2.6 - (0.28+0.05+0.35) cf design
-    const Double_t kPbExtenderHeight        =   1.92 * fgkmm;
-    const Double_t kPbExtenderWidthY        =  11.0  * fgkmm;
-    //design=?? 70 deg. seems OK
-    const Double_t kMcmExtenderSlopeAngle   =  70.0  * TMath::Pi()/180.;
-    const Double_t kMcmExtenderThickness    =   0.10 * fgkmm;
-    const Double_t kMcmExtenderHeight       =   1.8  * fgkmm;
-    const Double_t kMcmExtenderWidthY       =   kPbExtenderWidthY;
-    //   const Double_t groundingThickness    =   0.07  * fgkmm;
-    //   const Double_t grounding2pixelBusDz  =   0.625 * fgkmm;
-    //   const Double_t pixelBusThickness     =   0.28  * fgkmm;
-    //   const Double_t groundingWidthX       = 170.501 * fgkmm;
-    //   const Double_t pixelBusContactDx     =   1.099 * fgkmm;
-    //   const Double_t pixelBusWidthY        =  13.8   * fgkmm;
-    //design=20 deg.
-    //   const Double_t pixelBusContactPhi    =  20.0   * TMath::Pi()/180.
-    //   const Double_t pbExtenderTopZ        =   2.72  * fgkmm;
-    //   const Double_t mcmThickness          =   0.35  * fgkmm;
-    //   const Double_t halfStaveTotalLength  = 247.64  * fgkmm;
-    //   const Double_t deltaYOrigin          =  15.95/2.* fgkmm;
-    //   const Double_t deltaXOrigin          =   1.1    * fgkmm;
-    //   const Double_t deltaZOrigin          = halfStaveTotalLength / 2.;
-    //   const Double_t grounding2pixelBusDz2 = grounding2pixelBusDz+
-    //                           groundingThickness/2. + pixelBusThickness/2.;
-    //   const Double_t pixelBusWidthX        = groundingWidthX;
-    //   const Double_t pixelBusRaiseLength   = (pixelBusContactDx-
-    //                  pixelBusThickness*TMath::Sin(pixelBusContactPhi))/
-    //                                       TMath::Cos(pixelBusContactPhi);
-    //   const Double_t pbExtenderBaseZ       = grounding2pixelBusDz2 +
-    //        pixelBusRaiseLength*TMath::Sin(pixelBusContactPhi) +
-    //        2*pixelBusThickness*TMath::Sin(pixelBusContactPhi)*
-    //        TMath::Tan(pixelBusContactPhi);
-    //   const Double_t pbExtenderDeltaZ      = pbExtenderTopZ-pbExtenderBaseZ;
-    //   const Double_t pbExtenderEndPointX   = 2*deltaZOrigin - 
-    //    groundingWidthX - 2*pixelBusThickness*TMath::Sin(pixelBusContactPhi);
-    //   const Double_t pbExtenderXtru3L   = 1.5 * fgkmm; //arbitrary ?
-    //   const Double_t pbExtenderXtru4L   = (pbExtenderDeltaZ + 
-    //             pixelBusThickness*(TMath::Cos(extenderSlope)-2))/
-    //                                      TMath::Sin(extenderSlope);
-    //   const Double_t kMcmExtenderEndPointX  = deltaZOrigin - 48.2 * fgkmm;
-    //   const Double_t kMcmExtenderXtru3L     = 1.5  * fgkmm;
-    //   //=====  end constants  =====
-    const Double_t kPbExtenderInnerLength    = 10. * fgkmm;
-    const Double_t kPbExtenderOuterLength    = 15. * fgkmm;
-    const Double_t kMcmExtenderInnerLength   = 10. * fgkmm;
-    const Double_t kMcmExtenderOuterLength   = 15. * fgkmm;
-    Double_t pbExtenderParams[6]  = {kPbExtenderInnerLength,  //0
-                                     kPbextenderThickness,    //1
-                                     kPbExtenderSlopeAngle,   //2
-                                     kPbExtenderHeight,       //3
-                                     kPbExtenderOuterLength,  //4
-                                     kPbExtenderWidthY};      //5
-    
-    Double_t mcmExtenderParams[6] = {kMcmExtenderInnerLength, //0
-                                     kMcmExtenderThickness,   //1
-                                     kMcmExtenderSlopeAngle,  //2
-                                     kMcmExtenderHeight,      //3
-                                     kMcmExtenderOuterLength, //4
-                                     kMcmExtenderWidthY};     //5
-    
-    TArrayD sizes(3);
-    TGeoVolume* pbExtender  = CreateExtender(pbExtenderParams,medPBExtender,
-                                             sizes);
-    if(GetDebug(1))printf("CREATED AN EXTENDER : THICKNESS = %5.5f cm\t"
-              "LENGTH=%5.5f cm\tWIDTH=%5.5f cm\n",sizes[0],sizes[1],sizes[2]);
-    TGeoVolume* mcmExtender = CreateExtender(mcmExtenderParams,medMCMExtender,
-                                             sizes);
-    if(GetDebug(1))printf("CREATED AN EXTENDER : THICKNESS = %5.5f cm\t"
-             "LENGTH=%5.5f cm\tWIDTH=%5.5f cm\n",sizes[0],sizes[1],sizes[2]);
-    //   Double_t pixelBusValues[5]    = {pixelBusWidthX,        //0
-    //                     pixelBusThickness,     //1
-    //                     pixelBusContactPhi,    //2
-    //                     pixelBusRaiseLength,   //3
-    //                     pixelBusWidthY};      //4
-    
-    //   Double_t pbExtenderValues[8]  = {pixelBusRaiseLength,   //0
-    //                     pixelBusContactPhi,     //1
-    //                     pbExtenderXtru3L,       //2
-    //                     pixelBusThickness,      //3
-    //                     extenderSlope,     //4
-    //                     pbExtenderXtru4L,      //5
-    //                     pbExtenderEndPointX,   //6
-    //                     kPbExtenderWidthY};    //7
-    
-    //   Double_t mcmExtenderValues[6] = {mcmExtenderXtru3L,     //0
-    //                     mcmExtenderThickness,  //1
-    //                     extenderSlope,     //2
-    //                     deltaMcmMcmExtender,    //3
-    //                     mcmExtenderEndPointX,  //4
-    //                     mcmExtenderWidthY};    //5
-    //   TGeoVolumeAssembly *pixelBus=new TGeoVolumeAssembly("ITSSPDpixelBus");
-    //   CreatePixelBus(pixelBus,pixelBusValues,medPixelBus); 
-    //   TGeoVolumeAssembly *pbExtender = new TGeoVolumeAssembly(
-    //                                              "ITSSPDpixelBusExtender");
-    //   CreatePixelBusExtender(pbExtender,pbExtenderValues,medPBExtender);
-    //   TGeoVolumeAssembly *mcmExtender = new TGeoVolumeAssembly(
-    //                                                 "ITSSPDmcmExtender");
-    //   CreateMCMExtender(mcmExtender,mcmExtenderValues,medMCMExtender);
-    //--------------   DEFINITION OF GEOMETRICAL TRANSFORMATIONS --------
-    //   TGeoRotation    * commonRot  = new TGeoRotation("commonRot",0,90,0);
-    //   commonRot->MultiplyBy(new TGeoRotation("rot",-90,0,0));
-    //   TGeoTranslation * pixelBusTrans   = new TGeoTranslation(
-    //                      pixelBusThickness/2. - deltaXOrigin + 0.52*fgkmm ,
-    //                                   -pixelBusWidthY/2.   + deltaYOrigin , 
-    //                                   -groundingWidthX/2.  + deltaZOrigin);
-    //   TGeoRotation    *pixelBusRot     = new TGeoRotation(*commonRot);
-    //   TGeoTranslation *pbExtenderTrans =new TGeoTranslation(*pixelBusTrans);
-    //   TGeoRotation    *pbExtenderRot   = new TGeoRotation(*pixelBusRot);
-    //   pbExtenderTrans->SetDz(*(pbExtenderTrans->GetTranslation()+2) -
-    //                          pixelBusWidthX/2. - 2*pixelBusThickness*
-    //                                    TMath::Sin(pixelBusContactPhi));  
-    //   if (!zpos) {
-    //     pbExtenderTrans->SetDy(*(pbExtenderTrans->GetTranslation()+1) -
-    //                               (pixelBusWidthY - kPbExtenderWidthY)/2.);
-    //   } else {
-    //     pbExtenderTrans->SetDy(*(pbExtenderTrans->GetTranslation()+1) +
-    //                            (pixelBusWidthY - kPbExtenderWidthY)/2.);
-    //   }
-    //   pbExtenderTrans->SetDx(*(pbExtenderTrans->GetTranslation()) +
-    //                      pixelBusThickness/2 + 2*pixelBusThickness*
-    //                      TMath::Sin(pixelBusContactPhi)*
-    //                      TMath::Tan(pixelBusContactPhi));
-    //   TGeoTranslation * mcmExtenderTrans = new TGeoTranslation(0.12*fgkmm +
-    //                                    mcmThickness - deltaXOrigin,
-    //                                    pbExtenderTrans->GetTranslation()[1],
-    //                                    -4.82);
-    //   TGeoRotation    * mcmExtenderRot   = new TGeoRotation(*pbExtenderRot);
-    //   // add pt1000 components
-    //   Double_t pt1000Z = fgkmm * 64400. * 1E-4;
-    //   //Double_t pt1000X[10] = {319700.,  459700.,  599700.,  739700., 
-    //                             879700., 1029700., 1169700., 1309700.,
-    //                            1449700., 1589700.};
-    //   Double_t pt1000X[10] ={66160., 206200.,  346200.,  486200.,  626200.,
-    //                         776200., 916200., 1056200., 1196200., 1336200.};
-    //   Double_t pt1000size[3] = {fgkmm*1.5, fgkmm*0.6, fgkmm*3.1};
-    //   Int_t i;
-    //   for (i = 0; i < 10; i++) {
-    //     pt1000X[i] *= fgkmm * 1E-4;
-    //   }
-    //   TGeoVolume *pt1000 = mgr->MakeBox("ITSSPDpt1000",0,0.5*pt1000size[0],
-    //                              0.5*pt1000size[1], 0.5*pt1000size[2]);
-    //   pt1000->SetLineColor(kGray);
-    //   Double_t refThickness = - pixelBusThickness;
-    //   for (i = 0; i < 10; i++) {
-    //     TGeoTranslation *tr = new TGeoTranslation(pt1000X[i]-
-    //          0.5*pixelBusWidthX, 0.002+0.5*(-3.*refThickness+pt1000size[3]),
-    //                                            pt1000Z -0.5*pixelBusWidthY);
-    //     pixelBus->AddNode(pt1000, i+1, tr);
-    //   }
-    
-    //CREATE FINAL VOLUME ASSEMBLY AND ROTATE IT
-    TGeoVolumeAssembly *assembly = new TGeoVolumeAssembly("ITSSPDextenders");
-    //   assembly->AddNode((TGeoVolume*)pixelBus,1,
-    //          new TGeoCombiTrans(*pixelBusTrans,*pixelBusRot));
-    //   assembly->AddNode((TGeoVolume*)pbExtender,1,
-    //           new TGeoCombiTrans(*pbExtenderTrans,*pbExtenderRot));
-    //   assembly->AddNode((TGeoVolume*)mcmExtender,1,
-    //         new TGeoCombiTrans(*mcmExtenderTrans,*mcmExtenderRot));
-    //   assembly->AddNode(mcmExtender,1,new TGeoIdentity());
-    assembly->AddNode(pbExtender,1);
-    assembly->AddNode(mcmExtender,1);
-    //   assembly->SetTransparency(50);
-    
-    return assembly;
-}
+
 //______________________________________________________________________
 TGeoVolumeAssembly* AliITSv11GeometrySPD::CreateHalfStave(Bool_t isRight,
 Int_t layer,Int_t idxCentral,Int_t idxSide,TArrayD &sizes,TGeoManager *mgr)
 {
     //
-    // Implementation of an half-stave, which depends on the side where 
-    // we are on the stave. The convention for "left" and "right" is the 
-    // same as for the MCM. The return value is a TGeoAssembly which is 
-    // structured in such a way that the origin of its local reference 
+    // Implementation of an half-stave, which depends on the side where
+    // we are on the stave. The convention for "left" and "right" is the
+    // same as for the MCM. The return value is a TGeoAssembly which is
+    // structured in such a way that the origin of its local reference
     // frame coincides with the origin of the whole stave.
     // The TArrayD passed by reference will contain details of the shape:
     //  - sizes[0] = thickness
@@ -3398,11 +3137,11 @@ Int_t layer,Int_t idxCentral,Int_t idxSide,TArrayD &sizes,TGeoManager *mgr)
 
     // define the separations along Z direction between the objects
     Double_t sepLadderLadder = fgkmm * 0.2; // sep. btw the 2 ladders
-    Double_t sepLadderCenter = fgkmm * 0.4; // sep. btw the "central" ladder 
+    Double_t sepLadderCenter = fgkmm * 0.4; // sep. btw the "central" ladder
                                             // and the Z=0 plane in stave ref.
     Double_t sepLadderMCM    = fgkmm * 0.3; // sep. btw the "external" ladder
                                             // and MCM
-    Double_t sepBusCenter    = fgkmm * 0.3; // sep. btw the bus central edge 
+    Double_t sepBusCenter    = fgkmm * 0.3; // sep. btw the bus central edge
                                             // and the Z=0 plane in stave ref.
 
     // ** VOLUMES **
@@ -3427,7 +3166,7 @@ Int_t layer,Int_t idxCentral,Int_t idxSide,TArrayD &sizes,TGeoManager *mgr)
     Double_t mcmThickness = mcmSize[0];
     Double_t mcmLength = mcmSize[1];
     Double_t mcmWidth = mcmSize[2];
-    
+
     // bus
     TArrayD busSize(6);
     TGeoVolumeAssembly *bus = CreatePixelBus(isRight, layer, busSize, mgr);
@@ -3447,7 +3186,7 @@ Int_t layer,Int_t idxCentral,Int_t idxSide,TArrayD &sizes,TGeoManager *mgr)
     Double_t &fullThickness = sizes[0];
     Double_t &fullLength = sizes[1];
     Double_t &fullWidth = sizes[2];
-    
+
     // compute the full size of the container
     fullLength    = sepLadderCenter+2.0*ladderLength+sepLadderMCM+
                        sepLadderLadder+mcmLength;
@@ -3464,13 +3203,13 @@ Int_t layer,Int_t idxCentral,Int_t idxSide,TArrayD &sizes,TGeoManager *mgr)
     TGeoTranslation *grndTrans = new TGeoTranslation(xGrnd, 0.0, zGrnd);
 
     // ladders (translations along thickness and length)
-    // layers must be sorted going from the one at largest Z to the 
+    // layers must be sorted going from the one at largest Z to the
     // one at smallest Z:
     // -|Zmax| ------> |Zmax|
     //      3   2   1   0
-    // then, for layer 1 ladders they must be placed exactly this way, 
-    // and in layer 2 at the opposite. In order to remember the placements, 
-    // we define as "inner" and "outer" ladder respectively the one close 
+    // then, for layer 1 ladders they must be placed exactly this way,
+    // and in layer 2 at the opposite. In order to remember the placements,
+    // we define as "inner" and "outer" ladder respectively the one close
     // to barrel center, and the one closer to MCM, respectively.
     Double_t xLad, zLadIn, zLadOut;
     xLad    = xGrnd + 0.5*(grndThickness + ladderThickness) +
@@ -3493,15 +3232,15 @@ Int_t layer,Int_t idxCentral,Int_t idxSide,TArrayD &sizes,TGeoManager *mgr)
     TGeoCombiTrans *trLadOut = new TGeoCombiTrans(xLad,ladderShift,zLadOut,
                                                   rotLad);
 
-    // MCM (length and thickness direction, placing at same level as the 
-    // ladder, which implies to recompute the position of center, because 
-    // ladder and MCM have NOT the same thickness) the two copies of the 
+    // MCM (length and thickness direction, placing at same level as the
+    // ladder, which implies to recompute the position of center, because
+    // ladder and MCM have NOT the same thickness) the two copies of the
     // MCM are placed at the same distance from the center, on both sides
-    Double_t xMCM = xGrnd + 0.5*grndThickness + 0.5*mcmThickness + 
+    Double_t xMCM = xGrnd + 0.5*grndThickness + 0.5*mcmThickness +
                     0.01175 - fgkGapLadder;
     Double_t yMCM = 0.5*(fullWidth - mcmWidth);
     Double_t zMCM = zLadOut - 0.5*ladderLength - 0.5*mcmLength - sepLadderMCM;
-    if (!isRight) zMCM = zLadOut + 0.5*ladderLength + 0.5*mcmLength + 
+    if (!isRight) zMCM = zLadOut + 0.5*ladderLength + 0.5*mcmLength +
                          sepLadderMCM;
 
     // create the correction rotations
@@ -3510,7 +3249,7 @@ Int_t layer,Int_t idxCentral,Int_t idxSide,TArrayD &sizes,TGeoManager *mgr)
     TGeoCombiTrans *trMCM = new TGeoCombiTrans(xMCM, yMCM, zMCM, rotMCM);
 
     // glue between ladders and pixel bus
-    Double_t xLadGlue = xLad + 0.5*ladderThickness + 0.01175 - 
+    Double_t xLadGlue = xLad + 0.5*ladderThickness + 0.01175 -
                         fgkGapLadder + 0.5*ladGlueThickness;
 
     // bus (length and thickness direction)
@@ -3568,18 +3307,18 @@ TGeoVolumeAssembly* AliITSv11GeometrySPD::CreateStave(Int_t layer,
     // when put on the sector.
     // This assembly contains, going from bottom to top in the thickness
     // direction:
-    //   - the complete grounding foil, defined by the "CreateGroundingFoil" 
-    //     method which already joins some glue and real groudning foil 
+    //   - the complete grounding foil, defined by the "CreateGroundingFoil"
+    //     method which already joins some glue and real groudning foil
     //     layers for the whole stave (left + right);
-    //   - 4 ladders, which are sorted according to the ALICE numbering 
+    //   - 4 ladders, which are sorted according to the ALICE numbering
     //     scheme, which depends on the layer we are building this stave for;
     //   - 2 MCMs (a left and a right one);
     //   - 2 pixel buses (a left and a right one);
     // ---
     // Arguments:
-    //   - the layer number, which determines the displacement and naming 
+    //   - the layer number, which determines the displacement and naming
     //     of sensitive volumes
-    //   - a TArrayD passed by reference which will contain the size 
+    //   - a TArrayD passed by reference which will contain the size
     //     of virtual box containing the stave
     //   - the TGeoManager
     //
@@ -3588,13 +3327,13 @@ TGeoVolumeAssembly* AliITSv11GeometrySPD::CreateStave(Int_t layer,
     TGeoVolumeAssembly *container = new TGeoVolumeAssembly(Form(
                                                  "ITSSPDlay%d-Stave",layer));
     // define the indexes of the ladders in order to have the correct order
-    // keeping in mind that the staves will be inserted as they are on layer 
-    // 2, while they are rotated around their local Y axis when inserted 
-    // on layer 1, so in this case they must be put in the "wrong" order 
-    // to turn out to be right at the end. The convention is:    
+    // keeping in mind that the staves will be inserted as they are on layer
+    // 2, while they are rotated around their local Y axis when inserted
+    // on layer 1, so in this case they must be put in the "wrong" order
+    // to turn out to be right at the end. The convention is:
     //   -|Zmax| ------> |Zmax|
     //      3   2   1   0
-    // with respect to the "native" stave reference frame, "left" is in 
+    // with respect to the "native" stave reference frame, "left" is in
     // the positive Z this leads the definition of these indexes:
     Int_t idxCentralL, idxSideL, idxCentralR, idxSideR;
 
@@ -3609,7 +3348,7 @@ TGeoVolumeAssembly* AliITSv11GeometrySPD::CreateStave(Int_t layer,
         idxCentralR = 2;
         idxSideR = 3;
     } // end if layer ==1
-    
+
      // create the two half-staves
     TArrayD sizeL, sizeR;
     TGeoVolumeAssembly *hstaveL = CreateHalfStave(kFALSE, layer, idxCentralL,
@@ -3639,9 +3378,9 @@ void AliITSv11GeometrySPD::SetAddStave(Bool_t *mask)
 {
     //
     // Define a mask which states qhich staves must be placed.
-    // It is a string which must contain '0' or '1' depending if 
+    // It is a string which must contain '0' or '1' depending if
     // a stave must be placed or not.
-    // Each place is referred to one of the staves, so the first 
+    // Each place is referred to one of the staves, so the first
     // six characters of the string will be checked.
     //
      Int_t i;
@@ -3662,7 +3401,7 @@ void AliITSv11GeometrySPD::StavesInSector(TGeoVolume *moth, TGeoManager *mgr)
     // but it inserts in the mother volume (argument 'moth') all the stuff
     // which composes the complete SPD sector.
     // ---
-    // In the following, the stave numbering order used for arrays is the 
+    // In the following, the stave numbering order used for arrays is the
     // same as defined in the GetSectorMountingPoints():
     //                         /5
     //                        /\/4
@@ -3672,9 +3411,9 @@ void AliITSv11GeometrySPD::StavesInSector(TGeoVolume *moth, TGeoManager *mgr)
     // Arguments: see description of "CarbonFiberSector" method.
     //
 
-    Double_t shift[6];  // shift from the innermost position in the 
-                        // sector placement plane (where the stave 
-                        // edge is in the point where the rounded 
+    Double_t shift[6];  // shift from the innermost position in the
+                        // sector placement plane (where the stave
+                        // edge is in the point where the rounded
                         // corner begins)
 
     shift[0] = fgkmm * -0.691;
@@ -3683,37 +3422,37 @@ void AliITSv11GeometrySPD::StavesInSector(TGeoVolume *moth, TGeoManager *mgr)
     shift[3] = fgkmm * -0.610;
     shift[4] = fgkmm * -0.610;
     shift[5] = fgkmm * -0.610;
-    
+
     // corrections after interaction with Andrea and CAD
     Double_t corrX[6] = {0.0, 0.0, 0.0, 0.0, 0.0, 0.0};
     Double_t corrY[6] = {0.0, 0.0, 0.0, 0.0, 0.0, 0.0};
-    
+
     corrX[0] =  0.0046;
     corrX[1] = -0.0041;
     corrX[2] = corrX[3] = corrX[4] = corrX[5] = -0.0016;
-    
+
     corrY[0] = -0.0007;
     corrY[1] = -0.0009;
     corrY[2] = corrY[3] = corrY[4] = corrY[5] = -0.0003;
-    
+
     corrX[0] +=  0.00026;
     corrY[0] += -0.00080;
-    
+
     corrX[1] +=  0.00018;
     corrY[1] += -0.00086;
-    
+
     corrX[2] +=  0.00020;
     corrY[2] += -0.00062;
-    
+
     corrX[3] +=  0.00017;
     corrY[3] += -0.00076;
-    
+
     corrX[4] +=  0.00016;
     corrY[4] += -0.00096;
-    
+
     corrX[5] +=  0.00018;
     corrY[5] += -0.00107;
-    
+
     // create stave volumes (different for layer 1 and 2)
     TArrayD staveSizes1(9), staveSizes2(9), clipSize(5);
     Double_t &staveHeight = staveSizes1[2], &staveThickness = staveSizes1[0];
@@ -3727,13 +3466,13 @@ void AliITSv11GeometrySPD::StavesInSector(TGeoVolume *moth, TGeoManager *mgr)
     Double_t dx, dy;      // (xL - xR) and (yL - yR)
     Double_t widthLR;     // width of the segment L-R
     Double_t angle;       // stave rotation angle in degrees
-    Double_t diffWidth;   // difference between mounting plane width and 
+    Double_t diffWidth;   // difference between mounting plane width and
                           // stave width (smaller)
     Double_t xPos, yPos;  // final translation of the stave
     Double_t parMovement; // translation in the LR plane direction
-    
+
     staveThickness += fgkGapHalfStave;
-    
+
     // loop on staves
     Int_t i, iclip = 1;
     for (i = 0; i < 6; i++) {
@@ -3752,20 +3491,20 @@ void AliITSv11GeometrySPD::StavesInSector(TGeoVolume *moth, TGeoManager *mgr)
         // by an amount equal to the width difference
         // and then the fixed shift must also be added
         parMovement = diffWidth + shift[i];
-        // due to stave thickness, another movement must be done 
+        // due to stave thickness, another movement must be done
         // in the direction normal to the mounting plane
-        // which is computed using an internal method, in a reference 
-        // frame where the LR segment has its middle point in the origin 
+        // which is computed using an internal method, in a reference
+        // frame where the LR segment has its middle point in the origin
         // and axes parallel to the master reference frame
         if (i == 0) {
-            ParallelPosition(-0.5*staveThickness, -parMovement, angle, 
+            ParallelPosition(-0.5*staveThickness, -parMovement, angle,
                                   xPos, yPos);
         } // end if i==0
         if (i == 1) {
-            ParallelPosition( 0.5*staveThickness, -parMovement, angle, 
+            ParallelPosition( 0.5*staveThickness, -parMovement, angle,
                                   xPos, yPos);
         }else {
-            ParallelPosition( 0.5*staveThickness,  parMovement, angle, 
+            ParallelPosition( 0.5*staveThickness,  parMovement, angle,
                                   xPos, yPos);
         } // end if i==1
         // then we go into the true reference frame
@@ -3773,7 +3512,7 @@ void AliITSv11GeometrySPD::StavesInSector(TGeoVolume *moth, TGeoManager *mgr)
         yPos += yM;
         xPos += corrX[i];
         yPos += corrY[i];
-        // using the parameters found here, compute the 
+        // using the parameters found here, compute the
         // translation and rotation of this stave:
         TGeoRotation *rot = new TGeoRotation(*gGeoIdentity);
         if (i == 0 || i == 1) rot->RotateX(180.0);
@@ -3793,7 +3532,7 @@ void AliITSv11GeometrySPD::StavesInSector(TGeoVolume *moth, TGeoManager *mgr)
                 rotClip->RotateX(180.0);
                 Double_t x = staveSizes2[3] + fgkGapHalfStave;
                 Double_t y = staveSizes2[4];
-                Double_t z[4] = { staveSizes2[5], staveSizes2[6], 
+                Double_t z[4] = { staveSizes2[5], staveSizes2[6],
                                   staveSizes2[7], staveSizes2[8] };
                 for (j = 0; j < 4; j++) {
                     TGeoCombiTrans *trClip = new TGeoCombiTrans(x, y, z[j],
@@ -3804,6 +3543,31 @@ void AliITSv11GeometrySPD::StavesInSector(TGeoVolume *moth, TGeoManager *mgr)
             } // end if i!=2
         } // end if i==0||i==1 else
     } // end for i
+    
+    
+    // Add a box representing the collector for cooling tubes
+    Double_t collWidth     = fgkmm * 22.0;
+    Double_t collLength    = fgkmm * 50.0;
+    Double_t collThickness = fgkmm *  7.0;
+    Double_t collInSize    = fgkmm * 10.5;
+    
+    TGeoMedium *medColl   = GetMedium("INOX$");
+    TGeoMedium *medCollIn = GetMedium("COPPER$");
+    TGeoVolume *vColl     = mgr->MakeBox("ITSSPDSectorTubeColl"  , medColl, 0.5*collWidth, 0.5*collThickness, 0.5*collLength);
+    TGeoVolume *vCollIn   = mgr->MakeBox("ITSSPDSectorTubeCollIn", medCollIn, 0.5*collInSize, 0.5*collInSize, 0.5*collInSize);
+    vColl->SetLineColor(kGreen+2);
+    vCollIn->SetLineColor(kYellow);
+    
+    TGeoTranslation *tr1 = new TGeoTranslation( 0.1, 1.2,  35.0);
+    TGeoTranslation *tr2 = new TGeoTranslation(-0.1, 1.2, -35.0);
+    TGeoTranslation *tr3 = new TGeoTranslation( 0.1, 1.2 - 0.5*(collThickness+collInSize),  35.0 + 0.5*(collLength - collInSize));
+    TGeoTranslation *tr4 = new TGeoTranslation(-0.1, 1.2 - 0.5*(collThickness+collInSize), -35.0 - 0.5*(collLength - collInSize));
+    
+    moth->AddNode(vColl, 0, tr1);
+    moth->AddNode(vColl, 1, tr2);
+    moth->AddNode(vCollIn, 0, tr3);
+    moth->AddNode(vCollIn, 1, tr4);
+    
 }
 //______________________________________________________________________
 void AliITSv11GeometrySPD::ParallelPosition(Double_t dist1, Double_t dist2,
@@ -3811,16 +3575,16 @@ void AliITSv11GeometrySPD::ParallelPosition(Double_t dist1, Double_t dist2,
 {
     //
     // Performs the following steps:
-    // 1 - finds a straight line parallel to the one passing through 
+    // 1 - finds a straight line parallel to the one passing through
     //     the origin and with angle 'phi' with X axis(phi in RADIANS);
-    // 2 - finds another line parallel to the previous one, with a 
+    // 2 - finds another line parallel to the previous one, with a
     //     distance 'dist1' from it
-    // 3 - takes a reference point in the second line in the intersection 
+    // 3 - takes a reference point in the second line in the intersection
     //     between the normal to both lines  passing through the origin
-    // 4 - finds a point whith has distance 'dist2' from this reference, 
+    // 4 - finds a point whith has distance 'dist2' from this reference,
     //     in the second line (point 2)
     // ----
-    // According to the signs given to dist1 and dist2, the point is 
+    // According to the signs given to dist1 and dist2, the point is
     // found in different position w.r. to the origin
     // compute the point
     //
@@ -3872,109 +3636,7 @@ Double_t AliITSv11GeometrySPD::GetSPDSectorTranslation(
     c = y0-a*x0-r*b;
     return -c;
 }
-//______________________________________________________________________
-void AliITSv11GeometrySPD::CreateFigure0(const Char_t *filepath,
-                                         const Char_t *type,
-                                         TGeoManager *mgr) const
-{
-    //
-    // Creates Figure 0 for the documentation of this class. In this
-    // specific case, it creates the X,Y cross section of the SPD suport
-    // section, center and ends. The output is written to a standard
-    // file name to the path specificed.
-    // Inputs:
-    //   const Char_t *filepath  Path where the figure is to be drawn
-    //   const Char_t *type      The type of file, default is gif.
-    //   TGeoManager  *mgr       The TGeoManager default gGeoManager
-    // Output:
-    //   none.
-    // Return:
-    //   none.
-    //
-    TGeoXtru *sA0,*sA1,*sB0,*sB1;
-    //TPolyMarker *pmA,*pmB;
-    TPolyLine plA0,plA1,plB0,plB1;
-    TCanvas *canvas;
-    TLatex txt;
-    Double_t x=0.0,y=0.0;
-    Int_t i,kNRadii=6;
-
-    if(strcmp(filepath,"")){
-        Error("CreateFigure0","filepath=%s type=%s",filepath,type);
-    } // end if
-    //
-    sA0 = (TGeoXtru*) mgr->GetVolume("ITSSPDCarbonFiberSupportSectorA0_1")->
-              GetShape();
-    sA1 = (TGeoXtru*) mgr->GetVolume("ITSSPDCarbonFiberSupportSectorAirA1_1")->
-              GetShape();
-    sB0 = (TGeoXtru*) mgr->GetVolume("ITSSPDCarbonFiberSupportSectorEndB0_1")->
-             GetShape();
-    sB1 = (TGeoXtru*) mgr->GetVolume("ITSSPDCarbonFiberSupportSectorEndAirB1_1"
-           )->GetShape();
-    //pmA = new TPolyMarker();
-    //pmA.SetMarkerStyle(2); // +
-    //pmA.SetMarkerColor(7); // light blue
-    //pmB = new TPolyMarker();
-    //pmB.SetMarkerStyle(5); // X
-    //pmB.SetMarkerColor(6); // purple
-    plA0.SetPolyLine(sA0->GetNvert());
-    plA0.SetLineColor(1); // black
-    plA0.SetLineStyle(1);
-    plA1.SetPolyLine(sA1->GetNvert());
-    plA1.SetLineColor(2); // red
-    plA1.SetLineStyle(1);
-    plB0.SetPolyLine(sB0->GetNvert());
-    plB0.SetLineColor(3); // Green
-    plB0.SetLineStyle(2);
-    plB1.SetPolyLine(sB1->GetNvert());
-    plB1.SetLineColor(4); // Blue
-    plB1.SetLineStyle(2);
-    //for(i=0;i<kNRadii;i++) pmA.SetPoint(i,xyB1p[i][0],xyB1p[i][1]);
-    //for(i=0;i<kNRadii;i++) pmB.SetPoint(i,xyB1p[i][0],xyB1p[i][1]);
-    for(i=0;i<sA0->GetNvert();i++) plA0.SetPoint(i,sA0->GetX(i),sA0->GetY(i));
-    for(i=0;i<sA1->GetNvert();i++) plA1.SetPoint(i,sA1->GetX(i),sA1->GetY(i));
-    for(i=0;i<sB0->GetNvert();i++) plB0.SetPoint(i,sB0->GetX(i),sB0->GetY(i));
-    for(i=0;i<sB1->GetNvert();i++) plB1.SetPoint(i,sB1->GetX(i),sB1->GetY(i));
-    canvas = new TCanvas("AliITSv11GeometrySPDFig0","",1000,1000);
-    canvas->Range(-3.,-3.,3.,3.);
-    txt.SetTextSize(0.05);
-    txt.SetTextAlign(33);
-    txt.SetTextColor(1);
-    txt.DrawLatex(2.9,2.9,"Section A-A outer Carbon Fiber surface");
-    txt.SetTextColor(2);
-    txt.DrawLatex(2.9,2.5,"Section A-A Inner Carbon Fiber surface");
-    txt.SetTextColor(3);
-    txt.DrawLatex(2.9,2.1,"Section E-E outer Carbon Fiber surface");
-    txt.SetTextColor(4);
-    txt.DrawLatex(2.9,1.7,"Section E-E Inner Carbon Fiber surface");
-    plA0.Draw();
-    plA1.Draw();
-    plB0.Draw();
-    plB1.Draw();
-    //pmA.Draw();
-    //pmB.Draw();
-    //
-    x = 1.0;
-    y = -2.5;
-    Char_t chr[3];
-    for(i=0;i<kNRadii;i++){
-        sprintf(chr,"%2d",i);txt.DrawLatex(x-0.1,y,chr);
-        sprintf(chr,"%8.4f",5.000);txt.DrawLatex(x,y,chr);
-        sprintf(chr,"%8.4f",5.000);txt.DrawLatex(x+0.5,y,chr);
-        sprintf(chr,"%8.4f",5.000);txt.DrawLatex(x+1.0,y,chr);
-        sprintf(chr,"%8.4f",5.000);txt.DrawLatex(x+1.5,y,chr);
-        sprintf(chr,"%8.4f",5.000);txt.DrawLatex(x+2.0,y,chr);
-        if(kTRUE) txt.DrawLatex(x+2.5,y,"A-A/E-E");
-        else txt.DrawLatex(x+2.5,y,"E-E");
-    } // end for i
-    txt.DrawLatex(x,y,"x_{c} mm");
-    txt.DrawLatex(x+0.5,y,"y_{c} mm");
-    txt.DrawLatex(x+1.0,y,"R mm");
-    txt.DrawLatex(x+1.5,y,"#theta_{start}^{#circle}");
-    txt.DrawLatex(x+2.0,y,"#theta_{end}^{#circle}");
-    txt.DrawLatex(x+2.5,y,"Section");
-    //
-}
+
 //______________________________________________________________________
 void AliITSv11GeometrySPD::PrintAscii(ostream *os) const
 {
@@ -4009,9 +3671,9 @@ void AliITSv11GeometrySPD::PrintAscii(ostream *os) const
     for(i=0;i<fSPDsectorX1.GetSize();i++) *os<< fSPDsectorX1.GetAt(i) << " ";
     for(i=0;i<fSPDsectorX1.GetSize();i++) *os<< fSPDsectorY1.GetAt(i) << " ";
     *os<<10<<" "<< 2 <<" " << 6 << " "<< 3 <<" ";
-    for(k=0;k<10;k++)for(i=0;i<6;i++)for(j=0;j<3;j++) 
+    for(k=0;k<10;k++)for(i=0;i<6;i++)for(j=0;j<3;j++)
         *os<<fTubeEndSector[k][0][i][j]<<" ";
-    for(k=0;k<10;k++)for(i=0;i<6;i++)for(j=0;j<3;j++) 
+    for(k=0;k<10;k++)for(i=0;i<6;i++)for(j=0;j<3;j++)
         *os<<fTubeEndSector[k][1][i][j]<<" ";
     os->flags(fmt); // reset back to old Formating.
     return;
@@ -4031,14 +3693,18 @@ void AliITSv11GeometrySPD::ReadAscii(istream* is)
     //
     Int_t i,j,k,n;
     Double_t gapLadder,GapHalfStave;
-
+    const Int_t kLimits = 100;
     *is>>gapLadder>>GapHalfStave>>n;
     if(n!=6){
-        Warning("ReadAscii","fAddStave Array !=6 n=%d",n);
+      AliError(Form("fAddStave Array !=6 n=%d",n));
         return;
     } // end if
     for(i=0;i<n;i++) *is>>fAddStave[i];
     *is>>n;
+    if(n<0 || n> kLimits){
+      AliError("Anomalous value for parameter n");
+      return;
+    } 
     fSPDsectorX0.Set(n);
     fSPDsectorY0.Set(n);
     fSPDsectorX1.Set(n);
@@ -4053,9 +3719,9 @@ void AliITSv11GeometrySPD::ReadAscii(istream* is)
                 "found [%d][%d][%d]",i,j,n);
         return;
     } // end if
-    for(k=0;k<10;k++)for(i=0;i<6;i++)for(j=0;j<3;j++) 
+    for(k=0;k<10;k++)for(i=0;i<6;i++)for(j=0;j<3;j++)
         *is>>fTubeEndSector[k][0][i][j];
-    for(k=0;k<10;k++)for(i=0;i<6;i++)for(j=0;j<3;j++) 
+    for(k=0;k<10;k++)for(i=0;i<6;i++)for(j=0;j<3;j++)
         *is>>fTubeEndSector[k][1][i][j];
     return;
 }
@@ -4093,90 +3759,4 @@ istream &operator>>(istream &is,AliITSv11GeometrySPD &s)
     s.ReadAscii(&is);
     return is;
 }
-//
-//______________________________________________________________________
-Bool_t AliITSv11GeometrySPD::Make2DCrossSections(TPolyLine &a0,TPolyLine &a1,
-                             TPolyLine &b0,TPolyLine &b1,TPolyMarker &p)const
-{
-    //
-    // Fill the objects with the points representing
-    // a0 the outer carbon fiber SPD sector shape Cross Section A
-    // a1 the inner carbon fiber SPD sector shape Cross Section A
-    // b0 the outer carbon fiber SPD sector shape Cross Section B
-    // b1 the inner carbon fiber SPD sector shape Cross Section B
-    //
-    // Inputs:
-    //   TPolyLine &a0   The outer carbon fiber SPD sector shape
-    //   TPolyLine &a1   The Inner carbon fiber SPD sector shape
-    //   TPolyLine &b0   The outer carbon fiber SPD sector shape
-    //   TPolyLine &b1   The Inner carbon fiber SPD sector shape
-    //   TPolyMarker &p  The points where the ladders are to be placed
-    // Outputs:
-    //   TPolyLine &a0   The shape filled with the points
-    //   TPolyLine &a1   The shape filled with the points
-    //   TPolyLine &b0   The shape filled with the points
-    //   TPolyLine &b1   The shape filled with the points
-    //   TPolyMarker &p  The filled array of points
-    // Return:
-    //     An error flag.
-    //
-    Int_t n0,n1,i;
-    Double_t x,y;
-    TGeoVolume *a0V,*a1V,*b0V,*b1V;
-    TGeoXtru *a0S,*a1S,*b0S,*b1S;
-    TGeoManager *mgr = gGeoManager;
-
-    a0V = mgr->GetVolume("ITS SPD Carbon fiber support Sector A0");
-    a0S = dynamic_cast<TGeoXtru*>(a0V->GetShape());
-    n0 = a0S->GetNvert();
-    a0.SetPolyLine(n0+1);
-    //for(i=0;i<fSPDsectorPoints0.GetSize();i++) 
-    //  printf("%d %d %d\n",i,fSPDsectorPoints0[i],fSPDsectorPoints1[i]);
-    for(i=0;i<n0;i++){
-        x = a0S->GetX(i);
-          y = a0S->GetY(i);
-          //printf("%d %g %g\n",i,x,y);
-        a0.SetPoint(i,x,y);
-          if(i==0) a0.SetPoint(n0,x,y);
-    } // end for i
-    a1V = mgr->GetVolume("ITSSPDCarbonFiberSupportSectorAirA1");
-    a1S = dynamic_cast<TGeoXtru*>(a1V->GetShape());
-    n1 = a1S->GetNvert();
-    a1.SetPolyLine(n1+1);
-    for(i=0;i<n1;i++){
-        x = a1S->GetX(i);
-          y = a1S->GetY(i);
-        a1.SetPoint(i,x,y);
-          if(i==0) a1.SetPoint(n1,x,y);
-    } // end for i
-    // Cross Section B
-    b0V = mgr->GetVolume("ITSSPDCarbonFiberSupportSectorEndB0");
-    b0S = dynamic_cast<TGeoXtru*>(b0V->GetShape());
-    n0 = b0S->GetNvert();
-    b0.SetPolyLine(n0+1);
-    for(i=0;i<n0;i++){
-        x = b0S->GetX(i);
-          y = b0S->GetY(i);
-        b0.SetPoint(i,x,y);
-          if(i==0) b0.SetPoint(n0,x,y);
-    } // end for i
-    b1V = mgr->GetVolume("ITSSPDCarbonFiberSupportSectorEndAirB1");
-    b1S = dynamic_cast<TGeoXtru*>(b1V->GetShape());
-    n1 = b1S->GetNvert();
-    b1.SetPolyLine(n1+1);
-    for(i=0;i<n1;i++){
-        x = b1S->GetX(i);
-          y = b1S->GetY(i);
-        b1.SetPoint(i,x,y);
-          if(i==0) b1.SetPoint(n1,x,y);
-    } // end for i
-    //
-    Double_t x0,y0,x1,y1;
-    p.SetPolyMarker(2*fSPDsectorX0.GetSize());
-    for(i=0;i<fSPDsectorX0.GetSize();i++){
-          GetSectorMountingPoints(i,x0,y0,x1,y1);
-          p.SetPoint(2*i,x0,y0);
-          p.SetPoint(2*i+1,x1,y1);
-    } // end for i
-    return kTRUE;
-}
+