TPC/MakeTPCResMisAlignment.C

   1 void MakeTPCResMisAlignment(){
   2   // Create TClonesArray of residual misalignment objects for TPC
   3   //
   4   if(!AliGeomManager::GetGeometry()){
   5     if(!(AliCDBManager::Instance())->IsDefaultStorageSet())
   6       AliCDBManager::Instance()->SetDefaultStorage("local://$ALICE_ROOT");
   7       AliCDBManager::Instance()->SetRun(0);
   8     AliGeomManager::LoadGeometry();
   9   }
  10   // needed for the constructors with local coordinates not to fail
  11
  12   TClonesArray *array = new TClonesArray("AliAlignObjAngles",100);
  13   TClonesArray &alobj = *array;
  14
  15   TRandom *rnd   = new TRandom(4357);
  16   AliAlignObjAngles o;
  17   Double_t dx, dy, dz, dpsi, dtheta, dphi;
  18   Int_t j = 0;
  19
  20   // RS = local
  21   // sigma translation = 0.1 mm
  22   // sigma rotation = 0.1 mrad
  23   Float_t sigmatr=0.01;
  24   Float_t sigmarot = 0.006;
  25   for (Int_t iLayer = AliGeomManager::kTPC1; iLayer <= AliGeomManager::kTPC2; iLayer++) {
  26     for (Int_t iModule = 0; iModule < AliGeomManager::LayerSize(iLayer); iModule++) {
  27
  28       dx = rnd->Gaus(0,sigmatr);
  29       dy = rnd->Gaus(0,sigmatr);
  30       dz = rnd->Gaus(0,sigmatr);
  31       dpsi = rnd->Gaus(0,sigmarot);
  32       dtheta = rnd->Gaus(0,sigmarot);
  33       dphi = rnd->Gaus(0,sigmarot);
  34
  35       UShort_t volid = AliGeomManager::LayerToVolUID(iLayer,iModule);
  36       const char *symname = AliGeomManager::SymName(volid);
  37       new(alobj[j]) AliAlignObjAngles(symname, volid, dx, dy, dz, dpsi, dtheta, dphi, kFALSE);
  38       j++;
  39     }
  40   }
  41
  42
  43   const char* macroname = "MakeTPCResMisAlignment.C";
  44   if( gSystem->Getenv("TOCDB") != TString("kTRUE") ){
  45     // save on file
  46     const char* filename = "TPCresidualMisalignment.root";
  47     TFile f(filename,"RECREATE");
  48     if(!f){
  49       Error(macroname,"cannot open file for output\n");
  50       return;
  51     }
  52     Info(macroname,"Saving alignment objects to the file %s", filename);
  53     f.cd();
  54     f.WriteObject(array,"TPCAlignObjs","kSingleKey");
  55     f.Close();
  56   }else{
  57     // save in CDB storage
  58     TString Storage = gSystem->Getenv("STORAGE");
  59     if(!Storage.BeginsWith("local://") && !Storage.BeginsWith("alien://")) {
  60       Error(macroname,"STORAGE variable set to %s is not valid. Exiting\n",Storage.Data());
  61       return;
  62     }
  63     Info(macroname,"Saving alignment objects in CDB storage %s",
  64          Storage.Data());
  65     AliCDBManager* cdb = AliCDBManager::Instance();
  66     AliCDBStorage* storage = cdb->GetStorage(Storage.Data());
  67     if(!storage){
  68       Error(macroname,"Unable to open storage %s\n",Storage.Data());
  69       return;
  70     }
  71     AliCDBMetaData* md = new AliCDBMetaData();
  72     md->SetResponsible("Marian Ivanov");
  73     md->SetComment("Residual misalignment for TPC, sigmatr=0.01 and sigmarot=0.6 in the local RS");
  74     md->SetAliRootVersion(gSystem->Getenv("ARVERSION"));
  75     AliCDBId id("TPC/Align/Data",0,AliCDBRunRange::Infinity());
  76     storage->Put(array,id,md);
  77   }
  78
  79   array->Delete();
  80
  81 }
  82