X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=PHOS%2FMakePHOSResMisAlignment.C;h=cdf40c122506693e6a3630938135d2f59c5bfae3;hb=1ee5d9a3f69a886521b4d35c2ed1b34c6f4c72d8;hp=16ecbdf24e8599e53ccac50205d6b9391ea16702;hpb=3c759505f362cb280f68c4cb5c9566c65977b632;p=u%2Fmrichter%2FAliRoot.git diff --git a/PHOS/MakePHOSResMisAlignment.C b/PHOS/MakePHOSResMisAlignment.C index 16ecbdf24e8..cdf40c12250 100644 --- a/PHOS/MakePHOSResMisAlignment.C +++ b/PHOS/MakePHOSResMisAlignment.C @@ -8,79 +8,95 @@ void MakePHOSResMisAlignment(){ } AliPHOSEMCAGeometry *emca = phosGeom->GetEMCAGeometry(); - TClonesArray *array = new TClonesArray("AliAlignObjAngles", 16 + phosGeom->GetNModules() * + TClonesArray *array = new TClonesArray("AliAlignObjParams", 16 + phosGeom->GetNModules() * emca->GetNStripX() * emca->GetNStripZ()); TClonesArray &alobj = *array; Double_t dpsi=0., dtheta=0., dphi=0.; Double_t displacement = 0.2; Int_t iIndex=0; // let all modules have index=0 in a layer with no LUT - AliAlignObj::ELayerID iLayer = AliAlignObj::kInvalidLayer; - UShort_t volid = AliAlignObj::LayerToVolUID(iLayer,iIndex); + AliGeomManager::ELayerID iLayer = AliGeomManager::kInvalidLayer; + UShort_t volid = AliGeomManager::LayerToVolUID(iLayer,iIndex); Int_t i=0 ; // Alignment for 5 PHOS modules - new(alobj[i++]) AliAlignObjAngles("PHOS/Module1", + new(alobj[i++]) AliAlignObjParams("PHOS/Module1", volid, -0.20, -0.1, +0.0, dpsi, dtheta, 0.2, kTRUE); - new(alobj[i++]) AliAlignObjAngles("PHOS/Module2", + new(alobj[i++]) AliAlignObjParams("PHOS/Module2", volid, -0.10, +0.0, -0.2, dpsi, dtheta, 0.2, kTRUE); - new(alobj[i++]) AliAlignObjAngles("PHOS/Module3", + new(alobj[i++]) AliAlignObjParams("PHOS/Module3", volid, 0.05, -0.1, 0.2, dpsi, dtheta, 0.0, kTRUE); - new(alobj[i++]) AliAlignObjAngles("PHOS/Module4", + new(alobj[i++]) AliAlignObjParams("PHOS/Module4", volid, +0.10, -0.0, -0.1, dpsi, dtheta, 0.1, kTRUE); - new(alobj[i++]) AliAlignObjAngles("PHOS/Module5", + new(alobj[i++]) AliAlignObjParams("PHOS/Module5", volid, +0.20, -0.1, 0.1, dpsi, dtheta, 0.2, kTRUE); Double_t dx=0., dy=0., dz=0. ; // Alignment of CPV modules - new(alobj[i++]) AliAlignObjAngles("PHOS/Module1/CPV", + new(alobj[i++]) AliAlignObjParams("PHOS/Module1/CPV", volid, dx, dy, dz, dpsi, dtheta, dphi, kTRUE); - new(alobj[i++]) AliAlignObjAngles("PHOS/Module2/CPV", + new(alobj[i++]) AliAlignObjParams("PHOS/Module2/CPV", volid, dx, dy, dz, dpsi, dtheta, dphi, kTRUE); - new(alobj[i++]) AliAlignObjAngles("PHOS/Module3/CPV", + new(alobj[i++]) AliAlignObjParams("PHOS/Module3/CPV", volid, dx, dy, dz, dpsi, dtheta, dphi, kTRUE); - new(alobj[i++]) AliAlignObjAngles("PHOS/Module4/CPV", + new(alobj[i++]) AliAlignObjParams("PHOS/Module4/CPV", volid, dx, dy, dz, dpsi, dtheta, dphi, kTRUE); - new(alobj[i++]) AliAlignObjAngles("PHOS/Module5/CPV", + new(alobj[i++]) AliAlignObjParams("PHOS/Module5/CPV", volid, dx, dy, dz, dpsi, dtheta, dphi, kTRUE); // Alignment for PHOS cradle - new(alobj[i++]) AliAlignObjAngles("PHOS/Cradle0", + new(alobj[i++]) AliAlignObjParams("PHOS/Cradle0", volid, 0., 0., -displacement, dpsi, dtheta, dphi, kTRUE); - new(alobj[i++]) AliAlignObjAngles("PHOS/Cradle1", + new(alobj[i++]) AliAlignObjParams("PHOS/Cradle1", volid, 0., 0., +displacement, dpsi, dtheta, dphi, kTRUE); // Alignment for cradle wheels - new(alobj[i++]) AliAlignObjAngles("PHOS/Wheel0", + new(alobj[i++]) AliAlignObjParams("PHOS/Wheel0", volid, 0., 0., -displacement, dpsi, dtheta, dphi, kTRUE); - new(alobj[i++]) AliAlignObjAngles("PHOS/Wheel1", + new(alobj[i++]) AliAlignObjParams("PHOS/Wheel1", volid, 0., 0., -displacement, dpsi, dtheta, dphi, kTRUE); - new(alobj[i++]) AliAlignObjAngles("PHOS/Wheel2", + new(alobj[i++]) AliAlignObjParams("PHOS/Wheel2", volid, 0., 0., +displacement, dpsi, dtheta, dphi, kTRUE); - new(alobj[i++]) AliAlignObjAngles("PHOS/Wheel3", + new(alobj[i++]) AliAlignObjParams("PHOS/Wheel3", volid, 0., 0., +displacement, dpsi, dtheta, dphi, kTRUE); AliPHOSSurvey geodesicData("phos_mod3_survey.txt"); - geodesicData.CreateAliAlignObjAngles(alobj); + geodesicData.CreateAliAlignObjParams(alobj); - if(!gSystem->Getenv("TOCDB")){ + const char* macroname = "MakePHOSResMisAlignment.C"; + if( TString(gSystem->Getenv("TOCDB")) != TString("kTRUE") ){ // save on file - TFile f("PHOSresidualMisalignment.root","RECREATE"); - if(!f) cerr<<"cannot open file for output\n"; + const char* filename = "PHOSresidualMisalignment.root"; + TFile f(filename,"RECREATE"); + if(!f){ + Error(macroname,"cannot open file for output\n"); + return; + } + Info(macroname,"Saving alignment objects to the file %s", filename); f.cd(); - f.WriteObject(array,"PHOSResidualObjs ","kSingleKey"); + f.WriteObject(array,"PHOSAlignObjs","kSingleKey"); f.Close(); }else{ // save in CDB storage - const char* Storage = gSystem->Getenv("STORAGE"); - AliCDBManager *CDB = AliCDBManager::Instance(); - AliCDBStorage* storage = CDB->GetStorage(Storage); + TString Storage = gSystem->Getenv("STORAGE"); + if(!Storage.BeginsWith("local://") && !Storage.BeginsWith("alien://")) { + Error(macroname,"STORAGE variable set to %s is not valid. Exiting\n",Storage.Data()); + return; + } + Info(macroname,"Saving alignment objects in CDB storage %s", + Storage.Data()); + AliCDBManager* cdb = AliCDBManager::Instance(); + AliCDBStorage* storage = cdb->GetStorage(Storage.Data()); + if(!storage){ + Error(macroname,"Unable to open storage %s\n",Storage.Data()); + return; + } AliCDBMetaData *md= new AliCDBMetaData(); md->SetResponsible("Yuri Kharlov"); md->SetComment("Alignment objects for slightly misaligned geometry (residual misalignment"); md->SetAliRootVersion(gSystem->Getenv("ARVERSION")); - AliCDBId id("PHOS/Align/Data",0,9999999); + AliCDBId id("PHOS/Align/Data",0,AliCDBRunRange::Infinity()); storage->Put(array,id, md); }