- Return from Gstpar if material is not used.

[u/mrichter/AliRoot.git] / TRD / AliTRDReconstructor.cxx

diff --git a/TRD/AliTRDReconstructor.cxx b/TRD/AliTRDReconstructor.cxx
index 85d2dcee0978f096d4db33d80edd9655e3dde36d..d48f229f9b1e7f20db09fb73eebcf2ff59e224bf 100644 (file)
--- a/TRD/AliTRDReconstructor.cxx
+++ b/TRD/AliTRDReconstructor.cxx
@@ -27,6 +27,7 @@
 #include "AliTRDparameter.h"
 #include "AliTRDclusterizerV1.h"
 #include "AliTRDtracker.h"
+#include "AliTRDpidESD.h"
 #include <TFile.h>
 
 
@@ -58,6 +59,8 @@ void AliTRDReconstructor::Reconstruct(AliRunLoader* runLoader) const
     clusterer.MakeClusters();
     clusterer.WriteClusters(-1);
   }
+
+  loader->UnloadRecPoints();
 }
 
 //_____________________________________________________________________________
@@ -71,8 +74,17 @@ AliTracker* AliTRDReconstructor::CreateTracker(AliRunLoader* runLoader) const
 
 //_____________________________________________________________________________
 void AliTRDReconstructor::FillESD(AliRunLoader* /*runLoader*/, 
-                                 AliESD* /*esd*/) const
+                                 AliESD* esd) const
 {
+// make PID
+
+  Double_t parTRD[] = {
+    187., // Min. Ionizing Particle signal.  Check it !!!
+    0.23, // relative resolution             Check it !!!
+    10.   // PID range (in sigmas)
+  };
+  AliTRDpidESD trdPID(parTRD);
+  trdPID.MakePID(esd);
 }