Added functionality to generate module misalignments

[u/mrichter/AliRoot.git] / MUON / MUONCheckMisAligner.C

diff --git a/MUON/MUONCheckMisAligner.C b/MUON/MUONCheckMisAligner.C
index 21fc8f53d680d3f3f8ead99060fa2e0f4d7a7e7a..99301f39c9327dc0063aa0d990c13753b4396e95 100644 (file)
--- a/MUON/MUONCheckMisAligner.C
+++ b/MUON/MUONCheckMisAligner.C
@@ -38,6 +38,11 @@
   in this plane. Misalignments in the XZ and YZ plane will be very small 
   compared to those in the XY plane, which are small already - of the order 
   of microns)
+  --> Default behavior generates a "residual" misalignment using gaussian
+  distributions. Uniform distributions can still be used, see 
+  AliMUONGeometryAligner.
+  --> User can also generate module misalignments using SetModuleCartMisAlig
+  and SetModuleAngMisAlig
   Note : If the detection elements are allowed to be misaligned in all
   directions, this has consequences for the alignment algorithm, which 
   needs to know the number of free parameters. Eric only allowed 3 : 
@@ -48,22 +53,33 @@
 
 */
 
-void MUONCheckMisAligner(Double_t cartmisalig = 1.0, Double_t angmisalig = 10.)
+void MUONCheckMisAligner(Double_t xcartmisaligm = 0.0, Double_t xcartmisaligw = 0.004, 
+                        Double_t ycartmisaligm = 0.0, Double_t ycartmisaligw = 0.003, 
+                        Double_t angmisaligm = 0.0, Double_t angmisaligw = 0.0023,
+                        TString nameCDB = "ResMisAlignCDB")
 {
   
   AliMUONGeometryTransformer *transform = new AliMUONGeometryTransformer(true);
   transform->ReadGeometryData("volpath.dat", "transform.dat");
-
-  AliMUONGeometryMisAligner misAligner(cartmisalig,angmisalig);
+  if (gSystem->AccessPathName("geometry.root",kFileExists))
+    gGeoManager->Export("geometry.root");
+  AliMUONGeometryMisAligner misAligner(xcartmisaligm,xcartmisaligw,
+                                       ycartmisaligm,ycartmisaligw,
+                                      angmisaligm,angmisaligw);
 
   // Generate mis alignment data
+  // Uncomment lines below if you would like to generate module misalignments
+//  misAligner.SetModuleCartMisAlig(0.0,0.1,0.0,0.1,0.0,0.1); // Full
+//   misAligner.SetModuleAngMisAlig(0.0,0.02,0.0,0.04,0.0,0.02); // Full
+//   misAligner.SetModuleCartMisAlig(0.0,0.003,0.0,0.003,0.0,0.003); // Res
+//   misAligner.SetModuleAngMisAlig(0.0,0.0006,0.0,0.001,0.0,0.0005); // Res
   AliMUONGeometryTransformer *newTransform = misAligner.MisAlign(transform,true); 
   newTransform->WriteTransformations("transform2.dat");
   newTransform->WriteMisAlignmentData("misalign.root");
 
   // Apply misAlignment via AliRoot framework
   TGeoManager::Import("geometry.root");
-  AliRun::ApplyDisplacements(
+  AliSimulation::ApplyAlignObjsToGeom(
      const_cast<TClonesArray*>(newTransform->GetMisAlignmentData()));
   // Save new geometry file
   gGeoManager->Export("geometry2.root");
@@ -74,6 +90,21 @@ void MUONCheckMisAligner(Double_t cartmisalig = 1.0, Double_t angmisalig = 10.)
   transform3->WriteTransformations("transform3.dat");
                // Check that transform3.dat is equal to transform2.dat
 
+  // Generate misaligned data in local cdb
+  TClonesArray* array = newTransform->GetMisAlignmentData();
+   
+  TString sLocCDB("local://");
+  sLocCDB += nameCDB;
+  // CDB manager
+  AliCDBManager* cdbManager = AliCDBManager::Instance();
+  cdbManager->SetDefaultStorage(sLocCDB.Data());
+  
+  AliCDBMetaData* cdbData = new AliCDBMetaData();
+  cdbData->SetResponsible("Dimuon Offline project");
+  cdbData->SetComment("MUON alignment objects with residual misalignment");
+  AliCDBId id("MUON/Align/Data", 0, 0); 
+  cdbManager->Put(array, id, cdbData);
+
   // To run simulation with misaligned geometry, you have to set
   // the Align option in Config.C:
   // MUON->SetAlign("transform2.dat");