load AddTaskEmcalJet.C

[u/mrichter/AliRoot.git] / ITS / AliITSDigits2RecPoints.C

diff --git a/ITS/AliITSDigits2RecPoints.C b/ITS/AliITSDigits2RecPoints.C
index 962cccc716279f85c2cfc6af0f69b190e273e13f..404b024f3025721d118663ae7de19e9c1fd6c5b0 100644 (file)
--- a/ITS/AliITSDigits2RecPoints.C
+++ b/ITS/AliITSDigits2RecPoints.C
@@ -1,15 +1,48 @@
-void AliITSDigits2RecPoints(Int_t evNumber1=0,Int_t evNumber2=0, const char *filename="galice.root"){
-    TStopwatch timer;
-
-    if(gAlice){
-       delete gAlice;
-       gAlice = 0;
-    } // end if gAlice
-    cout << "Creating reconstructed points from digits for the ITS..." << endl;
-    AliITSreconstruction *itsr = new AliITSreconstruction(filename);
-    timer.Start();
-    itsr->Init();
-    itsr->Exec(); 
-    timer.Stop(); timer.Print();    
-    delete itsr;
+#if !defined(__CINT__) || defined(__MAKECINT__)
+
+#include <Riostream.h>
+#include <TClassTable.h>
+#include <TDatime.h>
+#include <TGeoManager.h>
+#include <TROOT.h>
+#include <TStopwatch.h>
+#include "AliRun.h"
+#include "AliRunDigitizer.h"
+#include "AliITSDetTypeRec.h"
+#include "AliITS.h"
+#include "AliITSresponseSDD.h"
+#include "AliITSreconstruction.h"
+
+#endif
+#define DEBUG
+
+Int_t AliITSDigits2RecPoints(TString filename="galice.root",TString fileRP=""){
+  // Standard ITS Digits to RecPoints.
+
+  // Get geometry
+  TGeoManager::Import("geometry.root");
+
+  // Dynamically link some shared libs
+  if (gClassTable->GetID("AliRun") < 0) {
+    gROOT->ProcessLine(".x $(ALICE_ROOT)/macros/loadlibs.C");
+  }else if (gAlice){
+    delete AliRunLoader::Instance();
+    delete gAlice;
+    gAlice=0;
+  } // end if
+
+  TStopwatch timer;
+#ifdef DEBUG
+  cout << "Creating reconstructed points from digits for the ITS..." << endl;
+#endif
+  AliITSreconstruction *itsr = new AliITSreconstruction(filename);
+
+  timer.Start();
+  if(!(fileRP.IsNull()))itsr->SetOutputFile(fileRP);
+  itsr->Init();
+  itsr->Exec(); 
+  timer.Stop(); 
+  timer.Print();
+  delete itsr;
+  return 0;
 }