From cba109a0896bf8c661fc886994ab068cf120445f Mon Sep 17 00:00:00 2001 From: kleinb Date: Fri, 16 Jul 2010 15:37:35 +0000 Subject: [PATCH] Added protection (QAMC), increased range for unfolding histograms (Spectrum) --- PWG4/JetTasks/AliAnalysisTaskJetSpectrum2.cxx | 6 +++--- PWG4/JetTasks/AliPWG4HighPtQAMC.cxx | 3 +++ 2 files changed, 6 insertions(+), 3 deletions(-) diff --git a/PWG4/JetTasks/AliAnalysisTaskJetSpectrum2.cxx b/PWG4/JetTasks/AliAnalysisTaskJetSpectrum2.cxx index 48484278b94..20734d90f8c 100644 --- a/PWG4/JetTasks/AliAnalysisTaskJetSpectrum2.cxx +++ b/PWG4/JetTasks/AliAnalysisTaskJetSpectrum2.cxx @@ -282,7 +282,7 @@ void AliAnalysisTaskJetSpectrum2::UserCreateOutputObjects() // // Histogram - const Int_t nBinPt = 240; + const Int_t nBinPt = 320; Double_t binLimitsPt[nBinPt+1]; for(Int_t iPt = 0;iPt <= nBinPt;iPt++){ if(iPt == 0){ @@ -1092,7 +1092,7 @@ void AliAnalysisTaskJetSpectrum2::MakeJetContainer(){ // link it // const Int_t kNvar = 3 ; //number of variables on the grid:pt,eta, phi - const Double_t kPtmin = 0.0, kPtmax = 240.; // we do not want to have empty bins at the beginning... + const Double_t kPtmin = 0.0, kPtmax = 320.; // we do not want to have empty bins at the beginning... const Double_t kEtamin = -3.0, kEtamax = 3.0; const Double_t kPhimin = 0., kPhimax = 2. * TMath::Pi(); const Double_t kZmin = 0., kZmax = 1; @@ -1104,7 +1104,7 @@ void AliAnalysisTaskJetSpectrum2::MakeJetContainer(){ //arrays for the number of bins in each dimension Int_t iBin[kNvar]; - iBin[0] = 240; //bins in pt + iBin[0] = 320; //bins in pt iBin[1] = 1; //bins in eta iBin[2] = 1; // bins in phi diff --git a/PWG4/JetTasks/AliPWG4HighPtQAMC.cxx b/PWG4/JetTasks/AliPWG4HighPtQAMC.cxx index c7fbdb2cbc2..8a6684fc341 100644 --- a/PWG4/JetTasks/AliPWG4HighPtQAMC.cxx +++ b/PWG4/JetTasks/AliPWG4HighPtQAMC.cxx @@ -394,11 +394,14 @@ void AliPWG4HighPtQAMC::Exec(Option_t *) { Int_t nTracks = fESD->GetNumberOfTracks(); AliDebug(2,Form("nTracks %d", nTracks)); + int nMCtracks = stack->GetNtrack(); + for (Int_t iTrack = 0; iTrack < nTracks; iTrack++) { AliESDtrack *track = fESD->GetTrack(iTrack); if(!track) continue; Int_t label = TMath::Abs(track->GetLabel()); + if(label>=nMCtracks)continue; TParticle *particle = stack->Particle(label) ; if(!particle) continue; -- 2.43.0