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 **************************************************************************/
16 /* $Id: runReconstruction.C 23207 2007-12-20 09:59:20Z ivana $ */
19 /// \file runDataReconstruction.C
20 /// \brief Macro for running reconstruction
22 /// Macro for running reconstruction on the cosmics run data.
24 /// \author Laurent Aphecetche, Nicole Bastid, Bogdan Vulpescu, ...
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"
36 //#include <TObjectTable.h>
39 // Data file, OCDB on Grid
40 #TString input="alien:///alice/data/2007/LHC07w/000014493/raw/07000014493001.10.root";
41 TString input="alien:///alice/data/2008/LHC08a/000021931/raw/08000021931001.50.root";
42 TString ocdbPath = "alien://folder=/alice/data/2008/LHC08a/OCDB";
44 // Data file, OCDB locally
45 //TString input="$ALICE_ROOT/MUON/test_out.100/raw.root";
46 //TString ocdbPath = "local://$ALICE_ROOT";
48 TString caliboption1 = "NOGAIN";
49 TString caliboption2 = "GAINCONSTANTCAPA";
50 TString recoptions = "SAVEDIGITS";
53 void runDataReconstruction(Int_t calib = 1)
55 TGrid::Connect("alien://");
57 AliCDBManager* man = AliCDBManager::Instance();
58 man->SetDefaultStorage(ocdbPath.Data());
59 //man->SetSpecificStorage("MUON/Calib/Mapping","local://$ALICE_ROOT");
60 //man->SetSpecificStorage("MUON/Calib/DDLStore","local://$ALICE_ROOT");
61 //man->SetSpecificStorage("MUON/Calib/Gains","local://$ALICE_ROOT");
63 gRandom->SetSeed(seed);
65 AliMagFMaps* field = new AliMagFMaps("Maps","Maps", 1, 1., 10., AliMagFMaps::k5kG);
66 AliTracker::SetFieldMap(field, kFALSE);
67 AliReconstruction *MuonRec = new AliReconstruction();
70 MuonRec->SetInput(input.Data());
71 MuonRec->SetRunVertexFinder(kFALSE);
72 MuonRec->SetRunLocalReconstruction("MUON");
73 MuonRec->SetRunTracking("MUON");
74 MuonRec->SetFillESD(" ");
75 MuonRec->SetLoadAlignData("MUON");
76 MuonRec->SetNumberOfEventsPerFile(1000);
77 MuonRec->SetOption("MUON",recoptions.Data());
78 // MuonRec->SetEventRange(319,319);
79 MuonRec->SetWriteAOD();
80 //MuonRec.SetEventRange(0,100);
81 AliMUONRecoParam *muonRecoParam = AliMUONRecoParam::GetLowFluxParam();
82 muonRecoParam->CombineClusterTrackReco(kFALSE);
83 TString caliboption = caliboption1;
84 if ( calib == 2 ) caliboption = caliboption2;
85 muonRecoParam->SetCalibrationMode(caliboption.Data());
86 //muonRecoParam->SetClusteringMode("PEAKFIT");
87 //muonRecoParam->SetClusteringMode("PEAKCOG");
88 muonRecoParam->Print("FULL");
89 AliRecoParam::Instance()->RegisterRecoParam(muonRecoParam);
95 //gObjectTable->Print();