20-sep-2005 NvE Dummy values 0 set for KAPPA-ADC and KAPPA-TOT for OMs in IceCal2Root
[u/mrichter/AliRoot.git] / RALICE / icepack / macros / icecalib.cc
1 //////////////////////////////////////////////////////////
2 // Example macro to demonstrate usage of the IceCalibrate
3 // processor for various calibrations as a subtask of
4 // an F2K conversion job. The macro also shows how one
5 // can (interactively) invoke one or more subtasks
6 // (i.e. IceXtalk and EvtAna.cxx) to be executed.
7 // The latter is very convenient in developing/testing
8 // new reconstruction/analysis algorithms.
9 //
10 // To run this macro in batch, just do
11 //
12 // root -b -q icecalib.cc
13 //
14 // For more details see the docs of class IceCalibrate
15 //
16 // NvE 20-sep-2005 Utrecht University
17 //////////////////////////////////////////////////////////
18 {
19  gSystem->Load("ralice");
20  gSystem->Load("icepack");
21  gSystem->Load("iceconvert");
22
23  // Interactively compile and load the EvtAna.cxx code
24  gROOT->LoadMacro("EvtAna.cxx+");
25
26  // The database loading job (needed for the Xtalk constants)
27  IceCal2Root cal("IceCal2Root","Calibration format conversion");
28  cal.SetAmacalibFile("amacalib_amanda2_2003.txt");
29  cal.ExecuteJob();
30
31  AliObjMatrix* omdb=cal.GetOMdbase();
32
33  // The calibration processor task
34  IceCalibrate calib("IceCalibrate","Signal calibration");
35  calib.SetOMdbase(omdb);
36
37  // The Xtalk correction processor task
38  IceXtalk xtalk("IceXtalk","Cross talk correction");
39  xtalk.SetOMdbase(omdb);
40  xtalk.SetMinProb(0.5);
41  xtalk.SetXtalkPE(1);
42
43  // The event analysis task
44  EvtAna evtana("evtana","Event analysis");
45
46  // The F2K event data processing job
47  IceF2k q("IceF2k","Processing of the F2K event data");
48  q.SetMaxEvents(1);
49  q.SetPrintFreq(0);
50  q.SetInputFile("real-reco.f2k");
51
52  q.Add(&calib);
53  q.Add(&xtalk);
54  q.Add(&evtana);
55
56  // Perform the conversion
57  q.ExecuteJob();
58 }