Added configurable eta range for multiplicity estimator for pp vs multiplicity analysis.

Will be tunable at train configuration level.

(+ bug fix in unrelated pPb testing code)

---David

diff --git a/PWGLF/STRANGENESS/LambdaK0/AliAnalysisTaskExtractPerformanceV0.cxx b/PWGLF/STRANGENESS/LambdaK0/AliAnalysisTaskExtractPerformanceV0.cxx
index cbeefb6c8b1fbdb6827f828248a5a9eb6f1f0fd9..b938d927f687a7ea0cd6231933d9f138a33a8e7f 100644 (file)
--- a/PWGLF/STRANGENESS/LambdaK0/AliAnalysisTaskExtractPerformanceV0.cxx
+++ b/PWGLF/STRANGENESS/LambdaK0/AliAnalysisTaskExtractPerformanceV0.cxx
@@ -116,6 +116,7 @@ AliAnalysisTaskExtractPerformanceV0::AliAnalysisTaskExtractPerformanceV0()
   fkSaveAssociatedOnly (kFALSE),
   fkDoNotCallTPCdEdx( kFALSE ),
   fDiffractiveOnly(kFALSE),
+  fEtaRefMult(0.5),
 //------------------------------------------------
 // Tree Variables 
 
@@ -478,6 +479,7 @@ AliAnalysisTaskExtractPerformanceV0::AliAnalysisTaskExtractPerformanceV0(const c
   fkSaveAssociatedOnly (kFALSE),
   fkDoNotCallTPCdEdx( kFALSE ),
   fDiffractiveOnly(kFALSE),
+  fEtaRefMult (0.5),
 //------------------------------------------------
 // Tree Variables 
 
@@ -2146,7 +2148,7 @@ void AliAnalysisTaskExtractPerformanceV0::UserExec(Option_t *)
     Int_t lMultiplicitySPD = -100;
   
    //testing purposes
-   if(fkIsNuclear == kFALSE) lMultiplicity = fESDtrackCuts->GetReferenceMultiplicity(lESDevent, AliESDtrackCuts::kTrackletsITSTPC,0.5);
+   if(fkIsNuclear == kFALSE) lMultiplicity = fESDtrackCuts->GetReferenceMultiplicity(lESDevent, AliESDtrackCuts::kTrackletsITSTPC,fEtaRefMult);
 
    //---> If this is a nuclear collision, then go nuclear on "multiplicity" variable...
    //---> Warning: Experimental

diff --git a/PWGLF/STRANGENESS/LambdaK0/AliAnalysisTaskExtractPerformanceV0.h b/PWGLF/STRANGENESS/LambdaK0/AliAnalysisTaskExtractPerformanceV0.h
index 3c1cd3f4b4c786e81e4cf9794146dca89320a7f7..9dd7b948f2a4f824e9113a3a70457587de07e0d6 100644 (file)
--- a/PWGLF/STRANGENESS/LambdaK0/AliAnalysisTaskExtractPerformanceV0.h
+++ b/PWGLF/STRANGENESS/LambdaK0/AliAnalysisTaskExtractPerformanceV0.h
@@ -82,6 +82,7 @@ class AliAnalysisTaskExtractPerformanceV0 : public AliAnalysisTaskSE {
   void SetV0VertexerMaxRadius      ( Double_t lParameter ){ fV0Sels[6] = lParameter; }
   //---------------------------------------------------------------------------------------
   void SetDiffractiveOnly ( Bool_t lDiffractiveOnly = kTRUE ) { fDiffractiveOnly = lDiffractiveOnly; }
+  void SetEtaRefMult ( Double_t lEtaRefMult = 0.5 ) { fEtaRefMult = lEtaRefMult; }
 
   
  private:
@@ -118,6 +119,7 @@ class AliAnalysisTaskExtractPerformanceV0 : public AliAnalysisTaskSE {
   Double_t        fV0Sels[7];                     // Array to store the 7 values for the different selections V0 related
 
   Bool_t fDiffractiveOnly; //Only look at diffractive generated events
+  Bool_t fEtaRefMult; // Ref Mult eta range 
 
 //===========================================================================================
 //   Variables for Tree, V0s

diff --git a/PWGLF/STRANGENESS/LambdaK0/AliAnalysisTaskExtractPerformanceV0pPb.cxx b/PWGLF/STRANGENESS/LambdaK0/AliAnalysisTaskExtractPerformanceV0pPb.cxx
index a957a9b6bb3a57c4cec05b6c13db035913b8fc46..6eccc7c176616690da04e4cf3fc278cc81d2199e 100644 (file)
--- a/PWGLF/STRANGENESS/LambdaK0/AliAnalysisTaskExtractPerformanceV0pPb.cxx
+++ b/PWGLF/STRANGENESS/LambdaK0/AliAnalysisTaskExtractPerformanceV0pPb.cxx
@@ -949,24 +949,24 @@ void AliAnalysisTaskExtractPerformanceV0pPb::UserExec(Option_t *)
       if ( lMCstack->IsPhysicalPrimary(iMc)!=kTRUE ) continue;
 
          if( lPdgcodeCurrentPart == 3122 ){
-            f3dHist_Generated_PtVsYVsV0A_Lambda->Fill(lPtCurrentPart, lRapCurrentPart, lMultiplicity);
-            f3dHist_Generated_PtVsYCMSVsV0A_Lambda->Fill(lPtCurrentPart, lRapCurrentPart+lpARapidityShift, lMultiplicity);
+            f3dHist_Analysis_PtVsYVsV0A_Lambda->Fill(lPtCurrentPart, lRapCurrentPart, lMultiplicity);
+            f3dHist_Analysis_PtVsYCMSVsV0A_Lambda->Fill(lPtCurrentPart, lRapCurrentPart+lpARapidityShift, lMultiplicity);
          }
          if( lPdgcodeCurrentPart == -3122 ){
-            f3dHist_Generated_PtVsYVsV0A_AntiLambda->Fill(lPtCurrentPart, lRapCurrentPart, lMultiplicity);
-            f3dHist_Generated_PtVsYCMSVsV0A_AntiLambda->Fill(lPtCurrentPart, lRapCurrentPart+lpARapidityShift, lMultiplicity);
+            f3dHist_Analysis_PtVsYVsV0A_AntiLambda->Fill(lPtCurrentPart, lRapCurrentPart, lMultiplicity);
+            f3dHist_Analysis_PtVsYCMSVsV0A_AntiLambda->Fill(lPtCurrentPart, lRapCurrentPart+lpARapidityShift, lMultiplicity);
          }
          if( lPdgcodeCurrentPart == 310 ){
-            f3dHist_Generated_PtVsYVsV0A_K0Short->Fill(lPtCurrentPart, lRapCurrentPart, lMultiplicity);
-            f3dHist_Generated_PtVsYCMSVsV0A_K0Short->Fill(lPtCurrentPart, lRapCurrentPart+lpARapidityShift, lMultiplicity);
+            f3dHist_Analysis_PtVsYVsV0A_K0Short->Fill(lPtCurrentPart, lRapCurrentPart, lMultiplicity);
+            f3dHist_Analysis_PtVsYCMSVsV0A_K0Short->Fill(lPtCurrentPart, lRapCurrentPart+lpARapidityShift, lMultiplicity);
          }
          if( lPdgcodeCurrentPart == 321 ){
-            f3dHist_Generated_PtVsYVsV0A_KPlus->Fill(lPtCurrentPart, lRapCurrentPart, lMultiplicity);
-            f3dHist_Generated_PtVsYCMSVsV0A_KPlus->Fill(lPtCurrentPart, lRapCurrentPart+lpARapidityShift, lMultiplicity);
+            f3dHist_Analysis_PtVsYVsV0A_KPlus->Fill(lPtCurrentPart, lRapCurrentPart, lMultiplicity);
+            f3dHist_Analysis_PtVsYCMSVsV0A_KPlus->Fill(lPtCurrentPart, lRapCurrentPart+lpARapidityShift, lMultiplicity);
          }
          if( lPdgcodeCurrentPart == -321 ){
-            f3dHist_Generated_PtVsYVsV0A_KMinus->Fill(lPtCurrentPart, lRapCurrentPart, lMultiplicity);
-            f3dHist_Generated_PtVsYCMSVsV0A_KMinus->Fill(lPtCurrentPart, lRapCurrentPart+lpARapidityShift, lMultiplicity);
+            f3dHist_Analysis_PtVsYVsV0A_KMinus->Fill(lPtCurrentPart, lRapCurrentPart, lMultiplicity);
+            f3dHist_Analysis_PtVsYCMSVsV0A_KMinus->Fill(lPtCurrentPart, lRapCurrentPart+lpARapidityShift, lMultiplicity);
          }
    }
 

diff --git a/PWGLF/STRANGENESS/LambdaK0/AliAnalysisTaskExtractV0.cxx b/PWGLF/STRANGENESS/LambdaK0/AliAnalysisTaskExtractV0.cxx
index 8a9f5d99ce237b38a7a49f2004e88ea6be40dca3..b43235a9bc26ef2f3b44d5187c5dd66b7ecf78f4 100644 (file)
--- a/PWGLF/STRANGENESS/LambdaK0/AliAnalysisTaskExtractV0.cxx
+++ b/PWGLF/STRANGENESS/LambdaK0/AliAnalysisTaskExtractV0.cxx
@@ -105,6 +105,7 @@ AliAnalysisTaskExtractV0::AliAnalysisTaskExtractV0()
   fkSpecialExecution( kFALSE ),
   fkSkipTrigger   ( kFALSE ),
   fTPCdEdxSelection ( kTRUE ),
+  fEtaRefMult ( 0.5 ),
 //------------------------------------------------
 // Initialize 
        fTreeVariableChi2V0(0),
@@ -251,6 +252,7 @@ AliAnalysisTaskExtractV0::AliAnalysisTaskExtractV0(const char *name)
   fkSpecialExecution( kFALSE ),
   fkSkipTrigger   ( kFALSE ),
   fTPCdEdxSelection ( kTRUE ),
+  fEtaRefMult ( 0.5 ),
 //------------------------------------------------
 // Initialize 
        fTreeVariableChi2V0(0),
@@ -852,7 +854,7 @@ void AliAnalysisTaskExtractV0::UserExec(Option_t *)
    Int_t lMultiplicityTRK = -100;
    Int_t lMultiplicitySPD = -100;  
 
-   if(fkIsNuclear == kFALSE) lMultiplicity = fESDtrackCuts->GetReferenceMultiplicity(lESDevent, AliESDtrackCuts::kTrackletsITSTPC,0.5);
+   if(fkIsNuclear == kFALSE) lMultiplicity = fESDtrackCuts->GetReferenceMultiplicity(lESDevent, AliESDtrackCuts::kTrackletsITSTPC,fEtaRefMult);
 
    //---> If this is a nuclear collision, then go nuclear on "multiplicity" variable...
    //---> Warning: Experimental

diff --git a/PWGLF/STRANGENESS/LambdaK0/AliAnalysisTaskExtractV0.h b/PWGLF/STRANGENESS/LambdaK0/AliAnalysisTaskExtractV0.h
index 875ead5782abf6236f2960e72bd2e857a933c7a2..9078350659171b0e821d248c6807fc2f18892932 100644 (file)
--- a/PWGLF/STRANGENESS/LambdaK0/AliAnalysisTaskExtractV0.h
+++ b/PWGLF/STRANGENESS/LambdaK0/AliAnalysisTaskExtractV0.h
@@ -86,6 +86,7 @@ class AliAnalysisTaskExtractV0 : public AliAnalysisTaskSE {
 //---------------------------------------------------------------------------------------
   
   void SetTPCdEdxSelection ( Bool_t lTPCdEdxSelection = kTRUE ) { fTPCdEdxSelection = lTPCdEdxSelection; }
+  void SetEtaRefMult ( Double_t lEtaRefMult = 0.5 ) { fEtaRefMult = lEtaRefMult; }
   
  private:
                                // Note : In ROOT, "//!" means "do not stream the data from Master node to Worker node" ...
@@ -115,9 +116,11 @@ class AliAnalysisTaskExtractV0 : public AliAnalysisTaskSE {
   Bool_t fkSkipTrigger; //To be used with ::SetCollisionCandidates
   
   Double_t        fV0Sels[7];                     // Array to store the 7 values for the different selections V0 related
+//  Double_t fEtaRefMult; //Reference multiplicity eta
 
   //Extra selections
   Bool_t fTPCdEdxSelection; //Configuration to apply extra TPC dE/dx selection for better filling of tree
+  Double_t fEtaRefMult; //Reference multiplicity eta
 
   //Variables for Tree
        Float_t fTreeVariableChi2V0;         //!