fDebugReconstruction and SetDebugReconstruction for AliPHOSReconstructioner. Debug...

diff --git a/PHOS/AliPHOSReconstructioner.cxx b/PHOS/AliPHOSReconstructioner.cxx
index 08ff29f14f40d447acfa954eb232846b09aa77d0..365fdfea6f7c6025ae679bffa60137503db962ec 100644 (file)
--- a/PHOS/AliPHOSReconstructioner.cxx
+++ b/PHOS/AliPHOSReconstructioner.cxx
@@ -29,6 +29,8 @@
 
 // --- Standard library ---
 
+#include <iomanip>
+
 // --- AliRoot header files ---
 
 #include "AliPHOSReconstructioner.h"
@@ -46,6 +48,7 @@ AliPHOSReconstructioner::AliPHOSReconstructioner(AliPHOSClusterizer * Clusterize
   fClusterizer        = Clusterizer ;
   fTrackSegmentMaker  = Tracker ;
   fPID                = Pid ; 
+  fDebugReconstruction = kFALSE ;
 } 
 
 
@@ -57,6 +60,7 @@ AliPHOSReconstructioner::AliPHOSReconstructioner(AliPHOSClusterizer * Clusterize
   fClusterizer        = Clusterizer ;
   fTrackSegmentMaker  = Tracker ;
   fPID                = Pid ; 
+  fDebugReconstruction = kFALSE ;
 } 
 
 //____________________________________________________________________________
@@ -68,10 +72,46 @@ AliPHOSReconstructioner::AliPHOSReconstructioner(AliPHOSClusterizer * Clusterize
   //                                                      Make the reconstructed particles
 
   Int_t index ; 
-  
-  cout << "Start making reconstructed points (clusterizing)" << endl;
+  // Digit Debuging
+
+
+  if  (fDebugReconstruction) 
+    {
+      cout << ">>>>>>>>>>>>>>>>>>>>>> DebugReconstruction  <<<<<<<<<<<<<<<<<<<<<<<<<<"  << endl ;
+      cout << "DebugReconstruction>>> Digit list entries is " <<    dl->GetEntries() << endl ;
+      AliPHOSDigit * Digit;
+      cout << "DebugReconstruction>>>    Vol Id " << 
+         " Energy (MeV) "             <<                         
+         " Index "                    << 
+         " Nprim "                     << 
+         " Prim1 "                      << 
+         " Prim2 "                      << 
+         " Prim3 "                      <<  endl;  
+
+      for (index = 0 ; index < dl->GetEntries() ; index++) {
+       Digit = (AliPHOSDigit * )  dl->At(index) ;
+       cout << "DebugReconstruction>>>  " << 
+         setw(8) <<Digit->GetId() << "  "  << 
+         setw(10) << 1000.*fClusterizer->Calibrate(Digit->GetAmp()) <<       "  "  <<                   
+         setw(6) <<  Digit->GetIndexInList() << "  "  << 
+         setw(5) <<  Digit->GetNprimary() <<"  "  << 
+         setw(5) <<  Digit->GetPrimary(1) <<"  "  << 
+         setw(5) <<  Digit->GetPrimary(2) <<"  "  << 
+         setw(5) <<  Digit->GetPrimary(3) << endl;      
+  }
+
+    }
+
+
+  if  (fDebugReconstruction)  cout << "DebugReconstruction > Start making reconstructed points (clusterizing)" << endl;
   fClusterizer->MakeClusters(dl, emccl, ppsdl);
 
+  cout << "Emccl list entries is " << emccl->GetEntries() << endl ;
+  cout << "Ppsdl list entries is " << emccl->GetEntries() << endl ;
+
+
+
+
   // mark the position of the RecPoints in the array
   AliPHOSEmcRecPoint * emcrp ; 
   for (index = 0 ; index < emccl->GetEntries() ; index++) {

diff --git a/PHOS/AliPHOSReconstructioner.h b/PHOS/AliPHOSReconstructioner.h
index 217687a6a524ea6db26be51f6840b9225c0cf905..ccb37a1dfa609965718cc47a4e35551067c065c6 100644 (file)
--- a/PHOS/AliPHOSReconstructioner.h
+++ b/PHOS/AliPHOSReconstructioner.h
@@ -39,9 +39,12 @@ public:
   void Make(TClonesArray * DL, RecPointsList * emccl, RecPointsList * ppsdl, 
            TrackSegmentsList * trsl, RecParticlesList * rpl) ; // does the job
 
+  void SetDebugReconstruction(Bool_t deb) { fDebugReconstruction = deb; }
+
 
 private:
   
+  Bool_t               fDebugReconstruction;      // For debuging of the Reconstruction procedure
   AliPHOSClusterizer * fClusterizer ;             // Method for clusterization 
   AliPHOSTrackSegmentMaker * fTrackSegmentMaker ; // Method for track segments finding
   AliPHOSPID * fPID ;                             // Method for identifying the type of particle