include also negative pT for the jet spectrum

diff --git a/PWG4/JetTasks/AliAnalysisTaskJetSpectrum2.cxx b/PWG4/JetTasks/AliAnalysisTaskJetSpectrum2.cxx
index 86b3b64682905dfecfb04f5c1deca211ddd27dd5..e96acdb3f9694570a7ff12fe873c643606c06ba6 100644 (file)
--- a/PWG4/JetTasks/AliAnalysisTaskJetSpectrum2.cxx
+++ b/PWG4/JetTasks/AliAnalysisTaskJetSpectrum2.cxx
@@ -349,11 +349,11 @@ void AliAnalysisTaskJetSpectrum2::UserCreateOutputObjects()
     
 
 
-  const Int_t nBinPt = 100;
+  const Int_t nBinPt = 120;
   Double_t binLimitsPt[nBinPt+1];
   for(Int_t iPt = 0;iPt <= nBinPt;iPt++){
     if(iPt == 0){
-      binLimitsPt[iPt] = 0.0;
+      binLimitsPt[iPt] = -50.0;
     }
     else {// 1.0
       binLimitsPt[iPt] =  binLimitsPt[iPt-1] + 2.5;
@@ -475,8 +475,8 @@ void AliAnalysisTaskJetSpectrum2::UserCreateOutputObjects()
 
     // Bins:  Jet number: pTJet, cent, mult, RP, Area.   total bins = 4.5M
     const Int_t nBinsSparse1 = 6;
-    const Int_t nBins1[nBinsSparse1] = {     kMaxJets+1,100, 10,  25,    fNRPBins, 10};
-    const Double_t xmin1[nBinsSparse1]  = {        -0.5,  0,  0,   0,        -0.5, 0.};
+    const Int_t nBins1[nBinsSparse1] = {     kMaxJets+1,120, 10,  25,    fNRPBins, 10};
+    const Double_t xmin1[nBinsSparse1]  = {        -0.5,-50,  0,   0,        -0.5, 0.};
     const Double_t xmax1[nBinsSparse1]  = {kMaxJets+0.5,250,100,5000,fNRPBins-0.5,1.0};
     
     fhnJetPt[ij] = new THnSparseF(Form("fhnJetPt%s",cAdd.Data()),";jet number;p_{T,jet};cent;# tracks;RP;area",nBinsSparse1,nBins1,xmin1,xmax1);
@@ -894,6 +894,7 @@ void AliAnalysisTaskJetSpectrum2::FillJetHistos(TList &jetsList,TList &particles
   for(int ij = 0;ij < nJets;ij++){
     AliAODJet *jet = (AliAODJet*)jetsList.At(ij);
     Float_t ptJet = jet->Pt();
+    if(ptJet<1.)jet->GetPtSubtracted(0);
     if(jet->Trigger()&fJetTriggerExcludeMask){
       fh1PtJetsInRej[iType]->Fill(ptJet);
       continue;

diff --git a/PWG4/macros/AddTaskPWG4HighPtSpectra.C b/PWG4/macros/AddTaskPWG4HighPtSpectra.C
index 5772ca29c6903b593a0dd2f8071529404096ca66..5e231372b54b5c39b66b1ba9c0fe5c7cb69fa5b3 100644 (file)
--- a/PWG4/macros/AddTaskPWG4HighPtSpectra.C
+++ b/PWG4/macros/AddTaskPWG4HighPtSpectra.C
@@ -187,7 +187,7 @@ AliPWG4HighPtSpectra* AddTaskPWG4HighPtSpectra(char *prodType = "LHC10e14", Bool
     // no requirements on SPD and ITSrefit failed
     trackCuts = CreateTrackCutsPWG4(10041005);   //no ITSrefit requirement
     trackCutsReject = CreateTrackCutsPWG4(1005); //ITSrefit requirement
-    trackCutsReject->SetEtaRange(-etamin,etamax);
+    trackCutsReject->SetEtaRange(etamin,etamax);
     trackCutsReject->SetPtRange(0.15, 1e10);
   }
   if(trackType==7 && cuts==1) {