Implementation to reduce output size in Pb-Pb central collisions
authorddobrigk <david.dobrigkeit.chinellato@cern.ch>
Tue, 25 Nov 2014 09:51:18 +0000 (07:51 -0200)
committerddobrigk <david.dobrigkeit.chinellato@cern.ch>
Tue, 25 Nov 2014 09:52:33 +0000 (07:52 -0200)
PWGLF/STRANGENESS/Cascades/AliAnalysisTaskExtractCascade.cxx
PWGLF/STRANGENESS/Cascades/AliAnalysisTaskExtractCascade.h
PWGLF/STRANGENESS/Cascades/AliAnalysisTaskExtractPerformanceCascade.cxx
PWGLF/STRANGENESS/Cascades/AliAnalysisTaskExtractPerformanceCascade.h

index 96645e2..4ada043 100644 (file)
@@ -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();
+          }
+      }
+      
+
   }
 
 //------------------------------------------------
index c76431b..be54906 100644 (file)
@@ -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
index 8ec4f0c..48e4dcb 100644 (file)
@@ -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.
index 1de8230..108ec93 100644 (file)
@@ -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