add cluster occupancy to the ESD friend

[u/mrichter/AliRoot.git] / JETAN / AliAnalysisTaskKMeans.cxx

diff --git a/JETAN/AliAnalysisTaskKMeans.cxx b/JETAN/AliAnalysisTaskKMeans.cxx
index fb03b20e529c20f73bc519b0ade8f3a5dbd39467..c47df6f0f5cdd85940b8519c1bab542b29a1a81c 100644 (file)
--- a/JETAN/AliAnalysisTaskKMeans.cxx
+++ b/JETAN/AliAnalysisTaskKMeans.cxx
@@ -105,8 +105,13 @@ AliAnalysisTaskKMeans::AliAnalysisTaskKMeans()
 
 {
   //
+
   // Constructor
   //
+    for (Int_t i=0; i< 10; i++) {
+       fA[i] = 0;
+       fB[i] = 0;
+    }
 }
 
 //________________________________________________________________________
@@ -153,6 +158,10 @@ AliAnalysisTaskKMeans::AliAnalysisTaskKMeans(const char *name)
   //
   // Constructor
   //
+    for (Int_t i=0; i< 10; i++) {
+       fA[i] = 0;
+       fB[i] = 0;
+    }
   DefineOutput(1,  TList::Class());
 }
 
@@ -196,6 +205,10 @@ AliAnalysisTaskKMeans::AliAnalysisTaskKMeans(const AliAnalysisTaskKMeans &res)
       ,fCuts(0)
 {
     // Dummy copy constructor
+    for (Int_t i=0; i< 10; i++) {
+       fA[i] = 0;
+       fB[i] = 0;
+    }
 }
 
 AliAnalysisTaskKMeans& AliAnalysisTaskKMeans::operator=(const AliAnalysisTaskKMeans& /*trclass*/)
@@ -294,11 +307,11 @@ void AliAnalysisTaskKMeans::UserExec(Option_t *)
 {
   // Main loop
   // Called for each event
-    Double_t   phi [500];
-    Double_t   phiR[500];
-    Double_t    eta[500];
-    Double_t   etaR[500];
-    Double_t    pt [500];
+    Double_t   phi [500] = {0};
+    Double_t   phiR[500] = {0};
+    Double_t    eta[500] = {0};
+    Double_t   etaR[500] = {0};
+    Double_t    pt [500] = {0};
 
   if (!fInputEvent) {
     Printf("ERROR: fESD not available");
@@ -492,7 +505,6 @@ void AliAnalysisTaskKMeans::UserExec(Option_t *)
   
     mPhi    = best.GetMx();
     mEta    = best.GetMy();
-    sigma2  = best.GetSigma2();
     nk      = best.GetK();
     im      = (best.GetInd())[0];
     etaC    = mEta[im];    
@@ -501,12 +513,12 @@ void AliAnalysisTaskKMeans::UserExec(Option_t *)
     // Cluster Multiplicity
     if (rmaxG > 0. && TMath::Abs(etaC) < 0.4) {
       for (Int_t i = 0; i < 1; i++) {
-       Int_t im = (best.GetInd())[i];
+       im = (best.GetInd())[i];
        fH1CEtaR->Fill(mEta[im]);
       }
       
       for (Int_t i = 0; i < 1; i++) {
-       Int_t im = (best.GetInd())[i];
+       im = (best.GetInd())[i];
        for (Int_t j = 0; j < ic; j++) {
          Double_t dphi = DeltaPhi(mPhi[im], phiR[j]);
          Double_t deta = mEta[im] - etaR[j];