X-Git-Url: http://git.uio.no/git/?p=u%2Fmrichter%2FAliRoot.git;a=blobdiff_plain;f=PWG4%2FJetTasks%2FAliAnalysisTaskJetSpectrum.h;h=c5256f31f7352976e5e59924f9467e739e53b023;hp=86790d54ebf6cc8a3327238daec43c1f9d0bd97b;hb=d40dc1ba4e5766eba31b1c3e79ae440f13c64403;hpb=22ceb537b68b7f8654fc8ab204f2003b9ddb2036 diff --git a/PWG4/JetTasks/AliAnalysisTaskJetSpectrum.h b/PWG4/JetTasks/AliAnalysisTaskJetSpectrum.h index 86790d54ebf..c5256f31f73 100644 --- a/PWG4/JetTasks/AliAnalysisTaskJetSpectrum.h +++ b/PWG4/JetTasks/AliAnalysisTaskJetSpectrum.h @@ -3,9 +3,15 @@ /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * * See cxx source for full Copyright notice */ - + +// ************************************** +// Task used for the correction of determiantion of reconstructed jet spectra +// Compares input (gen) and output (rec) jets +// ******************************************* + #include "AliAnalysisTaskSE.h" -#include "THnSparse.h" +#include "THnSparse.h" // cannot forward declare ThnSparseF + //////////////// class AliJetHeader; class AliESDEvent; @@ -20,6 +26,7 @@ class TH3F; class TProfile; + class AliAnalysisTaskJetSpectrum : public AliAnalysisTaskSE { public: @@ -39,19 +46,20 @@ class AliAnalysisTaskJetSpectrum : public AliAnalysisTaskSE virtual void SetUseExternalWeightOnly(Bool_t b){fUseExternalWeightOnly = b;} virtual void SetAODInput(Bool_t b){fUseAODInput = b;} virtual void SetLimitGenJetEta(Bool_t b){fLimitGenJetEta = b;} + virtual void SetRecEtaWindow(Float_t f){fRecEtaWindow = f;} virtual void SetAnalysisType(Int_t i){fAnalysisType = i;} - virtual void SetBranchGen(char* c){fBranchGen = c;} - virtual void SetBranchRec(char* c){fBranchRec = c;} + virtual void SetBranchGen(const char* c){fBranchGen = c;} + virtual void SetBranchRec(const char* c){fBranchRec = c;} // Helper - static void GetClosestJets(AliAODJet *genJets,Int_t &nGenJets, - AliAODJet *recJets,Int_t &nRecJets, + static void GetClosestJets(AliAODJet *genJets,const Int_t &nGenJets, + AliAODJet *recJets,const Int_t &nRecJets, Int_t *iGenIndex,Int_t *iRecIndex,Int_t iDebug = 0); // enum {kAnaMC = 0x1}; - enum {kMaxJets = 5}; + enum {kMaxJets = 4}; enum {kMaxCorrelation = 3}; private: @@ -71,41 +79,42 @@ class AliAnalysisTaskJetSpectrum : public AliAnalysisTaskSE TString fBranchGen; // AOD brnach for genereated TString fConfigGen; // Name of the Config file (if any) - Bool_t fUseAODInput; - Bool_t fUseExternalWeightOnly; - Bool_t fLimitGenJetEta; - Int_t fAnalysisType; - Float_t fExternalWeight; + Bool_t fUseAODInput; // use AOD input + Bool_t fUseExternalWeightOnly; // use only external weight + Bool_t fLimitGenJetEta; // Limit the eta of the generated jets + Int_t fAnalysisType; // Analysis type + Float_t fExternalWeight; // external weight + Float_t fRecEtaWindow; // eta window used for corraltion plots between rec and gen TProfile* fh1Xsec; // pythia cross section and trials TH1F* fh1Trials; // trials are added TH1F* fh1PtHard; // Pt har of the event... - TH1F* fh1PtHard_NoW; // Pt har of the event... - TH1F* fh1PtHard_Trials; // Number of trials + TH1F* fh1PtHardNoW; // Pt har of the event... + TH1F* fh1PtHardTrials; // Number of trials TH1F* fh1NGenJets; TH1F* fh1NRecJets; TH1F* fh1E[kMaxJets]; // Jet Energy TH1F* fh1PtRecIn[kMaxJets]; // Jet pt for all TH1F* fh1PtRecOut[kMaxJets]; // Jet pt with corellated generated jet TH1F* fh1PtGenIn[kMaxJets]; // Detection efficiency for given p_T.gen - TH1F* fh1PtGenOut[kMaxJets]; // - + TH1F* fh1PtGenOut[kMaxJets]; // gen pT of found jets - - TH2F* fh2PtFGen[kMaxJets]; // - TH2F* fh2PhiFGen[kMaxJets]; // - TH2F* fh2EtaFGen[kMaxJets]; // + TH2F* fh2PtFGen[kMaxJets]; // correlation betwen genreated and found jet pT + TH2F* fh2PhiFGen[kMaxJets]; // correlation betwen genreated and found jet phi + TH2F* fh2EtaFGen[kMaxJets]; // correlation betwen genreated and found jet eta TH2F* fh2Frag[kMaxJets]; // fragmentation function - TH2F* fh2FragLn[kMaxJets]; // - TH2F* fh2PtGenDeltaPhi[kMaxJets]; - TH2F* fh2PtGenDeltaEta[kMaxJets]; - - TH3F* fh3PtRecGenHard[kMaxJets]; // - TH3F* fh3PtRecGenHard_NoW[kMaxJets]; // - TH3F* fh3RecEtaPhiPt[kMaxJets]; // - TH3F* fh3RecEtaPhiPt_NoGen[kMaxJets]; // - TH3F* fh3GenEtaPhiPt_NoFound[kMaxJets]; // - TH3F* fh3GenEtaPhiPt[kMaxJets]; // + TH2F* fh2FragLn[kMaxJets]; // fragmetation in xi + TH2F* fh2PtGenDeltaPhi[kMaxJets]; // difference between generated and found jet phi + TH2F* fh2PtGenDeltaEta[kMaxJets]; // difference between generated and found jet eta + TH2F* fh2PtRecDeltaR[kMaxJets]; // distance of rec jet i to j > i p_T,j + TH2F* fh2PtGenDeltaR[kMaxJets]; // distance of jet i to j > i vs p_T,j + + TH3F* fh3PtRecGenHard[kMaxJets]; // correlation beetwen pt hard, rec and gen + TH3F* fh3PtRecGenHardNoW[kMaxJets]; // correlation beetwen pt hard, rec and gen no weight + TH3F* fh3RecEtaPhiPt[kMaxJets]; // correlation between eta phi and rec pt + TH3F* fh3RecEtaPhiPtNoGen[kMaxJets]; // correlation between eta phi and rec pt of jets without a partner + TH3F* fh3GenEtaPhiPtNoFound[kMaxJets]; // correlation between eta phi and gen pt of jets without a partner + TH3F* fh3GenEtaPhiPt[kMaxJets]; // correlation between eta phi and gen pt of jets without a partner // ========= Multiplicity dependence ====== // ==========TODO , flavaor dependence ======== @@ -118,12 +127,12 @@ class AliAnalysisTaskJetSpectrum : public AliAnalysisTaskSE TList *fHistList; // Output list ///////// For 2 dimensional unfolding ////////////////// - TH1F* fh1JetMultiplicity; - TH2F* fh2ERecZRec; - TH2F* fh2EGenZGen; - TH2F* fh2Efficiency; - TH3F* fh3EGenERecN; - THnSparseF* fhnCorrelation[kMaxCorrelation]; + TH1F* fh1JetMultiplicity; // JetMultiplicity + TH2F* fh2ERecZRec; // rec e and z + TH2F* fh2EGenZGen; // gen e and z + TH2F* fh2Efficiency; // 2 dimensional efficiency + TH3F* fh3EGenERecN; // correlation rec energy, gen energy N particles + THnSparseF* fhnCorrelation[kMaxCorrelation]; // 4d Correllation for unfolding ///////////////////////////////////////////////////////