Reconstruction of raw data (T.Kuhr)

diff --git a/TPC/AliTPCReconstructor.cxx b/TPC/AliTPCReconstructor.cxx
index 8e284e2d290a20fa6d0bacca4092192b526c2f9f..a491a174359bb40e9549e515aa8f1c787cd872de 100644 (file)
--- a/TPC/AliTPCReconstructor.cxx
+++ b/TPC/AliTPCReconstructor.cxx
@@ -25,6 +25,7 @@
 #include "AliTPCReconstructor.h"
 #include "AliRunLoader.h"
 #include "AliRun.h"
 #include "AliTPCReconstructor.h"
 #include "AliRunLoader.h"
 #include "AliRun.h"

+#include "AliRawReader.h"

 #include "AliTPCclustererMI.h"
 #include "AliTPCtrackerMI.h"
 #include "AliTPCpidESD.h"
 #include "AliTPCclustererMI.h"
 #include "AliTPCtrackerMI.h"
 #include "AliTPCpidESD.h"


@@ -77,6 +78,42 @@ void AliTPCReconstructor::Reconstruct(AliRunLoader* runLoader) const
   loader->UnloadDigits();
 }
 
   loader->UnloadDigits();
 }
 

+//_____________________________________________________________________________
+void AliTPCReconstructor::Reconstruct(AliRunLoader* runLoader,
+                                     AliRawReader* rawReader) const
+{
+// reconstruct clusters from raw data
+
+  AliLoader* loader = runLoader->GetLoader("TPCLoader");
+  if (!loader) {
+    Error("Reconstruct", "TPC loader not found");
+    return;
+  }
+  loader->LoadRecPoints("recreate");
+
+  AliTPCParam* param = GetTPCParam(runLoader);
+  if (!param) return;
+  AliTPCclustererMI clusterer(param);
+
+  Int_t iEvent = 0;
+  while (rawReader->NextEvent()) {
+    runLoader->GetEvent(iEvent++);
+
+    TTree* treeClusters = loader->TreeR();
+    if (!treeClusters) {
+      loader->MakeTree("R");
+      treeClusters = loader->TreeR();
+    }
+
+    clusterer.SetOutput(treeClusters);
+    clusterer.Digits2Clusters(rawReader);
+         
+    loader->WriteRecPoints("OVERWRITE");
+  }
+
+  loader->UnloadRecPoints();
+}
+

 //_____________________________________________________________________________
 AliTracker* AliTPCReconstructor::CreateTracker(AliRunLoader* runLoader) const
 {
 //_____________________________________________________________________________
 AliTracker* AliTPCReconstructor::CreateTracker(AliRunLoader* runLoader) const
 {

diff --git a/TPC/AliTPCReconstructor.h b/TPC/AliTPCReconstructor.h
index babcf6aa9f2b4e70e3bad61002a82fe14fb5be81..52baafc331b0d3447b6aad62710b2db31c126dd1 100644 (file)
--- a/TPC/AliTPCReconstructor.h
+++ b/TPC/AliTPCReconstructor.h
@@ -16,6 +16,8 @@ public:
   virtual ~AliTPCReconstructor() {};
 
   virtual void         Reconstruct(AliRunLoader* runLoader) const;
   virtual ~AliTPCReconstructor() {};
 
   virtual void         Reconstruct(AliRunLoader* runLoader) const;

+  virtual void         Reconstruct(AliRunLoader* runLoader,
+                                  AliRawReader* rawReader) const;

   virtual AliTracker*  CreateTracker(AliRunLoader* runLoader) const;
   virtual void         FillESD(AliRunLoader* runLoader, AliESD* esd) const;
 
   virtual AliTracker*  CreateTracker(AliRunLoader* runLoader) const;
   virtual void         FillESD(AliRunLoader* runLoader, AliESD* esd) const;