]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - ITS/AliITSvPPRasymmFMD.cxx
Completely reengineered version of CMake build system (Johny)
[u/mrichter/AliRoot.git] / ITS / AliITSvPPRasymmFMD.cxx
index edad6079665f91bb969985d5c4c42e1aa16774ad..ba3267e8a4ec664cfa95c75d5c99bf997ad75a44 100644 (file)
 
 // See AliITSvPPRasymmFMD::StepManager().
 
-#include <TClonesArray.h>
-#include <TGeometry.h>
-#include <TLorentzVector.h>
-#include <TGeoMatrix.h>
-#include <TGeoPhysicalNode.h>
 #include <TArrayD.h>
 #include <TArrayF.h>
-#include <TString.h>
-#include <TNode.h>
-#include <TTUBE.h>
+#include <TClonesArray.h>
+#include <TGeoGlobalMagField.h>
 #include <TGeoManager.h>
+#include <TGeoMatrix.h>
+#include <TGeoPhysicalNode.h>
 #include <TGeoVolume.h>
+#include <TLorentzVector.h>
+#include <TString.h>
 #include <TVirtualMC.h>
 
 #include "AliITS.h"
@@ -63,6 +61,7 @@
 #include "AliMagF.h"
 #include "AliRun.h"
 #include "AliTrackReference.h"
+#include "AliGeomManager.h"
 
 
 #define GEANTGEOMETRY kTRUE
@@ -72,15 +71,10 @@ ClassImp(AliITSvPPRasymmFMD)
 //______________________________________________________________________
 AliITSvPPRasymmFMD::AliITSvPPRasymmFMD():
 AliITS(),                    // Default AliITS Constructor
-fGeomDetOut(kFALSE),         // Flag to write .det file out
-fGeomDetIn(kFALSE),          // Flag to read .det file or directly from Geat.
 fByThick(kTRUE),             // Flag to use services materials by thickness
                              // ture, or mass false.
 fMajorVersion(IsVersion()),  // Major version number == IsVersion
 fMinorVersion(2),           // Minor version number
-fEuclidGeomDet("$ALICE_ROOT/ITS/ITSgeometry_vPPRasymm2.det"),// file where detector transormation are define.
-fRead("$ALICE_ROOT/ITS/ITSgeometry_vPPRasymm.det"),//! file name to read .det file
-fWrite("$ALICE_ROOT/ITS/ITSgeometry_vPPRasymm2.det"),//! file name to write .det file
 fDet1(200.0),     // thickness of detector in SPD layer 1 [microns]
 fDet2(200.0),     // thickness of detector in SPD layer 2 [microns]
 fChip1(150.0),    // thickness of chip in SPD layer 1 [microns]
@@ -100,18 +94,14 @@ fIgm((AliITSVersion_t)fMajorVersion,fMinorVersion)//! Get access to decoding
     //   none.
     fIgm.SetGeometryName("");
 }
+
 //______________________________________________________________________
 AliITSvPPRasymmFMD::AliITSvPPRasymmFMD(const Char_t *title):
 AliITS(title),               // Standard AliITS Constructor
-fGeomDetOut(kFALSE),         // Flag to write .det file out
-fGeomDetIn(kFALSE),          // Flag to read .det file or directly from Geat.
 fByThick(kTRUE),             // Flag to use services materials by thickness
                              // ture, or mass false.
 fMajorVersion(IsVersion()),  // Major version number == IsVersion
 fMinorVersion(2),            // Minor version number
-fEuclidGeomDet("$ALICE_ROOT/ITS/ITSgeometry_vPPRasymm2.det"),// file where detector transormation are define.
-fRead("$ALICE_ROOT/ITS/ITSgeometry_vPPRasymm.det"),//! file name to read .det file
-fWrite("$ALICE_ROOT/ITS/ITSgeometry_vPPRasymm2.det"),//! file name to write .det file
 fDet1(200.0),     // thickness of detector in SPD layer 1 [microns]
 fDet2(200.0),     // thickness of detector in SPD layer 2 [microns]
 fChip1(150.0),    // thickness of chip in SPD layer 1 [microns]
@@ -149,18 +139,14 @@ fIgm((AliITSVersion_t)fMajorVersion,fMinorVersion)//! Get access to decoding
     fIdSens    = new Int_t[fIdN];
     for(i=0;i<fIdN;i++) fIdSens[i] = 0;
 }
+
 //______________________________________________________________________
 AliITSvPPRasymmFMD::AliITSvPPRasymmFMD(const char *name, const char *title):
 AliITS(name,title),          // Extended AliITS Constructor
-fGeomDetOut(kFALSE),         // Flag to write .det file out
-fGeomDetIn(kFALSE),          // Flag to read .det file or directly from Geat.
 fByThick(kTRUE),             // Flag to use services materials by thickness
                              // ture, or mass false.
 fMajorVersion(IsVersion()),  // Major version number == IsVersion
 fMinorVersion(2),            // Minor version number
-fEuclidGeomDet("$ALICE_ROOT/ITS/ITSgeometry_vPPRasymm2.det"),// file where detector transormation are define.
-fRead("$ALICE_ROOT/ITS/ITSgeometry_vPPRasymm.det"),//! file name to read .det file
-fWrite("$ALICE_ROOT/ITS/ITSgeometry_vPPRasymm2.det"),//! file name to write .det file
 fDet1(200.0),     // thickness of detector in SPD layer 1 [microns]
 fDet2(200.0),     // thickness of detector in SPD layer 2 [microns]
 fChip1(150.0),    // thickness of chip in SPD layer 1 [microns]
@@ -206,21 +192,17 @@ fIgm((AliITSVersion_t)fMajorVersion,fMinorVersion)//! Get access to decoding
     SetThicknessChip2();
     SetDensityServicesByThickness();
 }
+
 //______________________________________________________________________
 AliITSvPPRasymmFMD::AliITSvPPRasymmFMD(const AliITSvPPRasymmFMD &source):
 AliITS(source.GetName(),source.GetTitle()),// Extended AliITS Constructor
-fGeomDetOut(source.GetWriteDet()),// Flag to write .det file out
-fGeomDetIn(source.GetReadDet()),  // Flag to read .det file or directly 
                                   // from Geat.
 fByThick(source.IsDensityServicesByThickness()),// Flag to use services 
                                                 // materials by thickness
                                                 // ture, or mass false.
 fMajorVersion(source.IsVersion()),  // Major version number == IsVersion
 fMinorVersion(source.GetMinorVersion()),// Minor version number
-fEuclidGeomDet(source.GetEuclidFile()),// file where detector transormation 
                                        // are define.
-fRead(source.GetReadDetFileName()),  //! file name to read .det file
-fWrite(source.GetWriteDetFileName()),//! file name to write .det file
 fDet1(source.GetThicknessDet1()),//thickness of detector in SPD layer 1[micron]
 fDet2(source.GetThicknessDet2()),//thickness of detector in SPD layer 2[micron]
 fChip1(source.GetThicknessChip1()),//thickness of chip in SPD layer 1 [microns]
@@ -245,6 +227,7 @@ fIgm(source.GetGeomInit())//! Get access to decoding
     Warning("Copy Constructor","Not allowed to copy AliITSvPPRasymmFMD");
     return;
 }
+
 //______________________________________________________________________
 AliITSvPPRasymmFMD& AliITSvPPRasymmFMD::operator=(const AliITSvPPRasymmFMD 
                                                  &source){
@@ -262,6 +245,7 @@ AliITSvPPRasymmFMD& AliITSvPPRasymmFMD::operator=(const AliITSvPPRasymmFMD
     Warning("= operator","Not allowed to copy AliITSvPPRasymmFMD");
     return *this; // return null pointer, copy not allowed.
 }
+
 //______________________________________________________________________
 AliITSvPPRasymmFMD::~AliITSvPPRasymmFMD() {
     //    Standard destructor for the ITS version 10.
@@ -272,6 +256,7 @@ AliITSvPPRasymmFMD::~AliITSvPPRasymmFMD() {
     // Return:
     //   none.
 }
+
 //______________________________________________________________________
 void AliITSvPPRasymmFMD::AddAlignableVolumes() const
 {
@@ -292,6 +277,9 @@ void AliITSvPPRasymmFMD::AddAlignableVolumes() const
     return;
   }
 
+  AliGeomManager::ELayerID layerId;
+  Int_t modUID, modnum;
+
   if( !gGeoManager->SetAlignableEntry("ITS","ALIC_1/ITSV_1") )
     AliFatal("Unable to set alignable entry!!");    
 
@@ -312,6 +300,8 @@ void AliITSvPPRasymmFMD::AddAlignableVolumes() const
 
   //===== SPD layer1 =====
   {
+    layerId = AliGeomManager::kSPD1;
+    modnum = 0;
     TString str0 = "ALIC_1/ITSV_1/ITSD_1/IT12_1/I12B_";
     TString str1 = "/I10B_";
     TString str1Bis = "/L1H-STAVE";
@@ -323,8 +313,8 @@ void AliITSvPPRasymmFMD::AddAlignableVolumes() const
     TString halfStave;
     TString module;
     
-    for(Int_t cSect = 0; cSect<10; cSect++) {
-      
+    for(Int_t cSect = 0; cSect<10; cSect++)
+    {
       sector = str0;
       sector += cSect+1; // this is one full sector
       strEntryName1 = strSPD;
@@ -335,8 +325,8 @@ void AliITSvPPRasymmFMD::AddAlignableVolumes() const
       if(!gGeoManager->SetAlignableEntry(strEntryName1.Data(),sector.Data()))
        AliFatal("Unable to set alignable entry!!");    
       
-      for(Int_t cStave = 0; cStave<2; cStave++) {
-       
+      for(Int_t cStave = 0; cStave<2; cStave++)
+      {
        stave = sector;
        stave += str1;
        stave += cStave+1;
@@ -345,8 +335,8 @@ void AliITSvPPRasymmFMD::AddAlignableVolumes() const
        strEntryName2 += cStave;
        //printf("%s   ==   %s\n",strEntryName2.Data(),stave.Data()); // this is a stave
 
-       for(Int_t cHS=0; cHS<2; cHS++) {
-
+       for(Int_t cHS=0; cHS<2; cHS++)
+       {
          halfStave = stave;
          halfStave += str1Bis;
          halfStave += cHS;
@@ -358,8 +348,9 @@ void AliITSvPPRasymmFMD::AddAlignableVolumes() const
          if(!gGeoManager->SetAlignableEntry(strEntryName3.Data(),halfStave.Data()))
            AliFatal("Unable to set alignable entry!!");    
 
-         for(Int_t cLadder = 0; cLadder<2; cLadder++) {
-           
+         for(Int_t cLadder = 0; cLadder<2; cLadder++)
+         {
+           modUID = AliGeomManager::LayerToVolUID(layerId,modnum++);
            module = halfStave;
            module += str2;
            module += cLadder+cHS*2+1;
@@ -367,15 +358,12 @@ void AliITSvPPRasymmFMD::AddAlignableVolumes() const
            strEntryName4 += strLadder;
            strEntryName4 += cLadder+cHS*2;
            //printf("%s   ==   %s\n",strEntryName4.Data(),module.Data());
-           if(!gGeoManager->SetAlignableEntry(strEntryName4.Data(),module.Data()))
+           if(!gGeoManager->SetAlignableEntry(strEntryName4.Data(),module.Data(),modUID))
              AliFatal("Unable to set alignable entry!!");    
 
            // Creates the TGeo Local to Tracking transformation matrix ...
-           TGeoPNEntry *alignableEntry = gGeoManager->GetAlignableEntry(strEntryName4.Data());
-           const char *path = alignableEntry->GetTitle();
-           if (!gGeoManager->cd(path))
-             AliFatal(Form("Volume path %s not valid!",path));
-           TGeoHMatrix* globMatrix = gGeoManager->GetCurrentMatrix();
+           TGeoPNEntry *alignableEntry = gGeoManager->GetAlignableEntryByUID(modUID);
+           TGeoHMatrix *globMatrix = alignableEntry->GetGlobalOrig();
            gtrans = globMatrix->GetTranslation();
            memcpy(&rotMatrix[0], globMatrix->GetRotationMatrix(), 9*sizeof(Double_t));
            al = TMath::ATan2(rotMatrix[1],rotMatrix[0]);
@@ -385,7 +373,7 @@ void AliITSvPPRasymmFMD::AddAlignableVolumes() const
 
            matLtoT->SetDy(-fChip1*0.0001/2.);
            matLtoT->SetDz(-gtrans[2]);
-           rotMatrix[0]= 0;  rotMatrix[1]= 1;  rotMatrix[2]= 0; // + rotation
+           rotMatrix[0]= 0;  rotMatrix[1]=-1;  rotMatrix[2]= 0; // + rotation
            rotMatrix[3]=-1;  rotMatrix[4]= 0;  rotMatrix[5]= 0; // ! flip in y for the SPD1 only
            rotMatrix[6]= 0;  rotMatrix[7]= 0;  rotMatrix[8]=-1;
            TGeoRotation rot;
@@ -400,12 +388,14 @@ void AliITSvPPRasymmFMD::AddAlignableVolumes() const
     }
     
     //===== SPD layer2 =====
+    layerId = AliGeomManager::kSPD2;
+    modnum = 0;
     str1Bis = "/L2H-STAVE";
     str1 = "/I20B_";
     str2 = "/I1D7_";
     
-    for(Int_t cSect = 0; cSect<10; cSect++) {
-      
+    for(Int_t cSect = 0; cSect<10; cSect++)
+    {
       sector = str0;
       sector += cSect+1; // this is one full sector
       strEntryName1 = strSPD;
@@ -417,8 +407,8 @@ void AliITSvPPRasymmFMD::AddAlignableVolumes() const
       // we don't need the previous lines because the whole sector is already define
       // with first layer ...
       
-      for(Int_t cStave =0; cStave<4; cStave++) {
-       
+      for(Int_t cStave =0; cStave<4; cStave++)
+      {
        stave = sector;
        stave += str1;
        stave += cStave+1;
@@ -426,8 +416,8 @@ void AliITSvPPRasymmFMD::AddAlignableVolumes() const
        strEntryName2 += strStave;
        strEntryName2 += cStave;
        
-       for(Int_t cHS=0; cHS<2; cHS++) {
-         
+       for(Int_t cHS=0; cHS<2; cHS++)
+       {
          halfStave = stave;
          halfStave += str1Bis;
          halfStave += cHS;
@@ -439,8 +429,9 @@ void AliITSvPPRasymmFMD::AddAlignableVolumes() const
          if(!gGeoManager->SetAlignableEntry(strEntryName3.Data(),halfStave.Data()))
            AliFatal("Unable to set alignable entry!!");    
          
-         for(Int_t cLad =0; cLad<2; cLad++) {
-           
+         for(Int_t cLad =0; cLad<2; cLad++)
+         {
+           modUID = AliGeomManager::LayerToVolUID(layerId,modnum++);
            module = halfStave;
            module += str2;
            module += cLad+cHS*2+1;
@@ -448,15 +439,12 @@ void AliITSvPPRasymmFMD::AddAlignableVolumes() const
            strEntryName4 += strLadder;
            strEntryName4 += cLad+cHS*2;
            //printf("%s   ==   %s\n",strEntryName4.Data(),module.Data());
-           if(!gGeoManager->SetAlignableEntry(strEntryName4.Data(),module.Data()))
+           if(!gGeoManager->SetAlignableEntry(strEntryName4.Data(),module.Data(),modUID))
              AliFatal("Unable to set alignable entry!!");
            
            // Creates the TGeo Local to Tracking transformation matrix ...
-           TGeoPNEntry *alignableEntry = gGeoManager->GetAlignableEntry(strEntryName4.Data());
-           const char *path = alignableEntry->GetTitle();
-           if (!gGeoManager->cd(path))
-             AliFatal(Form("Volume path %s not valid!",path));
-           TGeoHMatrix* globMatrix = gGeoManager->GetCurrentMatrix();
+           TGeoPNEntry *alignableEntry = gGeoManager->GetAlignableEntryByUID(modUID);
+           TGeoHMatrix *globMatrix = alignableEntry->GetGlobalOrig();
            gtrans = globMatrix->GetTranslation();
            memcpy(&rotMatrix[0], globMatrix->GetRotationMatrix(), 9*sizeof(Double_t));
            al = TMath::ATan2(rotMatrix[1],rotMatrix[0]) + TMath::Pi();
@@ -483,14 +471,16 @@ void AliITSvPPRasymmFMD::AddAlignableVolumes() const
   
   //===== SDD layer1 =====
   {
+    layerId = AliGeomManager::kSDD1;
+    modnum = 0;
     TString str0 = "ALIC_1/ITSV_1/ITSD_1/IT34_1/I004_";
     TString str1 = "/I302_";
 
     TString ladder;
     TString wafer;
 
-    for(Int_t c1 = 1; c1<=14; c1++){
-
+    for(Int_t c1 = 1; c1<=14; c1++)
+    {
       ladder = str0;
       ladder += c1; // the set of wafers from one ladder
       strEntryName1 = strSDD;
@@ -501,24 +491,22 @@ void AliITSvPPRasymmFMD::AddAlignableVolumes() const
        AliFatal("Unable to set alignable entry!!");    
       //printf("%s    ==    %s\n",strEntryName1.Data(),ladder.Data());
 
-      for(Int_t c2 =1; c2<=6; c2++){
-
+      for(Int_t c2 =1; c2<=6; c2++)
+      {
+       modUID = AliGeomManager::LayerToVolUID(layerId,modnum++);
        wafer = ladder;
        wafer += str1;
        wafer += c2;    // one wafer
        strEntryName2 = strEntryName1;
        strEntryName2 += strSensor;
        strEntryName2 += (c2-1);
-       if(!gGeoManager->SetAlignableEntry(strEntryName2.Data(),wafer.Data()))
+       if(!gGeoManager->SetAlignableEntry(strEntryName2.Data(),wafer.Data(),modUID))
          AliFatal("Unable to set alignable entry!!");    
        //printf("%s    ==    %s\n",strEntryName2.Data(),wafer.Data());
 
        // Creates the TGeo Local to Tracking transformation matrix ...
-       TGeoPNEntry *alignableEntry = gGeoManager->GetAlignableEntry(strEntryName2.Data());
-       const char *path = alignableEntry->GetTitle();
-       if (!gGeoManager->cd(path))
-         AliFatal(Form("Volume path %s not valid!",path));
-       TGeoHMatrix* globMatrix = gGeoManager->GetCurrentMatrix();
+       TGeoPNEntry *alignableEntry = gGeoManager->GetAlignableEntryByUID(modUID);
+       TGeoHMatrix *globMatrix = alignableEntry->GetGlobalOrig();
        gtrans = globMatrix->GetTranslation();
        memcpy(&rotMatrix[0], globMatrix->GetRotationMatrix(), 9*sizeof(Double_t));
        al = TMath::ATan2(rotMatrix[1],rotMatrix[0]) + TMath::Pi();
@@ -543,14 +531,16 @@ void AliITSvPPRasymmFMD::AddAlignableVolumes() const
 
   //===== SDD layer2 =====
   {
+    layerId = AliGeomManager::kSDD2;
+    modnum = 0;
     TString str0 = "ALIC_1/ITSV_1/ITSD_1/IT34_1/I005_";
     TString str1 = "/I402_";
 
     TString ladder;
     TString wafer;
 
-    for(Int_t c1 = 1; c1<=22; c1++){
-
+    for(Int_t c1 = 1; c1<=22; c1++)
+    {
       ladder = str0;
       ladder += c1; // the set of wafers from one ladder
       strEntryName1 = strSDD;
@@ -561,24 +551,22 @@ void AliITSvPPRasymmFMD::AddAlignableVolumes() const
        AliFatal("Unable to set alignable entry!!");    
       //printf("%s    ==    %s\n",strEntryName1.Data(),ladder.Data());
 
-      for(Int_t c2 =1; c2<=8; c2++){
-
+      for(Int_t c2 =1; c2<=8; c2++)
+      {
+       modUID = AliGeomManager::LayerToVolUID(layerId,modnum++);
        wafer = ladder;
        wafer += str1;
        wafer += c2;    // one wafer
        strEntryName2 = strEntryName1;
        strEntryName2 += strSensor;
        strEntryName2 += (c2-1);
-       if(!gGeoManager->SetAlignableEntry(strEntryName2.Data(),wafer.Data()))
+       if(!gGeoManager->SetAlignableEntry(strEntryName2.Data(),wafer.Data(),modUID))
          AliFatal("Unable to set alignable entry!!");    
        //printf("%s    ==    %s\n",strEntryName2.Data(),wafer.Data());
 
        // Creates the TGeo Local to Tracking transformation matrix ...
-       TGeoPNEntry *alignableEntry = gGeoManager->GetAlignableEntry(strEntryName2.Data());
-       const char *path = alignableEntry->GetTitle();
-       if (!gGeoManager->cd(path))
-         AliFatal(Form("Volume path %s not valid!",path));
-       TGeoHMatrix* globMatrix = gGeoManager->GetCurrentMatrix();
+       TGeoPNEntry *alignableEntry = gGeoManager->GetAlignableEntryByUID(modUID);
+       TGeoHMatrix *globMatrix = alignableEntry->GetGlobalOrig();
        gtrans = globMatrix->GetTranslation();
        memcpy(&rotMatrix[0], globMatrix->GetRotationMatrix(), 9*sizeof(Double_t));
        al = TMath::ATan2(rotMatrix[1],rotMatrix[0]) + TMath::Pi();
@@ -603,14 +591,16 @@ void AliITSvPPRasymmFMD::AddAlignableVolumes() const
 
   //===== SSD layer1 =====
   {
+    layerId = AliGeomManager::kSSD1;
+    modnum = 0;
     TString str0 = "ALIC_1/ITSV_1/ITSD_1/IT56_1/I565_";
     TString str1 = "/I562_";
 
     TString ladder;
     TString wafer;
 
-    for(Int_t c1 = 1; c1<=34; c1++){
-
+    for(Int_t c1 = 1; c1<=34; c1++)
+    {
       ladder = str0;
       ladder += c1; // the set of wafers from one ladder
       strEntryName1 = strSSD;
@@ -621,24 +611,22 @@ void AliITSvPPRasymmFMD::AddAlignableVolumes() const
        AliFatal("Unable to set alignable entry!!");    
       //printf("%s    ==    %s\n",strEntryName1.Data(),ladder.Data());
 
-      for(Int_t c2 = 1; c2<=22; c2++){
-
+      for(Int_t c2 = 1; c2<=22; c2++)
+      {
+       modUID = AliGeomManager::LayerToVolUID(layerId,modnum++);
        wafer = ladder;
        wafer += str1;
        wafer += c2;     // one wafer
        strEntryName2 = strEntryName1;
        strEntryName2 += strSensor;
        strEntryName2 += (c2-1);
-       if(!gGeoManager->SetAlignableEntry(strEntryName2.Data(),wafer.Data()))
+       if(!gGeoManager->SetAlignableEntry(strEntryName2.Data(),wafer.Data(),modUID))
          AliFatal("Unable to set alignable entry!!");    
        //printf("%s    ==    %s\n",strEntryName2.Data(),wafer.Data());
 
        // Creates the TGeo Local to Tracking transformation matrix ...
-       TGeoPNEntry *alignableEntry = gGeoManager->GetAlignableEntry(strEntryName2.Data());
-       const char *path = alignableEntry->GetTitle();
-       if (!gGeoManager->cd(path))
-         AliFatal(Form("Volume path %s not valid!",path));
-       TGeoHMatrix* globMatrix = gGeoManager->GetCurrentMatrix();
+       TGeoPNEntry *alignableEntry = gGeoManager->GetAlignableEntryByUID(modUID);
+       TGeoHMatrix *globMatrix = alignableEntry->GetGlobalOrig();
        gtrans = globMatrix->GetTranslation();
        memcpy(&rotMatrix[0], globMatrix->GetRotationMatrix(), 9*sizeof(Double_t));
        al = TMath::ATan2(rotMatrix[1],rotMatrix[0]) + TMath::Pi();
@@ -663,14 +651,16 @@ void AliITSvPPRasymmFMD::AddAlignableVolumes() const
 
   //===== SSD layer2 =====
   {
+    layerId = AliGeomManager::kSSD2;
+    modnum = 0;
     TString str0 = "ALIC_1/ITSV_1/ITSD_1/IT56_1/I569_";
     TString str1 = "/I566_";
 
     TString ladder;
     TString wafer;
 
-    for(Int_t c1 = 1; c1<=38; c1++){
-
+    for(Int_t c1 = 1; c1<=38; c1++)
+    {
       ladder = str0;
       ladder += c1; // the set of wafers from one ladder
       strEntryName1 = strSSD;
@@ -681,24 +671,22 @@ void AliITSvPPRasymmFMD::AddAlignableVolumes() const
        AliFatal("Unable to set alignable entry!!");    
       //printf("%s    ==    %s\n",strEntryName1.Data(),ladder.Data());
 
-      for(Int_t c2 = 1; c2<=25; c2++){
-
+      for(Int_t c2 = 1; c2<=25; c2++)
+      {
+       modUID = AliGeomManager::LayerToVolUID(layerId,modnum++);
        wafer = ladder;
        wafer += str1;
        wafer += c2;     // one wafer
        strEntryName2 = strEntryName1;
        strEntryName2 += strSensor;
        strEntryName2 += (c2-1);
-       if(!gGeoManager->SetAlignableEntry(strEntryName2.Data(),wafer.Data()))
+       if(!gGeoManager->SetAlignableEntry(strEntryName2.Data(),wafer.Data(),modUID))
          AliFatal("Unable to set alignable entry!!");    
        //printf("%s    ==    %s\n",strEntryName2.Data(),wafer.Data());
 
        // Creates the TGeo Local to Tracking transformation matrix ...
-       TGeoPNEntry *alignableEntry = gGeoManager->GetAlignableEntry(strEntryName2.Data());
-       const char *path = alignableEntry->GetTitle();
-       if (!gGeoManager->cd(path))
-         AliFatal(Form("Volume path %s not valid!",path));
-       TGeoHMatrix* globMatrix = gGeoManager->GetCurrentMatrix();
+       TGeoPNEntry *alignableEntry = gGeoManager->GetAlignableEntryByUID(modUID);
+       TGeoHMatrix *globMatrix = alignableEntry->GetGlobalOrig();
        gtrans = globMatrix->GetTranslation();
        memcpy(&rotMatrix[0], globMatrix->GetRotationMatrix(), 9*sizeof(Double_t));
        al = TMath::ATan2(rotMatrix[1],rotMatrix[0]) + TMath::Pi();
@@ -721,62 +709,7 @@ void AliITSvPPRasymmFMD::AddAlignableVolumes() const
     }
   }
 }
-//______________________________________________________________________
-void AliITSvPPRasymmFMD::BuildGeometry(){
-    //    Geometry builder for the ITS version 10. Event Display geometry.
-    // Inputs:
-    //   none.
-    // Outputs:
-    //   none.
-    // Return:
-    //   none.
 
-    TNode *node, *top;
-    
-    const Int_t kColorITS=kYellow;
-    //
-    top = gAlice->GetGeometry()->GetNode("alice");
-
-
-    new TTUBE("S_layer1","Layer1 of ITS","void",
-             3.8095,3.8095+1.03*9.36/100.,14.35);
-    top->cd();
-    node = new TNode("Layer1","Layer1","S_layer1",0,0,0,"");
-    node->SetLineColor(kColorITS);
-    fNodes->Add(node);
-
-    new TTUBE("S_layer2","Layer2 of ITS","void",7.,7.+1.03*9.36/100.,14.35);
-    top->cd();
-    node = new TNode("Layer2","Layer2","S_layer2",0,0,0,"");
-    node->SetLineColor(kColorITS);
-    fNodes->Add(node);
-
-    new TTUBE("S_layer3","Layer3 of ITS","void",15.,15.+0.94*9.36/100.,25.1);
-    top->cd();
-    node = new TNode("Layer3","Layer3","S_layer3",0,0,0,"");
-    node->SetLineColor(kColorITS);
-    fNodes->Add(node);
-
-    new TTUBE("S_layer4","Layer4 of ITS","void",24.1,24.1+0.95*9.36/100.,32.1);
-    top->cd();
-    node = new TNode("Layer4","Layer4","S_layer4",0,0,0,"");
-    node->SetLineColor(kColorITS);
-    fNodes->Add(node);
-
-    new TTUBE("S_layer5","Layer5 of ITS","void",
-             38.5,38.5+0.91*9.36/100.,49.405);
-    top->cd();
-    node = new TNode("Layer5","Layer5","S_layer5",0,0,0,"");
-    node->SetLineColor(kColorITS);
-    fNodes->Add(node);
-
-    new TTUBE("S_layer6","Layer6 of ITS","void",
-             43.5765,43.5765+0.87*9.36/100.,55.27);
-    top->cd();
-    node = new TNode("Layer6","Layer6","S_layer6",0,0,0,"");
-    node->SetLineColor(kColorITS);
-    fNodes->Add(node);
-}
 //______________________________________________________________________
 void AliITSvPPRasymmFMD::CreateGeometry(){
     //    This routine defines and Creates the geometry for version 10 of 
@@ -5118,11 +5051,6 @@ void AliITSvPPRasymmFMD::CreateGeometry(){
 
   // --- Outputs the geometry tree in the EUCLID/CAD format 
   
-  if (fEuclidOut) {
-    gMC->WriteEuclid("ITSgeometry", "ITSV", 1, 5);
-  }
-
-
   {
     if (!gGeoManager) {
       AliError("TGeoManager doesn't exist !");
@@ -5339,6 +5267,7 @@ void AliITSvPPRasymmFMD::CreateGeometry(){
   }
 
 }
+
 //______________________________________________________________________
 void AliITSvPPRasymmFMD::CreateMaterials(){
     // Create ITS materials
@@ -5357,8 +5286,8 @@ void AliITSvPPRasymmFMD::CreateMaterials(){
     // Return:
     //   none.
 
-    Int_t   ifield = gAlice->Field()->Integ();
-    Float_t fieldm = gAlice->Field()->Max();
+    Int_t   ifield = ((AliMagF*)TGeoGlobalMagField::Instance()->GetField())->Integ();
+    Float_t fieldm = ((AliMagF*)TGeoGlobalMagField::Instance()->GetField())->Max();
 
     Float_t tmaxfd = 0.1; // 1.0; // Degree
     Float_t stemax = 1.0; // cm
@@ -5854,6 +5783,7 @@ void AliITSvPPRasymmFMD::CreateMaterials(){
     AliMaterial(96, "SSD cone$",63.546, 29., 1.15, 1.265, 999);
     AliMedium(96,"SSD cone$",96,0,ifield,fieldm,tmaxfdServ,stemaxServ,deemaxServ,epsilServ,stminServ);
 }
+
 //______________________________________________________________________
 void AliITSvPPRasymmFMD::InitAliITSgeom(){
     //     Based on the geometry tree defined in Geant 3.21, this
@@ -5971,6 +5901,7 @@ void AliITSvPPRasymmFMD::InitAliITSgeom(){
     } // end for lay
     return;
 }
+
 //______________________________________________________________________
 void AliITSvPPRasymmFMD::Init(){
     //     Initialise the ITS after it has been created.
@@ -5985,20 +5916,17 @@ void AliITSvPPRasymmFMD::Init(){
                 fMinorVersion));
     //
     /*  obsolete initialization of AliITSgeom from external "det" file
-    if(fRead[0]=='\0') strncpy(fRead,fEuclidGeomDet,60);
-    if(fWrite[0]=='\0') strncpy(fWrite,fEuclidGeomDet,60);
     AliITSgeom* geom = new AliITSgeom();
     SetITSgeom(geom);
-    if(fGeomDetIn) GetITSgeom()->ReadNewFile(fRead);
     else this->InitAliITSgeom();
     */
     UpdateInternalGeometry();
     AliITS::Init();
-    if(fGeomDetOut) GetITSgeom()->WriteNewFile(fWrite);
 
     //
     fIDMother = gMC->VolId("ITSV"); // ITS Mother Volume ID.
 }
+
 //______________________________________________________________________
 void AliITSvPPRasymmFMD::SetDefaults(){
     // sets the default segmentation, response, digit and raw cluster classes
@@ -6009,98 +5937,20 @@ void AliITSvPPRasymmFMD::SetDefaults(){
     // Return:
     //   none.
 
-    const Float_t kconv = 1.0e+04; // convert cm to microns
-
     if(!fDetTypeSim){
        Warning("SetDefaults","Error fDetTypeSim not defined");
        return;
     }
 
-    AliITSgeomSPD  *s0;
-    AliITSgeomSDD  *s1;
-    AliITSgeomSSD  *s2;
-    Int_t i;
-    Float_t bx[256],bz[280];
-
     fDetTypeSim->SetDefaults();
     
-    //SPD
-    s0 = (AliITSgeomSPD*) GetITSgeom()->GetShape(kSPD);
-    // Get shape info. Do it this way for now.
-    //AliITSCalibrationSPD* resp0=new AliITSCalibrationSPD();
-    AliITSsegmentationSPD* seg0 = 
-       (AliITSsegmentationSPD*)fDetTypeSim->GetSegmentationModel(0);
-    seg0->SetDetSize(s0->GetDx()*2.*kconv, // base this on AliITSgeomSPD
-                    s0->GetDz()*2.*kconv, // for now.
-                    s0->GetDy()*2.*kconv); // x,z,y full width in microns.
-    seg0->SetNPads(256,160);// Number of Bins in x and z
-    for(i=000;i<256;i++) bx[i] =  50.0; // in x all are 50 microns.
-    for(i=000;i<160;i++) bz[i] = 425.0; // most are 425 microns except below
-    for(i=160;i<280;i++) bz[i] =   0.0; // Outside of detector.
-    bz[ 31] = bz[ 32] = 625.0; // first chip boundry
-    bz[ 63] = bz[ 64] = 625.0; // first chip boundry
-    bz[ 95] = bz[ 96] = 625.0; // first chip boundry
-    bz[127] = bz[128] = 625.0; // first chip boundry
-    bz[160] = 425.0; // Set so that there is no zero pixel size for fNz.
-    seg0->SetBinSize(bx,bz); // Based on AliITSgeomSPD for now.
-    SetSegmentationModel(kSPD,seg0);
-    // set digit and raw cluster classes to be used
-    const char *kData0=(fDetTypeSim->GetCalibrationModel(
-                           GetITSgeom()->GetStartSPD()))->DataType();
-    if (strstr(kData0,"real")) fDetTypeSim->SetDigitClassName(kSPD,
-                                                             "AliITSdigit");
-    else fDetTypeSim->SetDigitClassName(kSPD,"AliITSdigitSPD");
-    // SDD
-    s1 = (AliITSgeomSDD*) GetITSgeom()->GetShape(kSDD);
-    // Get shape info. Do it this way for now.
-
-    //AliITSCalibrationSDD* resp1=new AliITSCalibrationSDD("simulated");
-    AliITSsegmentationSDD* seg1 = 
-       (AliITSsegmentationSDD*)fDetTypeSim->GetSegmentationModel(1);
-    seg1->SetDetSize(s1->GetDx()*kconv, // base this on AliITSgeomSDD
-                    s1->GetDz()*2.*kconv, // for now.
-                    s1->GetDy()*2.*kconv); // x,z,y full width in microns.
-    seg1->SetNPads(256,256);// Use AliITSgeomSDD for now
-    SetSegmentationModel(kSDD,seg1);
-    const char *kData1=(fDetTypeSim->GetCalibrationModel(
-                           GetITSgeom()->GetStartSDD()))->DataType();
-    AliITSCalibrationSDD* rsp = 
-       (AliITSCalibrationSDD*)fDetTypeSim->GetCalibrationModel(
-           GetITSgeom()->GetStartSDD());
-    const char *kopt=rsp->GetZeroSuppOption();
-    if((!strstr(kopt,"2D")) && (!strstr(kopt,"1D")) || strstr(kData1,"real") ){
-       fDetTypeSim->SetDigitClassName(kSDD,"AliITSdigit");
-    } else fDetTypeSim->SetDigitClassName(kSDD,"AliITSdigitSDD");
-    // SSD  Layer 5
-
-    s2 = (AliITSgeomSSD*) GetITSgeom()->GetShape(kSSD);
-    // Get shape info. Do it this way for now.
-
-
-    //SetCalibrationModel(GetITSgeom()->GetStartSSD(),
-    // new AliITSCalibrationSSD("simulated"));
-    AliITSsegmentationSSD* seg2 = 
-       (AliITSsegmentationSSD*)fDetTypeSim->GetSegmentationModel(2);
-    seg2->SetDetSize(s2->GetDx()*2.*kconv, // base this on AliITSgeomSSD
-                    s2->GetDz()*2.*kconv, // for now.
-                    s2->GetDy()*2.*kconv); // x,z,y full width in microns.
-    seg2->SetPadSize(95.,0.); // strip x pitch in microns
-    seg2->SetNPads(768,0); // number of strips on each side.
-    seg2->SetAngles(0.0075,0.0275); // strip angels rad P and N side.
-    seg2->SetAnglesLay5(0.0075,0.0275); // strip angels rad P and N side.
-    seg2->SetAnglesLay6(0.0275,0.0075); // strip angels rad P and N side.
-    SetSegmentationModel(kSSD,seg2); 
-        const char *kData2=(fDetTypeSim->GetCalibrationModel(
-                               GetITSgeom()->GetStartSSD()))->DataType();
-    if(strstr(kData2,"real") ) fDetTypeSim->SetDigitClassName(kSSD,
-                                                             "AliITSdigit");
-    else fDetTypeSim->SetDigitClassName(kSSD,"AliITSdigitSSD");
     if(fgkNTYPES>3){
        Warning("SetDefaults",
                "Only the four basic detector types are initialised!");
     }// end if
     return;
 }
+
 //______________________________________________________________________
 void AliITSvPPRasymmFMD::DrawModule() const{
     //     Draw a shaded view of the FMD version 10.
@@ -6153,6 +6003,7 @@ void AliITSvPPRasymmFMD::DrawModule() const{
     gMC->Gdhead(1111, "Inner Tracking System Version 1");
     gMC->Gdman(17, 6, "MAN");
 }
+
 /*
 //______________________________________________________________________
 void AliITSvPPRasymmFMD::StepManager(){
@@ -6289,6 +6140,7 @@ void AliITSvPPRasymmFMD::StepManager(){
 }
 
 */
+
 //______________________________________________________________________
 void AliITSvPPRasymmFMD::StepManager(){
     //    Called for every step in the ITS, then calles the AliITShit class
@@ -6453,6 +6305,7 @@ void AliITSvPPRasymmFMD::StepManager(){
     */
     return;
 }
+
 //______________________________________________________________________
 void AliITSvPPRasymmFMD::PrintAscii(ostream *os)const{
     // Print out class data values in Ascii Form to output stream
@@ -6476,17 +6329,15 @@ void AliITSvPPRasymmFMD::PrintAscii(ostream *os)const{
 #endif
 #endif
 
-    *os << fGeomDetOut << " " << fGeomDetIn << " " << fByThick << " ";
+    *os << fByThick << " ";
     *os << fMajorVersion << " " << fMinorVersion << " ";
-    *os << "\"" << fEuclidGeomDet.Data() << "\"" << " ";
-    *os << "\"" << fRead.Data() << "\"" << " ";
-    *os << "\"" << fWrite.Data() << "\"" << " ";
     fmt = os->setf(ios::scientific); // set scientific floating point output
     *os << fDet1 << " "  << fDet2 << " " << fChip1 << " " << fChip2 << " ";
     *os << fRails << " " << fFluid << " " << fIDMother;
     os->flags(fmt); // reset back to old Formating.
     return;
 }
+
 //______________________________________________________________________
 void AliITSvPPRasymmFMD::ReadAscii(istream *is){
     // Read in class data values in Ascii Form to output stream
@@ -6496,22 +6347,16 @@ void AliITSvPPRasymmFMD::ReadAscii(istream *is){
     //   none.
     // Return:
     //   none.
-    Char_t name[120];
 
-    *is >> fGeomDetOut >> fGeomDetIn >> fByThick;
+    *is >> fByThick;
     *is >> fMajorVersion >> fMinorVersion;
-    *is >> name;
-    fEuclidGeomDet = name;
-    *is >> name;
-    fRead = name;
-    *is >> name;
-    fWrite = name;
     *is >> fDet1 >> fDet2 >> fChip1 >> fChip2;
     *is >> fRails >> fFluid >> fIDMother;
     fIgm.SetVersion((AliITSVersion_t)fMajorVersion,fMinorVersion);
     fIgm.SetGeometryName("ITS PPR aymmetric services with course"
                          " cables on cones");
 }
+
 //______________________________________________________________________
 ostream &operator<<(ostream &os,const AliITSvPPRasymmFMD &s){
     // Standard output streaming function
@@ -6526,6 +6371,7 @@ ostream &operator<<(ostream &os,const AliITSvPPRasymmFMD &s){
     s.PrintAscii(&os);
     return os;
 }
+
 //______________________________________________________________________
 istream &operator>>(istream &is,AliITSvPPRasymmFMD &s){
     // Standard inputput streaming function