+void AliTPCv2::AddAlignableVolumes() const
+{
+ //
+ // Create entries for alignable volumes associating the symbolic volume
+ // name with the corresponding volume path. Needs to be syncronized with
+ // eventual changes in the geometry.
+ //
+ SetInnerChambersAlignable();
+ SetOuterChambersAlignable();
+}
+
+//_____________________________________________________________________________
+void AliTPCv2::SetInnerChambersAlignable() const
+{
+ //
+ Int_t modnum = 0;
+ TString vpstr1 = "ALIC_1/TPC_M_1/TPC_Drift_1/TPC_ENDCAP_1/TPC_SECT_";
+ TString vpstr2 = "ALIC_1/TPC_M_1/TPC_Drift_1/TPC_ENDCAP_2/TPC_SECT_";
+ TString vpappend = "/TPC_IROC_1";
+ TString snstr1="TPC/EndcapA/Sector";
+ TString snstr2="TPC/EndcapC/Sector";
+ TString snappend="/InnerChamber";
+ TString volpath, symname;
+
+ for(Int_t cnt=1; cnt<=18; cnt++){
+ volpath = vpstr1;
+ volpath += cnt;
+ volpath += vpappend;
+ symname = snstr1;
+ symname += cnt;
+ symname += snappend;
+ if(!gGeoManager->SetAlignableEntry(symname.Data(),volpath.Data()))
+ AliFatal(Form("Alignable entry %s not created. Volume path %s not valid", symname.Data(),volpath.Data()));
+ //
+ TGeoPNEntry *alignableEntry = gGeoManager->GetAlignableEntry(symname.Data());
+ const char *path = alignableEntry->GetTitle();
+ if (!gGeoManager->cd(path))
+ AliFatal(Form("Volume path %s not valid!",path));
+ TGeoHMatrix* globMatrix = gGeoManager->GetCurrentMatrix();
+ TGeoHMatrix* matTtoL = fTPCParam->Tracking2LocalMatrix(globMatrix,cnt-1);
+ alignableEntry->SetMatrix(matTtoL);
+ modnum++;
+ }
+
+ for(Int_t cnt=1; cnt<=18; cnt++){
+ volpath = vpstr2;
+ volpath += cnt;
+ volpath += vpappend;
+ symname = snstr2;
+ symname += cnt;
+ symname += snappend;
+ if(!gGeoManager->SetAlignableEntry(symname.Data(),volpath.Data()))
+ AliFatal(Form("Alignable entry %s not created. Volume path %s not valid", symname.Data(),volpath.Data()));
+ TGeoPNEntry *alignableEntry = gGeoManager->GetAlignableEntry(symname.Data());
+ const char *path = alignableEntry->GetTitle();
+ if (!gGeoManager->cd(path))
+ AliFatal(Form("Volume path %s not valid!",path));
+ TGeoHMatrix* globMatrix = gGeoManager->GetCurrentMatrix();
+ TGeoHMatrix* matTtoL = fTPCParam->Tracking2LocalMatrix(globMatrix,18+cnt-1);
+ alignableEntry->SetMatrix(matTtoL);
+ modnum++;
+ }
+}
+
+//_____________________________________________________________________________
+void AliTPCv2::SetOuterChambersAlignable() const
+{
+ //
+ Int_t modnum = 0;
+ TString vpstr1 = "ALIC_1/TPC_M_1/TPC_Drift_1/TPC_ENDCAP_1/TPC_SECT_";
+ TString vpstr2 = "ALIC_1/TPC_M_1/TPC_Drift_1/TPC_ENDCAP_2/TPC_SECT_";
+ TString vpappend = "/TPC_OROC_1";
+ TString snstr1="TPC/EndcapA/Sector";
+ TString snstr2="TPC/EndcapC/Sector";
+ TString snappend="/OuterChamber";
+ TString volpath, symname;
+
+ for(Int_t cnt=1; cnt<=18; cnt++){
+ volpath = vpstr1;
+ volpath += cnt;
+ volpath += vpappend;
+ symname = snstr1;
+ symname += cnt;
+ symname += snappend;
+ if(!gGeoManager->SetAlignableEntry(symname.Data(),volpath.Data()))
+ AliFatal(Form("Alignable entry %s not created. Volume path %s not valid", symname.Data(),volpath.Data()));
+ TGeoPNEntry *alignableEntry = gGeoManager->GetAlignableEntry(symname.Data());
+ const char *path = alignableEntry->GetTitle();
+ if (!gGeoManager->cd(path))
+ AliFatal(Form("Volume path %s not valid!",path));
+ TGeoHMatrix* globMatrix = gGeoManager->GetCurrentMatrix();
+ TGeoHMatrix* matTtoL = fTPCParam->Tracking2LocalMatrix(globMatrix,36+cnt-1);
+ alignableEntry->SetMatrix(matTtoL);
+ modnum++;
+ }
+
+ for(Int_t cnt=1; cnt<=18; cnt++){
+ volpath = vpstr2;
+ volpath += cnt;
+ volpath += vpappend;
+ symname = snstr2;
+ symname += cnt;
+ symname += snappend;
+ if(!gGeoManager->SetAlignableEntry(symname.Data(),volpath.Data()))
+ AliFatal(Form("Alignable entry %s not created. Volume path %s not valid", symname.Data(),volpath.Data()));
+ TGeoPNEntry *alignableEntry = gGeoManager->GetAlignableEntry(symname.Data());
+ const char *path = alignableEntry->GetTitle();
+ if (!gGeoManager->cd(path))
+ AliFatal(Form("Volume path %s not valid!",path));
+ TGeoHMatrix* globMatrix = gGeoManager->GetCurrentMatrix();
+ TGeoHMatrix* matTtoL = fTPCParam->Tracking2LocalMatrix(globMatrix,36+18+cnt-1);
+ alignableEntry->SetMatrix(matTtoL);
+ modnum++;
+ }
+}
+
+//_____________________________________________________________________________
+void AliTPCv2::DrawDetector() const