/**************************************************************************
 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/

// $Id$
//
// Macro for generating the geometry data files:
// (volpath.dat, transform.dat, svmap.dat)
// and local CDB storage with zero, residual and full misalignment
// To be run from aliroot:
// .x MUONGenerateGeometryData.C
//
// The generated files do not replace the existing ones
// but have different names (with extension ".out").
//
//  Author: I. Hrivnacova, IPN Orsay

void MUONGenerateGeometryData(Bool_t volpaths = true,
                              Bool_t transforms = true, 
                              Bool_t svmaps = true,
			      Bool_t zeroAlign = true,
			      Bool_t resMisAlign = true,
			      Bool_t fullMisAlign = true)
{
  // Initialize
  gAlice->Init("$ALICE_ROOT/MUON/Config.C");
  cout << "Init done " << endl;

  // Get MUON detector
  AliMUON* muon = (AliMUON*)gAlice->GetModule("MUON");
  if (!muon) {
    AliFatal("MUON detector not defined.");
    return 0;
  }  

  // Get geometry builder
  AliMUONGeometryBuilder* builder = muon ->GetGeometryBuilder();
  
  if (volpaths) {
    cout << "Generating volpath file ..." << endl;
    builder->GetTransformer()->WriteVolumePaths("volpath.dat.out");
  }  

  if (transforms) {
    cout << "Generating transformation file ..." << endl;
    builder->GetTransformer()->WriteTransformations("transform.dat.out");
  }  

  if (svmaps) {
    cout << "Generating svmaps file ..." << endl;
    builder->WriteSVMaps();
  }  

  if (zeroAlign) {
    cout << "Generating zero misalignment data in  MUON/Align/Data..." << endl;
   
    // Create zero alignment data
    TClonesArray* array 
      = builder->GetTransformer()->CreateZeroAlignmentData();
   
    // CDB manager
    AliCDBManager* cdbManager = AliCDBManager::Instance();
    cdbManager->SetDefaultStorage("local://$ALICE_ROOT");
  
    AliCDBMetaData* cdbData = new AliCDBMetaData();
    cdbData->SetResponsible("Dimuon Offline project");
    cdbData->SetComment("MUON alignment objects for ideal geometry");
    AliCDBId id("MUON/Align/Data", 0, 0); 
    cdbManager->Put(array, id, cdbData);
    
    delete array;
  } 
  
  if (resMisAlign) {
    cout << "Generating residual misalignment data in  MUON/ResMisAlignCDB/Data..." << endl;
   
    AliMUONGeometryMisAligner misAligner(0.0, 0.004, 0.0, 0.003, 0.0, 0.0023);
    AliMUONGeometryTransformer* newTransform 
      = misAligner.MisAlign(builder->GetTransformer(), true);
    TClonesArray* array = newTransform->GetMisAlignmentData();

    AliCDBManager* cdbManager = AliCDBManager::Instance();
    cdbManager->SetDefaultStorage("local://ResMisAlignCDB");
  
    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);
    
    delete newTransform;
  }   
           
  if (fullMisAlign) {
    cout << "Generating residual misalignment data in  MUON/FullMisAlignCDB/Data..." << endl;
   
    AliMUONGeometryMisAligner misAligner(0.0, 0.03, 0.0, 0.03, 0.0, 0.03);
    AliMUONGeometryTransformer* newTransform 
      = misAligner.MisAlign(builder->GetTransformer(), true);
    TClonesArray* array = newTransform->GetMisAlignmentData();

    AliCDBManager* cdbManager = AliCDBManager::Instance();
    cdbManager->SetDefaultStorage("local://FullMisAlignCDB");
  
    AliCDBMetaData* cdbData = new AliCDBMetaData();
    cdbData->SetResponsible("Dimuon Offline project");
    cdbData->SetComment("MUON alignment objects with full misalignment");
    AliCDBId id("MUON/Align/Data", 0, 0); 
    cdbManager->Put(array, id, cdbData);
    
    delete newTransform;
  }   
           

}