]>
Commit | Line | Data |
---|---|---|
0531c107 | 1 | #ifndef ALIANALYSISTASKEMCALJETTRIGGERQA_H |
2 | #define ALIANALYSISTASKEMCALJETTRIGGERQA_H | |
3 | ||
4 | class TH1; | |
5 | class TH2; | |
6 | class TH3; | |
7 | class TH3F; | |
8 | class TClonesArray; | |
9 | class TArrayI; | |
0531c107 | 10 | |
11 | #include <TRef.h> | |
12 | #include <TBits.h> | |
13 | ||
9239b066 | 14 | #include "AliAnalysisTaskEmcalJet.h" |
0531c107 | 15 | |
9239b066 | 16 | class AliAnalysisTaskEmcalJetTriggerQA : public AliAnalysisTaskEmcalJet { |
0531c107 | 17 | public: |
18 | ||
19 | AliAnalysisTaskEmcalJetTriggerQA(); | |
20 | AliAnalysisTaskEmcalJetTriggerQA(const char *name); | |
21 | virtual ~AliAnalysisTaskEmcalJetTriggerQA(); | |
22 | ||
23 | void UserCreateOutputObjects(); | |
24 | void Terminate(Option_t *option); | |
25 | ||
0531c107 | 26 | Bool_t SelectEvent(); //decides if event is used for analysis |
27 | void FindTriggerPatch(); | |
28 | ||
0531c107 | 29 | //Setters |
30 | void SetDebug(Int_t d) { fDebug = d;} | |
31 | void SetTriggerClass(const char *n) { fTriggerClass = n; } | |
32 | void SetNFastorPatch(Int_t i) { fNFastOR = i;} | |
33 | ||
9e5eee5d | 34 | void SetContainerFull(Int_t c) { fContainerFull = c;} |
35 | void SetContainerCharged(Int_t c) { fContainerCharged = c;} | |
36 | ||
0714a353 | 37 | Int_t GetLeadingCellId(const AliVCluster *clus) const; |
38 | Double_t GetEnergyLeadingCell(const AliVCluster *clus) const; | |
698f34d8 | 39 | Double_t GetECross(Int_t absID) const; |
0531c107 | 40 | |
41 | Double_t GetZ(const AliVParticle *trk, const AliEmcalJet *jet) const; | |
42 | Double_t GetZ(const Double_t trkPx, const Double_t trkPy, const Double_t trkPz, const Double_t jetPx, const Double_t jetPy, const Double_t jetPz) const; | |
43 | ||
44 | protected: | |
45 | Bool_t FillHistograms() ; | |
46 | Bool_t Run() ; | |
47 | ||
48 | Bool_t TestFilterBit(Int_t trigBit, UInt_t bitJetTrig) const {return (Bool_t) ((trigBit & bitJetTrig) != 0);} | |
49 | ||
50 | ||
0531c107 | 51 | private: |
0714a353 | 52 | Bool_t fDebug; // debug level |
0714a353 | 53 | TString fTriggerClass; // trigger class to analyze EJ1 or EJ2 |
54 | UInt_t fBitJ1; // trigger bit of EJE1 | |
55 | UInt_t fBitJ2; // trigger bit of EJE2 | |
56 | Int_t fContainerFull; // number of container with full jets DET | |
57 | Int_t fContainerCharged; // number of container with charged jets DET | |
58 | Double_t fMaxPatchEnergy; // energy of patch with largest energy | |
59 | Int_t fTriggerType; // trigger type | |
60 | Int_t fNFastOR; // size of trigger patch fNFastORxfNFastOR | |
0531c107 | 61 | |
62 | TH1F *fhNEvents; //! Histo number of events | |
2fa783ab | 63 | TH2F *fHistRhovsCentFull; //!rho vs. centrality |
64 | TH2F *fHistRhovsCentCharged; //!rho vs. centrality | |
698f34d8 | 65 | TH3F *fh3PtEtaPhiTracks; //! pt,eta,phi of tracks |
a29bc208 CL |
66 | TH3F *fh3PtEtaPhiTracksOnEmcal; //! pt,eta,phi of tracks |
67 | TH3F *fh3PtEtaPhiTracksProp; //! pt,eta,phi of tracks | |
0531c107 | 68 | TH3F *fh3PtEtaPhiJetFull; //! pt,eta,phi of full jets |
69 | TH3F *fh3PtEtaPhiJetCharged; //! pt,eta,phi of charged jets | |
70 | TH2F *fh2NJetsPtFull; //! NJets per event vs pT,jet | |
71 | TH2F *fh2NJetsPtCharged; //! NJets per event vs pT,jet | |
72 | TH3F *fh3PtEtaAreaJetFull; //! pt,eta,area of full jet | |
73 | TH3F *fh3PtEtaAreaJetCharged; //! pt,eta,area of charged jets | |
74 | TH2F *fh2PtNConstituentsCharged; //! pt, # charged jet constituents | |
75 | TH2F *fh2PtNConstituents; //! pt, # jet constituents | |
76 | TH2F *fh2PtMeanPtConstituentsCharged; //! pt, <pt> charged constituents | |
77 | TH2F *fh2PtMeanPtConstituentsNeutral; //! pt, <pt> neutral constituents | |
78 | TH2F *fh2PtNEF; //! pt, NEF (neutral energy fraction) | |
3c709670 | 79 | TH3F *fh3NEFEtaPhi; //! NEF, eta, phi |
80 | TH2F *fh2NEFNConstituentsCharged; //! NEF, # charged jet constituents | |
81 | TH2F *fh2NEFNConstituentsNeutral; //! NEF, # neutral jet constituents | |
82 | TH2F *fh2Ptz; //! pt, z=pT,h,proj/p,jet full jet | |
83 | TH2F *fh2PtzCharged; //! pt, z=pT,h,proj/p,jet charged jet | |
0531c107 | 84 | TH2F *fh2PtLeadJet1VsLeadJet2; //! correlation between leading jet of the two branches |
85 | TH3F *fh3EEtaPhiCluster; //! cluster E, eta, phi | |
86 | TH3F *fh3PtLeadJet1VsPatchEnergy; //! leading jet energy vs leading patch energy vs jet trigger (J1/J2) | |
87 | TH3F *fh3PtLeadJet2VsPatchEnergy; //! leading jet energy vs leading patch energy vs jet trigger (J1/J2) | |
584095b1 | 88 | TH3F *fh3PatchEnergyEtaPhiCenterJ1; //! patch energy vs eta, phi at center of patch, high threshold |
89 | TH3F *fh3PatchEnergyEtaPhiCenterJ2; //! patch energy vs eta, phi at center of patch, low threshold | |
698f34d8 CL |
90 | TH3F *fh3PatchEnergyEtaPhiCenterJ1J2; //! patch energy vs eta, phi at center of patch, low + high threshold |
91 | TH3F *fh3PatchADCEnergyEtaPhiCenterJ1; //! patch ADC energy vs eta, phi at center of patch, high threshold | |
92 | TH3F *fh3PatchADCEnergyEtaPhiCenterJ2; //! patch ADC energy vs eta, phi at center of patch, low threshold | |
93 | TH3F *fh3PatchADCEnergyEtaPhiCenterJ1J2; //! patch ADC energy vs eta, phi at center of patch, low + high threshold | |
992189ac | 94 | TH2F *fh2CellEnergyVsTime; //! emcal cell energy vs time |
0714a353 | 95 | TH3F *fh3EClusELeadingCellVsTime; //! cluster energy vs energy of leading cell in cluster vs time of the leading cell |
0531c107 | 96 | |
97 | ||
98 | AliAnalysisTaskEmcalJetTriggerQA(const AliAnalysisTaskEmcalJetTriggerQA&); // not implemented | |
99 | AliAnalysisTaskEmcalJetTriggerQA &operator=(const AliAnalysisTaskEmcalJetTriggerQA&); // not implemented | |
100 | ||
2fa783ab | 101 | ClassDef(AliAnalysisTaskEmcalJetTriggerQA, 9) |
0531c107 | 102 | }; |
103 | #endif |