changed binning, vertex range, leading, subleading
[u/mrichter/AliRoot.git] / PWGJE / AliAnalysisTaskJetCore.cxx
index c418e3d4d5c73d13e1e3fb964d803fa90b04189e..ef551ab170d44484d7f2d058871a47ef7f4c02a3 100644 (file)
@@ -1,3 +1,29 @@
+
+// ******************************************
+// This task computes several jet observables like 
+// the fraction of energy in inner and outer coronnas,
+// the distance from track to jet axis and a 
+// correlation strength distribution of particles inside jets.    
+// Author: lcunquei@cern.ch
+// *******************************************
+
+
+/**************************************************************************
+ * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ *                                                                        *
+ * Author: The ALICE Off-line Project.                                    *
+ * Contributors are mentioned in the code where appropriate.              *
+ *                                                                        *
+ * Permission to use, copy, modify and distribute this software and its   *
+ * documentation strictly for non-commercial purposes is hereby granted   *
+ * without fee, provided that the above copyright notice appears in all   *
+ * copies and that both the copyright notice and this permission notice   *
+ * appear in the supporting documentation. The authors make no claims     *
+ * about the suitability of this software for any purpose. It is          *
+ * provided "as is" without express or implied warranty.                  *
+ **************************************************************************/
+
+
 #include "TChain.h"
 #include "TTree.h"
 #include "TMath.h"
@@ -39,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),
@@ -77,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),
@@ -142,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),
@@ -178,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),
@@ -216,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),
@@ -281,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),
@@ -395,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.); 
 
 
 
@@ -433,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); 
@@ -467,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);
@@ -499,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);
@@ -533,7 +1020,6 @@ void AliAnalysisTaskJetCore::UserCreateOutputObjects()
       fOutputList->Add(fh2DeltaRC20pt4);
       fOutputList->Add(fh2DeltaRC30pt4);
       fOutputList->Add(fh2DeltaRC60pt4);
-      
       fOutputList->Add(fh2DeltaEtaC10pt1);
       fOutputList->Add(fh2DeltaEtaC20pt1);
       fOutputList->Add(fh2DeltaEtaC30pt1);
@@ -566,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);
@@ -582,7 +1173,7 @@ void AliAnalysisTaskJetCore::UserCreateOutputObjects()
        fOutputList->Add(fh2AngStructpt1C60);
        fOutputList->Add(fh2AngStructpt2C60);
        fOutputList->Add(fh2AngStructpt3C60);
-       fOutputList->Add(fh2AngStructpt4C60);  
+       fOutputList->Add(fh2AngStructpt4C60);}  
 
 
 
@@ -689,7 +1280,8 @@ void AliAnalysisTaskJetCore::UserExec(Option_t *)
    // }
 
    
-   fHistEvtSelection->Fill(0); // accepted events  
+   fHistEvtSelection->Fill(0); 
+   // accepted events  
    // -- end event selection --
   
    // get background
@@ -711,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++) {
@@ -725,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;
@@ -741,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.};
@@ -782,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){
@@ -831,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);}
@@ -842,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);
@@ -856,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.)){
 
@@ -871,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.){
@@ -886,9 +1612,40 @@ 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
+     } 
+     //end of track loop
      Double_t coronain=rho*TMath::Pi()*(1.-0.8*0.8);
      Double_t coronaout=rho*TMath::Pi()*(0.6*0.6-0.4*0.4);
      if(centValue<10.){  
@@ -896,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//////////////////////////////////////
 
@@ -944,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++){
@@ -955,7 +1753,7 @@ void AliAnalysisTaskJetCore::UserExec(Option_t *)
       AliAODTrack* part2;
           for(Int_t it=0; it<nTracksGenJet; ++it){
              part1 = (AliAODTrack*)(genTrackList->At(it));
-           for(Int_t itu=it+1; itu<nTracksGenJet; ++itu){
+           for(Int_t itu=0; itu<nTracksGenJet; ++itu){
              part2 = (AliAODTrack*)(genTrackList->At(itu));
            Double_t ptm=part1->Pt();
            Double_t ptn=part2->Pt();   
@@ -990,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){  
@@ -1012,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]);}}}
 
     
 
@@ -1021,7 +1819,7 @@ void AliAnalysisTaskJetCore::UserExec(Option_t *)
 
 
    PostData(1, fOutputList);
-}
+     }
 
 void AliAnalysisTaskJetCore::Terminate(const Option_t *)
 {