Using symbolic names instead of volume paths (Raffaele)
authorhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 17 Oct 2006 15:59:42 +0000 (15:59 +0000)
committerhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 17 Oct 2006 15:59:42 +0000 (15:59 +0000)
12 files changed:
EMCAL/MakeEMCALFullMisAlignment.C
EMCAL/MakeEMCALResMisAlignment.C
ITS/MakeITSFullMisAlignment.C
ITS/MakeITSResMisAlignment.C
PHOS/MakePHOSFullMisAlignment.C
PHOS/MakePHOSResMisAlignment.C
STEER/AliAlignObj.cxx
TOF/MakeTOFFullMisAlignment.C
TPC/MakeTPCFullMisAlignment.C
TPC/MakeTPCResMisAlignment.C
TRD/MakeTRDFullMisAlignment.C
TRD/MakeTRDResMisAlignment.C

index 0d19e90..ca50227 100644 (file)
@@ -16,7 +16,7 @@ void MakeEMCALFullMisAlignment(){
 
   // null shifts and rotations
 
-  const TString basepath = "ALIC_1/XEN1_1/SMOD_";
+  const TString basepath = "EMCAL/FullSupermodule";
   TString pathstr;
 
   Int_t iIndex=0; //let all modules have index=0 in a layer with no LUT
index 7c49a2f..7fe39f1 100644 (file)
@@ -11,7 +11,7 @@ void MakeEMCALResMisAlignment(){
 
   Double_t dx, dy, dz, dpsi, dtheta, dphi;
 
-  const TString basepath = "ALIC_1/XEN1_1/SMOD_";
+  const TString basepath = "EMCAL/FullSupermodule";
   TString pathstr;
 
   Int_t iIndex=0; // let all modules have index=0 in a layer with no LUT
index 4563364..32d3639 100644 (file)
@@ -23,7 +23,7 @@ void MakeITSFullMisAlignment(){
   Double_t dx=0., dy=0., dz=0., dpsi=0., dtheta=0., dphi=0.;
 
   Int_t j = 0;
-  new(alobj[j]) AliAlignObjAngles("ALIC_1/ITSV_1", 0, dx, dy, globalZ, dpsi, dtheta, dphi, kTRUE);
+  new(alobj[j]) AliAlignObjAngles("ITS", 0, dx, dy, globalZ, dpsi, dtheta, dphi, kTRUE);
   j++;
 
   for ( Int_t l = AliAlignObj::kSPD1; l <= AliAlignObj::kSSD2; l++) {
@@ -76,9 +76,9 @@ void MakeITSFullMisAlignment(){
       default: Printf("Wrong layer index in ITS (%d) !",l);
       };
       UShort_t volid = AliAlignObj::LayerToVolUID(iLayer,iModule);
-      const char *path = AliAlignObj::SymName(volid);
+      const char *symname = AliAlignObj::SymName(volid);
 
-      new(alobj[j]) AliAlignObjAngles(path, volid, dx, dy, dz, dpsi, dtheta, dphi, kFALSE);
+      new(alobj[j]) AliAlignObjAngles(symname, volid, dx, dy, dz, dpsi, dtheta, dphi, kFALSE);
       j++;
 
     }
index cfbd0d8..bd647ba 100644 (file)
@@ -24,7 +24,8 @@ void MakeITSResMisAlignment(){
   AliAlignObj::ELayerID iLayer = AliAlignObj::kInvalidLayer; 
 
   Int_t j = 0;
-  new(alobj[j]) AliAlignObjAngles("ALIC_1/ITSV_1", 0, dx, dy, globalZ, dpsi, dtheta, dphi, kTRUE);
+
+  new(alobj[j]) AliAlignObjAngles("ITS", 0, dx, dy, globalZ, dpsi, dtheta, dphi, kTRUE);
   j++;
 
   for ( Int_t l = AliAlignObj::kSPD1; l <= AliAlignObj::kSSD2; l++) {
@@ -78,9 +79,9 @@ void MakeITSResMisAlignment(){
       default: Printf("Wrong layer index in ITS (%d) !",l);
       };
       UShort_t volid = AliAlignObj::LayerToVolUID(iLayer,iModule);
-      const char *path = AliAlignObj::SymName(volid);
+      const char *symname = AliAlignObj::SymName(volid);
 
-      new(alobj[j]) AliAlignObjAngles(path, volid, dx, dy, dz, dpsi, dtheta, dphi, kFALSE);
+      new(alobj[j]) AliAlignObjAngles(symname, volid, dx, dy, dz, dpsi, dtheta, dphi, kFALSE);
       j++;
 
     }
index 9e07379..723b417 100644 (file)
@@ -15,31 +15,31 @@ void MakePHOSFullMisAlignment(){
 
 
   // Alignment for 5 PHOS modules
-  new(alobj[0]) AliAlignObjAngles("/ALIC_1/PHOS_1",
+  new(alobj[0]) AliAlignObjAngles("PHOS/Module1",
          volid, -20., -10.,   0., dpsi, dtheta, 5, kTRUE);
-  new(alobj[1]) AliAlignObjAngles("/ALIC_1/PHOS_2",
+  new(alobj[1]) AliAlignObjAngles("PHOS/Module2",
          volid, -10.,   0., -10., dpsi, dtheta, 2, kTRUE);
-  new(alobj[2]) AliAlignObjAngles("/ALIC_1/PHOS_3",
+  new(alobj[2]) AliAlignObjAngles("PHOS/Module3",
          volid,   5., -10.,  10., dpsi, dtheta, 0, kTRUE);
-  new(alobj[3]) AliAlignObjAngles("/ALIC_1/PHOS_4",
+  new(alobj[3]) AliAlignObjAngles("PHOS/Module4",
          volid, +10.,  -0., -10., dpsi, dtheta, 2, kTRUE);
-  new(alobj[4]) AliAlignObjAngles("/ALIC_1/PHOS_5",
+  new(alobj[4]) AliAlignObjAngles("PHOS/Module5",
          volid, +20., -10.,   0., dpsi, dtheta, 5, kTRUE);
 
   // Alignment for PHOS cradle
-  new(alobj[5]) AliAlignObjAngles("/ALIC_1/PCRA_0",
+  new(alobj[5]) AliAlignObjAngles("PHOS/Cradle0",
          volid, 0., 0., -displacement, dpsi, dtheta, dphi, kTRUE);
-  new(alobj[6]) AliAlignObjAngles("/ALIC_1/PCRA_1",
+  new(alobj[6]) AliAlignObjAngles("PHOS/Cradle1",
          volid, 0., 0., +displacement, dpsi, dtheta, dphi, kTRUE);
 
   // Alignment for cradle wheels
-  new(alobj[7]) AliAlignObjAngles("/ALIC_1/PWHE_0",
+  new(alobj[7]) AliAlignObjAngles("PHOS/Wheel0",
          volid, 0., 0., -displacement, dpsi, dtheta, dphi, kTRUE);
-  new(alobj[8]) AliAlignObjAngles("/ALIC_1/PWHE_1",
+  new(alobj[8]) AliAlignObjAngles("PHOS/Wheel1",
          volid, 0., 0., -displacement, dpsi, dtheta, dphi, kTRUE);
-  new(alobj[9]) AliAlignObjAngles("/ALIC_1/PWHE_2",
+  new(alobj[9]) AliAlignObjAngles("PHOS/Wheel2",
          volid, 0., 0., +displacement, dpsi, dtheta, dphi, kTRUE);
-  new(alobj[10]) AliAlignObjAngles("/ALIC_1/PWHE_3",
+  new(alobj[10]) AliAlignObjAngles("PHOS/Wheel3",
          volid, 0., 0., +displacement, dpsi, dtheta, dphi, kTRUE);
 
   // *************************    2nd step    ***************
index 01aad31..27b0208 100644 (file)
@@ -14,31 +14,31 @@ void MakePHOSResMisAlignment(){
   UShort_t volid = AliAlignObj::LayerToVolUID(iLayer,iIndex);
 
   // Alignment for 5 PHOS modules
-  new(alobj[0]) AliAlignObjAngles("/ALIC_1/PHOS_1",
+  new(alobj[0]) AliAlignObjAngles("PHOS/Module1",
         volid, -0.20, -0.1, +0.0, dpsi, dtheta, 0.2, kTRUE);
-  new(alobj[1]) AliAlignObjAngles("/ALIC_1/PHOS_2",
+  new(alobj[1]) AliAlignObjAngles("PHOS/Module2",
        volid, -0.10, +0.0, -0.2, dpsi, dtheta, 0.2, kTRUE);
-  new(alobj[2]) AliAlignObjAngles("/ALIC_1/PHOS_3",
+  new(alobj[2]) AliAlignObjAngles("PHOS/Module3",
        volid,  0.05, -0.1,  0.2, dpsi, dtheta, 0.0, kTRUE);
-  new(alobj[3]) AliAlignObjAngles("/ALIC_1/PHOS_4",
+  new(alobj[3]) AliAlignObjAngles("PHOS/Module4",
        volid, +0.10, -0.0, -0.1, dpsi, dtheta, 0.1, kTRUE);
-  new(alobj[4]) AliAlignObjAngles("/ALIC_1/PHOS_5",
+  new(alobj[4]) AliAlignObjAngles("PHOS/Module5",
        volid, +0.20, -0.1,  0.1, dpsi, dtheta, 0.2, kTRUE);
 
   // Alignment for PHOS cradle
-  new(alobj[5]) AliAlignObjAngles("/ALIC_1/PCRA_0",
+  new(alobj[5]) AliAlignObjAngles("PHOS/Cradle0",
        volid, 0., 0., -displacement, dpsi, dtheta, dphi, kTRUE);
-  new(alobj[6]) AliAlignObjAngles("/ALIC_1/PCRA_1",
+  new(alobj[6]) AliAlignObjAngles("PHOS/Cradle1",
        volid, 0., 0., +displacement, dpsi, dtheta, dphi, kTRUE);
 
   // Alignment for cradle wheels
-  new(alobj[7])  AliAlignObjAngles("/ALIC_1/PWHE_0",
+  new(alobj[7])  AliAlignObjAngles("PHOS/Wheel0",
        volid, 0., 0., -displacement, dpsi, dtheta, dphi, kTRUE);
-  new(alobj[8])  AliAlignObjAngles("/ALIC_1/PWHE_1",
+  new(alobj[8])  AliAlignObjAngles("PHOS/Wheel1",
        volid, 0., 0., -displacement, dpsi, dtheta, dphi, kTRUE);
-  new(alobj[9])  AliAlignObjAngles("/ALIC_1/PWHE_2",
+  new(alobj[9])  AliAlignObjAngles("PHOS/Wheel2",
        volid, 0., 0., +displacement, dpsi, dtheta, dphi, kTRUE);
-  new(alobj[10]) AliAlignObjAngles("/ALIC_1/PWHE_3",
+  new(alobj[10]) AliAlignObjAngles("PHOS/Wheel3",
        volid, 0., 0., +displacement, dpsi, dtheta, dphi, kTRUE);
 
 
index 7ca7f2a..a6790ef 100644 (file)
  **************************************************************************/
 
 //-----------------------------------------------------------------
-//   Implementation of the alignment object class through the abstract
-//  class AliAlignObj. From it two derived concrete representation of
-//  alignment object class (AliAlignObjAngles, AliAlignObjMatrix) are
-//  derived in separate files.
+//  Implementation of the alignment object class, holding the alignment
+//  constants for a single volume, through the abstract class AliAlignObj.
+//  From it two derived concrete representation of alignment object class
+//  (AliAlignObjAngles, AliAlignObjMatrix) are derived in separate files.
 //-----------------------------------------------------------------
 #include <TGeoManager.h>
 #include <TGeoPhysicalNode.h>
@@ -153,10 +153,12 @@ void AliAlignObj::SetVolUID(ELayerID detId, Int_t modId)
 //_____________________________________________________________________________
 void AliAlignObj::GetVolUID(ELayerID &layerId, Int_t &modId) const
 {
-  // From detector name and module number (according to detector numbering)
-  // build fVolUID, unique numerical identity of that volume inside ALICE
-  // fVolUID is 16 bits, first 5 reserved for detID (32 possible values),
-  // remaining 11 for module ID inside det (2048 possible values).
+  // From the fVolUID, unique numerical identity of that volume inside ALICE,
+  // (voluid is 16 bits, first 5 reserved for layerID (32 possible values),
+  // remaining 11 for module ID inside det (2048 possible values)), sets
+  // the argument layerId to the identity of the layer to which that volume
+  // belongs and sets the argument modId to the identity of that volume
+  // internally to the layer.
   //
   layerId = VolUIDToLayer(fVolUID,modId);
 }
@@ -201,6 +203,7 @@ Int_t AliAlignObj::Compare(const TObject *obj) const
   // Used in the sorting during
   // the application of alignment
   // objects to the geometry
+  //
   Int_t level = GetLevel();
   Int_t level2 = ((AliAlignObj *)obj)->GetLevel();
   if (level == level2)
@@ -214,6 +217,7 @@ void AliAlignObj::AnglesToMatrix(const Double_t *angles, Double_t *rot) const
 {
   // Calculates the rotation matrix using the 
   // Euler angles in "x y z" notation
+  //
   Double_t degrad = TMath::DegToRad();
   Double_t sinpsi = TMath::Sin(degrad*angles[0]);
   Double_t cospsi = TMath::Cos(degrad*angles[0]);
@@ -240,6 +244,7 @@ Bool_t AliAlignObj::MatrixToAngles(const Double_t *rot, Double_t *angles) const
   // using the rotation matrix
   // Returns false in case the rotation angles can not be
   // extracted from the matrix
+  //
   if(TMath::Abs(rot[0])<1e-7 || TMath::Abs(rot[8])<1e-7) {
     AliError("Failed to extract roll-pitch-yall angles!");
     return kFALSE;
@@ -258,7 +263,7 @@ void AliAlignObj::Transform(AliTrackPoint &p) const
   // transformation matrix provided by the AliAlignObj
   // The covariance matrix is not affected since we assume
   // that the transformations are sufficiently small
-
+  //
   if (fVolUID != p.GetVolumeID())
     AliWarning(Form("Alignment object ID is not equal to the space-point ID (%d != %d)",fVolUID,p.GetVolumeID())); 
 
@@ -283,6 +288,7 @@ void AliAlignObj::Transform(AliTrackPointArray &array) const
 {
   // This method is used to transform all the track points
   // from the input AliTrackPointArray
+  // 
   AliTrackPoint p;
   for (Int_t i = 0; i < array.GetNPoints(); i++) {
     array.GetPoint(p,i);
@@ -297,6 +303,7 @@ void AliAlignObj::Print(Option_t *) const
   // Print the contents of the
   // alignment object in angles and
   // matrix representations
+  //
   Double_t tr[3];
   GetTranslation(tr);
   Double_t angles[3];
@@ -321,8 +328,9 @@ void AliAlignObj::Print(Option_t *) const
 //_____________________________________________________________________________
 Int_t AliAlignObj::LayerSize(Int_t layerId)
 {
-  // Get the corresponding layer size.
+  // Get the layer size for layer corresponding to layerId.
   // Implemented only for ITS,TPC,TRD,TOF and RICH
+  //
   if (layerId < kFirstLayer || layerId >= kLastLayer) {
     AliErrorClass(Form("Invalid layer index %d ! Layer range is (%d -> %d) !",layerId,kFirstLayer,kLastLayer));
     return 0;
@@ -335,8 +343,9 @@ Int_t AliAlignObj::LayerSize(Int_t layerId)
 //_____________________________________________________________________________
 const char* AliAlignObj::LayerName(Int_t layerId)
 {
-  // Get the corresponding layer name.
+  // Get the layer name corresponding to layerId.
   // Implemented only for ITS,TPC,TRD,TOF and RICH
+  //
   if (layerId < kFirstLayer || layerId >= kLastLayer) {
     AliErrorClass(Form("Invalid layer index %d ! Layer range is (%d -> %d) !",layerId,kFirstLayer,kLastLayer));
     return "Invalid Layer!";
@@ -349,8 +358,9 @@ const char* AliAlignObj::LayerName(Int_t layerId)
 //_____________________________________________________________________________
 UShort_t AliAlignObj::LayerToVolUID(ELayerID layerId, Int_t modId)
 {
-  // From detector (layer) name and module number (according to detector numbering)
-  // build fVolUID, unique numerical identity of that volume inside ALICE
+  // From detector (layer) name and module number (according to detector
+  // internal numbering) build the unique numerical identity of that volume
+  // inside ALICE
   // fVolUID is 16 bits, first 5 reserved for layerID (32 possible values),
   // remaining 11 for module ID inside det (2048 possible values).
   //
@@ -360,8 +370,9 @@ UShort_t AliAlignObj::LayerToVolUID(ELayerID layerId, Int_t modId)
 //_____________________________________________________________________________
 UShort_t AliAlignObj::LayerToVolUID(Int_t   layerId, Int_t modId)
 {
-  // From detector (layer) index and module number (according to detector numbering)
-  // build fVolUID, unique numerical identity of that volume inside ALICE
+  // From detector (layer) name and module number (according to detector
+  // internal numbering) build the unique numerical identity of that volume
+  // inside ALICE
   // fVolUID is 16 bits, first 5 reserved for layerID (32 possible values),
   // remaining 11 for module ID inside det (2048 possible values).
   //
@@ -371,10 +382,11 @@ UShort_t AliAlignObj::LayerToVolUID(Int_t   layerId, Int_t modId)
 //_____________________________________________________________________________
 AliAlignObj::ELayerID AliAlignObj::VolUIDToLayer(UShort_t voluid, Int_t &modId)
 {
-  // From detector (layer) name and module number (according to detector numbering)
-  // build fVolUID, unique numerical identity of that volume inside ALICE
-  // fVolUID is 16 bits, first 5 reserved for layerID (32 possible values),
-  // remaining 11 for module ID inside det (2048 possible values).
+  // From voluid, unique numerical identity of that volume inside ALICE,
+  // (voluid is 16 bits, first 5 reserved for layerID (32 possible values),
+  // remaining 11 for module ID inside det (2048 possible values)), return
+  // the identity of the layer to which that volume belongs and sets the
+  // argument modId to the identity of that volume internally to the layer.
   //
   modId = voluid & 0x7ff;
 
@@ -384,10 +396,10 @@ AliAlignObj::ELayerID AliAlignObj::VolUIDToLayer(UShort_t voluid, Int_t &modId)
 //_____________________________________________________________________________
 AliAlignObj::ELayerID AliAlignObj::VolUIDToLayer(UShort_t voluid)
 {
-  // From detector (layer) name and module number (according to detector numbering)
-  // build fVolUID, unique numerical identity of that volume inside ALICE
-  // fVolUID is 16 bits, first 5 reserved for layerID (32 possible values),
-  // remaining 11 for module ID inside det (2048 possible values).
+  // From voluid, unique numerical identity of that volume inside ALICE,
+  // (voluid is 16 bits, first 5 reserved for layerID (32 possible values),
+  // remaining 11 for module ID inside det (2048 possible values)), return
+  // the identity of the layer to which that volume belongs
   //
   return ELayerID((voluid >> 11) & 0x1f);
 }
@@ -396,8 +408,9 @@ AliAlignObj::ELayerID AliAlignObj::VolUIDToLayer(UShort_t voluid)
 void AliAlignObj::SetPars(Double_t x, Double_t y, Double_t z,
                          Double_t psi, Double_t theta, Double_t phi)
 {
-  // Set rotation matrix and translation
-  // using 3 angles and 3 translations
+  // Set rotation matrix and translation using 3 angles and 3 translations
+  // The three angles are expressed in degrees
+  // 
   SetTranslation(x,y,z);
   SetRotation(psi,theta,phi);
 }
@@ -406,11 +419,12 @@ void AliAlignObj::SetPars(Double_t x, Double_t y, Double_t z,
 Bool_t AliAlignObj::SetLocalPars(Double_t x, Double_t y, Double_t z,
                                 Double_t psi, Double_t theta, Double_t phi)
 {
-  // Set the translations and angles by using parameters
-  // defined in the local (in TGeo means) coordinate system
-  // of the alignable volume. In case that the TGeo was
-  // initialized, returns false and the object parameters are
-  // not set.
+  // Set the translations and angles (in degrees) by considering the
+  // parameters passed as arguments as expressed in the local reference
+  // system of the alignable volume (known by TGeo geometry).
+  // In case that the TGeo was not initialized or not closed,
+  // returns false and the object parameters are not set.
+  //
   TGeoHMatrix m;
   Double_t tr[3] = {x, y, z};
   m.SetTranslation(tr);
@@ -426,12 +440,12 @@ Bool_t AliAlignObj::SetLocalPars(Double_t x, Double_t y, Double_t z,
 //_____________________________________________________________________________
 Bool_t AliAlignObj::SetLocalMatrix(const TGeoMatrix& m)
 {
-  // Set the translations and angles by using TGeo matrix
-  // defined in the local (in TGeo means) coordinate system
-  // of the alignable volume. In case that the TGeo was
-  // initialized, returns false and the object parameters are
-  // not set.
-
+  // Set the translations and angles by considering the TGeo matrix
+  // passed as argument as expressing the transformation in the local
+  // reference system of the alignable volume (known by TGeo geometry).
+  // In case that the TGeo was not initialized or not closed,
+  // returns false and the object parameters are not set.
+  //
   if (!gGeoManager || !gGeoManager->IsClosed()) {
     AliError("Can't set the alignment object parameters! gGeoManager doesn't exist or it is still opened!");
     return kFALSE;
@@ -472,8 +486,10 @@ Bool_t AliAlignObj::SetLocalMatrix(const TGeoMatrix& m)
 //_____________________________________________________________________________
 Bool_t AliAlignObj::SetMatrix(const TGeoMatrix& m)
 {
-  // Set rotation matrix and translation
-  // using TGeoMatrix
+  // Set rotation matrix and translation using the TGeoMatrix passed
+  // as argument considering it as relative to the global reference
+  // system
+  //
   SetTranslation(m);
   return SetRotation(m);
 }
@@ -481,9 +497,10 @@ Bool_t AliAlignObj::SetMatrix(const TGeoMatrix& m)
 //_____________________________________________________________________________
 Bool_t AliAlignObj::ApplyToGeometry()
 {
-  // Apply the current alignment object
-  // to the TGeo geometry
-
+  // Apply the current alignment object to the TGeo geometry
+  // This method returns FALSE if the symname of the object was not
+  // valid neither to get a TGeoPEntry nor as a volume path
+  //
   if (!gGeoManager || !gGeoManager->IsClosed()) {
     AliError("Can't apply the alignment object! gGeoManager doesn't exist or it is still opened!");
     return kFALSE;
@@ -494,18 +511,22 @@ Bool_t AliAlignObj::ApplyToGeometry()
   TGeoPhysicalNode* node;
   TGeoPNEntry* pne = gGeoManager->GetAlignableEntry(symname);
   if(pne){
-    node = gGeoManager->MakeAlignablePN(pne);
-    if(!node) return kFALSE;
     path = pne->GetTitle();
+    if(!gGeoManager->CheckPath(path)){
+      AliDebug(1,Form("Valid PNEntry but invalid volume path %s!",path));
+          // this should happen only for volumes in disactivated branches
+      return kTRUE;
+    }
+    node = gGeoManager->MakeAlignablePN(pne);
   }else{
     AliWarning(Form("The symbolic volume name %s does not correspond to a physical entry. Using it as a volume path!",symname));
     path=symname;
-    if (gGeoManager->GetListOfPhysicalNodes()->FindObject(path)) {
-      AliError(Form("Volume %s has already been misaligned!",path));
+    if (!gGeoManager->CheckPath(path)) {
+      AliError(Form("Volume path %s not valid!",path));
       return kFALSE;
     }
-    if (!gGeoManager->cd(path)) {
-      AliError(Form("Volume path %s not valid!",path));
+    if (gGeoManager->GetListOfPhysicalNodes()->FindObject(path)) {
+      AliError(Form("Volume %s has already been misaligned!",path));
       return kFALSE;
     }
     node = (TGeoPhysicalNode*) gGeoManager->MakePhysicalNode(path);
@@ -540,6 +561,9 @@ Bool_t AliAlignObj::GetFromGeometry(const char *symname, AliAlignObj &alobj)
   // symname (in case equal to the TGeo volume path)
   // The method is extremely slow due to the searching by string.
   // Therefore it should be used with great care!!
+  // This method returns FALSE if the symname of the object was not
+  // valid neither to get a TGeoPEntry nor as a volume path, or if the path
+  // associated to the TGeoPNEntry was not valid.
   //
 
   // Reset the alignment object
@@ -624,6 +648,7 @@ void  AliAlignObj::InitAlignObjFromGeometry()
 //_____________________________________________________________________________
 AliAlignObj* AliAlignObj::GetAlignObj(UShort_t voluid) {
   // Returns the alignment object for given volume ID
+  //
   Int_t modId;
   ELayerID layerId = VolUIDToLayer(voluid,modId);
   return GetAlignObj(layerId,modId);
@@ -632,7 +657,8 @@ AliAlignObj* AliAlignObj::GetAlignObj(UShort_t voluid) {
 //_____________________________________________________________________________
 AliAlignObj* AliAlignObj::GetAlignObj(ELayerID layerId, Int_t modId)
 {
-  // Returns pointer to alignment object givent its layer and module ID
+  // Returns pointer to alignment object given its layer and module ID
+  //
   if(modId<0 || modId>=fgLayerSize[layerId-kFirstLayer]){
     AliWarningClass(Form("Module number %d not in the valid range (0->%d) !",modId,fgLayerSize[layerId-kFirstLayer]-1));
     return NULL;
@@ -644,7 +670,8 @@ AliAlignObj* AliAlignObj::GetAlignObj(ELayerID layerId, Int_t modId)
 
 //_____________________________________________________________________________
 const char* AliAlignObj::SymName(UShort_t voluid) {
-  // Returns the volume path for given volume ID
+  // Returns the symbolic volume name for given volume ID
+  //
   Int_t modId;
   ELayerID layerId = VolUIDToLayer(voluid,modId);
   return SymName(layerId,modId);
@@ -653,7 +680,9 @@ const char* AliAlignObj::SymName(UShort_t voluid) {
 //_____________________________________________________________________________
 const char* AliAlignObj::SymName(ELayerID layerId, Int_t modId)
 {
-  // Returns volume path to alignment object givent its layer and module ID
+  // Returns the symbolic volume name given for a given layer
+  // and module ID
+  //
   if(modId<0 || modId>=fgLayerSize[layerId-kFirstLayer]){
     AliWarningClass(Form("Module number %d not in the valid range (0->%d) !",modId,fgLayerSize[layerId-kFirstLayer]-1));
     return NULL;
@@ -666,40 +695,48 @@ const char* AliAlignObj::SymName(ELayerID layerId, Int_t modId)
 //_____________________________________________________________________________
 void AliAlignObj::InitSymNames()
 {
-  // Initialize the LUTs which contain
-  // the TGeo volume paths for each
-  // alignable volume. The LUTs are
-  // static, so they are created during
-  // the creation of the first intance
-  // of AliAlignObj
-
+  // Initialize the LUTs which associate the symbolic volume names
+  // for each alignable volume with their unique numerical identity.
+  // The LUTs are static, so they are created during the instantiation
+  // of the first intance of AliAlignObj
+  //
   if (fgVolPath[0]) return;
 
   for (Int_t iLayer = 0; iLayer < (kLastLayer - kFirstLayer); iLayer++)
     fgVolPath[iLayer] = new TString[fgLayerSize[iLayer]];
 
+  TString symname;
+  Int_t modnum; // in the following, set it to 0 at the start of each layer
+
+  /*********************       ITS layers  ***********************/
+  TString strSPD = "ITS/SPD";
+  TString strSDD = "ITS/SDD";
+  TString strSSD = "ITS/SSD";
+  TString strStave = "/Stave";
+  TString strLadder = "/Ladder";
+  TString strSector = "/Sector";
+  TString strSensor = "/Sensor";
+  TString strEntryName1;
+  TString strEntryName2;
+
   /*********************       SPD layer1  ***********************/
   {
-    Int_t modnum = 0;
-    TString str0 = "ALIC_1/ITSV_1/ITSD_1/IT12_1/I12B_"; //".../I12A_"
-    TString str1 = "/I10B_";    //"/I10A_";
-    TString str2 = "/I107_";    //"/I103_"
-    //    TString str3 = "/I101_1/ITS1_1";
-    TString symname, symname1, symname2;
+    modnum = 0;
 
     for(Int_t c1 = 1; c1<=10; c1++){
-      symname = str0;
-      symname += c1;
-      symname += str1;
+      strEntryName1 = strSPD;
+      strEntryName1 += 0;
+      strEntryName1 += strSector;
+      strEntryName1 += (c1-1);
       for(Int_t c2 =1; c2<=2; c2++){
-       symname1 = symname;
-       symname1 += c2;
-       symname1 += str2;
+       strEntryName2 = strEntryName1;
+       strEntryName2 += strStave;
+       strEntryName2 += (c2-1);
        for(Int_t c3 =1; c3<=4; c3++){
-         symname2 = symname1;
-         symname2 += c3;
-         //      symname2 += str3;
-         fgVolPath[kSPD1-kFirstLayer][modnum] = symname2.Data();
+         symname = strEntryName2;
+         symname += strLadder;
+         symname += (c3-1);
+         fgVolPath[kSPD1-kFirstLayer][modnum] = symname.Data();
          modnum++;
        }
       }
@@ -708,26 +745,22 @@ void AliAlignObj::InitSymNames()
   
   /*********************       SPD layer2  ***********************/
   {
-    Int_t modnum = 0;
-    TString str0 = "ALIC_1/ITSV_1/ITSD_1/IT12_1/I12B_";  //".../I12A_"
-    TString str1 = "/I20B_";  //"/I20A"
-    TString str2 = "/I1D7_";  //"/I1D3"
-    //    TString str3 = "/I1D1_1/ITS2_1";
-    TString symname, symname1, symname2;
+    modnum = 0;
 
     for(Int_t c1 = 1; c1<=10; c1++){
-      symname = str0;
-      symname += c1;
-      symname += str1;
+      strEntryName1 = strSPD;
+      strEntryName1 += 1;
+      strEntryName1 += strSector;
+      strEntryName1 += (c1-1);
       for(Int_t c2 =1; c2<=4; c2++){
-       symname1 = symname;
-       symname1 += c2;
-       symname1 += str2;
+       strEntryName2 = strEntryName1;
+       strEntryName2 += strStave;
+       strEntryName2 += (c2-1);
        for(Int_t c3 =1; c3<=4; c3++){
-         symname2 = symname1;
-         symname2 += c3;
-         //      symname2 += str3;
-         fgVolPath[kSPD2-kFirstLayer][modnum] = symname2.Data();
+         symname = strEntryName2;
+         symname += strLadder;
+         symname += (c3-1);
+         fgVolPath[kSPD2-kFirstLayer][modnum] = symname.Data();
          modnum++;
        }
       }
@@ -736,21 +769,18 @@ void AliAlignObj::InitSymNames()
 
   /*********************       SDD layer1  ***********************/
   {
-    Int_t modnum=0;
-    TString str0 = "ALIC_1/ITSV_1/ITSD_1/IT34_1/I004_";
-    TString str1 = "/I302_";
-    //    TString str2 = "/ITS3_1";
-    TString symname, symname1;
+    modnum=0;
 
     for(Int_t c1 = 1; c1<=14; c1++){
-      symname = str0;
-      symname += c1;
-      symname += str1;
+      strEntryName1 = strSDD;
+      strEntryName1 += 2;
+      strEntryName1 +=strLadder;
+      strEntryName1 += (c1-1);
       for(Int_t c2 =1; c2<=6; c2++){
-       symname1 = symname;
-       symname1 += c2;
-       //      symname1 += str2;
-       fgVolPath[kSDD1-kFirstLayer][modnum] = symname1.Data();
+       symname = strEntryName1;
+       symname += strSensor;
+       symname += (c2-1);
+       fgVolPath[kSDD1-kFirstLayer][modnum] = symname.Data();
        modnum++;
       }
     }
@@ -758,21 +788,18 @@ void AliAlignObj::InitSymNames()
 
   /*********************       SDD layer2  ***********************/
   {
-    Int_t modnum=0;
-    TString str0 = "ALIC_1/ITSV_1/ITSD_1/IT34_1/I005_";
-    TString str1 = "/I402_";
-    //    TString str2 = "/ITS4_1";
-    TString symname, symname1;
+    modnum=0;
 
     for(Int_t c1 = 1; c1<=22; c1++){
-      symname = str0;
-      symname += c1;
-      symname += str1;
+      strEntryName1 = strSDD;
+      strEntryName1 += 3;
+      strEntryName1 +=strLadder;
+      strEntryName1 += (c1-1);
       for(Int_t c2 = 1; c2<=8; c2++){
-       symname1 = symname;
-       symname1 += c2;
-       //      symname1 += str2;
-       fgVolPath[kSDD2-kFirstLayer][modnum] = symname1.Data();
+       symname = strEntryName1;
+       symname += strSensor;
+       symname += (c2-1);
+       fgVolPath[kSDD2-kFirstLayer][modnum] = symname.Data();
        modnum++;
       }
     }
@@ -780,67 +807,64 @@ void AliAlignObj::InitSymNames()
 
   /*********************       SSD layer1  ***********************/
   {
-    Int_t modnum=0;
-    TString str0 = "ALIC_1/ITSV_1/ITSD_1/IT56_1/I565_";
-    TString str1 = "/I562_";
-    //    TString str2 = "/ITS5_1";
-    TString symname, symname1;
+    modnum=0;
 
     for(Int_t c1 = 1; c1<=34; c1++){
-      symname = str0;
-      symname += c1;
-      symname += str1;
+      strEntryName1 = strSSD;
+      strEntryName1 += 4;
+      strEntryName1 +=strLadder;
+      strEntryName1 += (c1-1);
       for(Int_t c2 = 1; c2<=22; c2++){
-       symname1 = symname;
-       symname1 += c2;
-       //      symname1 += str2;
-       fgVolPath[kSSD1-kFirstLayer][modnum] = symname1.Data();
+       symname = strEntryName1;
+       symname += strSensor;
+       symname += (c2-1);
+       fgVolPath[kSSD1-kFirstLayer][modnum] = symname.Data();
        modnum++;
       }
     }
   }
 
-  /*********************       SSD layer1  ***********************/
+  /*********************       SSD layer2  ***********************/
   {
-    Int_t modnum=0;
-    TString str0 = "ALIC_1/ITSV_1/ITSD_1/IT56_1/I569_";
-    TString str1 = "/I566_";
-    //    TString str2 = "/ITS6_1";
-    TString symname, symname1;
+    modnum=0;
 
     for(Int_t c1 = 1; c1<=38; c1++){
-      symname = str0;
-      symname += c1;
-      symname += str1;
+      strEntryName1 = strSSD;
+      strEntryName1 += 5;
+      strEntryName1 +=strLadder;
+      strEntryName1 += (c1-1);
       for(Int_t c2 = 1; c2<=25; c2++){
-       symname1 = symname;
-       symname1 += c2;
-       //      symname1 += str2;
-       fgVolPath[kSSD2-kFirstLayer][modnum] = symname1.Data();
+       symname = strEntryName1;
+       symname += strSensor;
+       symname += (c2-1);
+       fgVolPath[kSSD2-kFirstLayer][modnum] = symname.Data();
        modnum++;
       }
     }
   }
 
+
+  /***************    TPC inner and outer layers    ****************/
+  TString sAsector="TPC/EndcapA/Sector";
+  TString sCsector="TPC/EndcapC/Sector";
+  TString sInner="/InnerChamber";
+  TString sOuter="/OuterChamber";
+  
   /***************    TPC inner chambers' layer    ****************/
   {
-    Int_t modnum = 0;
-    TString str1 = "ALIC_1/TPC_M_1/TPC_Drift_1/TPC_ENDCAP_1/TPC_SECT_";
-    TString str2 = "ALIC_1/TPC_M_1/TPC_Drift_1/TPC_ENDCAP_2/TPC_SECT_";
-    TString strIn = "/TPC_IROC_1";
-    TString symname;
+    modnum = 0;
     
     for(Int_t cnt=1; cnt<=18; cnt++){
-      symname = str1;
+      symname = sAsector;
       symname += cnt;
-      symname += strIn;
+      symname += sInner;
       fgVolPath[kTPC1-kFirstLayer][modnum] = symname.Data();
       modnum++;
     }
     for(Int_t cnt=1; cnt<=18; cnt++){
-      symname = str2;
+      symname = sCsector;
       symname += cnt;
-      symname += strIn;
+      symname += sInner;
       fgVolPath[kTPC1-kFirstLayer][modnum] = symname.Data();
       modnum++;
     }
@@ -848,23 +872,19 @@ void AliAlignObj::InitSymNames()
 
   /***************    TPC outer chambers' layer    ****************/
   {
-    Int_t modnum = 0;
-    TString str1 = "ALIC_1/TPC_M_1/TPC_Drift_1/TPC_ENDCAP_1/TPC_SECT_";
-    TString str2 = "ALIC_1/TPC_M_1/TPC_Drift_1/TPC_ENDCAP_2/TPC_SECT_";
-    TString strOut = "/TPC_OROC_1";
-    TString symname;
+    modnum = 0;
     
     for(Int_t cnt=1; cnt<=18; cnt++){
-      symname = str1;
+      symname = sAsector;
       symname += cnt;
-      symname += strOut;
+      symname += sOuter;
       fgVolPath[kTPC2-kFirstLayer][modnum] = symname.Data();
       modnum++;
     }
     for(Int_t cnt=1; cnt<=18; cnt++){
-      symname = str2;
+      symname = sCsector;
       symname += cnt;
-      symname += strOut;
+      symname += sOuter;
       fgVolPath[kTPC2-kFirstLayer][modnum] = symname.Data();
       modnum++;
     }
@@ -872,36 +892,26 @@ void AliAlignObj::InitSymNames()
 
   /*********************       TOF layer   ***********************/
   {
+    modnum=0;
+    
     Int_t nstrA=15;
     Int_t nstrB=19;
     Int_t nstrC=19;
-    Int_t nsec=18;
-    Int_t nStripSec=nstrA+2*nstrB+2*nstrC;
-    Int_t nStrip=nStripSec*nsec;
-
-    for (Int_t modnum=0; modnum < nStrip; modnum++) {
-
-      Int_t sector = modnum/nStripSec;
-      Char_t  string1[100];
-      Char_t  string2[100];
-      Int_t icopy=-1;
-      icopy=sector; 
-
-      // Old 6h convention
-      // if(sector<13){
-      //       icopy=sector+5;}  
-      // else{ icopy=sector-13;}
-
-      sprintf(string1,"/ALIC_1/B077_1/BSEGMO%i_1/BTOF%i_1/FTOA_0/FLTA_0",sector,sector);
-      
-      Int_t strInSec=modnum%nStripSec;
-      icopy= strInSec;
-      icopy++;
-      sprintf(string2,"FSTR_%i",icopy);      
-      Char_t  path[100];
-      sprintf(path,"%s/%s",string1,string2); 
-      //      printf("%d  %s\n",modnum,path);
-      fgVolPath[kTOF-kFirstLayer][modnum] = path;
+    Int_t nSectors=18;
+    Int_t nStrips=nstrA+2*nstrB+2*nstrC;
+    
+    TString snSM  = "TOF/sm";
+    TString snSTRIP = "/strip";
+    
+    for (Int_t isect = 0; isect < nSectors; isect++) {
+      for (Int_t istr = 1; istr <= nStrips; istr++) {  
+       symname  = snSM;
+       symname += Form("%02d",isect);
+       symname += snSTRIP;
+       symname += Form("%02d",istr);
+       fgVolPath[kTOF-kFirstLayer][modnum] = symname.Data();   
+       modnum++;
+      }
     }
   } 
 
@@ -910,53 +920,37 @@ void AliAlignObj::InitSymNames()
     TString str = "ALIC_1/RICH_";
     TString symname;
 
-    for (Int_t modnum=0; modnum < 7; modnum++) {
+    for (modnum=0; modnum < 7; modnum++) {
       symname = str;
       symname += (modnum+1);
       fgVolPath[kRICH-kFirstLayer][modnum] = symname.Data();
     }
   }
 
-  /*********************      TRD layers 0-6   *******************/
+  /*********************      TRD layers 1-6   *******************/
+  //!! 6 layers with index increasing in outwards direction
   {
-    TString strSM[18]={"ALIC_1/B077_1/BSEGMO0_1/BTRD0_1/UTR1_1/UTS1_1/UTI1_1/UT",
-                      "ALIC_1/B077_1/BSEGMO1_1/BTRD1_1/UTR1_1/UTS1_1/UTI1_1/UT",
-                      "ALIC_1/B077_1/BSEGMO2_1/BTRD2_1/UTR1_1/UTS1_1/UTI1_1/UT",
-                      "ALIC_1/B077_1/BSEGMO3_1/BTRD3_1/UTR1_1/UTS1_1/UTI1_1/UT",
-                      "ALIC_1/B077_1/BSEGMO4_1/BTRD4_1/UTR1_1/UTS1_1/UTI1_1/UT",
-                      "ALIC_1/B077_1/BSEGMO5_1/BTRD5_1/UTR1_1/UTS1_1/UTI1_1/UT",
-                      "ALIC_1/B077_1/BSEGMO6_1/BTRD6_1/UTR1_1/UTS1_1/UTI1_1/UT",
-                      "ALIC_1/B077_1/BSEGMO7_1/BTRD7_1/UTR1_1/UTS1_1/UTI1_1/UT",
-                      "ALIC_1/B077_1/BSEGMO8_1/BTRD8_1/UTR1_1/UTS1_1/UTI1_1/UT",
-                      "ALIC_1/B077_1/BSEGMO9_1/BTRD9_1/UTR1_1/UTS1_1/UTI1_1/UT",
-                      "ALIC_1/B077_1/BSEGMO10_1/BTRD10_1/UTR1_1/UTS1_1/UTI1_1/UT",
-                      "ALIC_1/B077_1/BSEGMO11_1/BTRD11_1/UTR1_1/UTS1_1/UTI1_1/UT",
-                      "ALIC_1/B077_1/BSEGMO12_1/BTRD12_1/UTR1_1/UTS1_1/UTI1_1/UT",
-                      "ALIC_1/B077_1/BSEGMO13_1/BTRD13_1/UTR1_1/UTS1_1/UTI1_1/UT",
-                      "ALIC_1/B077_1/BSEGMO14_1/BTRD14_1/UTR1_1/UTS1_1/UTI1_1/UT",
-                      "ALIC_1/B077_1/BSEGMO15_1/BTRD15_1/UTR1_1/UTS1_1/UTI1_1/UT",
-                      "ALIC_1/B077_1/BSEGMO16_1/BTRD16_1/UTR1_1/UTS1_1/UTI1_1/UT",
-                      "ALIC_1/B077_1/BSEGMO17_1/BTRD17_1/UTR1_1/UTS1_1/UTI1_1/UT"};
-    TString strPost = "_1";
-    TString zeroStr = "0";
-    TString symname;
-
     Int_t arTRDlayId[6] = {kTRD1, kTRD2, kTRD3, kTRD4, kTRD5, kTRD6};
 
+    TString snStr  = "TRD/sm";
+    TString snApp1 = "/st";
+    TString snApp2 = "/pl";
+    
     for(Int_t layer=0; layer<6; layer++){
-      Int_t modnum=0;
-      for(Int_t sm = 0; sm < 18; sm++){
-       for(Int_t stacknum = 0; stacknum < 5; stacknum++){
-         Int_t chnum = layer + stacknum*6;
-         symname = strSM[sm];
-         if(chnum<10) symname += zeroStr;
-         symname += chnum;
-         symname += strPost;
+      modnum=0;
+      for (Int_t isect = 0; isect < 18; isect++) {
+       for (Int_t icham = 0; icham < 5; icham++) {
+         symname  = snStr;
+         symname += Form("%02d",isect);
+         symname += snApp1;
+         symname += icham;
+         symname += snApp2;
+         symname += layer;
          fgVolPath[arTRDlayId[layer]-kFirstLayer][modnum] = symname.Data();
          modnum++;
        }
       }
     }
   }
-
 }
+
index fe2fa98..85d71fd 100644 (file)
@@ -23,7 +23,7 @@ void MakeTOFFullMisAlignment(){
 
   for(i=0; i<18; i++) {
     Char_t  path[100];
-    sprintf(path,"/ALIC_1/B077_1/BSEGMO%i_1/BTOF%i_1/FTOA_0",i,i);
+    sprintf(path,"TOF/sm%02d",i);
     dx = rnd->Gaus(0.,sigmatr);
     dy = 0;
     dz = rnd->Gaus(0.,sigmatr);
index f0b423b..db6c149 100644 (file)
@@ -18,16 +18,16 @@ void MakeTPCFullMisAlignment(){
   for (Int_t iLayer = AliAlignObj::kTPC1; iLayer <= AliAlignObj::kTPC2; iLayer++) {
     for (Int_t iModule = 0; iModule < AliAlignObj::LayerSize(iLayer); iModule++) {
 
-      Float_t dx = (rnd->Uniform()-0.5)*sigmatr;
-      Float_t dy = (rnd->Uniform()-0.5)*sigmatr;
-      Float_t dz = (rnd->Uniform()-0.5)*sigmatr;
-      Float_t dpsi = (rnd->Uniform()-0.5)*sigmarot;
-      Float_t dtheta = (rnd->Uniform()-0.5)*sigmarot;
-      Float_t dphi = (rnd->Uniform()-0.5)*sigmarot;
+      Double_t dx = (rnd->Uniform()-0.5)*sigmatr;
+      Double_t dy = (rnd->Uniform()-0.5)*sigmatr;
+      Double_t dz = (rnd->Uniform()-0.5)*sigmatr;
+      Double_t dpsi = (rnd->Uniform()-0.5)*sigmarot;
+      Double_t dtheta = (rnd->Uniform()-0.5)*sigmarot;
+      Double_t dphi = (rnd->Uniform()-0.5)*sigmarot;
 
       UShort_t volid = AliAlignObj::LayerToVolUID(iLayer,iModule);
-      const char *path = AliAlignObj::SymName(volid);
-      new(alobj[j]) AliAlignObjAngles(path, volid, dx, dy, dz, dpsi, dtheta, dphi, kFALSE);
+      const char *symname = AliAlignObj::SymName(volid);
+      new(alobj[j]) AliAlignObjAngles(symname, volid, dx, dy, dz, dpsi, dtheta, dphi, kFALSE);
       j++;
     }
   }
index 8d170c1..658fb17 100644 (file)
@@ -28,8 +28,8 @@ void MakeTPCResMisAlignment(){
       dphi = (rnd->Uniform()-0.5)*sigmarot;
 
       UShort_t volid = AliAlignObj::LayerToVolUID(iLayer,iModule);
-      const char *path = AliAlignObj::SymName(volid);
-      new(alobj[j]) AliAlignObjAngles(path, volid, dx, dy, dz, dpsi, dtheta, dphi, kFALSE);
+      const char *symname = AliAlignObj::SymName(volid);
+      new(alobj[j]) AliAlignObjAngles(symname, volid, dx, dy, dz, dpsi, dtheta, dphi, kFALSE);
       j++;
     }
   }
index 6a05178..5fc6403 100644 (file)
@@ -24,14 +24,6 @@ void MakeTRDFullMisAlignment(){
   UShort_t volid;
   const char *path;
 
-  //generate the paths for the 18 supermodules
-  Double_t r=325.;
-  Double_t phi;
-  Int_t isec;
-  TString ts0;
-  TString* sm_path = new TString[18];
-  TRegexp regexp(".*BTRD[0-9][^/]*");
-
   // create the chambers' alignment objects
   for (Int_t iLayer = AliAlignObj::kTRD1; iLayer <= AliAlignObj::kTRD6; iLayer++) {
     for (Int_t iModule = 0; iModule < AliAlignObj::LayerSize(iLayer); iModule++) {
@@ -45,8 +37,8 @@ void MakeTRDFullMisAlignment(){
       ry*=chry;
       rz*=chrz;
       volid = AliAlignObj::LayerToVolUID(iLayer,iModule);
-      path = AliAlignObj::SymName(volid);
-      new(alobj[j++]) AliAlignObjAngles(path,volid,dx,dy,dz,rx,ry,rz,kFALSE);
+      symname = AliAlignObj::SymName(volid);
+      new(alobj[j++]) AliAlignObjAngles(symname,volid,dx,dy,dz,rx,ry,rz,kFALSE);
     }
   }
 
index a4db0f9..507fe5b 100644 (file)
@@ -9,14 +9,6 @@ void MakeTRDResMisAlignment(){
 
   AliAlignObjAngles a;
 
-  // sigmas for the supermodules
-  Double_t smdx=3; // 300 microns
-  Double_t smdy=3; // 300 microns
-  Double_t smdz=3; // 300 microns
-  Double_t smrx=0.4/1000/TMath::Pi()*180; // 0.4 mrad
-  Double_t smry=2.0/1000/TMath::Pi()*180; // 2 mrad
-  Double_t smrz=0.4/1000/TMath::Pi()*180; // 0.4 mrad
-
   // sigmas for the chambers
   Double_t chdx=0.02; // 200 microns
   Double_t chdy=0.03; // 300 microns
@@ -32,36 +24,6 @@ void MakeTRDResMisAlignment(){
   UShort_t volid;
   const char *path;
 
-  //generate the paths for the 18 supermodules
-  Double_t r=325.;
-  Double_t phi;
-  Int_t isec;
-  TString ts0;
-  TString* sm_path = new TString[18];
-  TRegexp regexp(".*BTRD[0-9][^/]*");
-
-  for (isec=0; isec<18; isec++) {
-    phi=(isec+0.5)/18*2*TMath::Pi();
-    gGeoManager->FindNode(r*cos(phi),r*sin(phi),0.);
-    ts0 = gGeoManager->GetPath();
-    sm_path[isec] = ts0(regexp);
-    sm_path[isec] +='\0';
-  }
-
-  // create the supermodules' alignment objects
-  for (int i; i<18; i++) {
-    ran.Rannor(dx,rx);
-    ran.Rannor(dy,ry);
-    ran.Rannor(dz,rz);
-    dx*=smdx;
-    dy*=smdy;
-    dz*=smdz;
-    rx*=smrx;
-    ry*=smry;
-    rz*=smrz;
-    new(alobj[j++]) AliAlignObjAngles(sm_path[i].Data(),0,dx,dy,dz,rx,ry,rz,kFALSE);
-  }
-
   // create the chambers' alignment objects
   for (Int_t iLayer = AliAlignObj::kTRD1; iLayer <= AliAlignObj::kTRD6; iLayer++) {
     for (Int_t iModule = 0; iModule < AliAlignObj::LayerSize(iLayer); iModule++) {
@@ -75,8 +37,8 @@ void MakeTRDResMisAlignment(){
       ry*=chry;
       rz*=chrz;
       volid = AliAlignObj::LayerToVolUID(iLayer,iModule);
-      path = AliAlignObj::SymName(volid);
-      new(alobj[j++]) AliAlignObjAngles(path,volid,dx,dy,dz,rx,ry,rz,kFALSE);
+      symname = AliAlignObj::SymName(volid);
+      new(alobj[j++]) AliAlignObjAngles(symname,volid,dx,dy,dz,rx,ry,rz,kFALSE);
     }
   }