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 }