]>
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;} | |
38 | TH1F* GetJetJT(){return fhJetJT;} | |
39 | TH1F* GetPartonJT(){return fhPartonJT;} | |
40 | TH1F* GetJetPL(){return fhJetPL;} | |
41 | TH1F* GetPartonPL(){return fhPartonPL;} | |
42 | TH1F* GetJetEt(){return fhJetEt;} | |
43 | TH1F* GetJetEta(){return fhJetEta;} | |
44 | TH1F* GetPartonEta(){return fhPartonEta;} | |
45 | TH1F* GetPartonPhi(){return fhPartonPhi;} | |
46 | TH1F* GetJetPhi(){return fhJetPhi;} | |
47 | TH1F* GetEtaDiff(){return fhEtaDiff;} | |
48 | TH1F* GetPhiDiff(){return fhPhiDiff;} | |
49 | TH2F* GetEtaPhiSpread(){return fhEtaPhiSpread;} | |
50 | TH1F* GetNJets(){return fhNJets;} | |
73aa79d4 | 51 | |
52 | //========================== CASE 2 ======================== | |
53 | // Only consider events with at least 2 jets | |
ab01dff2 | 54 | TH1F* GetFragmFcn2(){return fhFragmFcn2;} |
55 | TH1F* GetPartonFragmFcn2(){return fhPartonFragmFcn2;} | |
56 | TH1F* GetJetJT2(){return fhJetJT2;} | |
57 | TH1F* GetPartonJT2(){return fhPartonJT2;} | |
58 | TH1F* GetJetPL2(){return fhJetPL2;} | |
59 | TH1F* GetPartonPL2(){return fhPartonPL2;} | |
60 | TH1F* GetJetEt2(){return fhJetEt2;} | |
61 | TH1F* GetJetEta2(){return fhJetEta2;} | |
62 | TH1F* GetPartonEta2(){return fhPartonEta2;} | |
63 | TH1F* GetPartonPhi2(){return fhPartonPhi2;} | |
64 | TH1F* GetJetPhi2(){return fhJetPhi2;} | |
65 | TH1F* GetEtaDiff2(){return fhEtaDiff2;} | |
66 | TH1F* GetPhiDiff2(){return fhPhiDiff2;} | |
67 | TH2F* GetEtaPhiSpread2(){return fhEtaPhiSpread2;} | |
68 | TH1F* GetNJets2(){return fhNJets2;} | |
69 | TH1F* GetJetEtSecond2(){return fhJetEtSecond2;} | |
70 | TH1F* GetJetEtRatio2(){return fhJetEtRatio2;} | |
71 | TH1F* GetEtaPhiDist2(){return fhEtaPhiDist2;} | |
73aa79d4 | 72 | |
f7d5860b | 73 | private: |
74 | void InitPlots(); | |
75 | Int_t fDebug; // Debug value | |
76 | Float_t fConeRadius; // Cone radius to be used in filling | |
77 | Float_t fNominalEnergy; // Force a nominal energy - specifically for 80+20 jets | |
78 | AliEMCALJetFinderOutput* fOutput; // Output object to be analysed | |
73aa79d4 | 79 | //===================== CASE 1 =========================================== |
ab01dff2 | 80 | TH1F *fhFragmFcn; // ("hFragmFcn","Fragmentation Function",100,0,1); |
81 | TH1F *fhPartonFragmFcn;// ("hFragmFcn","Parton Fragmentation Function",100,0,1); | |
82 | TH1F *fhPartonJT; // ("hPartonJT","Track Momentum Perpendicular to Parton Axis",100,0.,10.); | |
83 | TH1F *fhPartonPL; // ("hPartonPL","Track Momentum Parallel to Parton Axis ",100,0.,100.); | |
84 | TH1F *fhJetJT; // ("hJetJT","Track Momentum Perpendicular to Jet Axis",100,0.,10.); | |
85 | TH1F *fhJetPL; // ("hJetPL","Track Momentum Parallel to Jet Axis ",100,0.,100.); | |
86 | TH1F *fhJetEt; // ("hJetEt","E_{T}^{reco}",250,0.,250.); | |
87 | TH1F *fhJetEta; // ("hJetEta","#eta_{jet}^{reco}",180,-0.9,0.9); | |
88 | TH1F *fhJetPhi; // ("hJetPhi","#phi_{jet}^{reco}",62,0.,3.1); | |
89 | TH1F *fhPartonEta; // ("hPartonEta","#eta_{Parton}",180,-0.9,0.9); | |
90 | TH1F *fhPartonPhi; // ("hPartonPhi","#phi_{Parton}",62,0.,3.1); | |
91 | TH1F *fhEtaDiff; // ("hEtaDiff","#eta_{jet}^{reco}-#eta_{jet}^{input}",100,-0.5,0.5); | |
92 | TH1F *fhPhiDiff; // ("hPhiDiff","#phi_{jet}^{reco}-#phi_{jet}^{input}",100,-0.5,0.5); | |
93 | TH2F *fhEtaPhiSpread; // ("hEtaPhiSpread","#eta - #phi Distribution | |
f7d5860b | 94 | //of Reconstructed Jets",192,-0.7,0.7,288,pi/3,pi); |
ab01dff2 | 95 | TH1F *fhNJets; // ("hNJets","N Reconstructed jets",11,-0.5,10.5); |
73aa79d4 | 96 | |
97 | //============================== CASE 2 ============================================ | |
98 | ||
ab01dff2 | 99 | TH1F *fhFragmFcn2; // ("hFragmFcn2","Fragmentation Function",100,0,1); |
100 | TH1F *fhPartonFragmFcn2;// ("hFragmFcn2","Parton Fragmentation Function",100,0,1); | |
101 | TH1F *fhPartonJT2; // ("hPartonJT2","Track Momentum Perpendicular to Parton Axis",100,0.,10.); | |
102 | TH1F *fhPartonPL2; // ("hPartonPL2","Track Momentum Parallel to Parton Axis ",100,0.,100.); | |
103 | TH1F *fhJetJT2; // ("hJetJT2","Track Momentum Perpendicular to Jet Axis",100,0.,10.); | |
104 | TH1F *fhJetPL2; // ("hJetPL2","Track Momentum Parallel to Jet Axis ",100,0.,100.); | |
105 | TH1F *fhJetEt2; // ("hJetEt2","E_{T}^{reco}",250,0.,250.); | |
106 | TH1F *fhJetEta2; // ("hJetEta2","#eta_{jet}^{reco}",180,-0.9,0.9); | |
107 | TH1F *fhJetPhi2; // ("hJetPhi2","#phi_{jet}^{reco}",62,0.,3.1); | |
108 | TH1F *fhPartonEta2; // ("hPartonEta2","#eta_{Parton}",180,-0.9,0.9); | |
109 | TH1F *fhPartonPhi2; // ("hPartonPhi2","#phi_{Parton}",62,0.,3.1); | |
110 | TH1F *fhEtaDiff2; // ("hEtaDiff2","#eta_{jet}^{reco}-#eta_{jet}^{input}",100,-0.5,0.5); | |
111 | TH1F *fhPhiDiff2; // ("hPhiDiff2","#phi_{jet}^{reco}-#phi_{jet}^{input}",100,-0.5,0.5); | |
112 | TH2F *fhEtaPhiSpread2; // ("hEtaPhiSpread2","#eta - #phi Distribution | |
73aa79d4 | 113 | //of Reconstructed Jets",192,-0.7,0.7,288,pi/3,pi); |
ab01dff2 | 114 | TH1F *fhNJets2; // ("hNJets2","N Reconstructed jets",11,-0.5,10.5); |
115 | TH1F *fhJetEtSecond2; //("hJetEtSecond2","E_{T}^{reco}",250,0.,250.); | |
116 | TH1F *fhJetEtRatio2; //("hJetEtRatio2","Ratio of Second Highest to Highest",100,0,1); | |
117 | TH1F *fhEtaPhiDist2; //("hEtaPhiDist2","Angular Distance Between First and Second",100,0,3); | |
73aa79d4 | 118 | |
119 | ||
120 | ||
121 | ||
44f59d68 | 122 | Bool_t fInitialised; // have histograms been initialised |
f7d5860b | 123 | |
124 | ||
44f59d68 | 125 | ClassDef(AliEMCALJetFinderPlots,4) |
f7d5860b | 126 | |
127 | }; | |
128 | #endif | |
129 |