changed binning, vertex range, leading, subleading

diff --git a/PWGJE/AliAnalysisTaskJetCore.cxx b/PWGJE/AliAnalysisTaskJetCore.cxx
index edc2451c42b95a7a289746539ed1ee24f5e9112b..ef551ab170d44484d7f2d058871a47ef7f4c02a3 100644 (file)
--- 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<TClonesArray*>(fAOD->FindListObject(fJetBranchName[0].Data())); // 
+   aodJets[0] = dynamic_cast<TClonesArray*>(fAOD->FindListObject(fJetBranchName[0].Data())); 
    aodJets[1] = dynamic_cast<TClonesArray*>(fAOD->FindListObject(fJetBranchName[1].Data()));  }
    if(fAODExtension && !aodJets[0]){ 
-     aodJets[0] = dynamic_cast<TClonesArray*>(fAODExtension->GetAOD()->FindListObject(fJetBranchName[0].Data())); // 
+     aodJets[0] = dynamic_cast<TClonesArray*>(fAODExtension->GetAOD()->FindListObject(fJetBranchName[0].Data())); 
      aodJets[1] = dynamic_cast<TClonesArray*>(fAODExtension->GetAOD()->FindListObject(fJetBranchName[1].Data()));  }
 
+   Double_t ptsub[aodJets[0]->GetEntriesFast()];
+   Int_t inord[aodJets[0]->GetEntriesFast()];
+   for(Int_t n=0;n<aodJets[0]->GetEntriesFast();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<AliAODJet*>((*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;jj<fListJets[0]->GetEntries();jj++){
+   AliAODJet* jetlead = (AliAODJet*)(fListJets[0]->At(inord[jj]));
+   if(jetlead->Pt()-rho*jetlead->EffectiveAreaCharged()<=0) continue;
+   if((jetlead->Eta()<fJetEtaMin)||(jetlead->Eta()>fJetEtaMax)) continue;
+   indexlead=inord[jj];
+   indexstop=jj;
+   break;}
+   if((indexstop>-1)&&(indexstop+1<fListJets[0]->GetEntries()-1)){
+   for(Int_t k=indexstop+1;k<fListJets[0]->GetEntries();k++){
+   AliAODJet* jetsublead = (AliAODJet*)(fListJets[0]->At(inord[k]));
+   if(jetsublead->Pt()-rho*jetsublead->EffectiveAreaCharged()<=0) continue;
+   if((jetsublead->Eta()<fJetEtaMin)||(jetsublead->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; j<fListJets[1]->GetEntries(); ++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 9596e0c68fb4f17587b9ee5b27d5b58771b17de0..efcdcf351f497989cb194947b6aef80ce856660a 100644 (file)
--- 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 b19dda885721fe6f4de052b568c8287fec8dbd21..28f6d335562391dc8273580e1a1604ec0827d6b1 100644 (file)
--- 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));