//Init geometry
if(!fPHOSGeo){
AliOADBContainer geomContainer("phosGeo");
- geomContainer.InitFromFile("$ALICE_ROOT/OADB/PHOS/PHOSGeometry.root","PHOSRotationMatrixes");
+ if(fIsMC) //use excatly the same geometry as in simulation
+ geomContainer.InitFromFile("$ALICE_ROOT/OADB/PHOS/PHOSMCGeometry.root","PHOSMCRotationMatrixes");
+ else //Use best approaximation to real geometry
+ geomContainer.InitFromFile("$ALICE_ROOT/OADB/PHOS/PHOSGeometry.root","PHOSRotationMatrixes");
TObjArray *matrixes = (TObjArray*)geomContainer.GetObject(runNumber,"PHOSRotationMatrixes");
fPHOSGeo = AliPHOSGeometry::GetInstance("IHEP") ;
for(Int_t mod=0; mod<5; mod++) {
fPHOSGeo->SetMisalMatrix(((TGeoHMatrix*)matrixes->At(mod)),mod) ;
printf(".........Adding Matrix(%d), geo=%p\n",mod,fPHOSGeo) ;
((TGeoHMatrix*)matrixes->At(mod))->Print() ;
- }
+ }
}
//Init Bad channels map