Configurable eta cut for stack particles to avoid large outputs
authorddobrigk <david.dobrigkeit.chinellato@cern.ch>
Fri, 10 Oct 2014 11:39:32 +0000 (08:39 -0300)
committerddobrigk <david.dobrigkeit.chinellato@cern.ch>
Fri, 10 Oct 2014 11:40:13 +0000 (08:40 -0300)
PWGLF/SPECTRA/IdentifiedHighPt/train/AliAnalysisTaskHighPtDeDx.cxx
PWGLF/SPECTRA/IdentifiedHighPt/train/AliAnalysisTaskHighPtDeDx.h

index 8e2f3cd..0c76fe0 100755 (executable)
@@ -124,6 +124,7 @@ AliAnalysisTaskHighPtDeDx::AliAnalysisTaskHighPtDeDx():
   fVZEROArray(0x0),
   fVtxCut(10.0),  
   fEtaCut(0.9),  
+  fEtaCutStack(1.4),  
   fMinPt(0.1),
   fMinPtV0(0.1),
   fMinCent(0.0),
@@ -181,7 +182,8 @@ AliAnalysisTaskHighPtDeDx::AliAnalysisTaskHighPtDeDx(const char *name):
   fTrackArrayMC(0x0),
   fVZEROArray(0x0),
   fVtxCut(10.0),  
-  fEtaCut(0.9),  
+  fEtaCut(0.9),
+  fEtaCutStack(1.4),    
   fMinPt(0.1),
   fMinPtV0(0.1),
   fMinCent(0.0),
@@ -836,7 +838,7 @@ void AliAnalysisTaskHighPtDeDx::ProcessMCTruthESD()
     Double_t chargeMC = pdgPart->Charge();
 
     
-    if (TMath::Abs(trackMC->Eta()) > 2.4 )
+    if (TMath::Abs(trackMC->Eta()) > fEtaCutStack )
       continue;
     
     trackmult++;
@@ -925,7 +927,7 @@ void AliAnalysisTaskHighPtDeDx::ProcessMCTruthAOD()
     Double_t chargeMC = trackMC->Charge();
 
     
-    if (TMath::Abs(trackMC->Eta()) > 2.4 )
+    if (TMath::Abs(trackMC->Eta()) > fEtaCutStack )
       continue;
     //cout << " debug 34 " << endl;
     trackmult++;
index b723150..757593f 100755 (executable)
@@ -37,6 +37,7 @@ class AliAnalysisTaskHighPtDeDx : public AliAnalysisTaskSE {
   Bool_t   GetAnalysisMC() { return fAnalysisMC; }   
   Double_t GetVtxCut() { return fVtxCut; }   
   Double_t GetEtaCut() { return fEtaCut; }     
+  Double_t GetEtaCutStack() { return fEtaCutStack; }   
   Double_t GetMinPt() { return fMinPt; }   
   Double_t GetMinPtV0() { return fMinPtV0; }
   Int_t    GetTreeOption() { return fTreeOption; }
@@ -52,6 +53,7 @@ class AliAnalysisTaskHighPtDeDx : public AliAnalysisTaskSE {
   virtual void  SetAnalysisMC(Bool_t isMC) {fAnalysisMC = isMC;}
   virtual void  SetVtxCut(Double_t vtxCut){fVtxCut = vtxCut;}
   virtual void  SetEtaCut(Double_t etaCut){fEtaCut = etaCut;}
+  virtual void  SetEtaCutStack(Double_t etaCutStack){fEtaCutStack = etaCutStack;}
   virtual void  SetMinPt(Double_t value) {fMinPt = value;}   
   virtual void  SetMinPtV0(Double_t value) {fMinPtV0 = value;}
   virtual void  SetMinCent(Float_t minvalc) {fMinCent = minvalc;}
@@ -128,8 +130,9 @@ class AliAnalysisTaskHighPtDeDx : public AliAnalysisTaskSE {
   UInt_t       ftrigBit2;
   Double_t     fVtxCut;             // Vtx cut on z position in cm
   Double_t     fEtaCut;             // Eta cut used to select particles
+  Double_t     fEtaCutStack;        // Eta cut used to select particles - reduce saved stack size
   Double_t     fMinPt;              // Min pt - for histogram limits
-  Double_t     fMinPtV0;              // Min pt - for histogram limits
+  Double_t     fMinPtV0;            // Min pt - for histogram limits - V0s / strangeness part
   Double_t     fLowPtFraction;      // Fraction of tracks below min pt to keep
   Double_t     fMassCut;            // Reject all v0 with all dmass > masscut!
   Int_t        fTreeOption;         // 0: no tree, >0: enable debug tree