]>
Commit | Line | Data |
---|---|---|
6379d2b0 | 1 | |
2 | void MakePipeMisalignment(Int_t inputRun=0, | |
3 | const char* ocdbPath="local://$ALICE_ROOT/OCDB", | |
4 | Bool_t verbose=true, | |
5 | const char* outFileName=0) | |
6 | { | |
7 | // Creates misalignment of the beam pipe, frame and FMD2&3. | |
8 | // | |
9 | // The misalignment is based on the ITS global alignment object with | |
10 | // the exception of the beam pipe where only the translations are | |
11 | // used ignoring the rotation matrix. | |
12 | AliCDBManager* cdb = AliCDBManager::Instance(); | |
b1ce649d | 13 | cdb->SetDefaultStorage("alien://Folder=/alice/data/2010/OCDB"); |
6379d2b0 | 14 | cdb->SetRun(inputRun); |
15 | ||
16 | ||
17 | // Get the geometry | |
18 | AliGeomManager::LoadGeometry(); | |
19 | ||
20 | // Get ITS alignment | |
21 | AliCDBEntry* itsEntry = cdb->Get("ITS/Align/Data"); | |
22 | if (!itsEntry) { | |
23 | Error("MakePipeMisalignment", "Couldn't get ITS alignment data"); | |
24 | return; | |
25 | } | |
26 | TClonesArray* itsArray = static_cast<TClonesArray*>(itsEntry->GetObject()); | |
27 | if (!itsArray) { | |
28 | Error("MakePipeMisalignment", "No array in CDB entry"); | |
29 | return; | |
30 | } | |
31 | AliAlignObjParams* itsAlign = | |
32 | static_cast<AliAlignObjParams*>(itsArray->At(0)); | |
33 | if (!itsAlign) { | |
34 | Error("MakePipeMisalignment", "No alignment object in CDB entry"); | |
35 | return; | |
36 | } | |
37 | ||
38 | if (verbose) { | |
39 | Info("MakePipeMisalignment", "ITS alignment:"); | |
40 | itsAlign->Print(); | |
41 | } | |
42 | ||
43 | Double_t itsTrans[3]; | |
44 | Double_t itsRot[3]; | |
45 | itsAlign->GetTranslation(itsTrans); | |
46 | itsAlign->GetAngles(itsRot); | |
47 | ||
48 | Double_t local[] = { 0, 0, 0 }; | |
49 | ||
50 | if (!gGeoManager->cd("/ALIC_1/CP_1/Cp1_1/CpSupC_1")) return; | |
51 | Double_t pCpSupCA[3]; | |
52 | gGeoManager->LocalToMaster(local, pCpSupCA); | |
53 | ||
54 | if (!gGeoManager->cd("/ALIC_1/CP_1/Cp3_1/CpSupC_3")) return; | |
55 | Double_t pCpSupCC[3]; | |
56 | gGeoManager->LocalToMaster(local, pCpSupCC); | |
57 | ||
58 | if (!gGeoManager->cd("/ALIC_1/CP_1/Cp1_1")) return; | |
59 | Double_t pCp1[3]; | |
60 | gGeoManager->LocalToMaster(local, pCp1); | |
61 | ||
62 | if (!gGeoManager->cd("/ALIC_1/CP_1/Cp1_1/CP1MO_1/CP1AT_1")) return; | |
63 | do { | |
64 | TGeoNode* node = gGeoManager->GetCurrentNode(); | |
65 | if (!node) break; | |
66 | TGeoVolume* volume = node->GetVolume(); | |
67 | if (!volume) break; | |
68 | TGeoBBox* bbox = static_cast<TGeoBBox*>(volume->GetShape()); | |
69 | if (!bbox) break; | |
70 | local[2] = -bbox->GetDZ(); | |
71 | } while (false); | |
72 | if (local[2] == 0) { | |
73 | Error("MakePipeMisalignement", "Failed to find end of CP1 section"); | |
74 | return; | |
75 | } | |
76 | Double_t pCP1AT[3]; | |
77 | gGeoManager->LocalToMaster(local, pCP1AT); | |
78 | ||
79 | ||
80 | ||
81 | // Now, we need to misalign the C-side collar according to the ITS | |
82 | // misalignment. | |
83 | TGeoHMatrix* itsMat = new TGeoHMatrix(); | |
84 | itsAlign->GetMatrix(*itsMat); | |
85 | // itsTrans[2] = 0; | |
86 | // itsMat->SetTranslation(itsTrans); | |
87 | Double_t tCpSupCA[3]; | |
88 | itsMat->LocalToMaster(pCpSupCA, tCpSupCA); | |
89 | Double_t tCpSupCC[3]; | |
90 | itsMat->LocalToMaster(pCpSupCC, tCpSupCC); | |
91 | ||
92 | Info("MakePipeMisalignment", "\n" | |
93 | "\tCollar A-side (%f,%f,%f)\n" | |
94 | "\tCollar C-side (%f,%f,%f)\n" | |
95 | "\tBeam pipe center (%f,%f,%f)\n" | |
96 | "\tEnd of CP1 A-side (%f,%f,%f)\n" | |
97 | "\tCollar A-side after (%f,%f,%f)\n" | |
98 | "\tCollar C-side after (%f,%f,%f)", | |
99 | pCpSupCA[0], pCpSupCA[1], pCpSupCA[2], | |
100 | pCpSupCC[0], pCpSupCC[1], pCpSupCC[2], | |
101 | pCp1[0], pCp1[1], pCp1[2], | |
102 | pCP1AT[0], pCP1AT[1], pCP1AT[2], | |
103 | tCpSupCA[0], tCpSupCA[1], tCpSupCA[2], | |
104 | tCpSupCC[0], tCpSupCC[1], tCpSupCC[2]); | |
105 | ||
106 | // Distance from A-side end to mis-aligned C-side Collar | |
107 | Double_t dX = - tCpSupCC[0] + pCP1AT[0] ; | |
108 | Double_t dY = - tCpSupCC[1] + pCP1AT[1] ; | |
109 | Double_t dZ = - tCpSupCC[2] + pCP1AT[2] ; | |
110 | ||
111 | // Rotation angle around X and Y | |
112 | Double_t rotX = TMath::ATan2(dY, dZ); | |
113 | Double_t rotY = TMath::ATan2(dX, dZ); | |
114 | ||
115 | ||
116 | Double_t dZc = pCp1[2] - pCP1AT[2]; | |
117 | Double_t dXc = TMath::Tan(rotY) * dZc; | |
118 | Double_t dYc = TMath::Tan(rotX) * dZc; | |
119 | ||
120 | ||
121 | Double_t dZt = pCpSupCC[2] - pCP1AT[2] + itsTrans[2]; | |
122 | ||
123 | Info("MakePipeMisalignment", "\n" | |
124 | "\tDistances (%f,%f,%f)\n" | |
125 | "\trotX=%f rotY=%f\n" | |
126 | "\tCP1 displacement (x,y)=(%f,%f) dZ=%f\n" | |
127 | "\tCpSubCA displacement (x,y)=(%f,%f) dZ=%f", | |
128 | dX, dY, dZ, | |
129 | rotX * TMath::RadToDeg(), rotY * TMath::RadToDeg(), | |
130 | dXc, dYc, dZc, | |
131 | TMath::Tan(rotY)*dZt, TMath::Tan(rotX)*dZt, dZt); | |
132 | ||
133 | ||
134 | TGeoRotation* rot = new TGeoRotation(); | |
135 | rot->RotateX(rotX*TMath::RadToDeg()); | |
136 | rot->RotateY(rotY*TMath::RadToDeg()); | |
137 | TGeoCombiTrans* trans = new TGeoCombiTrans(dXc, dYc, 0, rot); | |
138 | trans->Print(); | |
139 | ||
140 | // Output array | |
141 | TClonesArray* structArray = new TClonesArray("AliAlignObjParams", 20); | |
142 | ||
143 | //dummy vol id | |
144 | UShort_t dvoluid = | |
145 | AliGeomManager::LayerToVolUID(AliGeomManager::kInvalidLayer,0); | |
146 | ||
147 | //base of symbolic name corresponding to base of path "ALIC_1/B077_1/BSEGMO"; | |
148 | const char* baseSymName = "FRAME/Sector"; | |
149 | for(Int_t sm=0; sm<18; sm++){ | |
150 | TString symname = baseSymName; | |
151 | symname += sm; | |
152 | new((*structArray)[sm]) AliAlignObjParams(symname.Data(),dvoluid, | |
153 | 0.,0.,0.,0.,0.,0.,kTRUE); | |
154 | } | |
155 | ||
156 | // --- Beam pipe --------------------------------------------------- | |
157 | AliAlignObjParams* c1Align = | |
158 | new((*structArray)[18]) AliAlignObjParams("CP1",dvoluid, | |
159 | dXc, dYc, 0, rotX, rotY, 0,kTRUE); | |
160 | c1Align->Print(); | |
161 | ||
162 | // --- C side Collar ----------------------------------------------- | |
163 | AliAlignObjParams* cSupCAlign = | |
164 | new((*structArray)[19]) AliAlignObjParams("CPSUPC", | |
165 | dvoluid, | |
166 | tCpSupCC[0]-pCpSupCC[0], | |
167 | tCpSupCC[1]-pCpSupCC[1], | |
168 | tCpSupCC[2]-pCpSupCC[2], | |
169 | 0, 0, 0, kTRUE); | |
170 | cSupCAlign->Print(); | |
171 | ||
172 | // --- save in CDB storage ----------------------------------------- | |
173 | Info("MakeForwardMisAlignment", | |
174 | "Saving alignment objects in CDB storage %s", ocdbPath); | |
175 | AliCDBStorage* storage = cdb->GetStorage(ocdbPath); | |
176 | if (!storage){ | |
177 | Error("MakeForwardMisAlignment","Unable to open storage %s\n",ocdbPath); | |
178 | return; | |
179 | } | |
180 | ||
181 | ||
182 | // --- Make objects and store -------------------------------------- | |
183 | AliCDBMetaData* md = new AliCDBMetaData(); | |
184 | md->SetResponsible("Grosso Raffaele"); | |
185 | md->SetComment("Misalignment for FRAME and beam pipe"); | |
186 | md->SetAliRootVersion(gSystem->Getenv("ARVERSION")); | |
187 | AliCDBId id("GRP/Align/Data",0,AliCDBRunRange::Infinity()); | |
188 | storage->Put(structArray,id,md); | |
189 | ||
190 | gGeoManager->Export("myGeom.root"); | |
191 | } |