update for bug #100726

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

diff --git a/JETAN/AliUA1JetFinderV2.cxx b/JETAN/AliUA1JetFinderV2.cxx
index 30980248105ae944567c16949e98d07e8ce5b12e..0e73b96a28d330625d4f67f0980b5caa628cdd2b 100644 (file)
--- a/JETAN/AliUA1JetFinderV2.cxx
+++ b/JETAN/AliUA1JetFinderV2.cxx
@@ -335,10 +335,27 @@ void AliUA1JetFinderV2::FindJets()
   Int_t      loop1     = 0;
   Int_t*     injet     = new Int_t[nCandidate];
   Int_t*     sflag     = new Int_t[nCandidate];
+  
+  for(Int_t i = 0; i < nCandidate; i++) {
+      ptT[i]       = 0;
+      en2T[i]      = 0;
+      pt2T[i]      = 0;
+      detT[i]      = 0;
+      etaT[i]      = 0;
+      phiT[i]      = 0;
+      cFlagT[i]    = 0;
+      cFlag2T[i]   = 0;
+      sFlagT[i]    = 0;
+      cClusterT[i] = 0;
+      vectT[i]     = 0;
+      injet[i]     = 0;
+      sflag[i]     = 0;
+}
+
   TRefArray* trackRef  = new TRefArray();
 
   //total energy in array
-  Float_t  etbgTotal = 0.0;
+  Float_t etbgTotal = 0.0;
   TH1F*    hPtTotal  = new TH1F("hPt","Pt distribution of all particles ",100,0.0,15.0);
 
   // Input cell info
@@ -354,11 +371,11 @@ void AliUA1JetFinderV2::FindJets()
     etCell[i]    = 0.;
     etaCell[i]   = 0.;
     phiCell[i]   = 0.;
-    flagCell[i]  = 0.;
+    flagCell[i]  = 0;
     etCell2[i]   = 0.;
     etaCell2[i]  = 0.;
     phiCell2[i]  = 0.;
-    flagCell2[i] = 0.;
+    flagCell2[i] = 0;
   }
   // Information extracted from fUnitArray
   // Load input vectors and calculate total energy in array
@@ -763,7 +780,7 @@ void AliUA1JetFinderV2::RunAlgoritm(Int_t nIn, Float_t* etCell, Float_t* const e
   Float_t etaAlgoJet[30]    = {0.};
   Float_t phiAlgoJet[30]    = {0.};
   Float_t etAlgoJet[30]     = {0.};
-  Int_t   ncellsAlgoJet[30] = {0.};
+  Int_t   ncellsAlgoJet[30] = {0};
 
   // Run algorithm//
 
@@ -966,10 +983,10 @@ void AliUA1JetFinderV2::RunAlgoritmC(Float_t etbgTotal, Double_t dEtTotal, Int_t
   Float_t etseed  = header->GetEtSeed();
 
   // Tmp array of jets form algoritm
-  Float_t etaAlgoJet[30];
-  Float_t phiAlgoJet[30];
-  Float_t etAlgoJet[30];
-  Int_t   ncellsAlgoJet[30];
+  Float_t etaAlgoJet[30]    = {0.};
+  Float_t phiAlgoJet[30]    = {0.};
+  Float_t etAlgoJet[30]     = {0.};
+  Int_t   ncellsAlgoJet[30] = {0};
 
   // Run algorithm//
 
@@ -1345,7 +1362,7 @@ void AliUA1JetFinderV2::SubtractBackgC(const Int_t& nIn, const Int_t&nJ, Float_t
   //calculate energy inside and outside cones
   AliUA1JetHeaderV1* header = (AliUA1JetHeaderV1*) fHeader;
   Float_t rc= header->GetRadius();
-  Float_t etIn[30];
+  Float_t etIn[30] = {0.};
   Float_t etOut = 0;
   for(Int_t jpart = 0; jpart < nIn; jpart++){ // loop for all particles in array
     // if((fReader->GetCutFlag(jpart)) != 1) continue; // pt cut
@@ -1467,7 +1484,7 @@ void AliUA1JetFinderV2::SubtractBackgStat(const Int_t& nIn, const Int_t&nJ,Float
 
 ////////////////////////////////////////////////////////////////////////
 void AliUA1JetFinderV2::SubtractBackgCone(const Int_t& nIn, const Int_t&nJ,Float_t& etbgTotalN,
-                      Float_t* ptT, Float_t* etaT, const Float_t* phiT, const Int_t* cFlagT, const Int_t* sFlagT,
+                      const Float_t* ptT, const Float_t* etaT, const Float_t* phiT, const Int_t* cFlagT, const Int_t* sFlagT,
                       Float_t* const etJet, const Float_t* etaJet, const Float_t* phiJet,
                       Float_t* const etsigJet, Int_t* const multJet, Int_t* const injet)
 {
@@ -1484,7 +1501,7 @@ void AliUA1JetFinderV2::SubtractBackgCone(const Int_t& nIn, const Int_t&nJ,Float
   TH1F* hAreaJet[30];
   for(Int_t mjet=0; mjet<nJ; mjet++){
     char hEtname[256]; char hAreaname[256];
-    sprintf(hEtname, "hEtJet%d", mjet); sprintf(hAreaname, "hAreaJet%d", mjet);
+    snprintf(hEtname, 256, "hEtJet%d", mjet); snprintf(hAreaname, 256, "hAreaJet%d", mjet);
     hEtJet[mjet] = new TH1F(hEtname,"et dist in eta ",ndiv,etamin,etamax);
     hAreaJet[mjet] = new TH1F(hAreaname,"area dist in eta ",ndiv,etamin,etamax);
   }
@@ -1581,7 +1598,7 @@ void AliUA1JetFinderV2::SubtractBackgCone(const Int_t& nIn, const Int_t&nJ,Float
 
 ////////////////////////////////////////////////////////////////////////
 void AliUA1JetFinderV2::SubtractBackgRatio(const Int_t& nIn, const Int_t&nJ,Float_t& etbgTotalN,
-                      Float_t* ptT, Float_t* etaT, const Float_t* phiT, const Int_t* cFlagT, const Int_t* sFlagT,
+                      const Float_t* ptT, const Float_t* etaT, const Float_t* phiT, const Int_t* cFlagT, const Int_t* sFlagT,
                       Float_t* const etJet, const Float_t* etaJet, const Float_t* phiJet,
                       Float_t* const etsigJet, Int_t* const multJet, Int_t* const injet)
 {
@@ -1601,7 +1618,7 @@ void AliUA1JetFinderV2::SubtractBackgRatio(const Int_t& nIn, const Int_t&nJ,Floa
   TH1F* hAreaJet[30];
   for(Int_t mjet=0; mjet<nJ; mjet++){
     char hEtname[256]; char hAreaname[256];
-    sprintf(hEtname, "hEtJet%d", mjet); sprintf(hAreaname, "hAreaJet%d", mjet);
+    snprintf(hEtname, 256, "hEtJet%d", mjet); snprintf(hAreaname, 256, "hAreaJet%d", mjet);
     hEtJet[mjet] = new TH1F(hEtname,"et dist in eta ",ndiv,etamin,etamax);        // change range
     hAreaJet[mjet] = new TH1F(hAreaname,"area dist in eta ",ndiv,etamin,etamax);  // change range
   }