[u/mrichter/AliRoot.git] / PWGGA / EMCALTasks / AliAnalysisTaskPi0V2.h

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

/* $Id: AliAnalysisTaskPi0V2.h 45956 2010-12-10 12:55:37Z agheata $ */
/* AliAnalysisTaskPi0V2.h
 *
 * Template task producing a P_t spectrum and pseudorapidity distribution.
 * Includes explanations of physics and primary track selections
 *
 * Based on tutorial example from offline pages
 * Edited by Arvinder Palaha
 */
#ifndef ALIANALYSISTASKEX01_H
#define ALIANALYSISTASKEX01_H

class TH1F;
class TH2F;
class TH3F;
class TList;
class AliESDCaloCluster;
class AliESDtrackCuts;
class AliESDEvent;
class THnSparse;

#ifndef ALIANALYSISTASKSE_H
#include "AliAnalysisTaskSE.h"
#endif

class AliAnalysisTaskPi0V2 : public AliAnalysisTaskSE {
 public:
    AliAnalysisTaskPi0V2();
    AliAnalysisTaskPi0V2(const char *name);
    virtual ~AliAnalysisTaskPi0V2();
    
    virtual void     UserCreateOutputObjects();
    virtual void     UserExec(Option_t *option);
    virtual void     Terminate(Option_t *);
    Double_t		GetCrossEnergy(const AliVCluster *cluster, Short_t &idmax) const;
    Double_t		GetMaxCellEnergy(const AliVCluster *cluster, Short_t &id) const;
    Bool_t		IsWithinFiducialVolume(Short_t id) const;
    Bool_t		IsGoodCluster(const AliESDCaloCluster *c) const;
    Bool_t		IsGoodPion(const TLorentzVector& p1, const TLorentzVector& p2) const;
    void		FillPion(const TLorentzVector& p1, const TLorentzVector& p2, Double_t EPV0r, Double_t EPV0A, Double_t EPV0C, Double_t EPTPC);
    void 		GetMom(TLorentzVector& p, const AliESDCaloCluster *c, Double_t *vertex);		
    void		SetsubEventMethod(Bool_t b )	{ fcheckEP2sub =b ;}

    
 private:
    TList           		*fOutput;        // Output list
    AliESDtrackCuts 		*fTrackCuts;     // Track cuts
    AliESDEvent			*fESD;          //!ESD object

    Bool_t 			fcheckEP2sub;	//! do 2 sub event method
    // NEW HISTO to be declared here
    Double_t			fCentrality;	//! Centrality
    Double_t			fEPTPC;	//! Evt plane TPC
    Double_t			fEPTPCreso;	//! resolution of TPC method
    Double_t			fEPV0;		//! EP V0
    Double_t			fEPV0A;	//! EP V0A
    Double_t			fEPV0C;	//! EP V0C
    Double_t			fEPV0Ar;	//! EP V0A reduced
    Double_t			fEPV0Cr;	//! EP V0C reduced
    Double_t			fEPV0r;	//! EP V0 reduced
    Double_t			fEPV0AR4;	//! EP V0A ring4 only
    Double_t			fEPV0AR5;	//! EP V0A ring5 only
    Double_t			fEPV0AR6;	//! EP V0A ring6 only
    Double_t			fEPV0AR7;	//! EP V0A ring7 only
    Double_t			fEPV0CR0;	//! EP V0C ring0 only	
    Double_t			fEPV0CR1;	//! EP V0C ring1 only	
    Double_t			fEPV0CR2;	//! EP V0C ring2 only	
    Double_t			fEPV0CR3;	//! EP V0C ring3 only	

    TH2F			*hEPTPC;	//! 2-D histo EPTPC  vs cent
    TH2F			*hresoTPC;	//! 2-D histo TPC resolution vs cent
    TH2F			*hEPV0;		//! 2-D histo EPV0   vs cent
    TH2F			*hEPV0A;	//! 2-D histo EPV0A  vs cent
    TH2F			*hEPV0C;	//! 2-D histo EPV0C  vs cent
    TH2F			*hEPV0Ar;	//! 2-D histo EPV0Ar vs cent
    TH2F			*hEPV0Cr;	//! 2-D histo EPV0Cr vs cent
    TH2F			*hEPV0r;	//! 2-D histo EPV0r  vs cent
    TH2F			*hEPV0AR4;	//! 2-D histo EPV0AR4 vs cent
    TH2F			*hEPV0AR7;	//! 2-D histo EPV0AR7 vs cent
    TH2F			*hEPV0CR0;	//! 2-D histo EPV0AR0 vs cent
    TH2F			*hEPV0CR3;	//! 2-D histo EPV0AR3 vs cent

    TH2F			*hdifV0A_V0CR0;		//! 2-D histo diff V0A, V0CR0 vs cent
    TH2F			*hdifV0A_V0CR3;		//! 2-D histo diff V0A, V0CR3 vs cent
    TH2F			*hdifV0ACR0_V0CR3;	//! 2-D histo diff V0CR0, V0CR3 vs cent
    TH2F			*hdifV0C_V0AR4;		//! 2-D histo diff V0C, V0AR4 vs cent
    TH2F			*hdifV0C_V0AR7;		//! 2-D histo diff V0C, V0AR7 vs cent
    TH2F			*hdifV0AR4_V0AR7;	//! 2-D histo diff V0AR7, V0AR4 vs cent

    TH2F			*hdifV0A_V0C;		//! 2-D histo diff V0A - V0C
    TH2F			*hdifV0A_TPC;		//! 2-D histo diff V0A - TPC
    TH2F			*hdifV0C_TPC;		//! 2-D histo diff V0C - TPC
    TH2F			*hdifV0C_V0A;		//! 2-D histo diff V0C - V0A

    TH3F			*hdifEMC_EP;		//! 2-D histo dif phi in EMC with EP
    TH3F			*hdifful_EP;		//! 2-D histo dif phi in full with EP
    TH3F			*hdifout_EP;		//! 2-D histo dif phi NOT in EMC with EP

    THnSparse                   *fHEPV0r;         //! Flow 4-D Histo
    THnSparse                   *fHEPV0A;        //! Flow 4-D Histo
    THnSparse                   *fHEPV0C;        //! Flow 4-D Histo
    THnSparse                   *fHEPTPC;        //! Flow 4-D Histo

    
    
    AliAnalysisTaskPi0V2(const AliAnalysisTaskPi0V2&); // not implemented
    AliAnalysisTaskPi0V2& operator=(const AliAnalysisTaskPi0V2&); // not implemented
    
    ClassDef(AliAnalysisTaskPi0V2, 1); // example of analysis
};

#endif
Commit	Line	Data
04b116e8	1	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
04b116e8	2	* See cxx source for full Copyright notice */
3c40321c	3
04b116e8	4	/* $Id: AliAnalysisTaskPi0V2.h 45956 2010-12-10 12:55:37Z agheata $ */
	5	/* AliAnalysisTaskPi0V2.h
	6	*
	7	* Template task producing a P_t spectrum and pseudorapidity distribution.
	8	* Includes explanations of physics and primary track selections
	9	*
	10	* Based on tutorial example from offline pages
	11	* Edited by Arvinder Palaha
	12	*/
	13	#ifndef ALIANALYSISTASKEX01_H
	14	#define ALIANALYSISTASKEX01_H
3c40321c	15
	16	class TH1F;
	17	class TH2F;
04b116e8	18	class TH3F;
3c40321c	19	class TList;
	20	class AliESDCaloCluster;
	21	class AliESDtrackCuts;
	22	class AliESDEvent;
04b116e8	23	class THnSparse;
3c40321c	24
04b116e8	25	#ifndef ALIANALYSISTASKSE_H
3c40321c	26	#include "AliAnalysisTaskSE.h"
04b116e8	27	#endif
3c40321c	28
	29	class AliAnalysisTaskPi0V2 : public AliAnalysisTaskSE {
	30	public:
04b116e8	31	AliAnalysisTaskPi0V2();
	32	AliAnalysisTaskPi0V2(const char *name);
	33	virtual ~AliAnalysisTaskPi0V2();
3c40321c	34
04b116e8	35	virtual void UserCreateOutputObjects();
	36	virtual void UserExec(Option_t *option);
	37	virtual void Terminate(Option_t *);
	38	Double_t GetCrossEnergy(const AliVCluster *cluster, Short_t &idmax) const;
	39	Double_t GetMaxCellEnergy(const AliVCluster *cluster, Short_t &id) const;
	40	Bool_t IsWithinFiducialVolume(Short_t id) const;
	41	Bool_t IsGoodCluster(const AliESDCaloCluster *c) const;
	42	Bool_t IsGoodPion(const TLorentzVector& p1, const TLorentzVector& p2) const;
	43	void FillPion(const TLorentzVector& p1, const TLorentzVector& p2, Double_t EPV0r, Double_t EPV0A, Double_t EPV0C, Double_t EPTPC);
	44	void GetMom(TLorentzVector& p, const AliESDCaloCluster c, Double_t vertex);
	45	void SetsubEventMethod(Bool_t b ) { fcheckEP2sub =b ;}
	46
3c40321c	47
3c40321c	48	private:
04b116e8	49	TList *fOutput; // Output list
	50	AliESDtrackCuts *fTrackCuts; // Track cuts
	51	AliESDEvent *fESD; //!ESD object
	52
	53	Bool_t fcheckEP2sub; //! do 2 sub event method
	54	// NEW HISTO to be declared here
	55	Double_t fCentrality; //! Centrality
	56	Double_t fEPTPC; //! Evt plane TPC
	57	Double_t fEPTPCreso; //! resolution of TPC method
	58	Double_t fEPV0; //! EP V0
	59	Double_t fEPV0A; //! EP V0A
	60	Double_t fEPV0C; //! EP V0C
	61	Double_t fEPV0Ar; //! EP V0A reduced
	62	Double_t fEPV0Cr; //! EP V0C reduced
	63	Double_t fEPV0r; //! EP V0 reduced
	64	Double_t fEPV0AR4; //! EP V0A ring4 only
	65	Double_t fEPV0AR5; //! EP V0A ring5 only
	66	Double_t fEPV0AR6; //! EP V0A ring6 only
	67	Double_t fEPV0AR7; //! EP V0A ring7 only
	68	Double_t fEPV0CR0; //! EP V0C ring0 only
	69	Double_t fEPV0CR1; //! EP V0C ring1 only
	70	Double_t fEPV0CR2; //! EP V0C ring2 only
	71	Double_t fEPV0CR3; //! EP V0C ring3 only
	72
	73	TH2F *hEPTPC; //! 2-D histo EPTPC vs cent
	74	TH2F *hresoTPC; //! 2-D histo TPC resolution vs cent
	75	TH2F *hEPV0; //! 2-D histo EPV0 vs cent
	76	TH2F *hEPV0A; //! 2-D histo EPV0A vs cent
	77	TH2F *hEPV0C; //! 2-D histo EPV0C vs cent
	78	TH2F *hEPV0Ar; //! 2-D histo EPV0Ar vs cent
	79	TH2F *hEPV0Cr; //! 2-D histo EPV0Cr vs cent
	80	TH2F *hEPV0r; //! 2-D histo EPV0r vs cent
	81	TH2F *hEPV0AR4; //! 2-D histo EPV0AR4 vs cent
	82	TH2F *hEPV0AR7; //! 2-D histo EPV0AR7 vs cent
	83	TH2F *hEPV0CR0; //! 2-D histo EPV0AR0 vs cent
	84	TH2F *hEPV0CR3; //! 2-D histo EPV0AR3 vs cent
	85
	86	TH2F *hdifV0A_V0CR0; //! 2-D histo diff V0A, V0CR0 vs cent
	87	TH2F *hdifV0A_V0CR3; //! 2-D histo diff V0A, V0CR3 vs cent
	88	TH2F *hdifV0ACR0_V0CR3; //! 2-D histo diff V0CR0, V0CR3 vs cent
	89	TH2F *hdifV0C_V0AR4; //! 2-D histo diff V0C, V0AR4 vs cent
	90	TH2F *hdifV0C_V0AR7; //! 2-D histo diff V0C, V0AR7 vs cent
	91	TH2F *hdifV0AR4_V0AR7; //! 2-D histo diff V0AR7, V0AR4 vs cent
	92
	93	TH2F *hdifV0A_V0C; //! 2-D histo diff V0A - V0C
	94	TH2F *hdifV0A_TPC; //! 2-D histo diff V0A - TPC
	95	TH2F *hdifV0C_TPC; //! 2-D histo diff V0C - TPC
	96	TH2F *hdifV0C_V0A; //! 2-D histo diff V0C - V0A
	97
	98	TH3F *hdifEMC_EP; //! 2-D histo dif phi in EMC with EP
	99	TH3F *hdifful_EP; //! 2-D histo dif phi in full with EP
	100	TH3F *hdifout_EP; //! 2-D histo dif phi NOT in EMC with EP
	101
	102	THnSparse *fHEPV0r; //! Flow 4-D Histo
	103	THnSparse *fHEPV0A; //! Flow 4-D Histo
	104	THnSparse *fHEPV0C; //! Flow 4-D Histo
	105	THnSparse *fHEPTPC; //! Flow 4-D Histo
	106
	107
3c40321c	108
04b116e8	109	AliAnalysisTaskPi0V2(const AliAnalysisTaskPi0V2&); // not implemented
	110	AliAnalysisTaskPi0V2& operator=(const AliAnalysisTaskPi0V2&); // not implemented
	111
	112	ClassDef(AliAnalysisTaskPi0V2, 1); // example of analysis
3c40321c	113	};
04b116e8	114
3c40321c	115	#endif
04b116e8	116