PHOS/MakePHOSZeroMisAlignment.C

   1 void MakePHOSZeroMisAlignment(){
   2   // Create TClonesArray of zero misalignment objects for PHOS
   3   //
   4   const char* macroname = "MakePHOSZeroMisAlignment.C";
   5   const AliPHOSGeometry *phosGeom = AliPHOSGeometry::GetInstance("IHEP", "IHEP");
   6   if (!phosGeom) {
   7     Error("MakePHOSFullMisAlignment", "Cannot obtain AliPHOSGeometry singleton\n");
   8     return;
   9   }
  10
  11   AliPHOSEMCAGeometry *emca = phosGeom->GetEMCAGeometry();
  12   TClonesArray *array = new TClonesArray("AliAlignObjParams", 16 + phosGeom->GetNModules() *
  13                                          emca->GetNStripX() * emca->GetNStripZ());
  14   TClonesArray &alobj = *array;
  15
  16   Double_t dx=0., dy=0., dz=0., dpsi=0., dtheta=0., dphi=0.;
  17   Int_t iIndex=0; // let all modules have index=0 in a layer with no LUT
  18   AliGeomManager::ELayerID iLayer = AliGeomManager::kInvalidLayer;
  19   UShort_t volid = AliGeomManager::LayerToVolUID(iLayer,iIndex);
  20
  21   Int_t i=0 ;
  22   // Alignment for 5 PHOS modules
  23   new(alobj[i++]) AliAlignObjParams("PHOS/Module1",
  24         volid, dx, dy, dz, dpsi, dtheta, dphi, kTRUE);
  25   new(alobj[i++]) AliAlignObjParams("PHOS/Module2",
  26         volid, dx, dy, dz, dpsi, dtheta, dphi, kTRUE);
  27   new(alobj[i++]) AliAlignObjParams("PHOS/Module3",
  28         volid, dx, dy, dz, dpsi, dtheta, dphi, kTRUE);
  29   new(alobj[i++]) AliAlignObjParams("PHOS/Module4",
  30         volid, dx, dy, dz, dpsi, dtheta, dphi, kTRUE);
  31   new(alobj[i++]) AliAlignObjParams("PHOS/Module5",
  32         volid, dx, dy, dz, dpsi, dtheta, dphi, kTRUE);
  33
  34   // Alignment of CPV modules
  35   new(alobj[i++]) AliAlignObjParams("PHOS/Module1/CPV",
  36         volid, dx, dy, dz, dpsi, dtheta, dphi, kTRUE);
  37   new(alobj[i++]) AliAlignObjParams("PHOS/Module2/CPV",
  38         volid, dx, dy, dz, dpsi, dtheta, dphi, kTRUE);
  39   new(alobj[i++]) AliAlignObjParams("PHOS/Module3/CPV",
  40         volid, dx, dy, dz, dpsi, dtheta, dphi, kTRUE);
  41   new(alobj[i++]) AliAlignObjParams("PHOS/Module4/CPV",
  42         volid, dx, dy, dz, dpsi, dtheta, dphi, kTRUE);
  43   new(alobj[i++]) AliAlignObjParams("PHOS/Module5/CPV",
  44         volid, dx, dy, dz, dpsi, dtheta, dphi, kTRUE);
  45
  46
  47   // Alignment for PHOS cradle
  48   new(alobj[i++]) AliAlignObjParams("PHOS/Cradle0",
  49         volid, dx, dy, dz, dpsi, dtheta, dphi, kTRUE);
  50   new(alobj[i++]) AliAlignObjParams("PHOS/Cradle1",
  51         volid, dx, dy, dz, dpsi, dtheta, dphi, kTRUE);
  52
  53   // Alignment for cradle wheels
  54   new(alobj[i++])  AliAlignObjParams("PHOS/Wheel0",
  55         volid, dx, dy, dz, dpsi, dtheta, dphi, kTRUE);
  56   new(alobj[i++])  AliAlignObjParams("PHOS/Wheel1",
  57         volid, dx, dy, dz, dpsi, dtheta, dphi, kTRUE);
  58   new(alobj[i++])  AliAlignObjParams("PHOS/Wheel2",
  59         volid, dx, dy, dz, dpsi, dtheta, dphi, kTRUE);
  60   new(alobj[i++]) AliAlignObjParams("PHOS/Wheel3",
  61         volid, dx, dy, dz, dpsi, dtheta, dphi, kTRUE);
  62
  63   AliPHOSSurvey geodesicData;
  64   geodesicData.CreateNullObjects(alobj, phosGeom);
  65
  66   if( TString(gSystem->Getenv("TOCDB")) != TString("kTRUE") ){
  67     // save on file
  68     const char* filename = "PHOSzeroMisalignment.root";
  69     TFile f(filename,"RECREATE");
  70     if(!f){
  71       Error(macroname,"cannot open file for output\n");
  72       return;
  73     }
  74     Info(macroname,"Saving alignment objects to the file %s", filename);
  75     f.cd();
  76     f.WriteObject(array,"PHOSAlignObjs","kSingleKey");
  77     f.Close();
  78   }else{
  79     // save in CDB storage
  80     TString Storage = gSystem->Getenv("STORAGE");
  81     if(!Storage.BeginsWith("local://") && !Storage.BeginsWith("alien://")) {
  82       Error(macroname,"STORAGE variable set to %s is not valid. Exiting\n",Storage.Data());
  83       return;
  84     }
  85     Info(macroname,"Saving alignment objects in CDB storage %s",
  86          Storage.Data());
  87     AliCDBManager* cdb = AliCDBManager::Instance();
  88     AliCDBStorage* storage = cdb->GetStorage(Storage.Data());
  89     if(!storage){
  90       Error(macroname,"Unable to open storage %s\n",Storage.Data());
  91       return;
  92     }
  93     AliCDBMetaData *md= new AliCDBMetaData();
  94     md->SetResponsible("Yuri Kharlov");
  95     md->SetComment("Zero misalignment objects");
  96     md->SetAliRootVersion(gSystem->Getenv("ARVERSION"));
  97     AliCDBId id("PHOS/Align/Data",0,AliCDBRunRange::Infinity());
  98     storage->Put(array,id, md);
  99   }
 100
 101   array->Delete();
 102
 103 }