coverity fix,
authorpkalinak <pkalinak@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 15 Apr 2013 10:36:59 +0000 (10:36 +0000)
committerpkalinak <pkalinak@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 15 Apr 2013 10:36:59 +0000 (10:36 +0000)
cleaning of the task,
histos for Xi and Omegas

PWGLF/STRANGENESS/LambdaK0PbPb/AliAnalysisTaskPerformanceStrange.cxx

index 949bf1b..a4613c2 100644 (file)
@@ -43,6 +43,7 @@
 #include "AliESDInputHandler.h"
 #include "AliESDtrack.h"
 #include "AliESDv0.h"
+#include "AliESDcascade.h"
 #include "AliESDtrackCuts.h"
 #include "AliESDpid.h"
 #include "AliMultiplicity.h"
@@ -82,779 +83,185 @@ ClassImp(AliAnalysisTaskPerformanceStrange)
 AliAnalysisTaskPerformanceStrange::AliAnalysisTaskPerformanceStrange()
 : AliAnalysisTaskSE(), fAnalysisMC(0), fAnalysisType("infoType"),  fCollidingSystems(0), fUsePID("infoPID"), fUseCut("infoCut"),fDown(0),fUp(0), fESD(0), fListHist(0),fCentrSelector(0),fTracksCuts(0),fPIDResponse(0),fQASelector(0), 
 
-  fHistMCPrimaryVertexX(0),
-  fHistMCPrimaryVertexY(0),
-  fHistMCPrimaryVertexZ(0),
-  fHistPtTracksITSRefit(0),
-  fHistPtTracks(0),
-  fHistPtTracksPITSRefit(0),
-  fHistPtTracksP(0),
-  fHistMCMultiplicityPrimary(0),
-  fHistMCMultiplicityTracks(0),
-  fHistTPCTracks(0),
-  fHistMCtracksProdRadiusK0s(0),
-  fHistMCtracksProdRadiusLambda(0),
-  fHistMCtracksProdRadiusAntiLambda(0),
-  fHistMCtracksDecayRadiusK0s(0),
-  fHistMCtracksDecayRadiusLambda(0),
-  fHistMCtracksDecayRadiusAntiLambda(0),
-
-  fHistMCPtAllK0s(0),
-  fHistMCPtAllLambda(0),
-  fHistMCPtAllAntiLambda(0),
-
-//Rap3
-  fHistMCPtAllK0sRap3(0),
-  fHistMCPtAllLambdaRap3(0),
-  fHistMCPtAllAntiLambdaRap3(0),
-
-
-  fHistMCProdRadiusK0s(0),
-  fHistMCProdRadiusLambda(0),
-  fHistMCProdRadiusAntiLambda(0),
-
-//Decay radius of primary V0s
-  fHistMCPrimDecayRadiusK0s(0),
-  fHistMCPrimDecayRadiusLambda(0),
-  fHistMCPrimDecayRadiusAntiLambda(0),
-
-  fHistMCRapK0s(0),
-  fHistMCRapInPtRangeK0s(0),
-  fHistMCRapLambda(0),
-  fHistMCRapInPtRangeLambda(0),
-  fHistMCRapAntiLambda(0),
-  fHistMCRapInPtRangeAntiLambda(0),
-  fHistMCRapXi(0),
-  fHistMCRapInPtRangeXi(0),
-  fHistMCRapPhi(0),
-  fHistMCRapInPtRangePhi(0),
-////////////////////////////////////////
-
-  fHistMCPtK0s(0),
-  fHistMCPtLambda(0),
-  fHistMCPtAntiLambda(0),
-
-//Rap3
-  fHistMCPtK0sRap3(0),
-  fHistMCPtLambdaRap3(0),
-  fHistMCPtAntiLambdaRap3(0),
-
-
-///////////////////////////////////////////
-
-  fHistMCPtLambdaFromSigma(0),
-  fHistMCPtAntiLambdaFromSigma(0),
-  fHistNTimesRecK0s(0),
-  fHistNTimesRecLambda(0),
-  fHistNTimesRecAntiLambda(0),
-  fHistNTimesRecK0sVsPt(0),
-  fHistNTimesRecLambdaVsPt(0),
-  fHistNTimesRecAntiLambdaVsPt(0),
-  fHistNumberEvents(0),
-  fHistTrackPerEvent(0),
-  fHistTPCMult(0),
-  fHistTrackletPerEvent(0),
-  fHistMCDaughterTrack(0),
-  fHistSPDPrimaryVertexZ(0),
-  fHistPrimaryVertexX(0),
-  fHistPrimaryVertexY(0),
-  fHistPrimaryVertexZ(0),
-  fHistPrimaryVertexResX(0),
-  fHistPrimaryVertexResY(0),
-  fHistPrimaryVertexResZ(0),
-  fHistPrimaryVertexPosXV0events(0), 
-  fHistPrimaryVertexPosYV0events(0), 
-  fHistPrimaryVertexPosZV0events(0),
-  fHistDaughterPt(0),
-////////////////////////////////////////
-
-  fHistDcaPosToPrimVertexK0(0),
-  fHistDcaNegToPrimVertexK0(0),
-  fHistRadiusV0K0(0),
-  fHistDecayLengthV0K0(0),
-  fHistDcaV0DaughtersK0(0),
-  fHistChi2K0(0),
-  fHistCosPointAngleK0(0),
-
-
-  fHistDcaPosToPrimVertexK0vsMassK0(0),
-  fHistDcaNegToPrimVertexK0vsMassK0(0),
-  fHistRadiusV0K0vsMassK0(0),
-  fHistDecayLengthV0K0vsMassK0(0),
-  fHistDcaV0DaughtersK0vsMassK0(0),
-  fHistCosPointAngleK0vsMassK0(0),
-  fHistDcaPosToPrimVertexK0vsMassK0pt1(0),
-  fHistDcaNegToPrimVertexK0vsMassK0pt1(0),
-  fHistRadiusV0K0vsMassK0pt1(0),
-  fHistDecayLengthV0K0vsMassK0pt1(0),
-  fHistDcaV0DaughtersK0vsMassK0pt1(0),
-  fHistCosPointAngleK0vsMassK0pt1(0),
-
-  fHistDcaPosToPrimVertexK0vsMassK0pt2(0),
-  fHistDcaNegToPrimVertexK0vsMassK0pt2(0),
-  fHistRadiusV0K0vsMassK0pt2(0),
-  fHistDecayLengthV0K0vsMassK0pt2(0),
-  fHistDcaV0DaughtersK0vsMassK0pt2(0),
-  fHistCosPointAngleK0vsMassK0pt2(0),
-
-  fHistDcaPosToPrimVertexK0vsMassK0pt3(0),
-  fHistDcaNegToPrimVertexK0vsMassK0pt3(0),
-  fHistRadiusV0K0vsMassK0pt3(0),
-  fHistDecayLengthV0K0vsMassK0pt3(0),
-  fHistDcaV0DaughtersK0vsMassK0pt3(0),
-  fHistCosPointAngleK0vsMassK0pt3(0),
-
-/////////// Lambda ///////////////////////////
-
-  fHistDcaPosToPrimVertexL(0),
-  fHistDcaNegToPrimVertexL(0),
-  fHistRadiusV0L(0),
-  fHistDecayLengthV0L(0),
-  fHistDcaV0DaughtersL(0),
-  fHistChi2L(0),
-  fHistCosPointAngleL(0),
-
-
-  fHistDcaPosToPrimVertexLvsMassL(0),
-  fHistDcaNegToPrimVertexLvsMassL(0),
-  fHistRadiusV0LvsMassL(0),
-  fHistDecayLengthV0LvsMassL(0),
-  fHistDcaV0DaughtersLvsMassL(0),
-  fHistCosPointAngleLvsMassL(0),
-  fHistCosPointAngleLvsMassVsPtsigL(0),
-  fHistCosPointAngleLvsMassVsPtbackL(0),
 
-    
-  fHistDcaPosToPrimVertexLambdaVsMasspt1(0),
-  fHistDcaNegToPrimVertexLambdaVsMasspt1(0),
-  fHistRadiusV0LambdaVsMasspt1(0),
-  fHistDecayLengthV0LambdaVsMasspt1(0),
-  fHistDcaV0DaughtersLambdaVsMasspt1(0),
-  fHistCosPointAngleLambdaVsMasspt1(0),
-
-  fHistDcaPosToPrimVertexLambdaVsMasspt2(0),
-  fHistDcaNegToPrimVertexLambdaVsMasspt2(0),
-  fHistRadiusV0LambdaVsMasspt2(0),
-  fHistDecayLengthV0LambdaVsMasspt2(0),
-  fHistDcaV0DaughtersLambdaVsMasspt2(0),
-  fHistCosPointAngleLambdaVsMasspt2(0),
-
-  fHistDcaPosToPrimVertexLambdaVsMasspt3(0),
-  fHistDcaNegToPrimVertexLambdaVsMasspt3(0),
-  fHistRadiusV0LambdaVsMasspt3(0),
-  fHistDecayLengthV0LambdaVsMasspt3(0),
-  fHistDcaV0DaughtersLambdaVsMasspt3(0),
-  fHistCosPointAngleLambdaVsMasspt3(0),
-
-
-
-/////////Antilambda ///////////////////
-  fHistDcaPosToPrimVertexAntiL(0),
-  fHistDcaNegToPrimVertexAntiL(0),
-  fHistRadiusV0AntiL(0),
-  fHistDecayLengthV0AntiL(0),
-  fHistDcaV0DaughtersAntiL(0),
-  fHistChi2AntiL(0),
-  fHistCosPointAngleAntiL(0),
-
-
-  fHistDcaPosToPrimVertexAntiLvsMass(0),
-  fHistDcaNegToPrimVertexAntiLvsMass(0),
-  fHistRadiusV0AntiLvsMass(0),
-  fHistDecayLengthV0AntiLvsMass(0),
-  fHistDcaV0DaughtersAntiLvsMass(0),
-  fHistCosPointAngleAntiLvsMass(0),
+    /////// primary vertex
+    fHistMCPrimaryVertexX(0),   fHistMCPrimaryVertexY(0),   fHistMCPrimaryVertexZ(0),
+    ////// tracks & multiplicity
+    fHistPtTracks(0), 
+    fHistMCMultiplicityPrimary(0),  fHistMCMultiplicityTracks(0),  fHistTPCTracks(0), 
+    ///// Transverse Momentum
+    fHistMCPtAllK0s(0),
+    fHistMCPtAllLambda(0),   fHistMCPtAllAntiLambda(0),
+    fHistMCPtAllXi(0),  fHistMCPtAllAntiXi(0),
+    fHistMCPtAllOmega(0), fHistMCPtAllAntiOmega(0),
+    /// Rapidity
+    fHistMCRapK0s(0),
+    fHistMCRapLambda(0), fHistMCRapAntiLambda(0),
+    fHistMCRapXi(0), 
 
-    
-    
-  fHistDcaPosToPrimVertexAntiLVsMasspt1(0),
-  fHistDcaNegToPrimVertexAntiLVsMasspt1(0),
-  fHistRadiusV0AntiLVsMasspt1(0),
-  fHistDecayLengthV0AntiLVsMasspt1(0),
-  fHistDcaV0DaughtersAntiLVsMasspt1(0),
-  fHistCosPointAngleAntiLVsMasspt1(0),
-
-  fHistDcaPosToPrimVertexAntiLVsMasspt2(0),
-  fHistDcaNegToPrimVertexAntiLVsMasspt2(0),
-  fHistRadiusV0AntiLVsMasspt2(0),
-  fHistDecayLengthV0AntiLVsMasspt2(0),
-  fHistDcaV0DaughtersAntiLVsMasspt2(0),
-  fHistCosPointAngleAntiLVsMasspt2(0),
-
-  fHistDcaPosToPrimVertexAntiLVsMasspt3(0),
-  fHistDcaNegToPrimVertexAntiLVsMasspt3(0),
-  fHistRadiusV0AntiLVsMasspt3(0),
-  fHistDecayLengthV0AntiLVsMasspt3(0),
-  fHistDcaV0DaughtersAntiLVsMasspt3(0),
-  fHistCosPointAngleAntiLVsMasspt3(0),
-
-/////////////////////////////////////////
-  fHistV0Multiplicity(0),
-  fHistMassK0(0),
-  fHistMassLambda(0),
-  fHistMassAntiLambda(0),
-  fHistMassVsRadiusK0(0),
-  fHistMassVsRadiusLambda(0),
-  fHistMassVsRadiusAntiLambda(0),
-
-///////////////////////////////////////
-  fHistPtVsMassK0(0),
-  fHistPtVsMassLambda(0),
-  fHistPtVsMassAntiLambda(0),
-
-/// Rap3
-  fHistPtVsMassK0Rap3(0),
-  fHistPtVsMassLambdaRap3(0),
-  fHistPtVsMassAntiLambdaRap3(0),
-
-    // cTau dist
-    fHistTranscTauVsMassL(0),
-    fHistTranscTauVsMassAntiL(0),
-    fHistTranscTauVsMassK0s(0),
-
-    // cTau Trans Vs. Mass Rap3
-    fHistTranscTauVsMassLRap3(0),
-    fHistTranscTauVsMassAntiLRap3(0),
-    fHistTranscTauVsMassK0sRap3(0),
-
-    // cTau Trans Vs Mass Low pt
-    fHistTranscTauVsMassLptLow(0),
-    fHistTranscTauVsMassAntiLptLow(0),
-    fHistTranscTauVsMassK0sptLow(0),
-
-    // cTau Trans Vs Mass Low pt Rap3
-    fHistTranscTauVsMassLptLowRap3(0),
-    fHistTranscTauVsMassAntiLptLowRap3(0),
-    fHistTranscTauVsMassK0sptLowRap3(0),
-
-
-////////////////////////////////////////
-
-  fHistArmenterosPodolanski(0),
-  fHistK0sMassVsLambdaMass(0),
-  fHistTPCsigPLambda(0),
-  fHistTPCsigPAntiLambda(0),
-  fHistNSigmaProton(0),    
-  fHistNsigmaPosPionAntiLambda(0),
-  fHistNsigmaNegProtonAntiLambda(0),
-  fHistNsigmaPosProtonLambda(0),
-  fHistNsigmaNegPionLambda(0),
-  fHistNsigmaPosProtonAntiLambda(0),
-  fHistNsigmaNegPionAntiLambda(0),
-  fHistNsigmaPosPionK0(0),
-  fHistNsigmaNegPionK0(0),
-  fHistAsMcRapK0(0),
-  fHistAsMcRapLambda(0),
-  fHistAsMcRapAntiLambda(0),
-
-// pt distribution    /////////////////////
-  fHistAsMcPtK0(0),
-  fHistAsMcPtLambda(0),
-  fHistAsMcPtAntiLambda(0),
-
-
-// Transverse cTau distribution
-  fHistAsMcTranscTauL(0),
-  fHistAsMcTranscTauAntiL(0),
-  fHistAsMcTranscTauK0s(0),
-
-    fHistAsMcTranscTauLptLow(0),
-    fHistAsMcTranscTauAntiLptLow(0),
-    fHistAsMcTranscTauK0sptLow(0),
-
-//Rap3
-    fHistAsMcTranscTauLptLowRap3(0),
-    fHistAsMcTranscTauAntiLptLowRap3(0),
-    fHistAsMcTranscTauK0sptLowRap3(0),
-
-//Rap3
-  fHistAsMcPtK0Rap3(0),
-  fHistAsMcPtLambdaRap3(0),
-  fHistAsMcPtAntiLambdaRap3(0),
-
-
-/////////////////////////////////////
-
-  fHistAsMcPtZoomK0(0),
-  fHistAsMcPtZoomLambda(0),
-  fHistAsMcPtZoomAntiLambda(0),
-  fHistAsMcProdRadiusK0(0),
-  fHistAsMcProdRadiusLambda(0),
-  fHistAsMcProdRadiusAntiLambda(0),
-  fHistAsMcProdRadiusXvsYK0s(0),
-  fHistAsMcProdRadiusXvsYLambda(0),
-  fHistAsMcProdRadiusXvsYAntiLambda(0),
-  fHistPidMcMassK0(0),
-  fHistPidMcMassLambda(0),
-  fHistPidMcMassAntiLambda(0),
-
-
-///////////////////////////////
-
-  fHistAsMcMassK0(0),
-  fHistAsMcMassLambda(0),
-  fHistAsMcMassAntiLambda(0),
-
-/// Rap3
-  fHistAsMcMassK0Rap3(0),
-  fHistAsMcMassLambdaRap3(0),
-  fHistAsMcMassAntiLambdaRap3(0),
-
-
-  fHistAsMcPtVsMassK0(0),
-  fHistAsMcPtVsMassLambda(0),
-  fHistAsMcPtVsMassAntiLambda(0),
-
-/// Rap3
-  fHistAsMcPtVsMassK0Rap3(0),
-  fHistAsMcPtVsMassLambdaRap3(0),
-  fHistAsMcPtVsMassAntiLambdaRap3(0),
-
-///////////////////////////////////////
-
-  fHistAsMcMassVsRadiusK0(0),
-  fHistAsMcMassVsRadiusLambda(0),
-  fHistAsMcMassVsRadiusAntiLambda(0),
-  fHistAsMcResxK0(0),
-  fHistAsMcResyK0(0),
-  fHistAsMcReszK0(0),
-  fHistAsMcResrVsRadiusK0(0),
-  fHistAsMcReszVsRadiusK0(0),
-  fHistAsMcResxLambda(0),
-  fHistAsMcResyLambda(0),
-  fHistAsMcReszLambda(0),
-  fHistAsMcResrVsRadiusLambda(0),
-  fHistAsMcReszVsRadiusLambda(0),
-  fHistAsMcResxAntiLambda(0),
-  fHistAsMcResyAntiLambda(0),
-  fHistAsMcReszAntiLambda(0),
-  fHistAsMcResrVsRadiusAntiLambda(0),
-  fHistAsMcReszVsRadiusAntiLambda(0),
-  fHistAsMcResPtK0(0),
-  fHistAsMcResPtLambda(0),
-  fHistAsMcResPtAntiLambda(0),
-  fHistAsMcResPtVsRapK0(0),
-  fHistAsMcResPtVsRapLambda(0),
-  fHistAsMcResPtVsRapAntiLambda(0),
-  fHistAsMcResPtVsPtK0(0),
-  fHistAsMcResPtVsPtLambda(0),
-  fHistAsMcResPtVsPtAntiLambda(0),
-  fHistAsMcMotherPdgCodeK0s(0),
-  fHistAsMcMotherPdgCodeLambda(0),
-  fHistAsMcMotherPdgCodeAntiLambda(0),
-  fHistAsMcPtLambdaFromSigma(0),
-  fHistAsMcPtAntiLambdaFromSigma(0),
-  fHistAsMcSecondaryPtVsRapK0s(0),
-  fHistAsMcSecondaryPtVsRapLambda(0),
-  fHistAsMcSecondaryPtVsRapAntiLambda(0),
-  fHistAsMcSecondaryProdRadiusK0s(0),
-  fHistAsMcSecondaryProdRadiusLambda(0),
-  fHistAsMcSecondaryProdRadiusAntiLambda(0),
-  fHistAsMcSecondaryProdRadiusXvsYK0s(0),
-  fHistAsMcSecondaryProdRadiusXvsYLambda(0),
-  fHistAsMcSecondaryProdRadiusXvsYAntiLambda(0),
-  fHistAsMcSecondaryMotherPdgCodeK0s(0),
-  fHistAsMcSecondaryMotherPdgCodeLambda(0),
-  fHistAsMcSecondaryMotherPdgCodeAntiLambda(0),
-  fHistAsMcSecondaryPtLambdaFromSigma(0),
-  fHistAsMcSecondaryPtAntiLambdaFromSigma(0)
+    ///// Transverse Momentum primary
+    fHistMCPtK0s(0),
+    fHistMCPtLambda(0),  fHistMCPtAntiLambda(0),
+
+    ///////////////////////////////////////////
+    ////   ESD
+
+    fHistNumberEvents(0),
+    fHistTrackPerEvent(0),
+    fHistTPCMult(0),
+    fHistTrackletPerEvent(0),
+    fHistSPDPrimaryVertexZ(0),
+    fHistPrimaryVertexX(0),
+    fHistPrimaryVertexY(0),
+    fHistPrimaryVertexZ(0),
+
+    fHistV0Multiplicity(0),
+    ///// inv. mass 
+    fHistMassK0(0),
+    fHistMassLambda(0),
+    fHistMassAntiLambda(0),
+    fHistMassXi(0),
+    fHistMassAntiXi(0),
+    fHistMassOmega(0),
+    fHistMassAntiOmega(0),
+
+    //inv mass vs PID
+    fHistMassXiVsPID(0),
  
-{
-  // Constructor
-}
+    ///////////////////////////////////////
+    fHistPtVsMassK0(0),
+    fHistPtVsMassLambda(0),
+    fHistPtVsMassAntiLambda(0),
 
+    ////////////////////////////////////////
 
+    fHistArmenterosPodolanski(0),
+    fHistK0sMassVsLambdaMass(0),
+    fHistTPCsigPLambda(0),
+    fHistTPCsigPAntiLambda(0),
+    fHistNSigmaProton(0),  
+  
+    /// Associated histos 
+    ///rapidity
+    fHistAsMcRapK0(0),
+    fHistAsMcRapLambda(0),   fHistAsMcRapAntiLambda(0),
 
+    // pt distribution    /////////////////////
+    fHistAsMcPtK0(0),
+    fHistAsMcPtLambda(0),   fHistAsMcPtAntiLambda(0),
 
+    /////////////////////////////////////
+    fHistPidMcMassK0(0),
+    fHistPidMcMassLambda(0),
+    fHistPidMcMassAntiLambda(0),
+
+    ///inv. mass
+    fHistAsMcMassK0(0),
+    fHistAsMcMassLambda(0), fHistAsMcMassAntiLambda(0),
+    ///transv. momentum
+    fHistAsMcPtVsMassK0(0),
+    fHistAsMcPtVsMassLambda(0), fHistAsMcPtVsMassAntiLambda(0),
+
+    /////background composition
+    fHistCompositionXi(0),
+    fHistCompositionAntiXi(0),
+    fHistCompositionOmega(0),
+  fHistCompositionAntiOmega(0),
+
+  fHistMCIndexes(0)
+{
+  // Constructor
+}
 
 //________________________________________________________________________
 AliAnalysisTaskPerformanceStrange::AliAnalysisTaskPerformanceStrange(const char *name)
   : AliAnalysisTaskSE(name), fAnalysisMC(0), fAnalysisType("infoType"), fCollidingSystems(0), fUsePID("infoPID"), fUseCut("infocut"),fDown(0),fUp(0), fESD(0), fListHist(),fCentrSelector(0), fTracksCuts(0),fPIDResponse(0),fQASelector(0),
 
-    fHistMCPrimaryVertexX(0),
-    fHistMCPrimaryVertexY(0),
-    fHistMCPrimaryVertexZ(0),
-    fHistPtTracksITSRefit(0),
-    fHistPtTracks(0),
-    fHistPtTracksPITSRefit(0),
-    fHistPtTracksP(0),
-    fHistMCMultiplicityPrimary(0),
-    fHistMCMultiplicityTracks(0),
-    fHistTPCTracks(0),
-    fHistMCtracksProdRadiusK0s(0),
-    fHistMCtracksProdRadiusLambda(0),
-    fHistMCtracksProdRadiusAntiLambda(0),
-    fHistMCtracksDecayRadiusK0s(0),
-    fHistMCtracksDecayRadiusLambda(0),
-    fHistMCtracksDecayRadiusAntiLambda(0),
-    // pt distribution
+    /////// primary vertex
+    fHistMCPrimaryVertexX(0),   fHistMCPrimaryVertexY(0),   fHistMCPrimaryVertexZ(0),
+    ////// tracks & multiplicity
+    fHistPtTracks(0), 
+    fHistMCMultiplicityPrimary(0),  fHistMCMultiplicityTracks(0),  fHistTPCTracks(0), 
+    ///// Transverse Momentum
     fHistMCPtAllK0s(0),
-    fHistMCPtAllLambda(0),
-    fHistMCPtAllAntiLambda(0),
-
-    //Rap3
-    fHistMCPtAllK0sRap3(0),
-    fHistMCPtAllLambdaRap3(0),
-    fHistMCPtAllAntiLambdaRap3(0),
-
-
-    fHistMCProdRadiusK0s(0),
-    fHistMCProdRadiusLambda(0),
-    fHistMCProdRadiusAntiLambda(0),
-  
-    // Decay radius of Primary V0s
-    fHistMCPrimDecayRadiusK0s(0),
-    fHistMCPrimDecayRadiusLambda(0),
-    fHistMCPrimDecayRadiusAntiLambda(0),
-
+    fHistMCPtAllLambda(0),   fHistMCPtAllAntiLambda(0),
+    fHistMCPtAllXi(0),  fHistMCPtAllAntiXi(0),
+    fHistMCPtAllOmega(0), fHistMCPtAllAntiOmega(0),
+    /// Rapidity
     fHistMCRapK0s(0),
-    fHistMCRapInPtRangeK0s(0),
-    fHistMCRapLambda(0),
-    fHistMCRapInPtRangeLambda(0),
-    fHistMCRapAntiLambda(0),
-    fHistMCRapInPtRangeAntiLambda(0),
-    fHistMCRapXi(0),
-    fHistMCRapInPtRangeXi(0),
-    fHistMCRapPhi(0),
-    fHistMCRapInPtRangePhi(0),
-    ////////////////////////////////////////
+    fHistMCRapLambda(0), fHistMCRapAntiLambda(0),
+    fHistMCRapXi(0), 
 
+    ///// Transverse Momentum primary
     fHistMCPtK0s(0),
-    fHistMCPtLambda(0),
-    fHistMCPtAntiLambda(0),
-
-    //Rap3
-    fHistMCPtK0sRap3(0),
-    fHistMCPtLambdaRap3(0),
-    fHistMCPtAntiLambdaRap3(0),
-
+    fHistMCPtLambda(0),  fHistMCPtAntiLambda(0),
 
     ///////////////////////////////////////////
+    ////   ESD
 
-    fHistMCPtLambdaFromSigma(0),
-    fHistMCPtAntiLambdaFromSigma(0),
-    fHistNTimesRecK0s(0),
-    fHistNTimesRecLambda(0),
-    fHistNTimesRecAntiLambda(0),
-    fHistNTimesRecK0sVsPt(0),
-    fHistNTimesRecLambdaVsPt(0),
-    fHistNTimesRecAntiLambdaVsPt(0),
     fHistNumberEvents(0),
     fHistTrackPerEvent(0),
     fHistTPCMult(0),
     fHistTrackletPerEvent(0),
-    fHistMCDaughterTrack(0),
     fHistSPDPrimaryVertexZ(0),
     fHistPrimaryVertexX(0),
     fHistPrimaryVertexY(0),
     fHistPrimaryVertexZ(0),
-    fHistPrimaryVertexResX(0),
-    fHistPrimaryVertexResY(0),
-    fHistPrimaryVertexResZ(0),
-    fHistPrimaryVertexPosXV0events(0), 
-    fHistPrimaryVertexPosYV0events(0), 
-    fHistPrimaryVertexPosZV0events(0),
-    fHistDaughterPt(0),
-    ////////////////////////////////////////
-
-    fHistDcaPosToPrimVertexK0(0),
-    fHistDcaNegToPrimVertexK0(0),
-    fHistRadiusV0K0(0),
-    fHistDecayLengthV0K0(0),
-    fHistDcaV0DaughtersK0(0),
-    fHistChi2K0(0),
-    fHistCosPointAngleK0(0),
-
-
-    fHistDcaPosToPrimVertexK0vsMassK0(0),
-    fHistDcaNegToPrimVertexK0vsMassK0(0),
-    fHistRadiusV0K0vsMassK0(0),
-    fHistDecayLengthV0K0vsMassK0(0),
-    fHistDcaV0DaughtersK0vsMassK0(0),
-    fHistCosPointAngleK0vsMassK0(0),
-    fHistDcaPosToPrimVertexK0vsMassK0pt1(0),
-    fHistDcaNegToPrimVertexK0vsMassK0pt1(0),
-    fHistRadiusV0K0vsMassK0pt1(0),
-    fHistDecayLengthV0K0vsMassK0pt1(0),
-    fHistDcaV0DaughtersK0vsMassK0pt1(0),
-    fHistCosPointAngleK0vsMassK0pt1(0),
-
-    fHistDcaPosToPrimVertexK0vsMassK0pt2(0),
-    fHistDcaNegToPrimVertexK0vsMassK0pt2(0),
-    fHistRadiusV0K0vsMassK0pt2(0),
-    fHistDecayLengthV0K0vsMassK0pt2(0),
-    fHistDcaV0DaughtersK0vsMassK0pt2(0),
-    fHistCosPointAngleK0vsMassK0pt2(0),
-
-    fHistDcaPosToPrimVertexK0vsMassK0pt3(0),
-    fHistDcaNegToPrimVertexK0vsMassK0pt3(0),
-    fHistRadiusV0K0vsMassK0pt3(0),
-    fHistDecayLengthV0K0vsMassK0pt3(0),
-    fHistDcaV0DaughtersK0vsMassK0pt3(0),
-    fHistCosPointAngleK0vsMassK0pt3(0),
-
-    /////////// Lambda ///////////////////////////
-
-    fHistDcaPosToPrimVertexL(0),
-    fHistDcaNegToPrimVertexL(0),
-    fHistRadiusV0L(0),
-    fHistDecayLengthV0L(0),
-    fHistDcaV0DaughtersL(0),
-    fHistChi2L(0),
-    fHistCosPointAngleL(0),
-
-    
-
-    fHistDcaPosToPrimVertexLvsMassL(0),
-    fHistDcaNegToPrimVertexLvsMassL(0),
-    fHistRadiusV0LvsMassL(0),
-    fHistDecayLengthV0LvsMassL(0),
-    fHistDcaV0DaughtersLvsMassL(0),
-    fHistCosPointAngleLvsMassL(0),
-    fHistCosPointAngleLvsMassVsPtsigL(0),
-    fHistCosPointAngleLvsMassVsPtbackL(0),
 
-    
-    fHistDcaPosToPrimVertexLambdaVsMasspt1(0),
-    fHistDcaNegToPrimVertexLambdaVsMasspt1(0),
-    fHistRadiusV0LambdaVsMasspt1(0),
-    fHistDecayLengthV0LambdaVsMasspt1(0),
-    fHistDcaV0DaughtersLambdaVsMasspt1(0),
-    fHistCosPointAngleLambdaVsMasspt1(0),
-
-    fHistDcaPosToPrimVertexLambdaVsMasspt2(0),
-    fHistDcaNegToPrimVertexLambdaVsMasspt2(0),
-    fHistRadiusV0LambdaVsMasspt2(0),
-    fHistDecayLengthV0LambdaVsMasspt2(0),
-    fHistDcaV0DaughtersLambdaVsMasspt2(0),
-    fHistCosPointAngleLambdaVsMasspt2(0),
-
-    fHistDcaPosToPrimVertexLambdaVsMasspt3(0),
-    fHistDcaNegToPrimVertexLambdaVsMasspt3(0),
-    fHistRadiusV0LambdaVsMasspt3(0),
-    fHistDecayLengthV0LambdaVsMasspt3(0),
-    fHistDcaV0DaughtersLambdaVsMasspt3(0),
-    fHistCosPointAngleLambdaVsMasspt3(0),
-
-
-
-    /////////Antilambda ///////////////////
-    fHistDcaPosToPrimVertexAntiL(0),
-    fHistDcaNegToPrimVertexAntiL(0),
-    fHistRadiusV0AntiL(0),
-    fHistDecayLengthV0AntiL(0),
-    fHistDcaV0DaughtersAntiL(0),
-    fHistChi2AntiL(0),
-    fHistCosPointAngleAntiL(0),
-
-
-    fHistDcaPosToPrimVertexAntiLvsMass(0),
-    fHistDcaNegToPrimVertexAntiLvsMass(0),
-    fHistRadiusV0AntiLvsMass(0),
-    fHistDecayLengthV0AntiLvsMass(0),
-    fHistDcaV0DaughtersAntiLvsMass(0),
-    fHistCosPointAngleAntiLvsMass(0),
-
-    
-    
-    fHistDcaPosToPrimVertexAntiLVsMasspt1(0),
-    fHistDcaNegToPrimVertexAntiLVsMasspt1(0),
-    fHistRadiusV0AntiLVsMasspt1(0),
-    fHistDecayLengthV0AntiLVsMasspt1(0),
-    fHistDcaV0DaughtersAntiLVsMasspt1(0),
-    fHistCosPointAngleAntiLVsMasspt1(0),
-
-    fHistDcaPosToPrimVertexAntiLVsMasspt2(0),
-    fHistDcaNegToPrimVertexAntiLVsMasspt2(0),
-    fHistRadiusV0AntiLVsMasspt2(0),
-    fHistDecayLengthV0AntiLVsMasspt2(0),
-    fHistDcaV0DaughtersAntiLVsMasspt2(0),
-    fHistCosPointAngleAntiLVsMasspt2(0),
-
-    fHistDcaPosToPrimVertexAntiLVsMasspt3(0),
-    fHistDcaNegToPrimVertexAntiLVsMasspt3(0),
-    fHistRadiusV0AntiLVsMasspt3(0),
-    fHistDecayLengthV0AntiLVsMasspt3(0),
-    fHistDcaV0DaughtersAntiLVsMasspt3(0),
-    fHistCosPointAngleAntiLVsMasspt3(0),
-
-    /////////////////////////////////////////
     fHistV0Multiplicity(0),
+    ///// inv. mass 
     fHistMassK0(0),
     fHistMassLambda(0),
     fHistMassAntiLambda(0),
-    fHistMassVsRadiusK0(0),
-    fHistMassVsRadiusLambda(0),
-    fHistMassVsRadiusAntiLambda(0),
+    fHistMassXi(0),
+    fHistMassAntiXi(0),
+    fHistMassOmega(0),
+    fHistMassAntiOmega(0),
 
+    //inv mass vs PID
+    fHistMassXiVsPID(0),
     ///////////////////////////////////////
     fHistPtVsMassK0(0),
     fHistPtVsMassLambda(0),
     fHistPtVsMassAntiLambda(0),
 
-    /// Rap3
-    fHistPtVsMassK0Rap3(0),
-    fHistPtVsMassLambdaRap3(0),
-    fHistPtVsMassAntiLambdaRap3(0),
-
-    // cTau dist
-    fHistTranscTauVsMassL(0),
-    fHistTranscTauVsMassAntiL(0),
-    fHistTranscTauVsMassK0s(0),
-
-    // cTau Trans Vs. Mass Rap3
-    fHistTranscTauVsMassLRap3(0),
-    fHistTranscTauVsMassAntiLRap3(0),
-    fHistTranscTauVsMassK0sRap3(0),
-
-    // cTau Trans Vs Mass Low pt
-    fHistTranscTauVsMassLptLow(0),
-    fHistTranscTauVsMassAntiLptLow(0),
-    fHistTranscTauVsMassK0sptLow(0),
-
-    // cTau Trans Vs Mass Low pt Rap3
-    fHistTranscTauVsMassLptLowRap3(0),
-    fHistTranscTauVsMassAntiLptLowRap3(0),
-    fHistTranscTauVsMassK0sptLowRap3(0),
-
     ////////////////////////////////////////
+
     fHistArmenterosPodolanski(0),
     fHistK0sMassVsLambdaMass(0),
     fHistTPCsigPLambda(0),
     fHistTPCsigPAntiLambda(0),
-    fHistNSigmaProton(0),    
-    fHistNsigmaPosPionAntiLambda(0),
-    fHistNsigmaNegProtonAntiLambda(0),
-    fHistNsigmaPosProtonLambda(0),
-    fHistNsigmaNegPionLambda(0),
-    fHistNsigmaPosProtonAntiLambda(0),
-    fHistNsigmaNegPionAntiLambda(0),
-    fHistNsigmaPosPionK0(0),
-    fHistNsigmaNegPionK0(0),
+    fHistNSigmaProton(0),  
+  
+    /// Associated histos 
+    ///rapidity
     fHistAsMcRapK0(0),
-    fHistAsMcRapLambda(0),
-    fHistAsMcRapAntiLambda(0),
+    fHistAsMcRapLambda(0),   fHistAsMcRapAntiLambda(0),
 
-    ///  Pt distribution  /////////////////////
+    // pt distribution    /////////////////////
     fHistAsMcPtK0(0),
-    fHistAsMcPtLambda(0),
-    fHistAsMcPtAntiLambda(0),
-
-    //    cTau distribution
-    fHistAsMcTranscTauL(0),
-    fHistAsMcTranscTauAntiL(0),
-    fHistAsMcTranscTauK0s(0),
-
-    //   cTau Rap3
-    fHistAsMcTranscTauLRap3(0),
-    fHistAsMcTranscTauAntiLRap3(0),
-    fHistAsMcTranscTauK0sRap3(0),
-
-
-    fHistAsMcTranscTauLptLow(0),
-    fHistAsMcTranscTauAntiLptLow(0),
-    fHistAsMcTranscTauK0sptLow(0),
-
-    //Rap3 
-    fHistAsMcTranscTauLptLowRap3(0),
-    fHistAsMcTranscTauAntiLptLowRap3(0),
-    fHistAsMcTranscTauK0sptLowRap3(0),
-    
-
-    ///Rap3
-    fHistAsMcPtK0Rap3(0),
-    fHistAsMcPtLambdaRap3(0),
-    fHistAsMcPtAntiLambdaRap3(0),
-
+    fHistAsMcPtLambda(0),   fHistAsMcPtAntiLambda(0),
 
     /////////////////////////////////////
-
-    fHistAsMcPtZoomK0(0),
-    fHistAsMcPtZoomLambda(0),
-    fHistAsMcPtZoomAntiLambda(0),
-    fHistAsMcProdRadiusK0(0),
-    fHistAsMcProdRadiusLambda(0),
-    fHistAsMcProdRadiusAntiLambda(0),
-    fHistAsMcProdRadiusXvsYK0s(0),
-    fHistAsMcProdRadiusXvsYLambda(0),
-    fHistAsMcProdRadiusXvsYAntiLambda(0),
     fHistPidMcMassK0(0),
     fHistPidMcMassLambda(0),
     fHistPidMcMassAntiLambda(0),
 
-
-    ///////////////////////////////
-
+    ///inv. mass
     fHistAsMcMassK0(0),
-    fHistAsMcMassLambda(0),
-    fHistAsMcMassAntiLambda(0),
-
-    /// Rap3
-    fHistAsMcMassK0Rap3(0),
-    fHistAsMcMassLambdaRap3(0),
-    fHistAsMcMassAntiLambdaRap3(0),
-
-
+    fHistAsMcMassLambda(0), fHistAsMcMassAntiLambda(0),
+    ///transv. momentum
     fHistAsMcPtVsMassK0(0),
-    fHistAsMcPtVsMassLambda(0),
-    fHistAsMcPtVsMassAntiLambda(0),
-
-    /// Rap3
-    fHistAsMcPtVsMassK0Rap3(0),
-    fHistAsMcPtVsMassLambdaRap3(0),
-    fHistAsMcPtVsMassAntiLambdaRap3(0),
+    fHistAsMcPtVsMassLambda(0), fHistAsMcPtVsMassAntiLambda(0),
 
-    ///////////////////////////////////////
+    /////background composition
+    fHistCompositionXi(0),
+    fHistCompositionAntiXi(0),
+    fHistCompositionOmega(0),
+    fHistCompositionAntiOmega(0),
 
-    fHistAsMcMassVsRadiusK0(0),
-    fHistAsMcMassVsRadiusLambda(0),
-    fHistAsMcMassVsRadiusAntiLambda(0),
-    fHistAsMcResxK0(0),
-    fHistAsMcResyK0(0),
-    fHistAsMcReszK0(0),
-    fHistAsMcResrVsRadiusK0(0),
-    fHistAsMcReszVsRadiusK0(0),
-    fHistAsMcResxLambda(0),
-    fHistAsMcResyLambda(0),
-    fHistAsMcReszLambda(0),
-    fHistAsMcResrVsRadiusLambda(0),
-    fHistAsMcReszVsRadiusLambda(0),
-    fHistAsMcResxAntiLambda(0),
-    fHistAsMcResyAntiLambda(0),
-    fHistAsMcReszAntiLambda(0),
-    fHistAsMcResrVsRadiusAntiLambda(0),
-    fHistAsMcReszVsRadiusAntiLambda(0),
-    fHistAsMcResPtK0(0),
-    fHistAsMcResPtLambda(0),
-    fHistAsMcResPtAntiLambda(0),
-    fHistAsMcResPtVsRapK0(0),
-    fHistAsMcResPtVsRapLambda(0),
-    fHistAsMcResPtVsRapAntiLambda(0),
-    fHistAsMcResPtVsPtK0(0),
-    fHistAsMcResPtVsPtLambda(0),
-    fHistAsMcResPtVsPtAntiLambda(0),
-    fHistAsMcMotherPdgCodeK0s(0),
-    fHistAsMcMotherPdgCodeLambda(0),
-    fHistAsMcMotherPdgCodeAntiLambda(0),
-    fHistAsMcPtLambdaFromSigma(0),
-    fHistAsMcPtAntiLambdaFromSigma(0),
-    fHistAsMcSecondaryPtVsRapK0s(0),
-    fHistAsMcSecondaryPtVsRapLambda(0),
-    fHistAsMcSecondaryPtVsRapAntiLambda(0),
-    fHistAsMcSecondaryProdRadiusK0s(0),
-    fHistAsMcSecondaryProdRadiusLambda(0),
-    fHistAsMcSecondaryProdRadiusAntiLambda(0),
-    fHistAsMcSecondaryProdRadiusXvsYK0s(0),
-    fHistAsMcSecondaryProdRadiusXvsYLambda(0),
-    fHistAsMcSecondaryProdRadiusXvsYAntiLambda(0),
-    fHistAsMcSecondaryMotherPdgCodeK0s(0),
-    fHistAsMcSecondaryMotherPdgCodeLambda(0),
-    fHistAsMcSecondaryMotherPdgCodeAntiLambda(0),
-    fHistAsMcSecondaryPtLambdaFromSigma(0),
-    fHistAsMcSecondaryPtAntiLambdaFromSigma(0)
-    
+    fHistMCIndexes(0)
 {
   // Constructor
 
@@ -903,19 +310,10 @@ void AliAnalysisTaskPerformanceStrange::UserCreateOutputObjects()
   fListHist->Add(fHistMCPrimaryVertexZ);
 
   // Pt track distribution
-  fHistPtTracksITSRefit            = new TH1F("h1PtTracksITSRefit", "Pt tracks;p_{t} (GeV/c);Counts", 240,0,12);
-  fListHist->Add(fHistPtTracksITSRefit);
 
   fHistPtTracks            = new TH1F("h1PtTracks", "Pt tracks;p_{t} (GeV/c);Counts", 240,0,12);
   fListHist->Add(fHistPtTracks);
 
-  // Pt track distribution
-  fHistPtTracksPITSRefit            = new TH1F("h1PtTracksPITSRefit", "Pt tracks;p_{t} (GeV/c);Counts", 240,0,12);
-  fListHist->Add(fHistPtTracksPITSRefit);
-
-  fHistPtTracksP            = new TH1F("h1PtTracksP", "Pt tracks;p_{t} (GeV/c);Counts", 240,0,12);
-  fListHist->Add(fHistPtTracksP);
-  
   // Multiplicity
   fHistMCMultiplicityPrimary           = new TH1F("h1MCMultiplicityPrimary", "MC Primary Particles;NPrimary;Count", 201, -0.5, 200.5);
   fListHist->Add(fHistMCMultiplicityPrimary);
@@ -923,26 +321,6 @@ void AliAnalysisTaskPerformanceStrange::UserCreateOutputObjects()
   fHistMCMultiplicityTracks            = new TH1F("h1MCMultiplicityTracks", "MC Tracks;Ntracks;Count", 201, -0.5, 200.5);
   fListHist->Add(fHistMCMultiplicityTracks);
 
-  // Production Radius of non-primary particles:
-  fHistMCtracksProdRadiusK0s           = new TH2F("h2MCtracksProdRadiusK0s","Non-primary MC K^{0} Production Radius;x (cm); y (cm)",200,-50,50,200,-50,50);
-  fListHist->Add(fHistMCtracksProdRadiusK0s);
-
-  fHistMCtracksProdRadiusLambda        = new TH2F("h2MCtracksProdRadiusLambda","Non-primary MC #Lambda^{0} Production Radius;x (cm); y (cm)",200,-50,50,200,-50,50);
-  fListHist->Add(fHistMCtracksProdRadiusLambda);
-
-  fHistMCtracksProdRadiusAntiLambda    = new TH2F("h2MCtracksProdRadiusAntiLambda","Non-primary MC #bar{#Lambda}^{0} Production Radius;x (cm); y (cm)",200,-50,50,200,-50,50);
-  fListHist->Add(fHistMCtracksProdRadiusAntiLambda);
-
-  // Decay Radius of non-primary particles:
-  fHistMCtracksDecayRadiusK0s          = new TH1F("h1MCtracksDecayRadiusK0s","Non-primary MC K^{0} Decay Radius;r (cm)",100,0,100);
-  fListHist->Add(fHistMCtracksDecayRadiusK0s);
-
-  fHistMCtracksDecayRadiusLambda       = new TH1F("h1MCtracksDecayRadiusLambda","Non-primary MC #Lambda^{0} Decay Radius;r (cm)",100,0,100);
-  fListHist->Add(fHistMCtracksDecayRadiusLambda);
-
-  fHistMCtracksDecayRadiusAntiLambda   = new TH1F("h1MCtracksDecayRadiusAntiLambda","Non-primary #bar{#Lambda}^{0} Decay Radius;r (cm)",100,0,100);
-  fListHist->Add(fHistMCtracksDecayRadiusAntiLambda);
-
   // Pt Distribution of non-primary particles:
   fHistMCPtAllK0s                      = new TH1F("h1MCPtAllK0s", "Non-primary MC K^{0};p_{t} (GeV/c);Counts",240,0,12);
   fListHist->Add(fHistMCPtAllK0s);
@@ -953,68 +331,30 @@ void AliAnalysisTaskPerformanceStrange::UserCreateOutputObjects()
   fHistMCPtAllAntiLambda               = new TH1F("h1MCPtAllAntiLambda", "Non-primary MC #bar{#Lambda}^{0};p_{t} (GeV/c);Counts",240,0,12);
   fListHist->Add(fHistMCPtAllAntiLambda);
 
+  fHistMCPtAllXi               = new TH1F("h1MCPtAllXi", "Non-primary MC #Xi;p_{t} (GeV/c);Counts",240,0,12);
+  fListHist->Add(fHistMCPtAllXi);
 
-  // Pt Distribution of non-primary particles Rap3:
-  fHistMCPtAllK0sRap3                      = new TH1F("h1MCPtAllK0sRap3", "Non-primary MC K^{0};p_{t} (GeV/c);Counts",240,0,12);
-  fListHist->Add(fHistMCPtAllK0sRap3);
-
-  fHistMCPtAllLambdaRap3                   = new TH1F("h1MCPtAllLambdaRap3", "Non-primary MC #LambdaRap3^{0};p_{t} (GeV/c);Counts",240,0,12);
-  fListHist->Add(fHistMCPtAllLambdaRap3);
-
-  fHistMCPtAllAntiLambdaRap3               = new TH1F("h1MCPtAllAntiLambdaRap3", "Non-primary MC #bar{#LambdaRap3}^{0};p_{t} (GeV/c);Counts",240,0,12);
-  fListHist->Add(fHistMCPtAllAntiLambdaRap3);
+  fHistMCPtAllAntiXi               = new TH1F("h1MCPtAllAntiXi", "Non-primary MC #bar{#Xi};p_{t} (GeV/c);Counts",240,0,12);
+  fListHist->Add(fHistMCPtAllAntiXi);
 
+  fHistMCPtAllOmega               = new TH1F("h1MCPtAllOmega", "Non-primary MC #Omega;p_{t} (GeV/c);Counts",240,0,12);
+  fListHist->Add(fHistMCPtAllOmega);
 
-  // Production Radius
-  fHistMCProdRadiusK0s                 = new TH1F("h1MCProdRadiusK0s", "MC K^{0} Production Radius;r (cm);Count", 400, -2, 2);
-  fListHist->Add(fHistMCProdRadiusK0s);
-
-  fHistMCProdRadiusLambda              = new TH1F("h1MCProdRadiusLambda", "MC #Lambda^{0} Production Radius;r (cm);Count", 400, -2, 2);
-  fListHist->Add(fHistMCProdRadiusLambda);
+  fHistMCPtAllAntiOmega               = new TH1F("h1MCPtAllAntiOmega", "Non-primary MC #bar{#Omega};p_{t} (GeV/c);Counts",240,0,12);
+  fListHist->Add(fHistMCPtAllAntiOmega);
 
-  fHistMCProdRadiusAntiLambda         = new TH1F("h1MCProdRadiusAntiLambda", "MC #bar{#Lambda}^{0} Production Radius;r (cm);Count", 400, -2, 2);
-  fListHist->Add(fHistMCProdRadiusAntiLambda);
-
-  // Decay radius of primary V0s
-  fHistMCPrimDecayRadiusK0s    = new TH1F("fHistMCPrimDecayRadiusK0s", "MC  K^{0} Decay Radius;r (cm);Count", 100, 0, 100);
-  fListHist->Add(fHistMCPrimDecayRadiusK0s);
-
-  fHistMCPrimDecayRadiusLambda    = new TH1F("fHistMCPrimDecayRadiusLambda", "MC #Lambda^{0} Decay Radius;r (cm);Count", 100, 0, 100);
-  fListHist->Add(fHistMCPrimDecayRadiusLambda);
-
-  fHistMCPrimDecayRadiusAntiLambda    = new TH1F("fHistMCPrimDecayRadiusAntiLambda", "MC #bar{#Lambda}^{0} Decay Radius;r (cm);Count", 100, 0, 100);
-  fListHist->Add(fHistMCPrimDecayRadiusAntiLambda);
 
   // Rapidity distribution:
   fHistMCRapK0s                 = new TH1F("h1MCRapK0s", "K^{0};y",160,-4,4);
   fListHist->Add(fHistMCRapK0s);
 
-  fHistMCRapInPtRangeK0s        = new TH1F("h1MCRapInPtRangeK0s", "K^{0};y",160,-4,4);
-  fListHist->Add(fHistMCRapInPtRangeK0s);
-
   fHistMCRapLambda              = new TH1F("h1MCRapLambda", "#Lambda;y",160,-4,4);
   fListHist->Add(fHistMCRapLambda);
 
-  fHistMCRapInPtRangeLambda     = new TH1F("h1MCRapInPtRangeLambda", "#Lambda;y",160,-4,4);
-  fListHist->Add(fHistMCRapInPtRangeLambda);
 
   fHistMCRapAntiLambda          = new TH1F("h1MCRapAntiLambda", "#bar{#Lambda};y",160,-4,4);
   fListHist->Add(fHistMCRapAntiLambda);
 
-  fHistMCRapInPtRangeAntiLambda = new TH1F("h1MCRapInPtRangeAntiLambda", "#bar{#Lambda};y",160,-4,4);
-  fListHist->Add(fHistMCRapInPtRangeAntiLambda);
-
-  fHistMCRapXi                  = new TH1F("h1MCRapXi", "Xi;y",160,-4,4);
-  fListHist->Add(fHistMCRapXi);
-
-  fHistMCRapInPtRangeXi         = new TH1F("h1MCRapInPtRangeXi", "Xi;y",160,-4,4);
-  fListHist->Add(fHistMCRapInPtRangeXi);
-
-  fHistMCRapPhi                  = new TH1F("h1MCRapPhi", "Phi;y",160,-4,4);
-  fListHist->Add(fHistMCRapPhi);
-
-  fHistMCRapInPtRangePhi         = new TH1F("h1MCRapInPtRangePhi", "Phi;y",160,-4,4);
-  fListHist->Add(fHistMCRapInPtRangePhi);
 
   // Pt distribution:
   fHistMCPtK0s               = new TH1F("h1MCPtK0s", "K^{0};p_{t} (GeV/c)",240,0,12);
@@ -1026,44 +366,6 @@ void AliAnalysisTaskPerformanceStrange::UserCreateOutputObjects()
   fHistMCPtAntiLambda            = new TH1F("h1MCPtAntiLambda", "#AntiLambda^{0};p_{t} (GeV/c)",240,0,12);
   fListHist->Add(fHistMCPtAntiLambda);
 
-
-  // Pt distribution Rap3:
-  fHistMCPtK0sRap3               = new TH1F("h1MCPtK0sRap3", "K^{0};p_{t} (GeV/c)",240,0,12);
-  fListHist->Add(fHistMCPtK0sRap3);
-
-  fHistMCPtLambdaRap3            = new TH1F("h1MCPtLambdaRap3", "#Lambda^{0};p_{t} (GeV/c)",240,0,12);
-  fListHist->Add(fHistMCPtLambdaRap3);
-
-  fHistMCPtAntiLambdaRap3            = new TH1F("h1MCPtAntiLambdaRap3", "#AntiLambda^{0};p_{t} (GeV/c)",240,0,12);
-  fListHist->Add(fHistMCPtAntiLambdaRap3);
-
-
-  // Pt distribution of Lambda coming from Sigma decay
-  fHistMCPtLambdaFromSigma      = new TH1F("h1MCPtLambdaFromSigma", "#Lambda^{0};p_{t} (GeV/c)",240,0,12);
-  fListHist->Add(fHistMCPtLambdaFromSigma);
-
-  fHistMCPtAntiLambdaFromSigma  = new TH1F("h1MCPtAntiLambdaFromSigma", "#Lambda^{0};p_{t} (GeV/c)",240,0,12);
-  fListHist->Add(fHistMCPtAntiLambdaFromSigma);
-  // Multiple reconstruction studies
-  fHistNTimesRecK0s             = new TH1F("h1NTimesRecK0s","number of times a K0s is reconstructed in -1<y<1;number of times;counts",500,-0.5,4.5);
-  fListHist->Add(fHistNTimesRecK0s);
-
-  fHistNTimesRecLambda          = new TH1F("h1NTimesRecLambda","number of times a Lambda is reconstructed in -1<y<1;number of times;counts",500,-0.5,4.5);
-  fListHist->Add(fHistNTimesRecLambda);
-
-  fHistNTimesRecAntiLambda      = new TH1F("h1NTimesRecAntiLambda","number of times an AntiLambda is reconstructed in -1<y<1;number of times;counts",500,-0.5,4.5);
-  fListHist->Add(fHistNTimesRecAntiLambda);
-
-  fHistNTimesRecK0sVsPt         = new TH2F("h2NTimesRecK0sVsPt","NTimes versus Pt, K^{0} in -1<y<1;p_{t} (GeV/c);number of times",75,0,15,5,-0.5,4.5);
-  fListHist->Add(fHistNTimesRecK0sVsPt);
-
-  fHistNTimesRecLambdaVsPt      = new TH2F("h2NTimesRecLambdaVsPt","NTimes versus Pt, #Lambda^{0} in -1<y<1;p_{t} (GeV/c);number of times",75,0,15,5,-0.5,4.5);
-  fListHist->Add(fHistNTimesRecLambdaVsPt);
-
-  fHistNTimesRecAntiLambdaVsPt  = new TH2F("h2NTimesRecAntiLambdaVsPt","NTimes versus Pt, #bar{#Lambda}^{0} in -1<y<1;p_{t} (GeV/c);number of times",75,0,15,5,-0.5,4.5);
-  fListHist->Add(fHistNTimesRecAntiLambdaVsPt);
-
   //***********************************
   // Reconstructed particles histograms
   //***********************************
@@ -1083,9 +385,6 @@ void AliAnalysisTaskPerformanceStrange::UserCreateOutputObjects()
   fHistTrackletPerEvent       = new TH1F("h1TrackletPerEvent", "Number of tracklets;Number of tracklets per events;Number of events",1000,0,1000);
   fListHist->Add(fHistTrackletPerEvent);
 
-  fHistMCDaughterTrack         = new TH1F("h1MCDaughterTrack","Distribution of mc id for daughters;id tags;Counts",15,0,15);
-  fListHist->Add(fHistMCDaughterTrack);
-
   fHistTPCTracks               = new TH1F("h1TPCTracks","Distribution of TPC tracks;Number of TPC tracks:Number of events",1000,0,10000);
   fListHist->Add(fHistTPCTracks);
 
@@ -1102,345 +401,8 @@ void AliAnalysisTaskPerformanceStrange::UserCreateOutputObjects()
   fHistPrimaryVertexZ          = new TH1F("h1PrimaryVertexZ", "Primary Vertex Position Z;Primary Vertex Position Z (cm);Events",200,-20,20);
   fListHist->Add(fHistPrimaryVertexZ);
 
-
-  // Primary vertex resolution
-  fHistPrimaryVertexResX          = new TH1F("h1PrimaryVertexResX", "Primary Vertex Resolution X;Primary Vertex Resolution X (cm);Events",100,-0.25,0.25);
-  fListHist->Add(fHistPrimaryVertexResX);
-
-  fHistPrimaryVertexResY          = new TH1F("h1PrimaryVertexResY", "Primary Vertex Resolution Y;Primary Vertex Resolution Y (cm);Events",100,-0.25,0.25);
-  fListHist->Add(fHistPrimaryVertexResY);
-
-  fHistPrimaryVertexResZ          = new TH1F("h1PrimaryVertexResZ", "Primary Vertex Resolution Z;Primary Vertex Resolution Z (cm);Events",200,-0.25,0.25);
-  fListHist->Add(fHistPrimaryVertexResZ);
-  
-
-  // Primary Vertex in events with V0 candidates:
-  fHistPrimaryVertexPosXV0events       = new TH1F("h1PrimaryVertexPosXV0events", "Primary Vertex Position X;Primary Vertex Position X (cm);Events",100,-0.5,0.5);
-  fListHist->Add(fHistPrimaryVertexPosXV0events);
-  fHistPrimaryVertexPosYV0events       = new TH1F("h1PrimaryVertexPosYV0events", "Primary Vertex Position Y;Primary Vertex Position Y (cm);Events",100,-0.5,0.5);
-  fListHist->Add(fHistPrimaryVertexPosYV0events);
-  fHistPrimaryVertexPosZV0events       = new TH1F("h1PrimaryVertexPosZV0events", "Primary Vertex Position Z;Primary Vertex Position Z (cm);Events",200,-20.0,20.0);
-  fListHist->Add(fHistPrimaryVertexPosZV0events);
-
-  // Daughters Pt:
-  fHistDaughterPt              = new TH2F("h2DaughterPt", "Daughter Pt;Positive Daughter Pt; Negative Daughter Pt",200,0,2,200,0,2);
-  fListHist->Add(fHistDaughterPt);
-
   //////K0s///////////////// 2D histos: cut vs on fly status////
 
-  fHistDcaPosToPrimVertexK0      = new TH2F("h2DcaPosToPrimVertexK0", "Positive V0 daughter;dca(cm);Status",100,0,10,2,-0.5,1.5);
-  fListHist->Add(fHistDcaPosToPrimVertexK0);
-
-  fHistDcaNegToPrimVertexK0      = new TH2F("h2DcaNegToPrimVertexK0", "Negative V0 daughter;dca(cm);Status",100,0,10,2,-0.5,1.5);
-  fListHist->Add(fHistDcaNegToPrimVertexK0);
-
-
-  fHistRadiusV0K0                = new TH2F("h2RadiusV0K0", "Radius;Radius(cm);Status",500,0,500,2,-0.5,1.5);
-  fListHist->Add(fHistRadiusV0K0);
-
-  fHistDecayLengthV0K0           = new TH2F("h2DecayLengthV0K0", "V0s decay Length;decay length(cm);Status", 500, 0, 500,2,-0.5,1.5);
-  fListHist->Add(fHistDecayLengthV0K0);
-
-  fHistDcaV0DaughtersK0          = new TH2F("h2DcaV0DaughtersK0", "DCA between daughters;dca(cm);Status", 300, 0, 3.0,2,-0.5,1.5);
-  fListHist->Add(fHistDcaV0DaughtersK0);
-
-  fHistChi2K0                    = new TH2F("h2Chi2K0", "V0s chi2;chi2;Status", 1000, 0, 0.1,2,-0.5,1.5);
-  fListHist->Add(fHistChi2K0);
-
-  fHistCosPointAngleK0           = new TH2F("h2CosPointAngleK0", "Cosine of V0's pointing angle", 200,0.99,1.01,2,-0.5,1.5);
-  fListHist->Add(fHistCosPointAngleK0);
-
-
-  ////////////K0s///////////////// 2D histos: cut vs mass////
-
-
-  fHistDcaPosToPrimVertexK0vsMassK0 = new TH2F("h2DcaPosToPrimVertexK0vsMassK0", "Positive V0 daughter;dca(cm);K0s inv. mass",500,0,10,200,0.4,0.6);
-  fListHist->Add(fHistDcaPosToPrimVertexK0vsMassK0);
-
-  fHistDcaNegToPrimVertexK0vsMassK0 = new TH2F("h2DcaNegToPrimVertexK0vsMassK0", "Negative V0 daughter;dca(cm);K0s inv. mass",500,0,10,200,0.4,0.6);
-  fListHist->Add(fHistDcaNegToPrimVertexK0vsMassK0);
-
-
-  fHistRadiusV0K0vsMassK0           = new TH2F("h2RadiusV0K0vsMassK0", "Radius;Radius(cm);K0s inv. mass",110,0,110,200,0.4,0.6);
-  fListHist->Add(fHistRadiusV0K0vsMassK0);
-
-  fHistDecayLengthV0K0vsMassK0      = new TH2F("h2DecayLengthV0K0vsMassK0", "V0s decay Length;decay length(cm);K0s inv. mass", 100, 0, 100,200,0.4,0.6);
-  fListHist->Add(fHistDecayLengthV0K0vsMassK0);
-
-  fHistDcaV0DaughtersK0vsMassK0     = new TH2F("h2DcaV0DaughtersK0vsMassK0", "DCA between daughters;dca(cm);K0s inv. mass", 110, 0, 1.1,200,0.4,0.6);
-  fListHist->Add(fHistDcaV0DaughtersK0vsMassK0);
-
-
-  fHistCosPointAngleK0vsMassK0      = new TH2F("h2CosPointAngleK0vsMassK0", "Cosine of V0's pointing angle", 200,0.997,1.007,200,0.4,0.6);
-  fListHist->Add(fHistCosPointAngleK0vsMassK0);
-  if(fQASelector) {
-    //// pt1
-    fHistDcaPosToPrimVertexK0vsMassK0pt1 = new TH2F("h2DcaPosToPrimVertexK0vsMassK0pt1", "Positive V0 daughter;dca(cm);K0s inv. mass",500,0,10,200,0.4,0.6);
-    fListHist->Add(fHistDcaPosToPrimVertexK0vsMassK0pt1);
-
-    fHistDcaNegToPrimVertexK0vsMassK0pt1 = new TH2F("h2DcaNegToPrimVertexK0vsMassK0pt1", "Negative V0 daughter;dca(cm);K0s inv. mass",500,0,10,200,0.4,0.6);
-    fListHist->Add(fHistDcaNegToPrimVertexK0vsMassK0pt1);
-
-    fHistRadiusV0K0vsMassK0pt1           = new TH2F("h2RadiusV0K0vsMassK0pt1", "Radius;Radius(cm);K0s inv. mass",110,0,110,200,0.4,0.6);
-    fListHist->Add(fHistRadiusV0K0vsMassK0pt1);
-
-    fHistDecayLengthV0K0vsMassK0pt1      = new TH2F("h2DecayLengthV0K0vsMassK0pt1", "V0s decay Length;decay length(cm);K0s inv. mass",100,0,100,200,0.4,0.6);
-    fListHist->Add(fHistDecayLengthV0K0vsMassK0pt1);
-
-    fHistDcaV0DaughtersK0vsMassK0pt1     = new TH2F("h2DcaV0DaughtersK0vsMassK0pt1", "DCA between daughters;dca(cm);K0s inv. mass", 110, 0, 1.1,200,0.4,0.6);
-    fListHist->Add(fHistDcaV0DaughtersK0vsMassK0pt1);
-
-    fHistCosPointAngleK0vsMassK0pt1      = new TH2F("h2CosPointAngleK0vsMassK0pt1", "Cosine of V0's pointing angle", 200,0.997,1.007,200,0.4,0.6);
-    fListHist->Add(fHistCosPointAngleK0vsMassK0pt1);
-
-    /// pt2
-    fHistDcaPosToPrimVertexK0vsMassK0pt2 = new TH2F("h2DcaPosToPrimVertexK0vsMassK0pt2", "Positive V0 daughter;dca(cm);K0s inv. mass",500,0,10,200,0.4,0.6);
-    fListHist->Add(fHistDcaPosToPrimVertexK0vsMassK0pt2);
-
-    fHistDcaNegToPrimVertexK0vsMassK0pt2 = new TH2F("h2DcaNegToPrimVertexK0vsMassK0pt2", "Negative V0 daughter;dca(cm);K0s inv. mass",500,0,10,200,0.4,0.6);
-    fListHist->Add(fHistDcaNegToPrimVertexK0vsMassK0pt2);
-
-    fHistRadiusV0K0vsMassK0pt2           = new TH2F("h2RadiusV0K0vsMassK0pt2", "Radius;Radius(cm);K0s inv. mass",110,0,110,200,0.4,0.6);
-    fListHist->Add(fHistRadiusV0K0vsMassK0pt2);
-
-    fHistDecayLengthV0K0vsMassK0pt2      = new TH2F("h2DecayLengthV0K0vsMassK0pt2", "V0s decay Length;decay length(cm);K0s inv. mass", 100, 0, 100,200,0.4,0.6);
-    fListHist->Add(fHistDecayLengthV0K0vsMassK0pt2);
-
-    fHistDcaV0DaughtersK0vsMassK0pt2     = new TH2F("h2DcaV0DaughtersK0vsMassK0pt2", "DCA between daughters;dca(cm);K0s inv. mass", 110, 0, 1.1,200,0.4,0.6);
-    fListHist->Add(fHistDcaV0DaughtersK0vsMassK0pt2);
-  
-    fHistCosPointAngleK0vsMassK0pt2      = new TH2F("h2CosPointAngleK0vsMassK0pt2", "Cosine of V0's pointing angle", 200,0.997,1.007,200,0.4,0.6);
-    fListHist->Add(fHistCosPointAngleK0vsMassK0pt2);
-
-    /// pt3
-    fHistDcaPosToPrimVertexK0vsMassK0pt3 = new TH2F("h2DcaPosToPrimVertexK0vsMassK0pt3", "Positive V0 daughter;dca(cm);K0s inv. mass",500,0,10,200,0.4,0.6);
-    fListHist->Add(fHistDcaPosToPrimVertexK0vsMassK0pt3);
-
-    fHistDcaNegToPrimVertexK0vsMassK0pt3 = new TH2F("h2DcaNegToPrimVertexK0vsMassK0pt3", "Negative V0 daughter;dca(cm);K0s inv. mass",500,0,10,200,0.4,0.6);
-    fListHist->Add(fHistDcaNegToPrimVertexK0vsMassK0pt3);
-
-    fHistRadiusV0K0vsMassK0pt3           = new TH2F("h2RadiusV0K0vsMassK0pt3", "Radius;Radius(cm);K0s inv. mass",110,0,110,200,0.4,0.6);
-    fListHist->Add(fHistRadiusV0K0vsMassK0pt3);
-
-    fHistDecayLengthV0K0vsMassK0pt3      = new TH2F("h2DecayLengthV0K0vsMassK0pt3", "V0s decay Length;decay length(cm);K0s inv. mass", 100, 0, 100,200,0.4,0.6);
-    fListHist->Add(fHistDecayLengthV0K0vsMassK0pt3);
-
-    fHistDcaV0DaughtersK0vsMassK0pt3     = new TH2F("h2DcaV0DaughtersK0vsMassK0pt3", "DCA between daughters;dca(cm);K0s inv. mass", 110, 0, 1.1,200,0.4,0.6);
-    fListHist->Add(fHistDcaV0DaughtersK0vsMassK0pt3);
-  
-    fHistCosPointAngleK0vsMassK0pt3      = new TH2F("h2CosPointAngleK0vsMassK0pt3", "Cosine of V0's pointing angle", 200,0.997,1.007,200,0.4,0.6);
-    fListHist->Add(fHistCosPointAngleK0vsMassK0pt3);
-  } //end QA condition
-  //////////Lambda////////////// 2D histos: cut vs on fly status////
-
-  fHistDcaPosToPrimVertexL      = new TH2F("h2DcaPosToPrimVertexL", "Positive V0 daughter;dca(cm);Status",100,0,10,2,-0.5,1.5);
-  fListHist->Add(fHistDcaPosToPrimVertexL);
-
-  fHistDcaNegToPrimVertexL      = new TH2F("h2DcaNegToPrimVertexL", "Negative V0 daughter;dca(cm);Status",100,0,10,2,-0.5,1.5);
-  fListHist->Add(fHistDcaNegToPrimVertexL);
-
-  fHistRadiusV0L                = new TH2F("h2RadiusV0L", "Radius;Radius(cm);Status",100,0,110,2,-0.5,1.5);
-  fListHist->Add(fHistRadiusV0L);
-
-  fHistDecayLengthV0L           = new TH2F("h2DecayLengthV0L", "V0s decay Length;decay length(cm);Status", 500, 0, 500,2,-0.5,1.5);
-  fListHist->Add(fHistDecayLengthV0L);
-
-  fHistDcaV0DaughtersL          = new TH2F("h2DcaV0DaughtersL", "DCA between daughters;dca(cm);Status", 300, 0, 3.0,2,-0.5,1.5);
-  fListHist->Add(fHistDcaV0DaughtersL);
-
-  fHistChi2L                    = new TH2F("h2Chi2L", "V0s chi2;chi2;Status", 100, 0, 0.10,2,-0.5,1.5);
-  fListHist->Add(fHistChi2L);
-
-  fHistCosPointAngleL           = new TH2F("h2CosPointAngleL", "Cosine of V0's pointing angle", 200,0.99,1.01,2,-0.5,1.5);
-  fListHist->Add(fHistCosPointAngleL);
-
-  //////////Lambda////////////// 2D histos: cut vs mass////
-  fHistDcaPosToPrimVertexLvsMassL      = new TH2F("h2DcaPosToPrimVertexLvsMassL", "Positive V0 daughter;dca(cm);Status",100,0,10,140, 1.06, 1.2);
-  fListHist->Add(fHistDcaPosToPrimVertexLvsMassL);
-
-  fHistDcaNegToPrimVertexLvsMassL      = new TH2F("h2DcaNegToPrimVertexLvsMassL", "Negative V0 daughter;dca(cm);Status",100,0,10,140, 1.06, 1.2);
-  fListHist->Add(fHistDcaNegToPrimVertexLvsMassL);
-
-
-  fHistRadiusV0LvsMassL                = new TH2F("h2RadiusV0LvsMassL", "Radius;Radius(cm);Status",110,0,110,140, 1.06, 1.2);
-  fListHist->Add(fHistRadiusV0LvsMassL);
-
-  fHistDecayLengthV0LvsMassL           = new TH2F("h2DecayLengthV0LvsMassL", "V0s decay Length;decay length(cm);Status", 120, 0, 120,140, 1.06, 1.2);
-  fListHist->Add(fHistDecayLengthV0LvsMassL);
-
-  fHistDcaV0DaughtersLvsMassL          = new TH2F("h2DcaV0DaughtersLvsMassL", "DCA between daughters;dca(cm);Status", 110, 0, 1.1,140, 1.06, 1.2);
-  fListHist->Add(fHistDcaV0DaughtersLvsMassL);
-
-  fHistCosPointAngleLvsMassL           = new TH2F("h2CosPointAngleLvsMassL", "Cosine of V0's pointing angle", 200,0.997,1.007,140, 1.06, 1.2);
-  fListHist->Add(fHistCosPointAngleLvsMassL);
-  /*
-    fHistCosPointAngleLvsMassVsPtsigL           = new TH3F("h3McCosPointAngleLvsMassVsPtsigL", "Cosine of V0's pointing angle",3,0,12, 2,00.997,1.007,140, 1.06, 1.2);
-    fListHist->Add(fHistCosPointAngleLvsMassVsPtsigL);
-    fHistCosPointAngleLvsMassVsPtbackL           = new TH3F("h3McCosPointAngleLvsMassVsPtbackL", "Cosine of V0's pointing angle",3,0,12, 20,0.997,1.007,140, 1.06, 1.2);
-    fListHist->Add(fHistCosPointAngleLvsMassVsPtbackL);
-  */
-
-  if(fQASelector){
-    //// pt1
-    fHistDcaPosToPrimVertexLambdaVsMasspt1 = new TH2F("h2DcaPosToPrimVertexLambdaVsMasspt1", "Positive V0 daughter;dca(cm);K0s inv. mass",500,0,10,140,1.06,1.2);
-    fListHist->Add(fHistDcaPosToPrimVertexLambdaVsMasspt1);
-
-    fHistDcaNegToPrimVertexLambdaVsMasspt1 = new TH2F("h2DcaNegToPrimVertexLambdaVsMasspt1", "Negative V0 daughter;dca(cm);K0s inv. mass",500,0,10,140,1.06,1.2);
-    fListHist->Add(fHistDcaNegToPrimVertexLambdaVsMasspt1);
-
-    fHistRadiusV0LambdaVsMasspt1           = new TH2F("h2RadiusV0LambdaVsMasspt1", "Radius;Radius(cm);K0s inv. mass",110,0,110,140,1.06,1.2);
-    fListHist->Add(fHistRadiusV0LambdaVsMasspt1);
-
-    fHistDecayLengthV0LambdaVsMasspt1      = new TH2F("h2DecayLengthV0LambdaVsMasspt1", "V0s decay Length;decay length(cm);K0s inv. mass",100,0,100,140,1.06,1.2);
-    fListHist->Add(fHistDecayLengthV0LambdaVsMasspt1);
-
-    fHistDcaV0DaughtersLambdaVsMasspt1     = new TH2F("h2DcaV0DaughtersLambdaVsMasspt1", "DCA between daughters;dca(cm);K0s inv. mass", 110, 0, 1.1,140,1.06,1.2);
-    fListHist->Add(fHistDcaV0DaughtersLambdaVsMasspt1);
-
-    fHistCosPointAngleLambdaVsMasspt1      = new TH2F("h2CosPointAngleLambdaVsMasspt1", "Cosine of V0's pointing angle", 200,0.997,1.007,140,1.06,1.2);
-    fListHist->Add(fHistCosPointAngleLambdaVsMasspt1);
-
-    /// pt2
-    fHistDcaPosToPrimVertexLambdaVsMasspt2 = new TH2F("h2DcaPosToPrimVertexLambdaVsMasspt2", "Positive V0 daughter;dca(cm);K0s inv. mass",500,0,10,140,1.06,1.2);
-    fListHist->Add(fHistDcaPosToPrimVertexLambdaVsMasspt2);
-
-    fHistDcaNegToPrimVertexLambdaVsMasspt2 = new TH2F("h2DcaNegToPrimVertexLambdaVsMasspt2", "Negative V0 daughter;dca(cm);K0s inv. mass",500,0,10,140,1.06,1.2);
-    fListHist->Add(fHistDcaNegToPrimVertexLambdaVsMasspt2);
-
-    fHistRadiusV0LambdaVsMasspt2           = new TH2F("h2RadiusV0LambdaVsMasspt2", "Radius;Radius(cm);K0s inv. mass",110,0,110,140,1.06,1.2);
-    fListHist->Add(fHistRadiusV0LambdaVsMasspt2);
-
-    fHistDecayLengthV0LambdaVsMasspt2      = new TH2F("h2DecayLengthV0LambdaVsMasspt2", "V0s decay Length;decay length(cm);K0s inv. mass", 100, 0, 100,140,1.06,1.2);
-    fListHist->Add(fHistDecayLengthV0LambdaVsMasspt2);
-
-    fHistDcaV0DaughtersLambdaVsMasspt2     = new TH2F("h2DcaV0DaughtersLambdaVsMasspt2", "DCA between daughters;dca(cm);K0s inv. mass", 110, 0, 1.1,140,1.06,1.2);
-    fListHist->Add(fHistDcaV0DaughtersLambdaVsMasspt2);
-  
-    fHistCosPointAngleLambdaVsMasspt2      = new TH2F("h2CosPointAngleLambdaVsMasspt2", "Cosine of V0's pointing angle", 200,0.997,1.007,140,1.06,1.2);
-    fListHist->Add(fHistCosPointAngleLambdaVsMasspt2);
-
-    /// pt3
-    fHistDcaPosToPrimVertexLambdaVsMasspt3 = new TH2F("h2DcaPosToPrimVertexLambdaVsMasspt3", "Positive V0 daughter;dca(cm);K0s inv. mass",500,0,10,140,1.06,1.2);
-    fListHist->Add(fHistDcaPosToPrimVertexLambdaVsMasspt3);
-
-    fHistDcaNegToPrimVertexLambdaVsMasspt3 = new TH2F("h2DcaNegToPrimVertexLambdaVsMasspt3", "Negative V0 daughter;dca(cm);K0s inv. mass",500,0,10,140,1.06,1.2);
-    fListHist->Add(fHistDcaNegToPrimVertexLambdaVsMasspt3);
-
-    fHistRadiusV0LambdaVsMasspt3           = new TH2F("h2RadiusV0LambdaVsMasspt3", "Radius;Radius(cm);K0s inv. mass",110,0,110,140,1.06,1.2);
-    fListHist->Add(fHistRadiusV0LambdaVsMasspt3);
-
-    fHistDecayLengthV0LambdaVsMasspt3      = new TH2F("h2DecayLengthV0LambdaVsMasspt3", "V0s decay Length;decay length(cm);K0s inv. mass", 100, 0, 100,140,1.06,1.2);
-    fListHist->Add(fHistDecayLengthV0LambdaVsMasspt3);
-
-    fHistDcaV0DaughtersLambdaVsMasspt3     = new TH2F("h2DcaV0DaughtersLambdaVsMasspt3", "DCA between daughters;dca(cm);K0s inv. mass", 110, 0, 1.1,140,1.06,1.2);
-    fListHist->Add(fHistDcaV0DaughtersLambdaVsMasspt3);
-  
-    fHistCosPointAngleLambdaVsMasspt3      = new TH2F("h2CosPointAngleLambdaVsMasspt3", "Cosine of V0's pointing angle", 200,0.997,1.007,140,1.06,1.2);
-    fListHist->Add(fHistCosPointAngleLambdaVsMasspt3);
-  } // end QA condition
-  //////////AntiLambda////////////// 2D histos: cut vs on fly status////
-
-  fHistDcaPosToPrimVertexAntiL      = new TH2F("h2DcaPosToPrimVertexAntiL", "Positive V0 daughter;dca(cm);Status",100,0,10,2,-0.5,1.5);
-  fListHist->Add(fHistDcaPosToPrimVertexAntiL);
-
-  fHistDcaNegToPrimVertexAntiL      = new TH2F("h2DcaNegToPrimVertexAntiL", "Negative V0 daughter;dca(cm);Status",100,0,10,2,-0.5,1.5);
-  fListHist->Add(fHistDcaNegToPrimVertexAntiL);
-
-
-  fHistRadiusV0AntiL                = new TH2F("h2RadiusV0AntiL", "Radius;Radius(cm);Status",100,0,110,2,-0.5,1.5);
-  fListHist->Add(fHistRadiusV0AntiL);
-
-  fHistDecayLengthV0AntiL           = new TH2F("h2DecayLengthV0AntiL", "V0s decay Length;decay length(cm);Status", 500, 0, 500,2,-0.5,1.5);
-  fListHist->Add(fHistDecayLengthV0AntiL);
-
-  fHistDcaV0DaughtersAntiL          = new TH2F("h2DcaV0DaughtersAntiL", "DCA between daughters;dca(cm);Status", 300, 0, 3.0,2,-0.5,1.5);
-  fListHist->Add(fHistDcaV0DaughtersAntiL);
-
-  fHistChi2AntiL                    = new TH2F("h2Chi2AntiL", "V0s chi2;chi2;Status", 100, 0, 0.10,2,-0.5,1.5);
-  fListHist->Add(fHistChi2AntiL);
-
-  fHistCosPointAngleAntiL           = new TH2F("h2CosPointAngleAntiL", "Cosine of V0's pointing angle", 200,0.99,1.01,2,-0.5,1.5);
-  fListHist->Add(fHistCosPointAngleAntiL);
-
-  //////////AntiLambda////////////// 2D histos: cut vs mass////
-
-  fHistDcaPosToPrimVertexAntiLvsMass      = new TH2F("h2DcaPosToPrimVertexAntiLvsMass", "Positive V0 daughter;dca(cm);Status",100,0,10,140, 1.06, 1.2);
-  fListHist->Add(fHistDcaPosToPrimVertexAntiLvsMass);
-
-  fHistDcaNegToPrimVertexAntiLvsMass      = new TH2F("h2DcaNegToPrimVertexAntiLvsMass", "Negative V0 daughter;dca(cm);Status",100,0,10,140, 1.06, 1.2);
-  fListHist->Add(fHistDcaNegToPrimVertexAntiLvsMass);
-
-
-  fHistRadiusV0AntiLvsMass                = new TH2F("h2RadiusV0AntiLvsMass", "Radius;Radius(cm);Status",110,0,110,140, 1.06, 1.2);
-  fListHist->Add(fHistRadiusV0AntiLvsMass);
-
-  fHistDecayLengthV0AntiLvsMass           = new TH2F("h2DecayLengthV0AntiLvsMass", "V0s decay Length;decay length(cm);Status", 120, 0, 120,140, 1.06, 1.2);
-  fListHist->Add(fHistDecayLengthV0AntiLvsMass);
-
-  fHistDcaV0DaughtersAntiLvsMass          = new TH2F("h2DcaV0DaughtersAntiLvsMass", "DCA between daughters;dca(cm);Status", 110, 0, 1.1,140, 1.06, 1.2);
-  fListHist->Add(fHistDcaV0DaughtersAntiLvsMass);
-
-  fHistCosPointAngleAntiLvsMass           = new TH2F("h2CosPointAngleAntiLvsMass", "Cosine of V0's pointing angle", 200,0.997,1.007,140, 1.06, 1.2);
-  fListHist->Add(fHistCosPointAngleAntiLvsMass);
-  if(fQASelector){
-    //// pt1
-    fHistDcaPosToPrimVertexAntiLVsMasspt1 = new TH2F("h2DcaPosToPrimVertexAntiLVsMasspt1", "Positive V0 daughter;dca(cm);K0s inv. mass",500,0,10,140,1.06,1.2);
-    fListHist->Add(fHistDcaPosToPrimVertexAntiLVsMasspt1);
-
-    fHistDcaNegToPrimVertexAntiLVsMasspt1 = new TH2F("h2DcaNegToPrimVertexAntiLVsMasspt1", "Negative V0 daughter;dca(cm);K0s inv. mass",500,0,10,140,1.06,1.2);
-    fListHist->Add(fHistDcaNegToPrimVertexAntiLVsMasspt1);
-
-    fHistRadiusV0AntiLVsMasspt1           = new TH2F("h2RadiusV0AntiLVsMasspt1", "Radius;Radius(cm);K0s inv. mass",110,0,110,140,1.06,1.2);
-    fListHist->Add(fHistRadiusV0AntiLVsMasspt1);
-
-    fHistDecayLengthV0AntiLVsMasspt1      = new TH2F("h2DecayLengthV0AntiLVsMasspt1", "V0s decay Length;decay length(cm);K0s inv. mass",100,0,100,140,1.06,1.2);
-    fListHist->Add(fHistDecayLengthV0AntiLVsMasspt1);
-
-    fHistDcaV0DaughtersAntiLVsMasspt1     = new TH2F("h2DcaV0DaughtersAntiLVsMasspt1", "DCA between daughters;dca(cm);K0s inv. mass", 110, 0, 1.1,140,1.06,1.2);
-    fListHist->Add(fHistDcaV0DaughtersAntiLVsMasspt1);
-
-    fHistCosPointAngleAntiLVsMasspt1      = new TH2F("h2CosPointAngleAntiLVsMasspt1", "Cosine of V0's pointing angle", 200,0.997,1.007,140,1.06,1.2);
-    fListHist->Add(fHistCosPointAngleAntiLVsMasspt1);
-
-    /// pt2
-    fHistDcaPosToPrimVertexAntiLVsMasspt2 = new TH2F("h2DcaPosToPrimVertexAntiLVsMasspt2", "Positive V0 daughter;dca(cm);K0s inv. mass",500,0,10,140,1.06,1.2);
-    fListHist->Add(fHistDcaPosToPrimVertexAntiLVsMasspt2);
-
-    fHistDcaNegToPrimVertexAntiLVsMasspt2 = new TH2F("h2DcaNegToPrimVertexAntiLVsMasspt2", "Negative V0 daughter;dca(cm);K0s inv. mass",500,0,10,140,1.06,1.2);
-    fListHist->Add(fHistDcaNegToPrimVertexAntiLVsMasspt2);
-
-    fHistRadiusV0AntiLVsMasspt2           = new TH2F("h2RadiusV0AntiLVsMasspt2", "Radius;Radius(cm);K0s inv. mass",110,0,110,140,1.06,1.2);
-    fListHist->Add(fHistRadiusV0AntiLVsMasspt2);
-
-    fHistDecayLengthV0AntiLVsMasspt2      = new TH2F("h2DecayLengthV0AntiLVsMasspt2", "V0s decay Length;decay length(cm);K0s inv. mass", 100, 0, 100,140,1.06,1.2);
-    fListHist->Add(fHistDecayLengthV0AntiLVsMasspt2);
-
-    fHistDcaV0DaughtersAntiLVsMasspt2     = new TH2F("h2DcaV0DaughtersAntiLVsMasspt2", "DCA between daughters;dca(cm);K0s inv. mass", 110, 0, 1.1,140,1.06,1.2);
-    fListHist->Add(fHistDcaV0DaughtersAntiLVsMasspt2);
-  
-    fHistCosPointAngleAntiLVsMasspt2      = new TH2F("h2CosPointAngleAntiLVsMasspt2", "Cosine of V0's pointing angle", 200,0.997,1.007,140,1.06,1.2);
-    fListHist->Add(fHistCosPointAngleAntiLVsMasspt2);
-
-    /// pt3
-    fHistDcaPosToPrimVertexAntiLVsMasspt3 = new TH2F("h2DcaPosToPrimVertexAntiLVsMasspt3", "Positive V0 daughter;dca(cm);K0s inv. mass",500,0,10,140,1.06,1.2);
-    fListHist->Add(fHistDcaPosToPrimVertexAntiLVsMasspt3);
-
-    fHistDcaNegToPrimVertexAntiLVsMasspt3 = new TH2F("h2DcaNegToPrimVertexAntiLVsMasspt3", "Negative V0 daughter;dca(cm);K0s inv. mass",500,0,10,140,1.06,1.2);
-    fListHist->Add(fHistDcaNegToPrimVertexAntiLVsMasspt3);
-
-    fHistRadiusV0AntiLVsMasspt3           = new TH2F("h2RadiusV0AntiLVsMasspt3", "Radius;Radius(cm);K0s inv. mass",110,0,110,140,1.06,1.2);
-    fListHist->Add(fHistRadiusV0AntiLVsMasspt3);
-
-    fHistDecayLengthV0AntiLVsMasspt3      = new TH2F("h2DecayLengthV0AntiLVsMasspt3", "V0s decay Length;decay length(cm);K0s inv. mass", 100, 0, 100,140,1.06,1.2);
-    fListHist->Add(fHistDecayLengthV0AntiLVsMasspt3);
-
-    fHistDcaV0DaughtersAntiLVsMasspt3     = new TH2F("h2DcaV0DaughtersAntiLVsMasspt3", "DCA between daughters;dca(cm);K0s inv. mass", 110, 0, 1.1,140,1.06,1.2);
-    fListHist->Add(fHistDcaV0DaughtersAntiLVsMasspt3);
-  
-    fHistCosPointAngleAntiLVsMasspt3      = new TH2F("h2CosPointAngleAntiLVsMasspt3", "Cosine of V0's pointing angle", 200,0.997,1.007,140,1.06,1.2);
-    fListHist->Add(fHistCosPointAngleAntiLVsMasspt3);
-  } //end QA condition for AntiLambda
-
   // V0 Multiplicity
   if (!fHistV0Multiplicity) {
     if (fCollidingSystems)
@@ -1461,22 +423,22 @@ void AliAnalysisTaskPerformanceStrange::UserCreateOutputObjects()
   fHistMassAntiLambda           = new TH1F("h1MassAntiLambda", "#bar{#Lambda}^{0} candidates;M(#bar{p}#pi^{+}) (GeV/c^{2});Counts", 150, 1.05, 1.2);
   fListHist->Add(fHistMassAntiLambda);
 
-  /*  // invariant mass vs radius
-      const Double_t radius[10] = {0.0,2.5,2.9,3.9,7.6,15.0,23.9,37.8,42.8,100.0};
-      Int_t lNbinRadius        = 9;
-      Int_t lNbinInvMassLambda = 300;
-  */
-  fHistMassVsRadiusK0           = new TH2F("h2MassVsRadiusK0", "K^{0} candidates;radius (cm);M(#pi^{+}#pi^{-}) (GeV/c^{2})",200,0,200, 200, 0.4, 0.6);
-  fListHist->Add(fHistMassVsRadiusK0);
+  fHistMassXi           = new TH1F("h1MassXi", "#Xi candidates;M(#bar{p}#pi^{+}) (GeV/c^{2});Counts", 150, 1.25, 1.4);  // 1321.71
+  fListHist->Add(fHistMassXi);
+
+  fHistMassAntiXi           = new TH1F("h1MassAntiXi", "#bar{#Xi} candidates;M(#bar{p}#pi^{+}) (GeV/c^{2});Counts", 150, 1.25, 1.4);
+  fListHist->Add(fHistMassAntiXi);
+
+  fHistMassOmega           = new TH1F("h1MassOmega", "#Omega candidates;M(#bar{p}#pi^{+}) (GeV/c^{2});Counts", 150, 1.6, 1.75); //1672.45
+  fListHist->Add(fHistMassOmega);
 
-  //fHistMassVsRadiusLambda       = new TH2F("h2MassVsRadiusLambda", "#Lambda candidates;radius (cm);M(p#pi^{-}) (GeV/c^{2})",lNbinRadius,radius, 140, 1.06, 1.2);
-  fHistMassVsRadiusLambda       = new TH2F("h2MassVsRadiusLambda", "#Lambda candidates;radius (cm);M(p#pi^{-}) (GeV/c^{2})",200,0,200, 140, 1.06, 1.2);
-  fListHist->Add(fHistMassVsRadiusLambda);
+  fHistMassAntiOmega           = new TH1F("h1MassAntiOmega", "#bar{#Omega} candidates;M(#bar{p}#pi^{+}) (GeV/c^{2});Counts", 150, 1.6, 1.75);
+  fListHist->Add(fHistMassAntiOmega);
 
+  //PID vs Mass
+  fHistMassXiVsPID           = new TH2F("h1MassXiVsPID", "#Xi candidates;M(#bar{p}#pi^{+}) (GeV/c^{2});PID", 150, 1.25, 1.4,5,0,5);  // 1321.71
+  fListHist->Add(fHistMassXiVsPID);
 
-  //fHistMassVsRadiusAntiLambda   = new TH2F("h2MassVsRadiusAntiLambda", "#bar{#Lambda} candidates;radius (cm);M(#bar{p}#pi^{+}) (GeV/c^{2})",lNbinRadius,radius, 140, 1.06, 1.2);
-  fHistMassVsRadiusAntiLambda   = new TH2F("h2MassVsRadiusAntiLambda", "#bar{#Lambda} candidates;radius (cm);M(#bar{p}#pi^{+}) (GeV/c^{2})",200,0,200, 140, 1.06, 1.2);
-  fListHist->Add(fHistMassVsRadiusAntiLambda);
 
 
   // Pt Vs Mass
@@ -1489,60 +451,6 @@ void AliAnalysisTaskPerformanceStrange::UserCreateOutputObjects()
   fHistPtVsMassAntiLambda           = new TH2F("h2PtVsMassAntiLambda","#AntiLambda^{0} candidates;M(p#pi^{-}) (GeV/c^{2});p_{t} (GeV/c)",280, 1.06, 1.2,240,0,12);
   fListHist->Add(fHistPtVsMassAntiLambda);
 
-
-  ////////// Transverse CTau distributions Vs Mass ////////////////
-  fHistTranscTauVsMassL                    = new TH2F("h2TranscTauVsMassL","cTauVsMass Lambdas",100,0,100,280, 1.06, 1.2);
-  fListHist->Add(fHistTranscTauVsMassL);
-
-  fHistTranscTauVsMassAntiL                    = new TH2F("h2TranscTauVsMassAntiL","cTauVsMass AntiLambdas",100,0,100,280, 1.06, 1.2);
-  fListHist->Add(fHistTranscTauVsMassAntiL);
-
-  fHistTranscTauVsMassK0s                    = new TH2F("h2TranscTauVsMassK0s","Transverse K0s cTau vs Mass ",100,0,100,400, 0.4, 0.6);
-  fListHist->Add(fHistTranscTauVsMassK0s);
-
-  // cTauVsMass Rap3
-  fHistTranscTauVsMassLRap3                    = new TH2F("h2TranscTauVsMassLRap3","cTauVsMass Lambdas",100,0,100,280, 1.06, 1.2);
-  fListHist->Add(fHistTranscTauVsMassLRap3);
-
-  fHistTranscTauVsMassAntiLRap3                    = new TH2F("h2TranscTauVsMassAntiLRap3","cTauVsMass AntiLambdas",100,0,100,280, 1.06, 1.2);
-  fListHist->Add(fHistTranscTauVsMassAntiLRap3);
-
-  fHistTranscTauVsMassK0sRap3                    = new TH2F("h2TranscTauVsMassK0sRap3","Transverse K0s cTau vs Mass ",100,0,100,400, 0.4, 0.6);
-  fListHist->Add(fHistTranscTauVsMassK0sRap3);
-
-  //cTauVsMass Low pt
-  fHistTranscTauVsMassLptLow                    = new TH2F("h2TranscTauVsMassLptLow","cTauVsMass Lambdas",100,0,100,280, 1.06, 1.2);
-  fListHist->Add(fHistTranscTauVsMassLptLow);
-
-  fHistTranscTauVsMassAntiLptLow                    = new TH2F("h2TranscTauVsMassAntiLptLow","cTauVsMass AntiLambdas",100,0,100,280, 1.06, 1.2);
-  fListHist->Add(fHistTranscTauVsMassAntiLptLow);
-
-  fHistTranscTauVsMassK0sptLow                    = new TH2F("h2TranscTauVsMassK0sptLow","Transverse K0s cTau vs Mass ",100,0,100,400, 0.4, 0.6);
-  fListHist->Add(fHistTranscTauVsMassK0sptLow);
-
-  //cTauVsMass Low pt Rap3
-  fHistTranscTauVsMassLptLowRap3                    = new TH2F("h2TranscTauVsMassLptLowRap3","cTauVsMass Lambdas",100,0,100,280, 1.06, 1.2);
-  fListHist->Add(fHistTranscTauVsMassLptLowRap3);
-
-  fHistTranscTauVsMassAntiLptLowRap3                    = new TH2F("h2TranscTauVsMassAntiLptLowRap3","cTauVsMass AntiLambdas",100,0,100,280, 1.06, 1.2);
-  fListHist->Add(fHistTranscTauVsMassAntiLptLowRap3);
-
-  fHistTranscTauVsMassK0sptLowRap3                    = new TH2F("h2TranscTauVsMassK0sptLowRap3","Transverse K0s cTau vs Mass ",100,0,100,400, 0.4, 0.6);
-  fListHist->Add(fHistTranscTauVsMassK0sptLowRap3);
-
-
-
-  // Pt Vs Mass Rap3
-  fHistPtVsMassK0Rap3               = new TH2F("h2PtVsMassK0Rap3","K^{0} candidates;M(#pi^{+}#pi^{-}) (GeV/c^{2});p_{t} (GeV/c)",400, 0.4, 0.6,240,0,12);
-  fListHist->Add(fHistPtVsMassK0Rap3);
-
-  fHistPtVsMassLambdaRap3           = new TH2F("h2PtVsMassLambdaRap3","#Lambda^{0} candidates;M(p#pi^{-}) (GeV/c^{2});p_{t} (GeV/c)",280, 1.06, 1.2,240,0,12);
-  fListHist->Add(fHistPtVsMassLambdaRap3);
-  
-  fHistPtVsMassAntiLambdaRap3           = new TH2F("h2PtVsMassAntiLambdaRap3","#AntiLambda^{0} candidates;M(p#pi^{-}) (GeV/c^{2});p_{t} (GeV/c)",280, 1.06, 1.2,240,0,12);
-  fListHist->Add(fHistPtVsMassAntiLambdaRap3);
   
   ////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
 
@@ -1567,32 +475,6 @@ void AliAnalysisTaskPerformanceStrange::UserCreateOutputObjects()
   fListHist->Add(fHistNSigmaProton);
 
 
-  //PID
-  fHistNsigmaPosPionAntiLambda   = new TH1F("h1NsigmaPosPionAntiLambda", "Positive daughter of Antilambda;NsigmaPion;Counts",25,0,5);
-  fListHist->Add(fHistNsigmaPosPionAntiLambda);
-
-  fHistNsigmaNegProtonAntiLambda = new TH1F("h1NsigmaNegProtonAntiLambda", "Negative daughter of Antilambda;NsigmaProton;Counts",25,0,5);
-  fListHist->Add(fHistNsigmaNegProtonAntiLambda);
-  
-  fHistNsigmaPosProtonLambda     = new TH1F("h1NsigmaPosProtonLambda", "Positive daughter of Lambda;NsigmaProton;Counts",25,0,5); 
-  fListHist->Add(fHistNsigmaPosProtonLambda);
-  
-  fHistNsigmaNegPionLambda       = new TH1F("h1NsigmaNegPionLambda", "Negative daughter of Lambda;NsigmaPion;Counts",25,0,5);
-  fListHist->Add(fHistNsigmaNegPionLambda);
-  
-  fHistNsigmaPosProtonAntiLambda     = new TH1F("h1NsigmaPosProtonAntiLambda", "Positive daughter of AntiLambda;NsigmaProton;Counts",25,0,5); 
-  fListHist->Add(fHistNsigmaPosProtonAntiLambda);
-  
-  fHistNsigmaNegPionAntiLambda       = new TH1F("h1NsigmaNegPionAntiLambda", "Negative daughter of AntiLambda;NsigmaPion;Counts",25,0,5);
-  fListHist->Add(fHistNsigmaNegPionAntiLambda);
-  
-  fHistNsigmaPosPionK0           = new TH1F("h1NsigmaPosPionK0", "Positive daughter of K0s;NsigmaPion;Counts",25,0,5);
-  fListHist->Add(fHistNsigmaPosPionK0);
-  
-  fHistNsigmaNegPionK0           = new TH1F("h1NsigmaNegPionK0", "Negative daughter of K0s;NsigmaPion;Counts",25,0,5);
-  fListHist->Add(fHistNsigmaNegPionK0);
-
-
   //********************************
   // Associated particles histograms
   //********************************
@@ -1617,82 +499,6 @@ void AliAnalysisTaskPerformanceStrange::UserCreateOutputObjects()
   fHistAsMcPtAntiLambda            = new TH1F("h1AsMcPtAntiLambda", "#AntiLambda^{0} associated;p_{t} (GeV/c);Counts", 240,0,12);
   fListHist->Add(fHistAsMcPtAntiLambda);
 
-
-  fHistAsMcPtZoomK0            = new TH1F("h1AsMcPtZoomK0", "K^{0} candidates in -1 <y<1;p_{t} (GeV/c);Counts",20,0,1);
-  fListHist->Add(fHistAsMcPtZoomK0);
-
-  fHistAsMcPtZoomLambda        = new TH1F("h1AsMcPtZoomLambda", "#Lambda^{0} candidates in -1 <y<1;p_{t} (GeV/c);Counts",20,0,1);
-  fListHist->Add(fHistAsMcPtZoomLambda);
-
-  fHistAsMcPtZoomAntiLambda        = new TH1F("h1AsMcPtZoomAntiLambda", "#AntiLambda^{0} candidates in -1 <y<1;p_{t} (GeV/c);Counts",20,0,1);
-  fListHist->Add(fHistAsMcPtZoomAntiLambda);
-
-  //Pt distribution Rap3
-  fHistAsMcPtK0Rap3                = new TH1F("h1AsMcPtK0Rap3", "K^{0} associated;p_{t} (GeV/c);Counts", 240,0,12);
-  fListHist->Add(fHistAsMcPtK0Rap3);
-
-  fHistAsMcPtLambdaRap3            = new TH1F("h1AsMcPtLambdaRap3", "#Lambda^{0} associated;p_{t} (GeV/c);Counts", 240,0,12);
-  fListHist->Add(fHistAsMcPtLambdaRap3);
-
-  fHistAsMcPtAntiLambdaRap3            = new TH1F("h1AsMcPtAntiLambdaRap3", "#AntiLambda^{0} associated;p_{t} (GeV/c);Counts", 240,0,12);
-  fListHist->Add(fHistAsMcPtAntiLambdaRap3);
-
-  //cTaud istribution
-  fHistAsMcTranscTauL                    = new TH1F("h1AsMcTranscTauL","cTau of AntiLambdas",100,0,100);
-  fListHist->Add(fHistAsMcTranscTauL);
-  fHistAsMcTranscTauAntiL                    = new TH1F("h1AsMcTranscTauAntiL","cTau of AntiLambdas",100,0,100);
-  fListHist->Add(fHistAsMcTranscTauAntiL);
-  fHistAsMcTranscTauK0s                    = new TH1F("h1AsMcTranscTauK0s","cTau of K0s",100,0,100);
-  fListHist->Add(fHistAsMcTranscTauK0s);
-
-  //cTau Rap3
-  fHistAsMcTranscTauLRap3                    = new TH1F("h1AsMcTranscTauLRap3","cTau of AntiLambdas",100,0,100);
-  fListHist->Add(fHistAsMcTranscTauLRap3);
-  fHistAsMcTranscTauAntiLRap3                    = new TH1F("h1AsMcTranscTauAntiLRap3","cTau of AntiLambdas",100,0,100);
-  fListHist->Add(fHistAsMcTranscTauAntiLRap3);
-  fHistAsMcTranscTauK0sRap3                    = new TH1F("h1AsMcTranscTauK0sRap3","cTau of K0s",100,0,100);
-  fListHist->Add(fHistAsMcTranscTauK0sRap3);
-
-
-  //cTau lowpt
-  fHistAsMcTranscTauLptLow                    = new TH1F("h1AsMcTranscTauLptLow","cTau of AntiLambdas",100,0,100);
-  fListHist->Add(fHistAsMcTranscTauLptLow);
-  fHistAsMcTranscTauAntiLptLow                    = new TH1F("h1AsMcTranscTauAntiLptLow","cTau of AntiLambdas",100,0,100);
-  fListHist->Add(fHistAsMcTranscTauAntiLptLow);
-  fHistAsMcTranscTauK0sptLow                    = new TH1F("h1AsMcTranscTauK0sptLow","cTau of K0s",100,0,100);
-  fListHist->Add(fHistAsMcTranscTauK0sptLow);
-
-  //cTau ptLowRap3 Low pt
-  fHistAsMcTranscTauLptLowRap3                    = new TH1F("h1AsMcTranscTauLptLowRap3","cTau of AntiLambdas",100,0,100);
-  fListHist->Add(fHistAsMcTranscTauLptLowRap3);
-  fHistAsMcTranscTauAntiLptLowRap3                    = new TH1F("h1AsMcTranscTauAntiLptLowRap3","cTau of AntiLambdas",100,0,100);
-  fListHist->Add(fHistAsMcTranscTauAntiLptLowRap3);
-  fHistAsMcTranscTauK0sptLowRap3                    = new TH1F("h1AsMcTranscTauK0sptLowRap3","cTau of K0s",100,0,100);
-  fListHist->Add(fHistAsMcTranscTauK0sptLowRap3);
-
-
-
-
-  // Radius distribution
-  fHistAsMcProdRadiusK0               = new TH1F("h1AsMcProdRadiusK0", "K^{0} associated;r (cm);Counts", 500, 0, 100);
-  fListHist->Add(fHistAsMcProdRadiusK0);
-
-  fHistAsMcProdRadiusLambda           = new TH1F("h1AsMcProdRadiusLambda", "#Lambda^{0} associated;r (cm);Counts", 500, 0, 100);
-  fListHist->Add(fHistAsMcProdRadiusLambda);
-
-  fHistAsMcProdRadiusAntiLambda       = new TH1F("h1AsMcProdRadiusAntiLambda", "#bar{#Lambda}^{0} associated;r (cm);Counts", 500, 0, 100);
-  fListHist->Add(fHistAsMcProdRadiusAntiLambda);
-
-  fHistAsMcProdRadiusXvsYK0s          = new TH2F("h2AsMcProdRadiusXvsYK0s","Associated Secondary K^{0} Production Radius;x (cm); y (cm)",200,-50,50,200,-50,50);
-  fListHist->Add(fHistAsMcProdRadiusXvsYK0s);
-
-  fHistAsMcProdRadiusXvsYLambda       = new TH2F("h2AsMcProdRadiusXvsYLambda","Associated Secondary #Lambda^{0} Production Radius;x (cm); y (cm)",200,-50,50,200,-50,50);
-  fListHist->Add(fHistAsMcProdRadiusXvsYLambda);
-
-  fHistAsMcProdRadiusXvsYAntiLambda   = new TH2F("h2AsMcProdRadiusXvsYAntiLambda","Associated Secondary #bar{#Lambda}^{0} Production Radius;x (cm); y (cm)",200,-50,50,200,-50,50);
-  fListHist->Add(fHistAsMcProdRadiusXvsYAntiLambda);
-
   // Mass
   fHistPidMcMassK0             = new TH1F("h1PidMcMassK0", "K^{0} MC PId checked;M(#pi^{+}#pi^{-}) (GeV/c^{2});Counts", 100, 0.4, 0.6);
   fListHist->Add(fHistPidMcMassK0);
@@ -1712,19 +518,6 @@ void AliAnalysisTaskPerformanceStrange::UserCreateOutputObjects()
   fHistAsMcMassAntiLambda      = new TH1F("h1AsMcMassAntiLambda", "#bar{#Lambda}^{0} associated;M(#bar{p}#pi^{+}) (GeV/c^{2});Counts", 75, 1.05, 1.2);
   fListHist->Add(fHistAsMcMassAntiLambda);
 
-
-  //Mass Rap3
-  fHistAsMcMassK0Rap3              = new TH1F("h1AsMcMassK0Rap3", "K^{0} associated;M(#pi^{+}#pi^{-}) (GeV/c^{2});Counts", 100, 0.4, 0.6);
-  fListHist->Add(fHistAsMcMassK0Rap3);
-  
-  fHistAsMcMassLambdaRap3          = new TH1F("h1AsMcMassLambdaRap3", "#Lambda^{0} associated;M(p#pi^{-}) (GeV/c^{2});Counts", 75, 1.05, 1.2);
-  fListHist->Add(fHistAsMcMassLambdaRap3);
-
-  fHistAsMcMassAntiLambdaRap3      = new TH1F("h1AsMcMassAntiLambdaRap3", "#bar{#Lambda}^{0} associated;M(#bar{p}#pi^{+}) (GeV/c^{2});Counts", 75, 1.05, 1.2);
-  fListHist->Add(fHistAsMcMassAntiLambdaRap3);
-
-
-
   //Pt versus Mass
   fHistAsMcPtVsMassK0               = new TH2F("h2AsMcPtVsMassK0","K^{0} associated;M(#pi^{+}#pi^{-}) (GeV/c^{2});p_{t} (GeV/c)",200, 0.4, 0.6,240,0,12);
   fListHist->Add(fHistAsMcPtVsMassK0);
@@ -1735,159 +528,21 @@ void AliAnalysisTaskPerformanceStrange::UserCreateOutputObjects()
   fHistAsMcPtVsMassAntiLambda       = new TH2F("h2AsMcPtVsMassAntiLambda","#bar{#Lambda}^{0} associated;M(#bar{p}#pi^{+}) (GeV/c^{2});p_{t} (GeV/c)",140, 1.06, 1.2,240,0,12);
   fListHist->Add(fHistAsMcPtVsMassAntiLambda);
 
+  fHistCompositionXi = new TH1F("h1CompisitionXi","background composition of Xi;Part ID;PID",8000,-4000,4000);
+  fListHist->Add(fHistCompositionXi);
 
-  //Pt versus Mass Rap3
-  fHistAsMcPtVsMassK0Rap3               = new TH2F("h2AsMcPtVsMassK0Rap3","K^{0} associated;M(#pi^{+}#pi^{-}) (GeV/c^{2});p_{t} (GeV/c)",200, 0.4, 0.6,240,0,12);
-  fListHist->Add(fHistAsMcPtVsMassK0Rap3);
-
-  fHistAsMcPtVsMassLambdaRap3           = new TH2F("h2AsMcPtVsMassLambdaRap3","#Lambda^{0} associated;M(p#pi^{-}) (GeV/c^{2});p_{t} (GeV/c)",140, 1.06, 1.2,240,0,12);
-  fListHist->Add(fHistAsMcPtVsMassLambdaRap3);
-
-  fHistAsMcPtVsMassAntiLambdaRap3       = new TH2F("h2AsMcPtVsMassAntiLambdaRap3","#bar{#Lambda}^{0} associated;M(#bar{p}#pi^{+}) (GeV/c^{2});p_{t} (GeV/c)",140, 1.06, 1.2,240,0,12);
-  fListHist->Add(fHistAsMcPtVsMassAntiLambdaRap3);
-
-
-
-
-
-  // invariant mass vs radius
-  //fHistAsMcMassVsRadiusK0             = new TH2F("h2AsMcMassVsRadiusK0", "K^{0} associated;radius (cm);M(#pi^{+}#pi^{-}) (GeV/c^{2})",lNbinRadius,radius, 500, 0.47, 0.52);
-  fHistAsMcMassVsRadiusK0             = new TH2F("h2AsMcMassVsRadiusK0", "K^{0} associated;radius (cm);M(#pi^{+}#pi^{-}) (GeV/c^{2})",200,0,200, 500, 0.47, 0.52);
-  fListHist->Add(fHistAsMcMassVsRadiusK0);
-
-  //fHistAsMcMassVsRadiusLambda         = new TH2F("h2AsMcMassVsRadiusLambda", "#Lambda associated;radius (cm);M(p#pi^{-}) (GeV/c^{2})",lNbinRadius,radius, lNbinInvMassLambda, 1.10, 1.13);
-  fHistAsMcMassVsRadiusLambda         = new TH2F("h2AsMcMassVsRadiusLambda", "#Lambda associated;radius (cm);M(p#pi^{-}) (GeV/c^{2})",200,0,200, 1.10, 1.13);
-  fListHist->Add(fHistAsMcMassVsRadiusLambda);
-
-  //fHistAsMcMassVsRadiusAntiLambda     = new TH2F("h2AsMcMassVsRadiusAntiLambda", "#bar{#Lambda} associated;radius (cm);M(#bar{p}#pi^{+}) (GeV/c^{2})",lNbinRadius,radius,lNbinInvMassLambda , 1.10, 1.13);
-  fHistAsMcMassVsRadiusAntiLambda     = new TH2F("h2AsMcMassVsRadiusAntiLambda", "#bar{#Lambda} associated;radius (cm);M(#bar{p}#pi^{+}) (GeV/c^{2})",200,0,200 , 1.10, 1.13);
-  fListHist->Add(fHistAsMcMassVsRadiusAntiLambda);
-  
-  // Position Resolution
-  fHistAsMcResxK0                     = new TH1F("h1AsMcResxK0", "K^{0} associated;#Delta x (cm);Counts", 50, -0.25, 0.25);
-  fListHist->Add(fHistAsMcResxK0);
-  fHistAsMcResyK0                     = new TH1F("h1AsMcResyK0", "K^{0} associated;#Delta y (cm);Counts", 50, -0.25, 0.25);
-  fListHist->Add(fHistAsMcResyK0);
-  fHistAsMcReszK0                     = new TH1F("h1AsMcReszK0", "K^{0} associated;#Delta z (cm);Counts", 50, -0.25, 0.25);
-  fListHist->Add(fHistAsMcReszK0);
-  fHistAsMcResrVsRadiusK0             = new TH2F("h2AsMcResrVsRadiusK0", "K^{0} associated;Radius (cm);#Delta r (cm)",200,0.0,50., 50, -0.25, 0.25);
-  fListHist->Add(fHistAsMcResrVsRadiusK0);
-  fHistAsMcReszVsRadiusK0             = new TH2F("h2AsMcReszVsRadiusK0", "K^{0} associated;Radius (cm);#Delta z (cm)",200,0.0,50.0, 50, -0.25, 0.25);
-  fListHist->Add(fHistAsMcReszVsRadiusK0);
-
-  fHistAsMcResxLambda                 = new TH1F("h1AsMcResxLambda", "#Lambda^{0} associated;#Delta x (cm);Counts", 50, -0.25, 0.25);
-  fListHist->Add(fHistAsMcResxLambda);
-  fHistAsMcResyLambda                 = new TH1F("h1AsMcResyLambda", "#Lambda^{0} associated;#Delta y (cm);Counts", 50, -0.25, 0.25);
-  fListHist->Add(fHistAsMcResyLambda);
-  fHistAsMcReszLambda                 = new TH1F("h1AsMcReszLambda", "#Lambda^{0} associated;#Delta z (cm);Counts", 50, -0.25, 0.25);
-  fListHist->Add(fHistAsMcReszLambda);
-  fHistAsMcResrVsRadiusLambda         = new TH2F("h2AsMcResrVsRadiusLambda", "#Lambda^{0} associated;Radius (cm);#Delta r (cm)",200,0.0,50.0, 50, -0.25, 0.25);
-  fListHist->Add(fHistAsMcResrVsRadiusLambda);
-  fHistAsMcReszVsRadiusLambda         = new TH2F("h2AsMcReszVsRadiusLambda", "#Lambda^{0} associated;Radius (cm);#Delta z (cm)",200,0.0,50.0, 50, -0.25, 0.25);
-  fListHist->Add(fHistAsMcReszVsRadiusLambda);
-
-  fHistAsMcResxAntiLambda             = new TH1F("h1AsMcResxAntiLambda", "#bar{#Lambda}^{0} associated;#Delta x (cm);Counts", 50, -0.25, 0.25);
-  fListHist->Add(fHistAsMcResxAntiLambda);
-  fHistAsMcResyAntiLambda             = new TH1F("h1AsMcResyAntiLambda", "#bar{#Lambda}^{0} associated;#Delta y (cm);Counts", 50, -0.25, 0.25);
-  fListHist->Add(fHistAsMcResyAntiLambda);
-  fHistAsMcReszAntiLambda             = new TH1F("h1AsMcReszAntiLambda", "#bar{#Lambda}^{0} associated;#Delta z (cm);Counts", 50, -0.25, 0.25);
-  fListHist->Add(fHistAsMcReszAntiLambda);
-  fHistAsMcResrVsRadiusAntiLambda     = new TH2F("h2AsMcResrVsRadiusAntiLambda", "#bar{#Lambda}^{0} associated;Radius (cm);#Delta r (cm)",200,0.0,50.0, 50, -0.25, 0.25);
-  fListHist->Add(fHistAsMcResrVsRadiusAntiLambda);
-  fHistAsMcReszVsRadiusAntiLambda     = new TH2F("h2AsMcReszVsRadiusAntiLambda", "#bar{#Lambda}^{0} associated;Radius (cm);#Delta z (cm)",200,0.0,50.0, 50, -0.25, 0.25);
-  fListHist->Add(fHistAsMcReszVsRadiusAntiLambda);
-
-  // Pt Resolution
-  fHistAsMcResPtK0                   = new TH1F("h1AsMcResPtK0","Pt Resolution K^{0};#Delta Pt;Counts",200,-1,1);
-  fListHist->Add(fHistAsMcResPtK0);
-  
-  fHistAsMcResPtLambda               = new TH1F("h1AsMcResPtLambda","Pt Resolution #Lambda^{0};#Delta Pt;Counts",200,-1,1);
-  fListHist->Add(fHistAsMcResPtLambda);
-
-  fHistAsMcResPtAntiLambda           = new TH1F("h1AsMcResPtAntiLambda","Pt Resolution #bar{#Lambda}^{0};#Delta Pt;Counts",200,-1,1);
-  fListHist->Add(fHistAsMcResPtAntiLambda);
-
-
-  fHistAsMcResPtVsRapK0              = new TH2F("h2AsMcResPtVsRapK0","Pt Resolution K^{0};#Delta Pt;Rap",200,-1,1,20,-1,1);
-  fListHist->Add(fHistAsMcResPtVsRapK0);
-  
-  fHistAsMcResPtVsRapLambda          = new TH2F("h2AsMcResPtVsRapLambda","Pt Resolution #Lambda^{0};#Delta Pt;Rap",200,-1,1,20,-1,1);
-  fListHist->Add(fHistAsMcResPtVsRapLambda);
-
-  fHistAsMcResPtVsRapAntiLambda      = new TH2F("h2AsMcResPtVsRapAntiLambda","Pt Resolution #bar{#Lambda}^{0};#Delta Pt;Rap",200,-1,1,20,-1,1);
-  fListHist->Add(fHistAsMcResPtVsRapAntiLambda);
-
-  fHistAsMcResPtVsPtK0               = new TH2F("h2AsMcResPtVsPtK0","Pt Resolution K^{0};#Delta Pt;Pt",600,-0.15,0.15,240,0,12);
-  fListHist->Add(fHistAsMcResPtVsPtK0);
-    
-  fHistAsMcResPtVsPtLambda           = new TH2F("h2AsMcResPtVsPtLambda","Pt Resolution #Lambda^{0};#Delta Pt;Pt",600,-0.15,0.15,240,0,12);
-  fListHist->Add(fHistAsMcResPtVsPtLambda);
-
-  fHistAsMcResPtVsPtAntiLambda       = new TH2F("h2AsMcResPtVsPtAntiLambda","Pt Resolution #bar{#Lambda}^{0};#Delta Pt;Pt",300,-0.15,0.15,240,0,12);
-  fListHist->Add(fHistAsMcResPtVsPtAntiLambda);
-
-  // pdgcode of mother
-  fHistAsMcMotherPdgCodeK0s           = new TH1F("h1AsMcMotherPdgCodeK0s","Mother of Associated K^{0};mother;counts",11,0,11);
-  fListHist->Add(fHistAsMcMotherPdgCodeK0s);
-
-  fHistAsMcMotherPdgCodeLambda        = new TH1F("h1AsMcMotherPdgCodeLambda","Mother of Associated #Lambda^{0};mother;counts",11,0,11);
-  fListHist->Add(fHistAsMcMotherPdgCodeLambda);
-
-  fHistAsMcMotherPdgCodeAntiLambda    = new TH1F("h1AsMcMotherPdgCodeAntiLambda","Mother of Associated #bar{#Lambda}^{0};mother;counts",11,0,11);
-  fListHist->Add(fHistAsMcMotherPdgCodeAntiLambda);
-
-  // Pt distribution Lambda from Sigma
-  fHistAsMcPtLambdaFromSigma          = new TH1F("h1AsMcPtLambdaFromSigma","#Lambda}^{0} associated from Sigma;p_{t} (GeV/c);Count",240,0,12);
-  fListHist->Add(fHistAsMcPtLambdaFromSigma);
-
-  fHistAsMcPtAntiLambdaFromSigma      = new TH1F("h1AsMcPtAntiLambdaFromSigma","#bar{#Lambda}^{0} associated from Sigma;p_{t} (GeV/c);Count",240,0,12);
-  fListHist->Add(fHistAsMcPtAntiLambdaFromSigma);
-
-  // Associated secondary particles:
-  // Pt and rapidity distribution
-  fHistAsMcSecondaryPtVsRapK0s          = new TH2F("h2AsMcSecondaryPtVsRapK0s", "K^{0} associated secondary;p_{t} (GeV/c);rapidity",240,0,12,30,-1.5,1.5);
-  fListHist->Add(fHistAsMcSecondaryPtVsRapK0s);
-
-  fHistAsMcSecondaryPtVsRapLambda       = new TH2F("h2AsMcSecondaryPtVsRapLambda", "#Lambda^{0} associated secondary;p_{t} (GeV/c);rapidity",240,0,12,30,-1.5,1.5);
-  fListHist->Add(fHistAsMcSecondaryPtVsRapLambda);
-
-  fHistAsMcSecondaryPtVsRapAntiLambda   = new TH2F("h2AsMcSecondaryPtVsRapAntiLambda", "#bar{#Lambda}^{0} associated secondary;p_{t} (GeV/c);rapidity",240,0,12,30,-1.5,1.5);
-  fListHist->Add(fHistAsMcSecondaryPtVsRapAntiLambda);
-
-  // Production radius
-  fHistAsMcSecondaryProdRadiusK0s              = new TH1F("h1AsMcSecondaryProdRadiusK0s", "K^{0} Production Radius;r (cm);Count", 170, -2, 15);
-  fListHist->Add(fHistAsMcSecondaryProdRadiusK0s);
-
-  fHistAsMcSecondaryProdRadiusLambda           = new TH1F("h1AsMcSecondaryProdRadiusLambda", "#Lambda^{0} Production Radius;r (cm);Count", 170, -2, 15);
-  fListHist->Add(fHistAsMcSecondaryProdRadiusLambda);
-
-  fHistAsMcSecondaryProdRadiusAntiLambda       = new TH1F("h1AsMcSecondaryProdRadiusAntiLambda", "#bar{#Lambda}^{0} Production Radius;r (cm);Count", 170, -2, 15);
-  fListHist->Add(fHistAsMcSecondaryProdRadiusAntiLambda);  
-
-  fHistAsMcSecondaryProdRadiusXvsYK0s          = new TH2F("h2AsMcSecondaryProdRadiusXvsYK0s","Associated Secondary K^{0} Production Radius;x (cm); y (cm)",200,-20,20,200,-20,20);
-  fListHist->Add(fHistAsMcSecondaryProdRadiusXvsYK0s);
-
-  fHistAsMcSecondaryProdRadiusXvsYLambda       = new TH2F("h2AsMcSecondaryProdRadiusXvsYLambda","Associated Secondary #Lambda^{0} Production Radius;x (cm); y (cm)",200,-20,20,200,-20,20);
-  fListHist->Add(fHistAsMcSecondaryProdRadiusXvsYLambda);
+  fHistCompositionAntiXi = new TH1F("h1CompisitionAntiXi","background composition of AntiXi;Part ID;Counts",8000,-4000,4000);
+  fListHist->Add(fHistCompositionAntiXi);
 
-  fHistAsMcSecondaryProdRadiusXvsYAntiLambda   = new TH2F("h2AsMcSecondaryProdRadiusXvsYAntiLambda","Associated Secondary #bar{#Lambda}^{0} Production Radius;x (cm); y (cm)",200,-20,20,200,-20,20);
-  fListHist->Add(fHistAsMcSecondaryProdRadiusXvsYAntiLambda);
+  fHistCompositionOmega = new TH1F("h1CompisitionOmega","background composition of Omega;Part ID;Counts",8000,-4000,4000);
+  fListHist->Add(fHistCompositionOmega);
 
-  fHistAsMcSecondaryMotherPdgCodeK0s           = new TH1F("h1AsMcSecondaryMotherPdgCodeK0s","Mother of Associated Secondary K^{0};mother;counts",11,0,11);
-  fListHist->Add(fHistAsMcSecondaryMotherPdgCodeK0s);
+  fHistCompositionAntiOmega = new TH1F("h1CompisitionAntiOmega","background composition of AntiOmega;Part ID;Counts",8000,-4000,4000);
+  fListHist->Add(fHistCompositionAntiOmega);
 
-  fHistAsMcSecondaryMotherPdgCodeLambda        = new TH1F("h1AsMcSecondaryMotherPdgCodeLambda","Mother of Associated Secondary #Lambda^{0};mother;counts",11,0,11);
-  fListHist->Add(fHistAsMcSecondaryMotherPdgCodeLambda);
+  fHistMCIndexes = new TH1I("h1MCIndexes","MC Indexes;Index;Counts",210,-10,200);
+  fListHist->Add(fHistMCIndexes);
 
-  fHistAsMcSecondaryMotherPdgCodeAntiLambda    = new TH1F("h1AsMcSecondaryMotherPdgCodeAntiLambda","Mother of Associated Secondary #bar{#Lambda}^{0};mother;counts",11,0,11);
-  fListHist->Add(fHistAsMcSecondaryMotherPdgCodeAntiLambda);
-
-  // Pt distribution Lambda from Sigma
-  fHistAsMcSecondaryPtLambdaFromSigma          = new TH1F("h1AsMcSecondaryPtLambdaFromSigma","#Lambda}^{0} associated from Sigma;p_{t} (GeV/c);Count",240,0,12);
-  fListHist->Add(fHistAsMcSecondaryPtLambdaFromSigma);
-
-  fHistAsMcSecondaryPtAntiLambdaFromSigma      = new TH1F("h1AsMcSecondaryPtAntiLambdaFromSigma","#bar{#Lambda}^{0} associated from Sigma;p_{t} (GeV/c);Count",240,0,12);
-  fListHist->Add(fHistAsMcSecondaryPtAntiLambdaFromSigma);
   PostData(1, fListHist);
   PostData(2, fCentrSelector);
   PostData(3, fTracksCuts);
@@ -1928,54 +583,15 @@ void AliAnalysisTaskPerformanceStrange::UserExec(Option_t *)
   }
 
   fHistNumberEvents->Fill(0.5);
-
-  //******************
-  // Trigger Selection ! Warning Works only for ESD, add protection in case of AOD loop
-  //******************
-  Bool_t isSelected = kFALSE;
-  if(fCollidingSystems == 2)
-        isSelected = (((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected()& AliVEvent::kCINT5 && ((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected()& AliVEvent::kINT7); 
-    //    isSelected = (((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected()& AliVEvent::kCINT7);
-  else
-   isSelected = (((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected()& AliVEvent::kMB);
-
-  if (!isSelected) return;
-
   fHistNumberEvents->Fill(1.5);  
-
-
   // Centrality selection
 
   static AliESDtrackCuts * trackCuts = AliESDtrackCuts::GetStandardITSTPCTrackCuts2010(kTRUE,1); 
   Bool_t isCentralitySelected = fCentrSelector->IsCentralityBinSelected(fESD,trackCuts);
   if(!isCentralitySelected) return;
-  
-  // Done by the AliPhysicsSelection Task ! Only the selected events are passed to this task
-  
-  AliESDtrackCuts *TestTrackCuts = new AliESDtrackCuts();
-  TestTrackCuts->SetRequireTPCRefit(kTRUE);
-  TestTrackCuts->SetRequireITSRefit(kFALSE);
-  // mytracksCuts->SetMinNClustersTPC(nbMinTPCclusters);                                                                      
-  TestTrackCuts->SetMinNCrossedRowsTPC(70);
-  TestTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.8);
-
-  AliESDtrack *track = 0x0;
-  Double_t trackP[3] = {0,0,0};
-  for (int i = 0; i<fESD->GetNumberOfTracks();i++){
-    track = fESD->GetTrack(i);
-    track->GetInnerPxPyPz(trackP);
-    TestTrackCuts->SetRequireITSRefit(kFALSE);
-    if(TestTrackCuts->IsSelected(track)) fHistPtTracks->Fill(trackP[2]);
-    TestTrackCuts->SetRequireITSRefit(kTRUE);
-    if(TestTrackCuts->IsSelected(track)) fHistPtTracksITSRefit->Fill(trackP[2]);
-  }
-  
 
   // Remove Events with no tracks
-  //if (!(fESD->GetNumberOfTracks()))  return;
-
-  fHistNumberEvents->Fill(2.5);
-  //  fHistTrackPerEvent->Fill(fESD->GetNumberOfTracks());
+  if (!(fESD->GetNumberOfTracks()))  return;
 
   //*************************************
   // Cut used:
@@ -1988,19 +604,14 @@ void AliAnalysisTaskPerformanceStrange::UserExec(Option_t *)
 
   Double_t lCutPseudorap = 0.8;
 
-  // Cut AliKF Chi2 for Reconstructed particles
-  //  Double_t cutChi2KF  = 1E3;
-
   // If PID is used:
   Double_t lLimitPPID    = 0.7;
   Float_t cutNSigmaLowP  = 1E3;
   Float_t cutNSigmaHighP = 1E3;
   if (fUsePID.Contains("withPID")) {
-    cutNSigmaLowP  = 3.0;
-    cutNSigmaHighP = 3.0;
+    cutNSigmaLowP  = 4.0;
+    cutNSigmaHighP = 4.0;
   }
-
-
   // Cut Daughters pt (GeV/c):
   //Double_t cutMinPtDaughter = 0.160;
 
@@ -2013,13 +624,9 @@ void AliAnalysisTaskPerformanceStrange::UserExec(Option_t *)
 
   // Min number of TPC clusters:
   // Int_t nbMinTPCclusters = 80;
-
-  
-  //
   // PID flags:
   Int_t LambdaPID = 0;
   Int_t AntiLambdaPID = 0;
-      
 
   //////******************************************
   ////// Access MC: ******************************
@@ -2060,7 +667,6 @@ void AliAnalysisTaskPerformanceStrange::UserExec(Option_t *)
     */
   }
 
-
   //**********************************************
   // MC loop
   //**********************************************
@@ -2075,7 +681,6 @@ void AliAnalysisTaskPerformanceStrange::UserExec(Option_t *)
   Double_t lPtCurrentPart   = 0;
   
   Int_t lComeFromSigma      = 0;
-
   
   // Production Radius
   Double_t mcPosX     = 0.0,  mcPosY      = 0.0,  mcPosZ      = 0.0;
@@ -2088,40 +693,32 @@ void AliAnalysisTaskPerformanceStrange::UserExec(Option_t *)
   Int_t iCurrentMother  = 0, lPdgCurrentMother    = 0;
   //  Bool_t lCurrentMotherIsPrimary;
 
-  // current mc particles 's daughter:
-  //  Int_t lPdgCurrentDaughter0 = 0, lPdgCurrentDaughter1 = 0; 
-
   // variables for multiple reconstruction studies:
   Int_t id0           = 0, id1          = 0;
-  //Int_t lLabelTrackN  = 0, lLabelTrackP = 0;
-  //Int_t lPartNMother  = 0, lPartPMother = 0;
-  //Int_t lPartPMotherPDGcode      = 0;
   Int_t lNtimesReconstructedK0s   = 0, lNtimesReconstructedLambda   = 0, lNtimesReconstructedAntiLambda   = 0;
   // Int_t lNtimesReconstructedK0sMI = 0, lNtimesReconstructedLambdaMI = 0, lNtimesReconstructedAntiLambdaMI = 0;
 
   ////********************************
   ////Start loop over MC particles****
   ////********************************
-
+  
   if (fAnalysisMC) {
 
     // Primary vertex
     fHistMCPrimaryVertexX->Fill(mcPrimaryVtx.At(0));
     fHistMCPrimaryVertexY->Fill(mcPrimaryVtx.At(1));
     fHistMCPrimaryVertexZ->Fill(mcPrimaryVtx.At(2));
-    
-    //    lmcPrimVtxR = TMath::Sqrt(mcPrimaryVtx.At(0)*mcPrimaryVtx.At(0)+mcPrimaryVtx.At(1)*mcPrimaryVtx.At(1));
   
+    //    lmcPrimVtxR = TMath::Sqrt(mcPrimaryVtx.At(0)*mcPrimaryVtx.At(0)+mcPrimaryVtx.At(1)*mcPrimaryVtx.At(1));
 
     if(fAnalysisType == "ESD") {
-      
-      lNbMCPrimary = stack->GetNprimary();
-      lNbMCPart    = stack->GetNtrack();
+    
+      lNbMCPrimary = stack->GetNprimary(); 
+     lNbMCPart    = stack->GetNtrack();
       
       fHistMCMultiplicityPrimary->Fill(lNbMCPrimary);
       fHistMCMultiplicityTracks->Fill(lNbMCPart);
       
-      
       for (Int_t iMc = 0; iMc < (stack->GetNtrack()); iMc++) {  
        TParticle *p0 = stack->Particle(iMc);
        if (!p0) {
@@ -2130,8 +727,8 @@ void AliAnalysisTaskPerformanceStrange::UserExec(Option_t *)
        }
        lPdgcodeCurrentPart = p0->GetPdgCode();
        
-       // Keep only K0s, Lambda and AntiLambda, Xi and Phi:
-       if ( (lPdgcodeCurrentPart != 310 ) && (lPdgcodeCurrentPart != 3122 ) && (lPdgcodeCurrentPart != -3122 ) && (lPdgcodeCurrentPart != 3312 ) && (lPdgcodeCurrentPart != -3312) && (lPdgcodeCurrentPart != -333) ) continue;
+       // Keep only K0s, Lambda and AntiLambda, Xi-, antiXi+ and Omega-, antiOmega+:
+       if ( (lPdgcodeCurrentPart != 310 ) && (lPdgcodeCurrentPart != 3122 ) && (lPdgcodeCurrentPart != -3122 ) && (lPdgcodeCurrentPart != 3312 ) && (lPdgcodeCurrentPart != -3312) && (lPdgcodeCurrentPart != 3334)  && (lPdgcodeCurrentPart != -3334) ) continue;
        
        lRapCurrentPart   = MyRapidity(p0->Energy(),p0->Pz());
        //lEtaCurrentPart   = p0->Eta();
@@ -2166,30 +763,27 @@ void AliAnalysisTaskPerformanceStrange::UserExec(Option_t *)
        }
        
        if (lPdgcodeCurrentPart==310)   {
-         fHistMCtracksProdRadiusK0s->Fill(mcPosX,mcPosY);
-         fHistMCtracksDecayRadiusK0s->Fill(mcDecayPosR);
-
-         if (TMath::Abs(lRapCurrentPart) < 0.3) fHistMCPtAllK0sRap3->Fill(lPtCurrentPart);
-
          if (TMath::Abs(lRapCurrentPart) < lCutRap) fHistMCPtAllK0s->Fill(lPtCurrentPart);
        }
        else if (lPdgcodeCurrentPart==3122)  {
-         fHistMCtracksProdRadiusLambda->Fill(mcPosX,mcPosY);
-         fHistMCtracksDecayRadiusLambda->Fill(mcDecayPosR);
-
-         if (TMath::Abs(lRapCurrentPart) < 0.3) fHistMCPtAllLambdaRap3->Fill(lPtCurrentPart);
-
          if (TMath::Abs(lRapCurrentPart) < lCutRap) fHistMCPtAllLambda->Fill(lPtCurrentPart);
        }
        else if (lPdgcodeCurrentPart==-3122) {
-         fHistMCtracksProdRadiusAntiLambda->Fill(mcPosX,mcPosY);
-         fHistMCtracksDecayRadiusAntiLambda->Fill(mcDecayPosR);
-
-         if (TMath::Abs(lRapCurrentPart) < 0.3) fHistMCPtAllAntiLambdaRap3->Fill(lPtCurrentPart);
-
          if (TMath::Abs(lRapCurrentPart) < lCutRap) fHistMCPtAllAntiLambda->Fill(lPtCurrentPart);
        }
-       
+       else if (lPdgcodeCurrentPart==3312) { // Xi
+         if (TMath::Abs(lRapCurrentPart) < lCutRap) fHistMCPtAllXi->Fill(lPtCurrentPart);
+       }
+       else if (lPdgcodeCurrentPart==-3312) { //Anti Xi
+         if (TMath::Abs(lRapCurrentPart) < lCutRap) fHistMCPtAllAntiXi->Fill(lPtCurrentPart);
+       }       
+       else if (lPdgcodeCurrentPart==3334) {  //Omega
+         if (TMath::Abs(lRapCurrentPart) < lCutRap) fHistMCPtAllXi->Fill(lPtCurrentPart);
+       }
+       else if (lPdgcodeCurrentPart==-3334) { //Anti Omega
+         if (TMath::Abs(lRapCurrentPart) < lCutRap) fHistMCPtAllAntiXi->Fill(lPtCurrentPart);
+       }       
+
        if ( ( ( TMath::Abs(lPdgCurrentMother) == 3212)  ||
               ( TMath::Abs(lPdgCurrentMother) == 3224)  ||
               ( TMath::Abs(lPdgCurrentMother) == 3214)  ||
@@ -2197,21 +791,20 @@ void AliAnalysisTaskPerformanceStrange::UserExec(Option_t *)
             && ( iCurrentMother <= lNbMCPrimary )
             ) lComeFromSigma = 1;
        else lComeFromSigma = 0;
-       
+      
        //*********************************************
        // Now keep only primary particles   
        // if ( ( iMc > lNbMCPrimary ) && (!lComeFromSigma) ) continue;
-
+       
        //*************************************//
         // new definition of primary particles //
         //*************************************//
-
        Double_t dx = 0;
        Double_t dy = 0;
        Double_t dz = 0;
        Double_t ProdDistance = 0;
 
-
        dx = ( (mcPrimaryVtx.At(0)) - (p0->Vx()) ); 
        dy = ( (mcPrimaryVtx.At(1)) - (p0->Vy()) );
        dz = ( (mcPrimaryVtx.At(2)) - (p0->Vz()) );
@@ -2222,99 +815,42 @@ void AliAnalysisTaskPerformanceStrange::UserExec(Option_t *)
     
        lNtimesReconstructedK0s   = 0; lNtimesReconstructedLambda   = 0; lNtimesReconstructedAntiLambda   = 0;
 
         // Rap distribution
         if (lPdgcodeCurrentPart==310) {
          fHistMCRapK0s->Fill(lRapCurrentPart);
-         if (lPtCurrentPart < 0.2 && lPtCurrentPart < 3.0)
-           fHistMCRapInPtRangeK0s->Fill(lRapCurrentPart);
        }
 
        if (lPdgcodeCurrentPart==3122) {
          fHistMCRapLambda->Fill(lRapCurrentPart);
-         if (lPtCurrentPart < 0.6 && lPtCurrentPart < 3.5)
-           fHistMCRapInPtRangeLambda->Fill(lRapCurrentPart);
        }
 
        if (lPdgcodeCurrentPart==-3122) {
          fHistMCRapAntiLambda->Fill(lRapCurrentPart);
-         if (lPtCurrentPart < 0.6 && lPtCurrentPart < 3.5)
-           fHistMCRapInPtRangeAntiLambda->Fill(lRapCurrentPart);
        }
-
-       if (lPdgcodeCurrentPart==3312 || lPdgcodeCurrentPart==-3312) {
+       /*      
+       if (lPdgcodeCurrentPart==3312) {
          fHistMCRapXi->Fill(lRapCurrentPart);
-         if (lPtCurrentPart < 0.6 && lPtCurrentPart < 3.0)
-           fHistMCRapInPtRangeXi->Fill(lRapCurrentPart);
-       }
-
-       if (lPdgcodeCurrentPart==333) {
-         fHistMCRapPhi->Fill(lRapCurrentPart);
-         if (lPtCurrentPart < 0.7 && lPtCurrentPart < 3.0)
-           fHistMCRapInPtRangePhi->Fill(lRapCurrentPart);
        }
+       */
        // Rapidity Cut
        if (TMath::Abs(lRapCurrentPart) > lCutRap) continue;
  
        if (lPdgcodeCurrentPart==310) {
-         fHistMCProdRadiusK0s->Fill(mcPosR);
-         fHistMCPrimDecayRadiusK0s->Fill(mcDecayPosR);
          fHistMCPtK0s->Fill(lPtCurrentPart);
-
-         if (TMath::Abs(lRapCurrentPart) < 0.3) fHistMCPtK0sRap3->Fill(lPtCurrentPart);
-
-
-         fHistNTimesRecK0s->Fill(lNtimesReconstructedK0s);
-         fHistNTimesRecK0sVsPt->Fill(lPtCurrentPart,lNtimesReconstructedK0s);
        }
        else 
          if (lPdgcodeCurrentPart==3122) {
-           fHistMCProdRadiusLambda->Fill(mcPosR);
-           fHistMCPrimDecayRadiusLambda->Fill(mcDecayPosR);
            fHistMCPtLambda->Fill(lPtCurrentPart);        
-
-           if (TMath::Abs(lRapCurrentPart) < 0.3) fHistMCPtLambdaRap3->Fill(lPtCurrentPart);
-
-           fHistNTimesRecLambda->Fill(lNtimesReconstructedLambda);
-           fHistNTimesRecLambdaVsPt->Fill(lPtCurrentPart,lNtimesReconstructedLambda);
-           if (lComeFromSigma) fHistMCPtLambdaFromSigma->Fill(lPtCurrentPart);
-
-           //printf("found Lambda MC pT=%e\n",lPtCurrentPart);
-           //printf("found Lambda MC Plabel=%d PPDGcode=%d Nlabel=%d NPDGcode=%d\n\n",id0,lPdgCurrentDaughter0,id1,lPdgCurrentDaughter1); 
-         
          }
          else 
            if (lPdgcodeCurrentPart==-3122) {
-             fHistMCProdRadiusAntiLambda->Fill(mcPosR);
-             fHistMCPrimDecayRadiusAntiLambda->Fill(mcDecayPosR);
-
              fHistMCPtAntiLambda->Fill(lPtCurrentPart);          
-
-             if (TMath::Abs(lRapCurrentPart) < 0.3) fHistMCPtAntiLambdaRap3->Fill(lPtCurrentPart);
-
-
-             fHistNTimesRecAntiLambda->Fill(lNtimesReconstructedAntiLambda);
-             fHistNTimesRecAntiLambdaVsPt->Fill(lPtCurrentPart,lNtimesReconstructedAntiLambda);
-             if (lComeFromSigma) fHistMCPtAntiLambdaFromSigma->Fill(lPtCurrentPart);
-
-             //printf("found Lambda MC pT=%e\n",lPtCurrentPart);
-             //printf("found Lambda MC Plabel=%d PPDGcode=%d Nlabel=%d NPDGcode=%d\n\n",id0,lPdgCurrentDaughter0,id1,lPdgCurrentDaughter1); 
-         
            }
-
-       
       } // end loop ESD MC
-      
     } // end ESD condition
-
-
   } // End Loop over MC condition
-
   
-
-
-
   //************************************
   // ESD loop 
   //************************************
@@ -2327,6 +863,10 @@ void AliAnalysisTaskPerformanceStrange::UserExec(Option_t *)
 
   //Multiplcity:
   Int_t    nv0sTot= 0, nv0s = 0;
+  Int_t    nCasTot= 0;
+
+  // Qualities
+  Double_t v0q = 0, kinCasQual = 0;
   //  Int_t nv0sMI =0;   
   // Variables:
   Double_t  lV0Position[3];
@@ -2344,22 +884,15 @@ void AliAnalysisTaskPerformanceStrange::UserExec(Option_t *)
   Double_t lcTauAntiLambda     = 0;
   Double_t lcTauK0s            = 0;
   Int_t    lOnFlyStatus        = 0;
-  //Float_t   tdcaPosToPrimVertexXYZ[2], tdcaNegToPrimVertexXYZ[2]; // ..[0] = Impact parameter in XY plane and ..[1] = Impact parameter in Z            
-  //Double_t  tdcaDaughterToPrimVertex[2];                          // ..[0] = Pos and ..[1] = Neg
-
   
 
   Double_t lInvMassK0s = 0, lInvMassLambda = 0, lInvMassAntiLambda = 0;
+  Double_t lInvMassXi = 0,lInvMassAntiXi = 0,lInvMassOmega = 0,lInvMassAntiOmega = 0;
   Double_t lPtK0s      = 0, lPtLambda      = 0, lPtAntiLambda      = 0;
   Double_t lRapK0s     = 0, lRapLambda     = 0, lRapAntiLambda     = 0;
   //  Double_t lEtaK0s     = 0, lEtaLambda     = 0, lEtaAntiLambda     = 0;
   Double_t lAlphaV0      = 0, lPtArmV0       = 0;
 
-  //  Double_t lPzK0s      = 0, lPzLambda      = 0,  lPzAntiLambda      = 0;
-
-
-  //  Double_t lV0Eta = 999;
-  
   // to study Associated V0s:
   Int_t    lIndexTrackPos       = 0, lIndexTrackNeg         = 0;
   UInt_t   lLabelTrackPos       = 0, lLabelTrackNeg         = 0;
@@ -2394,31 +927,14 @@ void AliAnalysisTaskPerformanceStrange::UserExec(Option_t *)
   // Daughters
   AliESDtrack  *myTrackPos  = NULL;
   AliESDtrack  *myTrackNeg  = NULL;
-  //  AliVParticle *lVPartPos   = NULL;
-  // AliVParticle *lVPartNeg   = NULL;
-
-  //  V0 momentum      
-  //Double_t lPosMom = 0;
-  //Double_t lNegMom = 0;
-  //  Double_t lPmom = 0;
-  //  Double_t lPLambda = 0;
 
   // Daughters' momentum:
   Double_t  lMomPos[3] = {999,999,999};
   Double_t  lMomNeg[3] = {999,999,999};
-  //  Double_t  lPtPos = 999, lPtNeg = 999;
-  //  Double_t  lPPos = 999, lPNeg = 999;
 
   // Inner Wall parameters:
   Double_t  lMomInnerWallPos =999, lMomInnerWallNeg = 999;
 
-  // AliKF Chi2 and Armenteros variables
-  //  Double_t lChi2KFK0s  = 0, lChi2KFLambda = 0,  lChi2KFAntiLambda = 0;
-  //  Double_t lAlphaV0K0s = 0, lAlphaV0Lambda = 0,  lAlphaV0AntiLambda = 0;
-  //Double_t lPtArmV0K0s = 0, lPtArmV0Lambda = 0,  lPtArmV0AntiLambda = 0;
-  //  Double_t lQlPos   = 0, lQlNeg   = 0;
-
-
   // PID
   Float_t nSigmaPosPion   = 0;
   Float_t nSigmaNegPion   = 0;
@@ -2426,19 +942,38 @@ void AliAnalysisTaskPerformanceStrange::UserExec(Option_t *)
   Float_t nSigmaPosProton = 0;
   Float_t nSigmaNegProton = 0;
   
+// Bachelor                                                                                                                                                                  
+Bool_t lIsBachelorKaonForTPC = kFALSE;
+Bool_t lIsBachelorPionForTPC = kFALSE;
+
+// Negative V0 daughter                                                                                                                                                      
+Bool_t lIsNegPionForTPC   = kFALSE;
+Bool_t lIsNegProtonForTPC = kFALSE;
 
+// Positive V0 daughter                                                                                                                                                      
+Bool_t lIsPosPionForTPC   = kFALSE;
+Bool_t lIsPosProtonForTPC = kFALSE;
+      
+ Double_t lEffMassXi      = 0. ;
+ Double_t lChi2Xi         = -1. ;
+ Double_t lDcaXiDaughters = -1. ;
+ Double_t lXiCosineOfPointingAngle = -1. ;
+ Double_t lPosXi[3] = { -1000.0, -1000.0, -1000.0 };
+ Double_t lXiRadius2D = -1000. ;
+ Double_t lXiRadius3D = -1000. ;
+
+ Int_t lIDs = 0;
+ Int_t lIndexBachelorMother = 0;
+ Int_t lIndexBachelor = 0;
   //  Int_t lCheckPIDK0sPosDaughter        = 0, lCheckPIDK0sNegDaughter        = 0;
   Int_t lCheckPIDLambdaPosDaughter     = 0;// lCheckPIDLambdaNegDaughter     = 0;
   //  Int_t lCheckPIDAntiLambdaPosDaughter = 0;
   Int_t lCheckPIDAntiLambdaNegDaughter = 0;
-
-  
   
   //****************************************************
   // Primary Vertex cuts &
   // Magnetic field and Quality tracks cuts 
   //****************************************************
-
   Double_t  lPrimaryVtxPosition[3];
   Double_t  lPrimaryVtxCov[6];
   Double_t  lPrimaryVtxChi2 = 999;
@@ -2448,8 +983,7 @@ void AliAnalysisTaskPerformanceStrange::UserExec(Option_t *)
      
   AliAODVertex *myPrimaryVertex = NULL;
   //const AliVVertex *mySPDPrimaryVertex = NULL;
-
-     
+    
   const AliMultiplicity *myMultiplicty = ((AliESDEvent*)fESD)->GetMultiplicity();
 
   if(fAnalysisType == "ESD") {  
@@ -2498,22 +1032,15 @@ void AliAnalysisTaskPerformanceStrange::UserExec(Option_t *)
       Double_t zRes = TMath::Sqrt(cov[5]);
       if (vtxTyp.Contains("vertexer:Z") && (zRes>0.25)) return;
       if (TMath::Abs(spdVtx->GetZ() - trkVtx->GetZ())>0.5) return;
-
       if (TMath::Abs(zvtx) > cutPrimVertex) return;
-
+      trkVtx->GetXYZ(lPrimaryVtxPosition);
     }
 
-
     // Number of Tracklets:
-
     fHistTrackletPerEvent->Fill(myMultiplicty->GetNumberOfTracklets());
-
     lMagneticField = ((AliESDEvent*)fESD)->GetMagneticField();
-
     fHistTPCTracks->Fill(AliESDtrackCuts::GetReferenceMultiplicity((AliESDEvent*)fESD, kTRUE));
-
     //////simple chack for multiplicity////////////////////////////////////////////////////////
-
     Int_t i =0;
 
     for (Int_t jTracks=0;jTracks<fESD->GetNumberOfTracks();jTracks++){
@@ -2523,40 +1050,159 @@ void AliAnalysisTaskPerformanceStrange::UserExec(Option_t *)
       Float_t z=0;
       tPCtrack->GetImpactParameters(xy,z);
       if ((fTracksCuts->IsSelected(tPCtrack))&&(xy<1.0)&&(z<1.0)) {i=i+1;}
-                       
     }
-
     fHistTPCMult->Fill(i);
-
     /////////////////////////////////////////////////////////////////////////////////////////
+      }
 
-
-  }
   fHistPrimaryVertexX->Fill(lPrimaryVtxPosition[0]);
   fHistPrimaryVertexY->Fill(lPrimaryVtxPosition[1]);
   fHistPrimaryVertexZ->Fill(lPrimaryVtxPosition[2]);
   //Double_t lrcPrimVtxR = TMath::Sqrt(lPrimaryVtxPosition[0]*lPrimaryVtxPosition[0]+lPrimaryVtxPosition[0]*lPrimaryVtxPosition[0]);
+  ////*************************
+  //// Cascade loop ***********
+  ////*************************
+  nCasTot = fESD->GetNumberOfCascades();
+
+  for (Int_t iCas = 0; iCas < nCasTot; iCas++) {
+
+    AliESDcascade *Cas = ((AliESDEvent*)fESD)->GetCascade(iCas);
+      if (!Cas) continue;
+
+      lIndexBachelor = Cas->GetBindex();
+      fHistMCIndexes->Fill(lIndexBachelor);
+
 
-  fHistPrimaryVertexResX->Fill(lResPrimaryVtxX);
-  fHistPrimaryVertexResY->Fill(lResPrimaryVtxY);
-  fHistPrimaryVertexResZ->Fill(lResPrimaryVtxZ);
+      UInt_t lIdxPosXi        = (UInt_t) TMath::Abs( Cas->GetPindex() );
+      UInt_t lIdxNegXi        = (UInt_t) TMath::Abs( Cas->GetNindex() );
+      UInt_t lBachIdx         = (UInt_t) TMath::Abs( lIndexBachelor );
+      // Care track label can be negative in MC production (linked with the track quality)                                                                                 
+      // However = normally, not the case for track index ...                                                                                                              
+
+      // FIXME : rejection of a double use of a daughter track (nothing but just a crosscheck of what is done in the cascade vertexer)                                     
+      if(lBachIdx == lIdxNegXi) {
+       AliWarning("Pb / Idx(Bach. track) = Idx(Neg. track) ... continue!"); continue;
+      }
+      if(lBachIdx == lIdxPosXi) {
+       AliWarning("Pb / Idx(Bach. track) = Idx(Pos. track) ... continue!"); continue;
+      }
+
+      lEffMassXi                      = Cas->GetEffMassXi();
+      lChi2Xi                         = Cas->GetChi2Xi();
+      lDcaXiDaughters                 = Cas->GetDcaXiDaughters();
+      lXiCosineOfPointingAngle        = Cas->GetCascadeCosineOfPointingAngle( lPrimaryVtxPosition[0],
+                                                                            lPrimaryVtxPosition[1],
+                                                                            lPrimaryVtxPosition[2] );
+      // Take care : the best available vertex should be used (like in AliCascadeVertexer)                                                                                 
+
+      Cas->GetXYZcascade( lPosXi[0],  lPosXi[1], lPosXi[2] );
+      lXiRadius2D    = TMath::Sqrt( lPosXi[0]*lPosXi[0]  +  lPosXi[1]*lPosXi[1] );
+      lXiRadius3D    = TMath::Sqrt( lPosXi[0]*lPosXi[0]  +  lPosXi[1]*lPosXi[1] +  lPosXi[2]*lPosXi[2]);
 
-  //***********************
-  // AliKF Primary Vertex
 
-  AliKFVertex primaryVtxKF( *myPrimaryVertex );
-  AliKFParticle::SetField(lMagneticField);
+      AliESDtrack *pTrackXi           = fESD->GetTrack( lIdxPosXi );
+      AliESDtrack *nTrackXi           = fESD->GetTrack( lIdxNegXi );
+      AliESDtrack *bachTrackXi        = fESD->GetTrack( lBachIdx );
+      if (!pTrackXi || !nTrackXi || !bachTrackXi ) {
+       AliWarning("ERROR: Could not retrieve one of the 3 ESD daughter tracks of the cascade ...");
+       continue;
+      }
+            
+      // Bachelor                                                                                                                                                                  
+      if (TMath::Abs(fPIDResponse->NumberOfSigmasTPC( bachTrackXi,AliPID::kKaon)) < 4) lIsBachelorKaonForTPC = kTRUE;
+      if (TMath::Abs(fPIDResponse->NumberOfSigmasTPC( bachTrackXi,AliPID::kPion)) < 4) lIsBachelorPionForTPC = kTRUE;
+
+      // Negative V0 daughter                                                                                                                                                      
+      if (TMath::Abs(fPIDResponse->NumberOfSigmasTPC( nTrackXi,AliPID::kPion   )) < 4) lIsNegPionForTPC   = kTRUE;
+      if (TMath::Abs(fPIDResponse->NumberOfSigmasTPC( nTrackXi,AliPID::kProton )) < 4) lIsNegProtonForTPC = kTRUE;
+
+      // Positive V0 daughter                                                                                                                                                      
+      if (TMath::Abs(fPIDResponse->NumberOfSigmasTPC( pTrackXi,AliPID::kPion   )) < 4) lIsPosPionForTPC   = kTRUE;
+      if (TMath::Abs(fPIDResponse->NumberOfSigmasTPC( pTrackXi,AliPID::kProton )) < 4) lIsPosProtonForTPC = kTRUE;
+      
 
+      if( bachTrackXi->Charge() < 0 ) {
+       v0q = 0;
+      kinCasQual =  Cas->ChangeMassHypothesis(v0q,3312);
+      lInvMassXi = Cas->GetEffMassXi();
+       v0q = 0;
+      kinCasQual =  Cas->ChangeMassHypothesis(v0q,3334);
+      lInvMassOmega = Cas->GetEffMassXi();
+      }
 
+      if( bachTrackXi->Charge() > 0 ) {
+       v0q = 0;
+      kinCasQual =  Cas->ChangeMassHypothesis(v0q,-3312);
+      lInvMassAntiXi = Cas->GetEffMassXi();
+       v0q = 0;
+      kinCasQual =  Cas->ChangeMassHypothesis(v0q,-3334);
+      lInvMassAntiOmega = Cas->GetEffMassXi();
+      }
 
 
-  //***Rerun the V0 finder
+      if( lIsBachelorPionForTPC == kTRUE && lIsNegPionForTPC == kTRUE && lIsPosProtonForTPC == kTRUE){
+      fHistMassXi->Fill(lInvMassXi);
+      if(fAnalysisMC && (lInvMassXi > 1.25) && (lInvMassXi < 1.4)){
+       TParticle *Bachelor = stack->Particle(lIndexBachelor);        
+               //        lIndexBachelorMother = Bachelor->GetMother(1);
+       lIndexBachelorMother = Bachelor->GetFirstMother();
+       if(lIndexBachelorMother == -1) continue;
+        if(lIndexBachelorMother == 1){  fHistMassXiVsPID->Fill(lInvMassXi,0);fHistCompositionXi->Fill(0); continue;} //combinatorial
+       lIDs = stack->Particle(lIndexBachelorMother)->GetPdgCode();
+       if(lIDs == 111) { fHistMassXiVsPID->Fill(lInvMassXi,1); } //pi0
+       if(lIDs == 213) { fHistMassXiVsPID->Fill(lInvMassXi,2); } //rho+
+       if(lIDs == 2212) { fHistMassXiVsPID->Fill(lInvMassXi,3); } //Proton
+       if(lIDs == 2112) { fHistMassXiVsPID->Fill(lInvMassXi,4); } //Neutron
+        fHistCompositionXi->Fill(lIDs);      
+       }
+
+      }
+      if( lIsBachelorKaonForTPC == kTRUE && lIsNegPionForTPC == kTRUE && lIsPosProtonForTPC == kTRUE){
+      fHistMassOmega->Fill(lInvMassOmega);
+      if(fAnalysisMC && lInvMassOmega > 1.6 && lInvMassOmega < 1.75){
+       TParticle *Bachelor = stack->Particle(lIndexBachelor);        
+               //        lIndexBachelorMother = Bachelor->GetMother(1);
+       lIndexBachelorMother = Bachelor->GetFirstMother();
+       if(lIndexBachelorMother == -1) continue;
+       if(lIndexBachelorMother == 1) { fHistCompositionOmega->Fill(0); continue;}
+       lIDs = stack->Particle(lIndexBachelorMother)->GetPdgCode();
+      fHistCompositionOmega->Fill(lIDs);      
+       }
+
+      }
+      if( lIsBachelorPionForTPC == kTRUE && lIsPosPionForTPC == kTRUE && lIsNegProtonForTPC == kTRUE){
+      fHistMassAntiXi->Fill(lInvMassAntiXi);
+      if(fAnalysisMC && lInvMassAntiXi > 1.25 && lInvMassAntiXi < 1.4){
+       TParticle *Bachelor = stack->Particle(lIndexBachelor);        
+               //        lIndexBachelorMother = Bachelor->GetMother(1);
+       lIndexBachelorMother = Bachelor->GetFirstMother();
+       if(lIndexBachelorMother == -1) continue;
+       if(lIndexBachelorMother == 1) { fHistCompositionAntiXi->Fill(0); continue;}
+       lIDs = stack->Particle(lIndexBachelorMother)->GetPdgCode();
+      fHistCompositionAntiXi->Fill(lIDs);      
+       }
+
+      }
+      if( lIsBachelorKaonForTPC == kTRUE && lIsPosPionForTPC == kTRUE && lIsNegProtonForTPC == kTRUE){
+      fHistMassAntiOmega->Fill(lInvMassAntiOmega);
+      if(fAnalysisMC && lInvMassAntiOmega > 1.6 && lInvMassAntiOmega < 1.75){
+       TParticle *Bachelor = stack->Particle(lIndexBachelor);        
+               //        lIndexBachelorMother = Bachelor->GetMother(1);
+       lIndexBachelorMother = Bachelor->GetFirstMother();
+       if(lIndexBachelorMother == -1) continue;
+       if(lIndexBachelorMother == 1) { fHistCompositionAntiOmega->Fill(0); continue;}
+       lIDs = stack->Particle(lIndexBachelorMother)->GetPdgCode();
+      fHistCompositionAntiOmega->Fill(lIDs);      
+       }
+
+      }
+           
+  
+    
+  
+  } /// end of Cascade loop
 
-  //  fESD->ResetV0s();
-  //  AliV0vertexer v0Vertexer;
-  //  v0Vertexer.SetCuts(fCuts);
-  //  v0Vertexer.Tracks2V0vertices(fESD);
   
   ////*************************
   //// V0 loop ****************
@@ -2566,27 +1212,19 @@ void AliAnalysisTaskPerformanceStrange::UserExec(Option_t *)
   if (!nv0sTot) fHistNumberEvents->Fill(6.5);
 
   for (Int_t iV0 = 0; iV0 < nv0sTot; iV0++) {
-    
+     
     lIndexPosMother     = 0; lIndexNegMother     = 0; lIndexMotherOfMother       = 0;
     lCheckPIdK0Short    = 0; lCheckMcK0Short     = 0; lCheckSecondaryK0s         = 0;
     lCheckPIdLambda     = 0; lCheckMcLambda      = 0; lCheckSecondaryLambda      = 0;
     lCheckPIdAntiLambda = 0; lCheckMcAntiLambda  = 0; lCheckSecondaryAntiLambda  = 0;       
     lComeFromSigma      = 0; //lCheckGamma = 0;
     
-    
     if(fAnalysisType == "ESD") {
 
-
       AliESDv0 *v0 = ((AliESDEvent*)fESD)->GetV0(iV0);
       if (!v0) continue;
-
       //      if ((v0->Pt())<0.6) continue;
-      
-      // Primary vertex:
-      fHistPrimaryVertexPosXV0events->Fill(lPrimaryVtxPosition[0]);
-      fHistPrimaryVertexPosYV0events->Fill(lPrimaryVtxPosition[1]);
-      fHistPrimaryVertexPosZV0events->Fill(lPrimaryVtxPosition[2]);
-      
+     
       // V0's Daughters
       lIndexTrackPos = TMath::Abs(v0->GetPindex());
       lIndexTrackNeg = TMath::Abs(v0->GetNindex());
@@ -2617,24 +1255,17 @@ void AliAnalysisTaskPerformanceStrange::UserExec(Option_t *)
       lV0tDecayLength = TMath::Sqrt(TMath::Power(lV0Position[0] - lPrimaryVtxPosition[0],2) +
                                    TMath::Power(lV0Position[1] - lPrimaryVtxPosition[1],2));
 
-
-
       if( myTrackPosTest->GetSign() ==1){
-       
        myTrackPos = ((AliESDEvent*)fESD)->GetTrack(lIndexTrackPos);
        myTrackNeg = ((AliESDEvent*)fESD)->GetTrack(lIndexTrackNeg);
-
        // Daughters' momentum;
        v0->GetPPxPyPz(lMomPos[0],lMomPos[1],lMomPos[2]);
        v0->GetNPxPyPz(lMomNeg[0],lMomNeg[1],lMomNeg[2]);
       }
            
       if( myTrackPosTest->GetSign() ==-1){
-       
        myTrackPos = ((AliESDEvent*)fESD)->GetTrack(lIndexTrackNeg);
        myTrackNeg = ((AliESDEvent*)fESD)->GetTrack(lIndexTrackPos);
-
-
        // Daughters' momentum;
        v0->GetPPxPyPz(lMomNeg[0],lMomNeg[1],lMomNeg[2]);
        v0->GetNPxPyPz(lMomPos[0],lMomPos[1],lMomPos[2]);
@@ -2643,22 +1274,6 @@ void AliAnalysisTaskPerformanceStrange::UserExec(Option_t *)
       lLabelTrackPos = (UInt_t)TMath::Abs(myTrackPos->GetLabel());
       lLabelTrackNeg = (UInt_t)TMath::Abs(myTrackNeg->GetLabel());
 
-      // Daughters Pt and P:
-      //      lPtPos = TMath::Sqrt(lMomPos[0]*lMomPos[0] + lMomPos[1]*lMomPos[1]);
-      //      lPtNeg = TMath::Sqrt(lMomNeg[0]*lMomNeg[0] + lMomNeg[1]*lMomNeg[1]);
-
-      //      lPPos = TMath::Sqrt(lMomPos[0]*lMomPos[0] + lMomPos[1]*lMomPos[1] + lMomPos[2]*lMomPos[2]);
-      //      lPNeg = TMath::Sqrt(lMomNeg[0]*lMomNeg[0] + lMomNeg[1]*lMomNeg[1] + lMomNeg[2]*lMomNeg[2]);
-
-
-      //  lPosMom = myTrackPos->GetInnerParam()->GetP();
-
-      //  lNegMom = myTrackNeg->GetInnerParam()->GetP();
-
-      //        lPmom  = TMath::Sqrt(lMomPos[0]*lMomPos[0]+lMomPos[1]*lMomPos[1]+lMomPos[2]+lMomPos[2]);
-
-
-
       // Inner Wall parameter:
       const AliExternalTrackParam *myInnerWallTrackPos = myTrackPos->GetInnerParam(); 
 
@@ -2669,7 +1284,6 @@ void AliAnalysisTaskPerformanceStrange::UserExec(Option_t *)
 
       if (myInnerWallTrackNeg) {lMomInnerWallNeg = myInnerWallTrackNeg->GetP();}  
       else {continue;}
-             
 
       // DCA between daughter and Primary Vertex:
       if (myTrackPos) lDcaPosToPrimVertex = TMath::Abs(myTrackPos->GetD(lPrimaryVtxPosition[0],lPrimaryVtxPosition[1],lMagneticField) );
@@ -2700,12 +1314,10 @@ void AliAnalysisTaskPerformanceStrange::UserExec(Option_t *)
       lPtLambda = v0->Pt();
       //lPzLambda = v0->Pz();
 
-
       v0->ChangeMassHypothesis(-3122);
       lInvMassAntiLambda = v0->GetEffMass();
       lPtAntiLambda = v0->Pt();
       //lPzAntiLambda = v0->Pz();
-  
 
       // Rapidity:
       lRapK0s    = v0->Y(310);
@@ -2735,8 +1347,6 @@ void AliAnalysisTaskPerformanceStrange::UserExec(Option_t *)
        nSigmaPosPion = 0; nSigmaNegPion =0; nSigmaPosProton = 0; nSigmaNegProton= 0;
       }
       
-      
-      
       // Monte-Carlo particle associated to reconstructed particles: 
       if (fAnalysisMC) {
        //if (lLabelTrackPos < 0 || lLabelTrackNeg < 0) continue;
@@ -2752,18 +1362,7 @@ void AliAnalysisTaskPerformanceStrange::UserExec(Option_t *)
        lIndexPosMother = lMCESDPartPos->GetFirstMother();
        lIndexNegMother = lMCESDPartNeg->GetFirstMother();
 
-       //////////////////////////////
-       //      if (lIndexPosMother == -1)      {
-       //      if (negPiKF) delete negPiKF; negPiKF=NULL;
-       //      if (posPiKF) delete posPiKF; posPiKF=NULL;
-       //      if (posPKF)  delete posPKF;  posPKF=NULL;
-       //      if (negAPKF) delete negAPKF; negAPKF=NULL;
-
-       //      continue;}
-
        if (lIndexPosMother == -1) {
-
-
          lPdgcodeMother = 0;
          lIndexMotherOfMother = 0;
          mcPosX = 0;
@@ -2779,8 +1378,6 @@ void AliAnalysisTaskPerformanceStrange::UserExec(Option_t *)
 
        else {
 
-
-
          TParticle  *lMCESDMother    = stack->Particle(lIndexPosMother);
          if (!lMCESDMother)    { continue;}
          lPdgcodeMother         = lMCESDMother->GetPdgCode();
@@ -2805,11 +1402,6 @@ void AliAnalysisTaskPerformanceStrange::UserExec(Option_t *)
        }
       }
     } // end ESD condition
-
-
-    
-
-    
     
     // Multiplicity:
     if(!lOnFlyStatus) nv0s++;
@@ -2823,13 +1415,13 @@ void AliAnalysisTaskPerformanceStrange::UserExec(Option_t *)
     // Look for associated particles:
     if (fAnalysisMC) {
       if( (lIndexPosMother==-1) || (lIndexNegMother==-1) ) {
-       fHistMCDaughterTrack->Fill(1);
+       //      fHistMCDaughterTrack->Fill(1);
       }
       
       else if( ( (lPDGCodePosDaughter==+211) && (lPDGCodeNegDaughter==-211) ) ) 
        {
          lCheckPIdK0Short    = 1;
-         fHistMCDaughterTrack->Fill(3);
+         //      fHistMCDaughterTrack->Fill(3);
          if ( (lIndexPosMother==lIndexNegMother) &&
               (lPdgcodeMother==310) ) {
 
@@ -2856,7 +1448,7 @@ void AliAnalysisTaskPerformanceStrange::UserExec(Option_t *)
       else if( ( (lPDGCodePosDaughter==+2212) && (lPDGCodeNegDaughter==-211)  )  ) 
        {
          lCheckPIdLambda     = 1;
-         fHistMCDaughterTrack->Fill(5);
+         //      fHistMCDaughterTrack->Fill(5);
          if ( (lIndexPosMother==lIndexNegMother) &&
               (lPdgcodeMother==3122)  ){
            if ( ( TMath::Abs(lPdgcodeMotherOfMother) == 3212) ||
@@ -2885,14 +1477,14 @@ void AliAnalysisTaskPerformanceStrange::UserExec(Option_t *)
 
            if (ProdDistance < 0.001) lCheckMcLambda  = 1;
            else lCheckSecondaryLambda = 1;
-
+           
        
          }
        }
       else if( ( (lPDGCodePosDaughter==211)   && (lPDGCodeNegDaughter==-2212) )        ) 
        {
          lCheckPIdAntiLambda = 1;
-         fHistMCDaughterTrack->Fill(7);
+         //      fHistMCDaughterTrack->Fill(7);
          if ( (lIndexPosMother==lIndexNegMother) &&
               (lPdgcodeMother==-3122) ) {
            if ( ( TMath::Abs(lPdgcodeMotherOfMother) == 3212) ||
@@ -2981,322 +1573,66 @@ void AliAnalysisTaskPerformanceStrange::UserExec(Option_t *)
 
       {continue;}
        
-    /*  {
-       if (negPiKF) delete negPiKF; negPiKF=NULL;
-       if (posPiKF) delete posPiKF; posPiKF=NULL;
-       if (posPKF)  delete posPKF;  posPKF=NULL;
-       if (negAPKF) delete negAPKF; negAPKF=NULL;
-
-       continue;}
-       //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
-       */
-
 
     /////////////////////////////////
     //PID for Lambda and AntiLambda
     /////////////////////////////////
-    
     if(fUsePID.Contains("withPID") && (lCheckPIDAntiLambdaNegDaughter==0) && (lCheckPIDLambdaPosDaughter==1)) LambdaPID = 1;
-
     else LambdaPID =0;
-
     if(fUsePID.Contains("withPID") && (lCheckPIDLambdaPosDaughter==0) && (lCheckPIDAntiLambdaNegDaughter==1)) AntiLambdaPID = 1;
-
     else AntiLambdaPID =0;
 
-
-    ///////////////////////////////////////////////////////////////////////////////////////////
-    //    if ((LambdaPID==1 && lPosMom <=1) || (lPosMom>1) ||  !(fUsePID.Contains("withPID"))){  
-    //      if ((TMath::Abs(lRapLambda) < lCutRap) && lOnFlyStatus==0) {
-    //      fHistTPCsigPLambda->Fill(lPosMom,myTrackPos->GetTPCsignal());
-    //
-    //            TestTrackCuts->SetRequireITSRefit(kFALSE);
-    //    if(TestTrackCuts->IsSelected(myTrackPos)) fHistPtTracksP->Fill(V0mom[2]);
-    //    TestTrackCuts->SetRequireITSRefit(kTRUE);
-    //    if(TestTrackCuts->IsSelected(myTrackPos)) fHistPtTracksPITSRefit->Fill(V0mom[2]);
-    //        }
-    //
-    //    }
-    //
-    //    if ((AntiLambdaPID==1 && lNegMom <=1) || (lNegMom>1) ||  !(fUsePID.Contains("withPID"))){  
-    //      if ((TMath::Abs(lRapAntiLambda) < lCutRap) && lOnFlyStatus==0) {
-    //      fHistTPCsigPAntiLambda->Fill(lNegMom,myTrackNeg->GetTPCsignal());
-    //
-    //      }
-    //    }
-    ///////////////////////////////////////////////////////////////////////////////////////////////
-
-
-    //    lPLambda = TMath::Sqrt(lPzLambda*lPzLambda + lPtLambda*lPtLambda);
-    //    lPtLambda = TMath::Sqrt(lPzLambda*lPzLambda + lPtLambda*lPtLambda);
-
-    //    lPAntiLambda = TMath::Sqrt(lPzAntiLambda*lPzAntiLambda + lPtAntiLambda*lPtAntiLambda);
-    //    lPtAntiLambda = TMath::Sqrt(lPtAntiLambda*lPtAntiLambda);
-
-    //    lPK0s = TMath::Sqrt(lPzK0s*lPzK0s + lPtK0s*lPtK0s);
-    //   lPtK0s = TMath::Sqrt(lPtK0s*lPtK0s);
-
-    //    lcTau     = (lV0DecayLength*lLambdaMass)/lPLambda;
-
     //ctau for lambda
     lcTauLambda     = (lV0tDecayLength*lInvMassLambda)/lPtLambda;
-
     //ctau for antilambda
     lcTauAntiLambda = (lV0tDecayLength*lInvMassAntiLambda)/lPtAntiLambda;
-
     //ctau for K0s
     lcTauK0s        = (lV0tDecayLength*lInvMassK0s)/lPtK0s;
 
-
     //*****************************
     // filling histograms
     //*****************************
-
     if (lPLambda <1 && lOnFlyStatus==0 ){
       fHistArmenterosPodolanski->Fill(lAlphaV0,lPtArmV0);
     }
 
-
-
     ////////////////
     //K0s particle
     ////////////////
-
-
-
     if (TMath::Abs(lRapK0s) < lCutRap ) {
-      if (lOnFlyStatus==0){
-         fHistTranscTauVsMassK0s->Fill(lcTauK0s,lInvMassK0s);
-             if (TMath::Abs(lRapK0s) < 0.3)
-              fHistTranscTauVsMassK0sRap3->Fill(lcTauK0s,lInvMassK0s);
-        if (lPtK0s >= 3 && lPtK0s <= 6 )
-          {
-            fHistTranscTauVsMassK0sptLow->Fill(lcTauK0s,lInvMassK0s);
-             if (TMath::Abs(lRapK0s) < 0.3)
-              fHistTranscTauVsMassK0sptLowRap3->Fill(lcTauK0s,lInvMassK0s);
-          }
-       }   
-
       if (lcTauK0s< cutcTauK0s) {
 
-       //////2D histos: cut vs on fly status/////////////////////
-
-       fHistDcaPosToPrimVertexK0->Fill(lDcaPosToPrimVertex,lOnFlyStatus);
-       fHistDcaNegToPrimVertexK0->Fill(lDcaNegToPrimVertex,lOnFlyStatus);
-       fHistRadiusV0K0->Fill(lV0Radius,lOnFlyStatus);
-       fHistDecayLengthV0K0->Fill(lV0DecayLength,lOnFlyStatus);
-       fHistDcaV0DaughtersK0->Fill(lDcaV0Daughters,lOnFlyStatus);
-       fHistChi2K0->Fill(lChi2V0,lOnFlyStatus);
-       fHistCosPointAngleK0->Fill(lV0cosPointAngle,lOnFlyStatus);
-
-       //////2D histos: cut vs mass///////////////////// 
-
        if (lOnFlyStatus==0){
-        
-
-         fHistMassK0->Fill(lInvMassK0s);
-         fHistMassVsRadiusK0->Fill(rcPosRK0s,lInvMassK0s);
+          fHistMassK0->Fill(lInvMassK0s);
          fHistPtVsMassK0->Fill(lInvMassK0s,lPtK0s);
-
-
-         if (TMath::Abs(lRapK0s) < 0.3 )       fHistPtVsMassK0Rap3->Fill(lInvMassK0s,lPtK0s);
-
-         fHistDcaPosToPrimVertexK0vsMassK0->Fill(lDcaPosToPrimVertex,lInvMassK0s);
-         fHistDcaNegToPrimVertexK0vsMassK0->Fill(lDcaNegToPrimVertex,lInvMassK0s);
-         fHistRadiusV0K0vsMassK0->Fill(lV0Radius,lInvMassK0s);
-         fHistDecayLengthV0K0vsMassK0->Fill(lV0DecayLength,lInvMassK0s);
-         fHistDcaV0DaughtersK0vsMassK0->Fill(lDcaV0Daughters,lInvMassK0s);
-         fHistCosPointAngleK0vsMassK0->Fill(lV0cosPointAngle,lInvMassK0s);
-         if(fQASelector){
-           if (lPtK0s>0 && lPtK0s <3){ 
-             fHistDcaPosToPrimVertexK0vsMassK0pt1->Fill(lDcaPosToPrimVertex,lInvMassK0s);
-             fHistDcaNegToPrimVertexK0vsMassK0pt1->Fill(lDcaNegToPrimVertex,lInvMassK0s);
-             fHistRadiusV0K0vsMassK0pt1->Fill(lV0Radius,lInvMassK0s);
-             fHistDecayLengthV0K0vsMassK0pt1->Fill(lV0DecayLength,lInvMassK0s);
-             fHistDcaV0DaughtersK0vsMassK0pt1->Fill(lDcaV0Daughters,lInvMassK0s);
-             fHistCosPointAngleK0vsMassK0pt1->Fill(lV0cosPointAngle,lInvMassK0s);
-           }   
-           if (lPtK0s > 3 && lPtK0s < 6){ 
-             fHistDcaPosToPrimVertexK0vsMassK0pt2->Fill(lDcaPosToPrimVertex,lInvMassK0s);
-             fHistDcaNegToPrimVertexK0vsMassK0pt2->Fill(lDcaNegToPrimVertex,lInvMassK0s);
-             fHistRadiusV0K0vsMassK0pt2->Fill(lV0Radius,lInvMassK0s);
-             fHistDecayLengthV0K0vsMassK0pt2->Fill(lV0DecayLength,lInvMassK0s);
-             fHistDcaV0DaughtersK0vsMassK0pt2->Fill(lDcaV0Daughters,lInvMassK0s);
-             fHistCosPointAngleK0vsMassK0pt2->Fill(lV0cosPointAngle,lInvMassK0s);
-           }   
-           if (lPtK0s > 6 && lPtK0s < 10){ 
-             fHistDcaPosToPrimVertexK0vsMassK0pt3->Fill(lDcaPosToPrimVertex,lInvMassK0s);
-             fHistDcaNegToPrimVertexK0vsMassK0pt3->Fill(lDcaNegToPrimVertex,lInvMassK0s);
-             fHistRadiusV0K0vsMassK0pt3->Fill(lV0Radius,lInvMassK0s);
-             fHistDecayLengthV0K0vsMassK0pt3->Fill(lV0DecayLength,lInvMassK0s);
-             fHistDcaV0DaughtersK0vsMassK0pt3->Fill(lDcaV0Daughters,lInvMassK0s);
-             fHistCosPointAngleK0vsMassK0pt3->Fill(lV0cosPointAngle,lInvMassK0s);
-           }   
-         } //fQASelector
-       }
-      } // if rap. condition
-    } //end ctau cut
-
+            } //end !lOnFlystatus
+          } //end ctau cut
+    } //end Rap condition
 
     ///////////////////
     //Lambda particle    
     ///////////////////
-
-
-
-
     if ((LambdaPID==1 && lMomInnerWallPos <=1 ) || (lMomInnerWallPos > 1) ||  !(fUsePID.Contains("withPID")  )){  
-
       if (TMath::Abs(lRapLambda) < lCutRap) {
-       if (lOnFlyStatus==0) {
-         fHistTranscTauVsMassL->Fill(lcTauLambda,lInvMassLambda);
-             if (TMath::Abs(lRapLambda) < 0.3)
-              fHistTranscTauVsMassLRap3->Fill(lcTauLambda,lInvMassLambda);
-        if (lPtLambda >= 3 && lPtLambda <= 6 )
-          {
-            fHistTranscTauVsMassLptLow->Fill(lcTauLambda,lInvMassLambda);
-             if (TMath::Abs(lRapLambda) < 0.3)
-              fHistTranscTauVsMassLptLowRap3->Fill(lcTauLambda,lInvMassLambda);
-          }
-       }   
-
        if (lcTauLambda < cutcTauL){
-         //////2D histos: cut vs on fly status/////////////////////
-
-         fHistDcaPosToPrimVertexL->Fill(lDcaPosToPrimVertex,lOnFlyStatus);
-         fHistDcaNegToPrimVertexL->Fill(lDcaNegToPrimVertex,lOnFlyStatus);
-         fHistRadiusV0L->Fill(lV0Radius,lOnFlyStatus);
-         fHistDecayLengthV0L->Fill(lV0DecayLength,lOnFlyStatus);
-         fHistDcaV0DaughtersL->Fill(lDcaV0Daughters,lOnFlyStatus);
-         fHistChi2L->Fill(lChi2V0,lOnFlyStatus);
-         fHistCosPointAngleL->Fill(lV0cosPointAngle,lOnFlyStatus);
-
-         //////2D histos: cut vs mass/////////////////////
-
          if (lOnFlyStatus==0){
            fHistMassLambda->Fill(lInvMassLambda);
-           fHistMassVsRadiusLambda->Fill(rcPosRLambda,lInvMassLambda);
            fHistPtVsMassLambda->Fill(lInvMassLambda,lPtLambda);
-
-           if (TMath::Abs(lRapLambda) < 0.3 ) fHistPtVsMassLambdaRap3->Fill(lInvMassLambda,lPtLambda);
-
            if(lPtLambda <=1) fHistNSigmaProton->Fill(nSigmaPosProton);
-
-           fHistDcaPosToPrimVertexLvsMassL->Fill(lDcaPosToPrimVertex,lInvMassLambda);
-           fHistDcaNegToPrimVertexLvsMassL->Fill(lDcaNegToPrimVertex,lInvMassLambda);
-           fHistRadiusV0LvsMassL->Fill(lV0Radius,lInvMassLambda);
-           fHistDecayLengthV0LvsMassL->Fill(lV0DecayLength,lInvMassLambda);
-           fHistDcaV0DaughtersLvsMassL->Fill(lDcaV0Daughters,lInvMassLambda);
-           fHistCosPointAngleLvsMassL->Fill(lV0cosPointAngle,lInvMassLambda);
-
-           if(fQASelector){            
-             if (lPtLambda>0 && lPtLambda <3){ 
-               fHistDcaPosToPrimVertexLambdaVsMasspt1->Fill(lDcaPosToPrimVertex,lInvMassLambda);
-               fHistDcaNegToPrimVertexLambdaVsMasspt1->Fill(lDcaNegToPrimVertex,lInvMassLambda);
-               fHistRadiusV0LambdaVsMasspt1->Fill(lV0Radius,lInvMassLambda);
-               fHistDecayLengthV0LambdaVsMasspt1->Fill(lV0DecayLength,lInvMassLambda);
-               fHistDcaV0DaughtersLambdaVsMasspt1->Fill(lDcaV0Daughters,lInvMassLambda);
-               fHistCosPointAngleLambdaVsMasspt1->Fill(lV0cosPointAngle,lInvMassLambda);
-             }   
-             if (lPtLambda > 3 && lPtLambda < 6){ 
-               fHistDcaPosToPrimVertexLambdaVsMasspt2->Fill(lDcaPosToPrimVertex,lInvMassLambda);
-               fHistDcaNegToPrimVertexLambdaVsMasspt2->Fill(lDcaNegToPrimVertex,lInvMassLambda);
-               fHistRadiusV0LambdaVsMasspt2->Fill(lV0Radius,lInvMassLambda);
-               fHistDecayLengthV0LambdaVsMasspt2->Fill(lV0DecayLength,lInvMassLambda);
-               fHistDcaV0DaughtersLambdaVsMasspt2->Fill(lDcaV0Daughters,lInvMassLambda);
-               fHistCosPointAngleLambdaVsMasspt2->Fill(lV0cosPointAngle,lInvMassLambda);
-             }   
-             if (lPtLambda > 6 && lPtLambda < 10){ 
-               fHistDcaPosToPrimVertexLambdaVsMasspt3->Fill(lDcaPosToPrimVertex,lInvMassLambda);
-               fHistDcaNegToPrimVertexLambdaVsMasspt3->Fill(lDcaNegToPrimVertex,lInvMassLambda);
-               fHistRadiusV0LambdaVsMasspt3->Fill(lV0Radius,lInvMassLambda);
-               fHistDecayLengthV0LambdaVsMasspt3->Fill(lV0DecayLength,lInvMassLambda);
-               fHistDcaV0DaughtersLambdaVsMasspt3->Fill(lDcaV0Daughters,lInvMassLambda);
-               fHistCosPointAngleLambdaVsMasspt3->Fill(lV0cosPointAngle,lInvMassLambda);
-             }   
-           }//QA Selector
          }
        }// end ctau condition
       } //end of Rap condition
     }// end of PID condition
 
-
     //////////////////////////////
     // Anti Lambda ///////////////
     //////////////////////////////
-
-
-
     if ((AntiLambdaPID==1 && lMomInnerWallNeg <=1) || (lMomInnerWallNeg>1) ||  !(fUsePID.Contains("withPID"))){  
-
       if (TMath::Abs(lRapAntiLambda) < lCutRap) {
-       if (lOnFlyStatus==0) {
-         fHistTranscTauVsMassAntiL->Fill(lcTauAntiLambda,lInvMassAntiLambda);
-             if (TMath::Abs(lRapAntiLambda) < 0.3)
-              fHistTranscTauVsMassAntiLRap3->Fill(lcTauAntiLambda,lInvMassAntiLambda);
-        if (lPtAntiLambda >= 3 && lPtAntiLambda <= 6 )
-          {
-            fHistTranscTauVsMassAntiLptLow->Fill(lcTauAntiLambda,lInvMassAntiLambda);
-             if (TMath::Abs(lRapAntiLambda) < 0.3)
-              fHistTranscTauVsMassAntiLptLowRap3->Fill(lcTauAntiLambda,lInvMassAntiLambda);
-          }
-       }   
-
        if (lcTauAntiLambda < cutcTauL){
-         //////2D histos: cut vs on fly status/////////////////////
-
-         fHistDcaPosToPrimVertexAntiL->Fill(lDcaPosToPrimVertex,lOnFlyStatus);
-         fHistDcaNegToPrimVertexAntiL->Fill(lDcaNegToPrimVertex,lOnFlyStatus);
-         fHistRadiusV0AntiL->Fill(lV0Radius,lOnFlyStatus);
-         fHistDecayLengthV0AntiL->Fill(lV0DecayLength,lOnFlyStatus);
-         fHistDcaV0DaughtersAntiL->Fill(lDcaV0Daughters,lOnFlyStatus);
-         fHistChi2AntiL->Fill(lChi2V0,lOnFlyStatus);
-         fHistCosPointAngleAntiL->Fill(lV0cosPointAngle,lOnFlyStatus);
-
-         //////2D histos: cut vs mass/////////////////////
-
          if (lOnFlyStatus==0){
-
            fHistMassAntiLambda->Fill(lInvMassAntiLambda);
-           fHistMassVsRadiusAntiLambda->Fill(rcPosRAntiLambda,lInvMassAntiLambda);
            fHistPtVsMassAntiLambda->Fill(lInvMassAntiLambda,lPtAntiLambda);
-
-           if (TMath::Abs(lRapAntiLambda) < 0.3) fHistPtVsMassAntiLambdaRap3->Fill(lInvMassAntiLambda,lPtAntiLambda);
-
-
-           fHistDcaPosToPrimVertexAntiLvsMass->Fill(lDcaPosToPrimVertex,lInvMassAntiLambda);
-           fHistDcaNegToPrimVertexAntiLvsMass->Fill(lDcaNegToPrimVertex,lInvMassAntiLambda);
-           fHistRadiusV0AntiLvsMass->Fill(lV0Radius,lInvMassAntiLambda);
-           fHistDecayLengthV0AntiLvsMass->Fill(lV0DecayLength,lInvMassAntiLambda);
-           fHistDcaV0DaughtersAntiLvsMass->Fill(lDcaV0Daughters,lInvMassAntiLambda);
-           fHistCosPointAngleAntiLvsMass->Fill(lV0cosPointAngle,lInvMassAntiLambda);
-
-           if(fQASelector){            
-             if (lPtAntiLambda>0 && lPtAntiLambda <3){ 
-               fHistDcaPosToPrimVertexAntiLVsMasspt1->Fill(lDcaPosToPrimVertex,lInvMassAntiLambda);
-               fHistDcaNegToPrimVertexAntiLVsMasspt1->Fill(lDcaNegToPrimVertex,lInvMassAntiLambda);
-               fHistRadiusV0AntiLVsMasspt1->Fill(lV0Radius,lInvMassAntiLambda);
-               fHistDecayLengthV0AntiLVsMasspt1->Fill(lV0DecayLength,lInvMassAntiLambda);
-               fHistDcaV0DaughtersAntiLVsMasspt1->Fill(lDcaV0Daughters,lInvMassAntiLambda);
-               fHistCosPointAngleAntiLVsMasspt1->Fill(lV0cosPointAngle,lInvMassAntiLambda);
-             }   
-             if (lPtAntiLambda > 3 && lPtAntiLambda < 6){ 
-               fHistDcaPosToPrimVertexAntiLVsMasspt2->Fill(lDcaPosToPrimVertex,lInvMassAntiLambda);
-               fHistDcaNegToPrimVertexAntiLVsMasspt2->Fill(lDcaNegToPrimVertex,lInvMassAntiLambda);
-               fHistRadiusV0AntiLVsMasspt2->Fill(lV0Radius,lInvMassAntiLambda);
-               fHistDecayLengthV0AntiLVsMasspt2->Fill(lV0DecayLength,lInvMassAntiLambda);
-               fHistDcaV0DaughtersAntiLVsMasspt2->Fill(lDcaV0Daughters,lInvMassAntiLambda);
-               fHistCosPointAngleAntiLVsMasspt2->Fill(lV0cosPointAngle,lInvMassAntiLambda);
-             }   
-             if (lPtAntiLambda > 6 && lPtAntiLambda < 10){ 
-               fHistDcaPosToPrimVertexAntiLVsMasspt3->Fill(lDcaPosToPrimVertex,lInvMassAntiLambda);
-               fHistDcaNegToPrimVertexAntiLVsMasspt3->Fill(lDcaNegToPrimVertex,lInvMassAntiLambda);
-               fHistRadiusV0AntiLVsMasspt3->Fill(lV0Radius,lInvMassAntiLambda);
-               fHistDecayLengthV0AntiLVsMasspt3->Fill(lV0DecayLength,lInvMassAntiLambda);
-               fHistDcaV0DaughtersAntiLVsMasspt3->Fill(lDcaV0Daughters,lInvMassAntiLambda);
-               fHistCosPointAngleAntiLVsMasspt3->Fill(lV0cosPointAngle,lInvMassAntiLambda);
-             }   
-           } //QA Selector
          }
        } //end of Rap condition
       } // end of PID condition
@@ -3318,72 +1654,16 @@ void AliAnalysisTaskPerformanceStrange::UserExec(Option_t *)
     ///////////////////////////////////////////////////    
     // K0s associated histograms in |rap| < lCutRap:
     ///////////////////////////////////////////////////
-
-
-
-
     if (TMath::Abs(lRapK0s) < lCutRap) {
-       
       switch (lOnFlyStatus){
       case 0 : 
-       if(lCheckMcK0Short) {
-         fHistAsMcTranscTauK0s->Fill(lcTauK0s);
-             if (TMath::Abs(lRapK0s) < 0.3)
-              fHistAsMcTranscTauK0sRap3->Fill(lcTauK0s);
-        if (lPtK0s >= 3 && lPtK0s <= 6 )
-          {
-            fHistAsMcTranscTauK0sptLow->Fill(lcTauK0s);
-             if (TMath::Abs(lRapK0s) < 0.3)
-              fHistAsMcTranscTauK0sptLowRap3->Fill(lcTauK0s);
-          }
-       }   
        if (lcTauK0s< cutcTauK0s) {
-         fHistNsigmaPosPionK0->Fill(nSigmaPosPion);
-         fHistNsigmaNegPionK0->Fill(nSigmaNegPion);
-
-
-         //      fHistMultVsPtVsMassK0->Fill(multiplicity ,lInvMassK0s,lPtK0s);
          if(lCheckPIdK0Short) fHistPidMcMassK0->Fill(lInvMassK0s);
          if(lCheckMcK0Short) {
            fHistAsMcMassK0->Fill(lInvMassK0s);
            fHistAsMcPtK0->Fill(lPtK0s);
            fHistAsMcPtVsMassK0->Fill(lInvMassK0s,lPtK0s);
-
-
-           if (TMath::Abs(lRapK0s) < 0.3){
-             fHistAsMcMassK0Rap3->Fill(lInvMassK0s);
-             fHistAsMcPtK0Rap3->Fill(lPtK0s);
-             fHistAsMcPtVsMassK0Rap3->Fill(lInvMassK0s,lPtK0s);
-           }
-
-
-           if (lPtK0s <= 1) fHistAsMcPtZoomK0->Fill(lPtK0s);
-           fHistAsMcMassVsRadiusK0->Fill(rcPosRK0s,lInvMassK0s);
-           fHistAsMcResxK0->Fill(rcPosXK0s-mcPosX);
-           fHistAsMcResyK0->Fill(rcPosYK0s-mcPosY);
-           fHistAsMcReszK0->Fill(rcPosZK0s-mcPosZ);
-           fHistAsMcResrVsRadiusK0->Fill(rcPosRK0s,rcPosRK0s-mcPosR);
-           fHistAsMcReszVsRadiusK0->Fill(rcPosZK0s,rcPosZK0s-mcPosZ);
-           fHistAsMcProdRadiusK0->Fill(mcPosMotherR);
-           fHistAsMcProdRadiusXvsYK0s->Fill(mcPosMotherX,mcPosMotherY);
-           fHistAsMcResPtK0->Fill(deltaPtK0s);
-           fHistAsMcResPtVsRapK0->Fill(deltaPtK0s,lRapK0s);
-           fHistAsMcResPtVsPtK0->Fill(deltaPtK0s,lPtK0s);
          }
-         else if (lCheckSecondaryK0s) {
-           fHistAsMcSecondaryPtVsRapK0s->Fill(lPtK0s,lRapK0s);
-           fHistAsMcSecondaryProdRadiusK0s->Fill(mcPosMotherR);
-           fHistAsMcSecondaryProdRadiusXvsYK0s->Fill(mcPosMotherX,mcPosMotherY);
-           switch (lPdgcodeMotherOfMother) {
-           case 130   : fHistAsMcSecondaryMotherPdgCodeK0s->Fill(0.5);break; // K0L
-           case 321   : fHistAsMcSecondaryMotherPdgCodeK0s->Fill(1.5);break; // K+
-           case -321  : fHistAsMcSecondaryMotherPdgCodeK0s->Fill(2.5);break; // K-
-           case -3122 : fHistAsMcSecondaryMotherPdgCodeK0s->Fill(3.5);break; //AntiLambda
-           default    : fHistAsMcSecondaryMotherPdgCodeK0s->Fill(6.5);break;
-           }
-         }
-         break;
-         
        }
       } // end rapidity condition
     } // end ctau cut
@@ -3391,104 +1671,18 @@ void AliAnalysisTaskPerformanceStrange::UserExec(Option_t *)
     ///////////////////////////////////////////////////
     // Associated Lambda histograms in |rap| < lCutRap
     ////////////////////////////////////////////////////
-
-
-
      if ((LambdaPID==1 && lMomInnerWallPos <=1) || (lMomInnerWallPos>1) ||  !(fUsePID.Contains("withPID"))){  
-
     if (TMath::Abs(lRapLambda) < lCutRap) {
       switch (lOnFlyStatus){
       case 0 : 
-       if(lCheckMcLambda) {
-         fHistAsMcTranscTauL->Fill(lcTauLambda);
-             if (TMath::Abs(lRapLambda) < 0.3)
-              fHistAsMcTranscTauLRap3->Fill(lcTauLambda);
-        if (lPtLambda >= 3 && lPtLambda <= 6 )
-          {
-            fHistAsMcTranscTauLptLow->Fill(lcTauLambda);
-             if (TMath::Abs(lRapLambda) < 0.3)
-              fHistAsMcTranscTauLptLowRap3->Fill(lcTauLambda);
-          }
-       }   
-
        if (lcTauLambda < cutcTauL){    
-         fHistNsigmaPosProtonLambda->Fill(nSigmaPosProton);
-         fHistNsigmaNegPionLambda->Fill(nSigmaNegPion);
-
-         //          fHistMultVsPtVsMassLambda->Fill(multiplicity ,lInvMassLambda,lPtLambda);
          if(lCheckPIdLambda) fHistPidMcMassLambda->Fill(lInvMassLambda);
-         
          if(lCheckMcLambda) {
-
            fHistAsMcMassLambda->Fill(lInvMassLambda);
            fHistAsMcPtLambda->Fill(lPtLambda);
            //    fHistCosPointAngleLvsMassVsPtsigL->Fill(lPtLambda,lV0cosPointAngle,lInvMassLambda);
            fHistAsMcPtVsMassLambda->Fill(lInvMassLambda,lPtLambda);
-
-
-           if (TMath::Abs(lRapLambda) < 0.3) {
-             fHistAsMcMassLambdaRap3->Fill(lInvMassLambda);
-             fHistAsMcPtLambdaRap3->Fill(lPtLambda);
-             fHistAsMcPtVsMassLambdaRap3->Fill(lInvMassLambda,lPtLambda);
-           }
-
-
-           if (lPtLambda <= 1) fHistAsMcPtZoomLambda->Fill(lPtLambda);
-           fHistAsMcMassVsRadiusLambda->Fill(rcPosRLambda,lInvMassLambda);
-           fHistAsMcResxLambda->Fill(rcPosXLambda-mcPosX);
-           fHistAsMcResyLambda->Fill(rcPosYLambda-mcPosY);
-           fHistAsMcReszLambda->Fill(rcPosZLambda-mcPosZ);
-           fHistAsMcResrVsRadiusLambda->Fill(rcPosRLambda,rcPosRLambda-mcPosR);
-           fHistAsMcReszVsRadiusLambda->Fill(rcPosZLambda,rcPosZLambda-mcPosZ);
-           fHistAsMcProdRadiusLambda->Fill(mcPosMotherR);
-           fHistAsMcProdRadiusXvsYLambda->Fill(mcPosMotherX,mcPosMotherY);
-           fHistAsMcResPtLambda->Fill(deltaPtLambda);
-           fHistAsMcResPtVsRapLambda->Fill(deltaPtLambda,lRapLambda);
-           fHistAsMcResPtVsPtLambda->Fill(deltaPtLambda,lPtLambda);
-           if (lComeFromSigma) fHistAsMcPtLambdaFromSigma->Fill(lPtLambda);
-           switch (lPdgcodeMotherOfMother) {
-           case 3222 : fHistAsMcMotherPdgCodeLambda->Fill(0.5); break; // Sigma +
-           case 3212 : fHistAsMcMotherPdgCodeLambda->Fill(1.5); break; // Sigma 0
-           case 3112 : fHistAsMcMotherPdgCodeLambda->Fill(2.5); break;// Sigma -
-           case 3224 : fHistAsMcMotherPdgCodeLambda->Fill(3.5); break;// Sigma(1385) +
-           case 3214 : fHistAsMcMotherPdgCodeLambda->Fill(4.5); break;// Sigma(1385) 0
-           case 3114 : fHistAsMcMotherPdgCodeLambda->Fill(5.5); break;// Sigma(1385) -
-           case 3322 : fHistAsMcMotherPdgCodeLambda->Fill(6.5); break; // Xi 0
-           case 3312 : fHistAsMcMotherPdgCodeLambda->Fill(7.5); break; // Xi -
-           case 3334 : fHistAsMcMotherPdgCodeLambda->Fill(8.5); break; // Omega
-           case -1   : fHistAsMcMotherPdgCodeLambda->Fill(9.5); break;
-           default   : fHistAsMcMotherPdgCodeLambda->Fill(10.5);break; 
-           }
-   
-           //printf("found Lambda RC dcaPos=%e dcaNeg=%e dcaDau=%e cosP=%e pT=%e mass=%e\n",lDcaPosToPrimVertex ,lDcaNegToPrimVertex ,lDcaV0Daughters,lV0cosPointAngle,lPtLambda,lInvMassLambda);
-           //printf("found Lambda RC Pindex=%d  Nindex=%d  Plabel=%d  Nlabel=%d\n\n",lIndexTrackPos,lIndexTrackNeg,lLabelTrackPos,lLabelTrackNeg);
-           
          }
-         
-         else if (lCheckSecondaryLambda) {
-           fHistAsMcSecondaryPtVsRapLambda->Fill(lPtLambda,lRapLambda);
-           fHistAsMcSecondaryProdRadiusLambda->Fill(mcPosMotherR); 
-           fHistAsMcSecondaryProdRadiusXvsYLambda->Fill(mcPosMotherX,mcPosMotherY);
-           if (lComeFromSigma) fHistAsMcSecondaryPtLambdaFromSigma->Fill(lPtLambda);
-           printf(" lPdgcodeMotherOfMother= %d",lPdgcodeMotherOfMother);
-           switch (lPdgcodeMotherOfMother) {
-           case 3222 : fHistAsMcSecondaryMotherPdgCodeLambda->Fill(0.5); break;// Sigma +
-           case 3212 : fHistAsMcSecondaryMotherPdgCodeLambda->Fill(1.5); break;// Sigma 0
-           case 3112 : fHistAsMcSecondaryMotherPdgCodeLambda->Fill(2.5); break;// Sigma -
-           case 3224 : fHistAsMcSecondaryMotherPdgCodeLambda->Fill(3.5); break;// Sigma(1385) +
-           case 3214 : fHistAsMcSecondaryMotherPdgCodeLambda->Fill(4.5); break;// Sigma(1385) 0
-           case 3114 : fHistAsMcSecondaryMotherPdgCodeLambda->Fill(5.5); break;// Sigma(1385) -
-           case 3322 : fHistAsMcSecondaryMotherPdgCodeLambda->Fill(6.5); break; // Xi 0
-           case 3312 : fHistAsMcSecondaryMotherPdgCodeLambda->Fill(7.5); break; // Xi -
-           case 3334 : fHistAsMcSecondaryMotherPdgCodeLambda->Fill(8.5); break; // Omega
-           case -1   : fHistAsMcSecondaryMotherPdgCodeLambda->Fill(9.5); break;
-           default   : fHistAsMcSecondaryMotherPdgCodeLambda->Fill(10.5);break;
-           }
-         }
-         
-         //         if(!lCheckMcLambda)fHistCosPointAngleLvsMassVsPtbackL->Fill(lPtLambda,lV0cosPointAngle,lInvMassLambda);
-         break;
-         
        }
       } // end rapidity condition
      }// end PID condition
@@ -3497,100 +1691,17 @@ void AliAnalysisTaskPerformanceStrange::UserExec(Option_t *)
     ////////////////////////////////////////////////////////
     // Associated AntiLambda histograms in |rap| < lCutRap
     ////////////////////////////////////////////////////////
-
       if ((AntiLambdaPID==1 && lMomInnerWallNeg <=1) || (lMomInnerWallNeg>1) ||  !(fUsePID.Contains("withPID"))){          
-
     if (TMath::Abs(lRapAntiLambda) < lCutRap) {
       switch (lOnFlyStatus){
       case 0 : 
-       if(lCheckMcAntiLambda) {
-         fHistAsMcTranscTauAntiL->Fill(lcTauAntiLambda);
-             if (TMath::Abs(lRapAntiLambda) < 0.3)
-              fHistAsMcTranscTauAntiLRap3->Fill(lcTauAntiLambda);
-        if (lPtAntiLambda >= 3 && lPtAntiLambda <= 6 )
-          {
-            fHistAsMcTranscTauAntiLptLow->Fill(lcTauAntiLambda);
-             if (TMath::Abs(lRapAntiLambda) < 0.3)
-              fHistAsMcTranscTauAntiLptLowRap3->Fill(lcTauAntiLambda);
-          }
-       }   
-
        if (lcTauAntiLambda < cutcTauL){
-         fHistNsigmaPosProtonAntiLambda->Fill(nSigmaPosProton);
-         fHistNsigmaNegPionAntiLambda->Fill(nSigmaNegPion);
-
-         //          fHistMultVsPtVsMassAntiLambda->Fill(multiplicity ,lInvMassAntiLambda,lPtAntiLambda);
          if(lCheckPIdAntiLambda) fHistPidMcMassAntiLambda->Fill(lInvMassAntiLambda);
-         
          if(lCheckMcAntiLambda) {
-
            fHistAsMcMassAntiLambda->Fill(lInvMassAntiLambda);
            fHistAsMcPtAntiLambda->Fill(lPtAntiLambda);
            fHistAsMcPtVsMassAntiLambda->Fill(lInvMassAntiLambda,lPtAntiLambda);
-
-
-           if (TMath::Abs(lRapAntiLambda) < 0.3) {
-             fHistAsMcMassAntiLambdaRap3->Fill(lInvMassAntiLambda);
-             fHistAsMcPtAntiLambdaRap3->Fill(lPtAntiLambda);
-             fHistAsMcPtVsMassAntiLambdaRap3->Fill(lInvMassAntiLambda,lPtAntiLambda);
-           }
-
-
-           if (lPtAntiLambda <= 1) fHistAsMcPtZoomAntiLambda->Fill(lPtAntiLambda);
-           fHistAsMcMassVsRadiusAntiLambda->Fill(rcPosRAntiLambda,lInvMassAntiLambda);
-           fHistAsMcResxAntiLambda->Fill(rcPosXAntiLambda-mcPosX);
-           fHistAsMcResyAntiLambda->Fill(rcPosYAntiLambda-mcPosY);
-           fHistAsMcReszAntiLambda->Fill(rcPosZAntiLambda-mcPosZ);
-           fHistAsMcResrVsRadiusAntiLambda->Fill(rcPosRAntiLambda,rcPosRAntiLambda-mcPosR);
-           fHistAsMcReszVsRadiusAntiLambda->Fill(rcPosZAntiLambda,rcPosZAntiLambda-mcPosZ);
-           fHistAsMcProdRadiusAntiLambda->Fill(mcPosMotherR);
-           fHistAsMcProdRadiusXvsYAntiLambda->Fill(mcPosMotherX,mcPosMotherY);
-           fHistAsMcResPtAntiLambda->Fill(deltaPtAntiLambda);
-           fHistAsMcResPtVsRapAntiLambda->Fill(deltaPtAntiLambda,lRapAntiLambda);
-           fHistAsMcResPtVsPtAntiLambda->Fill(deltaPtAntiLambda,lPtAntiLambda);
-           if (lComeFromSigma) fHistAsMcPtAntiLambdaFromSigma->Fill(lPtAntiLambda);
-           switch (lPdgcodeMotherOfMother) {
-           case 3222 : fHistAsMcMotherPdgCodeAntiLambda->Fill(0.5); break; // Sigma +
-           case 3212 : fHistAsMcMotherPdgCodeAntiLambda->Fill(1.5); break; // Sigma 0
-           case 3112 : fHistAsMcMotherPdgCodeAntiLambda->Fill(2.5); break;// Sigma -
-           case 3224 : fHistAsMcMotherPdgCodeAntiLambda->Fill(3.5); break;// Sigma(1385) +
-           case 3214 : fHistAsMcMotherPdgCodeAntiLambda->Fill(4.5); break;// Sigma(1385) 0
-           case 3114 : fHistAsMcMotherPdgCodeAntiLambda->Fill(5.5); break;// Sigma(1385) -
-           case 3322 : fHistAsMcMotherPdgCodeAntiLambda->Fill(6.5); break; // Xi 0
-           case 3312 : fHistAsMcMotherPdgCodeAntiLambda->Fill(7.5); break; // Xi -
-           case 3334 : fHistAsMcMotherPdgCodeAntiLambda->Fill(8.5); break; // Omega
-
-           case -1   : fHistAsMcMotherPdgCodeAntiLambda->Fill(9.5); break;
-           default   : fHistAsMcMotherPdgCodeAntiLambda->Fill(10.5);break; 
-           }
-    
-           //printf("found AntiLambda RC dcaPos=%e dcaNeg=%e dcaDau=%e cosP=%e pT=%e mass=%e\n",lDcaPosToPrimVertex ,lDcaNegToPrimVertex ,lDcaV0Daughters,lV0cosPointAngle,lPtAntiLambda,lInvMassAntiLambda);
-           //printf("found AntiLambda RC Pindex=%d  Nindex=%d  Plabel=%d  Nlabel=%d\n\n",lIndexTrackPos,lIndexTrackNeg,lLabelTrackPos,lLabelTrackNeg);
-           
          }
-         
-         else if (lCheckSecondaryAntiLambda) {
-           fHistAsMcSecondaryPtVsRapAntiLambda->Fill(lPtAntiLambda,lRapAntiLambda);
-           fHistAsMcSecondaryProdRadiusAntiLambda->Fill(mcPosMotherR); 
-           fHistAsMcSecondaryProdRadiusXvsYAntiLambda->Fill(mcPosMotherX,mcPosMotherY);
-           if (lComeFromSigma) fHistAsMcSecondaryPtAntiLambdaFromSigma->Fill(lPtAntiLambda);
-           printf(" lPdgcodeMotherOfMother= %d",lPdgcodeMotherOfMother);
-           switch (lPdgcodeMotherOfMother) {
-           case 3222 : fHistAsMcSecondaryMotherPdgCodeAntiLambda->Fill(0.5); break;// Sigma +
-           case 3212 : fHistAsMcSecondaryMotherPdgCodeAntiLambda->Fill(1.5); break;// Sigma 0
-           case 3112 : fHistAsMcSecondaryMotherPdgCodeAntiLambda->Fill(2.5); break;// Sigma -
-           case 3224 : fHistAsMcSecondaryMotherPdgCodeAntiLambda->Fill(3.5); break;// Sigma(1385) +
-           case 3214 : fHistAsMcSecondaryMotherPdgCodeAntiLambda->Fill(4.5); break;// Sigma(1385) 0
-           case 3114 : fHistAsMcSecondaryMotherPdgCodeAntiLambda->Fill(5.5); break;// Sigma(1385) -
-           case 3322 : fHistAsMcSecondaryMotherPdgCodeAntiLambda->Fill(6.5); break; // Xi 0
-           case 3312 : fHistAsMcSecondaryMotherPdgCodeAntiLambda->Fill(7.5); break; // Xi -
-           case 3334 : fHistAsMcSecondaryMotherPdgCodeAntiLambda->Fill(8.5); break; // Omega
-           case -1   : fHistAsMcSecondaryMotherPdgCodeAntiLambda->Fill(9.5); break;
-           default   : fHistAsMcSecondaryMotherPdgCodeAntiLambda->Fill(10.5);break;
-           }
-         }
-         break;
-         
        }
       } // end rapidity condition
      }// end PID condition     
@@ -3598,14 +1709,11 @@ void AliAnalysisTaskPerformanceStrange::UserExec(Option_t *)
   
   } // end V0 loop
 
-  //  if (primaryVtxKF) delete primaryVtxKF;primaryVtxKF=NULL ;
-
-
   fHistV0Multiplicity->Fill(nv0s);
   //  fHistV0MultiplicityMI->Fill(nv0sMI);
 
   if (fAnalysisType == "ESD") { if(myPrimaryVertex) delete myPrimaryVertex; }
-  if (fAnalysisType == "ESD") { if(TestTrackCuts) delete TestTrackCuts; }
+  //  if (fAnalysisType == "ESD") { if(TestTrackCuts) delete TestTrackCuts; }
 
   
   // Post output data