1 #ifndef ALIANALYSISTASKSAQA_H
2 #define ALIANALYSISTASKSAQA_H
12 #include "AliAnalysisTaskEmcalJet.h"
14 class AliAnalysisTaskSAQA : public AliAnalysisTaskEmcalJet {
16 AliAnalysisTaskSAQA();
17 AliAnalysisTaskSAQA(const char *name);
18 virtual ~AliAnalysisTaskSAQA();
20 void UserCreateOutputObjects();
21 void Terminate(Option_t *option);
23 void SetCellEnergyCut(Float_t cut) { fCellEnergyCut = cut ; }
24 void SetParticleLevel(Bool_t s) { fParticleLevel = s ; }
28 Bool_t FillHistograms() ;
29 Bool_t RetrieveEventObjects() ;
30 Int_t DoCellLoop(Float_t &sum, Float_t &sum_cut) ;
31 Float_t DoTrackLoop() ;
32 Float_t DoClusterLoop() ;
34 Double_t GetFcross(AliVCluster *cluster, AliVCaloCells *cells) ;
36 Float_t fCellEnergyCut; // Energy cell cut
37 Bool_t fParticleLevel; // Set particle level analysis
38 Int_t fNclusters; //!Number of accepted clusters in the event
39 Int_t fNtracks; //!Number of accepted tracks in the event
40 Int_t fNjets; //!Number of accepted jets in the event
43 TH2F *fHistTracksCent; //!Number of tracks vs. centrality
44 TH2F *fHistClusCent; //!Number of clusters vs. centrality
45 TH2F *fHistJetsCent; //!Number of jets vs. centrality
46 TH2F *fHistClusTracks; //!Number of clusters vs. number of tracks
47 TH2F *fHistJetsParts; //!Number of jets vs. number of particles (tracks+clusters)
48 TH2F *fHistCellsCent; //!Number of cells vs. centrality
49 TH2F *fHistCellsTracks; //!Number of cells vs. number of tracks
52 TH3F *fHistTrPhiEtaPt[4][4]; //!Phi-Eta-Pt distribution of tracks
53 TH2F *fHistTrEmcPhiEta; //!Phi-Eta emcal propagated distribution of tracks
54 TH2F *fHistTrPhiEtaNonProp; //!Phi-Eta distribution of non emcal propagated tracks
55 TH2F *fHistDeltaEtaPt; //!Eta-EtaProp vs. Pt
56 TH2F *fHistDeltaPhiPt; //!Phi-PhiProp vs. Pt
59 TH3F *fHistClusPhiEtaEnergy[4]; //!Phi-Eta-Energy distribution of clusters
60 TH2F *fHistNCellsEnergy; //!Number of cells vs. energy of cluster
61 TH2F *fHistFcrossEnergy; //!Fcross vs. energy of cluster
62 TH2F *fHistClusTimeEnergy; //!Time vs. energy of cluster
65 TH3F *fHistJetsPhiEtaPt[4]; //!Phi-Eta distribution of jets
66 TH2F *fHistJetsPtArea[4]; //!Pt vs. area of jets
69 TH2F *fHistCellsAbsIdEnergy; //!Energy spectrum of cells
72 TH2F *fHistChVSneCells; //!Charged vs. neutral (cells) energy
73 TH2F *fHistChVSneClus; //!Charged vs. neutral (clusters) energy
74 TH2F *fHistChVSneCorrCells; //!Charged vs. neutral (corrected cells) energy
77 AliAnalysisTaskSAQA(const AliAnalysisTaskSAQA&); // not implemented
78 AliAnalysisTaskSAQA &operator=(const AliAnalysisTaskSAQA&); // not implemented
80 ClassDef(AliAnalysisTaskSAQA, 17) // Quality task for Emcal analysis