[u/mrichter/AliRoot.git] / PWGGA / GammaConv / AliAnalysisTaskGammaConvDalitzV1.h

#ifndef ALIANALYSISTASKGAMMACONVDALITZV1_H
#define ALIANALYSISTASKGAMMACONVDALITZV1_H

/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice                               */

// Analysis task for pi0->e+e-gamma (Dalitz decay)

#include "AliAnalysisTaskSE.h"
#include "AliV0ReaderV1.h"
#include "AliKFConversionPhoton.h"
#include "AliDalitzElectronSelector.h"
#include "AliConversionMesonCuts.h"
#include "AliGammaConversionAODBGHandler.h"
#include "TProfile2D.h"

class AliESDInputHandler;
class AliMCEventHandler;
class AliESDEvent;
class AliESDtrack;
class AliESDtrackCuts;
class AliESDpidCuts;
class AliV0Reader;
class AliGammaConversionHistograms;
class AliTriggerAnalysis;

class AliAnalysisTaskGammaConvDalitzV1: public AliAnalysisTaskSE
{
  public:

	AliAnalysisTaskGammaConvDalitzV1();
	AliAnalysisTaskGammaConvDalitzV1( const char* name );
	virtual ~AliAnalysisTaskGammaConvDalitzV1();

	virtual void UserExec(Option_t *);
	virtual void UserCreateOutputObjects();
	virtual Bool_t Notify();
	virtual void Terminate(const Option_t *);


	void SetMoveParticleAccordingToVertex(Bool_t flag){fMoveParticleAccordingToVertex = flag;}
         
	void SetIsHeavyIon(Int_t flag){
		if (flag == 1 || flag ==2 ){
			fIsHeavyIon = 1;    
		} else {
			fIsHeavyIon = 0;    
		}
	}
	
	void SetIsMC(Bool_t isMC){fIsMC=isMC;}
	void SetConversionCutList(Int_t nCuts, TList *CutArray){
		fnCuts= nCuts;
		fCutGammaArray = CutArray;
	}
	void SetElectronCutList(TList *CutArray){
		fCutElectronArray = CutArray;
	}
	void SetMesonCutList(TList *CutArray){
		fCutMesonArray = CutArray;
	}
	void SetDoChicAnalysis(Bool_t flag){ fDoChicAnalysis = flag; }
	void SetDoMesonQA(Bool_t flag){ fDoMesonQA = flag; }
	void SetProductionVertextoVGamma(Bool_t flag) { fSetProductionVertextoVGamma = flag; }
  

	private:

		void InitBack();
		void ProcessPhotonCandidates();
		void ProcessTruePhotonCandidates(AliAODConversionPhoton*);
                void ProcessTrueMesonCandidates(AliAODConversionMother *Pi0Candidate, AliAODConversionPhoton *TrueGammaCandidate, AliAODConversionPhoton *TrueVirtualGammaCandidate);
                void MoveParticleAccordingToVertex(AliAODConversionPhoton* particle,const AliGammaConversionAODBGHandler::GammaConversionVertex *vertex);
    		void ProcessElectronCandidates();
		void ProcessVirtualGammasCandidates();
		void ProcessMCParticles();
		void CountESDTracks();
                void CalculatePi0DalitzCandidates();
                void CalculateBackground();
                void UpdateEventByEventData();
		void FillElectronQAHistos(AliAODConversionMother *Pi0Candidate,AliAODConversionPhoton *Vgamma) const;
                Double_t GetPsiPair( const AliESDtrack *trackPos, const AliESDtrack *trackNeg ) const;
		Bool_t IsDalitz(TParticle *fMCMother) const;
                Bool_t IsPi0DalitzDaughter( Int_t label ) const;
		
		
    AliV0ReaderV1 *fV0Reader;
    AliDalitzElectronSelector* fElecSelector;
    AliGammaConversionAODBGHandler **fBGHandler;
    AliESDEvent *fESDEvent;
    AliMCEvent *fMCEvent;
    AliStack *fMCStack;
    TList **fCutFolder;
    TList **fESDList;
    TList **fBackList;
    TList **fMotherList;
    TList **fTrueList;
    TList **fMCList;
    TList **fQAFolder;
    TList *fOutputContainer;
    TClonesArray *fReaderGammas;
    vector<Int_t> fSelectorElectronIndex;
    vector<Int_t> fSelectorPositronIndex;
    TList *fGoodGammas;
    TList *fGoodVirtualGammas;
    TList *fGoodElectrons;
    TList *fGoodPositrons;
    TList *fCutGammaArray;
    TList *fCutElectronArray;
    TList *fCutMesonArray;
    TList **fGammasPool;
    AliConversionCuts *fConversionCuts;
    TH1F **hESDConvGammaPt;
    TH1F **hESDConvGammaEta;
    TH2F **hESDConvGammaZR;
    TH1F **hESDDalitzElectronPt;
    TH1F **hESDDalitzPositronPt;
    TH1F **hESDDalitzElectronPhi;
    TH1F **hESDDalitzPositronPhi;
    TH1F **hESDDalitzElectronAfterPt;
    TH1F **hESDDalitzPositronAfterPt;
    TH1F **hESDDalitzElectronAfterEta;
    TH1F **hESDDalitzPositronAfterEta;
    TH1F **hESDDalitzElectronAfterPhi;
    TH1F **hESDDalitzPositronAfterPhi;
    TH1F **hESDDalitzElectronAfterNClsITS;
    TH1F **hESDDalitzPositronAfterNClsITS;
    TH2F **hESDDalitzElectronAfterNFindClsTPC;
    TH2F **hESDDalitzPositronAfterNFindClsTPC;
    TH2F **hESDDalitzElectronAfterNClsTPC;
    TH2F **hESDDalitzPositronAfterNClsTPC;
    TH2F **hESDDalitzElectronAfterNCrossedRowsTPC;
    TH2F **hESDDalitzPositronAfterNCrossedRowsTPC;
    TH2F **hESDDalitzPosEleAfterDCAxy;
    TH2F **hESDDalitzPosEleAfterDCAz;
    TH2F **hESDDalitzElectronAfterTPCdEdxVsP;
    TH2F **hESDDalitzPositronAfterTPCdEdxVsP;
    TH2F **hESDDalitzElectronAfterTPCdEdxSignalVsP;
    TH2F **hESDDalitzPositronAfterTPCdEdxSignalVsP;
    TH2F **hESDDalitzElectronAfterTPCdEdxVsEta;
    TH2F **hESDDalitzPositronAfterTPCdEdxVsEta;
    TH2F **hESDDalitzElectronAfterTPCdEdxVsPhi;
    TH2F **hESDDalitzPositronAfterTPCdEdxVsPhi;
    TH1F **hESDMotherPhi;
    TH2F **hESDEposEnegPsiPairDPhi;
    TH2F **hESDEposEnegInvMassPt;
    TH2F **hESDEposEnegAfterMassCutInvMassPi0Pt;
    TH2F **hESDEposEnegLikeSignBackInvMassPt;
    TH2F **hESDMotherInvMassPt;
    TH2F **hESDPi0MotherInvMassPt;
    TH2F **hESDPi0MotherDiffInvMassPt;
    TH2F **hESDPi0MotherDiffLimInvMassPt;
    THnSparseF **sESDMotherInvMassPtZM;
    TH2F **hESDMotherBackInvMassPt;
    THnSparseF **sESDMotherBackInvMassPtZM;
    TH1F **hMCAllGammaPt;
    TH1F **hMCConvGammaPt;
    TH1F **hMCConvGammaRSPt;
    TH1F **hMCAllPositronsPt;
    TH1F **hMCAllElectronsPt;
    TH1F **hMCConvGammaEta;
    TH1F **hMCAllPositronsEta;
    TH1F **hMCAllElectronsEta;
    TH1F **hMCPi0DalitzGammaPt;
    TH1F **hMCPi0DalitzElectronPt;
    TH1F **hMCPi0DalitzPositronPt;
    TH1F **hMCPi0Pt;
    TH1F **hMCPi0GGPt;
    TH1F **hMCEtaPt;
    TH1F **hMCEtaGGPt;
    TH1F **hMCPi0InAccPt;
    TH1F **hMCEtaInAccPt;
    TH1F **hMCChiCPt;
    TH1F **hMCChiCInAccPt;
    TH2F **hMCPi0EposEnegInvMassPt;
    TH2F **hMCEtaEposEnegInvMassPt;
    TH2F **hESDEposEnegTruePi0DalitzInvMassPt;
    TH1F **hESDEposEnegTruePrimPi0DalitzInvMass;
    TH2F **hESDEposEnegTruePi0DalitzPsiPairDPhi;
    TH2F **hESDEposEnegTrueEtaDalitzInvMassPt;
    TH1F **hESDEposEnegTruePrimEtaDalitzInvMass;
    TH2F **hESDEposEnegTrueEtaDalitzPsiPairDPhi;
    TH2F **hESDEposEnegTruePhotonInvMassPt;
    TH2F **hESDEposEnegTrueInvMassPt;
    TH2F **hESDEposEnegTruePhotonPsiPairDPhi;
    TH2F **hESDEposEnegTruePhotonPsiPairDPhiPtCut;
    TH2F **hESDEposEnegTrueJPsiInvMassPt;
    TH2F **hESDTrueMotherChiCInvMassPt;
    TH2F **hESDTrueMotherChiCDiffInvMassPt;
    TH2F **hESDTrueMotherInvMassPt;
    TH2F **hESDTrueMotherDalitzInvMassPt;
    TH2F **hESDTrueMotherPi0GGInvMassPt;
    TH2F **hESDTruePrimaryMotherPi0GGInvMassPt;
    TH2F **hESDTrueSecondaryMotherPi0GGInvMassPt;
    TH2F **hESDTruePrimaryMotherInvMassMCPt;
    TH2F **hESDTruePrimaryMotherInvMassPt;
    TH2F **hESDTruePrimaryMotherW0WeightingInvMassPt;
    TH2F **hESDTruePrimaryPi0DalitzESDPtMCPt;
    TH2F **hESDTrueSecondaryMotherInvMassPt;
    TH2F **hESDTrueSecondaryMotherFromK0sInvMassPt;
    TH2F **hESDTrueBckGGInvMassPt;
    TH2F **hESDTrueBckContInvMassPt;
    TH2F **hESDTrueMotherGGInvMassPt;
    TH1F **hESDTrueConvGammaPt;
    TH1F **hESDTruePositronPt;
    TH1F **hESDTrueElectronPt;
    TH1F **hESDTrueSecConvGammaPt;
    TH1F **hESDTrueSecPositronPt;
    TH1F **hESDTrueSecElectronPt;
    TH1F **hESDTruePi0DalitzConvGammaPt;
    TH1F **hESDTruePi0DalitzPositronPt;
    TH1F **hESDTruePi0DalitzElectronPt;
    TH1F **hESDTruePi0DalitzSecConvGammaPt;
    TH1F **hESDTruePi0DalitzSecPositronPt;
    TH1F **hESDTruePi0DalitzSecElectronPt;
      //if(fDoMesonAnalysis){


    TH1I **hNEvents;
    TH1I **hNGoodESDTracks;
    TH2F **hNGoodESDTracksVsNGoodGammas;
    TH2F **hNGoodESDTracksVsNGoodVGammas;
    TH1I **hNV0Tracks;
    TProfile **hEtaShift;
        
    TRandom3 fRandom;
    Double_t *fUnsmearedPx;
    Double_t *fUnsmearedPy;
    Double_t *fUnsmearedPz;
    Double_t *fUnsmearedE;
    Double_t *fUnsmearedVPx;
    Double_t *fUnsmearedVPy;
    Double_t *fUnsmearedVPz;
    Double_t *fUnsmearedVE;
    
    
    Int_t fnCuts;
    Int_t fiCut;
    Int_t fNumberOfESDTracks;
    Int_t fNumberOfESDTrackskBoth;
    Bool_t fMoveParticleAccordingToVertex;
    Bool_t fIsHeavyIon;
    Bool_t fDoMesonAnalysis;
    Bool_t fDoChicAnalysis;
    Bool_t fDoMesonQA;
    Bool_t fSetProductionVertextoVGamma;
    Bool_t fIsFromMBHeader;
    Bool_t fIsMC;

	private:
		AliAnalysisTaskGammaConvDalitzV1( const AliAnalysisTaskGammaConvDalitzV1& ); // Not implemented
		AliAnalysisTaskGammaConvDalitzV1& operator=( const AliAnalysisTaskGammaConvDalitzV1& ); // Not implemented

		ClassDef( AliAnalysisTaskGammaConvDalitzV1, 3 );
};

#endif // ALIANALYSISTASKGAMMACONVDALITZV1_H
Commit	Line	Data
2bb2434e	1	#ifndef ALIANALYSISTASKGAMMACONVDALITZV1_H
	2	#define ALIANALYSISTASKGAMMACONVDALITZV1_H
	3
	4	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	5	* See cxx source for full Copyright notice */
	6
	7	// Analysis task for pi0->e+e-gamma (Dalitz decay)
	8
	9	#include "AliAnalysisTaskSE.h"
	10	#include "AliV0ReaderV1.h"
	11	#include "AliKFConversionPhoton.h"
	12	#include "AliDalitzElectronSelector.h"
	13	#include "AliConversionMesonCuts.h"
	14	#include "AliGammaConversionAODBGHandler.h"
8a52eed5	15	#include "TProfile2D.h"
2bb2434e	16
	17	class AliESDInputHandler;
	18	class AliMCEventHandler;
	19	class AliESDEvent;
	20	class AliESDtrack;
	21	class AliESDtrackCuts;
	22	class AliESDpidCuts;
	23	class AliV0Reader;
	24	class AliGammaConversionHistograms;
	25	class AliTriggerAnalysis;
	26
	27	class AliAnalysisTaskGammaConvDalitzV1: public AliAnalysisTaskSE
	28	{
	29	public:
	30
	31	AliAnalysisTaskGammaConvDalitzV1();
	32	AliAnalysisTaskGammaConvDalitzV1( const char* name );
	33	virtual ~AliAnalysisTaskGammaConvDalitzV1();
	34
	35	virtual void UserExec(Option_t *);
	36	virtual void UserCreateOutputObjects();
ae947965	37	virtual Bool_t Notify();
2bb2434e	38	virtual void Terminate(const Option_t *);
	39
	40
	41
	42
ae947965	43	void SetMoveParticleAccordingToVertex(Bool_t flag){fMoveParticleAccordingToVertex = flag;}
	44
	45	void SetIsHeavyIon(Int_t flag){
	46	if (flag == 1 \|\| flag ==2 ){
	47	fIsHeavyIon = 1;
	48	} else {
	49	fIsHeavyIon = 0;
	50	}
	51	}
	52
	53	void SetIsMC(Bool_t isMC){fIsMC=isMC;}
	54	void SetConversionCutList(Int_t nCuts, TList *CutArray){
	55	fnCuts= nCuts;
	56	fCutGammaArray = CutArray;
	57	}
	58	void SetElectronCutList(TList *CutArray){
	59	fCutElectronArray = CutArray;
	60	}
	61	void SetMesonCutList(TList *CutArray){
	62	fCutMesonArray = CutArray;
	63	}
fcc79cf5	64	void SetDoChicAnalysis(Bool_t flag){ fDoChicAnalysis = flag; }
fcc79cf5	65	void SetDoMesonQA(Bool_t flag){ fDoMesonQA = flag; }
2d49ea4c	66	void SetProductionVertextoVGamma(Bool_t flag) { fSetProductionVertextoVGamma = flag; }
ae947965	67
2bb2434e	68
	69	private:
	70
2bb2434e	71	void InitBack();
	72	void ProcessPhotonCandidates();
	73	void ProcessTruePhotonCandidates(AliAODConversionPhoton*);
	74	void ProcessTrueMesonCandidates(AliAODConversionMother Pi0Candidate, AliAODConversionPhoton TrueGammaCandidate, AliAODConversionPhoton *TrueVirtualGammaCandidate);
	75	void MoveParticleAccordingToVertex(AliAODConversionPhoton* particle,const AliGammaConversionAODBGHandler::GammaConversionVertex *vertex);
	76	void ProcessElectronCandidates();
30b9a400	77	void ProcessVirtualGammasCandidates();
2bb2434e	78	void ProcessMCParticles();
	79	void CountESDTracks();
	80	void CalculatePi0DalitzCandidates();
	81	void CalculateBackground();
	82	void UpdateEventByEventData();
30b9a400	83	void FillElectronQAHistos(AliAODConversionMother Pi0Candidate,AliAODConversionPhoton Vgamma) const;
a280ac15	84	Double_t GetPsiPair( const AliESDtrack trackPos, const AliESDtrack trackNeg ) const;
ae947965	85	Bool_t IsDalitz(TParticle *fMCMother) const;
ccfa8c0d	86	Bool_t IsPi0DalitzDaughter( Int_t label ) const;
2bb2434e	87
	88
	89
	90	AliV0ReaderV1 *fV0Reader;
	91	AliDalitzElectronSelector* fElecSelector;
	92	AliGammaConversionAODBGHandler **fBGHandler;
	93	AliESDEvent *fESDEvent;
	94	AliMCEvent *fMCEvent;
	95	AliStack *fMCStack;
	96	TList **fCutFolder;
	97	TList **fESDList;
	98	TList **fBackList;
fcc79cf5	99	TList **fMotherList;
2bb2434e	100	TList **fTrueList;
2bb2434e	101	TList **fMCList;
d9d6352b	102	TList **fQAFolder;
2bb2434e	103	TList *fOutputContainer;
	104	TClonesArray *fReaderGammas;
	105	vector<Int_t> fSelectorElectronIndex;
	106	vector<Int_t> fSelectorPositronIndex;
	107	TList *fGoodGammas;
	108	TList *fGoodVirtualGammas;
	109	TList *fGoodElectrons;
	110	TList *fGoodPositrons;
	111	TList *fCutGammaArray;
	112	TList *fCutElectronArray;
	113	TList *fCutMesonArray;
	114	TList **fGammasPool;
	115	AliConversionCuts *fConversionCuts;
	116	TH1F **hESDConvGammaPt;
2279d237	117	TH1F **hESDConvGammaEta;
2279d237	118	TH2F **hESDConvGammaZR;
2bb2434e	119	TH1F **hESDDalitzElectronPt;
2bb2434e	120	TH1F **hESDDalitzPositronPt;
4803eb1f	121	TH1F **hESDDalitzElectronPhi;
	122	TH1F **hESDDalitzPositronPhi;
	123	TH1F **hESDDalitzElectronAfterPt;
	124	TH1F **hESDDalitzPositronAfterPt;
2279d237	125	TH1F **hESDDalitzElectronAfterEta;
2279d237	126	TH1F **hESDDalitzPositronAfterEta;
4803eb1f	127	TH1F **hESDDalitzElectronAfterPhi;
4803eb1f	128	TH1F **hESDDalitzPositronAfterPhi;
d9d6352b	129	TH1F **hESDDalitzElectronAfterNClsITS;
d9d6352b	130	TH1F **hESDDalitzPositronAfterNClsITS;
4803eb1f	131	TH2F **hESDDalitzElectronAfterNFindClsTPC;
4803eb1f	132	TH2F **hESDDalitzPositronAfterNFindClsTPC;
d9d6352b	133	TH2F **hESDDalitzElectronAfterNClsTPC;
d9d6352b	134	TH2F **hESDDalitzPositronAfterNClsTPC;
f47b2bb5	135	TH2F **hESDDalitzElectronAfterNCrossedRowsTPC;
f47b2bb5	136	TH2F **hESDDalitzPositronAfterNCrossedRowsTPC;
4803eb1f	137	TH2F **hESDDalitzPosEleAfterDCAxy;
4803eb1f	138	TH2F **hESDDalitzPosEleAfterDCAz;
d9d6352b	139	TH2F **hESDDalitzElectronAfterTPCdEdxVsP;
	140	TH2F **hESDDalitzPositronAfterTPCdEdxVsP;
	141	TH2F **hESDDalitzElectronAfterTPCdEdxSignalVsP;
	142	TH2F **hESDDalitzPositronAfterTPCdEdxSignalVsP;
	143	TH2F **hESDDalitzElectronAfterTPCdEdxVsEta;
	144	TH2F **hESDDalitzPositronAfterTPCdEdxVsEta;
	145	TH2F **hESDDalitzElectronAfterTPCdEdxVsPhi;
	146	TH2F **hESDDalitzPositronAfterTPCdEdxVsPhi;
4803eb1f	147	TH1F **hESDMotherPhi;
2bb2434e	148	TH2F **hESDEposEnegPsiPairDPhi;
a280ac15	149	TH2F **hESDEposEnegInvMassPt;
30b9a400	150	TH2F **hESDEposEnegAfterMassCutInvMassPi0Pt;
0a2b2b4b	151	TH2F **hESDEposEnegLikeSignBackInvMassPt;
2bb2434e	152	TH2F **hESDMotherInvMassPt;
a280ac15	153	TH2F **hESDPi0MotherInvMassPt;
a280ac15	154	TH2F **hESDPi0MotherDiffInvMassPt;
ae947965	155	TH2F **hESDPi0MotherDiffLimInvMassPt;
2bb2434e	156	THnSparseF **sESDMotherInvMassPtZM;
	157	TH2F **hESDMotherBackInvMassPt;
	158	THnSparseF **sESDMotherBackInvMassPtZM;
ccfa8c0d	159	TH1F **hMCAllGammaPt;
ae947965	160	TH1F **hMCConvGammaPt;
ae947965	161	TH1F **hMCConvGammaRSPt;
ccfa8c0d	162	TH1F **hMCAllPositronsPt;
ccfa8c0d	163	TH1F **hMCAllElectronsPt;
2d49ea4c	164	TH1F **hMCConvGammaEta;
	165	TH1F **hMCAllPositronsEta;
	166	TH1F **hMCAllElectronsEta;
ccfa8c0d	167	TH1F **hMCPi0DalitzGammaPt;
	168	TH1F **hMCPi0DalitzElectronPt;
	169	TH1F **hMCPi0DalitzPositronPt;
2bb2434e	170	TH1F **hMCPi0Pt;
a280ac15	171	TH1F **hMCPi0GGPt;
2bb2434e	172	TH1F **hMCEtaPt;
a280ac15	173	TH1F **hMCEtaGGPt;
2bb2434e	174	TH1F **hMCPi0InAccPt;
2bb2434e	175	TH1F **hMCEtaInAccPt;
0f8c33c1	176	TH1F **hMCChiCPt;
0f8c33c1	177	TH1F **hMCChiCInAccPt;
f4d5bd23	178	TH2F **hMCPi0EposEnegInvMassPt;
f4d5bd23	179	TH2F **hMCEtaEposEnegInvMassPt;
ccfa8c0d	180	TH2F **hESDEposEnegTruePi0DalitzInvMassPt;
2d49ea4c	181	TH1F **hESDEposEnegTruePrimPi0DalitzInvMass;
a072aeaa	182	TH2F **hESDEposEnegTruePi0DalitzPsiPairDPhi;
ccfa8c0d	183	TH2F **hESDEposEnegTrueEtaDalitzInvMassPt;
2d49ea4c	184	TH1F **hESDEposEnegTruePrimEtaDalitzInvMass;
a072aeaa	185	TH2F **hESDEposEnegTrueEtaDalitzPsiPairDPhi;
ccfa8c0d	186	TH2F **hESDEposEnegTruePhotonInvMassPt;
f4d5bd23	187	TH2F **hESDEposEnegTrueInvMassPt;
a072aeaa	188	TH2F **hESDEposEnegTruePhotonPsiPairDPhi;
30b9a400	189	TH2F **hESDEposEnegTruePhotonPsiPairDPhiPtCut;
0f8c33c1	190	TH2F **hESDEposEnegTrueJPsiInvMassPt;
0f8c33c1	191	TH2F **hESDTrueMotherChiCInvMassPt;
ae947965	192	TH2F **hESDTrueMotherChiCDiffInvMassPt;
2bb2434e	193	TH2F **hESDTrueMotherInvMassPt;
4803eb1f	194	TH2F **hESDTrueMotherDalitzInvMassPt;
2bb2434e	195	TH2F **hESDTrueMotherPi0GGInvMassPt;
86eaf9ae	196	TH2F **hESDTruePrimaryMotherPi0GGInvMassPt;
86eaf9ae	197	TH2F **hESDTrueSecondaryMotherPi0GGInvMassPt;
2bb2434e	198	TH2F **hESDTruePrimaryMotherInvMassMCPt;
86eaf9ae	199	TH2F **hESDTruePrimaryMotherInvMassPt;
86eaf9ae	200	TH2F **hESDTruePrimaryMotherW0WeightingInvMassPt;
2bb2434e	201	TH2F **hESDTruePrimaryPi0DalitzESDPtMCPt;
	202	TH2F **hESDTrueSecondaryMotherInvMassPt;
	203	TH2F **hESDTrueSecondaryMotherFromK0sInvMassPt;
	204	TH2F **hESDTrueBckGGInvMassPt;
	205	TH2F **hESDTrueBckContInvMassPt;
	206	TH2F **hESDTrueMotherGGInvMassPt;
	207	TH1F **hESDTrueConvGammaPt;
ccfa8c0d	208	TH1F **hESDTruePositronPt;
ccfa8c0d	209	TH1F **hESDTrueElectronPt;
ae947965	210	TH1F **hESDTrueSecConvGammaPt;
	211	TH1F **hESDTrueSecPositronPt;
	212	TH1F **hESDTrueSecElectronPt;
ccfa8c0d	213	TH1F **hESDTruePi0DalitzConvGammaPt;
	214	TH1F **hESDTruePi0DalitzPositronPt;
	215	TH1F **hESDTruePi0DalitzElectronPt;
ae947965	216	TH1F **hESDTruePi0DalitzSecConvGammaPt;
	217	TH1F **hESDTruePi0DalitzSecPositronPt;
	218	TH1F **hESDTruePi0DalitzSecElectronPt;
ccfa8c0d	219	//if(fDoMesonAnalysis){
	220
	221
2bb2434e	222	TH1I **hNEvents;
2bb2434e	223	TH1I **hNGoodESDTracks;
58326f11	224	TH2F **hNGoodESDTracksVsNGoodGammas;
58326f11	225	TH2F **hNGoodESDTracksVsNGoodVGammas;
7b3733f6	226	TH1I **hNV0Tracks;
8a52eed5	227	TProfile **hEtaShift;
2bb2434e	228
	229	TRandom3 fRandom;
	230	Double_t *fUnsmearedPx;
	231	Double_t *fUnsmearedPy;
	232	Double_t *fUnsmearedPz;
	233	Double_t *fUnsmearedE;
30b9a400	234	Double_t *fUnsmearedVPx;
	235	Double_t *fUnsmearedVPy;
	236	Double_t *fUnsmearedVPz;
	237	Double_t *fUnsmearedVE;
	238
	239
2bb2434e	240	Int_t fnCuts;
	241	Int_t fiCut;
	242	Int_t fNumberOfESDTracks;
58326f11	243	Int_t fNumberOfESDTrackskBoth;
2bb2434e	244	Bool_t fMoveParticleAccordingToVertex;
	245	Bool_t fIsHeavyIon;
	246	Bool_t fDoMesonAnalysis;
fcc79cf5	247	Bool_t fDoChicAnalysis;
fcc79cf5	248	Bool_t fDoMesonQA;
2d49ea4c	249	Bool_t fSetProductionVertextoVGamma;
ae947965	250	Bool_t fIsFromMBHeader;
2d49ea4c	251	Bool_t fIsMC;
ae947965	252
2bb2434e	253	private:
	254	AliAnalysisTaskGammaConvDalitzV1( const AliAnalysisTaskGammaConvDalitzV1& ); // Not implemented
	255	AliAnalysisTaskGammaConvDalitzV1& operator=( const AliAnalysisTaskGammaConvDalitzV1& ); // Not implemented
	256
ae947965	257	ClassDef( AliAnalysisTaskGammaConvDalitzV1, 3 );
2bb2434e	258	};
	259
	260	#endif // ALIANALYSISTASKGAMMACONVDALITZV1_H
	261