Adding libEMCALjet

[u/mrichter/AliRoot.git] / TOF / AliTOFReconstructor.cxx

diff --git a/TOF/AliTOFReconstructor.cxx b/TOF/AliTOFReconstructor.cxx
index 6ff0e34a49549a1ccd9051a4dff9b77118dd4606..744180f32c760bb4d9ed1ea12baf8433624f4125 100644 (file)
--- a/TOF/AliTOFReconstructor.cxx
+++ b/TOF/AliTOFReconstructor.cxx
@@ -22,20 +22,63 @@
 ///////////////////////////////////////////////////////////////////////////////
 
 
-#include "AliTOFReconstructor.h"
-#include "AliRunLoader.h"
+#include <TFile.h>
+
+#include "AliLog.h"
 #include "AliRun.h"
-#include "AliTOF.h"
-#include "AliTOFtracker.h"
+#include "AliRunLoader.h"
 
+#include "AliTOFGeometry.h"
+#include "AliTOFGeometryV4.h"
+#include "AliTOFGeometryV5.h"
+#include "AliTOFtracker.h"
+#include "AliTOFtrackerMI.h"
+#include "AliTOFClusterFinder.h"
+#include "AliTOFReconstructor.h"
 
 ClassImp(AliTOFReconstructor)
 
+//_____________________________________________________________________________
+  void AliTOFReconstructor::Reconstruct(AliRunLoader* runLoader) const
+{
+// reconstruct clusters from digits
+
+  AliTOFClusterFinder *tofClus = new AliTOFClusterFinder(runLoader);
+  tofClus->Load();
+  for (Int_t iEvent = 0; iEvent < runLoader->GetNumberOfEvents(); iEvent++)
+    {
+      tofClus->Digits2RecPoints(iEvent);
+    }
+  tofClus->UnLoad();
+
+}
 
 //_____________________________________________________________________________
-void AliTOFReconstructor::Reconstruct(AliRunLoader* /*runLoader*/) const
+void AliTOFReconstructor::Reconstruct(AliRunLoader* runLoader,
+                                      AliRawReader *rawReader) const
 {
-// nothing to be done
+// reconstruct clusters from Raw Data
+
+  AliTOFClusterFinder tofClus(runLoader);
+  tofClus.LoadClusters();
+  Int_t iEvent = 0;
+  while (rawReader->NextEvent()) {
+    tofClus.Digits2RecPoints(iEvent,rawReader);
+    //tofClus.Raw2Digits(iEvent,rawReader); // temporary solution
+    iEvent++;
+  }
+  tofClus.UnLoadClusters();
+
+}
+
+//_____________________________________________________________________________
+void AliTOFReconstructor::Reconstruct(AliRawReader *rawReader,
+                                      TTree *clustersTree) const
+{
+// reconstruct clusters from Raw Data
+
+  AliTOFClusterFinder tofClus;
+  tofClus.Digits2RecPoints(rawReader, clustersTree);
 
 }
 
@@ -46,7 +89,12 @@ AliTracker* AliTOFReconstructor::CreateTracker(AliRunLoader* runLoader) const
 
   AliTOFGeometry* geom = GetTOFGeometry(runLoader);
   if (!geom) return NULL;
-  Double_t parPID[] = {130., 5.};
+  //  Double_t parPID[] = {130., 5.};
+  Double_t parPID[] = {80., 5.};
+  TString selectedTracker = GetOption();
+  // use MI tracker if selected
+  if (selectedTracker.Contains("MI")) return new AliTOFtrackerMI(geom,parPID);
+
   return new AliTOFtracker(geom, parPID);
 }
 
@@ -58,25 +106,30 @@ void AliTOFReconstructor::FillESD(AliRunLoader* /*runLoader*/,
 
 }
 
-
 //_____________________________________________________________________________
 AliTOFGeometry* AliTOFReconstructor::GetTOFGeometry(AliRunLoader* runLoader) const
 {
 // get the TOF parameters
 
-  if (!runLoader->GetAliRun()) runLoader->LoadgAlice();
-  if (!runLoader->GetAliRun()) {
-    Error("GetTOFGeometry", "couldn't get AliRun object");
-    return NULL;
+  AliTOFGeometry *tofGeom;
+
+  runLoader->CdGAFile();
+  TDirectory *savedir=gDirectory; 
+  TFile *in=(TFile*)gFile;  
+  if (!in->IsOpen()) {
+    AliWarning("Geometry file is not open default  TOF geometry will be used");
+    tofGeom = new AliTOFGeometry();
   }
-  AliTOF* tof = (AliTOF*) runLoader->GetAliRun()->GetDetector("TOF");
-  if (!tof) {
-    Error("GetTOFGeometry", "couldn't get TOF detector");
-    return NULL;
+  else {
+    in->cd();  
+    tofGeom = (AliTOFGeometry*) in->Get("TOFgeometry");
   }
-  if (!tof->GetGeometry()) {
-    Error("GetTOFGeometry", "no TOF geometry available");
+
+  savedir->cd();  
+
+  if (!tofGeom) {
+    AliError("no TOF geometry available");
     return NULL;
   }
-  return tof->GetGeometry();
+  return tofGeom;
 }