Implementation to reduce output size in Pb-Pb central collisions

diff --git a/PWGLF/STRANGENESS/Cascades/AliAnalysisTaskExtractCascade.cxx b/PWGLF/STRANGENESS/Cascades/AliAnalysisTaskExtractCascade.cxx
index 96645e2502675ad9328ac86a0d48d321c992e0e1..4ada043069299686ee169d63dd37389981ab1872 100644 (file)
--- a/PWGLF/STRANGENESS/Cascades/AliAnalysisTaskExtractCascade.cxx
+++ b/PWGLF/STRANGENESS/Cascades/AliAnalysisTaskExtractCascade.cxx
@@ -108,6 +108,7 @@ AliAnalysisTaskExtractCascade::AliAnalysisTaskExtractCascade()
 fkSelectCentrality (kFALSE),
 fCentSel_Low(0.0),
 fCentSel_High(0.0),
+fLowPtCutoff(0.0),
 //------------------------------------------------
 // Tree Variables
 //------------------------------------------------
@@ -262,6 +263,7 @@ AliAnalysisTaskExtractCascade::AliAnalysisTaskExtractCascade(const char *name)
    fkSelectCentrality (kFALSE),
     fCentSel_Low(0.0),
     fCentSel_High(0.0),
+fLowPtCutoff(0.0),
 //------------------------------------------------
 // Tree Variables
 //------------------------------------------------
@@ -1493,7 +1495,19 @@ void AliAnalysisTaskExtractCascade::UserExec(Option_t *)
 
   if( (fTreeCascVarMassAsXi<1.32+0.075&&fTreeCascVarMassAsXi>1.32-0.075) ||
       (fTreeCascVarMassAsOmega<1.68+0.075&&fTreeCascVarMassAsOmega>1.68-0.075) ){
-      fTreeCascade->Fill();
+      
+      if( !fkIsNuclear ) fTreeCascade->Fill();
+      if( fkIsNuclear  ){
+          //Extra selections in case this is a nuclear collision...
+          if (TMath::Abs(fTreeCascVarNegEta) < 0.8 &&
+              TMath::Abs(fTreeCascVarPosEta) < 0.8 &&
+              TMath::Abs(fTreeCascVarBachEta) < 0.8 &&
+              fTreeCascVarPt > fLowPtCutoff){ //beware ptMC and ptreco differences
+                fTreeCascade->Fill();
+          }
+      }
+      
+
   }
 
 //------------------------------------------------

diff --git a/PWGLF/STRANGENESS/Cascades/AliAnalysisTaskExtractCascade.h b/PWGLF/STRANGENESS/Cascades/AliAnalysisTaskExtractCascade.h
index c76431baa30a3d401890da7915378429ffc8d262..be549068e8189e9e886e4be072a7d7e5a30e272f 100644 (file)
--- a/PWGLF/STRANGENESS/Cascades/AliAnalysisTaskExtractCascade.h
+++ b/PWGLF/STRANGENESS/Cascades/AliAnalysisTaskExtractCascade.h
@@ -74,7 +74,11 @@ class AliAnalysisTaskExtractCascade : public AliAnalysisTaskSE {
         fCentSel_High = lCentSelHigh;
     }
     //---------------------------------------------------------------------------------------
-//Setters for the V0 Vertexer Parameters
+    void SetLowPtCutoff ( Double_t lLowPtCutoff = 1.0) {
+        fLowPtCutoff = lLowPtCutoff;
+    }
+    //---------------------------------------------------------------------------------------
+    //Setters for the V0 Vertexer Parameters
   void SetV0VertexerMaxChisquare   ( Double_t lParameter ){ fV0VertexerSels[0] = lParameter; }
   void SetV0VertexerDCAFirstToPV   ( Double_t lParameter ){ fV0VertexerSels[1] = lParameter; }
   void SetV0VertexerDCASecondtoPV  ( Double_t lParameter ){ fV0VertexerSels[2] = lParameter; }
@@ -122,6 +126,8 @@ class AliAnalysisTaskExtractCascade : public AliAnalysisTaskSE {
     
     Double_t fCentSel_Low;
     Double_t fCentSel_High;
+    
+    Double_t fLowPtCutoff; //Reduction of data volume
 
        //Double_t        fV0Sels[7];                     // Array to store the 7 values for the different selections V0 related
        //Double_t        fCascSels[8];                   // Array to store the 8 values for the different selections Casc. related

diff --git a/PWGLF/STRANGENESS/Cascades/AliAnalysisTaskExtractPerformanceCascade.cxx b/PWGLF/STRANGENESS/Cascades/AliAnalysisTaskExtractPerformanceCascade.cxx
index 8ec4f0c32458af26a4421e621c8ad7ac20279f7b..48e4dcbc033a32ba52b6381f31619de0f946e6ee 100644 (file)
--- a/PWGLF/STRANGENESS/Cascades/AliAnalysisTaskExtractPerformanceCascade.cxx
+++ b/PWGLF/STRANGENESS/Cascades/AliAnalysisTaskExtractPerformanceCascade.cxx
@@ -108,6 +108,7 @@ AliAnalysisTaskExtractPerformanceCascade::AliAnalysisTaskExtractPerformanceCasca
 fkSelectCentrality (kFALSE),
 fCentSel_Low(0.0),
 fCentSel_High(0.0),
+fLowPtCutoff(0.0),
 //------------------------------------------------
 // Tree Variables
 //------------------------------------------------
@@ -384,6 +385,7 @@ AliAnalysisTaskExtractPerformanceCascade::AliAnalysisTaskExtractPerformanceCasca
 fkSelectCentrality (kFALSE),
 fCentSel_Low(0.0),
 fCentSel_High(0.0),
+fLowPtCutoff(0.0),
 //------------------------------------------------
 // Tree Variables
 //------------------------------------------------
@@ -2884,10 +2886,20 @@ void AliAnalysisTaskExtractPerformanceCascade::UserExec(Option_t *)
   //Xi    Mass window: 150MeV wide
   //Omega mass window: 150MeV wide
 
-  if( (fTreeCascVarMassAsXi<1.32+0.075&&fTreeCascVarMassAsXi>1.32-0.075) ||
-      (fTreeCascVarMassAsOmega<1.68+0.075&&fTreeCascVarMassAsOmega>1.68-0.075) ){
-      fTreeCascade->Fill();
-  }
+      if( (fTreeCascVarMassAsXi<1.32+0.075&&fTreeCascVarMassAsXi>1.32-0.075) ||
+         (fTreeCascVarMassAsOmega<1.68+0.075&&fTreeCascVarMassAsOmega>1.68-0.075) ){
+          
+          if( !fkIsNuclear ) fTreeCascade->Fill();
+          if( fkIsNuclear  ){
+              //Extra selections in case this is a nuclear collision...
+              if (TMath::Abs(fTreeCascVarNegEta) < 0.8 &&
+                  TMath::Abs(fTreeCascVarPosEta) < 0.8 &&
+                  TMath::Abs(fTreeCascVarBachEta) < 0.8 &&
+                  fTreeCascVarPt > fLowPtCutoff){ //beware ptMC and ptreco differences
+                  fTreeCascade->Fill();
+              }
+          }
+      }
 
 //------------------------------------------------
 // Fill tree over.

diff --git a/PWGLF/STRANGENESS/Cascades/AliAnalysisTaskExtractPerformanceCascade.h b/PWGLF/STRANGENESS/Cascades/AliAnalysisTaskExtractPerformanceCascade.h
index 1de82306022defb3ea9b8aaef0ca1072b5f4b19d..108ec93c1fafd7601f6a1c3baaa726aa4888bfc8 100644 (file)
--- a/PWGLF/STRANGENESS/Cascades/AliAnalysisTaskExtractPerformanceCascade.h
+++ b/PWGLF/STRANGENESS/Cascades/AliAnalysisTaskExtractPerformanceCascade.h
@@ -96,6 +96,10 @@ class AliAnalysisTaskExtractPerformanceCascade : public AliAnalysisTaskSE {
         fCentSel_High = lCentSelHigh;
     }
     //---------------------------------------------------------------------------------------
+    void SetLowPtCutoff ( Double_t lLowPtCutoff = 1.0) {
+        fLowPtCutoff = lLowPtCutoff;
+    }
+    //---------------------------------------------------------------------------------------
     
  private:
         // Note : In ROOT, "//!" means "do not stream the data from Master node to Worker node" ...
@@ -127,6 +131,7 @@ class AliAnalysisTaskExtractPerformanceCascade : public AliAnalysisTaskSE {
     Bool_t    fkSelectCentrality; //Switch to skip anything other than 60-80% V0M
     Double_t fCentSel_Low;
     Double_t fCentSel_High;
+    Double_t fLowPtCutoff; //Reduction of data volume
     
        //Double_t        fV0Sels[7];                     // Array to store the 7 values for the different selections V0 related
        //Double_t        fCascSels[8];                   // Array to store the 8 values for the different selections Casc. related