fix separator

[u/mrichter/AliRoot.git] / PMD / AliPMDDiscriminator.cxx

diff --git a/PMD/AliPMDDiscriminator.cxx b/PMD/AliPMDDiscriminator.cxx
index 097124622d4915343059e814f252846855641e1d..fccb270d3568729e10118007d6a799e667e1063c 100644 (file)
--- a/PMD/AliPMDDiscriminator.cxx
+++ b/PMD/AliPMDDiscriminator.cxx
@@ -24,16 +24,14 @@
 
 #include <Riostream.h>
 #include <TMath.h>
-#include <TBRIK.h>
-#include <TNode.h>
 #include <TTree.h>
-#include <TGeometry.h>
 #include <TObjArray.h>
 #include <TClonesArray.h>
 #include <TFile.h>
 #include <TNtuple.h>
 #include <TParticle.h>
 
+#include "AliLog.h"
 #include "AliPMDcluster.h"
 #include "AliPMDclupid.h"
 #include "AliPMDDiscriminator.h"
@@ -41,8 +39,7 @@
 ClassImp(AliPMDDiscriminator)
 
 AliPMDDiscriminator::AliPMDDiscriminator() :
-  fDiscrim(0),
-  fDebug(0)
+  fDiscrim(0)
 {
 //
 // Default Constructor
@@ -79,11 +76,11 @@ void AliPMDDiscriminator::EmpDiscrimination(TObjArray *pmdcontin, TObjArray *pmd
   //
   const  Int_t kumperdet = 24;
   static Int_t neibx[6]={1,0,-1,-1,0,1}, neiby[6]={0,1,1,0,-1,-1}; 
-  Int_t   det,smn;
+  Int_t   det = 0, smn = 0;
   Int_t   iprecount[24], icpvcount[24];
-  Float_t xpos,ypos;
-  Float_t adc, ncell, rad;
-  Float_t clusdata[6];
+  Float_t xpos = 0., ypos = 0.;
+  Float_t adc = 0., ncell = 0., rad = 0.;
+  Float_t clusdata[6] = {0.,0.,0.,0.,0.,0.};
 
   for(Int_t i = 0; i < kumperdet; i++)
     {
@@ -94,7 +91,8 @@ void AliPMDDiscriminator::EmpDiscrimination(TObjArray *pmdcontin, TObjArray *pmd
   AliPMDclupid   *pmdclout = 0;
 
   Int_t nentries1 = pmdcontin->GetEntries();
-  cout << " nentries1 = " << nentries1 << endl;
+
+  AliDebug(1,Form("Number of total clusters from CPV PRE = %d",nentries1));
   for (Int_t ient1 = 0; ient1 < nentries1; ient1++)
     {
       pmdcl = (AliPMDcluster*)pmdcontin->UncheckedAt(ient1);
@@ -105,7 +103,7 @@ void AliPMDDiscriminator::EmpDiscrimination(TObjArray *pmdcontin, TObjArray *pmd
       if(det == 1) icpvcount[smn]++;
     } // Entries of TObjArray loop
 
-  Int_t   idet, ismn;
+  Int_t   idet = 0, ismn = 0;
   Float_t edepcpv[48][96];
   Int_t   statuscpv[48][96];
 
@@ -216,7 +214,7 @@ void AliPMDDiscriminator::EmpDiscrimination(TObjArray *pmdcontin, TObjArray *pmd
          clusdata[2] = adcpre[jjsort];
          clusdata[3] = ncellpre[jjsort];
          clusdata[4] = radpre[jjsort];
-         clusdata[5] = (Float_t) clupidpre[jjsort];
+         //PH    clusdata[5] = (Float_t) clupidpre[jjsort];
 
          // Temporary the cluster PID is set to 1 if the
          // adc > 3MIP units which will be changed later on.