ROOT logo
// $Id$
/*
 * Example macro to run the ITS tracker with the TPC reconstruction.
 * The reconstruction is done from the TPC and ITS raw data.
 *
 * Usage:
 * <pre>
 *   aliroot -b -q rec-tpc-its.C | tee rec-tpc-its.log
 *   aliroot -b -q rec-tpc-its.C'("./","spd")' | tee rec-tpc-its.log
 * </pre>
 *
 * The macro asumes raw data to be available in the rawx folders, either
 * simulated or real data. A different input can be specified as parameter
 * <pre>
 *   aliroot -b -q rec-tpc-its.C'("input.root")'
 * </pre>
 * 
 * In the second parameter you can change the clusterfinders to run only
 * SPD, SDD, SSD or All:
 *    - SPD runs only Silicon Pixels
 *    - SDD runs only Silicon Drift
 *    - SSD runs only Silicon Stips
 *    - All will run the full ITS. This is default
 *
 * In the first section, an analysis chain is defined. The scale of the
 * chain can be defined by choosing the range of sectors and partitions.
 *
 * The reconstruction is steered by the AliReconstruction object in the
 * usual way.
 *
 * @ingroup alihlt_tpc
 * @author Gaute.Ovrebekk@ift.uib.no
 */
void rec_tpc_its(const char* input="./", char* opt="All")
{
  
  if(!gSystem->AccessPathName("galice.root")){
    cerr << "please delete the galice.root or run at 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();
 
  ///////////////////////////////////////////////////////////////////////////////////////////////////
  //
  // Setting up which output to give
  //
  TString option="libAliHLTUtil.so libAliHLTRCU.so libAliHLTTPC.so libAliHLTITS.so libAliHLTGlobal.so loglevel=0x7c chains=";
  Bool_t runspd=kFALSE, runsdd=kFALSE, runssd=kFALSE;
  TString allArgs=opt;
  TString argument;
  TObjArray* pTokens=allArgs.Tokenize(" ");
  if (pTokens) {
    for (int i=0; i<pTokens->GetEntries(); i++) {
      argument=((TObjString*)pTokens->At(i))->GetString();
      if (argument.IsNull()) continue;
 
      if (argument.CompareTo("spd", TString::kIgnoreCase)==0) {
	runspd=kTRUE;
	continue;
      } 
      if (argument.CompareTo("sdd", TString::kIgnoreCase)==0) {
	runsdd=kTRUE;
	continue;
      }
      if (argument.CompareTo("ssd",TString::kIgnoreCase)==0) {
	runssd=kTRUE;
	continue;
      }
      if (argument.CompareTo("all",TString::kIgnoreCase)==0) {
	runspd=kTRUE;
	runsdd=kTRUE;
	runssd=kTRUE;
	continue;
      }
      else {
	cerr << "Unknown argument" << endl;
	break;
      }
    }
  }

  ///////////////////////////////////////////////////////////////////////////////////////////////////
  //
  // define the analysis chain to be run
  //
  int minddl=0;          //min ddl number for SPD
  int maxddl=19;         //max ddl number for SPD
  int spec=0x1;          //spec for ddl's
  int ddlno=0;
  TString cfout="";

  if(runspd){
    for(ddlno=minddl;ddlno<=maxddl;ddlno++){  
      TString arg, publisher, cf;
      
      arg.Form("-minid %d -datatype 'DDL_RAW ' 'ISPD ' -dataspec 0x%08x -verbose",ddlno, spec);
      publisher.Form("DP_%d", ddlno);
      AliHLTConfiguration pubconf(publisher.Data(), "AliRawReaderPublisher", NULL , arg.Data());
      
      cf.Form("CF_%d",ddlno);
      AliHLTConfiguration cfconf(cf.Data(), "ITSClusterFinderSPD", publisher.Data(), "");
      
      if (cfout.Length()>0) cfout+=" ";
      cfout+=cf;
      
      spec=spec<<1;
    }
  }
  
  if(runsdd){
    minddl=256;        //min ddl number for SDD    
    maxddl=279;        //max ddl number for SDD
    spec=0x1;          //spec for ddl's
    
    for(ddlno=minddl;ddlno<=maxddl;ddlno++){  
      TString arg, publisher, cf;
      
      arg.Form("-minid %d -datatype 'DDL_RAW ' 'ISDD ' -dataspec 0x%08x -verbose",ddlno, spec); 
      publisher.Form("DP_%d", ddlno);
      AliHLTConfiguration pubconf(publisher.Data(), "AliRawReaderPublisher", NULL , arg.Data());
      
      cf.Form("CF_%d",ddlno);
      AliHLTConfiguration cfconf(cf.Data(), "ITSClusterFinderSDD", publisher.Data(), "");
      
      if (cfout.Length()>0) cfout+=" ";
      cfout+=cf;
      
      spec=spec<<1;
    }
  }

  if(runssd){
    minddl=512;      //min ddl number for SSD     
    maxddl=527;      //max ddl number for SSD
    spec=0x1;        //spec for ddl's
    
    for(ddlno=minddl;ddlno<=maxddl;ddlno++){  
      TString arg, publisher, cf;
      
      arg.Form("-minid %d -datatype 'DDL_RAW ' 'ISSD ' -dataspec 0x%08x -verbose",ddlno, spec);
      publisher.Form("DP_%d", ddlno);
      AliHLTConfiguration pubconf(publisher.Data(), "AliRawReaderPublisher", NULL , arg.Data());
      
      cf.Form("CF_%d",ddlno);
      AliHLTConfiguration cfconf(cf.Data(), "ITSClusterFinderSSD", publisher.Data(), "");
      
      if (cfout.Length()>0) cfout+=" ";
      cfout+=cf;
      
      spec=spec<<1;
    }
  }

  TString ITSinput = "TPC-globalmerger ";
  ITSinput += cfout;
  AliHLTConfiguration itstrackerconf("itstracker","ITSTracker",ITSinput.Data(),""); 
  //option+="itstracker";
  
  AliHLTConfiguration globalConverter("globalConverter", "GlobalEsdConverter"   , "TPC-globalmerger itstracker", "");
  AliHLTConfiguration sink("esdfile", "EsdCollector"   , "globalConverter", "-directory hlt-tpc-esd");
  option+="esdfile";

  ///////////////////////////////////////////////////////////////////////////////////////////////////
  //
  // Init and run the reconstruction
  // All but HLT reconstructio is switched off
  //
  AliReconstruction rec;
  rec.SetInput(input);
  rec.SetRunVertexFinder(kFALSE);
  rec.SetRunReconstruction("HLT");
  rec.SetLoadAlignFromCDB(0);
  rec.SetRunQA(":");
  rec.SetDefaultStorage("local://$ALICE_ROOT/OCDB");
  rec.SetSpecificStorage("GRP/GRP/Data",Form("local://%s",gSystem->pwd())); 
  rec.SetOption("HLT", option);
  rec.Run();
}
 rec-tpc-its.C:1
 rec-tpc-its.C:2
 rec-tpc-its.C:3
 rec-tpc-its.C:4
 rec-tpc-its.C:5
 rec-tpc-its.C:6
 rec-tpc-its.C:7
 rec-tpc-its.C:8
 rec-tpc-its.C:9
 rec-tpc-its.C:10
 rec-tpc-its.C:11
 rec-tpc-its.C:12
 rec-tpc-its.C:13
 rec-tpc-its.C:14
 rec-tpc-its.C:15
 rec-tpc-its.C:16
 rec-tpc-its.C:17
 rec-tpc-its.C:18
 rec-tpc-its.C:19
 rec-tpc-its.C:20
 rec-tpc-its.C:21
 rec-tpc-its.C:22
 rec-tpc-its.C:23
 rec-tpc-its.C:24
 rec-tpc-its.C:25
 rec-tpc-its.C:26
 rec-tpc-its.C:27
 rec-tpc-its.C:28
 rec-tpc-its.C:29
 rec-tpc-its.C:30
 rec-tpc-its.C:31
 rec-tpc-its.C:32
 rec-tpc-its.C:33
 rec-tpc-its.C:34
 rec-tpc-its.C:35
 rec-tpc-its.C:36
 rec-tpc-its.C:37
 rec-tpc-its.C:38
 rec-tpc-its.C:39
 rec-tpc-its.C:40
 rec-tpc-its.C:41
 rec-tpc-its.C:42
 rec-tpc-its.C:43
 rec-tpc-its.C:44
 rec-tpc-its.C:45
 rec-tpc-its.C:46
 rec-tpc-its.C:47
 rec-tpc-its.C:48
 rec-tpc-its.C:49
 rec-tpc-its.C:50
 rec-tpc-its.C:51
 rec-tpc-its.C:52
 rec-tpc-its.C:53
 rec-tpc-its.C:54
 rec-tpc-its.C:55
 rec-tpc-its.C:56
 rec-tpc-its.C:57
 rec-tpc-its.C:58
 rec-tpc-its.C:59
 rec-tpc-its.C:60
 rec-tpc-its.C:61
 rec-tpc-its.C:62
 rec-tpc-its.C:63
 rec-tpc-its.C:64
 rec-tpc-its.C:65
 rec-tpc-its.C:66
 rec-tpc-its.C:67
 rec-tpc-its.C:68
 rec-tpc-its.C:69
 rec-tpc-its.C:70
 rec-tpc-its.C:71
 rec-tpc-its.C:72
 rec-tpc-its.C:73
 rec-tpc-its.C:74
 rec-tpc-its.C:75
 rec-tpc-its.C:76
 rec-tpc-its.C:77
 rec-tpc-its.C:78
 rec-tpc-its.C:79
 rec-tpc-its.C:80
 rec-tpc-its.C:81
 rec-tpc-its.C:82
 rec-tpc-its.C:83
 rec-tpc-its.C:84
 rec-tpc-its.C:85
 rec-tpc-its.C:86
 rec-tpc-its.C:87
 rec-tpc-its.C:88
 rec-tpc-its.C:89
 rec-tpc-its.C:90
 rec-tpc-its.C:91
 rec-tpc-its.C:92
 rec-tpc-its.C:93
 rec-tpc-its.C:94
 rec-tpc-its.C:95
 rec-tpc-its.C:96
 rec-tpc-its.C:97
 rec-tpc-its.C:98
 rec-tpc-its.C:99
 rec-tpc-its.C:100
 rec-tpc-its.C:101
 rec-tpc-its.C:102
 rec-tpc-its.C:103
 rec-tpc-its.C:104
 rec-tpc-its.C:105
 rec-tpc-its.C:106
 rec-tpc-its.C:107
 rec-tpc-its.C:108
 rec-tpc-its.C:109
 rec-tpc-its.C:110
 rec-tpc-its.C:111
 rec-tpc-its.C:112
 rec-tpc-its.C:113
 rec-tpc-its.C:114
 rec-tpc-its.C:115
 rec-tpc-its.C:116
 rec-tpc-its.C:117
 rec-tpc-its.C:118
 rec-tpc-its.C:119
 rec-tpc-its.C:120
 rec-tpc-its.C:121
 rec-tpc-its.C:122
 rec-tpc-its.C:123
 rec-tpc-its.C:124
 rec-tpc-its.C:125
 rec-tpc-its.C:126
 rec-tpc-its.C:127
 rec-tpc-its.C:128
 rec-tpc-its.C:129
 rec-tpc-its.C:130
 rec-tpc-its.C:131
 rec-tpc-its.C:132
 rec-tpc-its.C:133
 rec-tpc-its.C:134
 rec-tpc-its.C:135
 rec-tpc-its.C:136
 rec-tpc-its.C:137
 rec-tpc-its.C:138
 rec-tpc-its.C:139
 rec-tpc-its.C:140
 rec-tpc-its.C:141
 rec-tpc-its.C:142
 rec-tpc-its.C:143
 rec-tpc-its.C:144
 rec-tpc-its.C:145
 rec-tpc-its.C:146
 rec-tpc-its.C:147
 rec-tpc-its.C:148
 rec-tpc-its.C:149
 rec-tpc-its.C:150
 rec-tpc-its.C:151
 rec-tpc-its.C:152
 rec-tpc-its.C:153
 rec-tpc-its.C:154
 rec-tpc-its.C:155
 rec-tpc-its.C:156
 rec-tpc-its.C:157
 rec-tpc-its.C:158
 rec-tpc-its.C:159
 rec-tpc-its.C:160
 rec-tpc-its.C:161
 rec-tpc-its.C:162
 rec-tpc-its.C:163
 rec-tpc-its.C:164
 rec-tpc-its.C:165
 rec-tpc-its.C:166
 rec-tpc-its.C:167
 rec-tpc-its.C:168
 rec-tpc-its.C:169
 rec-tpc-its.C:170
 rec-tpc-its.C:171
 rec-tpc-its.C:172
 rec-tpc-its.C:173
 rec-tpc-its.C:174
 rec-tpc-its.C:175
 rec-tpc-its.C:176
 rec-tpc-its.C:177
 rec-tpc-its.C:178
 rec-tpc-its.C:179
 rec-tpc-its.C:180
 rec-tpc-its.C:181
 rec-tpc-its.C:182
 rec-tpc-its.C:183
 rec-tpc-its.C:184
 rec-tpc-its.C:185
 rec-tpc-its.C:186
 rec-tpc-its.C:187
 rec-tpc-its.C:188