]>
Commit | Line | Data |
---|---|---|
1895a097 | 1 | void MakeTOFFullMisAlignment(){ |
2 | // Create TClonesArray of full misalignment objects for TOF | |
6fce62af | 3 | // Expects to read objects for FRAME |
1895a097 | 4 | // |
90dbf5fb | 5 | TClonesArray *array = new TClonesArray("AliAlignObjParams",2000); |
6fce62af | 6 | const char* macroname = "MakeTOFFullMisAlignment.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; | |
a24be56b | 14 | TString Storage; |
6fce62af | 15 | |
a24be56b | 16 | if( TString(gSystem->Getenv("TOCDB")) == TString("kTRUE") ){ |
17 | Storage = gSystem->Getenv("STORAGE"); | |
6fce62af | 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 run in advance) and apply them to geometry | |
6fce62af | 39 | AliCDBPath fpath("GRP","Align","Data"); |
a24be56b | 40 | if( TString(gSystem->Getenv("TOCDB")) == TString("kTRUE") ){ |
41 | Info(macroname,"Loading FRAME alignment objects from CDB storage %s", | |
42 | Storage.Data()); | |
6fce62af | 43 | AliCDBEntry *eFrame = storage->Get(fpath.GetPath(),cdb->GetRun()); |
6fce62af | 44 | }else{ |
45 | AliCDBEntry *eFrame = cdb->Get(fpath.GetPath()); | |
dfe9c69d | 46 | } |
a24be56b | 47 | if(!eFrame) Fatal(macroname,"Could not get the specified CDB entry!"); |
48 | TClonesArray* arFrame = (TClonesArray*) eFrame->GetObject(); | |
49 | arFrame->Sort(); | |
50 | Int_t nvols = arFrame->GetEntriesFast(); | |
51 | Bool_t flag = kTRUE; | |
52 | for(Int_t j=0; j<nvols; j++) | |
53 | { | |
54 | AliAlignObj* alobj = (AliAlignObj*) arFrame->UncheckedAt(j); | |
55 | if (alobj->ApplyToGeometry() == kFALSE) flag = kFALSE; | |
56 | } | |
57 | if(!flag) Fatal(macroname,"Error in the application of FRAME alignment objects"); | |
1895a097 | 58 | |
6fce62af | 59 | //Produce objects for TOF supermodules |
1895a097 | 60 | Int_t iIndex=0; //let all modules have index=0 in a layer with no LUT |
ae079791 | 61 | AliGeomManager::ELayerID iLayer = AliGeomManager::kInvalidLayer; |
62 | UShort_t dvoluid = AliGeomManager::LayerToVolUID(iLayer,iIndex); //dummy vol id | |
d3a9c52a | 63 | |
d3a9c52a | 64 | Int_t nSMTOF = 18; |
a24be56b | 65 | Int_t sActive[18]={0,1,1,0,0,0,1,1,0,1,1,1,1,0,0,1,1,1}; |
6fce62af | 66 | Int_t j=0; |
a24be56b | 67 | Double_t smdx, smdy, smdz=0., dpsi=0., dtheta, dphi=0.; |
1895a097 | 68 | TRandom *rnd = new TRandom(2345); |
69 | Double_t sigmatr = 0.4; // max shift in cm w.r.t. local ideal RS | |
70 | Double_t sigmarot = 0.06; // max rot in deg w.r.t. local ideal RS (~ 1 mrad) | |
d3a9c52a | 71 | |
85ce4b22 | 72 | for(Int_t isect=0; isect<nSMTOF; isect++) { |
73 | TString symname(Form("TOF/sm%02d",isect)); | |
6fce62af | 74 | smdx = rnd->Gaus(0.,sigmatr); |
75 | smdy = rnd->Gaus(0.,sigmatr); | |
1895a097 | 76 | dtheta = rnd->Gaus(0.,sigmarot); |
85ce4b22 | 77 | if( (TString(gSystem->Getenv("PARTGEOM")) == TString("kTRUE")) && !sActive[isect] ) continue; |
6fce62af | 78 | new((*array)[j++]) AliAlignObjParams(symname.Data(), dvoluid, smdx, smdy, smdz, dpsi, dtheta, dphi, kFALSE); |
d3a9c52a | 79 | } |
6fce62af | 80 | // Apply objects for TOF supermodules |
a24be56b | 81 | Int_t smCounter=0; |
85ce4b22 | 82 | for(Int_t isect=0; isect<nSMTOF; isect++){ |
83 | if( (TString(gSystem->Getenv("PARTGEOM")) == TString("kTRUE")) && !sActive[isect] ) continue; | |
a24be56b | 84 | AliAlignObjParams* smobj = (AliAlignObjParams*)array->UncheckedAt(smCounter++); |
85ce4b22 | 85 | Info(macroname,Form("Applying object for sector %d ",isect)); |
d3a9c52a | 86 | if(!smobj->ApplyToGeometry()){ |
85ce4b22 | 87 | Fatal(macroname,Form("application of full misalignment object for sector %d failed!",isect)); |
d3a9c52a | 88 | return; |
89 | } | |
90 | } | |
d3a9c52a | 91 | |
6fce62af | 92 | //Produce objects for TOF strips (same sigmas as for residual misalignment) |
ae079791 | 93 | AliGeomManager::ELayerID idTOF = AliGeomManager::kTOF; |
a24be56b | 94 | Int_t strId=-1; |
d3a9c52a | 95 | |
a24be56b | 96 | Double_t sdx=0., sdy=0., sdz=0., sdpsi=0., sdtheta=0., sdphi=0.; |
d3a9c52a | 97 | TRandom *rnds = new TRandom(4357); |
98 | sigmatr = 0.1; // max shift in cm w.r.t. local ideal RS | |
99 | ||
a24be56b | 100 | Int_t nstrA=15; |
101 | Int_t nstrB=19; | |
102 | Int_t nstrC=19; | |
103 | Int_t nSectors=18; | |
104 | Int_t nStrips=nstrA+2*nstrB+2*nstrC; | |
105 | ||
85ce4b22 | 106 | for (Int_t isect = 0; isect < nSectors; isect++) { |
a24be56b | 107 | for (Int_t istr = 1; istr <= nStrips; istr++) { |
d3a9c52a | 108 | sdy = rnds->Gaus(0.,sigmatr); |
109 | sdz = rnds->Gaus(0.,sigmatr); | |
85ce4b22 | 110 | strId++; |
111 | if ((isect==13 || isect==14 || isect==15) && (istr >= 39 && istr <= 53)) continue; | |
112 | if( (TString(gSystem->Getenv("PARTGEOM")) == TString("kTRUE")) && !sActive[isect] ) continue; | |
113 | new((*array)[j++]) AliAlignObjParams(AliGeomManager::SymName(idTOF,strId),AliGeomManager::LayerToVolUID(idTOF,strId), sdx, sdy, sdz, sdpsi, sdtheta, sdphi, kFALSE); | |
a24be56b | 114 | } |
1895a097 | 115 | } |
116 | ||
a24be56b | 117 | if( TString(gSystem->Getenv("TOCDB")) != TString("kTRUE") ){ |
6fce62af | 118 | // save in file |
dfe9c69d | 119 | const char* filename = "TOFfullMisalignment.root"; |
120 | TFile f(filename,"RECREATE"); | |
121 | if(!f){ | |
122 | Error(macroname,"cannot open file for output\n"); | |
123 | return; | |
124 | } | |
6fce62af | 125 | Info(macroname,"Saving alignment objects in file %s", filename); |
1895a097 | 126 | f.cd(); |
127 | f.WriteObject(array,"TOFAlignObjs","kSingleKey"); | |
128 | f.Close(); | |
129 | }else{ | |
130 | // save in CDB storage | |
dfe9c69d | 131 | Info(macroname,"Saving alignment objects in CDB storage %s", |
132 | Storage.Data()); | |
1895a097 | 133 | AliCDBMetaData* md = new AliCDBMetaData(); |
134 | md->SetResponsible("Silvia Arcelli"); | |
6fce62af | 135 | md->SetComment("Full misalignment for TOF"); |
5bd470e1 | 136 | md->SetAliRootVersion(gSystem->Getenv("ARVERSION")); |
dfe9c69d | 137 | AliCDBId id("TOF/Align/Data",0,AliCDBRunRange::Infinity()); |
1895a097 | 138 | storage->Put(array,id,md); |
139 | } | |
140 | ||
141 | array->Delete(); | |
142 | ||
143 | } | |
144 | ||
145 |