[u/mrichter/AliRoot.git] / PWGJE / EMCALJetTasks / UserTasks / AliAnalysisTaskSAQA.h

#ifndef ALIANALYSISTASKSAQA_H
#define ALIANALYSISTASKSAQA_H

// $Id$

class TClonesArray;
class TString;
class TH1F;
class TH2F;
class TH3F;

#include "AliAnalysisTaskEmcalJet.h"

class AliAnalysisTaskSAQA : public AliAnalysisTaskEmcalJet {
 public:
  AliAnalysisTaskSAQA();
  AliAnalysisTaskSAQA(const char *name);
  virtual ~AliAnalysisTaskSAQA();

  void                        UserCreateOutputObjects();
  void                        Terminate(Option_t *option);

  void                        SetCellEnergyCut(Float_t cut)                        { fCellEnergyCut      = cut        ; }

 protected:

  Bool_t                      FillHistograms()                                              ;
  Bool_t                      RetrieveEventObjects()                                        ;
  Int_t                       DoCellLoop(Float_t &sum, Float_t &sum_cut)                    ;
  Float_t                     DoTrackLoop()                                                 ;
  Float_t                     DoClusterLoop()                                               ;
  void                        DoJetLoop()                                                   ;

  Float_t                     fCellEnergyCut;            // Energy cell cut
  Int_t                       fNclusters;                //!Number of accepted clusters in the event
  Int_t                       fNtracks;                  //!Number of accepted tracks in the event
  Int_t                       fNjets;                    //!Number of accepted jets in the event
 
  // General histograms
  TH2F                       *fHistTracksCent;           //!Number of tracks vs. centrality
  TH2F                       *fHistClusCent;             //!Number of clusters vs. centrality
  TH2F                       *fHistJetsCent;             //!Number of jets vs. centrality
  TH2F                       *fHistClusTracks;           //!Number of clusters vs. number of tracks
  TH2F                       *fHistJetsParts;            //!Number of jets vs. number of particles (tracks+clusters)
  TH2F                       *fHistCellsCent;            //!Number of cells vs. centrality
  TH2F                       *fHistCellsTracks;          //!Number of cells vs. number of tracks

  // Tracks
  TH3F                       *fHistTrPhiEtaPt[4];        //!Phi-Eta-Pt distribution of tracks
  TH2F                       *fHistTrEmcPhiEta;          //!Phi-Eta emcal propagated distribution of tracks
  TH2F                       *fHistTrPhiEtaNonProp;      //!Phi-Eta distribution of non emcal propagated tracks
  TH2F                       *fHistDeltaEtaPt;           //!Eta-EtaProp vs. Pt
  TH2F                       *fHistDeltaPhiPt;           //!Phi-PhiProp vs. Pt

  // Clusters
  TH3F                       *fHistClusPhiEtaEnergy[4];  //!Phi-Eta-Energy distribution of clusters
  TH2F                       *fHistNCellsEnergy;         //!Number of cells vs. energy of cluster
  TH2F                       *fHistClusTimeEnergy;       //!Time vs. energy of cluster

  //Jets
  TH3F                       *fHistJetsPhiEtaPt[4];      //!Phi-Eta distribution of jets
  TH2F                       *fHistJetsPtArea[4];        //!Pt vs. area of jets

  // EMCAL Cells
  TH2F                       *fHistCellsAbsIdEnergy;    //!Energy spectrum of cells

  // Had corr QA
  TH2F                       *fHistChVSneCells;          //!Charged vs. neutral (cells) energy
  TH2F                       *fHistChVSneClus;           //!Charged vs. neutral (clusters) energy
  TH2F                       *fHistChVSneCorrCells;      //!Charged vs. neutral (corrected cells) energy

  // Hybrid tracks
  TH1F                       *fHistTrackPhi[5];          //!Phi distribution of hybrid tracks
  TH1F                       *fHistTrackEta[5];          //!Eta distribution of hybrid tracks

 private:
  AliAnalysisTaskSAQA(const AliAnalysisTaskSAQA&);            // not implemented
  AliAnalysisTaskSAQA &operator=(const AliAnalysisTaskSAQA&); // not implemented

  ClassDef(AliAnalysisTaskSAQA, 14) // Quality task for Emcal analysis
};
#endif
Commit	Line	Data
c3ba2d3d	1	#ifndef ALIANALYSISTASKSAQA_H
	2	#define ALIANALYSISTASKSAQA_H
	3
6e8d91c9	4	// $Id$
c3ba2d3d	5
	6	class TClonesArray;
	7	class TString;
	8	class TH1F;
	9	class TH2F;
3000c095	10	class TH3F;
c3ba2d3d	11
6fd5039f	12	#include "AliAnalysisTaskEmcalJet.h"
c3ba2d3d	13
6fd5039f	14	class AliAnalysisTaskSAQA : public AliAnalysisTaskEmcalJet {
c3ba2d3d	15	public:
	16	AliAnalysisTaskSAQA();
	17	AliAnalysisTaskSAQA(const char *name);
	18	virtual ~AliAnalysisTaskSAQA();
	19
	20	void UserCreateOutputObjects();
	21	void Terminate(Option_t *option);
	22
11d18b51	23	void SetCellEnergyCut(Float_t cut) { fCellEnergyCut = cut ; }
c3ba2d3d	24
	25	protected:
	26
079b4732	27	Bool_t FillHistograms() ;
	28	Bool_t RetrieveEventObjects() ;
	29	Int_t DoCellLoop(Float_t &sum, Float_t &sum_cut) ;
079b4732	30	Float_t DoTrackLoop() ;
	31	Float_t DoClusterLoop() ;
	32	void DoJetLoop() ;
c3ba2d3d	33
3e951e32	34	Float_t fCellEnergyCut; // Energy cell cut
7c1d624c	35	Int_t fNclusters; //!Number of accepted clusters in the event
	36	Int_t fNtracks; //!Number of accepted tracks in the event
	37	Int_t fNjets; //!Number of accepted jets in the event
	38
3e951e32	39	// General histograms
3e951e32	40	TH2F *fHistTracksCent; //!Number of tracks vs. centrality
3e951e32	41	TH2F *fHistClusCent; //!Number of clusters vs. centrality
7c1d624c	42	TH2F *fHistJetsCent; //!Number of jets vs. centrality
3e951e32	43	TH2F *fHistClusTracks; //!Number of clusters vs. number of tracks
7c1d624c	44	TH2F *fHistJetsParts; //!Number of jets vs. number of particles (tracks+clusters)
3e951e32	45	TH2F *fHistCellsCent; //!Number of cells vs. centrality
3e951e32	46	TH2F *fHistCellsTracks; //!Number of cells vs. number of tracks
a487deae	47
3e951e32	48	// Tracks
a487deae	49	TH3F *fHistTrPhiEtaPt[4]; //!Phi-Eta-Pt distribution of tracks
3e951e32	50	TH2F *fHistTrEmcPhiEta; //!Phi-Eta emcal propagated distribution of tracks
079b4732	51	TH2F *fHistTrPhiEtaNonProp; //!Phi-Eta distribution of non emcal propagated tracks
3e951e32	52	TH2F *fHistDeltaEtaPt; //!Eta-EtaProp vs. Pt
3e951e32	53	TH2F *fHistDeltaPhiPt; //!Phi-PhiProp vs. Pt
a487deae	54
3e951e32	55	// Clusters
a487deae	56	TH3F *fHistClusPhiEtaEnergy[4]; //!Phi-Eta-Energy distribution of clusters
7cf4626b	57	TH2F *fHistNCellsEnergy; //!Number of cells vs. energy of cluster
090a0c3e	58	TH2F *fHistClusTimeEnergy; //!Time vs. energy of cluster
a487deae	59
3e951e32	60	//Jets
a487deae	61	TH3F *fHistJetsPhiEtaPt[4]; //!Phi-Eta distribution of jets
	62	TH2F *fHistJetsPtArea[4]; //!Pt vs. area of jets
	63
3e951e32	64	// EMCAL Cells
f483218e	65	TH2F *fHistCellsAbsIdEnergy; //!Energy spectrum of cells
a487deae	66
3e951e32	67	// Had corr QA
	68	TH2F *fHistChVSneCells; //!Charged vs. neutral (cells) energy
	69	TH2F *fHistChVSneClus; //!Charged vs. neutral (clusters) energy
	70	TH2F *fHistChVSneCorrCells; //!Charged vs. neutral (corrected cells) energy
a487deae	71
3e951e32	72	// Hybrid tracks
	73	TH1F *fHistTrackPhi[5]; //!Phi distribution of hybrid tracks
	74	TH1F *fHistTrackEta[5]; //!Eta distribution of hybrid tracks
c3ba2d3d	75
	76	private:
	77	AliAnalysisTaskSAQA(const AliAnalysisTaskSAQA&); // not implemented
	78	AliAnalysisTaskSAQA &operator=(const AliAnalysisTaskSAQA&); // not implemented
	79
a487deae	80	ClassDef(AliAnalysisTaskSAQA, 14) // Quality task for Emcal analysis
c3ba2d3d	81	};
c3ba2d3d	82	#endif