cut in on min eta (Leonardo)

fix in symmetrization

diff --git a/PWGCF/Correlations/Base/AliAnalyseLeadingTrackUE.cxx b/PWGCF/Correlations/Base/AliAnalyseLeadingTrackUE.cxx
index a8820a0fd8b0b531ad76de4705f48b66e6f69ada..680b0491c877e7be73710e0a685098a85dec03f6 100644 (file)
--- a/PWGCF/Correlations/Base/AliAnalyseLeadingTrackUE.cxx
+++ b/PWGCF/Correlations/Base/AliAnalyseLeadingTrackUE.cxx
@@ -79,6 +79,7 @@ AliAnalyseLeadingTrackUE::AliAnalyseLeadingTrackUE() :
   fOnlyHadrons(kFALSE),
   fCheckMotherPDG(kTRUE),
   fTrackEtaCut(0.8),
+  fTrackEtaCutMin(-1.),
   fTrackPtMin(0),
   fEventSelection(AliVEvent::kMB|AliVEvent::kUserDefined),
   fDCAXYCut(0),
@@ -219,7 +220,7 @@ TObjArray*  AliAnalyseLeadingTrackUE::FindLeadingObjects(TObject *obj)
        AliVParticle* part = ParticleWithCuts( obj, ipart );
         if (!part) continue;
        // Accept leading-tracks in a limited pseudo-rapidity range     
-       if( TMath::Abs(part->Eta()) > fTrackEtaCut ) continue;
+       if( TMath::Abs(part->Eta()) > fTrackEtaCut || TMath::Abs(part->Eta()) < fTrackEtaCutMin ) continue;
        tracks->AddLast( part );
        }
   // Order tracks by pT        
@@ -397,7 +398,7 @@ TObjArray* AliAnalyseLeadingTrackUE::GetAcceptedParticles(TObject* obj, TObject*
     if (!part) continue;
     
     if (useEtaPtCuts)
-      if (TMath::Abs(part->Eta()) > fTrackEtaCut || part->Pt() < fTrackPtMin)
+      if (TMath::Abs(part->Eta()) > fTrackEtaCut || TMath::Abs(part->Eta()) < fTrackEtaCutMin || part->Pt() < fTrackPtMin)
       {
        if (hasOwnership)
          delete part;
@@ -444,7 +445,7 @@ TObjArray* AliAnalyseLeadingTrackUE::GetFakeParticles(TObject* obj, TObject* arr
     if (!partReconstructed) continue;
 
     if (useEtaPtCuts)
-      if (TMath::Abs(partReconstructed->Eta()) > fTrackEtaCut || partReconstructed->Pt() < fTrackPtMin)
+      if (TMath::Abs(partReconstructed->Eta()) > fTrackEtaCut || TMath::Abs(partReconstructed->Eta()) < fTrackEtaCutMin || partReconstructed->Pt() < fTrackPtMin)
       {
         if (hasOwnership)
           delete partReconstructed;
@@ -905,7 +906,7 @@ TObjArray*  AliAnalyseLeadingTrackUE::SortRegions(const AliVParticle* leading, T
        TVector3 partVect(part->Px(), part->Py(), part->Pz());
  
        Int_t region = 0;
-       if( TMath::Abs(partVect.Eta()) > fTrackEtaCut ) continue;
+       if( TMath::Abs(partVect.Eta()) > fTrackEtaCut || TMath::Abs(partVect.Eta()) < fTrackEtaCutMin) continue;
        // transverse regions
        if (leadVect.DeltaPhi(partVect) < -k60rad && leadVect.DeltaPhi(partVect) > -k120rad )region = -1; //left
        if (leadVect.DeltaPhi(partVect) > k60rad && leadVect.DeltaPhi(partVect) < k120rad ) region = 1;   //right

diff --git a/PWGCF/Correlations/Base/AliAnalyseLeadingTrackUE.h b/PWGCF/Correlations/Base/AliAnalyseLeadingTrackUE.h
index 5faffa5c234d2b28e1a03c0f9d28d524af2e052a..593b1d71b6fdc314eebf66d313b702dcd407a922 100644 (file)
--- a/PWGCF/Correlations/Base/AliAnalyseLeadingTrackUE.h
+++ b/PWGCF/Correlations/Base/AliAnalyseLeadingTrackUE.h
@@ -43,9 +43,10 @@ class AliAnalyseLeadingTrackUE : public TObject {
   virtual ~AliAnalyseLeadingTrackUE();                                          //virtual destructor
 
   // Setters
-  void  SetParticleSelectionCriteria(Int_t filterbit, Bool_t onlyhadrons, Double_t etacut, Double_t ptMin = 0) { fFilterBit   = filterbit;
+  void  SetParticleSelectionCriteria(Int_t filterbit, Bool_t onlyhadrons, Double_t etacut, Double_t etacutMin=-1., Double_t ptMin = 0) { fFilterBit   = filterbit;
                                                                                             fOnlyHadrons = onlyhadrons;
                                                                                             fTrackEtaCut = etacut;     
+                                                                                            fTrackEtaCutMin = etacutMin;     
                                                                                             fTrackPtMin = ptMin;}
 
   void SetEventSelection(UInt_t bits) { fEventSelection = bits; }
@@ -83,6 +84,7 @@ private:
   Bool_t         fOnlyHadrons;       // consider only charged Pions, Protons and Kaons 
   Bool_t         fCheckMotherPDG;     // Check the PDG code of mother for secondaries 
   Double_t       fTrackEtaCut;       // pseudo-rapidity limit of transverse regions     
+  Double_t       fTrackEtaCutMin;       // minimum of the pseudo-rapidity limit of transverse regions     
   Double_t       fTrackPtMin;        // pt limit for selecting particles
   UInt_t         fEventSelection;    // bit for physics selection
   TFormula*      fDCAXYCut;          // additional pt dependent cut on DCA XY (only for AOD)

diff --git a/PWGCF/Correlations/Base/AliUEHist.cxx b/PWGCF/Correlations/Base/AliUEHist.cxx
index 620d58edd42364f5c1ced25491fe12ad759a79a5..6759c099d35836dfde8e0fb817ef0f1a31af9fd3 100644 (file)
--- a/PWGCF/Correlations/Base/AliUEHist.cxx
+++ b/PWGCF/Correlations/Base/AliUEHist.cxx
@@ -2594,9 +2594,13 @@ void AliUEHist::SymmetrizepTBins()
       // for symmetric bins
       THnSparse* source = (THnSparse*) target->Clone();
       
+      Int_t zVtxBins = 1;
+      if (target->GetNdimensions() > 5)
+       zVtxBins = target->GetAxis(5)->GetNbins();
+      
       // axes: 0 delta eta; 1 pT,a; 2 pT,t; 3 centrality; 4 delta phi; 5 vtx-z
       for (Int_t i3 = 1; i3 <= target->GetAxis(3)->GetNbins(); i3++)
-       for (Int_t i5 = 1; i5 <= target->GetAxis(5)->GetNbins(); i5++)
+       for (Int_t i5 = 1; i5 <= zVtxBins; i5++)
        {
          for (Int_t i1 = 1; i1 <= target->GetAxis(1)->GetNbins(); i1++)
            for (Int_t i2 = 1; i2 <= target->GetAxis(2)->GetNbins(); i2++)

diff --git a/PWGCF/Correlations/DPhi/AliAnalysisTaskPhiCorrelations.cxx b/PWGCF/Correlations/DPhi/AliAnalysisTaskPhiCorrelations.cxx
index a009c5ddf1c6052ab109132faa0a9665d793183c..5cb2e059a24eb5d52390286cac9b58b071991881 100644 (file)
--- a/PWGCF/Correlations/DPhi/AliAnalysisTaskPhiCorrelations.cxx
+++ b/PWGCF/Correlations/DPhi/AliAnalysisTaskPhiCorrelations.cxx
@@ -120,6 +120,7 @@ fZVertex(7.),
 fCentralityMethod("V0M"),
 // track cuts
 fTrackEtaCut(0.8),
+fTrackEtaCutMin(-1.),
 fOnlyOneEtaSide(0),
 fPtMin(0.5),
 fDCAXYCut(0),
@@ -218,7 +219,7 @@ void  AliAnalysisTaskPhiCorrelations::CreateOutputObjects()
    
   // Initialize class with main algorithms, event and track selection. 
   fAnalyseUE = new AliAnalyseLeadingTrackUE();
-  fAnalyseUE->SetParticleSelectionCriteria(fFilterBit, fUseChargeHadrons, fTrackEtaCut, fPtMin);
+  fAnalyseUE->SetParticleSelectionCriteria(fFilterBit, fUseChargeHadrons, fTrackEtaCut, fTrackEtaCutMin, fPtMin);
   fAnalyseUE->SetDCAXYCut(fDCAXYCut);
   fAnalyseUE->SetTrackStatus(fTrackStatus);
   fAnalyseUE->SetCheckMotherPDG(fCheckMotherPDG);
@@ -385,6 +386,7 @@ void  AliAnalysisTaskPhiCorrelations::AddSettingsTree()
   settingsTree->Branch("fZVertex", &fZVertex,"ZVertex/D");
   //settingsTree->Branch("fCentralityMethod", fCentralityMethod.Data(),"CentralityMethod/C");
   settingsTree->Branch("fTrackEtaCut", &fTrackEtaCut, "TrackEtaCut/D");
+  settingsTree->Branch("fTrackEtaCutMin", &fTrackEtaCutMin, "TrackEtaCutMin/D");
   settingsTree->Branch("fOnlyOneEtaSide", &fOnlyOneEtaSide,"OnlyOneEtaSide/I");
   settingsTree->Branch("fPtMin", &fPtMin, "PtMin/D");
   settingsTree->Branch("fFilterBit", &fFilterBit,"FilterBit/I");

diff --git a/PWGCF/Correlations/DPhi/AliAnalysisTaskPhiCorrelations.h b/PWGCF/Correlations/DPhi/AliAnalysisTaskPhiCorrelations.h
index afddb64c38f2211672a20cdf7a34c681367611a9..9897f023fa12f8a7c3f7e557bc1830ee4a22e372 100644 (file)
--- a/PWGCF/Correlations/DPhi/AliAnalysisTaskPhiCorrelations.h
+++ b/PWGCF/Correlations/DPhi/AliAnalysisTaskPhiCorrelations.h
@@ -85,6 +85,7 @@ class  AliAnalysisTaskPhiCorrelations : public AliAnalysisTask
     
     // track cuts
     void   SetTrackEtaCut( Double_t val )    { fTrackEtaCut = val; }
+    void   SetTrackEtaCutMin( Double_t val )    { fTrackEtaCutMin = val; }
     void   SetOnlyOneEtaSide(Int_t flag)     { fOnlyOneEtaSide = flag; }
     void   SetPtMin(Double_t val)            { fPtMin = val; }
     void   SetFilterBit( UInt_t val )        { fFilterBit = val;  }
@@ -174,7 +175,8 @@ class  AliAnalysisTaskPhiCorrelations : public AliAnalysisTask
     TString             fCentralityMethod;     // Method to determine centrality
     
     // Track cuts
-    Double_t           fTrackEtaCut;          // Eta cut on particles
+    Double_t           fTrackEtaCut;          // Maximum Eta cut on particles
+    Double_t           fTrackEtaCutMin;          // Minimum Eta cut on particles
     Int_t              fOnlyOneEtaSide;       // decides that only trigger particle from one eta side are considered (0 = all; -1 = negative, 1 = positive)
     Double_t            fPtMin;                // Min pT to start correlations
     TFormula*           fDCAXYCut;             // additional pt dependent cut on DCA XY (only for AOD)
@@ -206,7 +208,7 @@ class  AliAnalysisTaskPhiCorrelations : public AliAnalysisTask
     
     Bool_t fFillpT;                // fill sum pT instead of number density
     
-    ClassDef( AliAnalysisTaskPhiCorrelations, 33); // Analysis task for delta phi correlations
+    ClassDef( AliAnalysisTaskPhiCorrelations, 34); // Analysis task for delta phi correlations
   };
 
 class AliDPhiBasicParticle : public AliVParticle