PHOS/MakePHOSResMisAlignment.C

   1 void MakePHOSResMisAlignment(){
   2   // Create TClonesArray of residual misalignment objects for PHOS
   3   //
   4   TClonesArray *array = new TClonesArray("AliAlignObjAngles",11);
   5   TClonesArray &alobj = *array;
   6
   7   AliAlignObjAngles a;
   8
   9   Double_t dpsi=0., dtheta=0., dphi=0.;
  10   Double_t displacement = 0.2;
  11
  12   Int_t iIndex=0; // let all modules have index=0 in a layer with no LUT
  13   AliAlignObj::ELayerID iLayer = AliAlignObj::kInvalidLayer;
  14   UShort_t volid = AliAlignObj::LayerToVolUID(iLayer,iIndex);
  15
  16   Int_t i=0 ;
  17   // Alignment for 5 PHOS modules
  18   new(alobj[i++]) AliAlignObjAngles("PHOS/Module1",
  19         volid, -0.20, -0.1, +0.0, dpsi, dtheta, 0.2, kTRUE);
  20   new(alobj[i++]) AliAlignObjAngles("PHOS/Module2",
  21         volid, -0.10, +0.0, -0.2, dpsi, dtheta, 0.2, kTRUE);
  22   new(alobj[i++]) AliAlignObjAngles("PHOS/Module3",
  23         volid,  0.05, -0.1,  0.2, dpsi, dtheta, 0.0, kTRUE);
  24   new(alobj[i++]) AliAlignObjAngles("PHOS/Module4",
  25         volid, +0.10, -0.0, -0.1, dpsi, dtheta, 0.1, kTRUE);
  26   new(alobj[i++]) AliAlignObjAngles("PHOS/Module5",
  27         volid, +0.20, -0.1,  0.1, dpsi, dtheta, 0.2, kTRUE);
  28
  29   Double_t dx=0., dy=0., dz=0. ;
  30   // Alignment of CPV modules
  31   new(alobj[i++]) AliAlignObjAngles("PHOS/Module1/CPV",
  32         volid, dx, dy, dz, dpsi, dtheta, dphi, kTRUE);
  33   new(alobj[i++]) AliAlignObjAngles("PHOS/Module2/CPV",
  34         volid, dx, dy, dz, dpsi, dtheta, dphi, kTRUE);
  35   new(alobj[i++]) AliAlignObjAngles("PHOS/Module3/CPV",
  36         volid, dx, dy, dz, dpsi, dtheta, dphi, kTRUE);
  37   new(alobj[i++]) AliAlignObjAngles("PHOS/Module4/CPV",
  38         volid, dx, dy, dz, dpsi, dtheta, dphi, kTRUE);
  39   new(alobj[i++]) AliAlignObjAngles("PHOS/Module5/CPV",
  40         volid, dx, dy, dz, dpsi, dtheta, dphi, kTRUE);
  41
  42
  43   // Alignment for PHOS cradle
  44   new(alobj[i++]) AliAlignObjAngles("PHOS/Cradle0",
  45         volid, 0., 0., -displacement, dpsi, dtheta, dphi, kTRUE);
  46   new(alobj[i++]) AliAlignObjAngles("PHOS/Cradle1",
  47         volid, 0., 0., +displacement, dpsi, dtheta, dphi, kTRUE);
  48
  49   // Alignment for cradle wheels
  50   new(alobj[i++])  AliAlignObjAngles("PHOS/Wheel0",
  51         volid, 0., 0., -displacement, dpsi, dtheta, dphi, kTRUE);
  52   new(alobj[i++])  AliAlignObjAngles("PHOS/Wheel1",
  53         volid, 0., 0., -displacement, dpsi, dtheta, dphi, kTRUE);
  54   new(alobj[i++])  AliAlignObjAngles("PHOS/Wheel2",
  55         volid, 0., 0., +displacement, dpsi, dtheta, dphi, kTRUE);
  56   new(alobj[i++]) AliAlignObjAngles("PHOS/Wheel3",
  57         volid, 0., 0., +displacement, dpsi, dtheta, dphi, kTRUE);
  58
  59
  60   if(!gSystem->Getenv("$TOCDB")){
  61     // save on file
  62     TFile f("PHOSresidualMisalignment.root","RECREATE");
  63     if(!f) cerr<<"cannot open file for output\n";
  64     f.cd();
  65     f.WriteObject(array,"PHOSResidualObjs ","kSingleKey");
  66     f.Close();
  67   }else{
  68     // save in CDB storage
  69     const char* Storage = gSystem->Getenv("$STORAGE");
  70     AliCDBManager *CDB = AliCDBManager::Instance();
  71     AliCDBStorage* storage = CDB->GetStorage(Storage);
  72     AliCDBMetaData *md= new AliCDBMetaData();
  73     md->SetResponsible("Yuri Kharlov");
  74     md->SetComment("Alignment objects for slightly misaligned geometry (residual misalignment");
  75     md->SetAliRootVersion(gSystem->Getenv("$ARVERSION"));
  76     AliCDBId id("PHOS/Align/Data",0,9999999);
  77     storage->Put(array,id, md);
  78   }
  79
  80   array->Delete();
  81
  82 }