AliAnalysisTaskSELc2eleLambdafromAODtracks.h

Go to the documentation of this file.

#ifndef ALIANALYSISTASKSELC2ELELAMBDAFROMAODTRACKS_H
#define ALIANALYSISTASKSELC2ELELAMBDAFROMAODTRACKS_H

/**************************************************************************
 * Copyright(c) 1998-2009, 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.                  *
 **************************************************************************/

/* $Id$ */ 

#include "TROOT.h"
#include "TVector.h"
#include "TVector2.h"
#include "TSystem.h"
#include "TProfile.h"
#include <vector>

#include "AliAnalysisTaskSE.h"
#include "AliAODEvent.h"
#include "AliPID.h"
#include "AliRDHFCutsLctoeleLambdafromAODtracks.h"


class THnSparse;
class TH1F;
class TClonesArray;
class AliAODRecoCascadeHF;
class AliESDVertex;
class AliAODMCParticle;
class AliEventPoolManager;
class AliEventplane;
class AliNormalizationCounter;

class AliAnalysisTaskSELc2eleLambdafromAODtracks : public AliAnalysisTaskSE 
{
 public:
  AliAnalysisTaskSELc2eleLambdafromAODtracks();
  AliAnalysisTaskSELc2eleLambdafromAODtracks(const Char_t* name, AliRDHFCutsLctoeleLambdafromAODtracks* cuts, Bool_t writeVariableTree=kTRUE);
  virtual ~AliAnalysisTaskSELc2eleLambdafromAODtracks();

  virtual void UserCreateOutputObjects();
  virtual void Init();
  virtual void LocalInit() {Init();}
  virtual void UserExec(Option_t *option);
  virtual void Terminate(Option_t *option);

  void FillROOTObjects(AliAODRecoCascadeHF *elobj, AliAODv0 *v0, AliAODTrack *trk, AliAODTrack *trkpid, AliAODEvent *event, TClonesArray *mcArray);
  void FillMixROOTObjects(TLorentzVector *et, TLorentzVector *ev, TVector *tinfo, TVector *v0info2, Int_t charge);
  void FillElectronROOTObjects(AliAODTrack *trk, AliAODEvent *event, TClonesArray *mcArray);
  void FillV0ROOTObjects(AliAODv0 *v0, AliAODEvent *event, TClonesArray *mcArray);
  void FillMCROOTObjects(AliAODMCParticle *part, AliAODMCParticle *mcepart, AliAODMCParticle *mcv0part, Int_t decaytype);
  void FillMCEleROOTObjects(AliAODMCParticle *mcepart, TClonesArray *mcArray);
  void FillMCV0ROOTObjects(AliAODMCParticle *mcv0part, TClonesArray *mcArray);
  void FillMCGenPairROOTObjects(AliAODMCParticle *e, AliAODMCParticle *v, TClonesArray *array);
  Bool_t MakeMCAnalysis(TClonesArray *mcArray);
  void MakeAnalysis(AliAODEvent *aod, TClonesArray *mcArray);

  void SelectV0( const AliVEvent *event,Int_t nV0,Int_t &nSeleV0, Bool_t *seleV0Flags, TClonesArray *mcArray);
  void SelectTrack( const AliVEvent *event, Int_t trkEntries, Int_t &nSeleTrks,Bool_t *seleFlags, TClonesArray *mcArray);

  void SetMC(Bool_t theMCon) {fUseMCInfo = theMCon;}
  Bool_t GetMC() const {return fUseMCInfo;}
  void SetUseCentralityV0M(Bool_t centon) {fUseCentralityV0M = centon;}
  Bool_t GetUseCentralityV0M() const {return fUseCentralityV0M;}
  void SetUseCentralitySPDTracklet(Bool_t centon) {fUseCentralitySPDTracklet = centon;}
  Bool_t GetUseCentralitySPDTracklet() const {return fUseCentralitySPDTracklet;}
  void SetUseEventPlane(Int_t rpon) {fUseEventPlane = rpon;}
  Int_t GetUseEventPlane() const {return fUseEventPlane;}
  void SetWriteEachVariableTree(Bool_t a) {fWriteEachVariableTree = a;}
  Bool_t GetWriteEachVariableTree() const {return fWriteEachVariableTree;}
  void SetWriteMCVariableTree(Bool_t a) {fWriteMCVariableTree = a;}
  Bool_t GetWriteMCVariableTree() const {return fWriteMCVariableTree;}
  void SetMCEventType(Int_t theevt) {fMCEventType = theevt;}
  Int_t GetMCEventType() const {return fMCEventType;}
  void SetMCDoPairAnalysis(Bool_t a) {fMCDoPairAnalysis = a;}
  Bool_t GetMCDoPairAnalysis() const {return fMCDoPairAnalysis;}

  Float_t GetPhi0Pi(Float_t phi);
  Float_t GetEventPlaneForCandidate(AliAODTrack* d, AliAODv0 *v0, AliEventplane *pl,TVector2* q,TVector2* qsub1,TVector2* qsub2);

  void SetReconstructPrimVert(Bool_t a) { fReconstructPrimVert=a; }

  AliAODRecoCascadeHF* MakeCascadeHF(AliAODv0 *casc, AliAODTrack *trk, AliAODTrack *trkpid, AliAODEvent *aod, AliAODVertex *vert);
  AliAODVertex* ReconstructSecondaryVertex(AliAODv0 *casc, AliAODTrack *trk, AliAODEvent *aod);
  Int_t MatchToMC(AliAODRecoCascadeHF *elobj, TClonesArray *mcArray, Int_t *pdgele_array, Int_t *pdgv0_array, Int_t *labelele_array, Int_t *labelv0_array,  Int_t &ngen_ele, Int_t &ngen_v0);
  void  GetMCDecayHistory(AliAODMCParticle *mcpart, TClonesArray *mcArray, Int_t *pdgarray, Int_t *labelarray, Int_t &ngen);

  void StoreGlobalTrackReference(AliAODTrack *track, Int_t);
  void ResetGlobalTrackReference();

  // multiplicity dep analysis
  void SetMultiplVsZProfileLHC10b(TProfile* hprof){
    if(fMultEstimatorAvg[0]) delete fMultEstimatorAvg[0];
    fMultEstimatorAvg[0]=new TProfile(*hprof);
  }
  void SetMultiplVsZProfileLHC10c(TProfile* hprof){
    if(fMultEstimatorAvg[1]) delete fMultEstimatorAvg[1];
    fMultEstimatorAvg[1]=new TProfile(*hprof);
  }
  void SetMultiplVsZProfileLHC10d(TProfile* hprof){
    if(fMultEstimatorAvg[2]) delete fMultEstimatorAvg[2];
    fMultEstimatorAvg[2]=new TProfile(*hprof);
  }
  void SetMultiplVsZProfileLHC10e(TProfile* hprof){
    if(fMultEstimatorAvg[3]) delete fMultEstimatorAvg[3];
    fMultEstimatorAvg[3]=new TProfile(*hprof);
  }

  void SetReferenceMultiplcity(Double_t rmu){fRefMult=rmu;}

  void SetEventMixingWithPools(){fDoEventMixing=1;}
  void SetEventMixingOff(){fDoEventMixing=0;}
  void SetMixWithoutConversionFlag(Bool_t a){fMixWithoutConversionFlag=a;}
  void SetNumberOfEventsForMixing(Int_t events){fNumberOfEventsForMixing=events;}
  void SetPoolPVzBinLimits(Int_t Nzvtxbins,const Double_t *ZvtxBins){
    fNzVtxBins = Nzvtxbins;
    for(int ix = 0;ix<fNzVtxBins+1;ix++){fZvtxBins[ix] = ZvtxBins[ix];}
  }
  void SetPoolCentBinLimits(Int_t Ncentbins,const Double_t *CentBins){
    fNCentBins = Ncentbins;
    for(int ix = 0;ix<fNCentBins+1;ix++){fCentBins[ix] = CentBins[ix];}
  }
  void SetPoolRPBinLimits(Int_t Nrpbins,const Double_t *RPBins){
    fNRPBins = Nrpbins;
    for(int ix = 0;ix<fNRPBins+1;ix++){fRPBins[ix] = RPBins[ix];}
  }
  void DoEventMixingWithPools(Int_t index);
  void FillBackground(std::vector<TLorentzVector * > mixTypeE,std::vector<TVector * > mixTypeEVars, std::vector<TLorentzVector * > mixTypeL, std::vector<TVector * > mixTypeLVars, Int_t charge_v0pr);
  Int_t GetPoolIndex(Double_t zvert, Double_t mult, Double_t rp);


 private:

  AliAnalysisTaskSELc2eleLambdafromAODtracks(const AliAnalysisTaskSELc2eleLambdafromAODtracks &source);
  AliAnalysisTaskSELc2eleLambdafromAODtracks& operator=(const AliAnalysisTaskSELc2eleLambdafromAODtracks& source); 

  void DefineTreeVariables();
  void DefineEleTreeVariables();
  void DefineSingleTreeVariables();
  void DefineV0TreeVariables();
  void DefineMCTreeVariables();
  void DefineMCEleTreeVariables();
  void DefineMCV0TreeVariables();
  void DefineMCGenPairTreeVariables();
  void DefineCorrelationTreeVariables();
  void DefineGeneralHistograms();
  void DefineAnalysisHistograms();
  Bool_t HaveCharmInHistory(Int_t *history);
  Bool_t HaveBottomInHistory(Int_t *history);
  Int_t FromSemileptonicDecays(Int_t *history);

  AliAODVertex *CallPrimaryVertex(AliAODv0 *v0, AliAODTrack *trk, AliAODEvent *evt);
  AliAODVertex* PrimaryVertex(const TObjArray *trkArray,AliVEvent *event);
  TProfile* GetEstimatorHistogram(const AliVEvent *event);

  Bool_t fUseMCInfo;                 
  TList *fOutput;                    
  TList *fOutputAll;                 
  TList *fListCuts;                  
  TH1F *fCEvents;                    
  TH1F *fHTrigger;                   
  TH1F *fHCentrality;                
  TH1F *fHEventPlane;                
  TH2F *fHNTrackletvsZ;                
  TH2F *fHNTrackletCorrvsZ;                
  AliRDHFCutsLctoeleLambdafromAODtracks *fAnalCuts;// Cuts - sent to output slot 2
  Bool_t fIsEventSelected;          
  Bool_t    fWriteVariableTree;     
  Bool_t    fWriteEachVariableTree;     
  Bool_t    fWriteMCVariableTree;     
  TTree    *fVariablesTree;         
  TTree    *fEleVariablesTree;         
  TTree    *fV0VariablesTree;         
  TTree    *fSingleVariablesTree;         
  TTree    *fMCVariablesTree;         
  TTree    *fMCEleVariablesTree;         
  TTree    *fMCV0VariablesTree;         
  TTree    *fMCGenPairVariablesTree;         
  TTree* fCorrelationVariablesTree;         
  Bool_t fReconstructPrimVert;       
  Bool_t fIsMB;       
  Bool_t fIsSemi;     
  Bool_t fIsCent;     
  Bool_t fIsINT7;     
  Bool_t fIsEMC7;     
  Float_t *fCandidateVariables;   
  Float_t *fCandidateEleVariables;   
  Float_t *fCandidateV0Variables;   
  Float_t *fCandidateSingleVariables;   
  Float_t *fCandidateMCVariables;   
  Float_t *fCandidateMCEleVariables;   
  Float_t *fCandidateMCV0Variables;   
  Float_t *fCandidateMCGenPairVariables;   
  Float_t *fCorrelationVariables;   
  AliAODVertex *fVtx1;            
  AliESDVertex *fV1;              
  Float_t  fVtxZ;         
  Float_t  fEventPlane;         
  TVector2 *fQ;         
  TVector2 *fQSub1;         
  TVector2 *fQSub2;         
  Double_t fBzkG;                 
  Float_t  fCentrality;           
  Int_t  fRunNumber;           
  Float_t  fTriggerCheck;         
  Bool_t  fUseCentralityV0M;         
  Bool_t  fUseCentralitySPDTracklet;         
  Int_t  fUseEventPlane;         
  Int_t  fEvNumberCounter;         
  Int_t fMCEventType; 
  Bool_t fMCDoPairAnalysis; 

  //--------------------- My histograms ------------------
  THnSparse* fHistoEleLambdaMass;         
  THnSparse* fHistoEleLambdaMassRS;         
  THnSparse* fHistoEleLambdaMassRS1;         
  THnSparse* fHistoEleLambdaMassRS2;         
  THnSparse* fHistoEleLambdaMassWS;         
  THnSparse* fHistoEleLambdaMassWS1;         
  THnSparse* fHistoEleLambdaMassWS2;         
  THnSparse* fHistoEleLambdaMassRSMix;         
  THnSparse* fHistoEleLambdaMassRSMix1;         
  THnSparse* fHistoEleLambdaMassRSMix2;         
  THnSparse* fHistoEleLambdaMassWSMix;         
  THnSparse* fHistoEleLambdaMassWSMix1;         
  THnSparse* fHistoEleLambdaMassWSMix2;         
  THnSparse* fHistoEleLambdaMassRSSide;         
  THnSparse* fHistoEleLambdaMassRSSide1;         
  THnSparse* fHistoEleLambdaMassRSSide2;         
  THnSparse* fHistoEleLambdaMassWSSide;         
  THnSparse* fHistoEleLambdaMassWSSide1;         
  THnSparse* fHistoEleLambdaMassWSSide2;         

  THnSparse* fHistoEleLambdaMassRSAway;         
  THnSparse* fHistoEleLambdaMassRS1Away;         
  THnSparse* fHistoEleLambdaMassRS2Away;         
  THnSparse* fHistoEleLambdaMassWSAway;         
  THnSparse* fHistoEleLambdaMassWS1Away;         
  THnSparse* fHistoEleLambdaMassWS2Away;         
  THnSparse* fHistoEleLambdaMassRSMixAway;         
  THnSparse* fHistoEleLambdaMassRSMix1Away;         
  THnSparse* fHistoEleLambdaMassRSMix2Away;         
  THnSparse* fHistoEleLambdaMassWSMixAway;         
  THnSparse* fHistoEleLambdaMassWSMix1Away;         
  THnSparse* fHistoEleLambdaMassWSMix2Away;         
  THnSparse* fHistoEleLambdaMassRSSideAway;         
  THnSparse* fHistoEleLambdaMassRSSide1Away;         
  THnSparse* fHistoEleLambdaMassRSSide2Away;         
  THnSparse* fHistoEleLambdaMassWSSideAway;         
  THnSparse* fHistoEleLambdaMassWSSide1Away;         
  THnSparse* fHistoEleLambdaMassWSSide2Away;         

  THnSparse* fHistoEleLambdaMassvsElePtRS;         
  THnSparse* fHistoEleLambdaMassvsElePtWS;         
  THnSparse* fHistoEleLambdaMassvsElePtRSMix;         
  THnSparse* fHistoEleLambdaMassvsElePtWSMix;         
  THnSparse* fHistoEleLambdaMassvsElePtRSSide;         
  THnSparse* fHistoEleLambdaMassvsElePtWSSide;         
  THnSparse* fHistoEleLambdaMassvsElePtRS1;         
  THnSparse* fHistoEleLambdaMassvsElePtWS1;         
  THnSparse* fHistoEleLambdaMassvsElePtRSMix1;         
  THnSparse* fHistoEleLambdaMassvsElePtWSMix1;         
  THnSparse* fHistoEleLambdaMassvsElePtRSSide1;         
  THnSparse* fHistoEleLambdaMassvsElePtWSSide1;         
  THnSparse* fHistoEleLambdaMassvsElePtRS2;         
  THnSparse* fHistoEleLambdaMassvsElePtWS2;         
  THnSparse* fHistoEleLambdaMassvsElePtRSMix2;         
  THnSparse* fHistoEleLambdaMassvsElePtWSMix2;         
  THnSparse* fHistoEleLambdaMassvsElePtRSSide2;         
  THnSparse* fHistoEleLambdaMassvsElePtWSSide2;         
  TH2F* fHistoElePtRS;         
  TH2F* fHistoElePtWS;         
  TH2F* fHistoElePtRSMix;         
  TH2F* fHistoElePtWSMix;         

  //Efficiency calculation
  THnSparse* fHistoEleLambdaMassMCS;         
  THnSparse* fHistoEleLambdaMassLambdabMCS;         
  THnSparse* fHistoEleLambdaMassXibMCS;         
  THnSparse* fHistoEleLambdaMassMCGen;         
  THnSparse* fHistoEleLambdaMassvsElePtMCS;         
  THnSparse* fHistoEleLambdaMassvsElePtMCGen;         
  THnSparse* fHistoEleLambdaMassvsElePtMCS1;         
  THnSparse* fHistoEleLambdaMassvsElePtMCGen1;         
  THnSparse* fHistoEleLambdaMassvsElePtMCS2;         
  THnSparse* fHistoEleLambdaMassvsElePtMCGen2;         
  TH2F* fHistoElePtMCS;         
  TH2F* fHistoElePtMCGen;         

  THnSparse* fHistoElePtvsEtaRS;         
  THnSparse* fHistoElePtvsEtaWS;         
  THnSparse* fHistoElePtvsEtaRSMix;         
  THnSparse* fHistoElePtvsEtaWSMix;         
  THnSparse* fHistoElePtvsEtaMCS;         
  THnSparse* fHistoElePtvsEtaMCGen;         

  THnSparse* fHistoElePtvsLambdaPtRS;         
  THnSparse* fHistoElePtvsLambdaPtWS;         
  THnSparse* fHistoElePtvsLambdaPtRSMix;         
  THnSparse* fHistoElePtvsLambdaPtWSMix;         
  THnSparse* fHistoElePtvsLambdaPtMCS;         
  THnSparse* fHistoElePtvsLambdaPtvsLcPtMCS;         
  THnSparse* fHistoElePtvsLambdaPtMCGen;         
  THnSparse* fHistoElePtvsLambdaPtvsLcPtMCGen;         
  THnSparse* fHistoElePtvsLambdaPtMCLcGen;         

  THnSparse* fHistoElePtvsd0RS;         
  THnSparse* fHistoElePtvsd0WS;         
  THnSparse* fHistoElePtvsd0RSMix;         
  THnSparse* fHistoElePtvsd0WSMix;         
  THnSparse* fHistoElePtvsd0MCS;         
  THnSparse* fHistoElePtvsd0PromptMCS;         
  THnSparse* fHistoElePtvsd0BFeeddownMCS;         

  //Feeddown from Xic0
  THnSparse* fHistoEleLambdaMassFeeddownXic0MCS;         
  THnSparse* fHistoEleLambdaMassFeeddownXic0MCGen;         
  THnSparse* fHistoEleLambdaMassvsElePtFeeddownXic0MCS;         
  THnSparse* fHistoEleLambdaMassvsElePtFeeddownXic0MCGen;         
  THnSparse* fHistoEleLambdaMassvsElePtFeeddownXic0MCS1;         
  THnSparse* fHistoEleLambdaMassvsElePtFeeddownXic0MCGen1;         
  THnSparse* fHistoEleLambdaMassvsElePtFeeddownXic0MCS2;         
  THnSparse* fHistoEleLambdaMassvsElePtFeeddownXic0MCGen2;         
  TH2F* fHistoElePtFeeddownXic0MCS;         
  TH2F* fHistoElePtFeeddownXic0MCGen;         
  THnSparse* fHistoElePtvsEtaFeeddownXic0MCS;         
  THnSparse* fHistoElePtvsEtaFeeddownXic0MCGen;         
  THnSparse* fHistoElePtvsLambdaPtFeeddownXic0MCS;         
  THnSparse* fHistoElePtvsLambdaPtFeeddownXic0MCGen;         

  //Feeddown from XicPlus
  THnSparse* fHistoEleLambdaMassFeeddownXicPlusMCS;         
  THnSparse* fHistoEleLambdaMassFeeddownXicPlusMCGen;         
  THnSparse* fHistoEleLambdaMassvsElePtFeeddownXicPlusMCS;         
  THnSparse* fHistoEleLambdaMassvsElePtFeeddownXicPlusMCGen;         
  THnSparse* fHistoEleLambdaMassvsElePtFeeddownXicPlusMCS1;         
  THnSparse* fHistoEleLambdaMassvsElePtFeeddownXicPlusMCGen1;         
  THnSparse* fHistoEleLambdaMassvsElePtFeeddownXicPlusMCS2;         
  THnSparse* fHistoEleLambdaMassvsElePtFeeddownXicPlusMCGen2;         
  TH2F* fHistoElePtFeeddownXicPlusMCS;         
  TH2F* fHistoElePtFeeddownXicPlusMCGen;         
  THnSparse* fHistoElePtvsEtaFeeddownXicPlusMCS;         
  THnSparse* fHistoElePtvsEtaFeeddownXicPlusMCGen;         
  THnSparse* fHistoElePtvsLambdaPtFeeddownXicPlusMCS;         
  THnSparse* fHistoElePtvsLambdaPtFeeddownXicPlusMCGen;         


  //Checking histograms
  TH1F* fHistoBachPt;      
  TH1F* fHistoBachPtMCS;      
  TH1F* fHistoBachPtMCGen;      
  TH1F* fHistoBachFracSharedITS;      
  TH1F* fHistoBachChi2NclsITS;      
  TH1F* fHistod0Bach;      
  TH2F* fHistoLambdaMassvsPt;     
  TH2F* fHistoLambdaMassvsPtMCS;     
  TH2F* fHistoLambdaMassvsPtMCGen;     
  TH2F* fHistoLambdaPtvsEta;     
  TH2F* fHistoLambdaPtvsEtaMCS;     
  TH2F* fHistoLambdaPtvsEtaMCGen;     
  THnSparse* fHistoLambdaArmenteros;     
  THnSparse* fHistoLambdaArmenterosMCS;     
  TH2F* fHistoLambdaPtvsDl;     
  TH2F* fHistoLambdaPtvsDlSide;     
  TH2F* fHistoLambdaPtvsDlMCS;     
  TH2F* fHistoLambdaPtvsDlFeeddownXi0MCS;    
  TH2F* fHistoLambdaPtvsDlFeeddownXiMinusMCS;    
  TH2F* fHistoLambdaPtvsDlFeeddownOmegaMCS;    
  TH2F* fHistoK0sMassvsPt;     
  TH2F* fHistoElectronTPCPID;     
  TH2F* fHistoElectronTOFPID;     
  TH2F* fHistoElectronITSPID;     
  TH2F* fHistoElectronTPCSelPID;     
  TH2F* fHistoElectronTOFSelPID;     
  TH2F* fHistoElectronITSSelPID;     
  TH2F* fHistoElectronTPCPIDSelTOF;     
  TH2F* fHistoElectronTPCPIDSelTOFSmallEta;     
  TH2F* fHistoElectronTPCPIDSelTOFLargeEta;     
  TH2F* fHistoElectronTPCPIDSelTOFEtaDep[8];     
  TH2F* fHistoElectronITSPIDSelTPCTOF;     
  TH2F* fHistoElectronITSPIDSelTPCTOFSmallEta;     
  TH2F* fHistoElectronITSPIDSelTPCTOFLargeEta;     
  TH2F* fHistoElectronITSPIDSelTPCTOFEtaDep[8];     
  TH1F* fHistoMassConversionsMin; 
  TH1F* fHistoMassConversionsSameSignMin; 
  THnSparse* fHistoElePtvsCutVarsRS[17];         
  THnSparse* fHistoElePtvsCutVarsWS[17];         
  THnSparse* fHistoElePtvsCutVarsMCS[17];         
  TH2F* fHistoElectronQovPtvsPhi;     
  TH2F* fHistoLambdaQovPtvsPhi;     
  THnSparse* fHistoLcMCGen;         
  THnSparse* fHistoLcMCGen1;         
  THnSparse* fHistoLcMCGen2;         
  THnSparse* fHistoLbMCGen;         
  THnSparse* fHistoLbMCGen1;         
  THnSparse* fHistoLbMCGen2;         
  THnSparse* fHistoLcMCGenWithEleA;         
  TH2D* fHistoLcMCGenWithEleB;         
  THnSparse* fHistoLbMCGenWithEleA;         
  TH2D* fHistoLbMCGenWithEleB;         
  THnSparse* fHistoLbMCGenWithLambdac;         
  THnSparse* fHistoLcMCS;         
  THnSparse* fHistoLcMCS1;         
  THnSparse* fHistoLcMCS2;         
  THnSparse* fHistoLbMCS;         
  THnSparse* fHistoLbMCS1;         
  THnSparse* fHistoLbMCS2;         
  THnSparse* fHistoFeedDownXic0MCGen;     
  THnSparse* fHistoFeedDownXic0MCGen1;     
  THnSparse* fHistoFeedDownXic0MCGen2;     
  THnSparse* fHistoFeedDownXic0MCS;     
  THnSparse* fHistoFeedDownXic0MCS1;     
  THnSparse* fHistoFeedDownXic0MCS2;     
  THnSparse* fHistoFeedDownXicPlusMCGen;     
  THnSparse* fHistoFeedDownXicPlusMCGen1;     
  THnSparse* fHistoFeedDownXicPlusMCGen2;     
  THnSparse* fHistoFeedDownXicPlusMCS;     
  THnSparse* fHistoFeedDownXicPlusMCS1;     
  THnSparse* fHistoFeedDownXicPlusMCS2;     
  THnSparse* fHistoFeedDownSigmaMCGen;     
  THnSparse* fHistoFeedDownSigmaMCGen1;     
  THnSparse* fHistoFeedDownSigmaMCGen2;     
  THnSparse* fHistoFeedDownSigmaMCS;     
  THnSparse* fHistoFeedDownSigmaMCS1;     
  THnSparse* fHistoFeedDownSigmaMCS2;     
  THnSparse* fHistoLcElectronMCGen;         
  THnSparse* fHistoLcElectronMCGen1;         
  THnSparse* fHistoLcElectronMCGen2;         
  THnSparse* fHistoLcElectronMCS;         
  THnSparse* fHistoLcElectronMCS1;         
  THnSparse* fHistoLcElectronMCS2;         
  THnSparse* fHistoLbElectronMCGen;         
  THnSparse* fHistoLbElectronMCGen1;         
  THnSparse* fHistoLbElectronMCGen2;         
  THnSparse* fHistoLbElectronMCS;         
  THnSparse* fHistoLbElectronMCS1;         
  THnSparse* fHistoLbElectronMCS2;         
  THnSparse* fHistoElectronFeedDownXic0MCGen;     
  THnSparse* fHistoElectronFeedDownXic0MCGen1;     
  THnSparse* fHistoElectronFeedDownXic0MCGen2;     
  THnSparse* fHistoElectronFeedDownXic0MCS;     
  THnSparse* fHistoElectronFeedDownXic0MCS1;     
  THnSparse* fHistoElectronFeedDownXic0MCS2;     
  THnSparse* fHistoElectronFeedDownXicPlusMCGen;     
  THnSparse* fHistoElectronFeedDownXicPlusMCGen1;     
  THnSparse* fHistoElectronFeedDownXicPlusMCGen2;     
  THnSparse* fHistoElectronFeedDownXicPlusMCS;     
  THnSparse* fHistoElectronFeedDownXicPlusMCS1;     
  THnSparse* fHistoElectronFeedDownXicPlusMCS2;     
  THnSparse* fHistoElectronFeedDownXibMCGen;     
  THnSparse* fHistoElectronFeedDownXibMCGen1;     
  THnSparse* fHistoElectronFeedDownXibMCGen2;     
  THnSparse* fHistoElectronFeedDownXibMCS;     
  THnSparse* fHistoElectronFeedDownXibMCS1;     
  THnSparse* fHistoElectronFeedDownXibMCS2;     
  THnSparse* fHistoElectronMCGen;         
  THnSparse* fHistoBottomElectronMCGen;         
  THnSparse* fHistoCharmElectronMCGen;         
  THnSparse* fHistoLambdaMCGen;         

  THnSparse* fHistoElePtvsV0dlRS;         
  THnSparse* fHistoElePtvsV0dlRS1;         
  THnSparse* fHistoElePtvsV0dlRS2;         
  THnSparse* fHistoElePtvsV0dlRSSide;         
  THnSparse* fHistoElePtvsV0dlRSSide1;         
  THnSparse* fHistoElePtvsV0dlRSSide2;         
  THnSparse* fHistoElePtvsV0dlRSMix;         
  THnSparse* fHistoElePtvsV0dlRSMix1;         
  THnSparse* fHistoElePtvsV0dlRSMix2;         
  THnSparse* fHistoElePtvsV0dlWS;         
  THnSparse* fHistoElePtvsV0dlWS1;         
  THnSparse* fHistoElePtvsV0dlWS2;         
  THnSparse* fHistoElePtvsV0dlWSSide;         
  THnSparse* fHistoElePtvsV0dlWSSide1;         
  THnSparse* fHistoElePtvsV0dlWSSide2;         
  THnSparse* fHistoElePtvsV0dlWSMix;         
  THnSparse* fHistoElePtvsV0dlWSMix1;         
  THnSparse* fHistoElePtvsV0dlWSMix2;         
  THnSparse* fHistoElePtvsV0dlMCS;         
  THnSparse* fHistoElePtvsV0dlMCS1;         
  THnSparse* fHistoElePtvsV0dlMCS2;         
  THnSparse* fHistoElePtvsV0dlFeedDownXic0MCS;         
  THnSparse* fHistoElePtvsV0dlFeedDownXic0MCS1;         
  THnSparse* fHistoElePtvsV0dlFeedDownXic0MCS2;         
  THnSparse* fHistoElePtvsV0dlFeedDownXicPlusMCS;         
  THnSparse* fHistoElePtvsV0dlFeedDownXicPlusMCS1;         
  THnSparse* fHistoElePtvsV0dlFeedDownXicPlusMCS2;         

  THnSparse* fHistoElePtvsV0dcaRS;         
  THnSparse* fHistoElePtvsV0dcaRS1;         
  THnSparse* fHistoElePtvsV0dcaRS2;         
  THnSparse* fHistoElePtvsV0dcaRSSide;         
  THnSparse* fHistoElePtvsV0dcaRSSide1;         
  THnSparse* fHistoElePtvsV0dcaRSSide2;         
  THnSparse* fHistoElePtvsV0dcaRSMix;         
  THnSparse* fHistoElePtvsV0dcaRSMix1;         
  THnSparse* fHistoElePtvsV0dcaRSMix2;         
  THnSparse* fHistoElePtvsV0dcaWS;         
  THnSparse* fHistoElePtvsV0dcaWS1;         
  THnSparse* fHistoElePtvsV0dcaWS2;         
  THnSparse* fHistoElePtvsV0dcaWSSide;         
  THnSparse* fHistoElePtvsV0dcaWSSide1;         
  THnSparse* fHistoElePtvsV0dcaWSSide2;         
  THnSparse* fHistoElePtvsV0dcaWSMix;         
  THnSparse* fHistoElePtvsV0dcaWSMix1;         
  THnSparse* fHistoElePtvsV0dcaWSMix2;         
  THnSparse* fHistoElePtvsV0dcaMCS;         
  THnSparse* fHistoElePtvsV0dcaMCS1;         
  THnSparse* fHistoElePtvsV0dcaMCS2;         
  THnSparse* fHistoElePtvsV0dcaFeedDownXic0MCS;         
  THnSparse* fHistoElePtvsV0dcaFeedDownXic0MCS1;         
  THnSparse* fHistoElePtvsV0dcaFeedDownXic0MCS2;         
  THnSparse* fHistoElePtvsV0dcaFeedDownXicPlusMCS;         
  THnSparse* fHistoElePtvsV0dcaFeedDownXicPlusMCS1;         
  THnSparse* fHistoElePtvsV0dcaFeedDownXicPlusMCS2;         

  THnSparse* fHistoEleLambdaPtvsV0dlRS;         
  THnSparse* fHistoEleLambdaPtvsV0dlRS1;         
  THnSparse* fHistoEleLambdaPtvsV0dlRS2;         
  THnSparse* fHistoEleLambdaPtvsV0dlRSSide;         
  THnSparse* fHistoEleLambdaPtvsV0dlRSSide1;         
  THnSparse* fHistoEleLambdaPtvsV0dlRSSide2;         
  THnSparse* fHistoEleLambdaPtvsV0dlRSMix;         
  THnSparse* fHistoEleLambdaPtvsV0dlRSMix1;         
  THnSparse* fHistoEleLambdaPtvsV0dlRSMix2;         
  THnSparse* fHistoEleLambdaPtvsV0dlWS;         
  THnSparse* fHistoEleLambdaPtvsV0dlWS1;         
  THnSparse* fHistoEleLambdaPtvsV0dlWS2;         
  THnSparse* fHistoEleLambdaPtvsV0dlWSSide;         
  THnSparse* fHistoEleLambdaPtvsV0dlWSSide1;         
  THnSparse* fHistoEleLambdaPtvsV0dlWSSide2;         
  THnSparse* fHistoEleLambdaPtvsV0dlWSMix;         
  THnSparse* fHistoEleLambdaPtvsV0dlWSMix1;         
  THnSparse* fHistoEleLambdaPtvsV0dlWSMix2;         
  THnSparse* fHistoEleLambdaPtvsV0dlMCS;         
  THnSparse* fHistoEleLambdaPtvsV0dlMCS1;         
  THnSparse* fHistoEleLambdaPtvsV0dlMCS2;         
  THnSparse* fHistoEleLambdaPtvsV0dlPromptMCS;         
  THnSparse* fHistoEleLambdaPtvsV0dlPromptMCS1;         
  THnSparse* fHistoEleLambdaPtvsV0dlPromptMCS2;         
  THnSparse* fHistoEleLambdaPtvsV0dlBFeeddownMCS;         
  THnSparse* fHistoEleLambdaPtvsV0dlBFeeddownMCS1;         
  THnSparse* fHistoEleLambdaPtvsV0dlBFeeddownMCS2;         
  THnSparse* fHistoEleLambdaPtvsV0dlFeedDownXic0MCS;         
  THnSparse* fHistoEleLambdaPtvsV0dlFeedDownXic0MCS1;         
  THnSparse* fHistoEleLambdaPtvsV0dlFeedDownXic0MCS2;         
  THnSparse* fHistoEleLambdaPtvsV0dlFeedDownXic0PromptMCS;         
  THnSparse* fHistoEleLambdaPtvsV0dlFeedDownXic0PromptMCS1;         
  THnSparse* fHistoEleLambdaPtvsV0dlFeedDownXic0PromptMCS2;         
  THnSparse* fHistoEleLambdaPtvsV0dlFeedDownXic0BFeeddownMCS;         
  THnSparse* fHistoEleLambdaPtvsV0dlFeedDownXic0BFeeddownMCS1;         
  THnSparse* fHistoEleLambdaPtvsV0dlFeedDownXic0BFeeddownMCS2;         
  THnSparse* fHistoEleLambdaPtvsV0dlFeedDownXicPlusMCS;         
  THnSparse* fHistoEleLambdaPtvsV0dlFeedDownXicPlusMCS1;         
  THnSparse* fHistoEleLambdaPtvsV0dlFeedDownXicPlusMCS2;         
  THnSparse* fHistoEleLambdaPtvsV0dlFeedDownXicPlusPromptMCS;         
  THnSparse* fHistoEleLambdaPtvsV0dlFeedDownXicPlusPromptMCS1;         
  THnSparse* fHistoEleLambdaPtvsV0dlFeedDownXicPlusPromptMCS2;         
  THnSparse* fHistoEleLambdaPtvsV0dlFeedDownXicPlusBFeeddownMCS;         
  THnSparse* fHistoEleLambdaPtvsV0dlFeedDownXicPlusBFeeddownMCS1;         
  THnSparse* fHistoEleLambdaPtvsV0dlFeedDownXicPlusBFeeddownMCS2;         
  THnSparse* fHistoEleLambdaPtvsV0dlFeedDownSigmaMCS;         
  THnSparse* fHistoEleLambdaPtvsV0dlFeedDownSigmaMCS1;         
  THnSparse* fHistoEleLambdaPtvsV0dlFeedDownSigmaMCS2;         
  THnSparse* fHistoEleLambdaPtvsV0dlFeedDownSigmaPromptMCS;         
  THnSparse* fHistoEleLambdaPtvsV0dlFeedDownSigmaPromptMCS1;         
  THnSparse* fHistoEleLambdaPtvsV0dlFeedDownSigmaPromptMCS2;         
  THnSparse* fHistoEleLambdaPtvsV0dlFeedDownSigmaBFeeddownMCS;         
  THnSparse* fHistoEleLambdaPtvsV0dlFeedDownSigmaBFeeddownMCS1;         
  THnSparse* fHistoEleLambdaPtvsV0dlFeedDownSigmaBFeeddownMCS2;         
  THnSparse* fHistoEleLambdaPtvsV0dlLambdabMCS;         
  THnSparse* fHistoEleLambdaPtvsV0dlLambdabMCS1;         
  THnSparse* fHistoEleLambdaPtvsV0dlLambdabMCS2;         
  THnSparse* fHistoEleLambdaPtvsV0dlXibMCS;         
  THnSparse* fHistoEleLambdaPtvsV0dlXibMCS1;         
  THnSparse* fHistoEleLambdaPtvsV0dlXibMCS2;         

  THnSparse* fHistoEleLambdaPtvsV0dlRSAway;         
  THnSparse* fHistoEleLambdaPtvsV0dlRS1Away;         
  THnSparse* fHistoEleLambdaPtvsV0dlRS2Away;         
  THnSparse* fHistoEleLambdaPtvsV0dlRSSideAway;         
  THnSparse* fHistoEleLambdaPtvsV0dlRSSide1Away;         
  THnSparse* fHistoEleLambdaPtvsV0dlRSSide2Away;         
  THnSparse* fHistoEleLambdaPtvsV0dlRSMixAway;         
  THnSparse* fHistoEleLambdaPtvsV0dlRSMix1Away;         
  THnSparse* fHistoEleLambdaPtvsV0dlRSMix2Away;         
  THnSparse* fHistoEleLambdaPtvsV0dlWSAway;         
  THnSparse* fHistoEleLambdaPtvsV0dlWS1Away;         
  THnSparse* fHistoEleLambdaPtvsV0dlWS2Away;         
  THnSparse* fHistoEleLambdaPtvsV0dlWSSideAway;         
  THnSparse* fHistoEleLambdaPtvsV0dlWSSide1Away;         
  THnSparse* fHistoEleLambdaPtvsV0dlWSSide2Away;         
  THnSparse* fHistoEleLambdaPtvsV0dlWSMixAway;         
  THnSparse* fHistoEleLambdaPtvsV0dlWSMix1Away;         
  THnSparse* fHistoEleLambdaPtvsV0dlWSMix2Away;         

  THnSparse* fHistoCorrelationVariablesvsEleLambdaPt;         
  THnSparse* fHistoCorrelationVariablesvsEleLambdaPtMix;         
  THnSparse* fHistoCorrelationVariablesvsEleLambdaPtMC;         
  THnSparse* fHistoCorrelationVariablesvsElePt;         
  THnSparse* fHistoCorrelationVariablesvsElePtMix;         
  THnSparse* fHistoCorrelationVariablesvsElePtMC;         
  THnSparse* fHistoCorrelationVariablesvsLambdaPt;         
  THnSparse* fHistoCorrelationVariablesvsLambdaPtMix;         
  THnSparse* fHistoCorrelationVariablesvsLambdaPtMC;         

  THnSparse* fHistoMassVariablesvsEleLambdaPt;         
  THnSparse* fHistoMassVariablesvsEleLambdaPtMix;         
  THnSparse* fHistoMassVariablesvsEleLambdaPtMC;         
  THnSparse* fHistoMassVariablesvsElePt;         
  THnSparse* fHistoMassVariablesvsElePtMix;         
  THnSparse* fHistoMassVariablesvsElePtMC;         
  THnSparse* fHistoMassVariablesvsLambdaPt;         
  THnSparse* fHistoMassVariablesvsLambdaPtMix;         
  THnSparse* fHistoMassVariablesvsLambdaPtMC;         

  TH2D *fHistoResponseElePt; 
  TH2D *fHistoResponseElePt1; 
  TH2D *fHistoResponseElePt2; 
  TH2D *fHistoResponseEleLambdaPt; 
  TH2D *fHistoResponseEleLambdaPt1; 
  TH2D *fHistoResponseEleLambdaPt2; 
  TH2D *fHistoResponseEleLambdaPtFeeddownXic0; 
  TH2D *fHistoResponseEleLambdaPtFeeddownXic01; 
  TH2D *fHistoResponseEleLambdaPtFeeddownXic02; 
  TH2D *fHistoResponseEleLambdaPtFeeddownXicPlus; 
  TH2D *fHistoResponseEleLambdaPtFeeddownXicPlus1; 
  TH2D *fHistoResponseEleLambdaPtFeeddownXicPlus2; 
  TH2D *fHistoResponseEleLambdaPtFeeddownSigma; 
  TH2D *fHistoResponseEleLambdaPtFeeddownSigma1; 
  TH2D *fHistoResponseEleLambdaPtFeeddownSigma2; 
  TH2D *fHistoResponseEleLambdaPtLambdab; 
  TH2D *fHistoResponseElePtLambdab; 
  TH2D *fHistoResponseLambdaPtLambdab; 
  TH2D *fHistoResponseLambdaPt; 
  TH2D *fHistoResponseLambdaPtFeeddownXic0; 
  TH2D *fHistoResponseLambdaPtFeeddownXicPlus; 
  TH2D *fHistoResponseLambdaPtFeeddownSigma; 
  TH2D *fHistoResponseLambdaPtvsEleLambdaPt; 
  TH2D *fHistoResponseLambdaPtvsEleLambdaPtLambdab; 
  TH2D *fHistoResponseMcGenLambdabPtvsLambdacPt; 
  THnSparse* fHistoLcPtvseleLambdaPtvsElePtvsLambdaPt;         

  TH2F *fHistoEleLambdaPtvsRapidityRS; 
  TH2F *fHistoEleLambdaPtvsRapidityWS; 
  TH2F *fHistoEleLambdaPtvsRapidityMCS; 

  TH1F* fHistoPi0MCGen;         
  THnSparse* fHistoElectronPi0Total;         
  THnSparse* fHistoElectronPi0Tag;         
  TH1F* fHistoEtaMCGen;         
  THnSparse* fHistoElectronEtaTotal;         
  THnSparse* fHistoElectronEtaTag;         
  TH1F* fHistoKaonMCGen;         
  TH1F* fHistoD0MCGen;         

  AliNormalizationCounter *fCounter;
  TH1F *fHistonEvtvsRunNumber;
  TH1F *fHistonElevsRunNumber;
  TH1F *fHistonLambdavsRunNumber;
  TH1F *fHistoMCEventType;
  TH1F *fHistoMCDeltaPhiccbar;

  //Multiplicity dep analysis
  TProfile* fMultEstimatorAvg[4]; 
  Double_t fRefMult;   

  // Store pointers to global tracks for pid and dca
  AliAODTrack **fGTI;                
  Int_t *fGTIndex;                
  const UShort_t  fTrackBuffSize;          
  TH2D *fHistodPhiSdEtaSElectronProtonR125RS;
  TH2D *fHistodPhiSdEtaSElectronProtonR125WS;
  TH2D *fHistodPhiSdEtaSElectronProtonR125RSMix;
  TH2D *fHistodPhiSdEtaSElectronProtonR125WSMix;
  TH2D *fHistodPhiSdEtaSElectronPionR125RS;
  TH2D *fHistodPhiSdEtaSElectronPionR125WS;
  TH2D *fHistodPhiSdEtaSElectronPionR125RSMix;
  TH2D *fHistodPhiSdEtaSElectronPionR125WSMix;

  //Mixing
  Int_t fDoEventMixing; 
  Bool_t fMixWithoutConversionFlag; 
  Int_t  fNumberOfEventsForMixing; 
  Int_t fNzVtxBins;               
  Double_t fZvtxBins[100];            // [fNzVtxBinsDim]
  Int_t fNCentBins;               
  Double_t fCentBins[100];            // [fNCentBinsDim]
  Int_t fNRPBins;               
  Double_t fRPBins[100];            // [fRPBinsDim]
  Int_t  fNOfPools; 
  Int_t fPoolIndex; 
  TH2F *fHistoPoolMonitor;
  TH1F *fHistoPoolIDCounter;
  std::vector<Int_t> nextResVec; 
  std::vector<Bool_t> reservoirsReady; 
  std::vector<std::vector< std::vector< TLorentzVector * > > > m_ReservoirE; 
  std::vector<std::vector< std::vector< TLorentzVector * > > > m_ReservoirL1; 
  std::vector<std::vector< std::vector< TLorentzVector * > > > m_ReservoirL2; 
  std::vector<std::vector< std::vector< TVector * > > > m_ReservoirVarsE; 
  std::vector<std::vector< std::vector< TVector * > > > m_ReservoirVarsL1; 
  std::vector<std::vector< std::vector< TVector * > > > m_ReservoirVarsL2; 

  ClassDef(AliAnalysisTaskSELc2eleLambdafromAODtracks,42); 
};
#endif