Correct function Compare() for "pixels" from MLEM cluster finder.

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

/**************************************************************************
 * 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 residual misalignment data.
// The macro is triggered from AliRoot/macros/MakeAllDETsResMisAlignment.C
//
//  Author: I. Hrivnacova, IPN Orsay

#if !defined(__CINT__) || defined(__MAKECINT__)

#include "AliMUONGeometryTransformer.h"
#include "AliMUONGeometryMisAligner.h"

#include "AliCDBManager.h"
#include "AliCDBStorage.h"
#include "AliCDBId.h"

#include <TGeoManager.h>
#include <TClonesArray.h>
#include <TString.h>
#include <TFile.h>
#include <Riostream.h>

#endif

void MakeMUONResMisAlignment()
{
  // Check first if geometry is loaded,
  // if not loaded try to load it from galice.root file
  if ( ! gGeoManager && ! TGeoManager::Import("geometry.root") )  {
    cerr << "Loading geometry failed." << endl;
    return;
  }  

  AliMUONGeometryTransformer transformer(true);
  transformer.ReadGeometryData("volpath.dat", gGeoManager);
  
  AliMUONGeometryMisAligner misAligner(0.0, 0.004, 0.0, 0.003, 0.0, 0.0023);
  AliMUONGeometryTransformer* newTransform 
    = misAligner.MisAlign(&transformer, true);
  const TClonesArray* array = newTransform->GetMisAlignmentData();

  if ( TString(gSystem->Getenv("TOCDB")) != TString("kTRUE") ) {
    cout << "Generating residual misalignment data in a file" << endl;
  
    // Create a File to store the alignment data
    TFile f("MUONresidualMisalignment.root","RECREATE");
    if( !f.IsOpen() ) {
      cerr << "cannot open file for output" << endl;
    }
    
    f.cd();
    f.WriteObject(array,"MUONAlignObjs ","kSingleKey");
    f.Close();
  }
  else {
    cout << "Generating residual misalignment data in CDB" << endl;
    // save in CDB storage
    const char* Storage = gSystem->Getenv("STORAGE");
    AliCDBManager* cdbManager = AliCDBManager::Instance();
    AliCDBStorage* storage = cdbManager->GetStorage(Storage);
    AliCDBMetaData* cdbData = new AliCDBMetaData();
    cdbData->SetResponsible("Dimuon Offline project");
    cdbData->SetComment("MUON alignment objects with residual misalignment");
    cdbData->SetAliRootVersion(gSystem->Getenv("ARVERSION"));
    AliCDBId id("MUON/Align/Data", 0, 9999999); 
    storage->Put(const_cast<TClonesArray*>(array), id, cdbData);
  }
  delete newTransform;
 }