]>
Commit | Line | Data |
---|---|---|
1895a097 | 1 | void MakePMDFullMisAlignment(){ |
2 | // Create TClonesArray of full misalignment objects for PMD | |
3 | // | |
4 | ||
5 | // Macro to randomly displace the 4 sectors of the PMD | |
6 | // in each plane. Each sector (to be misaligned) | |
7 | // of PMD houses the following : | |
8 | // (a) 6 modules of preshower plane | |
9 | // (b) 6 modules of veto plane | |
10 | // (c) The FEE boards on back plane of each module | |
11 | // (d) 6 modules of convertor plates | |
12 | // The clustering is done module-wise | |
13 | // The actual amount displacement will be provided | |
14 | // by the survey data and has to be converted into | |
15 | // displacement in x,y,z,theta, phi and psi | |
16 | ||
17 | ||
18 | // Now specify the path of the module to be misaligned | |
19 | // as followed in the PMD geant | |
20 | ||
21 | /* | |
22 | _____________ | |
23 | | | | | |
24 | | 1 | 3 | | |
25 | | |________| | |
26 | |____|___| | | |
27 | | | 2 | | |
28 | | 4 | | | |
29 | |________|____| | |
30 | ||
31 | // Misalignment Matrix is expected to be | |
32 | // same for sectors 1 and 4 | |
33 | // and for the sectors 2 and 3 | |
34 | // As these will be mounted on the same | |
35 | // Steel plate | |
36 | */ | |
37 | ||
6fce62af | 38 | const char* macroname = "MakePMDFullMisAlignment.C"; |
39 | ||
40 | // Activate CDB storage and load geometry from CDB | |
41 | AliCDBManager* cdb = AliCDBManager::Instance(); | |
162637e4 | 42 | if(!cdb->IsDefaultStorageSet()) cdb->SetDefaultStorage("local://$ALICE_ROOT/OCDB"); |
6fce62af | 43 | cdb->SetRun(0); |
44 | ||
45 | AliCDBStorage* storage; | |
46 | ||
a24be56b | 47 | if( TString(gSystem->Getenv("TOCDB")) == TString("kTRUE") ){ |
6fce62af | 48 | TString Storage = gSystem->Getenv("STORAGE"); |
49 | if(!Storage.BeginsWith("local://") && !Storage.BeginsWith("alien://")) { | |
50 | Error(macroname,"STORAGE variable set to %s is not valid. Exiting\n",Storage.Data()); | |
51 | return; | |
52 | } | |
53 | storage = cdb->GetStorage(Storage.Data()); | |
54 | if(!storage){ | |
55 | Error(macroname,"Unable to open storage %s\n",Storage.Data()); | |
56 | return; | |
57 | } | |
58 | AliCDBPath path("GRP","Geometry","Data"); | |
59 | AliCDBEntry *entry = storage->Get(path.GetPath(),cdb->GetRun()); | |
60 | if(!entry) Fatal(macroname,"Could not get the specified CDB entry!"); | |
61 | entry->SetOwner(0); | |
62 | TGeoManager* geom = (TGeoManager*) entry->GetObject(); | |
63 | AliGeomManager::SetGeometry(geom); | |
64 | }else{ | |
65 | AliGeomManager::LoadGeometry(); //load geom from default CDB storage | |
dfe9c69d | 66 | } |
1895a097 | 67 | |
68 | Float_t max_trans=0.1; | |
69 | Float_t max_rot=0.1; | |
70 | ||
d3a9c52a | 71 | const char *Sector1="PMD/Sector1"; |
72 | const char *Sector2="PMD/Sector2"; | |
73 | const char *Sector3="PMD/Sector3"; | |
74 | const char *Sector4="PMD/Sector4"; | |
1895a097 | 75 | |
76 | ||
77 | //Sectors 1 and 4 | |
78 | Double_t dx14, dy14, dz14; // Misalignment in X,Y and Z | |
79 | Double_t dpsi14, dtheta14, dphi14; // Angular displacements | |
80 | //Sectors 2 and 3 | |
81 | Double_t dx23, dy23, dz23; // Misalignment in X,Y and Z | |
82 | Double_t dpsi23, dtheta23, dphi23; // Angular displacements | |
83 | ||
84 | ||
85 | ||
86 | // At present give some random shifts | |
87 | // generated by random numbers | |
88 | // max_trans : Maximun shifts in X,Y,Z in centimeters | |
89 | // max_rot : Maximum shifts in angles in degrees | |
90 | // Double_t max_trans, max_rot; | |
91 | ||
92 | ||
93 | TRandom *rnd = new TRandom(4357); | |
94 | ||
95 | // For sectors 1 and 4 | |
96 | // Translation displacement | |
97 | dx14 = (rnd->Uniform()-0.5)*max_trans; | |
98 | dy14 = (rnd->Uniform()-0.5)*max_trans; | |
99 | dz14 = (rnd->Uniform()-0.5)*max_trans; | |
100 | //Rotation angles | |
101 | dpsi14 = (rnd->Uniform()-0.5)*max_rot; | |
102 | dtheta14 = (rnd->Uniform()-0.5)*max_rot; | |
103 | dphi14 = (rnd->Uniform()-0.5)*max_rot; | |
104 | ||
105 | // For sectors 2 and 3 | |
106 | // Translation displacement | |
107 | dx23 = (rnd->Uniform()-0.5)*max_trans; | |
108 | dy23 = (rnd->Uniform()-0.5)*max_trans; | |
109 | dz23 = (rnd->Uniform()-0.5)*max_trans; | |
110 | //Rotation angles | |
111 | dpsi23 = (rnd->Uniform()-0.5)*max_rot; | |
112 | dtheta23 = (rnd->Uniform()-0.5)*max_rot; | |
113 | dphi23 = (rnd->Uniform()-0.5)*max_rot; | |
114 | ||
115 | ||
116 | //Create a TClonesArray of Align Object to store displacement Angles | |
90dbf5fb | 117 | TClonesArray *array = new TClonesArray("AliAlignObjParams",10); |
1895a097 | 118 | TClonesArray &alobj = *array; |
119 | ||
1895a097 | 120 | Int_t iIndex=0; //let all modules have index=0 in a layer with no LUT |
ae079791 | 121 | AliGeomManager::ELayerID iLayer = AliGeomManager::kInvalidLayer; |
122 | UShort_t volid = AliGeomManager::LayerToVolUID(iLayer,iIndex); | |
1895a097 | 123 | |
90dbf5fb | 124 | new(alobj[0]) AliAlignObjParams(Sector1, volid, dx14, dy14, dz14, dpsi14, dtheta14, dphi14, kFALSE); |
125 | new(alobj[1]) AliAlignObjParams(Sector2, volid, dx14, dy14, dz14, dpsi14, dtheta14, dphi14, kFALSE); | |
126 | new(alobj[2]) AliAlignObjParams(Sector3, volid, dx23, dy23, dz23, dpsi23, dtheta23, dphi23, kFALSE); | |
127 | new(alobj[3]) AliAlignObjParams(Sector4, volid, dx23, dy23, dz23, dpsi23, dtheta23, dphi23, kFALSE); | |
1895a097 | 128 | |
a24be56b | 129 | if( TString(gSystem->Getenv("TOCDB")) != TString("kTRUE") ){ |
1895a097 | 130 | // Create a File to store the alignment data |
dfe9c69d | 131 | const char* filename = "PMDfullMisalignment.root"; |
132 | TFile f(filename,"RECREATE"); | |
133 | if(!f){ | |
134 | Error(macroname,"cannot open file for output\n"); | |
135 | return; | |
136 | } | |
137 | Info(macroname,"Saving alignment objects to the file %s", filename); | |
1895a097 | 138 | f.cd(); |
dfe9c69d | 139 | f.WriteObject(array,"PMDAlignObjs","kSingleKey"); |
1895a097 | 140 | f.Close(); |
141 | }else{ | |
142 | // save in CDB storage | |
1895a097 | 143 | AliCDBMetaData* md = new AliCDBMetaData(); |
144 | md->SetResponsible(""); | |
145 | md->SetComment("Full misalignment for PMD, produced with sigmatr=0.1 and sigmarot=0.1 in the local RS"); | |
5bd470e1 | 146 | md->SetAliRootVersion(gSystem->Getenv("ARVERSION")); |
dfe9c69d | 147 | AliCDBId id("PMD/Align/Data",0,AliCDBRunRange::Infinity()); |
1895a097 | 148 | storage->Put(array,id,md); |
149 | } | |
150 | array->Delete(); | |
151 | ||
152 | } |