mPhi = best.GetMx();
mEta = best.GetMy();
- sigma2 = best.GetSigma2();
nk = best.GetK();
im = (best.GetInd())[0];
etaC = mEta[im];
pxT = new Double_t[nmctracks];
pyT = new Double_t[nmctracks];
ptT = new Double_t[nmctracks];
-
+ for (Int_t i = 0; i < nmctracks; i++)
+ {
+ pxT[i] = 0;
+ pyT[i] = 0;
+ ptT[i] = 0;
+ }
for (Int_t iMCTracks = 0; iMCTracks < nPrim; ++iMCTracks) {
TParticle* trackmc = stack->Particle(iMCTracks);
if (!trackmc) continue;
Int_t fOpt = fReader->GetReaderHeader()->GetDetector();
// check if we are reading AOD jets
- TRefArray *refs = 0;
Bool_t fromAod = !strcmp(fReader->ClassName(),"AliJetAODReader");
- if (fromAod) { refs = fReader->GetReferences(); }
// RUN ALGORITHM
// read input particles -----------------------------
Float_t BkgChargedFastJet();
Float_t BkgStat();
Float_t BkgFastJetCone(TClonesArray* fAODJets);
- Float_t BkgRemoveJetLeading(TClonesArray* fAODJets);
+// Float_t BkgRemoveJetLeading(TClonesArray* fAODJets);
Float_t BkgRemoveJetLeadingFromUArray(TClonesArray* fAODJets);
Float_t EtaToTheta(Float_t arg);
Bool_t EmcalAcceptance(const Float_t eta, const Float_t phi, const Float_t radius);
fSamplingFraction(0)
{
fgHadrCorr = this;
+ for (Int_t i = 0; i < 8; i++) fPar[i] = 0.;
}
AliJetHadronCorrectionv1*
void TrackPositionEMCal(AliAODTrack* track,Double_t &eta, Double_t &phi);
protected:
- AliJetHadronCorrectionv1():fSamplingFraction(0) {;}
+ AliJetHadronCorrectionv1():fSamplingFraction(0) {for (Int_t i = 0; i < 8; i++) fPar[i] = 0.;}
AliJetHadronCorrectionv1(const char *name, const char *title);
private:
}
//------------------------------------------------------------------------
-Bool_t AliJetUnitArray::operator>(AliJetUnitArray unit) const
+Bool_t AliJetUnitArray::operator>(AliJetUnitArray* unit) const
{
// Greater than operator used by sort
- if( fUnitEnergy > unit.GetUnitEnergy())
+ if( fUnitEnergy > unit->GetUnitEnergy())
return kTRUE;
else
return kFALSE;
}
//------------------------------------------------------------------------
-Bool_t AliJetUnitArray::operator<( AliJetUnitArray unit) const
+Bool_t AliJetUnitArray::operator<( AliJetUnitArray* unit) const
{
// Less than operator used by sort
- if( fUnitEnergy < unit.GetUnitEnergy())
+ if( fUnitEnergy < unit->GetUnitEnergy())
return kTRUE;
else
return kFALSE;
}
//------------------------------------------------------------------------
-Bool_t AliJetUnitArray::operator==( AliJetUnitArray unit) const
+Bool_t AliJetUnitArray::operator==( AliJetUnitArray* unit) const
{
// equality operator used by sort
- if( fUnitEnergy == unit.GetUnitEnergy())
+ if( fUnitEnergy == unit->GetUnitEnergy())
return kTRUE;
else
return kFALSE;
protected:
AliJetUnitArray(const AliJetUnitArray& rUnit);
AliJetUnitArray& operator = (const AliJetUnitArray& rhs);
- Bool_t operator> ( AliJetUnitArray unit1) const;
- Bool_t operator< ( AliJetUnitArray unit1) const;
- Bool_t operator== ( AliJetUnitArray unit1) const;
+ Bool_t operator> ( AliJetUnitArray* unit1) const;
+ Bool_t operator< ( AliJetUnitArray* unit1) const;
+ Bool_t operator== ( AliJetUnitArray* unit1) const;
Float_t fUnitEnergy; // Energy (Pt,et) of the unit
Float_t fUnitEta; // Eta of the unit
Int_t fOpt = fReader->GetReaderHeader()->GetDetector();
// Check if we are reading AOD jets
- TRefArray *refs = 0;
Bool_t fromAod = !strcmp(fReader->ClassName(),"AliJetAODReader");
- if (fromAod) { refs = fReader->GetReferences(); }
-
//******************************** SISCONE PLUGIN CONFIGURATION
// Here we look for SISCone parameters in the header and we define our plugin.
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);
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//
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
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