From a9e585a792827b46341ba0190fa4921450179132 Mon Sep 17 00:00:00 2001 From: kleinb Date: Thu, 16 Feb 2012 12:27:33 +0000 Subject: [PATCH] changed binning, vertex range, leading, subleading --- PWGJE/AliAnalysisTaskJetCore.cxx | 912 ++++++++++++++++++++++++++++--- PWGJE/AliAnalysisTaskJetCore.h | 148 ++++- PWGJE/macros/AddTaskJetCore.C | 6 +- 3 files changed, 989 insertions(+), 77 deletions(-) diff --git a/PWGJE/AliAnalysisTaskJetCore.cxx b/PWGJE/AliAnalysisTaskJetCore.cxx index edc2451c42b..ef551ab170d 100644 --- a/PWGJE/AliAnalysisTaskJetCore.cxx +++ b/PWGJE/AliAnalysisTaskJetCore.cxx @@ -1,3 +1,4 @@ + // ****************************************** // This task computes several jet observables like // the fraction of energy in inner and outer coronnas, @@ -26,7 +27,6 @@ #include "TChain.h" #include "TTree.h" #include "TMath.h" -#include "TH1I.h" #include "TH1F.h" #include "TH2F.h" #include "TH3F.h" @@ -65,8 +65,8 @@ fNonStdFile(""), fIsPbPb(kTRUE), fOfflineTrgMask(AliVEvent::kAny), fMinContribVtx(1), -fVtxZMin(-8.), -fVtxZMax(8.), +fVtxZMin(-10.), +fVtxZMax(10.), fEvtClassMin(0), fEvtClassMax(4), fRadioFrac(0.2), @@ -103,23 +103,62 @@ fh2JetCoreMethod3C30(0x0), fh2JetCoreMethod1C60(0x0), fh2JetCoreMethod2C60(0x0), fh2JetCoreMethod3C60(0x0), +fh2JetCoreMethod3C10lead(0x0), +fh2JetCoreMethod3C20lead(0x0), +fh2JetCoreMethod3C30lead(0x0), +fh2JetCoreMethod3C60lead(0x0), +fh2JetCoreMethod3C10sublead(0x0), +fh2JetCoreMethod3C20sublead(0x0), +fh2JetCoreMethod3C30sublead(0x0), +fh2JetCoreMethod3C60sublead(0x0), fh2SumPtInC10(0x0), fh2SumPtInC20(0x0), fh2SumPtInC30(0x0), fh2SumPtInC60(0x0), +fh2SumPtInC10lead(0x0), +fh2SumPtInC20lead(0x0), +fh2SumPtInC30lead(0x0), +fh2SumPtInC60lead(0x0), +fh2SumPtInC10sublead(0x0), +fh2SumPtInC20sublead(0x0), +fh2SumPtInC30sublead(0x0), +fh2SumPtInC60sublead(0x0), fh2SumPtOutC10(0x0), -fh2SumPtOutC10b(0x0), fh2SumPtOutC20(0x0), fh2SumPtOutC30(0x0), fh2SumPtOutC60(0x0), +fh2SumPtOutC10lead(0x0), +fh2SumPtOutC20lead(0x0), +fh2SumPtOutC30lead(0x0), +fh2SumPtOutC60lead(0x0), +fh2SumPtOutC10sublead(0x0), +fh2SumPtOutC20sublead(0x0), +fh2SumPtOutC30sublead(0x0), +fh2SumPtOutC60sublead(0x0), fh2SumPtInC10bkg(0x0), fh2SumPtInC20bkg(0x0), fh2SumPtInC30bkg(0x0), fh2SumPtInC60bkg(0x0), +fh2SumPtInC10bkglead(0x0), +fh2SumPtInC20bkglead(0x0), +fh2SumPtInC30bkglead(0x0), +fh2SumPtInC60bkglead(0x0), +fh2SumPtInC10bkgsublead(0x0), +fh2SumPtInC20bkgsublead(0x0), +fh2SumPtInC30bkgsublead(0x0), +fh2SumPtInC60bkgsublead(0x0), fh2SumPtOutC10bkg(0x0), fh2SumPtOutC20bkg(0x0), fh2SumPtOutC30bkg(0x0), fh2SumPtOutC60bkg(0x0), +fh2SumPtOutC10bkglead(0x0), +fh2SumPtOutC20bkglead(0x0), +fh2SumPtOutC30bkglead(0x0), +fh2SumPtOutC60bkglead(0x0), +fh2SumPtOutC10bkgsublead(0x0), +fh2SumPtOutC20bkgsublead(0x0), +fh2SumPtOutC30bkgsublead(0x0), +fh2SumPtOutC60bkgsublead(0x0), fh2DeltaRC10pt1(0x0), fh2DeltaRC20pt1(0x0), fh2DeltaRC30pt1(0x0), @@ -168,6 +207,102 @@ fh2DeltaPhiC10pt4(0x0), fh2DeltaPhiC20pt4(0x0), fh2DeltaPhiC30pt4(0x0), fh2DeltaPhiC60pt4(0x0), +fh2DeltaRC10pt1lead(0x0), +fh2DeltaRC20pt1lead(0x0), +fh2DeltaRC30pt1lead(0x0), +fh2DeltaRC60pt1lead(0x0), +fh2DeltaRC10pt2lead(0x0), +fh2DeltaRC20pt2lead(0x0), +fh2DeltaRC30pt2lead(0x0), +fh2DeltaRC60pt2lead(0x0), +fh2DeltaRC10pt3lead(0x0), +fh2DeltaRC20pt3lead(0x0), +fh2DeltaRC30pt3lead(0x0), +fh2DeltaRC60pt3lead(0x0), +fh2DeltaRC10pt4lead(0x0), +fh2DeltaRC20pt4lead(0x0), +fh2DeltaRC30pt4lead(0x0), +fh2DeltaRC60pt4lead(0x0), +fh2DeltaEtaC10pt1lead(0x0), +fh2DeltaEtaC20pt1lead(0x0), +fh2DeltaEtaC30pt1lead(0x0), +fh2DeltaEtaC60pt1lead(0x0), +fh2DeltaEtaC10pt2lead(0x0), +fh2DeltaEtaC20pt2lead(0x0), +fh2DeltaEtaC30pt2lead(0x0), +fh2DeltaEtaC60pt2lead(0x0), +fh2DeltaEtaC10pt3lead(0x0), +fh2DeltaEtaC20pt3lead(0x0), +fh2DeltaEtaC30pt3lead(0x0), +fh2DeltaEtaC60pt3lead(0x0), +fh2DeltaEtaC10pt4lead(0x0), +fh2DeltaEtaC20pt4lead(0x0), +fh2DeltaEtaC30pt4lead(0x0), +fh2DeltaEtaC60pt4lead(0x0), +fh2DeltaPhiC10pt1lead(0x0), +fh2DeltaPhiC20pt1lead(0x0), +fh2DeltaPhiC30pt1lead(0x0), +fh2DeltaPhiC60pt1lead(0x0), +fh2DeltaPhiC10pt2lead(0x0), +fh2DeltaPhiC20pt2lead(0x0), +fh2DeltaPhiC30pt2lead(0x0), +fh2DeltaPhiC60pt2lead(0x0), +fh2DeltaPhiC10pt3lead(0x0), +fh2DeltaPhiC20pt3lead(0x0), +fh2DeltaPhiC30pt3lead(0x0), +fh2DeltaPhiC60pt3lead(0x0), +fh2DeltaPhiC10pt4lead(0x0), +fh2DeltaPhiC20pt4lead(0x0), +fh2DeltaPhiC30pt4lead(0x0), +fh2DeltaPhiC60pt4lead(0x0), +fh2DeltaRC10pt1sublead(0x0), +fh2DeltaRC20pt1sublead(0x0), +fh2DeltaRC30pt1sublead(0x0), +fh2DeltaRC60pt1sublead(0x0), +fh2DeltaRC10pt2sublead(0x0), +fh2DeltaRC20pt2sublead(0x0), +fh2DeltaRC30pt2sublead(0x0), +fh2DeltaRC60pt2sublead(0x0), +fh2DeltaRC10pt3sublead(0x0), +fh2DeltaRC20pt3sublead(0x0), +fh2DeltaRC30pt3sublead(0x0), +fh2DeltaRC60pt3sublead(0x0), +fh2DeltaRC10pt4sublead(0x0), +fh2DeltaRC20pt4sublead(0x0), +fh2DeltaRC30pt4sublead(0x0), +fh2DeltaRC60pt4sublead(0x0), +fh2DeltaEtaC10pt1sublead(0x0), +fh2DeltaEtaC20pt1sublead(0x0), +fh2DeltaEtaC30pt1sublead(0x0), +fh2DeltaEtaC60pt1sublead(0x0), +fh2DeltaEtaC10pt2sublead(0x0), +fh2DeltaEtaC20pt2sublead(0x0), +fh2DeltaEtaC30pt2sublead(0x0), +fh2DeltaEtaC60pt2sublead(0x0), +fh2DeltaEtaC10pt3sublead(0x0), +fh2DeltaEtaC20pt3sublead(0x0), +fh2DeltaEtaC30pt3sublead(0x0), +fh2DeltaEtaC60pt3sublead(0x0), +fh2DeltaEtaC10pt4sublead(0x0), +fh2DeltaEtaC20pt4sublead(0x0), +fh2DeltaEtaC30pt4sublead(0x0), +fh2DeltaEtaC60pt4sublead(0x0), +fh2DeltaPhiC10pt1sublead(0x0), +fh2DeltaPhiC20pt1sublead(0x0), +fh2DeltaPhiC30pt1sublead(0x0), +fh2DeltaPhiC60pt1sublead(0x0), +fh2DeltaPhiC10pt2sublead(0x0), +fh2DeltaPhiC20pt2sublead(0x0), +fh2DeltaPhiC30pt2sublead(0x0), +fh2DeltaPhiC60pt2sublead(0x0), +fh2DeltaPhiC10pt3sublead(0x0), +fh2DeltaPhiC20pt3sublead(0x0), +fh2DeltaPhiC30pt3sublead(0x0), +fh2DeltaPhiC60pt3sublead(0x0), +fh2DeltaPhiC10pt4sublead(0x0), +fh2DeltaPhiC20pt4sublead(0x0), +fh2DeltaPhiC30pt4sublead(0x0), +fh2DeltaPhiC60pt4sublead(0x0), fh2AngStructpt1C10(0x0), fh2AngStructpt2C10(0x0), fh2AngStructpt3C10(0x0), @@ -204,8 +339,8 @@ fNonStdFile(""), fIsPbPb(kTRUE), fOfflineTrgMask(AliVEvent::kAny), fMinContribVtx(1), -fVtxZMin(-8.), -fVtxZMax(8.), +fVtxZMin(-10.), +fVtxZMax(10.), fEvtClassMin(0), fEvtClassMax(4), fRadioFrac(0.2), @@ -242,23 +377,62 @@ fh2JetCoreMethod3C30(0x0), fh2JetCoreMethod1C60(0x0), fh2JetCoreMethod2C60(0x0), fh2JetCoreMethod3C60(0x0), +fh2JetCoreMethod3C10lead(0x0), +fh2JetCoreMethod3C20lead(0x0), +fh2JetCoreMethod3C30lead(0x0), +fh2JetCoreMethod3C60lead(0x0), +fh2JetCoreMethod3C10sublead(0x0), +fh2JetCoreMethod3C20sublead(0x0), +fh2JetCoreMethod3C30sublead(0x0), +fh2JetCoreMethod3C60sublead(0x0), fh2SumPtInC10(0x0), fh2SumPtInC20(0x0), fh2SumPtInC30(0x0), fh2SumPtInC60(0x0), +fh2SumPtInC10lead(0x0), +fh2SumPtInC20lead(0x0), +fh2SumPtInC30lead(0x0), +fh2SumPtInC60lead(0x0), +fh2SumPtInC10sublead(0x0), +fh2SumPtInC20sublead(0x0), +fh2SumPtInC30sublead(0x0), +fh2SumPtInC60sublead(0x0), fh2SumPtOutC10(0x0), -fh2SumPtOutC10b(0x0), fh2SumPtOutC20(0x0), fh2SumPtOutC30(0x0), fh2SumPtOutC60(0x0), +fh2SumPtOutC10lead(0x0), +fh2SumPtOutC20lead(0x0), +fh2SumPtOutC30lead(0x0), +fh2SumPtOutC60lead(0x0), +fh2SumPtOutC10sublead(0x0), +fh2SumPtOutC20sublead(0x0), +fh2SumPtOutC30sublead(0x0), +fh2SumPtOutC60sublead(0x0), fh2SumPtInC10bkg(0x0), fh2SumPtInC20bkg(0x0), fh2SumPtInC30bkg(0x0), fh2SumPtInC60bkg(0x0), +fh2SumPtInC10bkglead(0x0), +fh2SumPtInC20bkglead(0x0), +fh2SumPtInC30bkglead(0x0), +fh2SumPtInC60bkglead(0x0), +fh2SumPtInC10bkgsublead(0x0), +fh2SumPtInC20bkgsublead(0x0), +fh2SumPtInC30bkgsublead(0x0), +fh2SumPtInC60bkgsublead(0x0), fh2SumPtOutC10bkg(0x0), fh2SumPtOutC20bkg(0x0), fh2SumPtOutC30bkg(0x0), fh2SumPtOutC60bkg(0x0), +fh2SumPtOutC10bkglead(0x0), +fh2SumPtOutC20bkglead(0x0), +fh2SumPtOutC30bkglead(0x0), +fh2SumPtOutC60bkglead(0x0), +fh2SumPtOutC10bkgsublead(0x0), +fh2SumPtOutC20bkgsublead(0x0), +fh2SumPtOutC30bkgsublead(0x0), +fh2SumPtOutC60bkgsublead(0x0), fh2DeltaRC10pt1(0x0), fh2DeltaRC20pt1(0x0), fh2DeltaRC30pt1(0x0), @@ -307,6 +481,102 @@ fh2DeltaPhiC10pt4(0x0), fh2DeltaPhiC20pt4(0x0), fh2DeltaPhiC30pt4(0x0), fh2DeltaPhiC60pt4(0x0), +fh2DeltaRC10pt1lead(0x0), +fh2DeltaRC20pt1lead(0x0), +fh2DeltaRC30pt1lead(0x0), +fh2DeltaRC60pt1lead(0x0), +fh2DeltaRC10pt2lead(0x0), +fh2DeltaRC20pt2lead(0x0), +fh2DeltaRC30pt2lead(0x0), +fh2DeltaRC60pt2lead(0x0), +fh2DeltaRC10pt3lead(0x0), +fh2DeltaRC20pt3lead(0x0), +fh2DeltaRC30pt3lead(0x0), +fh2DeltaRC60pt3lead(0x0), +fh2DeltaRC10pt4lead(0x0), +fh2DeltaRC20pt4lead(0x0), +fh2DeltaRC30pt4lead(0x0), +fh2DeltaRC60pt4lead(0x0), +fh2DeltaEtaC10pt1lead(0x0), +fh2DeltaEtaC20pt1lead(0x0), +fh2DeltaEtaC30pt1lead(0x0), +fh2DeltaEtaC60pt1lead(0x0), +fh2DeltaEtaC10pt2lead(0x0), +fh2DeltaEtaC20pt2lead(0x0), +fh2DeltaEtaC30pt2lead(0x0), +fh2DeltaEtaC60pt2lead(0x0), +fh2DeltaEtaC10pt3lead(0x0), +fh2DeltaEtaC20pt3lead(0x0), +fh2DeltaEtaC30pt3lead(0x0), +fh2DeltaEtaC60pt3lead(0x0), +fh2DeltaEtaC10pt4lead(0x0), +fh2DeltaEtaC20pt4lead(0x0), +fh2DeltaEtaC30pt4lead(0x0), +fh2DeltaEtaC60pt4lead(0x0), +fh2DeltaPhiC10pt1lead(0x0), +fh2DeltaPhiC20pt1lead(0x0), +fh2DeltaPhiC30pt1lead(0x0), +fh2DeltaPhiC60pt1lead(0x0), +fh2DeltaPhiC10pt2lead(0x0), +fh2DeltaPhiC20pt2lead(0x0), +fh2DeltaPhiC30pt2lead(0x0), +fh2DeltaPhiC60pt2lead(0x0), +fh2DeltaPhiC10pt3lead(0x0), +fh2DeltaPhiC20pt3lead(0x0), +fh2DeltaPhiC30pt3lead(0x0), +fh2DeltaPhiC60pt3lead(0x0), +fh2DeltaPhiC10pt4lead(0x0), +fh2DeltaPhiC20pt4lead(0x0), +fh2DeltaPhiC30pt4lead(0x0), +fh2DeltaPhiC60pt4lead(0x0), +fh2DeltaRC10pt1sublead(0x0), +fh2DeltaRC20pt1sublead(0x0), +fh2DeltaRC30pt1sublead(0x0), +fh2DeltaRC60pt1sublead(0x0), +fh2DeltaRC10pt2sublead(0x0), +fh2DeltaRC20pt2sublead(0x0), +fh2DeltaRC30pt2sublead(0x0), +fh2DeltaRC60pt2sublead(0x0), +fh2DeltaRC10pt3sublead(0x0), +fh2DeltaRC20pt3sublead(0x0), +fh2DeltaRC30pt3sublead(0x0), +fh2DeltaRC60pt3sublead(0x0), +fh2DeltaRC10pt4sublead(0x0), +fh2DeltaRC20pt4sublead(0x0), +fh2DeltaRC30pt4sublead(0x0), +fh2DeltaRC60pt4sublead(0x0), +fh2DeltaEtaC10pt1sublead(0x0), +fh2DeltaEtaC20pt1sublead(0x0), +fh2DeltaEtaC30pt1sublead(0x0), +fh2DeltaEtaC60pt1sublead(0x0), +fh2DeltaEtaC10pt2sublead(0x0), +fh2DeltaEtaC20pt2sublead(0x0), +fh2DeltaEtaC30pt2sublead(0x0), +fh2DeltaEtaC60pt2sublead(0x0), +fh2DeltaEtaC10pt3sublead(0x0), +fh2DeltaEtaC20pt3sublead(0x0), +fh2DeltaEtaC30pt3sublead(0x0), +fh2DeltaEtaC60pt3sublead(0x0), +fh2DeltaEtaC10pt4sublead(0x0), +fh2DeltaEtaC20pt4sublead(0x0), +fh2DeltaEtaC30pt4sublead(0x0), +fh2DeltaEtaC60pt4sublead(0x0), +fh2DeltaPhiC10pt1sublead(0x0), +fh2DeltaPhiC20pt1sublead(0x0), +fh2DeltaPhiC30pt1sublead(0x0), +fh2DeltaPhiC60pt1sublead(0x0), +fh2DeltaPhiC10pt2sublead(0x0), +fh2DeltaPhiC20pt2sublead(0x0), +fh2DeltaPhiC30pt2sublead(0x0), +fh2DeltaPhiC60pt2sublead(0x0), +fh2DeltaPhiC10pt3sublead(0x0), +fh2DeltaPhiC20pt3sublead(0x0), +fh2DeltaPhiC30pt3sublead(0x0), +fh2DeltaPhiC60pt3sublead(0x0), +fh2DeltaPhiC10pt4sublead(0x0), +fh2DeltaPhiC20pt4sublead(0x0), +fh2DeltaPhiC30pt4sublead(0x0), +fh2DeltaPhiC60pt4sublead(0x0), fh2AngStructpt1C10(0x0), fh2AngStructpt2C10(0x0), fh2AngStructpt3C10(0x0), @@ -421,25 +691,68 @@ void AliAnalysisTaskJetCore::UserCreateOutputObjects() fh2JetCoreMethod1C60 = new TH2F("JetCoreMethod1C60","",150, 0., 150.,100, 0., 1.5); fh2JetCoreMethod2C60 = new TH2F("JetCoreMethod2C60","",150, 0., 150.,100, 0., 1.5); fh2JetCoreMethod3C60 = new TH2F("JetCoreMethod3C60","",150, 0., 150.,100, 0., 1.5); + fh2JetCoreMethod3C10lead = new TH2F("JetCoreMethod3C10lead","",150, 0., 150.,100, 0., 1.5); + fh2JetCoreMethod3C20lead = new TH2F("JetCoreMethod3C20lead","",150, 0., 150.,100, 0., 1.5); + fh2JetCoreMethod3C30lead = new TH2F("JetCoreMethod3C30lead","",150, 0., 150.,100, 0., 1.5); + fh2JetCoreMethod3C60lead = new TH2F("JetCoreMethod3C60lead","",150, 0., 150.,100, 0., 1.5); + fh2JetCoreMethod3C10sublead = new TH2F("JetCoreMethod3C10sublead","",150, 0., 150.,100, 0., 1.5); + fh2JetCoreMethod3C20sublead = new TH2F("JetCoreMethod3C20sublead","",150, 0., 150.,100, 0., 1.5); + fh2JetCoreMethod3C30sublead = new TH2F("JetCoreMethod3C30sublead","",150, 0., 150.,100, 0., 1.5); + fh2JetCoreMethod3C60sublead = new TH2F("JetCoreMethod3C60sublead","",150, 0., 150.,100, 0., 1.5); - fh2SumPtInC10 = new TH2F("PtSumInC10","",150, 0., 150.,5000, 0., 50.); - fh2SumPtInC20 = new TH2F("PtSumInC20","",150, 0., 150.,5000, 0., 50.); - fh2SumPtInC30 = new TH2F("PtSumInC30","",150, 0., 150.,5000, 0., 50.); - fh2SumPtInC60 = new TH2F("PtSumInC60","",150, 0., 150.,5000, 0., 50.); - fh2SumPtOutC10 = new TH2F("PtSumOutC10","",150, 0., 150.,5000, 0., 50.); - fh2SumPtOutC10b= new TH2F("PtSumOutC10b","",150, 0., 150.,500, 0., 50.); - fh2SumPtOutC20 = new TH2F("PtSumOutC20","",150, 0., 150.,5000, 0., 50.); - fh2SumPtOutC30 = new TH2F("PtSumOutC30","",150, 0., 150.,5000, 0., 50.); - fh2SumPtOutC60 = new TH2F("PtSumOutC60","",150, 0., 150.,5000, 0., 50.); - - fh2SumPtInC10bkg = new TH2F("PtSumInC10bkg","",150, 0., 150.,5000, 0., 50.); - fh2SumPtInC20bkg = new TH2F("PtSumInC20bkg","",150, 0., 150.,5000, 0., 50.); - fh2SumPtInC30bkg = new TH2F("PtSumInC30bkg","",150, 0., 150.,5000, 0., 50.); - fh2SumPtInC60bkg = new TH2F("PtSumInC60bkg","",150, 0., 150.,5000, 0., 50.); - fh2SumPtOutC10bkg = new TH2F("PtSumOutC10bkg","",150, 0., 150.,5000, 0., 50.); - fh2SumPtOutC20bkg = new TH2F("PtSumOutC20bkg","",150, 0., 150.,5000, 0., 50.); - fh2SumPtOutC30bkg = new TH2F("PtSumOutC30bkg","",150, 0., 150.,5000, 0., 50.); - fh2SumPtOutC60bkg = new TH2F("PtSumOutC60bkg","",150, 0., 150.,5000, 0., 50.); + fh2SumPtInC10 = new TH2F("PtSumInC10","",150, 0., 150.,5000, 0., 25.); + fh2SumPtInC20 = new TH2F("PtSumInC20","",150, 0., 150.,5000, 0., 25.); + fh2SumPtInC30 = new TH2F("PtSumInC30","",150, 0., 150.,5000, 0., 25.); + fh2SumPtInC60 = new TH2F("PtSumInC60","",150, 0., 150.,5000, 0., 25.); + fh2SumPtInC10lead = new TH2F("PtSumInC10lead","",150, 0., 150.,5000, 0., 25.); + fh2SumPtInC20lead = new TH2F("PtSumInC20lead","",150, 0., 150.,5000, 0., 25.); + fh2SumPtInC30lead = new TH2F("PtSumInC30lead","",150, 0., 150.,5000, 0., 25.); + fh2SumPtInC60lead = new TH2F("PtSumInC60lead","",150, 0., 150.,5000, 0., 25.); + fh2SumPtInC10sublead = new TH2F("PtSumInC10sublead","",150, 0., 150.,5000, 0., 25.); + fh2SumPtInC20sublead = new TH2F("PtSumInC20sublead","",150, 0., 150.,5000, 0., 25.); + fh2SumPtInC30sublead = new TH2F("PtSumInC30sublead","",150, 0., 150.,5000, 0., 25.); + fh2SumPtInC60sublead = new TH2F("PtSumInC60sublead","",150, 0., 150.,5000, 0., 25.); + + fh2SumPtOutC10 = new TH2F("PtSumOutC10","",150, 0., 150.,5000, 0., 25.); + fh2SumPtOutC20 = new TH2F("PtSumOutC20","",150, 0., 150.,5000, 0., 25.); + fh2SumPtOutC30 = new TH2F("PtSumOutC30","",150, 0., 150.,5000, 0., 25.); + fh2SumPtOutC60 = new TH2F("PtSumOutC60","",150, 0., 150.,5000, 0., 25.); + fh2SumPtOutC10lead = new TH2F("PtSumOutC10lead","",150, 0., 150.,5000, 0., 25.); + fh2SumPtOutC20lead = new TH2F("PtSumOutC20lead","",150, 0., 150.,5000, 0., 25.); + fh2SumPtOutC30lead = new TH2F("PtSumOutC30lead","",150, 0., 150.,5000, 0., 25.); + fh2SumPtOutC60lead = new TH2F("PtSumOutC60lead","",150, 0., 150.,5000, 0., 25.); + fh2SumPtOutC10sublead = new TH2F("PtSumOutC10sublead","",150, 0., 150.,5000, 0., 25.); + fh2SumPtOutC20sublead = new TH2F("PtSumOutC20sublead","",150, 0., 150.,5000, 0., 25.); + fh2SumPtOutC30sublead = new TH2F("PtSumOutC30sublead","",150, 0., 150.,5000, 0., 25.); + fh2SumPtOutC60sublead = new TH2F("PtSumOutC60sublead","",150, 0., 150.,5000, 0., 25.); + + + fh2SumPtInC10bkg = new TH2F("PtSumInC10bkg","",150, 0., 150.,5000, 0., 25.); + fh2SumPtInC20bkg = new TH2F("PtSumInC20bkg","",150, 0., 150.,5000, 0., 25.); + fh2SumPtInC30bkg = new TH2F("PtSumInC30bkg","",150, 0., 150.,5000, 0., 25.); + fh2SumPtInC60bkg = new TH2F("PtSumInC60bkg","",150, 0., 150.,5000, 0., 25.); + fh2SumPtInC10bkglead = new TH2F("PtSumInC10bkglead","",150, 0., 150.,5000, 0., 25.); + fh2SumPtInC20bkglead = new TH2F("PtSumInC20bkglead","",150, 0., 150.,5000, 0., 25.); + fh2SumPtInC30bkglead = new TH2F("PtSumInC30bkglead","",150, 0., 150.,5000, 0., 25.); + fh2SumPtInC60bkglead = new TH2F("PtSumInC60bkglead","",150, 0., 150.,5000, 0., 25.); + fh2SumPtInC10bkgsublead = new TH2F("PtSumInC10bkgsublead","",150, 0., 150.,5000, 0., 25.); + fh2SumPtInC20bkgsublead = new TH2F("PtSumInC20bkgsublead","",150, 0., 150.,5000, 0., 25.); + fh2SumPtInC30bkgsublead = new TH2F("PtSumInC30bkgsublead","",150, 0., 150.,5000, 0., 25.); + fh2SumPtInC60bkgsublead = new TH2F("PtSumInC60bkgsublead","",150, 0., 150.,5000, 0., 25.); + + + fh2SumPtOutC10bkg = new TH2F("PtSumOutC10bkg","",150, 0., 150.,5000, 0., 25.); + fh2SumPtOutC20bkg = new TH2F("PtSumOutC20bkg","",150, 0., 150.,5000, 0., 25.); + fh2SumPtOutC30bkg = new TH2F("PtSumOutC30bkg","",150, 0., 150.,5000, 0., 25.); + fh2SumPtOutC60bkg = new TH2F("PtSumOutC60bkg","",150, 0., 150.,5000, 0., 25.); + fh2SumPtOutC10bkglead = new TH2F("PtSumOutC10bkglead","",150, 0., 150.,5000, 0., 25.); + fh2SumPtOutC20bkglead = new TH2F("PtSumOutC20bkglead","",150, 0., 150.,5000, 0., 25.); + fh2SumPtOutC30bkglead = new TH2F("PtSumOutC30bkglead","",150, 0., 150.,5000, 0., 25.); + fh2SumPtOutC60bkglead = new TH2F("PtSumOutC60bkglead","",150, 0., 150.,5000, 0., 25.); + fh2SumPtOutC10bkgsublead = new TH2F("PtSumOutC10bkgsublead","",150, 0., 150.,5000, 0., 25.); + fh2SumPtOutC20bkgsublead = new TH2F("PtSumOutC20bkgsublead","",150, 0., 150.,5000, 0., 25.); + fh2SumPtOutC30bkgsublead = new TH2F("PtSumOutC30bkgsublead","",150, 0., 150.,5000, 0., 25.); + fh2SumPtOutC60bkgsublead = new TH2F("PtSumOutC60bkgsublead","",150, 0., 150.,5000, 0., 25.); @@ -459,7 +772,6 @@ void AliAnalysisTaskJetCore::UserCreateOutputObjects() fh2DeltaRC20pt4 = new TH2F("DeltaRC20pt4","",150, 0., 150.,100,0.,1.5); fh2DeltaRC30pt4 = new TH2F("DeltaRC30pt4","",150, 0., 150.,100,0.,1.5); fh2DeltaRC60pt4 = new TH2F("DeltaRC60pt4","",150, 0., 150.,100,0.,1.5); - fh2DeltaEtaC10pt1 = new TH2F("DeltaEtaC10pt1","",150, 0., 150.,100,-1.5,1.5); fh2DeltaEtaC20pt1 = new TH2F("DeltaEtaC20pt1","",150, 0., 150.,100,-1.5,1.5); fh2DeltaEtaC30pt1 = new TH2F("DeltaEtaC30pt1","",150, 0., 150.,100,-1.5,1.5); @@ -493,26 +805,131 @@ void AliAnalysisTaskJetCore::UserCreateOutputObjects() fh2DeltaPhiC30pt4 = new TH2F("DeltaPhiC30pt4","",150, 0., 150.,100,-6.5,6.5); fh2DeltaPhiC60pt4 = new TH2F("DeltaPhiC60pt4","",150, 0., 150.,100,-6.5,6.5); - fh2AngStructpt1C10 = new TH2F("Ang struct pt1 C10","",15,0.,1.5,150,0.,5.); - fh2AngStructpt2C10 = new TH2F("Ang struct pt2 C10","",15,0.,1.5,150,0.,5.); - fh2AngStructpt3C10 = new TH2F("Ang struct pt3 C10","",15,0.,1.5,150,0.,5.); - fh2AngStructpt4C10 = new TH2F("Ang struct pt4 C10","",15,0.,1.5,150,0.,5.); - fh2AngStructpt1C20 = new TH2F("Ang struct pt1 C20","",15,0.,1.5,150,0.,5.); - fh2AngStructpt2C20 = new TH2F("Ang struct pt2 C20","",15,0.,1.5,150,0.,5.); - fh2AngStructpt3C20 = new TH2F("Ang struct pt3 C20","",15,0.,1.5,150,0.,5.); - fh2AngStructpt4C20 = new TH2F("Ang struct pt4 C20","",15,0.,1.5,150,0.,5.); - fh2AngStructpt1C30 = new TH2F("Ang struct pt1 C30","",15,0.,1.5,150,0.,5.); - fh2AngStructpt2C30 = new TH2F("Ang struct pt2 C30","",15,0.,1.5,150,0.,5.); - fh2AngStructpt3C30 = new TH2F("Ang struct pt3 C30","",15,0.,1.5,150,0.,5.); - fh2AngStructpt4C30 = new TH2F("Ang struct pt4 C30","",15,0.,1.5,150,0.,5.); - fh2AngStructpt1C60 = new TH2F("Ang struct pt1 C60","",15,0.,1.5,150,0.,5.); - fh2AngStructpt2C60 = new TH2F("Ang struct pt2 C60","",15,0.,1.5,150,0.,5.); - fh2AngStructpt3C60 = new TH2F("Ang struct pt3 C60","",15,0.,1.5,150,0.,5.); - fh2AngStructpt4C60 = new TH2F("Ang struct pt4 C60","",15,0.,1.5,150,0.,5.); + fh2DeltaRC10pt1lead = new TH2F("DeltaRC10pt1lead","",150, 0., 150.,100,0.,1.5); + fh2DeltaRC20pt1lead = new TH2F("DeltaRC20pt1lead","",150, 0., 150.,100,0.,1.5); + fh2DeltaRC30pt1lead = new TH2F("DeltaRC30pt1lead","",150, 0., 150.,100,0.,1.5); + fh2DeltaRC60pt1lead = new TH2F("DeltaRC60pt1lead","",150, 0., 150.,100,0.,1.5); + fh2DeltaRC10pt2lead = new TH2F("DeltaRC10pt2lead","",150, 0., 150.,100,0.,1.5); + fh2DeltaRC20pt2lead = new TH2F("DeltaRC20pt2lead","",150, 0., 150.,100,0.,1.5); + fh2DeltaRC30pt2lead = new TH2F("DeltaRC30pt2lead","",150, 0., 150.,100,0.,1.5); + fh2DeltaRC60pt2lead = new TH2F("DeltaRC60pt2lead","",150, 0., 150.,100,0.,1.5); + fh2DeltaRC10pt3lead = new TH2F("DeltaRC10pt3lead","",150, 0., 150.,100,0.,1.5); + fh2DeltaRC20pt3lead = new TH2F("DeltaRC20pt3lead","",150, 0., 150.,100,0.,1.5); + fh2DeltaRC30pt3lead = new TH2F("DeltaRC30pt3lead","",150, 0., 150.,100,0.,1.5); + fh2DeltaRC60pt3lead = new TH2F("DeltaRC60pt3lead","",150, 0., 150.,100,0.,1.5); + fh2DeltaRC10pt4lead = new TH2F("DeltaRC10pt4lead","",150, 0., 150.,100,0.,1.5); + fh2DeltaRC20pt4lead = new TH2F("DeltaRC20pt4lead","",150, 0., 150.,100,0.,1.5); + fh2DeltaRC30pt4lead = new TH2F("DeltaRC30pt4lead","",150, 0., 150.,100,0.,1.5); + fh2DeltaRC60pt4lead = new TH2F("DeltaRC60pt4lead","",150, 0., 150.,100,0.,1.5); + fh2DeltaEtaC10pt1lead = new TH2F("DeltaEtaC10pt1lead","",150, 0., 150.,100,-1.5,1.5); + fh2DeltaEtaC20pt1lead = new TH2F("DeltaEtaC20pt1lead","",150, 0., 150.,100,-1.5,1.5); + fh2DeltaEtaC30pt1lead = new TH2F("DeltaEtaC30pt1lead","",150, 0., 150.,100,-1.5,1.5); + fh2DeltaEtaC60pt1lead = new TH2F("DeltaEtaC60pt1lead","",150, 0., 150.,100,-1.5,1.5); + fh2DeltaEtaC10pt2lead = new TH2F("DeltaEtaC10pt2lead","",150, 0., 150.,100,-1.5,1.5); + fh2DeltaEtaC20pt2lead = new TH2F("DeltaEtaC20pt2lead","",150, 0., 150.,100,-1.5,1.5); + fh2DeltaEtaC30pt2lead = new TH2F("DeltaEtaC30pt2lead","",150, 0., 150.,100,-1.5,1.5); + fh2DeltaEtaC60pt2lead = new TH2F("DeltaEtaC60pt2lead","",150, 0., 150.,100,-1.5,1.5); + fh2DeltaEtaC10pt3lead = new TH2F("DeltaEtaC10pt3lead","",150, 0., 150.,100,-1.5,1.5); + fh2DeltaEtaC20pt3lead = new TH2F("DeltaEtaC20pt3lead","",150, 0., 150.,100,-1.5,1.5); + fh2DeltaEtaC30pt3lead = new TH2F("DeltaEtaC30pt3lead","",150, 0., 150.,100,-1.5,1.5); + fh2DeltaEtaC60pt3lead = new TH2F("DeltaEtaC60pt3lead","",150, 0., 150.,100,-1.5,1.5); + fh2DeltaEtaC10pt4lead = new TH2F("DeltaEtaC10pt4lead","",150, 0., 150.,100,-1.5,1.5); + fh2DeltaEtaC20pt4lead = new TH2F("DeltaEtaC20pt4lead","",150, 0., 150.,100,-1.5,1.5); + fh2DeltaEtaC30pt4lead = new TH2F("DeltaEtaC30pt4lead","",150, 0., 150.,100,-1.5,1.5); + fh2DeltaEtaC60pt4lead = new TH2F("DeltaEtaC60pt4lead","",150, 0., 150.,100,-1.5,1.5); + fh2DeltaPhiC10pt1lead = new TH2F("DeltaPhiC10pt1lead","",150, 0., 150.,100,-6.5,6.5); + fh2DeltaPhiC20pt1lead = new TH2F("DeltaPhiC20pt1lead","",150, 0., 150.,100,-6.5,6.5); + fh2DeltaPhiC30pt1lead = new TH2F("DeltaPhiC30pt1lead","",150, 0., 150.,100,-6.5,6.5); + fh2DeltaPhiC60pt1lead = new TH2F("DeltaPhiC60pt1lead","",150, 0., 150.,100,-6.5,6.5); + fh2DeltaPhiC10pt2lead = new TH2F("DeltaPhiC10pt2lead","",150, 0., 150.,100,-6.5,6.5); + fh2DeltaPhiC20pt2lead = new TH2F("DeltaPhiC20pt2lead","",150, 0., 150.,100,-6.5,6.5); + fh2DeltaPhiC30pt2lead = new TH2F("DeltaPhiC30pt2lead","",150, 0., 150.,100,-6.5,6.5); + fh2DeltaPhiC60pt2lead = new TH2F("DeltaPhiC60pt2lead","",150, 0., 150.,100,-6.5,6.5); + fh2DeltaPhiC10pt3lead = new TH2F("DeltaPhiC10pt3lead","",150, 0., 150.,100,-6.5,6.5); + fh2DeltaPhiC20pt3lead = new TH2F("DeltaPhiC20pt3lead","",150, 0., 150.,100,-6.5,6.5); + fh2DeltaPhiC30pt3lead = new TH2F("DeltaPhiC30pt3lead","",150, 0., 150.,100,-6.5,6.5); + fh2DeltaPhiC60pt3lead = new TH2F("DeltaPhiC60pt3lead","",150, 0., 150.,100,-6.5,6.5); + fh2DeltaPhiC10pt4lead = new TH2F("DeltaPhiC10pt4lead","",150, 0., 150.,100,-6.5,6.5); + fh2DeltaPhiC20pt4lead = new TH2F("DeltaPhiC20pt4lead","",150, 0., 150.,100,-6.5,6.5); + fh2DeltaPhiC30pt4lead = new TH2F("DeltaPhiC30pt4lead","",150, 0., 150.,100,-6.5,6.5); + fh2DeltaPhiC60pt4lead = new TH2F("DeltaPhiC60pt4lead","",150, 0., 150.,100,-6.5,6.5); + + + fh2DeltaRC10pt1sublead = new TH2F("DeltaRC10pt1sublead","",150, 0., 150.,100,0.,1.5); + fh2DeltaRC20pt1sublead = new TH2F("DeltaRC20pt1sublead","",150, 0., 150.,100,0.,1.5); + fh2DeltaRC30pt1sublead = new TH2F("DeltaRC30pt1sublead","",150, 0., 150.,100,0.,1.5); + fh2DeltaRC60pt1sublead = new TH2F("DeltaRC60pt1sublead","",150, 0., 150.,100,0.,1.5); + fh2DeltaRC10pt2sublead = new TH2F("DeltaRC10pt2sublead","",150, 0., 150.,100,0.,1.5); + fh2DeltaRC20pt2sublead = new TH2F("DeltaRC20pt2sublead","",150, 0., 150.,100,0.,1.5); + fh2DeltaRC30pt2sublead = new TH2F("DeltaRC30pt2sublead","",150, 0., 150.,100,0.,1.5); + fh2DeltaRC60pt2sublead = new TH2F("DeltaRC60pt2sublead","",150, 0., 150.,100,0.,1.5); + fh2DeltaRC10pt3sublead = new TH2F("DeltaRC10pt3sublead","",150, 0., 150.,100,0.,1.5); + fh2DeltaRC20pt3sublead = new TH2F("DeltaRC20pt3sublead","",150, 0., 150.,100,0.,1.5); + fh2DeltaRC30pt3sublead = new TH2F("DeltaRC30pt3sublead","",150, 0., 150.,100,0.,1.5); + fh2DeltaRC60pt3sublead = new TH2F("DeltaRC60pt3sublead","",150, 0., 150.,100,0.,1.5); + fh2DeltaRC10pt4sublead = new TH2F("DeltaRC10pt4sublead","",150, 0., 150.,100,0.,1.5); + fh2DeltaRC20pt4sublead = new TH2F("DeltaRC20pt4sublead","",150, 0., 150.,100,0.,1.5); + fh2DeltaRC30pt4sublead = new TH2F("DeltaRC30pt4sublead","",150, 0., 150.,100,0.,1.5); + fh2DeltaRC60pt4sublead = new TH2F("DeltaRC60pt4sublead","",150, 0., 150.,100,0.,1.5); + fh2DeltaEtaC10pt1sublead = new TH2F("DeltaEtaC10pt1sublead","",150, 0., 150.,100,-1.5,1.5); + fh2DeltaEtaC20pt1sublead = new TH2F("DeltaEtaC20pt1sublead","",150, 0., 150.,100,-1.5,1.5); + fh2DeltaEtaC30pt1sublead = new TH2F("DeltaEtaC30pt1sublead","",150, 0., 150.,100,-1.5,1.5); + fh2DeltaEtaC60pt1sublead = new TH2F("DeltaEtaC60pt1sublead","",150, 0., 150.,100,-1.5,1.5); + fh2DeltaEtaC10pt2sublead = new TH2F("DeltaEtaC10pt2sublead","",150, 0., 150.,100,-1.5,1.5); + fh2DeltaEtaC20pt2sublead = new TH2F("DeltaEtaC20pt2sublead","",150, 0., 150.,100,-1.5,1.5); + fh2DeltaEtaC30pt2sublead = new TH2F("DeltaEtaC30pt2sublead","",150, 0., 150.,100,-1.5,1.5); + fh2DeltaEtaC60pt2sublead = new TH2F("DeltaEtaC60pt2sublead","",150, 0., 150.,100,-1.5,1.5); + fh2DeltaEtaC10pt3sublead = new TH2F("DeltaEtaC10pt3sublead","",150, 0., 150.,100,-1.5,1.5); + fh2DeltaEtaC20pt3sublead = new TH2F("DeltaEtaC20pt3sublead","",150, 0., 150.,100,-1.5,1.5); + fh2DeltaEtaC30pt3sublead = new TH2F("DeltaEtaC30pt3sublead","",150, 0., 150.,100,-1.5,1.5); + fh2DeltaEtaC60pt3sublead = new TH2F("DeltaEtaC60pt3sublead","",150, 0., 150.,100,-1.5,1.5); + fh2DeltaEtaC10pt4sublead = new TH2F("DeltaEtaC10pt4sublead","",150, 0., 150.,100,-1.5,1.5); + fh2DeltaEtaC20pt4sublead = new TH2F("DeltaEtaC20pt4sublead","",150, 0., 150.,100,-1.5,1.5); + fh2DeltaEtaC30pt4sublead = new TH2F("DeltaEtaC30pt4sublead","",150, 0., 150.,100,-1.5,1.5); + fh2DeltaEtaC60pt4sublead = new TH2F("DeltaEtaC60pt4sublead","",150, 0., 150.,100,-1.5,1.5); + fh2DeltaPhiC10pt1sublead = new TH2F("DeltaPhiC10pt1sublead","",150, 0., 150.,100,-6.5,6.5); + fh2DeltaPhiC20pt1sublead = new TH2F("DeltaPhiC20pt1sublead","",150, 0., 150.,100,-6.5,6.5); + fh2DeltaPhiC30pt1sublead = new TH2F("DeltaPhiC30pt1sublead","",150, 0., 150.,100,-6.5,6.5); + fh2DeltaPhiC60pt1sublead = new TH2F("DeltaPhiC60pt1sublead","",150, 0., 150.,100,-6.5,6.5); + fh2DeltaPhiC10pt2sublead = new TH2F("DeltaPhiC10pt2sublead","",150, 0., 150.,100,-6.5,6.5); + fh2DeltaPhiC20pt2sublead = new TH2F("DeltaPhiC20pt2sublead","",150, 0., 150.,100,-6.5,6.5); + fh2DeltaPhiC30pt2sublead = new TH2F("DeltaPhiC30pt2sublead","",150, 0., 150.,100,-6.5,6.5); + fh2DeltaPhiC60pt2sublead = new TH2F("DeltaPhiC60pt2sublead","",150, 0., 150.,100,-6.5,6.5); + fh2DeltaPhiC10pt3sublead = new TH2F("DeltaPhiC10pt3sublead","",150, 0., 150.,100,-6.5,6.5); + fh2DeltaPhiC20pt3sublead = new TH2F("DeltaPhiC20pt3sublead","",150, 0., 150.,100,-6.5,6.5); + fh2DeltaPhiC30pt3sublead = new TH2F("DeltaPhiC30pt3sublead","",150, 0., 150.,100,-6.5,6.5); + fh2DeltaPhiC60pt3sublead = new TH2F("DeltaPhiC60pt3sublead","",150, 0., 150.,100,-6.5,6.5); + fh2DeltaPhiC10pt4sublead = new TH2F("DeltaPhiC10pt4sublead","",150, 0., 150.,100,-6.5,6.5); + fh2DeltaPhiC20pt4sublead = new TH2F("DeltaPhiC20pt4sublead","",150, 0., 150.,100,-6.5,6.5); + fh2DeltaPhiC30pt4sublead = new TH2F("DeltaPhiC30pt4sublead","",150, 0., 150.,100,-6.5,6.5); + fh2DeltaPhiC60pt4sublead = new TH2F("DeltaPhiC60pt4sublead","",150, 0., 150.,100,-6.5,6.5); + + + + + + if(fAngStructCloseTracks>0){ + fh2AngStructpt1C10 = new TH2F("Ang struct pt1 C10","",15,0.,1.5,150,0.,10.); + fh2AngStructpt2C10 = new TH2F("Ang struct pt2 C10","",15,0.,1.5,150,0.,10.); + fh2AngStructpt3C10 = new TH2F("Ang struct pt3 C10","",15,0.,1.5,150,0.,10.); + fh2AngStructpt4C10 = new TH2F("Ang struct pt4 C10","",15,0.,1.5,150,0.,10.); + fh2AngStructpt1C20 = new TH2F("Ang struct pt1 C20","",15,0.,1.5,150,0.,10.); + fh2AngStructpt2C20 = new TH2F("Ang struct pt2 C20","",15,0.,1.5,150,0.,10.); + fh2AngStructpt3C20 = new TH2F("Ang struct pt3 C20","",15,0.,1.5,150,0.,10.); + fh2AngStructpt4C20 = new TH2F("Ang struct pt4 C20","",15,0.,1.5,150,0.,10.); + fh2AngStructpt1C30 = new TH2F("Ang struct pt1 C30","",15,0.,1.5,150,0.,10.); + fh2AngStructpt2C30 = new TH2F("Ang struct pt2 C30","",15,0.,1.5,150,0.,10.); + fh2AngStructpt3C30 = new TH2F("Ang struct pt3 C30","",15,0.,1.5,150,0.,10.); + fh2AngStructpt4C30 = new TH2F("Ang struct pt4 C30","",15,0.,1.5,150,0.,10.); + fh2AngStructpt1C60 = new TH2F("Ang struct pt1 C60","",15,0.,1.5,150,0.,10.); + fh2AngStructpt2C60 = new TH2F("Ang struct pt2 C60","",15,0.,1.5,150,0.,10.); + fh2AngStructpt3C60 = new TH2F("Ang struct pt3 C60","",15,0.,1.5,150,0.,10.); + fh2AngStructpt4C60 = new TH2F("Ang struct pt4 C60","",15,0.,1.5,150,0.,10.); } fOutputList->Add(fHistEvtSelection); fOutputList->Add(fHistJetSelection); fOutputList->Add(fh2JetSelection); + fOutputList->Add(fh2JetCoreMethod1C10); fOutputList->Add(fh2JetCoreMethod2C10); fOutputList->Add(fh2JetCoreMethod3C10); @@ -525,24 +942,68 @@ void AliAnalysisTaskJetCore::UserCreateOutputObjects() fOutputList->Add(fh2JetCoreMethod1C60); fOutputList->Add(fh2JetCoreMethod2C60); fOutputList->Add(fh2JetCoreMethod3C60); + + fOutputList->Add(fh2JetCoreMethod3C10lead); + fOutputList->Add(fh2JetCoreMethod3C20lead); + fOutputList->Add(fh2JetCoreMethod3C30lead); + fOutputList->Add(fh2JetCoreMethod3C60lead); + fOutputList->Add(fh2JetCoreMethod3C10sublead); + fOutputList->Add(fh2JetCoreMethod3C20sublead); + fOutputList->Add(fh2JetCoreMethod3C30sublead); + fOutputList->Add(fh2JetCoreMethod3C60sublead); + fOutputList->Add(fh2SumPtInC10); fOutputList->Add(fh2SumPtInC20); fOutputList->Add(fh2SumPtInC30); fOutputList->Add(fh2SumPtInC60); + fOutputList->Add(fh2SumPtInC10lead); + fOutputList->Add(fh2SumPtInC20lead); + fOutputList->Add(fh2SumPtInC30lead); + fOutputList->Add(fh2SumPtInC60lead); + fOutputList->Add(fh2SumPtInC10sublead); + fOutputList->Add(fh2SumPtInC20sublead); + fOutputList->Add(fh2SumPtInC30sublead); + fOutputList->Add(fh2SumPtInC60sublead); + fOutputList->Add(fh2SumPtOutC10); - fOutputList->Add(fh2SumPtOutC10b); fOutputList->Add(fh2SumPtOutC20); fOutputList->Add(fh2SumPtOutC30); fOutputList->Add(fh2SumPtOutC60); + fOutputList->Add(fh2SumPtOutC10lead); + fOutputList->Add(fh2SumPtOutC20lead); + fOutputList->Add(fh2SumPtOutC30lead); + fOutputList->Add(fh2SumPtOutC60lead); + fOutputList->Add(fh2SumPtOutC10sublead); + fOutputList->Add(fh2SumPtOutC20sublead); + fOutputList->Add(fh2SumPtOutC30sublead); + fOutputList->Add(fh2SumPtOutC60sublead); fOutputList->Add(fh2SumPtInC10bkg); fOutputList->Add(fh2SumPtInC20bkg); fOutputList->Add(fh2SumPtInC30bkg); fOutputList->Add(fh2SumPtInC60bkg); + fOutputList->Add(fh2SumPtInC10bkglead); + fOutputList->Add(fh2SumPtInC20bkglead); + fOutputList->Add(fh2SumPtInC30bkglead); + fOutputList->Add(fh2SumPtInC60bkglead); + fOutputList->Add(fh2SumPtInC10bkgsublead); + fOutputList->Add(fh2SumPtInC20bkgsublead); + fOutputList->Add(fh2SumPtInC30bkgsublead); + fOutputList->Add(fh2SumPtInC60bkgsublead); + fOutputList->Add(fh2SumPtOutC10bkg); fOutputList->Add(fh2SumPtOutC20bkg); fOutputList->Add(fh2SumPtOutC30bkg); fOutputList->Add(fh2SumPtOutC60bkg); + fOutputList->Add(fh2SumPtOutC10bkglead); + fOutputList->Add(fh2SumPtOutC20bkglead); + fOutputList->Add(fh2SumPtOutC30bkglead); + fOutputList->Add(fh2SumPtOutC60bkglead); + fOutputList->Add(fh2SumPtOutC10bkgsublead); + fOutputList->Add(fh2SumPtOutC20bkgsublead); + fOutputList->Add(fh2SumPtOutC30bkgsublead); + fOutputList->Add(fh2SumPtOutC60bkgsublead); + fOutputList->Add(fh2DeltaRC10pt1); fOutputList->Add(fh2DeltaRC20pt1); fOutputList->Add(fh2DeltaRC30pt1); @@ -559,7 +1020,6 @@ void AliAnalysisTaskJetCore::UserCreateOutputObjects() fOutputList->Add(fh2DeltaRC20pt4); fOutputList->Add(fh2DeltaRC30pt4); fOutputList->Add(fh2DeltaRC60pt4); - fOutputList->Add(fh2DeltaEtaC10pt1); fOutputList->Add(fh2DeltaEtaC20pt1); fOutputList->Add(fh2DeltaEtaC30pt1); @@ -592,7 +1052,112 @@ void AliAnalysisTaskJetCore::UserCreateOutputObjects() fOutputList->Add(fh2DeltaPhiC20pt4); fOutputList->Add(fh2DeltaPhiC30pt4); fOutputList->Add(fh2DeltaPhiC60pt4); + + fOutputList->Add(fh2DeltaRC10pt1lead); + fOutputList->Add(fh2DeltaRC20pt1lead); + fOutputList->Add(fh2DeltaRC30pt1lead); + fOutputList->Add(fh2DeltaRC60pt1lead); + fOutputList->Add(fh2DeltaRC10pt2lead); + fOutputList->Add(fh2DeltaRC20pt2lead); + fOutputList->Add(fh2DeltaRC30pt2lead); + fOutputList->Add(fh2DeltaRC60pt2lead); + fOutputList->Add(fh2DeltaRC10pt3lead); + fOutputList->Add(fh2DeltaRC20pt3lead); + fOutputList->Add(fh2DeltaRC30pt3lead); + fOutputList->Add(fh2DeltaRC60pt3lead); + fOutputList->Add(fh2DeltaRC10pt4lead); + fOutputList->Add(fh2DeltaRC20pt4lead); + fOutputList->Add(fh2DeltaRC30pt4lead); + fOutputList->Add(fh2DeltaRC60pt4lead); + fOutputList->Add(fh2DeltaEtaC10pt1lead); + fOutputList->Add(fh2DeltaEtaC20pt1lead); + fOutputList->Add(fh2DeltaEtaC30pt1lead); + fOutputList->Add(fh2DeltaEtaC60pt1lead); + fOutputList->Add(fh2DeltaEtaC10pt2lead); + fOutputList->Add(fh2DeltaEtaC20pt2lead); + fOutputList->Add(fh2DeltaEtaC30pt2lead); + fOutputList->Add(fh2DeltaEtaC60pt2lead); + fOutputList->Add(fh2DeltaEtaC10pt3lead); + fOutputList->Add(fh2DeltaEtaC20pt3lead); + fOutputList->Add(fh2DeltaEtaC30pt3lead); + fOutputList->Add(fh2DeltaEtaC60pt3lead); + fOutputList->Add(fh2DeltaEtaC10pt4lead); + fOutputList->Add(fh2DeltaEtaC20pt4lead); + fOutputList->Add(fh2DeltaEtaC30pt4lead); + fOutputList->Add(fh2DeltaEtaC60pt4lead); + fOutputList->Add(fh2DeltaPhiC10pt1lead); + fOutputList->Add(fh2DeltaPhiC20pt1lead); + fOutputList->Add(fh2DeltaPhiC30pt1lead); + fOutputList->Add(fh2DeltaPhiC60pt1lead); + fOutputList->Add(fh2DeltaPhiC10pt2lead); + fOutputList->Add(fh2DeltaPhiC20pt2lead); + fOutputList->Add(fh2DeltaPhiC30pt2lead); + fOutputList->Add(fh2DeltaPhiC60pt2lead); + fOutputList->Add(fh2DeltaPhiC10pt3lead); + fOutputList->Add(fh2DeltaPhiC20pt3lead); + fOutputList->Add(fh2DeltaPhiC30pt3lead); + fOutputList->Add(fh2DeltaPhiC60pt3lead); + fOutputList->Add(fh2DeltaPhiC10pt4lead); + fOutputList->Add(fh2DeltaPhiC20pt4lead); + fOutputList->Add(fh2DeltaPhiC30pt4lead); + fOutputList->Add(fh2DeltaPhiC60pt4lead); + + + + fOutputList->Add(fh2DeltaRC10pt1sublead); + fOutputList->Add(fh2DeltaRC20pt1sublead); + fOutputList->Add(fh2DeltaRC30pt1sublead); + fOutputList->Add(fh2DeltaRC60pt1sublead); + fOutputList->Add(fh2DeltaRC10pt2sublead); + fOutputList->Add(fh2DeltaRC20pt2sublead); + fOutputList->Add(fh2DeltaRC30pt2sublead); + fOutputList->Add(fh2DeltaRC60pt2sublead); + fOutputList->Add(fh2DeltaRC10pt3sublead); + fOutputList->Add(fh2DeltaRC20pt3sublead); + fOutputList->Add(fh2DeltaRC30pt3sublead); + fOutputList->Add(fh2DeltaRC60pt3sublead); + fOutputList->Add(fh2DeltaRC10pt4sublead); + fOutputList->Add(fh2DeltaRC20pt4sublead); + fOutputList->Add(fh2DeltaRC30pt4sublead); + fOutputList->Add(fh2DeltaRC60pt4sublead); + fOutputList->Add(fh2DeltaEtaC10pt1sublead); + fOutputList->Add(fh2DeltaEtaC20pt1sublead); + fOutputList->Add(fh2DeltaEtaC30pt1sublead); + fOutputList->Add(fh2DeltaEtaC60pt1sublead); + fOutputList->Add(fh2DeltaEtaC10pt2sublead); + fOutputList->Add(fh2DeltaEtaC20pt2sublead); + fOutputList->Add(fh2DeltaEtaC30pt2sublead); + fOutputList->Add(fh2DeltaEtaC60pt2sublead); + fOutputList->Add(fh2DeltaEtaC10pt3sublead); + fOutputList->Add(fh2DeltaEtaC20pt3sublead); + fOutputList->Add(fh2DeltaEtaC30pt3sublead); + fOutputList->Add(fh2DeltaEtaC60pt3sublead); + fOutputList->Add(fh2DeltaEtaC10pt4sublead); + fOutputList->Add(fh2DeltaEtaC20pt4sublead); + fOutputList->Add(fh2DeltaEtaC30pt4sublead); + fOutputList->Add(fh2DeltaEtaC60pt4sublead); + fOutputList->Add(fh2DeltaPhiC10pt1sublead); + fOutputList->Add(fh2DeltaPhiC20pt1sublead); + fOutputList->Add(fh2DeltaPhiC30pt1sublead); + fOutputList->Add(fh2DeltaPhiC60pt1sublead); + fOutputList->Add(fh2DeltaPhiC10pt2sublead); + fOutputList->Add(fh2DeltaPhiC20pt2sublead); + fOutputList->Add(fh2DeltaPhiC30pt2sublead); + fOutputList->Add(fh2DeltaPhiC60pt2sublead); + fOutputList->Add(fh2DeltaPhiC10pt3sublead); + fOutputList->Add(fh2DeltaPhiC20pt3sublead); + fOutputList->Add(fh2DeltaPhiC30pt3sublead); + fOutputList->Add(fh2DeltaPhiC60pt3sublead); + fOutputList->Add(fh2DeltaPhiC10pt4sublead); + fOutputList->Add(fh2DeltaPhiC20pt4sublead); + fOutputList->Add(fh2DeltaPhiC30pt4sublead); + fOutputList->Add(fh2DeltaPhiC60pt4sublead); + + + + + if(fAngStructCloseTracks>0){ fOutputList->Add(fh2AngStructpt1C10); fOutputList->Add(fh2AngStructpt2C10); fOutputList->Add(fh2AngStructpt3C10); @@ -608,7 +1173,7 @@ void AliAnalysisTaskJetCore::UserCreateOutputObjects() fOutputList->Add(fh2AngStructpt1C60); fOutputList->Add(fh2AngStructpt2C60); fOutputList->Add(fh2AngStructpt3C60); - fOutputList->Add(fh2AngStructpt4C60); + fOutputList->Add(fh2AngStructpt4C60);} @@ -738,12 +1303,18 @@ void AliAnalysisTaskJetCore::UserExec(Option_t *) TClonesArray *aodJets[2]; aodJets[0]=0; if(fAOD&&!aodJets[0]){ - aodJets[0] = dynamic_cast(fAOD->FindListObject(fJetBranchName[0].Data())); // + aodJets[0] = dynamic_cast(fAOD->FindListObject(fJetBranchName[0].Data())); aodJets[1] = dynamic_cast(fAOD->FindListObject(fJetBranchName[1].Data())); } if(fAODExtension && !aodJets[0]){ - aodJets[0] = dynamic_cast(fAODExtension->GetAOD()->FindListObject(fJetBranchName[0].Data())); // + aodJets[0] = dynamic_cast(fAODExtension->GetAOD()->FindListObject(fJetBranchName[0].Data())); aodJets[1] = dynamic_cast(fAODExtension->GetAOD()->FindListObject(fJetBranchName[1].Data())); } + Double_t ptsub[aodJets[0]->GetEntriesFast()]; + Int_t inord[aodJets[0]->GetEntriesFast()]; + for(Int_t n=0;nGetEntriesFast();n++){ + ptsub[n]=0; + inord[n]=0;} + TList ParticleList; Int_t nT = GetListOfTracks(&ParticleList); for (Int_t iJetType = 0; iJetType < 2; iJetType++) { @@ -752,12 +1323,13 @@ void AliAnalysisTaskJetCore::UserExec(Option_t *) if(fDebug) Printf("%s: %d jets",fJetBranchName[iJetType].Data(),aodJets[iJetType]->GetEntriesFast()); + for (Int_t iJet = 0; iJet < aodJets[iJetType]->GetEntriesFast(); iJet++) { AliAODJet *jet = dynamic_cast((*aodJets[iJetType])[iJet]); if (jet) fListJets[iJetType]->Add(jet); - } - fListJets[iJetType]->Sort(); - } + if(iJetType==0){ + ptsub[iJet]=jet->Pt()-rho*jet->EffectiveAreaCharged(); + }}} Double_t etabig=0; Double_t ptbig=0; @@ -768,6 +1340,27 @@ void AliAnalysisTaskJetCore::UserExec(Option_t *) Double_t areasmall=0; Double_t distr=0.; Double_t phismall=0.; + Int_t indexlead=-1; + Int_t indexsublead=-1; + Int_t indexstop=-1; + + TMath::Sort(fListJets[0]->GetEntries(),ptsub,inord); + + for(Int_t jj=0;jjGetEntries();jj++){ + AliAODJet* jetlead = (AliAODJet*)(fListJets[0]->At(inord[jj])); + if(jetlead->Pt()-rho*jetlead->EffectiveAreaCharged()<=0) continue; + if((jetlead->Eta()Eta()>fJetEtaMax)) continue; + indexlead=inord[jj]; + indexstop=jj; + break;} + if((indexstop>-1)&&(indexstop+1GetEntries()-1)){ + for(Int_t k=indexstop+1;kGetEntries();k++){ + AliAODJet* jetsublead = (AliAODJet*)(fListJets[0]->At(inord[k])); + if(jetsublead->Pt()-rho*jetsublead->EffectiveAreaCharged()<=0) continue; + if((jetsublead->Eta()Eta()>fJetEtaMax)) continue; + indexsublead=inord[k]; + break;}} + Double_t up1[15]={0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.}; Double_t up2[15]={0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.}; @@ -809,15 +1402,21 @@ void AliAnalysisTaskJetCore::UserExec(Option_t *) distr=TMath::Sqrt(distr); if(distr<=fRadioFrac){ fracin=fracin+genTrack->Pt();}} - if(centValue<10) fh2JetCoreMethod3C10->Fill(ptcorr,fracin/ptbig); + if(centValue<10) {fh2JetCoreMethod3C10->Fill(ptcorr,fracin/ptbig); + if(i==indexlead) fh2JetCoreMethod3C10lead->Fill(ptcorr,fracin/ptbig); + if(i==indexsublead) fh2JetCoreMethod3C10sublead->Fill(ptcorr,fracin/ptbig);} + - if((centValue>20)&&(centValue<40)) fh2JetCoreMethod3C20->Fill(ptcorr,fracin/ptbig); - - if((centValue>30)&&(centValue<60)) fh2JetCoreMethod3C30->Fill(ptcorr,fracin/ptbig); - - if(centValue>60) fh2JetCoreMethod3C60->Fill(ptcorr,fracin/ptbig); - - + if((centValue>20)&&(centValue<40)) {fh2JetCoreMethod3C20->Fill(ptcorr,fracin/ptbig); + if(i==indexlead) fh2JetCoreMethod3C20lead->Fill(ptcorr,fracin/ptbig); + if(i==indexsublead) fh2JetCoreMethod3C20sublead->Fill(ptcorr,fracin/ptbig);} + if((centValue>30)&&(centValue<60)){ fh2JetCoreMethod3C30->Fill(ptcorr,fracin/ptbig); + if(i==indexlead) fh2JetCoreMethod3C30lead->Fill(ptcorr,fracin/ptbig); + if(i==indexsublead) fh2JetCoreMethod3C30sublead->Fill(ptcorr,fracin/ptbig);} + if(centValue>60){fh2JetCoreMethod3C60->Fill(ptcorr,fracin/ptbig); + if(i==indexlead) fh2JetCoreMethod3C60lead->Fill(ptcorr,fracin/ptbig); + if(i==indexsublead) fh2JetCoreMethod3C60sublead->Fill(ptcorr,fracin/ptbig);} + for(Int_t j=0; jGetEntries(); ++j){ @@ -858,6 +1457,7 @@ void AliAnalysisTaskJetCore::UserExec(Option_t *) if((deltaR>=0.4)&&(deltaR<=0.6))sumPtIn=sumPtIn+part->Pt(); if((deltaR>=0.8)&&(deltaR<=1.))sumPtOut=sumPtOut+part->Pt(); if(centValue<10.){ + if((ptcorr>=70.)&&(ptcorr<85.)) {fh2DeltaRC10pt1->Fill(part->Pt(),deltaR); if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R))fh2DeltaEtaC10pt1->Fill(part->Pt(),deltaEta); if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC10pt1->Fill(part->Pt(),deltaPhi);} @@ -869,7 +1469,41 @@ void AliAnalysisTaskJetCore::UserExec(Option_t *) if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC10pt3->Fill(part->Pt(),deltaPhi);} if((ptcorr>=120.)&&(ptcorr<140.)) {fh2DeltaRC10pt4->Fill(part->Pt(),deltaR); if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC10pt4->Fill(part->Pt(),deltaEta); - if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC10pt4->Fill(part->Pt(),deltaPhi); }} + if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC10pt4->Fill(part->Pt(),deltaPhi); } + + if(i==indexlead){ + if((ptcorr>=70.)&&(ptcorr<85.)) {fh2DeltaRC10pt1lead->Fill(part->Pt(),deltaR); + if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R))fh2DeltaEtaC10pt1lead->Fill(part->Pt(),deltaEta); + if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC10pt1lead->Fill(part->Pt(),deltaPhi);} + if((ptcorr>=85.)&&(ptcorr<100.)) {fh2DeltaRC10pt2lead->Fill(part->Pt(),deltaR); + if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC10pt2lead->Fill(part->Pt(),deltaEta); + if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC10pt2lead->Fill(part->Pt(),deltaPhi); } + if((ptcorr>=100.)&&(ptcorr<120.)) {fh2DeltaRC10pt3lead->Fill(part->Pt(),deltaR); + if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC10pt3lead->Fill(part->Pt(),deltaEta); + if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC10pt3lead->Fill(part->Pt(),deltaPhi);} + if((ptcorr>=120.)&&(ptcorr<140.)) {fh2DeltaRC10pt4lead->Fill(part->Pt(),deltaR); + if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC10pt4lead->Fill(part->Pt(),deltaEta); + if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC10pt4lead->Fill(part->Pt(),deltaPhi); }} + + + if(i==indexsublead){ + if((ptcorr>=70.)&&(ptcorr<85.)) {fh2DeltaRC10pt1sublead->Fill(part->Pt(),deltaR); + if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R))fh2DeltaEtaC10pt1sublead->Fill(part->Pt(),deltaEta); + if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC10pt1sublead->Fill(part->Pt(),deltaPhi);} + if((ptcorr>=85.)&&(ptcorr<100.)) {fh2DeltaRC10pt2sublead->Fill(part->Pt(),deltaR); + if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC10pt2sublead->Fill(part->Pt(),deltaEta); + if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC10pt2sublead->Fill(part->Pt(),deltaPhi); } + if((ptcorr>=100.)&&(ptcorr<120.)) {fh2DeltaRC10pt3sublead->Fill(part->Pt(),deltaR); + if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC10pt3sublead->Fill(part->Pt(),deltaEta); + if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC10pt3sublead->Fill(part->Pt(),deltaPhi);} + if((ptcorr>=120.)&&(ptcorr<140.)) {fh2DeltaRC10pt4sublead->Fill(part->Pt(),deltaR); + if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC10pt4sublead->Fill(part->Pt(),deltaEta); + if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC10pt4sublead->Fill(part->Pt(),deltaPhi); }} + + + + + } if((centValue>20.)&&(centValue<40.)){ if((ptcorr>=70.)&&(ptcorr<85.)) {fh2DeltaRC20pt1->Fill(part->Pt(),deltaR); @@ -883,7 +1517,37 @@ void AliAnalysisTaskJetCore::UserExec(Option_t *) if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC20pt3->Fill(part->Pt(),deltaPhi);} if((ptcorr>=120.)&&(ptcorr<140.)) {fh2DeltaRC20pt4->Fill(part->Pt(),deltaR); if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC20pt4->Fill(part->Pt(),deltaEta); - if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC20pt4->Fill(part->Pt(),deltaPhi); }} + if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC20pt4->Fill(part->Pt(),deltaPhi); } + + if(i==indexlead){ + if((ptcorr>=70.)&&(ptcorr<85.)) {fh2DeltaRC20pt1lead->Fill(part->Pt(),deltaR); + if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R))fh2DeltaEtaC20pt1lead->Fill(part->Pt(),deltaEta); + if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC20pt1lead->Fill(part->Pt(),deltaPhi);} + if((ptcorr>=85.)&&(ptcorr<100.)) {fh2DeltaRC20pt2lead->Fill(part->Pt(),deltaR); + if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC20pt2lead->Fill(part->Pt(),deltaEta); + if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC20pt2lead->Fill(part->Pt(),deltaPhi); } + if((ptcorr>=100.)&&(ptcorr<120.)) {fh2DeltaRC20pt3lead->Fill(part->Pt(),deltaR); + if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC20pt3lead->Fill(part->Pt(),deltaEta); + if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC20pt3lead->Fill(part->Pt(),deltaPhi);} + if((ptcorr>=120.)&&(ptcorr<140.)) {fh2DeltaRC20pt4lead->Fill(part->Pt(),deltaR); + if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC20pt4lead->Fill(part->Pt(),deltaEta); + if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC20pt4lead->Fill(part->Pt(),deltaPhi); }} + + if(i==indexsublead){ + if((ptcorr>=70.)&&(ptcorr<85.)) {fh2DeltaRC20pt1sublead->Fill(part->Pt(),deltaR); + if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R))fh2DeltaEtaC20pt1sublead->Fill(part->Pt(),deltaEta); + if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC20pt1sublead->Fill(part->Pt(),deltaPhi);} + if((ptcorr>=85.)&&(ptcorr<100.)) {fh2DeltaRC20pt2sublead->Fill(part->Pt(),deltaR); + if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC20pt2sublead->Fill(part->Pt(),deltaEta); + if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC20pt2sublead->Fill(part->Pt(),deltaPhi); } + if((ptcorr>=100.)&&(ptcorr<120.)) {fh2DeltaRC20pt3sublead->Fill(part->Pt(),deltaR); + if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC20pt3sublead->Fill(part->Pt(),deltaEta); + if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC20pt3sublead->Fill(part->Pt(),deltaPhi);} + if((ptcorr>=120.)&&(ptcorr<140.)) {fh2DeltaRC20pt4sublead->Fill(part->Pt(),deltaR); + if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC20pt4sublead->Fill(part->Pt(),deltaEta); + if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC20pt4sublead->Fill(part->Pt(),deltaPhi); }} + + } if((centValue>30.)&&(centValue<60.)){ @@ -898,7 +1562,42 @@ void AliAnalysisTaskJetCore::UserExec(Option_t *) if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC30pt3->Fill(part->Pt(),deltaPhi);} if((ptcorr>=120.)&&(ptcorr<140.)) {fh2DeltaRC30pt4->Fill(part->Pt(),deltaR); if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC30pt4->Fill(part->Pt(),deltaEta); - if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC30pt4->Fill(part->Pt(),deltaPhi); }} + if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC30pt4->Fill(part->Pt(),deltaPhi); } + + if(i==indexlead){ + if((ptcorr>=70.)&&(ptcorr<85.)) {fh2DeltaRC30pt1lead->Fill(part->Pt(),deltaR); + if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R))fh2DeltaEtaC30pt1lead->Fill(part->Pt(),deltaEta); + if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC30pt1lead->Fill(part->Pt(),deltaPhi);} + if((ptcorr>=85.)&&(ptcorr<100.)) {fh2DeltaRC30pt2lead->Fill(part->Pt(),deltaR); + if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC30pt2lead->Fill(part->Pt(),deltaEta); + if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC30pt2lead->Fill(part->Pt(),deltaPhi); } + if((ptcorr>=100.)&&(ptcorr<120.)) {fh2DeltaRC30pt3lead->Fill(part->Pt(),deltaR); + if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC30pt3lead->Fill(part->Pt(),deltaEta); + if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC30pt3lead->Fill(part->Pt(),deltaPhi);} + if((ptcorr>=120.)&&(ptcorr<140.)) {fh2DeltaRC30pt4lead->Fill(part->Pt(),deltaR); + if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC30pt4lead->Fill(part->Pt(),deltaEta); + if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC30pt4lead->Fill(part->Pt(),deltaPhi); }} + + + if(i==indexsublead){ + if((ptcorr>=70.)&&(ptcorr<85.)) {fh2DeltaRC30pt1sublead->Fill(part->Pt(),deltaR); + if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R))fh2DeltaEtaC30pt1sublead->Fill(part->Pt(),deltaEta); + if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC30pt1sublead->Fill(part->Pt(),deltaPhi);} + if((ptcorr>=85.)&&(ptcorr<100.)) {fh2DeltaRC30pt2sublead->Fill(part->Pt(),deltaR); + if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC30pt2sublead->Fill(part->Pt(),deltaEta); + if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC30pt2sublead->Fill(part->Pt(),deltaPhi); } + if((ptcorr>=100.)&&(ptcorr<120.)) {fh2DeltaRC30pt3sublead->Fill(part->Pt(),deltaR); + if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC30pt3sublead->Fill(part->Pt(),deltaEta); + if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC30pt3sublead->Fill(part->Pt(),deltaPhi);} + if((ptcorr>=120.)&&(ptcorr<140.)) {fh2DeltaRC30pt4sublead->Fill(part->Pt(),deltaR); + if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC30pt4sublead->Fill(part->Pt(),deltaEta); + if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC30pt4sublead->Fill(part->Pt(),deltaPhi); }} + + + + + + } if(centValue>60.){ @@ -913,7 +1612,37 @@ void AliAnalysisTaskJetCore::UserExec(Option_t *) if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC60pt3->Fill(part->Pt(),deltaPhi);} if((ptcorr>=120.)&&(ptcorr<140.)) {fh2DeltaRC60pt4->Fill(part->Pt(),deltaR); if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC60pt4->Fill(part->Pt(),deltaEta); - if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC60pt4->Fill(part->Pt(),deltaPhi);}} + if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC60pt4->Fill(part->Pt(),deltaPhi);} + + if(i==indexlead){ + if((ptcorr>=70.)&&(ptcorr<85.)) {fh2DeltaRC60pt1lead->Fill(part->Pt(),deltaR); + if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R))fh2DeltaEtaC60pt1lead->Fill(part->Pt(),deltaEta); + if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC60pt1lead->Fill(part->Pt(),deltaPhi);} + if((ptcorr>=85.)&&(ptcorr<100.)) {fh2DeltaRC60pt2lead->Fill(part->Pt(),deltaR); + if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC60pt2lead->Fill(part->Pt(),deltaEta); + if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC60pt2lead->Fill(part->Pt(),deltaPhi); } + if((ptcorr>=100.)&&(ptcorr<120.)) {fh2DeltaRC60pt3lead->Fill(part->Pt(),deltaR); + if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC60pt3lead->Fill(part->Pt(),deltaEta); + if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC60pt3lead->Fill(part->Pt(),deltaPhi);} + if((ptcorr>=120.)&&(ptcorr<140.)) {fh2DeltaRC60pt4lead->Fill(part->Pt(),deltaR); + if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC60pt4lead->Fill(part->Pt(),deltaEta); + if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC60pt4lead->Fill(part->Pt(),deltaPhi); }} + + if(i==indexsublead){ + if((ptcorr>=70.)&&(ptcorr<85.)) {fh2DeltaRC60pt1sublead->Fill(part->Pt(),deltaR); + if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R))fh2DeltaEtaC60pt1sublead->Fill(part->Pt(),deltaEta); + if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC60pt1sublead->Fill(part->Pt(),deltaPhi);} + if((ptcorr>=85.)&&(ptcorr<100.)) {fh2DeltaRC60pt2sublead->Fill(part->Pt(),deltaR); + if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC60pt2sublead->Fill(part->Pt(),deltaEta); + if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC60pt2sublead->Fill(part->Pt(),deltaPhi); } + if((ptcorr>=100.)&&(ptcorr<120.)) {fh2DeltaRC60pt3sublead->Fill(part->Pt(),deltaR); + if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC60pt3sublead->Fill(part->Pt(),deltaEta); + if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC60pt3sublead->Fill(part->Pt(),deltaPhi);} + if((ptcorr>=120.)&&(ptcorr<140.)) {fh2DeltaRC60pt4sublead->Fill(part->Pt(),deltaR); + if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC60pt4sublead->Fill(part->Pt(),deltaEta); + if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC60pt4sublead->Fill(part->Pt(),deltaPhi); }} + + } } //end of track loop @@ -924,23 +1653,64 @@ void AliAnalysisTaskJetCore::UserExec(Option_t *) fh2SumPtOutC10bkg->Fill(ptcorr,coronaout/ptbig); fh2SumPtInC10->Fill(ptcorr,sumPtIn/ptbig); fh2SumPtOutC10->Fill(ptcorr,sumPtOut/ptbig); - fh2SumPtOutC10b->Fill(ptcorr,sumPtOut/ptbig); -} + if(i==indexlead){ fh2SumPtInC10bkglead->Fill(ptcorr,coronain/ptbig); + fh2SumPtOutC10bkglead->Fill(ptcorr,coronaout/ptbig); + fh2SumPtInC10lead->Fill(ptcorr,sumPtIn/ptbig); + fh2SumPtOutC10lead->Fill(ptcorr,sumPtOut/ptbig);} + if(i==indexsublead){ fh2SumPtInC10bkgsublead->Fill(ptcorr,coronain/ptbig); + fh2SumPtOutC10bkgsublead->Fill(ptcorr,coronaout/ptbig); + fh2SumPtInC10sublead->Fill(ptcorr,sumPtIn/ptbig); + fh2SumPtOutC10sublead->Fill(ptcorr,sumPtOut/ptbig);} + } if((centValue>20.)&&(centValue<40.)){ fh2SumPtInC20bkg->Fill(ptcorr,coronain/ptbig); fh2SumPtOutC20bkg->Fill(ptcorr,coronaout/ptbig); fh2SumPtInC20->Fill(ptcorr,sumPtIn/ptbig); - fh2SumPtOutC20->Fill(ptcorr,sumPtOut/ptbig); } + fh2SumPtOutC20->Fill(ptcorr,sumPtOut/ptbig); + if(i==indexlead){ fh2SumPtInC20bkglead->Fill(ptcorr,coronain/ptbig); + fh2SumPtOutC20bkglead->Fill(ptcorr,coronaout/ptbig); + fh2SumPtInC20lead->Fill(ptcorr,sumPtIn/ptbig); + fh2SumPtOutC20lead->Fill(ptcorr,sumPtOut/ptbig);} + if(i==indexsublead){ fh2SumPtInC20bkgsublead->Fill(ptcorr,coronain/ptbig); + fh2SumPtOutC20bkgsublead->Fill(ptcorr,coronaout/ptbig); + fh2SumPtInC20sublead->Fill(ptcorr,sumPtIn/ptbig); + fh2SumPtOutC20sublead->Fill(ptcorr,sumPtOut/ptbig);} + + + + } if((centValue>30.)&&(centValue<60.)){ fh2SumPtInC30bkg->Fill(ptcorr,coronain/ptbig); fh2SumPtOutC30bkg->Fill(ptcorr,coronaout/ptbig); fh2SumPtInC30->Fill(ptcorr,sumPtIn/ptbig); - fh2SumPtOutC30->Fill(ptcorr,sumPtOut/ptbig); } + fh2SumPtOutC30->Fill(ptcorr,sumPtOut/ptbig); + + if(i==indexlead){ fh2SumPtInC30bkglead->Fill(ptcorr,coronain/ptbig); + fh2SumPtOutC30bkglead->Fill(ptcorr,coronaout/ptbig); + fh2SumPtInC30lead->Fill(ptcorr,sumPtIn/ptbig); + fh2SumPtOutC30lead->Fill(ptcorr,sumPtOut/ptbig);} + if(i==indexsublead){ fh2SumPtInC30bkgsublead->Fill(ptcorr,coronain/ptbig); + fh2SumPtOutC30bkgsublead->Fill(ptcorr,coronaout/ptbig); + fh2SumPtInC30sublead->Fill(ptcorr,sumPtIn/ptbig); + fh2SumPtOutC30sublead->Fill(ptcorr,sumPtOut/ptbig);} + + } if(centValue>60.){ fh2SumPtInC60bkg->Fill(ptcorr,coronain/ptbig); fh2SumPtOutC60bkg->Fill(ptcorr,coronaout/ptbig); fh2SumPtInC60->Fill(ptcorr,sumPtIn/ptbig); - fh2SumPtOutC60->Fill(ptcorr,sumPtOut/ptbig);} + fh2SumPtOutC60->Fill(ptcorr,sumPtOut/ptbig); + + if(i==indexlead){ fh2SumPtInC60bkglead->Fill(ptcorr,coronain/ptbig); + fh2SumPtOutC60bkglead->Fill(ptcorr,coronaout/ptbig); + fh2SumPtInC60lead->Fill(ptcorr,sumPtIn/ptbig); + fh2SumPtOutC60lead->Fill(ptcorr,sumPtOut/ptbig);} + if(i==indexsublead){ fh2SumPtInC60bkgsublead->Fill(ptcorr,coronain/ptbig); + fh2SumPtOutC60bkgsublead->Fill(ptcorr,coronaout/ptbig); + fh2SumPtInC60sublead->Fill(ptcorr,sumPtIn/ptbig); + fh2SumPtOutC60sublead->Fill(ptcorr,sumPtOut/ptbig);} + + } //////////////////ANGULAR STRUCTURE////////////////////////////////////// @@ -972,7 +1742,7 @@ void AliAnalysisTaskJetCore::UserExec(Option_t *) } //only jet constituents - if(fAngStructCloseTracks==0){ + if(fAngStructCloseTracks==2){ Double_t difR=0.04; for(Int_t l=0;l<15;l++){ @@ -1018,7 +1788,7 @@ void AliAnalysisTaskJetCore::UserExec(Option_t *) //end loop over R=0.4 jets - + if(fAngStructCloseTracks>0){ for(Int_t l=0;l<15;l++){ Double_t rr=l*0.1+0.1; if(down1[l]!=0){ @@ -1040,7 +1810,7 @@ void AliAnalysisTaskJetCore::UserExec(Option_t *) if(centValue<10.) fh2AngStructpt4C10->Fill(rr,rr*up4[l]/down4[l]); if(centValue>20. && centValue<40.) fh2AngStructpt4C20->Fill(rr,rr*up4[l]/down4[l]); if(centValue>30. && centValue<60.) fh2AngStructpt4C30->Fill(rr,rr*up4[l]/down4[l]); - if(centValue>60.) fh2AngStructpt4C60->Fill(rr,rr*up4[l]/down4[l]);}} + if(centValue>60.) fh2AngStructpt4C60->Fill(rr,rr*up4[l]/down4[l]);}}} @@ -1049,7 +1819,7 @@ void AliAnalysisTaskJetCore::UserExec(Option_t *) PostData(1, fOutputList); -} + } void AliAnalysisTaskJetCore::Terminate(const Option_t *) { diff --git a/PWGJE/AliAnalysisTaskJetCore.h b/PWGJE/AliAnalysisTaskJetCore.h index 9596e0c68fb..efcdcf351f4 100644 --- a/PWGJE/AliAnalysisTaskJetCore.h +++ b/PWGJE/AliAnalysisTaskJetCore.h @@ -144,25 +144,64 @@ private: TH2F *fh2JetCoreMethod1C60; //Energy fraction in the core C60 method 1 TH2F *fh2JetCoreMethod2C60; //Energy fraction in the core C60 method 2 TH2F *fh2JetCoreMethod3C60; //Energy fraction in the core C60 method 3 + TH2F *fh2JetCoreMethod3C10lead; //Energy fraction in the core C30 method 3 + TH2F *fh2JetCoreMethod3C20lead; //Energy fraction in the core C60 method 1 + TH2F *fh2JetCoreMethod3C30lead; //Energy fraction in the core C60 method 2 + TH2F *fh2JetCoreMethod3C60lead; //Energy fraction in the core C60 method 3 + TH2F *fh2JetCoreMethod3C10sublead; //Energy fraction in the core C30 method 3 + TH2F *fh2JetCoreMethod3C20sublead; //Energy fraction in the core C60 method 1 + TH2F *fh2JetCoreMethod3C30sublead; //Energy fraction in the core C60 method 2 + TH2F *fh2JetCoreMethod3C60sublead; //Energy fraction in the core C60 method 3 + TH2F *fh2SumPtInC10; //energy fraction in inner corona C10 TH2F *fh2SumPtInC20; //energy fraction in inner corona C20 TH2F *fh2SumPtInC30; //energy fraction in inner corona C30 TH2F *fh2SumPtInC60; //energy fraction in inner corona C60 + TH2F *fh2SumPtInC10lead; //energy fraction in inner corona C10 leading + TH2F *fh2SumPtInC20lead; //energy fraction in inner corona C20 leading + TH2F *fh2SumPtInC30lead; //energy fraction in inner corona C30 leading + TH2F *fh2SumPtInC60lead; //energy fraction in inner corona C60 leading + TH2F *fh2SumPtInC10sublead; //energy fraction in inner corona C10 subleading + TH2F *fh2SumPtInC20sublead; //energy fraction in inner corona C20 subleading + TH2F *fh2SumPtInC30sublead; //energy fraction in inner corona C30 subleading + TH2F *fh2SumPtInC60sublead; //energy fraction in inner corona C60 subleading TH2F *fh2SumPtOutC10; //energy fraction in outer corona C10 - TH2F *fh2SumPtOutC10b; //energy fraction in outer corona C10b TH2F *fh2SumPtOutC20; //energy fraction in outer corona C20 TH2F *fh2SumPtOutC30; //energy fraction in outer corona C30 TH2F *fh2SumPtOutC60; //energy fraction in outer corona C60 + TH2F *fh2SumPtOutC10lead; //energy fraction in outer corona C10 leading + TH2F *fh2SumPtOutC20lead; //energy fraction in outer corona C20 leading + TH2F *fh2SumPtOutC30lead; //energy fraction in outer corona C30 leading + TH2F *fh2SumPtOutC60lead; //energy fraction in outer corona C60 leading + TH2F *fh2SumPtOutC10sublead; //energy fraction in outer corona C10 subleading + TH2F *fh2SumPtOutC20sublead; //energy fraction in outer corona C20 subleading + TH2F *fh2SumPtOutC30sublead; //energy fraction in outer corona C30 subleading + TH2F *fh2SumPtOutC60sublead; //energy fraction in outer corona C60 subleading TH2F *fh2SumPtInC10bkg; //expected from background inner C10 TH2F *fh2SumPtInC20bkg; //expected from background inner C20 TH2F *fh2SumPtInC30bkg; //expected from background inner C30 TH2F *fh2SumPtInC60bkg; //expected from background inner C60 + TH2F *fh2SumPtInC10bkglead; //expected from background inner C10 lead + TH2F *fh2SumPtInC20bkglead; //expected from background inner C20 lead + TH2F *fh2SumPtInC30bkglead; //expected from background inner C30 lead + TH2F *fh2SumPtInC60bkglead; //expected from background inner C60 lead + TH2F *fh2SumPtInC10bkgsublead; //expected from background inner C10 sublead + TH2F *fh2SumPtInC20bkgsublead; //expected from background inner C20 sublead + TH2F *fh2SumPtInC30bkgsublead; //expected from background inner C30 sublead + TH2F *fh2SumPtInC60bkgsublead; //expected from background inner C60 sublead + TH2F *fh2SumPtOutC10bkg; //expected from background outer C10 TH2F *fh2SumPtOutC20bkg; //expected from background outer C10 TH2F *fh2SumPtOutC30bkg; //expected from background outer C10 TH2F *fh2SumPtOutC60bkg; //expected from background outer C10 - - + TH2F *fh2SumPtOutC10bkglead; //expected from background outer C10 lead + TH2F *fh2SumPtOutC20bkglead; //expected from background outer C10 lead + TH2F *fh2SumPtOutC30bkglead; //expected from background outer C10 lead + TH2F *fh2SumPtOutC60bkglead; //expected from background outer C10 lead + TH2F *fh2SumPtOutC10bkgsublead; //expected from background outer C10 sublead + TH2F *fh2SumPtOutC20bkgsublead; //expected from background outer C10 sublead + TH2F *fh2SumPtOutC30bkgsublead; //expected from background outer C10 sublead + TH2F *fh2SumPtOutC60bkgsublead; //expected from background outer C10 sublead TH2F* fh2DeltaRC10pt1; //Jet track R distance:C10 pt1 TH2F* fh2DeltaRC20pt1; //C20 pt1 @@ -212,6 +251,109 @@ private: TH2F* fh2DeltaPhiC20pt4; //C20 pt4 TH2F* fh2DeltaPhiC30pt4; //C30 pt4 TH2F* fh2DeltaPhiC60pt4; //C60 pt4 + + TH2F* fh2DeltaRC10pt1lead; //Jet track R distance:C10 pt1 + TH2F* fh2DeltaRC20pt1lead; //C20 pt1 + TH2F* fh2DeltaRC30pt1lead; //C30 pt1 + TH2F* fh2DeltaRC60pt1lead; //C60 pt1 + TH2F* fh2DeltaRC10pt2lead; //C10 pt2 + TH2F* fh2DeltaRC20pt2lead; //C20 pt2 + TH2F* fh2DeltaRC30pt2lead; //C30 pt2 + TH2F* fh2DeltaRC60pt2lead; //C60 pt2 + TH2F* fh2DeltaRC10pt3lead; //C10 pt3 + TH2F* fh2DeltaRC20pt3lead; //C20 pt3 + TH2F* fh2DeltaRC30pt3lead; //C30 pt3 + TH2F* fh2DeltaRC60pt3lead; //C60 pt3 + TH2F* fh2DeltaRC10pt4lead; //C10 pt4 + TH2F* fh2DeltaRC20pt4lead; //C20 pt4 + TH2F* fh2DeltaRC30pt4lead; //C30 pt4 + TH2F* fh2DeltaRC60pt4lead; //C60 pt4 + TH2F* fh2DeltaEtaC10pt1lead; //The same but eta distance:C10 pt1 + TH2F* fh2DeltaEtaC20pt1lead; //C20 pt1 + TH2F* fh2DeltaEtaC30pt1lead; //C30 pt1 + TH2F* fh2DeltaEtaC60pt1lead; //C60 pt1 + TH2F* fh2DeltaEtaC10pt2lead; //C10 pt2 + TH2F* fh2DeltaEtaC20pt2lead; //C20 pt2 + TH2F* fh2DeltaEtaC30pt2lead; //C30 pt2 + TH2F* fh2DeltaEtaC60pt2lead; //C60 pt2 + TH2F* fh2DeltaEtaC10pt3lead; //C10 pt3 + TH2F* fh2DeltaEtaC20pt3lead; //C20 pt3 + TH2F* fh2DeltaEtaC30pt3lead; //C30 pt3 + TH2F* fh2DeltaEtaC60pt3lead; //C60 pt3 + TH2F* fh2DeltaEtaC10pt4lead; //C10 pt4 + TH2F* fh2DeltaEtaC20pt4lead; //C20 pt4 + TH2F* fh2DeltaEtaC30pt4lead; //C30 pt4 + TH2F* fh2DeltaEtaC60pt4lead; //C60 pt4 + TH2F* fh2DeltaPhiC10pt1lead; //The same but phi distance:C10 pt1 + TH2F* fh2DeltaPhiC20pt1lead; //C20 pt1 + TH2F* fh2DeltaPhiC30pt1lead; //C30 pt1 + TH2F* fh2DeltaPhiC60pt1lead; //C60 pt1 + TH2F* fh2DeltaPhiC10pt2lead; //C10 pt2 + TH2F* fh2DeltaPhiC20pt2lead; //C20 pt2 + TH2F* fh2DeltaPhiC30pt2lead; //C30 pt2 + TH2F* fh2DeltaPhiC60pt2lead; //C60 pt2 + TH2F* fh2DeltaPhiC10pt3lead; //C10 pt3 + TH2F* fh2DeltaPhiC20pt3lead; //C20 pt3 + TH2F* fh2DeltaPhiC30pt3lead; //C30 pt3 + TH2F* fh2DeltaPhiC60pt3lead; //C60 pt3 + TH2F* fh2DeltaPhiC10pt4lead; //C10 pt4 + TH2F* fh2DeltaPhiC20pt4lead; //C20 pt4 + TH2F* fh2DeltaPhiC30pt4lead; //C30 pt4 + TH2F* fh2DeltaPhiC60pt4lead; //C60 pt4 + + TH2F* fh2DeltaRC10pt1sublead; //Jet track R distance:C10 pt1 + TH2F* fh2DeltaRC20pt1sublead; //C20 pt1 + TH2F* fh2DeltaRC30pt1sublead; //C30 pt1 + TH2F* fh2DeltaRC60pt1sublead; //C60 pt1 + TH2F* fh2DeltaRC10pt2sublead; //C10 pt2 + TH2F* fh2DeltaRC20pt2sublead; //C20 pt2 + TH2F* fh2DeltaRC30pt2sublead; //C30 pt2 + TH2F* fh2DeltaRC60pt2sublead; //C60 pt2 + TH2F* fh2DeltaRC10pt3sublead; //C10 pt3 + TH2F* fh2DeltaRC20pt3sublead; //C20 pt3 + TH2F* fh2DeltaRC30pt3sublead; //C30 pt3 + TH2F* fh2DeltaRC60pt3sublead; //C60 pt3 + TH2F* fh2DeltaRC10pt4sublead; //C10 pt4 + TH2F* fh2DeltaRC20pt4sublead; //C20 pt4 + TH2F* fh2DeltaRC30pt4sublead; //C30 pt4 + TH2F* fh2DeltaRC60pt4sublead; //C60 pt4 + TH2F* fh2DeltaEtaC10pt1sublead; //The same but eta distance:C10 pt1 + TH2F* fh2DeltaEtaC20pt1sublead; //C20 pt1 + TH2F* fh2DeltaEtaC30pt1sublead; //C30 pt1 + TH2F* fh2DeltaEtaC60pt1sublead; //C60 pt1 + TH2F* fh2DeltaEtaC10pt2sublead; //C10 pt2 + TH2F* fh2DeltaEtaC20pt2sublead; //C20 pt2 + TH2F* fh2DeltaEtaC30pt2sublead; //C30 pt2 + TH2F* fh2DeltaEtaC60pt2sublead; //C60 pt2 + TH2F* fh2DeltaEtaC10pt3sublead; //C10 pt3 + TH2F* fh2DeltaEtaC20pt3sublead; //C20 pt3 + TH2F* fh2DeltaEtaC30pt3sublead; //C30 pt3 + TH2F* fh2DeltaEtaC60pt3sublead; //C60 pt3 + TH2F* fh2DeltaEtaC10pt4sublead; //C10 pt4 + TH2F* fh2DeltaEtaC20pt4sublead; //C20 pt4 + TH2F* fh2DeltaEtaC30pt4sublead; //C30 pt4 + TH2F* fh2DeltaEtaC60pt4sublead; //C60 pt4 + TH2F* fh2DeltaPhiC10pt1sublead; //The same but phi distance:C10 pt1 + TH2F* fh2DeltaPhiC20pt1sublead; //C20 pt1 + TH2F* fh2DeltaPhiC30pt1sublead; //C30 pt1 + TH2F* fh2DeltaPhiC60pt1sublead; //C60 pt1 + TH2F* fh2DeltaPhiC10pt2sublead; //C10 pt2 + TH2F* fh2DeltaPhiC20pt2sublead; //C20 pt2 + TH2F* fh2DeltaPhiC30pt2sublead; //C30 pt2 + TH2F* fh2DeltaPhiC60pt2sublead; //C60 pt2 + TH2F* fh2DeltaPhiC10pt3sublead; //C10 pt3 + TH2F* fh2DeltaPhiC20pt3sublead; //C20 pt3 + TH2F* fh2DeltaPhiC30pt3sublead; //C30 pt3 + TH2F* fh2DeltaPhiC60pt3sublead; //C60 pt3 + TH2F* fh2DeltaPhiC10pt4sublead; //C10 pt4 + TH2F* fh2DeltaPhiC20pt4sublead; //C20 pt4 + TH2F* fh2DeltaPhiC30pt4sublead; //C30 pt4 + TH2F* fh2DeltaPhiC60pt4sublead; //C60 pt4 + + + + + TH2F* fh2AngStructpt1C10; //Average two particle correlation function:C10 pt1 TH2F* fh2AngStructpt2C10; //C10 pt2 TH2F* fh2AngStructpt3C10; //C10 pt3 diff --git a/PWGJE/macros/AddTaskJetCore.C b/PWGJE/macros/AddTaskJetCore.C index b19dda88572..28f6d335562 100644 --- a/PWGJE/macros/AddTaskJetCore.C +++ b/PWGJE/macros/AddTaskJetCore.C @@ -1,6 +1,6 @@ -AliAnalysisTaskJetCore* AddTaskJetCore(const char* bRec1,const char* bRec2, UInt_t filterMask = 272 , Float_t ptTrackMin = 0.15, Int_t eventClassMin = 0, Int_t eventClassMax = 4){ +AliAnalysisTaskJetCore* AddTaskJetCore(const char* bRec1,const char* bRec2, UInt_t filterMask = 272 , Float_t ptTrackMin = 0.15, Int_t kTriggerMask=0, Int_t eventClassMin = 0, Int_t eventClassMax = 4){ Printf("adding task jet response\n"); @@ -20,7 +20,7 @@ AliAnalysisTaskJetCore* AddTaskJetCore(const char* bRec1,const char* bRec2, UInt TString typeRec(bRec1); TString typeGen(bRec2); - AliAnalysisTaskJetCore *task = new AliAnalysisTaskJetCore(Form("JetCore_%s_%s",bRec1,bRec2)); + AliAnalysisTaskJetCore *task = new AliAnalysisTaskJetCore(Form("JetCore_%s_%s_%d",bRec1,bRec2,kTriggerMask)); @@ -42,7 +42,7 @@ AliAnalysisTaskJetCore* AddTaskJetCore(const char* bRec1,const char* bRec2, UInt mgr->AddTask(task); - AliAnalysisDataContainer *coutputJetCore = mgr->CreateContainer(Form("pwgjejetcore_%s_%s",bRec1,bRec2), TList::Class(),AliAnalysisManager::kOutputContainer,Form("%s:PWGJE_jetcore_%s_%s",AliAnalysisManager::GetCommonFileName(),bRec1,bRec2)); + AliAnalysisDataContainer *coutputJetCore = mgr->CreateContainer(Form("pwgjejetcore_%s_%s_%d",bRec1,bRec2,kTriggerMask), TList::Class(),AliAnalysisManager::kOutputContainer,Form("%s:PWGJE_jetcore_%s_%s_%d",AliAnalysisManager::GetCommonFileName(),bRec1,bRec2,kTriggerMask)); -- 2.31.1