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
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
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//
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//
//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
////////////////////////////////////////////////////////////////////////
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)
{
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);
}
////////////////////////////////////////////////////////////////////////
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)
{
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
}