Modifications to optionally release some track cuts

diff --git a/PWGLF/STRANGENESS/LambdaK0/AliAnalysisTaskExtractPerformanceV0.cxx b/PWGLF/STRANGENESS/LambdaK0/AliAnalysisTaskExtractPerformanceV0.cxx
index bdb46e9ec19729a0e76a09799e7a8f46b6dc78d9..6e1cf5397fec849f6031a7417812e2c8209f57dc 100644 (file)
--- a/PWGLF/STRANGENESS/LambdaK0/AliAnalysisTaskExtractPerformanceV0.cxx
+++ b/PWGLF/STRANGENESS/LambdaK0/AliAnalysisTaskExtractPerformanceV0.cxx
@@ -101,6 +101,88 @@ AliAnalysisTaskExtractPerformanceV0::AliAnalysisTaskExtractPerformanceV0()
    fkLowEnergyPP ( kFALSE ),
    fkUseOnTheFly ( kFALSE ),
 
+//------------------------------------------------
+// Tree Variables 
+
+  fTreeVariablePrimaryStatus(0),
+  fTreeVariablePrimaryStatusMother(0),
+  fTreeVariableChi2V0(0),
+       fTreeVariableDcaV0Daughters(0),
+       fTreeVariableDcaV0ToPrimVertex(0),
+       fTreeVariableDcaPosToPrimVertex(0),
+       fTreeVariableDcaNegToPrimVertex(0),
+       fTreeVariableV0CosineOfPointingAngle(0),
+       fTreeVariableV0Radius(0),
+       fTreeVariablePt(0),
+       fTreeVariablePtMC(0),
+       fTreeVariableRapK0Short(0),
+       fTreeVariableRapLambda(0),
+       fTreeVariableRapMC(0),
+       fTreeVariableInvMassK0s(0),
+       fTreeVariableInvMassLambda(0),
+       fTreeVariableInvMassAntiLambda(0),
+       fTreeVariableAlphaV0(0),
+       fTreeVariablePtArmV0(0),
+       fTreeVariableNegTotMomentum(0),
+       fTreeVariablePosTotMomentum(0),
+       fTreeVariableNegTransvMomentum(0),
+       fTreeVariablePosTransvMomentum(0),
+       fTreeVariableNegTransvMomentumMC(0),
+       fTreeVariablePosTransvMomentumMC(0),
+   
+       fTreeVariableNSigmasPosProton(0),
+       fTreeVariableNSigmasPosPion(0),
+       fTreeVariableNSigmasNegProton(0),
+       fTreeVariableNSigmasNegPion(0),
+
+       fTreeVariablePtMother(0),
+       fTreeVariableV0CreationRadius(0),
+  fTreeVariablePID(0),
+  fTreeVariablePIDPositive(0),
+  fTreeVariablePIDNegative(0),
+  fTreeVariablePIDMother(0),
+  fTreeVariableIndexStatus(0),
+  fTreeVariableIndexStatusMother(0),
+
+  fTreeVariableRunNumber(0),
+  fTreeVariableEventNumber(0),
+
+       fTreeVariableDistOverTotMom(0),
+
+       fTreeVariablePosEta(0),
+       fTreeVariableNegEta(0),
+
+       fTreeVariableVertexZ(0),
+
+  fTreeVariableLeastNbrCrossedRows(0),
+  fTreeVariableLeastRatioCrossedRowsOverFindable(0),
+  fTreeVariableMultiplicity(0),
+  fTreeVariableMultiplicityMC(0),
+
+  fTreeVariableV0x(0),
+  fTreeVariableV0y(0),
+  fTreeVariableV0z(0),
+
+  fTreeVariableV0Px(0),
+  fTreeVariableV0Py(0),
+  fTreeVariableV0Pz(0),
+
+  fTreeVariableMCV0x(0),
+  fTreeVariableMCV0y(0),
+  fTreeVariableMCV0z(0),
+
+  fTreeVariableMCV0Px(0),
+  fTreeVariableMCV0Py(0),
+  fTreeVariableMCV0Pz(0),
+
+  fTreeVariablePVx(0),
+  fTreeVariablePVy(0),
+  fTreeVariablePVz(0),
+
+  fTreeVariableMCPVx(0),
+  fTreeVariableMCPVy(0),
+  fTreeVariableMCPVz(0),
+
 //------------------------------------------------
 // HISTOGRAMS
 // --- Filled on an Event-by-event basis
@@ -185,6 +267,88 @@ AliAnalysisTaskExtractPerformanceV0::AliAnalysisTaskExtractPerformanceV0(const c
    fkLowEnergyPP ( kFALSE ),
    fkUseOnTheFly ( kFALSE ),
      
+//------------------------------------------------
+// Tree Variables 
+
+  fTreeVariablePrimaryStatus(0),
+  fTreeVariablePrimaryStatusMother(0),
+  fTreeVariableChi2V0(0),
+       fTreeVariableDcaV0Daughters(0),
+       fTreeVariableDcaV0ToPrimVertex(0),
+       fTreeVariableDcaPosToPrimVertex(0),
+       fTreeVariableDcaNegToPrimVertex(0),
+       fTreeVariableV0CosineOfPointingAngle(0),
+       fTreeVariableV0Radius(0),
+       fTreeVariablePt(0),
+       fTreeVariablePtMC(0),
+       fTreeVariableRapK0Short(0),
+       fTreeVariableRapLambda(0),
+       fTreeVariableRapMC(0),
+       fTreeVariableInvMassK0s(0),
+       fTreeVariableInvMassLambda(0),
+       fTreeVariableInvMassAntiLambda(0),
+       fTreeVariableAlphaV0(0),
+       fTreeVariablePtArmV0(0),
+       fTreeVariableNegTotMomentum(0),
+       fTreeVariablePosTotMomentum(0),
+       fTreeVariableNegTransvMomentum(0),
+       fTreeVariablePosTransvMomentum(0),
+       fTreeVariableNegTransvMomentumMC(0),
+       fTreeVariablePosTransvMomentumMC(0),
+   
+       fTreeVariableNSigmasPosProton(0),
+       fTreeVariableNSigmasPosPion(0),
+       fTreeVariableNSigmasNegProton(0),
+       fTreeVariableNSigmasNegPion(0),
+
+       fTreeVariablePtMother(0),
+       fTreeVariableV0CreationRadius(0),
+  fTreeVariablePID(0),
+  fTreeVariablePIDPositive(0),
+  fTreeVariablePIDNegative(0),
+  fTreeVariablePIDMother(0),
+  fTreeVariableIndexStatus(0),
+  fTreeVariableIndexStatusMother(0),
+
+  fTreeVariableRunNumber(0),
+  fTreeVariableEventNumber(0),
+
+       fTreeVariableDistOverTotMom(0),
+
+       fTreeVariablePosEta(0),
+       fTreeVariableNegEta(0),
+
+       fTreeVariableVertexZ(0),
+
+  fTreeVariableLeastNbrCrossedRows(0),
+  fTreeVariableLeastRatioCrossedRowsOverFindable(0),
+  fTreeVariableMultiplicity(0),
+  fTreeVariableMultiplicityMC(0),
+
+  fTreeVariableV0x(0),
+  fTreeVariableV0y(0),
+  fTreeVariableV0z(0),
+
+  fTreeVariableV0Px(0),
+  fTreeVariableV0Py(0),
+  fTreeVariableV0Pz(0),
+
+  fTreeVariableMCV0x(0),
+  fTreeVariableMCV0y(0),
+  fTreeVariableMCV0z(0),
+
+  fTreeVariableMCV0Px(0),
+  fTreeVariableMCV0Py(0),
+  fTreeVariableMCV0Pz(0),
+
+  fTreeVariablePVx(0),
+  fTreeVariablePVy(0),
+  fTreeVariablePVz(0),
+
+  fTreeVariableMCPVx(0),
+  fTreeVariableMCPVy(0),
+  fTreeVariableMCPVz(0),
+
 //------------------------------------------------
 // HISTOGRAMS
 // --- Filled on an Event-by-event basis

diff --git a/PWGLF/STRANGENESS/LambdaK0/AliAnalysisTaskExtractV0.cxx b/PWGLF/STRANGENESS/LambdaK0/AliAnalysisTaskExtractV0.cxx
index 6c90c0f3727a0f1d048c7208f8116dcd67811bb5..200cc8984e810417ecced782aa39a082fd322027 100644 (file)
--- a/PWGLF/STRANGENESS/LambdaK0/AliAnalysisTaskExtractV0.cxx
+++ b/PWGLF/STRANGENESS/LambdaK0/AliAnalysisTaskExtractV0.cxx
@@ -91,9 +91,58 @@ ClassImp(AliAnalysisTaskExtractV0)
 
 AliAnalysisTaskExtractV0::AliAnalysisTaskExtractV0() 
   : AliAnalysisTaskSE(), fListHistV0(0), fTree(0), fPIDResponse(0),fESDtrackCuts(0),
-   fkIsNuclear   ( kFALSE ), 
-   fkLowEnergyPP ( kFALSE ),
-   fkUseOnTheFly ( kFALSE ),
+   fkIsNuclear     ( kFALSE ), 
+   fkLowEnergyPP   ( kFALSE ),
+   fkUseOnTheFly   ( kFALSE ),
+   fkTakeAllTracks ( kFALSE ),
+//------------------------------------------------
+// Initialize 
+       fTreeVariableChi2V0(0),
+       fTreeVariableDcaV0Daughters(0),
+       fTreeVariableDcaV0ToPrimVertex(0),
+       fTreeVariableDcaPosToPrimVertex(0),
+       fTreeVariableDcaNegToPrimVertex(0),
+       fTreeVariableV0CosineOfPointingAngle(0),
+       fTreeVariableV0Radius(0),
+       fTreeVariablePt(0),
+       fTreeVariableRapK0Short(0),
+       fTreeVariableRapLambda(0),
+       fTreeVariableInvMassK0s(0),
+       fTreeVariableInvMassLambda(0),
+       fTreeVariableInvMassAntiLambda(0),
+       fTreeVariableAlphaV0(0),
+       fTreeVariablePtArmV0(0),
+       fTreeVariableNegTotMomentum(0),
+       fTreeVariablePosTotMomentum(0),
+       fTreeVariableNegdEdxSig(0),
+       fTreeVariablePosdEdxSig(0),
+       fTreeVariableNegEta(0),
+       fTreeVariablePosEta(0),
+
+       fTreeVariableNSigmasPosProton(0),
+       fTreeVariableNSigmasPosPion(0),
+       fTreeVariableNSigmasNegProton(0),
+       fTreeVariableNSigmasNegPion(0),
+       
+       fTreeVariableDistOverTotMom(0),
+       fTreeVariableLeastNbrCrossedRows(0),
+       fTreeVariableLeastRatioCrossedRowsOverFindable(0),
+       fTreeVariableMultiplicity(0),
+  
+  fTreeVariableRunNumber(0),
+  fTreeVariableEventNumber(0),
+  
+  fTreeVariableV0x(0),
+  fTreeVariableV0y(0),
+  fTreeVariableV0z(0),
+
+  fTreeVariableV0Px(0),
+  fTreeVariableV0Py(0),
+  fTreeVariableV0Pz(0),
+
+  fTreeVariablePVx(0),
+  fTreeVariablePVy(0),
+  fTreeVariablePVz(0),
 
 //------------------------------------------------
 // HISTOGRAMS
@@ -130,10 +179,60 @@ AliAnalysisTaskExtractV0::AliAnalysisTaskExtractV0()
 
 AliAnalysisTaskExtractV0::AliAnalysisTaskExtractV0(const char *name) 
   : AliAnalysisTaskSE(name), fListHistV0(0), fTree(0), fPIDResponse(0),fESDtrackCuts(0),
-   fkIsNuclear   ( kFALSE ), 
-   fkLowEnergyPP ( kFALSE ),
-   fkUseOnTheFly ( kFALSE ),
+   fkIsNuclear     ( kFALSE ), 
+   fkLowEnergyPP   ( kFALSE ),
+   fkUseOnTheFly   ( kFALSE ),
+   fkTakeAllTracks ( kFALSE ),
      
+//------------------------------------------------
+// Initialize 
+       fTreeVariableChi2V0(0),
+       fTreeVariableDcaV0Daughters(0),
+       fTreeVariableDcaV0ToPrimVertex(0),
+       fTreeVariableDcaPosToPrimVertex(0),
+       fTreeVariableDcaNegToPrimVertex(0),
+       fTreeVariableV0CosineOfPointingAngle(0),
+       fTreeVariableV0Radius(0),
+       fTreeVariablePt(0),
+       fTreeVariableRapK0Short(0),
+       fTreeVariableRapLambda(0),
+       fTreeVariableInvMassK0s(0),
+       fTreeVariableInvMassLambda(0),
+       fTreeVariableInvMassAntiLambda(0),
+       fTreeVariableAlphaV0(0),
+       fTreeVariablePtArmV0(0),
+       fTreeVariableNegTotMomentum(0),
+       fTreeVariablePosTotMomentum(0),
+       fTreeVariableNegdEdxSig(0),
+       fTreeVariablePosdEdxSig(0),
+       fTreeVariableNegEta(0),
+       fTreeVariablePosEta(0),
+
+       fTreeVariableNSigmasPosProton(0),
+       fTreeVariableNSigmasPosPion(0),
+       fTreeVariableNSigmasNegProton(0),
+       fTreeVariableNSigmasNegPion(0),
+       
+       fTreeVariableDistOverTotMom(0),
+       fTreeVariableLeastNbrCrossedRows(0),
+       fTreeVariableLeastRatioCrossedRowsOverFindable(0),
+       fTreeVariableMultiplicity(0),
+  
+  fTreeVariableRunNumber(0),
+  fTreeVariableEventNumber(0),
+  
+  fTreeVariableV0x(0),
+  fTreeVariableV0y(0),
+  fTreeVariableV0z(0),
+
+  fTreeVariableV0Px(0),
+  fTreeVariableV0Py(0),
+  fTreeVariableV0Pz(0),
+
+  fTreeVariablePVx(0),
+  fTreeVariablePVy(0),
+  fTreeVariablePVz(0),
+
 //------------------------------------------------
 // HISTOGRAMS
 // --- Filled on an Event-by-event basis
@@ -700,7 +799,7 @@ void AliAnalysisTaskExtractV0::UserExec(Option_t *)
       if( !(pTrack->GetStatus() & AliESDtrack::kTPCrefit)) continue;
       if( !(nTrack->GetStatus() & AliESDtrack::kTPCrefit)) continue;
 
-      if ( ( ( pTrack->GetTPCClusterInfo(2,1) ) < 70 ) || ( ( nTrack->GetTPCClusterInfo(2,1) ) < 70 ) ) continue;
+      if ( ( ( ( pTrack->GetTPCClusterInfo(2,1) ) < 70 ) || ( ( nTrack->GetTPCClusterInfo(2,1) ) < 70 ) )&&(fkTakeAllTracks==kFALSE) ) continue;
        
       //GetKinkIndex condition
       if( pTrack->GetKinkIndex(0)>0 || nTrack->GetKinkIndex(0)>0 ) continue;
@@ -718,7 +817,7 @@ void AliAnalysisTaskExtractV0::UserExec(Option_t *)
          fTreeVariableLeastRatioCrossedRowsOverFindable = lNegTrackCrossedRowsOverFindable;
 
       //Lowest Cut Level for Ratio Crossed Rows / Findable = 0.8, set here
-      if ( fTreeVariableLeastRatioCrossedRowsOverFindable < 0.8 ) continue;
+      if ( (fTreeVariableLeastRatioCrossedRowsOverFindable < 0.8)&&(fkTakeAllTracks==kFALSE) ) continue;
 
       //End track Quality Cuts
       //________________________________________________________________________

diff --git a/PWGLF/STRANGENESS/LambdaK0/AliAnalysisTaskExtractV0.h b/PWGLF/STRANGENESS/LambdaK0/AliAnalysisTaskExtractV0.h
index 74500b1de20303706a71fe810d37058d36984526..bafb8b415b623f59057f5d85fa5218be471a757e 100644 (file)
--- a/PWGLF/STRANGENESS/LambdaK0/AliAnalysisTaskExtractV0.h
+++ b/PWGLF/STRANGENESS/LambdaK0/AliAnalysisTaskExtractV0.h
@@ -61,6 +61,7 @@ class AliAnalysisTaskExtractV0 : public AliAnalysisTaskSE {
   void SetIsNuclear           (Bool_t lIsNuclear   = kTRUE ) { fkIsNuclear   = lIsNuclear;   }
   void SetIsLowEnergyPP       (Bool_t lLowEnergyPP = kTRUE ) { fkLowEnergyPP = lLowEnergyPP; }
   void SetUseOnTheFly         (Bool_t lUseOnTheFly = kTRUE ) { fkUseOnTheFly = lUseOnTheFly; }
+  void SetTakeAllTracks       (Bool_t lTakeAllTracks = kTRUE ) { fkTakeAllTracks = lTakeAllTracks; }
 
  private:
                                // Note : In ROOT, "//!" means "do not stream the data from Master node to Worker node" ...
@@ -77,6 +78,7 @@ class AliAnalysisTaskExtractV0 : public AliAnalysisTaskSE {
   Bool_t fkIsNuclear;   // if true, replace multiplicity est. by centrality (default FALSE) 
   Bool_t fkLowEnergyPP; // if true, skip FASTOnly (default FALSE)
   Bool_t fkUseOnTheFly; // if true, will use On-the-fly V0s instead of Offline V0s (default FALSE)
+  Bool_t fkTakeAllTracks; // if true, no TPC crossed rows and ratio cut
 
   //Variables for Tree
        Float_t fTreeVariableChi2V0;         //!