Updated alignment macros and related changes (Raffaele)
[u/mrichter/AliRoot.git] / TRD / MakeTRDFullMisAlignment.C
1 void MakeTRDFullMisAlignment(){
2   // Create TClonesArray of full misalignment objects for TRD
3   // Expects to read objects for FRAME
4   // 
5   TClonesArray *array = new TClonesArray("AliAlignObjParams",1000);
6   const char* macroname = "MakeTRDFullMisAlignment.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   TString Storage;
15   
16   if( TString(gSystem->Getenv("TOCDB")) == TString("kTRUE") ){
17     Storage = gSystem->Getenv("STORAGE");
18     if(!Storage.BeginsWith("local://") && !Storage.BeginsWith("alien://")) {
19       Error(macroname,"STORAGE variable set to %s is not valid. Exiting\n",Storage.Data());
20       return;
21     }
22     storage = cdb->GetStorage(Storage.Data());
23     if(!storage){
24       Error(macroname,"Unable to open storage %s\n",Storage.Data());
25       return;
26     }
27     AliCDBPath path("GRP","Geometry","Data");
28     AliCDBEntry *entry = storage->Get(path.GetPath(),cdb->GetRun());
29     if(!entry) Fatal(macroname,"Could not get the specified CDB entry!");
30     entry->SetOwner(0);
31     TGeoManager* geom = (TGeoManager*) entry->GetObject();
32     AliGeomManager::SetGeometry(geom);
33   }else{
34     AliGeomManager::LoadGeometry(); //load geom from default CDB storage
35   }    
36                   
37   // load FRAME full misalignment objects (if needed, the macro
38   // for FRAME has to be ran in advance) and apply them to geometry
39   AliCDBPath fpath("GRP","Align","Data");
40   AliCDBEntry *eFrame;
41   if( TString(gSystem->Getenv("TOCDB")) == TString("kTRUE") ){
42     Info(macroname,"Loading FRAME alignment objects from CDB storage %s",
43       Storage.Data());
44     eFrame = storage->Get(fpath.GetPath(),cdb->GetRun());
45   }else{
46     eFrame = cdb->Get(fpath.GetPath());
47   }
48   if(!eFrame) Fatal(macroname,"Could not get the specified CDB entry!");
49   TClonesArray* arFrame = (TClonesArray*) eFrame->GetObject();
50   arFrame->Sort();
51   Int_t nvols = arFrame->GetEntriesFast();
52   Bool_t flag = kTRUE;
53   for(Int_t j=0; j<nvols; j++)
54   {
55     AliAlignObj* alobj = (AliAlignObj*) arFrame->UncheckedAt(j);
56     if (alobj->ApplyToGeometry() == kFALSE) flag = kFALSE;
57   }
58   if(!flag) Fatal(macroname,"Error in the application of FRAME objects");
59
60    
61   // sigmas for the supermodules
62   Double_t smdx=0.3; // 3 mm
63   Double_t smdy=0.3; // 3 mm
64   Double_t smdz=0.3; // 3 mm
65   Double_t smrx=0.4/1000/TMath::Pi()*180; // 0.4 mrad
66   Double_t smry=2.0/1000/TMath::Pi()*180; // 2 mrad
67   Double_t smrz=0.4/1000/TMath::Pi()*180; // 0.4 mrad
68
69
70   // sigmas for the chambers
71   Double_t chdx=0.1; // 1 mm
72   Double_t chdy=0.1; // 1 mm
73   Double_t chdz=0.1; // 1 mm
74   Double_t chrx=1.0/1000/TMath::Pi()*180; // 1 mrad
75   Double_t chry=1.0/1000/TMath::Pi()*180; // 1 mrad
76   Double_t chrz=0.7/1000/TMath::Pi()*180; // 0.7 mrad
77
78   Int_t sActive[18]={0,0,1,1,1,0,1,0,0,0,0,1,1,0,1,1,0,0};
79   Double_t dx,dy,dz,rx,ry,rz;
80
81   Int_t j=0;
82   TRandom *ran = new TRandom(4357);
83   UShort_t volid;
84   const char *symname;
85
86   // create the supermodules' alignment objects
87   for (int iSect; iSect<18; iSect++) {
88     TString sm_symname(Form("TRD/sm%02d",iSect));
89     ran->Rannor(dx,rx);
90     ran->Rannor(dy,ry);
91     ran->Rannor(dz,rz);
92     dx*=smdx;
93     dy*=smdy;
94     dz*=smdz;
95     rx*=smrx;
96     ry*=smry;
97     rz*=smrz;
98     if( (TString(gSystem->Getenv("PARTGEOM")) == TString("kTRUE")) && !sActive[iSect] ) continue;
99     new((*array)[j++]) AliAlignObjParams(sm_symname.Data(),0,dx,dy,dz,rx,ry,rz,kFALSE);
100   }
101   // apply supermodules' alignment objects
102   Int_t smCounter=0;
103   for(Int_t iSect=0; iSect<18; iSect++){
104     if( (TString(gSystem->Getenv("PARTGEOM")) == TString("kTRUE")) && !sActive[iSect] ) continue;
105     AliAlignObjParams* smobj =
106       (AliAlignObjParams*)array->UncheckedAt(smCounter++);
107     if(!smobj->ApplyToGeometry()){
108       Fatal(macroname,Form("application of full misalignment object for sector %d failed!",iSect));
109       return;
110     }
111   }
112
113   // create the chambers' alignment objects
114   ran = new TRandom(4357);
115   Int_t chId;
116   for (Int_t iLayer = AliGeomManager::kTRD1; iLayer <= AliGeomManager::kTRD6; iLayer++) {
117     chId=-1;
118     for (Int_t iSect = 0; iSect < 18; iSect++){
119       for (Int_t iCh = 0; iCh < 5; iCh++) {
120       ran->Rannor(dx,rx);
121       ran->Rannor(dy,ry);
122       ran->Rannor(dz,rz);
123       dx*=chdx;
124       dy*=chdy;
125       dz*=chdz;
126       rx*=chrx;
127       ry*=chry;
128       rz*=chrz;
129       chId++;
130       volid = AliGeomManager::LayerToVolUID(iLayer,chId);
131       symname = AliGeomManager::SymName(volid);
132       if( (TString(gSystem->Getenv("PARTGEOM")) == TString("kTRUE")) && !sActive[iSect] ) continue;
133       new((*array)[j++]) AliAlignObjParams(symname,volid,dx,dy,dz,rx,ry,rz,kFALSE);
134     }
135   }
136   }
137
138   if( TString(gSystem->Getenv("TOCDB")) != TString("kTRUE") ){
139     // save on file
140     const char* filename = "TRDfullMisalignment.root";
141     TFile f(filename,"RECREATE");
142     if(!f){
143       Error(macroname,"cannot open file for output\n");
144       return;
145     }
146     Info(macroname,"Saving alignment objects to the file %s", filename);
147     f.cd();
148     f.WriteObject(array,"TRDAlignObjs","kSingleKey");
149     f.Close();
150   }else{
151     // save in CDB storage
152     Info(macroname,"Saving alignment objects in CDB storage %s",
153          Storage.Data());
154     AliCDBMetaData* md = new AliCDBMetaData();
155     md->SetResponsible("Dariusz Miskowiec");
156     md->SetComment("Full misalignment for TRD");
157     md->SetAliRootVersion(gSystem->Getenv("ARVERSION"));
158     AliCDBId id("TRD/Align/Data",0,AliCDBRunRange::Infinity());
159     storage->Put(array,id,md);
160   }
161
162   array->Delete();
163 }
164