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 **************************************************************************/
19 /// \file MakeMUONResMisAlignment.C
20 /// \brief Macro for generating the residual misalignment data.
22 /// The macro is triggered from AliRoot/macros/MakeAllDETsResMisAlignment.C
24 /// \author: I. Hrivnacova, IPN Orsay
26 #if !defined(__CINT__) || defined(__MAKECINT__)
28 #include "AliMUONGeometryTransformer.h"
29 #include "AliMUONGeometryMisAligner.h"
31 #include "AliGeomManager.h"
32 #include "AliCDBManager.h"
33 #include "AliCDBStorage.h"
34 #include "AliCDBEntry.h"
39 #include <TClonesArray.h>
42 #include <Riostream.h>
46 void MakeMUONResMisAlignment()
48 const char* macroname = "MakeMUONResMisAlignment.C";
49 // Activate CDB storage and load geometry from CDB
50 AliCDBManager* cdb = AliCDBManager::Instance();
51 if(!cdb->IsDefaultStorageSet()) cdb->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
54 AliCDBStorage* storage = 0;
56 if( TString(gSystem->Getenv("TOCDB")) == TString("kTRUE") ){
57 TString Storage = gSystem->Getenv("STORAGE");
58 if(!Storage.BeginsWith("local://") && !Storage.BeginsWith("alien://")) {
59 Error(macroname,"STORAGE variable set to %s is not valid. Exiting\n",Storage.Data());
62 storage = cdb->GetStorage(Storage.Data());
64 Error(macroname,"Unable to open storage %s\n",Storage.Data());
67 AliCDBPath path("GRP","Geometry","Data");
68 AliCDBEntry *entry = storage->Get(path.GetPath(),cdb->GetRun());
69 if(!entry) Fatal(macroname,"Could not get the specified CDB entry!");
71 TGeoManager* geom = (TGeoManager*) entry->GetObject();
72 AliGeomManager::SetGeometry(geom);
74 AliGeomManager::LoadGeometry(); //load geom from default CDB storage
77 AliMUONGeometryTransformer transformer;
78 transformer.LoadGeometryData();
80 AliMUONGeometryMisAligner misAligner(0.0, 0.004, 0.0, 0.003, 0.0, 0.0023);
81 AliMUONGeometryTransformer* newTransform
82 = misAligner.MisAlign(&transformer, true);
83 const TClonesArray* array = newTransform->GetMisAlignmentData();
85 // 100 mum residual resolution for chamber misalignments?
86 misAligner.SetAlignmentResolution(array,-1,0.01,0.01);
88 if ( TString(gSystem->Getenv("TOCDB")) != TString("kTRUE") ) {
90 const char* filename = "MUONresidualMisalignment.root";
91 TFile f(filename,"RECREATE");
93 Error(macroname,"cannot open file for output\n");
96 Info(macroname,"Saving alignment objects to the file %s", filename);
98 f.WriteObject(array,"MUONAlignObjs","kSingleKey");
101 // Save in CDB storage
102 AliCDBMetaData* cdbData = new AliCDBMetaData();
103 cdbData->SetResponsible("Dimuon Offline project");
104 cdbData->SetComment("MUON alignment objects with residual misalignment");
105 cdbData->SetAliRootVersion(gSystem->Getenv("ARVERSION"));
106 AliCDBId id("MUON/Align/Data", 0, AliCDBRunRange::Infinity());
107 storage->Put(const_cast<TClonesArray*>(array), id, cdbData);