AliAnalysisTaskSEXic2eleXifromAODtracks.h

Go to the documentation of this file.

#ifndef ALIANALYSISTASKSEXIC2ELEXIFROMAODTRACKS_H
#define ALIANALYSISTASKSEXIC2ELEXIFROMAODTRACKS_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 "AliRDHFCutsXictoeleXifromAODtracks.h"


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

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

  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, AliAODcascade *casc, AliAODTrack *trk, AliAODTrack *trkpid, AliAODEvent *event,  TClonesArray *mcArray);
  void FillMixROOTObjects(TLorentzVector *et, TLorentzVector *ev, TVector *tinfo, TVector *vinfo, Int_t charge);
  void FillElectronROOTObjects(AliAODTrack *trk, AliAODEvent *event, TClonesArray *mcArray);
  void FillCascROOTObjects(AliAODcascade *casc, AliAODEvent *event, TClonesArray *mcArray);
  void FillMCROOTObjects(AliAODMCParticle *part, AliAODMCParticle *mcepart, AliAODMCParticle *mcv0part, Int_t decaytype);
  void FillMCEleROOTObjects(AliAODMCParticle *mcepart, TClonesArray *mcArray);
  void FillMCCascROOTObjects(AliAODMCParticle *mccpart, TClonesArray *mcArray);
  void FillMCGenPairROOTObjects(AliAODMCParticle *mcparte, AliAODMCParticle* mcpartv, TClonesArray *mcArray);
  Bool_t MakeMCAnalysis(TClonesArray *mcArray);
  void MakeAnalysis(AliAODEvent *aod, TClonesArray *mcArray);

  void SelectCascade( const AliVEvent *event,Int_t nCasc,Int_t &nSeleCasc, Bool_t *seleCascFlags, 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, AliAODcascade *v0, AliEventplane *pl,TVector2* q,TVector2* qsub1,TVector2* qsub2);

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

  AliAODRecoCascadeHF* MakeCascadeHF(AliAODcascade *casc, AliAODTrack *trk, AliAODTrack *trkpid, AliAODEvent *aod, AliAODVertex *vert);
  AliAODVertex* ReconstructSecondaryVertex(AliAODcascade *casc, AliAODTrack *trk, AliAODEvent *aod);
  Int_t MatchToMC(AliAODRecoCascadeHF *elobj, TClonesArray *mcArray, Int_t *pdgarray_ele, Int_t *pdgarray_casc, Int_t *labelarray_ele, Int_t *labelarray_casc,  Int_t &ngen_ele, Int_t &ngen_casc);
  Int_t MatchToMCCascade(AliAODcascade *theCascade, Int_t pdgabscasc, Int_t *pdgDgcasc, Int_t *pdgDgv0, TClonesArray *mcArray) const;
  void  GetMCDecayHistory(AliAODMCParticle *mcpart, TClonesArray *mcArray, Int_t *pdgarray, Int_t *labelarray, Int_t &ngen);

  void StoreGlobalTrackReference(AliAODTrack *track, Int_t id);
  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 chargexi);
  Int_t GetPoolIndex(Double_t zvert, Double_t mult, Double_t rp);


 private:

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

  void DefineTreeVariables();
  void DefineEleTreeVariables();
  void DefineSingleTreeVariables();
  void DefineCascTreeVariables();
  void DefineMCTreeVariables();
  void DefineMCEleTreeVariables();
  void DefineMCCascTreeVariables();
  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(AliAODcascade *casc, 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;                
  AliRDHFCutsXictoeleXifromAODtracks *fAnalCuts;// Cuts - sent to output slot 2
  Bool_t fIsEventSelected;          
  Bool_t    fWriteVariableTree;     
  TTree    *fVariablesTree;         
  Bool_t    fWriteEachVariableTree;     
  Bool_t    fWriteMCVariableTree;     
  TTree    *fEleVariablesTree;         
  TTree    *fCascVariablesTree;         
  TTree    *fSingleVariablesTree;         
  TTree    *fMCVariablesTree;         
  TTree    *fMCEleVariablesTree;         
  TTree    *fMCCascVariablesTree;         
  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 *fCandidateCascVariables;   
  Float_t *fCandidateSingleVariables;   
  Float_t *fCandidateMCVariables;   
  Float_t *fCandidateMCEleVariables;   
  Float_t *fCandidateMCCascVariables;   
  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* fHistoEleXiMass;         
  THnSparse* fHistoEleXiMassRS;         
  THnSparse* fHistoEleXiMassWS;         
  THnSparse* fHistoEleXiMassRSMix;         
  THnSparse* fHistoEleXiMassWSMix;         
  THnSparse* fHistoEleXiMassRSSide;         
  THnSparse* fHistoEleXiMassWSSide;         
  THnSparse* fHistoEleXiMassRS1;         
  THnSparse* fHistoEleXiMassWS1;         
  THnSparse* fHistoEleXiMassRSMix1;         
  THnSparse* fHistoEleXiMassWSMix1;         
  THnSparse* fHistoEleXiMassRSSide1;         
  THnSparse* fHistoEleXiMassWSSide1;         
  THnSparse* fHistoEleXiMassRS2;         
  THnSparse* fHistoEleXiMassWS2;         
  THnSparse* fHistoEleXiMassRSMix2;         
  THnSparse* fHistoEleXiMassWSMix2;         
  THnSparse* fHistoEleXiMassRSSide2;         
  THnSparse* fHistoEleXiMassWSSide2;         
  THnSparse* fHistoEleXiMassAway;         
  THnSparse* fHistoEleXiMassRSAway;         
  THnSparse* fHistoEleXiMassWSAway;         
  THnSparse* fHistoEleXiMassRSMixAway;         
  THnSparse* fHistoEleXiMassWSMixAway;         
  THnSparse* fHistoEleXiMassRSSideAway;         
  THnSparse* fHistoEleXiMassWSSideAway;         
  THnSparse* fHistoEleXiMassRS1Away;         
  THnSparse* fHistoEleXiMassWS1Away;         
  THnSparse* fHistoEleXiMassRSMix1Away;         
  THnSparse* fHistoEleXiMassWSMix1Away;         
  THnSparse* fHistoEleXiMassRSSide1Away;         
  THnSparse* fHistoEleXiMassWSSide1Away;         
  THnSparse* fHistoEleXiMassRS2Away;         
  THnSparse* fHistoEleXiMassWS2Away;         
  THnSparse* fHistoEleXiMassRSMix2Away;         
  THnSparse* fHistoEleXiMassWSMix2Away;         
  THnSparse* fHistoEleXiMassRSSide2Away;         
  THnSparse* fHistoEleXiMassWSSide2Away;         

  THnSparse* fHistoEleXiMassvsElePtRS;         
  THnSparse* fHistoEleXiMassvsElePtWS;         
  THnSparse* fHistoEleXiMassvsElePtRSMix;         
  THnSparse* fHistoEleXiMassvsElePtWSMix;         
  THnSparse* fHistoEleXiMassvsElePtRSSide;         
  THnSparse* fHistoEleXiMassvsElePtWSSide;         
  THnSparse* fHistoEleXiMassvsElePtRS1;         
  THnSparse* fHistoEleXiMassvsElePtWS1;         
  THnSparse* fHistoEleXiMassvsElePtRSMix1;         
  THnSparse* fHistoEleXiMassvsElePtWSMix1;         
  THnSparse* fHistoEleXiMassvsElePtRSSide1;         
  THnSparse* fHistoEleXiMassvsElePtWSSide1;         
  THnSparse* fHistoEleXiMassvsElePtRS2;         
  THnSparse* fHistoEleXiMassvsElePtWS2;         
  THnSparse* fHistoEleXiMassvsElePtRSMix2;         
  THnSparse* fHistoEleXiMassvsElePtWSMix2;         
  THnSparse* fHistoEleXiMassvsElePtRSSide2;         
  THnSparse* fHistoEleXiMassvsElePtWSSide2;         
  TH2F* fHistoElePtRS;         
  TH2F* fHistoElePtWS;         
  TH2F* fHistoElePtRSMix;         
  TH2F* fHistoElePtWSMix;         
  THnSparse* fHistoEleXiMassMCS;         
  THnSparse* fHistoEleXiMassMCS1;         
  THnSparse* fHistoEleXiMassMCS2;         
  THnSparse* fHistoEleXiMassXibMCS;         
  THnSparse* fHistoEleXiMassXibMCS1;         
  THnSparse* fHistoEleXiMassXibMCS2;         
  THnSparse* fHistoEleXiMassPromptMCS;         
  THnSparse* fHistoEleXiMassPromptMCS1;         
  THnSparse* fHistoEleXiMassPromptMCS2;         
  THnSparse* fHistoEleXiMassBFeeddownMCS;         
  THnSparse* fHistoEleXiMassBFeeddownMCS1;         
  THnSparse* fHistoEleXiMassBFeeddownMCS2;         
  THnSparse* fHistoEleXiMassMCGen;         
  THnSparse* fHistoEleXiMassvsElePtMCS;         
  THnSparse* fHistoEleXiMassvsElePtMCGen;         
  THnSparse* fHistoEleXiMassvsElePtMCS1;         
  THnSparse* fHistoEleXiMassvsElePtMCGen1;         
  THnSparse* fHistoEleXiMassvsElePtMCS2;         
  THnSparse* fHistoEleXiMassvsElePtMCGen2;         
  TH2F* fHistoElePtMCS;         
  TH2F* fHistoElePtMCGen;         

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

  THnSparse* fHistoElePtvsXiPtRS;         
  THnSparse* fHistoElePtvsXiPtWS;         
  THnSparse* fHistoElePtvsXiPtRSMix;         
  THnSparse* fHistoElePtvsXiPtWSMix;         
  THnSparse* fHistoElePtvsXiPtMCS;         
  THnSparse* fHistoElePtvsXiPtvsXicPtMCS;         
  THnSparse* fHistoElePtvsXiPtMCGen;         
  THnSparse* fHistoElePtvsXiPtvsXicPtMCGen;         
  THnSparse* fHistoElePtvsXiPtMCXicGen;         

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

  TH1F* fHistoBachPt;      
  TH1F* fHistoBachPtMCS;      
  TH1F* fHistoBachPtMCGen;      
  TH1F* fHistod0Bach;      
  TH2F* fHistoXiMassvsPt;     
  TH2F* fHistoXiMassvsPtMCS;     
  TH2F* fHistoXiMassvsPtMCGen;     
  TH2F* fHistoOmegaMassvsPt;     
  TH2F* fHistoElectronTPCPID;     
  TH2F* fHistoElectronTOFPID;     
  TH2F* fHistoElectronTPCSelPID;     
  TH2F* fHistoElectronTOFSelPID;     
  TH2F* fHistoElectronTPCPIDSelTOF;     
  TH2F* fHistoElectronTOFPIDSelTPC;     
  TH2F* fHistoElectronTPCPIDSelTOFSmallEta;     
  TH2F* fHistoElectronTPCPIDSelTOFLargeEta;     
  TH2F* fHistoElectronTPCPIDSelTOFEtaDep[8];     
  TH1F* fHistoMassConversionsMin; 
  TH1F* fHistoMassConversionsSameSignMin; 
  THnSparse* fHistoElePtvsCutVarsRS[23];         
  THnSparse* fHistoElePtvsCutVarsWS[23];         
  THnSparse* fHistoElePtvsCutVarsMCS[23];         
  TH2F* fHistoElectronQovPtvsPhi;     
  TH2F* fHistoXiQovPtvsPhi;     
  THnSparse* fHistoXicMCGen;         
  THnSparse* fHistoXicMCGen1;         
  THnSparse* fHistoXicMCGen2;         
  THnSparse* fHistoXicMCS;         
  THnSparse* fHistoXicMCS1;         
  THnSparse* fHistoXicMCS2;         
  THnSparse* fHistoXibMCGen;         
  THnSparse* fHistoXibMCGenWithXic;         
  THnSparse* fHistoXibMCS;         
  THnSparse* fHistoXicElectronMCGen;         
  THnSparse* fHistoXicElectronMCGen1;         
  THnSparse* fHistoXicElectronMCGen2;         
  THnSparse* fHistoXicElectronMCS;         
  THnSparse* fHistoXicElectronMCS1;         
  THnSparse* fHistoXicElectronMCS2;         
  THnSparse* fHistoElectronMCGen;         
  THnSparse* fHistoBottomElectronMCGen;         
  THnSparse* fHistoCharmElectronMCGen;         
  THnSparse* fHistoXiMCGen;         

  TH2F *fHistoLambdaPtvsDl; 
  TH2F *fHistoLambdaPtvsDlSide; 
  TH2F *fHistoLambdaPtvsDlMCS; 
  TH2F *fHistoLambdaPtvsDR; 
  TH2F *fHistoLambdaPtvsDRSide; 
  TH2F *fHistoLambdaPtvsDRMCS; 

  TH2F *fHistoEleXiPtvsRapidityRS; 
  TH2F *fHistoEleXiPtvsRapidityWS; 
  TH2F *fHistoEleXiPtvsRapidityMCS; 

  THnSparse* fHistoCorrelationVariablesvsEleXiPt;         
  THnSparse* fHistoCorrelationVariablesvsEleXiPtMix;         
  THnSparse* fHistoCorrelationVariablesvsEleXiPtMC;         
  THnSparse* fHistoCorrelationVariablesvsElePt;         
  THnSparse* fHistoCorrelationVariablesvsElePtMix;         
  THnSparse* fHistoCorrelationVariablesvsElePtMC;         
  THnSparse* fHistoCorrelationVariablesvsXiPt;         
  THnSparse* fHistoCorrelationVariablesvsXiPtMix;         
  THnSparse* fHistoCorrelationVariablesvsXiPtMC;         

  THnSparse* fHistoMassVariablesvsEleXiPt;         
  THnSparse* fHistoMassVariablesvsEleXiPtMix;         
  THnSparse* fHistoMassVariablesvsEleXiPtMC;         
  THnSparse* fHistoMassVariablesvsElePt;         
  THnSparse* fHistoMassVariablesvsElePtMix;         
  THnSparse* fHistoMassVariablesvsElePtMC;         
  THnSparse* fHistoMassVariablesvsXiPt;         
  THnSparse* fHistoMassVariablesvsXiPtMix;         
  THnSparse* fHistoMassVariablesvsXiPtMC;         

  TH2D *fHistoResponseElePt; 
  TH2D *fHistoResponseXiPt; 
  TH2D *fHistoResponseEleXiPt; 
  TH2D *fHistoResponseXiPtvsEleXiPt; 
  TH2D *fHistoResponseXiPtXib; 
  TH2D *fHistoResponseEleXiPtXib; 
  TH2D *fHistoResponseMcGenXibPtvsXicPt; 

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


  AliNormalizationCounter *fCounter;
  TH1F *fHistonEvtvsRunNumber;
  TH1F *fHistonElevsRunNumber;
  TH1F *fHistonXivsRunNumber;
  TH1F *fHistoMCEventType;
  TH1F *fHistoMCXic0Decays;
  TH1F *fHistoMCDeltaPhiccbar;
  TH1F *fHistoMCNccbar;

  //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;
  TH2D *fHistodPhiSdEtaSElectronBachelorR125RS;
  TH2D *fHistodPhiSdEtaSElectronBachelorR125WS;
  TH2D *fHistodPhiSdEtaSElectronBachelorR125RSMix;
  TH2D *fHistodPhiSdEtaSElectronBachelorR125WSMix;

  //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; 
  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(AliAnalysisTaskSEXic2eleXifromAODtracks,34); 
};
#endif