]> git.uio.no Git - u/mrichter/AliRoot.git/blob - PWGGA/EMCALJetTasks/AliAnalysisTaskEmcalJetSpectra.h
adapted to one-to-many arrays from Salvatore. Needs checks.
[u/mrichter/AliRoot.git] / PWGGA / EMCALJetTasks / AliAnalysisTaskEmcalJetSpectra.h
1 #ifndef AliAnalysisTaskEmcalJetSpectra_cxx
2 #define AliAnalysisTaskEmcalJetSpectra_cxx
3
4
5 class TList;
6 class TH1F;
7 class TH2F;
8 class AliESDEvent;
9
10 #include "AliAnalysisTaskSE.h"
11
12 class AliAnalysisTaskEmcalJetSpectra : public AliAnalysisTaskSE {
13  public:
14   AliAnalysisTaskEmcalJetSpectra();
15   AliAnalysisTaskEmcalJetSpectra(const char *name);
16   virtual ~AliAnalysisTaskEmcalJetSpectra() {}
17   
18   virtual void   UserCreateOutputObjects();
19   virtual void   UserExec(Option_t *option);
20   virtual void   Terminate(Option_t *);
21   virtual void   SetTracksName(const char *n)             { fTracksName = n;     }
22   virtual void   SetJetsName(const char *n)               { fJetsName =   n;     }
23   virtual void   SetRhos1Name(const char *rn)             { fRhos1Name=   rn;    }
24   virtual void   SetRhos2Name(const char *rn)             { fRhos2Name=   rn;    }
25   virtual void   SetRhos3Name(const char *rn)             { fRhos3Name=   rn;    }
26   virtual void   SetJetPhi(Double_t pmin, Double_t pmax)  { phimin =      pmin;  phimax = pmax; }
27   virtual void   SetJetEta(Double_t emin, Double_t emax)  { etamin =      emin;  etamax = emax; }
28   virtual void   SetAreaCut(Double_t a)                   { areacut =     a;     }
29   
30  protected:
31   virtual Int_t GetCentBin(Double_t cent) const;
32  
33   
34  private:
35    AliESDEvent *fESD;    //! ESD object
36   TList       *fOutputList; //! Output list
37   TH1F        *fHistCentrality;
38   TH1F        *fHistJetArea; 
39   TH1F        *fHistJetMaxPt;
40   TH1F        *fHistJetZ;
41   TH1F        *fHistJetNEF;
42   TH2F        *fHistJetPtvsCent;  
43   TH2F        *fHistJetPtM3vsCent;  
44   TH2F        *fHistLeadingJetPtvsCent;  
45   TH2F        *fHistLeadingJetPtM3vsCent;  
46   TH2F        *fHistJetAreavsCent;  
47   TH2F        *fHistJetMaxPtvsCent;  
48   TH2F        *fHistJetZvsCent;  
49   TH2F        *fHistJetNEFvsCent;
50   TH2F        *fHistNjetvsCent;
51   
52   TH2F        *fHistJetPtvsNtrack;
53   TH2F        *fHistJetAreavsNtrack;
54   TH2F        *fHistJetMaxPtvsNtrack;
55   TH2F        *fHistJetZvsNtrack;
56   TH2F        *fHistJetNEFvsNtrack;
57   TH2F        *fHistNjetvsNtrack;
58
59   TH2F       *fHistDeltaRho12vsCent;
60   TH2F       *fHistDeltaRho13vsCent;
61   TH2F       *fHistDeltaRho23vsCent;
62
63   TH2F       *fHistDeltaJetPt12vsCent;
64   TH2F       *fHistDeltaJetPt13vsCent;
65   TH2F       *fHistDeltaJetPt23vsCent;
66   
67   TH2F       *fHistRho1vsCent;
68   TH2F       *fHistRho2vsCent;
69   TH2F       *fHistRho3vsCent;
70
71   TString                fTracksName;             // name of track collection
72   TString                fJetsName;             // name of jet collection
73   TString                fClustersName;             // name of clusters collection
74   TString                fRhos1Name;             // name of Rho1 array output
75   TString                fRhos2Name;             // name of Rho2 array output
76   TString                fRhos3Name;             // name of Rho2 array output
77
78   Double_t phimin;
79   Double_t phimax;
80   Double_t etamin;
81   Double_t etamax;
82   Double_t areacut;
83
84   TH2F       *fHistNEFvsPt[6][4];
85   TH2F       *fHistZvsPt[6][4];
86   TH2F       *fHistZchvsPt[6][4];
87   TH2F       *fHistZemvsPt[6][4];
88   TH2F       *fHistAreavsPt[6][4];
89   TH1F       *fHistJetPt[6][4];
90   TH1F       *fHistJetPt5[6][4];
91   TH1F       *fHistJetPt6[6][4];
92   TH1F       *fHistJetPt7[6][4];
93   TH1F       *fHistJetPt8[6][4];
94   TH2F       *fHistNconsvsPt[6][4];
95   TH1F       *fHistRawJetPt[6];
96   TH2F       *fHistAreavsRawPt[6];
97
98
99    AliAnalysisTaskEmcalJetSpectra(const AliAnalysisTaskEmcalJetSpectra&); // not implemented
100   AliAnalysisTaskEmcalJetSpectra& operator=(const AliAnalysisTaskEmcalJetSpectra&); // not implemented
101   
102   ClassDef(AliAnalysisTaskEmcalJetSpectra, 1); // example of analysis
103 };
104
105 #endif