MUON/runDataReconstruction.C

   1 /**************************************************************************
   2  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   3  *                                                                        *
   4  * Author: The ALICE Off-line Project.                                    *
   5  * Contributors are mentioned in the code where appropriate.              *
   6  *                                                                        *
   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  **************************************************************************/
  15
  16 /* $Id: runReconstruction.C 23207 2007-12-20 09:59:20Z ivana $ */
  17
  18 /// \ingroup macros
  19 /// \file runDataReconstruction.C
  20 /// \brief Macro for running reconstruction
  21 ///
  22 /// Macro for running reconstruction on the cosmics run data.
  23 ///
  24 /// \author Laurent Aphecetche, Nicole Bastid, Bogdan Vulpescu, ...
  25
  26 #if !defined(__CINT__) || defined(__MAKECINT__)
  27 #include "AliMUONReconstructor.h"
  28 #include "AliMUONRecoParam.h"
  29 #include "AliRecoParam.h"
  30 #include "AliCDBManager.h"
  31 #include "AliMagFMaps.h"
  32 #include "AliTracker.h"
  33 #include "AliReconstruction.h"
  34 #include <TRandom.h>
  35 #include <TGrid.h>
  36 //#include <TObjectTable.h>
  37 #endif
  38
  39 // Data file, OCDB on Grid
  40 TString input="alien:///alice/data/2008/LHC08a/000021931/raw/08000021931001.50.root";
  41 TString ocdbPath = "alien://folder=/alice/data/2008/LHC08a/OCDB";
  42
  43 // Data file, OCDB locally
  44 //TString input="$ALICE_ROOT/MUON/test_out.100/raw.root";
  45 //TString ocdbPath = "local://$ALICE_ROOT";
  46
  47 TString caliboption1 = "NOGAIN";
  48 TString caliboption2 = "GAINCONSTANTCAPA";
  49 TString recoptions = "SAVEDIGITS";
  50 Int_t seed = 1234567;
  51
  52 void runDataReconstruction(Int_t calib = 1)
  53 {
  54   TGrid::Connect("alien://");
  55
  56   AliCDBManager* man = AliCDBManager::Instance();
  57   man->SetDefaultStorage(ocdbPath.Data());
  58   //man->SetSpecificStorage("MUON/Calib/Mapping","local://$ALICE_ROOT");
  59   //man->SetSpecificStorage("MUON/Calib/DDLStore","local://$ALICE_ROOT");
  60   //man->SetSpecificStorage("MUON/Calib/Gains","local://$ALICE_ROOT");
  61
  62   gRandom->SetSeed(seed);
  63
  64   AliMagFMaps* field = new AliMagFMaps("Maps","Maps", 1, 1., 10., AliMagFMaps::k5kG);
  65   AliTracker::SetFieldMap(field, kFALSE);
  66   AliReconstruction *MuonRec = new AliReconstruction();
  67
  68
  69   MuonRec->SetInput(input.Data());
  70   MuonRec->SetRunVertexFinder(kFALSE);
  71   MuonRec->SetRunLocalReconstruction("MUON");
  72   MuonRec->SetRunTracking("MUON");
  73   MuonRec->SetFillESD(" ");
  74   MuonRec->SetLoadAlignData("MUON");
  75   MuonRec->SetNumberOfEventsPerFile(1000);
  76   MuonRec->SetOption("MUON",recoptions.Data());
  77   //  MuonRec->SetEventRange(319,319);
  78   MuonRec->SetWriteAOD();
  79   //MuonRec.SetEventRange(0,100);
  80   AliMUONRecoParam *muonRecoParam = AliMUONRecoParam::GetLowFluxParam();
  81   muonRecoParam->CombineClusterTrackReco(kFALSE);
  82   TString caliboption = caliboption1;
  83   if ( calib == 2 ) caliboption = caliboption2;
  84   muonRecoParam->SetCalibrationMode(caliboption.Data());
  85   //muonRecoParam->SetClusteringMode("PEAKFIT");
  86   //muonRecoParam->SetClusteringMode("PEAKCOG");
  87   muonRecoParam->Print("FULL");
  88   AliRecoParam::Instance()->RegisterRecoParam(muonRecoParam);
  89
  90   MuonRec->Run();
  91
  92   delete MuonRec;
  93
  94   //gObjectTable->Print();
  95 }
  96