Int_t* cFlagT = new Int_t[nIn]; // Temporarily added
Int_t* sFlagT = new Int_t[nIn]; // Temporarily added
Int_t* injet = new Int_t[nIn];
-
+
+ for (Int_t i = 0; i < nIn; i++) {
+ ptT[i] = 0.;
+ etaT[i] = 0.;
+ phiT[i] = 0.;
+ cFlagT[i] = 0;
+ sFlagT[i] = 0;
+ injet[i] = 0;
+ }
//total energy in array
Float_t etbgTotal = 0.0;
TH1F* hPtTotal = new TH1F("hPt","Pt distribution of all particles ",100,0.0,15.0);
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
Float_t *etaCell2 = new Float_t[nIn]; //! Cell eta - Extracted from UnitArray
Float_t *phiCell2 = new Float_t[nIn]; //! Cell phi - Extracted from UnitArray
Int_t *flagCell2 = new Int_t[nIn]; //! Cell phi - Extracted from UnitArray
-
+ for(Int_t i = 0; i < nIn; i++) {
+ etCell[i] = 0.;
+ etaCell[i] = 0.;
+ phiCell[i] = 0.;
+ flagCell[i] = 0;
+ etCell2[i] = 0.;
+ etaCell2[i] = 0.;
+ phiCell2[i] = 0.;
+ flagCell2[i] = 0;
+ }
// Information extracted from fUnitArray
// Load input vectors and calculate total energy in array
for(Int_t i=0; i<nIn; i++)
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//
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//
AliUA1JetHeaderV1* header = (AliUA1JetHeaderV1*) fHeader;
fOpt = fReader->GetReaderHeader()->GetDetector();
Float_t rc= header->GetRadius();
- Float_t etIn[30];
+ Float_t etIn[30] = {0.};
Float_t etOut = 0;
-
- for(Int_t j=0;j<30;j++){etIn[j]=0.;}
-
+
for(Int_t jpart = 0; jpart < nIn; jpart++){ // loop for all particles in array
for(Int_t ijet=0; ijet<nJ; ijet++){
//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
//calculate energy inside
Float_t rc= header->GetRadius();
- Float_t etIn[30];
+ Float_t etIn[30] = {0.};
for(Int_t jpart = 0; jpart < nIn; jpart++)
{ // loop for all particles in array
////////////////////////////////////////////////////////////////////////
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
}