1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
16 /** @file AliFMDAlignFaker.cxx
17 @author Christian Holm Christensen <cholm@nbi.dk>
18 @date Sun Mar 26 17:57:55 2006
19 @brief Implementation of AliFMDAlignFaker
21 //____________________________________________________________________
30 //____________________________________________________________________
32 // Forward Multiplicity Detector based on Silicon wafers.
34 // This task creates fake alignment. Which alignment, depends on
35 // the bit mask passed to the constructor, or added by `AddAlign'.
37 // The default is to write all alignment parameters to a local
38 // storage `local://cdb' which is a directory in the current
41 #include "AliLog.h" // ALILOG_H
42 #include "AliFMDAlignFaker.h" // ALIFMDALIGNFAKER_H
43 #include <AliCDBManager.h> // ALICDBMANAGER_H
44 #include <AliCDBEntry.h> // ALICDBMANAGER_H
45 // #include <AliAlignObj.h>
46 #include <AliAlignObjAngles.h>
47 // #include <Riostream.h>
48 // #include <TSystem.h>
51 #include <TClonesArray.h>
54 #include <TGeoManager.h>
56 // #include <TGeoVolume.h>
59 //====================================================================
60 ClassImp(AliFMDAlignFaker)
62 ; // This is here to keep Emacs for indenting the next line
65 //____________________________________________________________________
66 AliFMDAlignFaker::AliFMDAlignFaker(Int_t mask, const char* geo,
74 // Default constructor
75 SetSensorDisplacement();
77 SetHalfDisplacement();
82 //__________________________________________________________________
84 AliFMDAlignFaker::SetSensorDisplacement(Double_t x1, Double_t y1, Double_t z1,
85 Double_t x2, Double_t y2, Double_t z2)
87 // Set sensor displacement (unit is centimeters)
88 fSensorTransMin.SetXYZ(x1, y1, z1);
89 fSensorTransMax.SetXYZ(x2, y2, z2);
92 //__________________________________________________________________
94 AliFMDAlignFaker::SetSensorRotation(Double_t x1, Double_t y1, Double_t z1,
95 Double_t x2, Double_t y2, Double_t z2)
97 // Set sensor rotations (unit is degrees)
98 fSensorRotMin.SetXYZ(x1, y1, z1);
99 fSensorRotMax.SetXYZ(x2, y2, z2);
102 //__________________________________________________________________
104 AliFMDAlignFaker::SetHalfDisplacement(Double_t x1, Double_t y1, Double_t z1,
105 Double_t x2, Double_t y2, Double_t z2)
107 // Set half ring/cone displacement (unit is centimeters)
108 fHalfTransMin.SetXYZ(x1, y1, z1);
109 fHalfTransMax.SetXYZ(x2, y2, z2);
112 //__________________________________________________________________
114 AliFMDAlignFaker::SetHalfRotation(Double_t x1, Double_t y1, Double_t z1,
115 Double_t x2, Double_t y2, Double_t z2)
117 // Set half ring/cone rotations (unit is degrees)
118 fHalfRotMin.SetXYZ(x1, y1, z1);
119 fHalfRotMax.SetXYZ(x2, y2, z2);
122 //__________________________________________________________________
123 #define IS_NODE_HALF(name) \
124 (name[0] == 'F' && name[2] == 'M' && (name[3] == 'T' || name[3] == 'B'))
125 #define IS_NODE_SENSOR(name) \
126 (name[0] == 'F' && name[2] == 'S' && name[3] == 'E')
128 //__________________________________________________________________
130 AliFMDAlignFaker::Exec(Option_t*)
136 if (!TGeoManager::Import(GetName())) {
137 AliFatal(Form("Failed to import geometry from %s", GetName()));
142 TGeoVolume* topVolume = gGeoManager->GetTopVolume();
144 AliFatal("No top-level volume defined");
147 // Make container of transforms
148 if (!fArray) fArray = new TClonesArray("AliAlignObjAngles");
152 TGeoIterator next(topVolume);
155 // Loop over all entries in geometry to find our nodes.
156 while ((node = static_cast<TGeoNode*>(next()))) {
157 const char* name = node->GetName();
158 if (!(IS_NODE_HALF(name) && TESTBIT(fMask, kHalves)) &&
159 !(IS_NODE_SENSOR(name) && TESTBIT(fMask, kSensors)))
163 TString path(Form("/%s", gGeoManager->GetNode(0)->GetName()));
164 Int_t nLevel = next.GetLevel();
165 for (Int_t lvl = 0; lvl <= nLevel; lvl++) {
166 TGeoNode* p = next.GetNode(lvl);
169 AliWarning(Form("No node at level %d in path %s",lvl,path.Data()));
172 if (!path.IsNull()) path.Append("/");
173 path.Append(p->GetName());
175 Int_t id = node->GetVolume()->GetNumber();
176 if (IS_NODE_HALF(name)) MakeAlignHalf(path, id);
177 if (IS_NODE_SENSOR(name)) MakeAlignSensor(path, id);
180 TString t(GetTitle());
181 if (t.IsNull() || t.Contains("local://") || t.Contains("alien://"))
187 //__________________________________________________________________
189 AliFMDAlignFaker::MakeAlign(const TString& path, Int_t id,
190 Double_t transX, Double_t transY, Double_t transZ,
191 Double_t rotX, Double_t rotY, Double_t rotZ)
193 // make alignment for a path
197 // transX Translation in X
198 // transZ Translation in Y
199 // transZ Translation in Z
200 // rotX Rotation about X-axis
201 // rotY Rotation about Y-axis
202 // rotZ Rotation about Z-axis
203 AliDebug(1, Form("Make alignment for %s (volume %d): (%f,%f,%f) (%f,%f,%f)",
204 path.Data(), id, transX, transY, transZ, rotX, rotY, rotZ));
205 Int_t nAlign = fArray->GetEntries();
207 AliAlignObjAngles* obj =
208 new ((*fArray)[nAlign]) AliAlignObjAngles(path.Data(), id,0,0,0,0,0,0);
210 AliError(Form("Failed to create alignment object for %s", path.Data()));
213 if (!obj->SetLocalPars(transX, transY, transZ, rotX, rotY, rotZ)) {
214 AliError(Form("Failed to set local transforms on %s", path.Data()));
220 //__________________________________________________________________
222 AliFMDAlignFaker::MakeAlignHalf(const TString& path, Int_t id)
224 // Make alignment of a half ring/cone
225 AliDebug(15, Form("Make alignment for half-ring/cone %s", path.Data()));
226 Double_t transX = gRandom->Uniform(fHalfTransMin.X(), fHalfTransMax.X());
227 Double_t transY = gRandom->Uniform(fHalfTransMin.Y(), fHalfTransMax.Y());
228 Double_t transZ = gRandom->Uniform(fHalfTransMin.Z(), fHalfTransMax.Z());
229 Double_t rotX = gRandom->Uniform(fHalfRotMin.X(), fHalfRotMax.X());
230 Double_t rotY = gRandom->Uniform(fHalfRotMin.Y(), fHalfRotMax.Y());
231 Double_t rotZ = gRandom->Uniform(fHalfRotMin.Z(), fHalfRotMax.Z());
232 return MakeAlign(path, id, transX, transY, transZ, rotX, rotY, rotZ);
236 //__________________________________________________________________
238 AliFMDAlignFaker::MakeAlignSensor(const TString& path, Int_t id)
240 // Make alignment of a sensor
241 AliDebug(15, Form("Make alignment for sensor %s", path.Data()));
242 Double_t transX = gRandom->Uniform(fSensorTransMin.X(), fSensorTransMax.X());
243 Double_t transY = gRandom->Uniform(fSensorTransMin.Y(), fSensorTransMax.Y());
244 Double_t transZ = gRandom->Uniform(fSensorTransMin.Z(), fSensorTransMax.Z());
245 Double_t rotX = gRandom->Uniform(fSensorRotMin.X(), fSensorRotMax.X());
246 Double_t rotY = gRandom->Uniform(fSensorRotMin.Y(), fSensorRotMax.Y());
247 Double_t rotZ = gRandom->Uniform(fSensorRotMin.Z(), fSensorRotMax.Z());
248 return MakeAlign(path, id, transX, transY, transZ, rotX, rotY, rotZ);
251 //__________________________________________________________________
253 AliFMDAlignFaker::WriteToCDB()
256 AliCDBManager* cdb = AliCDBManager::Instance();
257 AliCDBMetaData* meta = new AliCDBMetaData;
258 meta->SetResponsible(gSystem->GetUserInfo()->fRealName.Data());
259 meta->SetAliRootVersion(gROOT->GetVersion());
260 meta->SetBeamPeriod(1);
261 meta->SetComment(fComment.Data());
263 AliCDBId id("FMD/Align/Data", fRunMin, fRunMax);
264 cdb->Put(fArray, id, meta);
268 //__________________________________________________________________
270 AliFMDAlignFaker::WriteToFile()
272 // Write to a local file
273 TFile* file = TFile::Open(GetTitle(), "RECREATE");
275 AliFatal(Form("Failed to open file '%s' for output", GetTitle()));
279 fArray->Write("FMDAlignment");
286 //____________________________________________________________________