The component path in CDB changed - probably it was wrong before.
[u/mrichter/AliRoot.git] / PWG4 / JetTasks / AliAnalysisTaskJetSpectrum.h
CommitLineData
f3050824 1#ifndef ALIANALYSISTASKJETSPECTRUM_H
2#define ALIANALYSISTASKJETSPECTRUM_H
3
4/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
5 * See cxx source for full Copyright notice */
77bbf113 6
7// **************************************
8// Task used for the correction of determiantion of reconstructed jet spectra
9// Compares input (gen) and output (rec) jets
10// *******************************************
11
f3050824 12#include "AliAnalysisTaskSE.h"
77bbf113 13#include "THnSparse.h" // cannot forward declare ThnSparseF
14
df65bddb 15////////////////
f3050824 16class AliJetHeader;
17class AliESDEvent;
18class AliAODEvent;
19class AliAODJet;
20class AliGenPythiaEventHeader;
21
22class TList;
23class TChain;
24class TH2F;
25class TH3F;
4dbfdecc 26class TProfile;
27
f3050824 28
77bbf113 29
f3050824 30class AliAnalysisTaskJetSpectrum : public AliAnalysisTaskSE
31{
32 public:
33 AliAnalysisTaskJetSpectrum();
34 AliAnalysisTaskJetSpectrum(const char* name);
35 virtual ~AliAnalysisTaskJetSpectrum() {;}
36 // Implementation of interface methods
37 virtual void UserCreateOutputObjects();
38 virtual void Init();
39 virtual void LocalInit() { Init(); }
40 virtual void UserExec(Option_t *option);
41 virtual void Terminate(Option_t *option);
4dbfdecc 42 virtual Bool_t Notify();
43
f3050824 44
45 virtual void SetExternalWeight(Float_t f){fExternalWeight = f;}
46 virtual void SetUseExternalWeightOnly(Bool_t b){fUseExternalWeightOnly = b;}
47 virtual void SetAODInput(Bool_t b){fUseAODInput = b;}
48 virtual void SetLimitGenJetEta(Bool_t b){fLimitGenJetEta = b;}
49 virtual void SetAnalysisType(Int_t i){fAnalysisType = i;}
77bbf113 50 virtual void SetBranchGen(const char* c){fBranchGen = c;}
51 virtual void SetBranchRec(const char* c){fBranchRec = c;}
f3050824 52
53 // Helper
77bbf113 54 static void GetClosestJets(AliAODJet *genJets,const Int_t &nGenJets,
55 AliAODJet *recJets,const Int_t &nRecJets,
f3050824 56 Int_t *iGenIndex,Int_t *iRecIndex,Int_t iDebug = 0);
57
58 //
59
60 enum {kAnaMC = 0x1};
df65bddb 61 enum {kMaxJets = 5};
62 enum {kMaxCorrelation = 3};
f3050824 63
64 private:
65
66 AliAnalysisTaskJetSpectrum(const AliAnalysisTaskJetSpectrum&);
67 AliAnalysisTaskJetSpectrum& operator=(const AliAnalysisTaskJetSpectrum&);
68
df65bddb 69 static const Float_t fgkJetNpartCut[kMaxCorrelation];
f3050824 70
f3050824 71
72 AliJetHeader *fJetHeaderRec;
73 AliJetHeader *fJetHeaderGen;
74 AliAODEvent *fAOD; // where we take the jets from can be input or output AOD
75
76 TString fBranchRec; // AOD branch name for reconstructed
77 TString fConfigRec; // Name of the Config file
78 TString fBranchGen; // AOD brnach for genereated
79 TString fConfigGen; // Name of the Config file (if any)
80
81 Bool_t fUseAODInput;
82 Bool_t fUseExternalWeightOnly;
83 Bool_t fLimitGenJetEta;
84 Int_t fAnalysisType;
85 Float_t fExternalWeight;
86
4dbfdecc 87 TProfile* fh1Xsec; // pythia cross section and trials
88 TH1F* fh1Trials; // trials are added
f3050824 89 TH1F* fh1PtHard; // Pt har of the event...
77bbf113 90 TH1F* fh1PtHardNoW; // Pt har of the event...
91 TH1F* fh1PtHardTrials; // Number of trials
f3050824 92 TH1F* fh1NGenJets;
93 TH1F* fh1NRecJets;
94 TH1F* fh1E[kMaxJets]; // Jet Energy
95 TH1F* fh1PtRecIn[kMaxJets]; // Jet pt for all
96 TH1F* fh1PtRecOut[kMaxJets]; // Jet pt with corellated generated jet
97 TH1F* fh1PtGenIn[kMaxJets]; // Detection efficiency for given p_T.gen
77bbf113 98 TH1F* fh1PtGenOut[kMaxJets]; // gen pT of found jets
f3050824 99
100
101
77bbf113 102 TH2F* fh2PtFGen[kMaxJets]; // correlation betwen genreated and found jet pT
103 TH2F* fh2PhiFGen[kMaxJets]; // correlation betwen genreated and found jet phi
104 TH2F* fh2EtaFGen[kMaxJets]; // correlation betwen genreated and found jet eta
f3050824 105 TH2F* fh2Frag[kMaxJets]; // fragmentation function
77bbf113 106 TH2F* fh2FragLn[kMaxJets]; // fragmetation in xi
107 TH2F* fh2PtGenDeltaPhi[kMaxJets]; // difference between generated and found jet phi
108 TH2F* fh2PtGenDeltaEta[kMaxJets]; // difference between generated and found jet eta
109
110 TH3F* fh3PtRecGenHard[kMaxJets]; // correlation beetwen pt hard, rec and gen
111 TH3F* fh3PtRecGenHardNoW[kMaxJets]; // correlation beetwen pt hard, rec and gen no weight
112 TH3F* fh3RecEtaPhiPt[kMaxJets]; // correlation between eta phi and rec pt
113 TH3F* fh3RecEtaPhiPtNoGen[kMaxJets]; // correlation between eta phi and rec pt of jets without a partner
114 TH3F* fh3GenEtaPhiPtNoFound[kMaxJets]; // correlation between eta phi and gen pt of jets without a partner
115 TH3F* fh3GenEtaPhiPt[kMaxJets]; // correlation between eta phi and gen pt of jets without a partner
f3050824 116 // ========= Multiplicity dependence ======
117
118 // ==========TODO , flavaor dependence ========
119 // ============================================
120
121
122 // ============= TODO , phi dependence ========
123 // ============================================
124
125 TList *fHistList; // Output list
df65bddb 126
127 ///////// For 2 dimensional unfolding //////////////////
77bbf113 128 TH1F* fh1JetMultiplicity; // JetMultiplicity
129 TH2F* fh2ERecZRec; // rec e and z
130 TH2F* fh2EGenZGen; // gen e and z
131 TH2F* fh2Efficiency; // 2 dimensional efficiency
132 TH3F* fh3EGenERecN; // correlation rec energy, gen energy N particles
133 THnSparseF* fhnCorrelation[kMaxCorrelation]; // 4d Correllation for unfolding
df65bddb 134 ///////////////////////////////////////////////////////
f3050824 135
136
137 ClassDef(AliAnalysisTaskJetSpectrum, 1) // Analysis task for standard jet analysis
138};
139
140#endif