/**************************************************************************
 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/

// +-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
//
// Modified version of AliAnalysisTaskCheckCascade.h
// Used bits of code from AliAnalysisTaskCheckPerformanceStrange
//
// --- David Dobrigkeit Chinellato
//
// +-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+

#ifndef ALIANALYSISTASKEXTRACTPERFORMANCEV0_H
#define ALIANALYSISTASKEXTRACTPERFORMANCEV0_H

class TList;
class TH1F;
class TH2F;
class TH3F;
class TVector3;
class THnSparse;

class AliESDpid;
class AliESDtrackCuts;
class AliAnalysisUtils;
class AliESDEvent;
class AliPhysicsSelection;
class AliCFContainer;

//#include "TString.h"
//#include "AliESDtrackCuts.h"
//#include "AliAnalysisTaskSE.h"

class AliAnalysisTaskExtractPerformanceV0 : public AliAnalysisTaskSE {
 public:
  AliAnalysisTaskExtractPerformanceV0();
  AliAnalysisTaskExtractPerformanceV0(const char *name);
  virtual ~AliAnalysisTaskExtractPerformanceV0();
  
  virtual void   UserCreateOutputObjects();
  virtual void   UserExec(Option_t *option);
  virtual void   Terminate(Option_t *);
  Double_t MyRapidity(Double_t rE, Double_t rPz) const;
  void CheckChargeV0(AliESDv0 *thisv0);

  void SetIsNuclear           (Bool_t lIsNuclear   = kTRUE ) { fkIsNuclear   = lIsNuclear;   }
  void SetINT7Trigger         (Bool_t lSwitchINT7  = kTRUE ) { fkSwitchINT7   = lSwitchINT7; }
  void SetUseOnTheFly         (Bool_t lUseOnTheFly = kTRUE ) { fkUseOnTheFly = lUseOnTheFly; }
  void SetTakeAllTracks       (Bool_t lTakeAllTracks = kTRUE ) { fkTakeAllTracks = lTakeAllTracks; }
  void SetpARapidityShift     (Double_t lRapShift = 0.465 ) { fpArapidityShift = lRapShift; }
  void SetCentralityEstimator (TString lCentralityEstimator = "V0M" ) { fCentralityEstimator = lCentralityEstimator; }
  void SetLightWeightAnalysis (Bool_t lLightWeight = kTRUE) {fkLightWeight = lLightWeight;  }
  void SetFastOnly (TString lFastOnly = "kFastOnly") {fkFastOnly = lFastOnly;  }
  void SetpAVertexSelection   (Bool_t lpAVertexSelection = kTRUE) {fkpAVertexSelection = lpAVertexSelection;  }
  void SetRunV0Vertexer       ( Bool_t lRunV0Vertexer = kTRUE) { fkRunV0Vertexer = lRunV0Vertexer; }
  void SetRejectPileup        ( Bool_t lRejectPileup = kTRUE) { fkRejectPileup = lRejectPileup; }
  void SetSpecialExecution    ( Bool_t lSpecialExecution = kTRUE) { fkSpecialExecution = lSpecialExecution; }
  void SetSaveAssociatedOnly  ( Bool_t lSaveAssocOnly = kTRUE) { fkSaveAssociatedOnly = lSaveAssocOnly; }
  void SetSkipTrigger         ( Bool_t lSkipTrigger = kTRUE ){ fkSkipTrigger = lSkipTrigger; }
  void SetDoNotCallTPCdEdx    ( Bool_t lDoNotCallTPCdEdx = kTRUE) { fkDoNotCallTPCdEdx = lDoNotCallTPCdEdx; }
  
  //---------------------------------------------------------------------------------------
  //Setters for the V0 Vertexer Parameters
  void SetV0VertexerMaxChisquare   ( Double_t lParameter ){ fV0Sels[0] = lParameter; }
  void SetV0VertexerDCAFirstToPV   ( Double_t lParameter ){ fV0Sels[1] = lParameter; }
  void SetV0VertexerDCASecondtoPV  ( Double_t lParameter ){ fV0Sels[2] = lParameter; }
  void SetV0VertexerDCAV0Daughters ( Double_t lParameter ){ fV0Sels[3] = lParameter; }
  void SetV0VertexerCosinePA       ( Double_t lParameter ){ fV0Sels[4] = lParameter; }
  void SetV0VertexerMinRadius      ( Double_t lParameter ){ fV0Sels[5] = lParameter; }
  void SetV0VertexerMaxRadius      ( Double_t lParameter ){ fV0Sels[6] = lParameter; }
  //---------------------------------------------------------------------------------------
  void SetDiffractiveOnly ( Bool_t lDiffractiveOnly = kTRUE ) { fDiffractiveOnly = lDiffractiveOnly; }
  void SetEtaRefMult ( Double_t lEtaRefMult = 0.5 ) { fEtaRefMult = lEtaRefMult; }

  
 private:
        // Note : In ROOT, "//!" means "do not stream the data from Master node to Worker node" ...
        // your data member object is created on the worker nodes and streaming is not needed.
        // http://root.cern.ch/download/doc/11InputOutput.pdf, page 14
  TList  *fListHistV0;  //! List of Cascade histograms
  TTree  *fTree;              //! Output Tree, V0

  AliPIDResponse *fPIDResponse;     // PID response object
  AliESDtrackCuts *fESDtrackCuts;   // ESD track cuts used for primary track definition
  AliAnalysisUtils *fUtils;         // analysis utils (for pA vertex selection)


  //Objects Controlling Task Behaviour 
  
  Bool_t fkIsNuclear;   //if true, replace multiplicity est. by centrality (default FALSE) 
  Bool_t fkSwitchINT7 ; //if true, skip FASTOnly (default FALSE)
  Bool_t fkUseOnTheFly; //if true, will use On-the-fly V0s instead of Offline V0s (default FALSE)
  Bool_t fkTakeAllTracks; // if true, no TPC crossed rows and ratio cut
  Double_t fpArapidityShift; //pA rapidity shift (should be 0.465, usually)
  TString fCentralityEstimator; //Centrality Estimator String value (default V0M, DEPRECATED)
  Bool_t fkLightWeight; //if true, skip a number of debugging information branches in TTree
                        //(to make resulting tree output significantly smaller!
  TString fkFastOnly; //"" if no extra selection, "kFastOnly" -> without SDD, "NotkFastOnly" -> With SDD
  Bool_t fkpAVertexSelection; //if true, select vertex with pPb Methods
  Bool_t fkRunV0Vertexer; //if true, re-run vertexer with loose cuts. CARE MUST BE TAKEN in PbPb!
  Bool_t fkRejectPileup; //Reject pileup switch
  Bool_t fkSpecialExecution; //Special Exploratory Execution
  Bool_t fkSkipTrigger; //Skip trigger selection (use SetCollisionCandidates)
  Bool_t fkSaveAssociatedOnly; //Save Only associated V0s for smaller output
  Bool_t fkDoNotCallTPCdEdx; //For running on datasets that don't have BB curves... 
  
  Double_t        fV0Sels[7];                     // Array to store the 7 values for the different selections V0 related

  Bool_t fDiffractiveOnly; //Only look at diffractive generated events
  Bool_t fEtaRefMult; // Ref Mult eta range 

//===========================================================================================
//   Variables for Tree, V0s
//===========================================================================================
   Int_t    fTreeVariablePrimaryStatus;      //!
   Int_t    fTreeVariablePrimaryStatusMother;      //!
   Float_t fTreeVariableChi2V0;             //!
   Float_t fTreeVariableDcaV0Daughters; //!
   Float_t fTreeVariableDcaV0ToPrimVertex; //!
   Float_t fTreeVariableDcaPosToPrimVertex; //!
   Float_t fTreeVariableDcaNegToPrimVertex; //!
   Float_t fTreeVariableDCAV0ToPrimVertex; //!
   Float_t fTreeVariableV0CosineOfPointingAngle; //!
   Float_t fTreeVariableV0Radius; //!
   Float_t fTreeVariablePt; //!
   Float_t fTreeVariablePtMC; //!
   Float_t fTreeVariableRapK0Short; //!
   Float_t fTreeVariableRapLambda; //!
   Float_t fTreeVariableRapMC; //!
   Float_t fTreeVariableInvMassK0s; //!
   Float_t fTreeVariableInvMassLambda; //!
   Float_t fTreeVariableInvMassAntiLambda; //!
   Float_t fTreeVariableAlphaV0; //!
   Float_t fTreeVariablePtArmV0;//!
   Float_t fTreeVariableNegTotMomentum; //!               
   Float_t fTreeVariablePosTotMomentum; //!
   Float_t fTreeVariableNegTransvMomentum; //!   
   Float_t fTreeVariablePosTransvMomentum; //!
   Float_t fTreeVariableNegTransvMomentumMC; //!   
   Float_t fTreeVariablePosTransvMomentumMC; //!
   
   Float_t fTreeVariableNSigmasPosProton; //!
   Float_t fTreeVariableNSigmasPosPion; //! 
   Float_t fTreeVariableNSigmasNegProton; //!
   Float_t fTreeVariableNSigmasNegPion; //! 

   Float_t fTreeVariablePtMother; //!
   Float_t fTreeVariableV0CreationRadius; //!
   Int_t fTreeVariablePID; //!
   Int_t fTreeVariablePIDPositive; //!
   Int_t fTreeVariablePIDNegative; //!
   Int_t fTreeVariablePIDMother; //!
   Int_t fTreeVariableIndexStatus; //!
   Int_t fTreeVariableIndexStatusMother; //!

   Int_t   fTreeVariableRunNumber; //! 
   ULong64_t fTreeVariableEventNumber; //!

   //Note: TDistOverTotMom needs a mass hypothesis to be converted to proper decaylength.
   Float_t fTreeVariableDistOverTotMom;//!

   Float_t fTreeVariablePosEta; //!
   Float_t fTreeVariableNegEta; //!

   Float_t fTreeVariableVertexZ; //!

   Int_t fTreeVariableLeastNbrCrossedRows;//!
   Float_t fTreeVariableLeastRatioCrossedRowsOverFindable;//!
  
  Int_t fTreeVariableMultiplicity;//!
  Int_t fTreeVariableMultiplicityV0A;//!
  Int_t fTreeVariableMultiplicityZNA;//!
  Int_t fTreeVariableMultiplicityTRK;//!
  Int_t fTreeVariableMultiplicitySPD;//!
  
   Int_t fTreeVariableMultiplicityMC;//!

  Float_t fTreeVariableV0x; //!
  Float_t fTreeVariableV0y; //!
  Float_t fTreeVariableV0z; //!

  Float_t fTreeVariableV0Px; //!
  Float_t fTreeVariableV0Py; //!
  Float_t fTreeVariableV0Pz; //!

  Float_t fTreeVariableMCV0x; //!
  Float_t fTreeVariableMCV0y; //!
  Float_t fTreeVariableMCV0z; //!

  Float_t fTreeVariableMCV0Px; //!
  Float_t fTreeVariableMCV0Py; //!
  Float_t fTreeVariableMCV0Pz; //!

  Float_t fTreeVariablePVx; //!
  Float_t fTreeVariablePVy; //!
  Float_t fTreeVariablePVz; //!

  Float_t fTreeVariableMCPVx; //!
  Float_t fTreeVariableMCPVy; //!
  Float_t fTreeVariableMCPVz; //!

  Bool_t fTreeVariableIsNonInjected; //!
  //Decay Length issue debugging: ULong_t with track status
  ULong64_t fTreeVariableNegTrackStatus;
  ULong64_t fTreeVariablePosTrackStatus;
  
  //Physical Primary, Sec-Weak, Sec-Material -- debug only
  Int_t fTreeVariableNegPhysicalStatus;
  Int_t fTreeVariablePosPhysicalStatus;
  
  Float_t fTreeVariableNegPx; //!
  Float_t fTreeVariableNegPy; //!
  Float_t fTreeVariableNegPz; //!
  Float_t fTreeVariablePosPx; //!
  Float_t fTreeVariablePosPy; //!
  Float_t fTreeVariablePosPz; //!
  
//===========================================================================================
//   Histograms
//===========================================================================================

   TH1F      *fHistV0MultiplicityBeforeTrigSel;              //! V0 multiplicity distribution
   TH1F      *fHistV0MultiplicityForTrigEvt;                 //! V0 multiplicity distribution
   TH1F      *fHistV0MultiplicityForSelEvt;                  //! V0 multiplicity distribution
   TH1F      *fHistV0MultiplicityForSelEvtNoTPCOnly;         //! V0 multiplicity distribution
   TH1F      *fHistV0MultiplicityForSelEvtNoTPCOnlyNoPileup; //! V0 multiplicity distribution

  //Default V0M Centrality
   TH1F      *fHistMultiplicityBeforeTrigSel;     //! multiplicity distribution      
   TH1F      *fHistMultiplicityForTrigEvt;        //! multiplicity distribution
   TH1F      *fHistMultiplicity;                  //! multiplicity distribution
   TH1F      *fHistMultiplicityNoTPCOnly;         //! multiplicity distribution
   TH1F      *fHistMultiplicityNoTPCOnlyNoPileup; //! multiplicity distribution

  //V0A Centrality
  TH1F    *fHistMultiplicityV0ABeforeTrigSel; 	        //! multiplicity distribution
	TH1F    *fHistMultiplicityV0AForTrigEvt;  		        //! multiplicity distribution
	TH1F    *fHistMultiplicityV0A;     					        //! multiplicity distribution
	TH1F    *fHistMultiplicityV0ANoTPCOnly;			        //! multiplicity distribution
	TH1F    *fHistMultiplicityV0ANoTPCOnlyNoPileup;			//! multiplicity distribution
  
  //ZNA Centrality
  TH1F    *fHistMultiplicityZNABeforeTrigSel; 	        //! multiplicity distribution
	TH1F    *fHistMultiplicityZNAForTrigEvt;  		        //! multiplicity distribution
	TH1F    *fHistMultiplicityZNA;     					        //! multiplicity distribution
	TH1F    *fHistMultiplicityZNANoTPCOnly;			        //! multiplicity distribution
	TH1F    *fHistMultiplicityZNANoTPCOnlyNoPileup;			//! multiplicity distribution
  
  //TRK Centrality
  TH1F    *fHistMultiplicityTRKBeforeTrigSel; 	        //! multiplicity distribution
	TH1F    *fHistMultiplicityTRKForTrigEvt;  		        //! multiplicity distribution
	TH1F    *fHistMultiplicityTRK;     					        //! multiplicity distribution
	TH1F    *fHistMultiplicityTRKNoTPCOnly;			        //! multiplicity distribution
	TH1F    *fHistMultiplicityTRKNoTPCOnlyNoPileup;			//! multiplicity distribution
  
  //SPD Centrality
  TH1F    *fHistMultiplicitySPDBeforeTrigSel; 	        //! multiplicity distribution
	TH1F    *fHistMultiplicitySPDForTrigEvt;  		        //! multiplicity distribution
	TH1F    *fHistMultiplicitySPD;     					        //! multiplicity distribution
	TH1F    *fHistMultiplicitySPDNoTPCOnly;			        //! multiplicity distribution
	TH1F    *fHistMultiplicitySPDNoTPCOnlyNoPileup;			//! multiplicity distribution
  
  //Raw Data for J/Psi paper Technique
	TH2F    *f2dHistMultiplicityVsTrueBeforeTrigSel; 	        //! multiplicity distribution    
	TH2F    *f2dHistMultiplicityVsTrueForTrigEvt;  		        //! multiplicity distribution
	TH2F    *f2dHistMultiplicityVsTrue;     					        //! multiplicity distribution
	TH2F    *f2dHistMultiplicityVsTrueNoTPCOnly;			        //! multiplicity distribution
	TH2F    *f2dHistMultiplicityVsTrueNoTPCOnlyNoPileup;			//! multiplicity distribution

  //Raw Data for Vertex Z position estimator change
	TH2F    *f2dHistMultiplicityVsVertexZBeforeTrigSel; 	        //! multiplicity distribution    
	TH2F    *f2dHistMultiplicityVsVertexZForTrigEvt;  		        //! multiplicity distribution
	TH2F    *f2dHistMultiplicityVsVertexZ;     					        //! multiplicity distribution
	TH2F    *f2dHistMultiplicityVsVertexZNoTPCOnly;			        //! multiplicity distribution
	TH2F    *f2dHistMultiplicityVsVertexZNoTPCOnlyNoPileup;			//! multiplicity distribution

   TH1F      *fHistGenVertexZBeforeTrigSel;     //! multiplicity distribution      
   TH1F      *fHistGenVertexZForTrigEvt;        //! multiplicity distribution
   TH1F      *fHistGenVertexZ;                  //! multiplicity distribution
   TH1F      *fHistGenVertexZNoTPCOnly;         //! multiplicity distribution
   TH1F      *fHistGenVertexZNoTPCOnlyNoPileup; //! multiplicity distribution

//---> Filled At Analysis Scope

   TH3F      *f3dHistPrimAnalysisPtVsYVsMultLambda;     //! Lambda
   TH3F      *f3dHistPrimAnalysisPtVsYVsMultAntiLambda; //! AntiLambda
   TH3F      *f3dHistPrimAnalysisPtVsYVsMultK0Short;    //! K0Short

   TH3F      *f3dHistPrimAnalysisPtVsYCMSVsMultLambda;     //! Lambda
   TH3F      *f3dHistPrimAnalysisPtVsYCMSVsMultAntiLambda; //! AntiLambda
   TH3F      *f3dHistPrimAnalysisPtVsYCMSVsMultK0Short;    //! K0Short
  
//---> TRUE Multiplicity Containers
  
  TH3F      *f3dHistPrimAnalysisPtVsYVsMultMCLambda;     //! Lambda
  TH3F      *f3dHistPrimAnalysisPtVsYVsMultMCAntiLambda; //! AntiLambda
  TH3F      *f3dHistPrimAnalysisPtVsYVsMultMCK0Short;    //! K0Short
  
  TH3F      *f3dHistPrimAnalysisPtVsYCMSVsMultMCLambda;     //! Lambda
  TH3F      *f3dHistPrimAnalysisPtVsYCMSVsMultMCAntiLambda; //! AntiLambda
  TH3F      *f3dHistPrimAnalysisPtVsYCMSVsMultMCK0Short;    //! K0Short

//V0A containers
  
  TH3F      *f3dHistPrimAnalysisPtVsYVsMultV0ALambda;     //! Lambda
  TH3F      *f3dHistPrimAnalysisPtVsYVsMultV0AAntiLambda; //! AntiLambda
  TH3F      *f3dHistPrimAnalysisPtVsYVsMultV0AK0Short;    //! K0Short
  
  TH3F      *f3dHistPrimAnalysisPtVsYCMSVsMultV0ALambda;     //! Lambda
  TH3F      *f3dHistPrimAnalysisPtVsYCMSVsMultV0AAntiLambda; //! AntiLambda
  TH3F      *f3dHistPrimAnalysisPtVsYCMSVsMultV0AK0Short;    //! K0Short

//Cross-checking histograms: Charged Kaons (to compare with neutral ones at generator level) 
  TH3F      *f3dHistPrimRawPtVsYCMSVsMultV0AKPlus;     //! Added for cross-check of any bias
  TH3F      *f3dHistPrimRawPtVsYCMSVsMultV0AKMinus;    //! Added for cross-check of any bias
  TH3F      *f3dHistPrimAnalysisPtVsYCMSVsMultV0AKPlus;     //! Added for cross-check of any bias
  TH3F      *f3dHistPrimAnalysisPtVsYCMSVsMultV0AKMinus;    //! Added for cross-check of any bias

//ZNA containers
  
  TH3F      *f3dHistPrimAnalysisPtVsYVsMultZNALambda;     //! Lambda
  TH3F      *f3dHistPrimAnalysisPtVsYVsMultZNAAntiLambda; //! AntiLambda
  TH3F      *f3dHistPrimAnalysisPtVsYVsMultZNAK0Short;    //! K0Short
  
  TH3F      *f3dHistPrimAnalysisPtVsYCMSVsMultZNALambda;     //! Lambda
  TH3F      *f3dHistPrimAnalysisPtVsYCMSVsMultZNAAntiLambda; //! AntiLambda
  TH3F      *f3dHistPrimAnalysisPtVsYCMSVsMultZNAK0Short;    //! K0Short
  
//TRK containers
  
  TH3F      *f3dHistPrimAnalysisPtVsYVsMultTRKLambda;     //! Lambda
  TH3F      *f3dHistPrimAnalysisPtVsYVsMultTRKAntiLambda; //! AntiLambda
  TH3F      *f3dHistPrimAnalysisPtVsYVsMultTRKK0Short;    //! K0Short
  
  TH3F      *f3dHistPrimAnalysisPtVsYCMSVsMultTRKLambda;     //! Lambda
  TH3F      *f3dHistPrimAnalysisPtVsYCMSVsMultTRKAntiLambda; //! AntiLambda
  TH3F      *f3dHistPrimAnalysisPtVsYCMSVsMultTRKK0Short;    //! K0Short
  
//SPD containers
  
  TH3F      *f3dHistPrimAnalysisPtVsYVsMultSPDLambda;     //! Lambda
  TH3F      *f3dHistPrimAnalysisPtVsYVsMultSPDAntiLambda; //! AntiLambda
  TH3F      *f3dHistPrimAnalysisPtVsYVsMultSPDK0Short;    //! K0Short
  
  TH3F      *f3dHistPrimAnalysisPtVsYCMSVsMultSPDLambda;     //! Lambda
  TH3F      *f3dHistPrimAnalysisPtVsYCMSVsMultSPDAntiLambda; //! AntiLambda
  TH3F      *f3dHistPrimAnalysisPtVsYCMSVsMultSPDK0Short;    //! K0Short
  
//---> Containers for monte carlo information for calculating efficiency!

   TH3F      *f3dHistPrimRawPtVsYVsMultLambda;     //! Lambda
   TH3F      *f3dHistPrimRawPtVsYVsMultAntiLambda; //! AntiLambda
   TH3F      *f3dHistPrimRawPtVsYVsMultK0Short;    //! K0Short

   TH3F      *f3dHistPrimRawPtVsYCMSVsMultLambda;     //! Lambda
   TH3F      *f3dHistPrimRawPtVsYCMSVsMultAntiLambda; //! AntiLambda
   TH3F      *f3dHistPrimRawPtVsYCMSVsMultK0Short;    //! K0Short

//V0A Containers
  
  TH3F      *f3dHistPrimRawPtVsYVsMultV0ALambda;     //! Lambda
  TH3F      *f3dHistPrimRawPtVsYVsMultV0AAntiLambda; //! AntiLambda
  TH3F      *f3dHistPrimRawPtVsYVsMultV0AK0Short;    //! K0Short
  
  TH3F      *f3dHistPrimRawPtVsYCMSVsMultV0ALambda;     //! Lambda
  TH3F      *f3dHistPrimRawPtVsYCMSVsMultV0AAntiLambda; //! AntiLambda
  TH3F      *f3dHistPrimRawPtVsYCMSVsMultV0AK0Short;    //! K0Short

//ZNA Containers
  
  TH3F      *f3dHistPrimRawPtVsYVsMultZNALambda;     //! Lambda
  TH3F      *f3dHistPrimRawPtVsYVsMultZNAAntiLambda; //! AntiLambda
  TH3F      *f3dHistPrimRawPtVsYVsMultZNAK0Short;    //! K0Short
  
  TH3F      *f3dHistPrimRawPtVsYCMSVsMultZNALambda;     //! Lambda
  TH3F      *f3dHistPrimRawPtVsYCMSVsMultZNAAntiLambda; //! AntiLambda
  TH3F      *f3dHistPrimRawPtVsYCMSVsMultZNAK0Short;    //! K0Short
  
//TRK Containers
  
  TH3F      *f3dHistPrimRawPtVsYVsMultTRKLambda;     //! Lambda
  TH3F      *f3dHistPrimRawPtVsYVsMultTRKAntiLambda; //! AntiLambda
  TH3F      *f3dHistPrimRawPtVsYVsMultTRKK0Short;    //! K0Short
  
  TH3F      *f3dHistPrimRawPtVsYCMSVsMultTRKLambda;     //! Lambda
  TH3F      *f3dHistPrimRawPtVsYCMSVsMultTRKAntiLambda; //! AntiLambda
  TH3F      *f3dHistPrimRawPtVsYCMSVsMultTRKK0Short;    //! K0Short
  
//SPD Containers
  
  TH3F      *f3dHistPrimRawPtVsYVsMultSPDLambda;     //! Lambda
  TH3F      *f3dHistPrimRawPtVsYVsMultSPDAntiLambda; //! AntiLambda
  TH3F      *f3dHistPrimRawPtVsYVsMultSPDK0Short;    //! K0Short
  
  TH3F      *f3dHistPrimRawPtVsYCMSVsMultSPDLambda;     //! Lambda
  TH3F      *f3dHistPrimRawPtVsYCMSVsMultSPDAntiLambda; //! AntiLambda
  TH3F      *f3dHistPrimRawPtVsYCMSVsMultSPDK0Short;    //! K0Short
  
//Miscellaneous checking containers
  
   TH3F      *f3dHistPrimRawPtVsYVsMultNonInjLambda;     //! Non-injected Lambda
   TH3F      *f3dHistPrimRawPtVsYVsMultNonInjAntiLambda; //! Non-injected AntiLambda
   TH3F      *f3dHistPrimRawPtVsYVsMultNonInjK0Short;    //! Non-injected K0Short

   TH3F      *f3dHistPrimRawPtVsYVsMultMCLambda;     //! Lambda
   TH3F      *f3dHistPrimRawPtVsYVsMultMCAntiLambda; //! AntiLambda
   TH3F      *f3dHistPrimRawPtVsYVsMultMCK0Short;    //! K0Short

   TH3F      *f3dHistPrimRawPtVsYVsVertexZLambda;     //! Lambda
   TH3F      *f3dHistPrimRawPtVsYVsVertexZAntiLambda; //! AntiLambda
   TH3F      *f3dHistPrimRawPtVsYVsVertexZK0Short;    //! K0Short

   TH3F      *f3dHistPrimCloseToPVPtVsYVsMultLambda;     //! Lambda
   TH3F      *f3dHistPrimCloseToPVPtVsYVsMultAntiLambda; //! AntiLambda
   TH3F      *f3dHistPrimCloseToPVPtVsYVsMultK0Short;    //! K0Short

//---> Filled vs Decay Length

   TH3F      *f3dHistPrimRawPtVsYVsDecayLengthLambda;     //! Lambda
   TH3F      *f3dHistPrimRawPtVsYVsDecayLengthAntiLambda; //! AntiLambda
   TH3F      *f3dHistPrimRawPtVsYVsDecayLengthK0Short;    //! K0Short

//---> Needed for FeedDown Corrections

   TH3F      *f3dHistGenPtVsYVsMultXiMinus;      //! Generated Xi- Distrib
   TH3F      *f3dHistGenPtVsYVsMultXiPlus;       //! Generated Xi+ Distrib

   TH3F      *f3dHistGenSelectedPtVsYVsMultXiMinus;      //! Generated Xi- Distrib, at event selection level
   TH3F      *f3dHistGenSelectedPtVsYVsMultXiPlus;       //! Generated Xi+ Distrib, at event selection level

   TH3F      *f3dHistGenPtVsYCMSVsMultXiMinus;      //! Generated Xi- Distrib
   TH3F      *f3dHistGenPtVsYCMSVsMultXiPlus;       //! Generated Xi+ Distrib

   TH3F      *f3dHistGenSelectedPtVsYCMSVsMultXiMinus;      //! Generated Xi- Distrib, at event selection level
   TH3F      *f3dHistGenSelectedPtVsYCMSVsMultXiPlus;       //! Generated Xi+ Distrib, at event selection level
  
//V0A Containers

  TH3F      *f3dHistGenPtVsYVsMultV0AXiMinus;      //! Generated Xi- Distrib
  TH3F      *f3dHistGenPtVsYVsMultV0AXiPlus;       //! Generated Xi+ Distrib
  
  TH3F      *f3dHistGenSelectedPtVsYVsMultV0AXiMinus;      //! Generated Xi- Distrib, at event selection level
  TH3F      *f3dHistGenSelectedPtVsYVsMultV0AXiPlus;       //! Generated Xi+ Distrib, at event selection level
  
  TH3F      *f3dHistGenPtVsYCMSVsMultV0AXiMinus;      //! Generated Xi- Distrib
  TH3F      *f3dHistGenPtVsYCMSVsMultV0AXiPlus;       //! Generated Xi+ Distrib
  
  TH3F      *f3dHistGenSelectedPtVsYCMSVsMultV0AXiMinus;      //! Generated Xi- Distrib, at event selection level
  TH3F      *f3dHistGenSelectedPtVsYCMSVsMultV0AXiPlus;       //! Generated Xi+ Distrib, at event selection level

//ZNA Containers
  
  TH3F      *f3dHistGenPtVsYVsMultZNAXiMinus;      //! Generated Xi- Distrib
  TH3F      *f3dHistGenPtVsYVsMultZNAXiPlus;       //! Generated Xi+ Distrib
  
  TH3F      *f3dHistGenSelectedPtVsYVsMultZNAXiMinus;      //! Generated Xi- Distrib, at event selection level
  TH3F      *f3dHistGenSelectedPtVsYVsMultZNAXiPlus;       //! Generated Xi+ Distrib, at event selection level
  
  TH3F      *f3dHistGenPtVsYCMSVsMultZNAXiMinus;      //! Generated Xi- Distrib
  TH3F      *f3dHistGenPtVsYCMSVsMultZNAXiPlus;       //! Generated Xi+ Distrib
  
  TH3F      *f3dHistGenSelectedPtVsYCMSVsMultZNAXiMinus;      //! Generated Xi- Distrib, at event selection level
  TH3F      *f3dHistGenSelectedPtVsYCMSVsMultZNAXiPlus;       //! Generated Xi+ Distrib, at event selection level

//TRK Containers
  
  TH3F      *f3dHistGenPtVsYVsMultTRKXiMinus;      //! Generated Xi- Distrib
  TH3F      *f3dHistGenPtVsYVsMultTRKXiPlus;       //! Generated Xi+ Distrib
  
  TH3F      *f3dHistGenSelectedPtVsYVsMultTRKXiMinus;      //! Generated Xi- Distrib, at event selection level
  TH3F      *f3dHistGenSelectedPtVsYVsMultTRKXiPlus;       //! Generated Xi+ Distrib, at event selection level
  
  TH3F      *f3dHistGenPtVsYCMSVsMultTRKXiMinus;      //! Generated Xi- Distrib
  TH3F      *f3dHistGenPtVsYCMSVsMultTRKXiPlus;       //! Generated Xi+ Distrib
  
  TH3F      *f3dHistGenSelectedPtVsYCMSVsMultTRKXiMinus;      //! Generated Xi- Distrib, at event selection level
  TH3F      *f3dHistGenSelectedPtVsYCMSVsMultTRKXiPlus;       //! Generated Xi+ Distrib, at event selection level
  
//SPD Containers
  
  TH3F      *f3dHistGenPtVsYVsMultSPDXiMinus;      //! Generated Xi- Distrib
  TH3F      *f3dHistGenPtVsYVsMultSPDXiPlus;       //! Generated Xi+ Distrib
  
  TH3F      *f3dHistGenSelectedPtVsYVsMultSPDXiMinus;      //! Generated Xi- Distrib, at event selection level
  TH3F      *f3dHistGenSelectedPtVsYVsMultSPDXiPlus;       //! Generated Xi+ Distrib, at event selection level
  
  TH3F      *f3dHistGenPtVsYCMSVsMultSPDXiMinus;      //! Generated Xi- Distrib
  TH3F      *f3dHistGenPtVsYCMSVsMultSPDXiPlus;       //! Generated Xi+ Distrib
  
  TH3F      *f3dHistGenSelectedPtVsYCMSVsMultSPDXiMinus;      //! Generated Xi- Distrib, at event selection level
  TH3F      *f3dHistGenSelectedPtVsYCMSVsMultSPDXiPlus;       //! Generated Xi+ Distrib, at event selection level
  
  
   TH1F      *fHistPVx;                      //! PVx distrib
   TH1F      *fHistPVy;                      //! PVy distrib
   TH1F      *fHistPVz;                      //! PVz distrib
   TH1F      *fHistPVxAnalysis;                      //! PVx distrib
   TH1F      *fHistPVyAnalysis;                      //! PVy distrib
   TH1F      *fHistPVzAnalysis;                      //! PVz distrib
   TH1F      *fHistPVxAnalysisHasHighPtLambda;                      //! PVx distrib
   TH1F      *fHistPVyAnalysisHasHighPtLambda;                      //! PVy distrib
   TH1F      *fHistPVzAnalysisHasHighPtLambda;                      //! PVz distrib

   TH1F      *fHistSwappedV0Counter;                      //! Swapped v0 counter

   AliAnalysisTaskExtractPerformanceV0(const AliAnalysisTaskExtractPerformanceV0&);            // not implemented
   AliAnalysisTaskExtractPerformanceV0& operator=(const AliAnalysisTaskExtractPerformanceV0&); // not implemented
   
   ClassDef(AliAnalysisTaskExtractPerformanceV0, 11);
};

#endif