TOF/MakeTOFZeroMisAlignment.C

   1 void MakeTOFZeroMisAlignment(){
   2   // Create TClonesArray of zero misalignment objects for TOF
   3   //
   4   TClonesArray *array = new TClonesArray("AliAlignObjParams",2000);
   5   TClonesArray &alobj = *array;
   6   const char* macroname = "MakeTOFZeroMisAlignment.C";
   7
   8   // Activate CDB storage and load geometry from CDB
   9   AliCDBManager* cdb = AliCDBManager::Instance();
  10   if(!cdb->IsDefaultStorageSet()) cdb->SetDefaultStorage("local://$ALICE_ROOT");
  11   cdb->SetRun(0);
  12
  13   AliCDBStorage* storage;
  14
  15   if( TString(gSystem->Getenv("TOCDB")) == TString("kTRUE") ){
  16     TString Storage = gSystem->Getenv("STORAGE");
  17     if(!Storage.BeginsWith("local://") && !Storage.BeginsWith("alien://")) {
  18       Error(macroname,"STORAGE variable set to %s is not valid. Exiting\n",Storage.Data());
  19       return;
  20     }
  21     storage = cdb->GetStorage(Storage.Data());
  22     if(!storage){
  23       Error(macroname,"Unable to open storage %s\n",Storage.Data());
  24       return;
  25     }
  26     AliCDBPath path("GRP","Geometry","Data");
  27     AliCDBEntry *entry = storage->Get(path.GetPath(),cdb->GetRun());
  28     if(!entry) Fatal(macroname,"Could not get the specified CDB entry!");
  29     entry->SetOwner(0);
  30     TGeoManager* geom = (TGeoManager*) entry->GetObject();
  31     AliGeomManager::SetGeometry(geom);
  32   }else{
  33     AliGeomManager::LoadGeometry(); //load geom from default CDB storage
  34   }
  35
  36
  37   AliGeomManager::ELayerID idTOF = AliGeomManager::kTOF;
  38   Int_t j=0;
  39   Int_t strId=-1;
  40   Double_t dx=0., dy=0., dz=0., dpsi=0., dtheta=0., dphi=0.;
  41
  42   Int_t nstrA=15;
  43   Int_t nstrB=19;
  44   Int_t nstrC=19;
  45   Int_t nSectors=18;
  46   Int_t nStrips=nstrA+2*nstrB+2*nstrC;
  47
  48   for (Int_t isect = 0; isect < nSectors; isect++) {
  49     for (Int_t istr = 1; istr <= nStrips; istr++) {
  50       strId++;
  51       if ((isect==13 || isect==14 || isect==15) && (istr >= 39 && istr <= 53)) continue;
  52       new(alobj[j++]) AliAlignObjParams(AliGeomManager::SymName(idTOF,strId), AliGeomManager::LayerToVolUID(idTOF,strId), dx, dy, dz, dpsi, dtheta, dphi, kTRUE);
  53     }
  54   }
  55
  56   if( TString(gSystem->Getenv("TOCDB")) != TString("kTRUE") ){
  57     // save on file
  58     const char* filename = "TOFzeroMisalignment.root";
  59     TFile f(filename,"RECREATE");
  60     if(!f){
  61       Error(macroname,"cannot open file for output\n");
  62       return;
  63     }
  64     Info(macroname,"Saving alignment objects to the file %s", filename);
  65     f.cd();
  66     f.WriteObject(array,"TOFAlignObjs","kSingleKey");
  67     f.Close();
  68   }else{
  69     // save in CDB storage
  70     AliCDBMetaData* md = new AliCDBMetaData();
  71     md->SetResponsible("Silvia Arcelli");
  72     md->SetComment("Zero misalignment for TOF");
  73     md->SetAliRootVersion("HEAD");
  74     AliCDBId id("TOF/Align/Data",0,AliCDBRunRange::Infinity());
  75     storage->Put(array,id,md);
  76   }
  77
  78   array->Delete();
  79
  80 }
  81
  82