]>
Commit | Line | Data |
---|---|---|
f7d5860b | 1 | #ifndef ALIEMCALJETFINDERPLOTS_H |
2 | #define ALIEMCALJETFINDERPLOTS_H | |
3 | ||
4 | /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * | |
5 | * * * See cxx source for full Copyright notice */ | |
6 | ||
73aa79d4 | 7 | |
f7d5860b | 8 | /* $Id$ */ |
9 | ||
10 | //_________________________________________________________________________ | |
11 | // Class for Filling jetfinder plots | |
12 | // | |
13 | //*-- Author: Mark Horner (LBL/UCT) | |
14 | // | |
15 | // | |
16 | ||
17 | ||
18 | ||
19 | #include "TObject.h" | |
20 | #include "TH1F.h" | |
21 | #include "TH2F.h" | |
22 | ||
23 | #include "AliEMCALJetFinderOutput.h" | |
24 | ||
25 | class AliEMCALJetFinderPlots : public TObject | |
26 | { | |
27 | public: | |
28 | AliEMCALJetFinderPlots(); | |
29 | ~AliEMCALJetFinderPlots(); | |
30 | void SetConeRadius(Float_t coneradius){fConeRadius = coneradius;} | |
31 | void SetNominalEnergy(Float_t energy){fNominalEnergy = energy;} | |
32 | void SetDebug(Int_t debug){fDebug = debug;} | |
33 | void FillFromOutput(AliEMCALJetFinderOutput* output); | |
73aa79d4 | 34 | //========================== CASE 1 ======================== |
35 | // Only consider events with only 1 jet | |
ab01dff2 | 36 | TH1F* GetFragmFcn(){return fhFragmFcn;} |
37 | TH1F* GetPartonFragmFcn(){return fhPartonFragmFcn;} | |
1515a8af | 38 | TH1F* GetPT(){return fhJetPT;} |
39 | TH1F* GetPartonPT(){return fhPartonPT;} | |
ab01dff2 | 40 | TH1F* GetJetJT(){return fhJetJT;} |
41 | TH1F* GetPartonJT(){return fhPartonJT;} | |
42 | TH1F* GetJetPL(){return fhJetPL;} | |
43 | TH1F* GetPartonPL(){return fhPartonPL;} | |
44 | TH1F* GetJetEt(){return fhJetEt;} | |
45 | TH1F* GetJetEta(){return fhJetEta;} | |
46 | TH1F* GetPartonEta(){return fhPartonEta;} | |
47 | TH1F* GetPartonPhi(){return fhPartonPhi;} | |
48 | TH1F* GetJetPhi(){return fhJetPhi;} | |
49 | TH1F* GetEtaDiff(){return fhEtaDiff;} | |
50 | TH1F* GetPhiDiff(){return fhPhiDiff;} | |
51 | TH2F* GetEtaPhiSpread(){return fhEtaPhiSpread;} | |
52 | TH1F* GetNJets(){return fhNJets;} | |
73aa79d4 | 53 | |
54 | //========================== CASE 2 ======================== | |
55 | // Only consider events with at least 2 jets | |
ab01dff2 | 56 | TH1F* GetFragmFcn2(){return fhFragmFcn2;} |
57 | TH1F* GetPartonFragmFcn2(){return fhPartonFragmFcn2;} | |
1515a8af | 58 | TH1F* GetPT2(){return fhJetPT2;} |
59 | TH1F* GetPartonPT2(){return fhPartonPT2;} | |
ab01dff2 | 60 | TH1F* GetJetJT2(){return fhJetJT2;} |
61 | TH1F* GetPartonJT2(){return fhPartonJT2;} | |
62 | TH1F* GetJetPL2(){return fhJetPL2;} | |
63 | TH1F* GetPartonPL2(){return fhPartonPL2;} | |
64 | TH1F* GetJetEt2(){return fhJetEt2;} | |
65 | TH1F* GetJetEta2(){return fhJetEta2;} | |
66 | TH1F* GetPartonEta2(){return fhPartonEta2;} | |
67 | TH1F* GetPartonPhi2(){return fhPartonPhi2;} | |
68 | TH1F* GetJetPhi2(){return fhJetPhi2;} | |
69 | TH1F* GetEtaDiff2(){return fhEtaDiff2;} | |
70 | TH1F* GetPhiDiff2(){return fhPhiDiff2;} | |
71 | TH2F* GetEtaPhiSpread2(){return fhEtaPhiSpread2;} | |
72 | TH1F* GetNJets2(){return fhNJets2;} | |
73 | TH1F* GetJetEtSecond2(){return fhJetEtSecond2;} | |
74 | TH1F* GetJetEtRatio2(){return fhJetEtRatio2;} | |
75 | TH1F* GetEtaPhiDist2(){return fhEtaPhiDist2;} | |
73aa79d4 | 76 | |
492f773f | 77 | |
78 | TH1F* GetJetPt(){return fhJetPT ;} | |
79 | TH1F* GetPartonPt(){return fhPartonPT ;} | |
80 | TH1F* GetJetPt2(){return fhJetPT2;} | |
81 | TH1F* GetPartonPt2(){return fhPartonPT2;} | |
82 | TH1F* GetRecoBinFragmFcn() {return fhRecoBinFragmFcn;} | |
83 | TH1F* GetRecoBinPartonFragmFcn() {return fhRecoBinPartonFragmFcn;} | |
84 | ||
63131144 | 85 | //============================== ALL CASES ============================================ |
86 | ||
87 | TH2F* GetInputOutput(){return fhInputOutput;} | |
88 | ||
89 | //============================== Reconstruction Bin Comparison ============================================ | |
90 | ||
91 | TH1F* GetRecoBinPt(){return fhRecoBinPt;} // ("fhRecoBinPt","Reconstructed Pt Distribution",100,0,1); | |
92 | TH1F* GetRecoBinPartonPt(){return fhRecoBinPartonPt;} // ("fhRecoBinPartonPt","Input Pt Distribution",100,0,1); | |
93 | TH1F* GetRecoBinJetEt(){return fhRecoBinJetEt;} // ("fhRecoJetEt","E_{T}^{reco}",250,0.,250.); | |
94 | TH1F* GetRecoBinInputJetEt(){return fhRecoBinInputJetEt;} // ("fhRecoInputJetEt","E_{T}^{reco}",250,0.,250.); | |
95 | ||
f7d5860b | 96 | private: |
97 | void InitPlots(); | |
98 | Int_t fDebug; // Debug value | |
99 | Float_t fConeRadius; // Cone radius to be used in filling | |
100 | Float_t fNominalEnergy; // Force a nominal energy - specifically for 80+20 jets | |
101 | AliEMCALJetFinderOutput* fOutput; // Output object to be analysed | |
73aa79d4 | 102 | //===================== CASE 1 =========================================== |
ab01dff2 | 103 | TH1F *fhFragmFcn; // ("hFragmFcn","Fragmentation Function",100,0,1); |
104 | TH1F *fhPartonFragmFcn;// ("hFragmFcn","Parton Fragmentation Function",100,0,1); | |
105 | TH1F *fhPartonJT; // ("hPartonJT","Track Momentum Perpendicular to Parton Axis",100,0.,10.); | |
106 | TH1F *fhPartonPL; // ("hPartonPL","Track Momentum Parallel to Parton Axis ",100,0.,100.); | |
107 | TH1F *fhJetJT; // ("hJetJT","Track Momentum Perpendicular to Jet Axis",100,0.,10.); | |
108 | TH1F *fhJetPL; // ("hJetPL","Track Momentum Parallel to Jet Axis ",100,0.,100.); | |
109 | TH1F *fhJetEt; // ("hJetEt","E_{T}^{reco}",250,0.,250.); | |
110 | TH1F *fhJetEta; // ("hJetEta","#eta_{jet}^{reco}",180,-0.9,0.9); | |
111 | TH1F *fhJetPhi; // ("hJetPhi","#phi_{jet}^{reco}",62,0.,3.1); | |
112 | TH1F *fhPartonEta; // ("hPartonEta","#eta_{Parton}",180,-0.9,0.9); | |
113 | TH1F *fhPartonPhi; // ("hPartonPhi","#phi_{Parton}",62,0.,3.1); | |
114 | TH1F *fhEtaDiff; // ("hEtaDiff","#eta_{jet}^{reco}-#eta_{jet}^{input}",100,-0.5,0.5); | |
115 | TH1F *fhPhiDiff; // ("hPhiDiff","#phi_{jet}^{reco}-#phi_{jet}^{input}",100,-0.5,0.5); | |
116 | TH2F *fhEtaPhiSpread; // ("hEtaPhiSpread","#eta - #phi Distribution | |
f7d5860b | 117 | //of Reconstructed Jets",192,-0.7,0.7,288,pi/3,pi); |
ab01dff2 | 118 | TH1F *fhNJets; // ("hNJets","N Reconstructed jets",11,-0.5,10.5); |
73aa79d4 | 119 | |
120 | //============================== CASE 2 ============================================ | |
121 | ||
ab01dff2 | 122 | TH1F *fhFragmFcn2; // ("hFragmFcn2","Fragmentation Function",100,0,1); |
123 | TH1F *fhPartonFragmFcn2;// ("hFragmFcn2","Parton Fragmentation Function",100,0,1); | |
124 | TH1F *fhPartonJT2; // ("hPartonJT2","Track Momentum Perpendicular to Parton Axis",100,0.,10.); | |
125 | TH1F *fhPartonPL2; // ("hPartonPL2","Track Momentum Parallel to Parton Axis ",100,0.,100.); | |
126 | TH1F *fhJetJT2; // ("hJetJT2","Track Momentum Perpendicular to Jet Axis",100,0.,10.); | |
127 | TH1F *fhJetPL2; // ("hJetPL2","Track Momentum Parallel to Jet Axis ",100,0.,100.); | |
128 | TH1F *fhJetEt2; // ("hJetEt2","E_{T}^{reco}",250,0.,250.); | |
129 | TH1F *fhJetEta2; // ("hJetEta2","#eta_{jet}^{reco}",180,-0.9,0.9); | |
130 | TH1F *fhJetPhi2; // ("hJetPhi2","#phi_{jet}^{reco}",62,0.,3.1); | |
131 | TH1F *fhPartonEta2; // ("hPartonEta2","#eta_{Parton}",180,-0.9,0.9); | |
132 | TH1F *fhPartonPhi2; // ("hPartonPhi2","#phi_{Parton}",62,0.,3.1); | |
133 | TH1F *fhEtaDiff2; // ("hEtaDiff2","#eta_{jet}^{reco}-#eta_{jet}^{input}",100,-0.5,0.5); | |
134 | TH1F *fhPhiDiff2; // ("hPhiDiff2","#phi_{jet}^{reco}-#phi_{jet}^{input}",100,-0.5,0.5); | |
135 | TH2F *fhEtaPhiSpread2; // ("hEtaPhiSpread2","#eta - #phi Distribution | |
73aa79d4 | 136 | //of Reconstructed Jets",192,-0.7,0.7,288,pi/3,pi); |
ab01dff2 | 137 | TH1F *fhNJets2; // ("hNJets2","N Reconstructed jets",11,-0.5,10.5); |
138 | TH1F *fhJetEtSecond2; //("hJetEtSecond2","E_{T}^{reco}",250,0.,250.); | |
139 | TH1F *fhJetEtRatio2; //("hJetEtRatio2","Ratio of Second Highest to Highest",100,0,1); | |
140 | TH1F *fhEtaPhiDist2; //("hEtaPhiDist2","Angular Distance Between First and Second",100,0,3); | |
73aa79d4 | 141 | |
63131144 | 142 | //============================== ALL CASES ============================================ |
73aa79d4 | 143 | |
63131144 | 144 | TH2F *fhInputOutput; //("hJetEtRatio2","Ratio of Second Highest to Highest",100,0,1); |
73aa79d4 | 145 | |
63131144 | 146 | //============================== Reconstruction Bin Comparison ============================================ |
147 | ||
148 | TH1F *fhRecoBinPt; // ("fhRecoBinPt","Reconstructed Pt Distribution",100,0,1); | |
149 | TH1F *fhRecoBinPartonPt; // ("fhRecoBinPartonPt","Input Pt Distribution",100,0,1); | |
150 | TH1F *fhRecoBinJetEt; // ("fhRecoJetEt","E_{T}^{reco}",250,0.,250.); | |
151 | TH1F *fhRecoBinInputJetEt; // ("fhRecoInputJetEt","E_{T}^{reco}",250,0.,250.); | |
1515a8af | 152 | TH1F* fhJetPT ;// new TH1F("hJetPT","P_{T} Distribution",200,0,200); |
153 | TH1F* fhPartonPT ;// new TH1F("hPartonPT","Parton P_{T} Distribution",200,0,1); | |
154 | TH1F* fhJetPT2 ;// new TH1F("hJetPT","P_{T} Distribution",200,0,200); | |
155 | TH1F* fhPartonPT2 ;// new TH1F("hPartonPT","Parton P_{T} Distribution",200,0,1); | |
156 | TH1F* fhRecoBinFragmFcn;//new TH1F("fhRecoBinFragmFcn","Reconstructed Frag. Fcn",100,0,1); | |
157 | TH1F* fhRecoBinPartonFragmFcn;// new TH1F("fhRecoBinPartonFragmFcn","Input Bin Fragm Fcn Distribution",100,0,1); | |
da2f0870 | 158 | |
159 | TH1F* fhJetInvE;// new TH1F("fhJetInvE","#frac{1}{E_{R}}",100,0,1); | |
160 | TH1F* fhJetInvE2;// new TH1F("fhJetInvE2","#frac{1}{E_{R}}",100,0,1); | |
63131144 | 161 | |
44f59d68 | 162 | Bool_t fInitialised; // have histograms been initialised |
f7d5860b | 163 | |
164 | ||
492f773f | 165 | ClassDef(AliEMCALJetFinderPlots,6) |
f7d5860b | 166 | |
167 | }; | |
168 | #endif | |
169 |