1 #ifndef ALIANALYSISTASKEMCALJETTRIGGERQA_H
2 #define ALIANALYSISTASKEMCALJETTRIGGERQA_H
14 #include "AliAnalysisTaskEmcalJet.h"
16 class AliAnalysisTaskEmcalJetTriggerQA : public AliAnalysisTaskEmcalJet {
19 AliAnalysisTaskEmcalJetTriggerQA();
20 AliAnalysisTaskEmcalJetTriggerQA(const char *name);
21 virtual ~AliAnalysisTaskEmcalJetTriggerQA();
23 void UserCreateOutputObjects();
24 void Terminate(Option_t *option);
26 Bool_t SelectEvent(); //decides if event is used for analysis
27 void FindTriggerPatch();
30 void SetDebug(Int_t d) { fDebug = d;}
31 void SetTriggerClass(const char *n) { fTriggerClass = n; }
32 void SetNFastorPatch(Int_t i) { fNFastOR = i;}
34 void SetContainerFull(Int_t c) { fContainerFull = c;}
35 void SetContainerCharged(Int_t c) { fContainerCharged = c;}
37 Int_t GetLeadingCellId(const AliVCluster *clus) const;
38 Double_t GetEnergyLeadingCell(const AliVCluster *clus) const;
40 Double_t GetZ(const AliVParticle *trk, const AliEmcalJet *jet) const;
41 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;
44 Bool_t FillHistograms() ;
47 Bool_t TestFilterBit(Int_t trigBit, UInt_t bitJetTrig) const {return (Bool_t) ((trigBit & bitJetTrig) != 0);}
51 Bool_t fDebug; // debug level
52 TString fTriggerClass; // trigger class to analyze EJ1 or EJ2
53 UInt_t fBitJ1; // trigger bit of EJE1
54 UInt_t fBitJ2; // trigger bit of EJE2
55 Int_t fContainerFull; // number of container with full jets DET
56 Int_t fContainerCharged; // number of container with charged jets DET
57 Double_t fMaxPatchEnergy; // energy of patch with largest energy
58 Int_t fTriggerType; // trigger type
59 Int_t fNFastOR; // size of trigger patch fNFastORxfNFastOR
61 TH1F *fhNEvents; //! Histo number of events
62 TH3F *fh3PtEtaPhiJetFull; //! pt,eta,phi of full jets
63 TH3F *fh3PtEtaPhiJetCharged; //! pt,eta,phi of charged jets
64 TH2F *fh2NJetsPtFull; //! NJets per event vs pT,jet
65 TH2F *fh2NJetsPtCharged; //! NJets per event vs pT,jet
66 TH3F *fh3PtEtaAreaJetFull; //! pt,eta,area of full jet
67 TH3F *fh3PtEtaAreaJetCharged; //! pt,eta,area of charged jets
68 TH2F *fh2PtNConstituentsCharged; //! pt, # charged jet constituents
69 TH2F *fh2PtNConstituents; //! pt, # jet constituents
70 TH2F *fh2PtMeanPtConstituentsCharged; //! pt, <pt> charged constituents
71 TH2F *fh2PtMeanPtConstituentsNeutral; //! pt, <pt> neutral constituents
72 TH2F *fh2PtNEF; //! pt, NEF (neutral energy fraction)
73 TH3F *fh3NEFEtaPhi; //! NEF, eta, phi
74 TH2F *fh2NEFNConstituentsCharged; //! NEF, # charged jet constituents
75 TH2F *fh2NEFNConstituentsNeutral; //! NEF, # neutral jet constituents
76 TH2F *fh2Ptz; //! pt, z=pT,h,proj/p,jet full jet
77 TH2F *fh2PtzCharged; //! pt, z=pT,h,proj/p,jet charged jet
78 TH2F *fh2PtLeadJet1VsLeadJet2; //! correlation between leading jet of the two branches
79 TH3F *fh3EEtaPhiCluster; //! cluster E, eta, phi
80 TH3F *fh3PtLeadJet1VsPatchEnergy; //! leading jet energy vs leading patch energy vs jet trigger (J1/J2)
81 TH3F *fh3PtLeadJet2VsPatchEnergy; //! leading jet energy vs leading patch energy vs jet trigger (J1/J2)
82 TH3F *fh3PatchEnergyEtaPhiCenterJ1; //! patch energy vs eta, phi at center of patch, high threshold
83 TH3F *fh3PatchEnergyEtaPhiCenterJ2; //! patch energy vs eta, phi at center of patch, low threshold
84 TH2F *fh2CellEnergyVsTime; //! emcal cell energy vs time
85 TH3F *fh3EClusELeadingCellVsTime; //! cluster energy vs energy of leading cell in cluster vs time of the leading cell
88 AliAnalysisTaskEmcalJetTriggerQA(const AliAnalysisTaskEmcalJetTriggerQA&); // not implemented
89 AliAnalysisTaskEmcalJetTriggerQA &operator=(const AliAnalysisTaskEmcalJetTriggerQA&); // not implemented
91 ClassDef(AliAnalysisTaskEmcalJetTriggerQA, 5)