ROOT logo
// $Id$
/*
 * Example macro to run the HLT tracker embedded into AliRoot reconstruction. 
 * The reconstruction is done from the TPC raw data. The calibration components 
 * have been added to the standard HLT TPC reconstruction.
 *
 * They take input from TPC-clusters, TPC-globalmerger and GLOBAL-esd-converter. 
 *
 * Usage:
 * <pre>
 * 
 * Usage: aliroot -b -q -l testTPCCalibration.C'("file", "cdb", minEvent, maxEvent, option)'
 *
 * Examples:
 *     testTPCCalibration.C'("raw://run12345", minEvent, MaxEvent,"calibtime")'
 *     testTPCCalibration.C'("raw.root", "local://$PWD", minEvent, MaxEvent, "calibtime")'
 * 
 * Defaults
 *     cdb="raw://"  -> takes OCDB from GRID
 *     minEvent=-1   -> no lower event selection
 *     maxEvent=-1   -> no upper event selection
 *     option="task" -> loads the analysis task that calls both the drift velocity and the gain calibration
 * 
 * </pre>
 *
 * To see usage examples run:
 * 
 * aliroot -q testTPCCalibration.C
 *
 * The argument option should be 'calibtime' for now, as the rest of the options 
 * call components that are not working yet.
 *
 * The reconstruction is steered by the AliReconstruction object in the
 * usual way.
 * 
 * @ingroup alihlt_tpc
 * @author Kalliopi.Kanaki@ift.uib.no
 */

// void testTPCCalibration(const char* filename, const char *cdbURI, int minEvent=-1, int maxEvent=-1, const char* option="task"){
// 
//   gSystem->Load("libANALYSIS");
//   gSystem->Load("libTPCcalib");  
//   
//   if(!gSystem->AccessPathName("galice.root")){
//     cerr << "Please delete file galice.root or run at a different place." << endl;
//     return;
//   }
// 
//   if (!filename) {
//     cerr << "Please specify input file name." << endl;
//     return;
//   }
// 
//   /////////////////////////////////////////////////////////////////////////
//   //
//   // init the HLT system in order to define the analysis chain below
//   //
//  
//   AliHLTSystem *gHLT = AliHLTPluginBase::GetInstance();
//   
//   TString calibInput;
//   if(calibInput.Length()>0) calibInput+=" ";
//   calibInput+="GLOBAL-esd-converter TPC-clusters TPC-globalmerger";
//     
//   TString calibOption = option;
//   if     (calibOption.CompareTo("task")==0)          AliHLTConfiguration calibtimeconf("TPCcalib",      "TPCCalibration",   calibInput.Data(), "");
//   else if(calibOption.CompareTo("calibtime")==0)     AliHLTConfiguration calibtimeconf("calibTime",     "TPCCalibTime",     calibInput.Data(), "");
//   else if(calibOption.CompareTo("calibtimegain")==0) AliHLTConfiguration calibtimeconf("calibTimeGain", "TPCCalibTimeGain", calibInput.Data(), "");
//   else 
//     {
//       cerr << "\nPlease specify an option for the calibration component you want to run, e.g. aliroot -b -q testTPCCalibration.C'(\"raw.root\",\"calibtime\")' \n"<< endl;
//       return;
//     }
// 
// 
//   /////////////////////////////////////////////////////////////////////////
//   //
//   // Init and run the reconstruction
//   // All but HLT reconstruction is switched off
//  
//   // connect to the GRID if we use a file or OCDB from the GRID
//   TString struri  = cdbURI;
//   TString strfile = filename;
//   if(struri.BeginsWith("raw://") ||  strfile.Contains("://") && !strfile.Contains("local://")){
//      TGrid::Connect("alien");
//   }
//  
//   // Set the CDB storage location
//   AliCDBManager *man = AliCDBManager::Instance();
//   man->SetDefaultStorage(cdbURI);
//   
//   AliReconstruction rec;
//   rec.SetInput(filename);
//   rec.SetRunVertexFinder(kFALSE);
//   //rec.SetWriteESDfriend(kTRUE);
//  
//   if(minEvent>=0 || maxEvent>minEvent){
//      if(minEvent<0) minEvent=0;
//      if(maxEvent<minEvent) maxEvent=minEvent;
//      rec.SetEventRange(minEvent,maxEvent);
//   } 
//   
//   if     (calibOption.CompareTo("task")==0)          rec.SetOption("HLT", "libAliHLTUtil.so libAliHLTTPC.so libAliHLTGlobal.so libAliHLTTPCCalibration.so loglevel=0x7c chains=TPCcalib"); 
//   if     (calibOption.CompareTo("calibtime")==0)     rec.SetOption("HLT", "libAliHLTUtil.so libAliHLTTPC.so libAliHLTGlobal.so libAliHLTTPCCalibration.so loglevel=0x7c chains=calibTime"); 
//   else if(calibOption.CompareTo("calibtimegain")==0) rec.SetOption("HLT", "libAliHLTUtil.so libAliHLTTPC.so libAliHLTGlobal.so libAliHLTTPCCalibration.so loglevel=0x7c chains=calibTimeGain");
//   
//   rec.SetRunPlaneEff(kFALSE);    
//   // switch off cleanESD
//   rec.SetCleanESD(kFALSE);
//   
//   rec.SetRunReconstruction("HLT");
//   rec.SetSpecificStorage("GRP/GRP/Data", Form("local://%s",gSystem->pwd()));
// 
// //   rec.SetRunLocalReconstruction("HLT");
// //   rec.SetRunTracking("");
// //   rec.SetFillESD("");
//  
//   rec.SetRunQA(":");
// 
//   AliLog::Flush();
//   rec.Run();
// }
// 
// void testTPCCalibration(const char *filename, int minEvent=-1, int maxEvent=-1, const char* option="task"){
//  
//   testTPCCalibration(filename, "raw://", minEvent, maxEvent, option);
// }
// 
// void testTPCCalibration(){
// 
//   cout << "testTPCCalibration: Run AliRoot reconstruction and TPC drift velocity calibration locally" << endl;
//   cout << " Usage: aliroot -b -q -l testTPCCalibration.C'(\"file\", \"cdb\", minEvent, maxEvent, option)'" << endl;
//   cout << "" << endl;
//   cout << " Examples:" << endl;
//   cout << "     testTPCCalibration.C'(\"raw://run12345\", minEvent, MaxEvent,\"calibtime\")'" << endl;
//   cout << "     testTPCCalibration.C'(\"raw.root\", \"local://$PWD\", minEvent, MaxEvent, \"calibtime\")'" << endl;
//   cout << "" << endl;
//   cout << " Defaults" << endl;
//   cout << "     cdb=\"raw://\"  -> takes OCDB from GRID" << endl;
//   cout << "     minEvent=-1   -> no lower event selection" << endl;
//   cout << "     maxEvent=-1   -> no upper event selection" << endl;
//   cout << "     option=\"task\" -> loads the analysis task that calls both the drift velocity and the gain calibration" << endl;
// }





void testTPCCalibration(const char* input="./", const char* option="task"){

  gSystem->Load("libANALYSIS");
  gSystem->Load("libTPCcalib");  
  
  if(!gSystem->AccessPathName("galice.root")){
    cerr << "Please delete file galice.root or run at a different place." << endl;
    return;
  }

  if (!input) {
    cerr << "Please specify input or run without arguments." << endl;
    return;
  }

  /////////////////////////////////////////////////////////////////////////
  //
  // init the HLT system in order to define the analysis chain below
  //
 
  AliHLTSystem *gHLT=AliHLTPluginBase::GetInstance();
  
  TString calibInput;
  if(calibInput.Length()>0) calibInput+=" ";
  calibInput+="GLOBAL-esd-converter TPC-globalmerger TPC-clusters";
    
  TString calibOption = option;
  if     (calibOption.CompareTo("task")==0)          AliHLTConfiguration calibtimeconf("TPCcalib",      "TPCCalibration",   calibInput.Data(), "");
  else if(calibOption.CompareTo("calibtime")==0)     AliHLTConfiguration calibtimeconf("calibTime",     "TPCCalibTime",     calibInput.Data(), "");
  else if(calibOption.CompareTo("calibtimegain")==0) AliHLTConfiguration calibtimeconf("calibTimeGain", "TPCCalibTimeGain", calibInput.Data(), "");
  else 
    {
      cerr << "\nPlease specify an option for the calibration component you want to run, e.g. aliroot -b -q testTPCCalibration.C'(\"raw.root\",\"calibtime\")' \n"<< endl;
      return;
    }


  /////////////////////////////////////////////////////////////////////////
  //
  // Init and run the reconstruction
  // All but HLT reconstruction is switched off
  //
  AliReconstruction rec;
  rec.SetInput(input);
  rec.SetRunVertexFinder(kFALSE);
  rec.SetRunLocalReconstruction("HLT");
  rec.SetRunTracking("");
  rec.SetLoadAlignFromCDB(0);
  rec.SetRunQA(":");
  rec.SetDefaultStorage("local://$ALICE_ROOT/OCDB");
  rec.SetSpecificStorage("GRP/GRP/Data", Form("local://%s",gSystem->pwd()));


  // NOTE: FillESD is a step in the AliReconstruction sequence and has
  // nothing to do with the fact that this macro writes ESD output
  // HLT processes the HLTOUT during FillESD and extracts data which
  // has already been prepared. This step is currently not necessary for
  // this macro
  rec.SetFillESD("");
  if     (calibOption.CompareTo("task")==0)          rec.SetOption("HLT", "libAliHLTUtil.so libAliHLTTPC.so libAliHLTGlobal.so libAliHLTTPCCalibration.so loglevel=0x7c chains=TPCcalib"); 
  if     (calibOption.CompareTo("calibtime")==0)     rec.SetOption("HLT", "libAliHLTUtil.so libAliHLTTPC.so libAliHLTGlobal.so libAliHLTTPCCalibration.so loglevel=0x7c chains=calibTime"); 
  else if(calibOption.CompareTo("calibtimegain")==0) rec.SetOption("HLT", "libAliHLTUtil.so libAliHLTTPC.so libAliHLTGlobal.so libAliHLTTPCCalibration.so loglevel=0x7c chains=calibTimeGain");

  rec.Run();
}
 testTPCCalibration.C:1
 testTPCCalibration.C:2
 testTPCCalibration.C:3
 testTPCCalibration.C:4
 testTPCCalibration.C:5
 testTPCCalibration.C:6
 testTPCCalibration.C:7
 testTPCCalibration.C:8
 testTPCCalibration.C:9
 testTPCCalibration.C:10
 testTPCCalibration.C:11
 testTPCCalibration.C:12
 testTPCCalibration.C:13
 testTPCCalibration.C:14
 testTPCCalibration.C:15
 testTPCCalibration.C:16
 testTPCCalibration.C:17
 testTPCCalibration.C:18
 testTPCCalibration.C:19
 testTPCCalibration.C:20
 testTPCCalibration.C:21
 testTPCCalibration.C:22
 testTPCCalibration.C:23
 testTPCCalibration.C:24
 testTPCCalibration.C:25
 testTPCCalibration.C:26
 testTPCCalibration.C:27
 testTPCCalibration.C:28
 testTPCCalibration.C:29
 testTPCCalibration.C:30
 testTPCCalibration.C:31
 testTPCCalibration.C:32
 testTPCCalibration.C:33
 testTPCCalibration.C:34
 testTPCCalibration.C:35
 testTPCCalibration.C:36
 testTPCCalibration.C:37
 testTPCCalibration.C:38
 testTPCCalibration.C:39
 testTPCCalibration.C:40
 testTPCCalibration.C:41
 testTPCCalibration.C:42
 testTPCCalibration.C:43
 testTPCCalibration.C:44
 testTPCCalibration.C:45
 testTPCCalibration.C:46
 testTPCCalibration.C:47
 testTPCCalibration.C:48
 testTPCCalibration.C:49
 testTPCCalibration.C:50
 testTPCCalibration.C:51
 testTPCCalibration.C:52
 testTPCCalibration.C:53
 testTPCCalibration.C:54
 testTPCCalibration.C:55
 testTPCCalibration.C:56
 testTPCCalibration.C:57
 testTPCCalibration.C:58
 testTPCCalibration.C:59
 testTPCCalibration.C:60
 testTPCCalibration.C:61
 testTPCCalibration.C:62
 testTPCCalibration.C:63
 testTPCCalibration.C:64
 testTPCCalibration.C:65
 testTPCCalibration.C:66
 testTPCCalibration.C:67
 testTPCCalibration.C:68
 testTPCCalibration.C:69
 testTPCCalibration.C:70
 testTPCCalibration.C:71
 testTPCCalibration.C:72
 testTPCCalibration.C:73
 testTPCCalibration.C:74
 testTPCCalibration.C:75
 testTPCCalibration.C:76
 testTPCCalibration.C:77
 testTPCCalibration.C:78
 testTPCCalibration.C:79
 testTPCCalibration.C:80
 testTPCCalibration.C:81
 testTPCCalibration.C:82
 testTPCCalibration.C:83
 testTPCCalibration.C:84
 testTPCCalibration.C:85
 testTPCCalibration.C:86
 testTPCCalibration.C:87
 testTPCCalibration.C:88
 testTPCCalibration.C:89
 testTPCCalibration.C:90
 testTPCCalibration.C:91
 testTPCCalibration.C:92
 testTPCCalibration.C:93
 testTPCCalibration.C:94
 testTPCCalibration.C:95
 testTPCCalibration.C:96
 testTPCCalibration.C:97
 testTPCCalibration.C:98
 testTPCCalibration.C:99
 testTPCCalibration.C:100
 testTPCCalibration.C:101
 testTPCCalibration.C:102
 testTPCCalibration.C:103
 testTPCCalibration.C:104
 testTPCCalibration.C:105
 testTPCCalibration.C:106
 testTPCCalibration.C:107
 testTPCCalibration.C:108
 testTPCCalibration.C:109
 testTPCCalibration.C:110
 testTPCCalibration.C:111
 testTPCCalibration.C:112
 testTPCCalibration.C:113
 testTPCCalibration.C:114
 testTPCCalibration.C:115
 testTPCCalibration.C:116
 testTPCCalibration.C:117
 testTPCCalibration.C:118
 testTPCCalibration.C:119
 testTPCCalibration.C:120
 testTPCCalibration.C:121
 testTPCCalibration.C:122
 testTPCCalibration.C:123
 testTPCCalibration.C:124
 testTPCCalibration.C:125
 testTPCCalibration.C:126
 testTPCCalibration.C:127
 testTPCCalibration.C:128
 testTPCCalibration.C:129
 testTPCCalibration.C:130
 testTPCCalibration.C:131
 testTPCCalibration.C:132
 testTPCCalibration.C:133
 testTPCCalibration.C:134
 testTPCCalibration.C:135
 testTPCCalibration.C:136
 testTPCCalibration.C:137
 testTPCCalibration.C:138
 testTPCCalibration.C:139
 testTPCCalibration.C:140
 testTPCCalibration.C:141
 testTPCCalibration.C:142
 testTPCCalibration.C:143
 testTPCCalibration.C:144
 testTPCCalibration.C:145
 testTPCCalibration.C:146
 testTPCCalibration.C:147
 testTPCCalibration.C:148
 testTPCCalibration.C:149
 testTPCCalibration.C:150
 testTPCCalibration.C:151
 testTPCCalibration.C:152
 testTPCCalibration.C:153
 testTPCCalibration.C:154
 testTPCCalibration.C:155
 testTPCCalibration.C:156
 testTPCCalibration.C:157
 testTPCCalibration.C:158
 testTPCCalibration.C:159
 testTPCCalibration.C:160
 testTPCCalibration.C:161
 testTPCCalibration.C:162
 testTPCCalibration.C:163
 testTPCCalibration.C:164
 testTPCCalibration.C:165
 testTPCCalibration.C:166
 testTPCCalibration.C:167
 testTPCCalibration.C:168
 testTPCCalibration.C:169
 testTPCCalibration.C:170
 testTPCCalibration.C:171
 testTPCCalibration.C:172
 testTPCCalibration.C:173
 testTPCCalibration.C:174
 testTPCCalibration.C:175
 testTPCCalibration.C:176
 testTPCCalibration.C:177
 testTPCCalibration.C:178
 testTPCCalibration.C:179
 testTPCCalibration.C:180
 testTPCCalibration.C:181
 testTPCCalibration.C:182
 testTPCCalibration.C:183
 testTPCCalibration.C:184
 testTPCCalibration.C:185
 testTPCCalibration.C:186
 testTPCCalibration.C:187
 testTPCCalibration.C:188
 testTPCCalibration.C:189
 testTPCCalibration.C:190
 testTPCCalibration.C:191
 testTPCCalibration.C:192
 testTPCCalibration.C:193
 testTPCCalibration.C:194
 testTPCCalibration.C:195
 testTPCCalibration.C:196
 testTPCCalibration.C:197
 testTPCCalibration.C:198
 testTPCCalibration.C:199
 testTPCCalibration.C:200
 testTPCCalibration.C:201
 testTPCCalibration.C:202
 testTPCCalibration.C:203
 testTPCCalibration.C:204
 testTPCCalibration.C:205
 testTPCCalibration.C:206
 testTPCCalibration.C:207
 testTPCCalibration.C:208
 testTPCCalibration.C:209
 testTPCCalibration.C:210
 testTPCCalibration.C:211
 testTPCCalibration.C:212
 testTPCCalibration.C:213
 testTPCCalibration.C:214
 testTPCCalibration.C:215