Loop only over MC true tracks in case of MC information

[u/mrichter/AliRoot.git] / PWGJE / AliAnalysisTaskJetCore.cxx

diff --git a/PWGJE/AliAnalysisTaskJetCore.cxx b/PWGJE/AliAnalysisTaskJetCore.cxx
index 59108af07bf2418ce18fbf04cdf861060085d6cd..a403ea616e132500c1fd18a9b870c4179bc6cccf 100644 (file)
--- a/PWGJE/AliAnalysisTaskJetCore.cxx
+++ b/PWGJE/AliAnalysisTaskJetCore.cxx
@@ -1,8 +1,7 @@
-
 // ******************************************
 // 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 
+// jet-track correlations,triggered jet shapes and 
 // correlation strength distribution of particles inside jets.    
 // Author: lcunquei@cern.ch
 // *******************************************
@@ -32,7 +31,7 @@
 #include "TH3F.h"
 #include "THnSparse.h"
 #include "TCanvas.h"
-
+#include "TRandom3.h"
 #include "AliLog.h"
 
 #include "AliAnalysisTask.h"
@@ -45,20 +44,24 @@
 #include "AliAnalysisHelperJetTasks.h"
 #include "AliInputEventHandler.h"
 #include "AliAODJetEventBackground.h"
+#include "AliAODMCParticle.h"
 #include "AliAnalysisTaskFastEmbedding.h"
-
 #include "AliAODEvent.h"
 #include "AliAODHandler.h"
 #include "AliAODJet.h"
 
 #include "AliAnalysisTaskJetCore.h"
 
+using std::cout;
+using std::endl;
+
 ClassImp(AliAnalysisTaskJetCore)
 
 AliAnalysisTaskJetCore::AliAnalysisTaskJetCore() :
 AliAnalysisTaskSE(),
 fESD(0x0),
-fAOD(0x0),
+fAODIn(0x0),
+fAODOut(0x0),
 fAODExtension(0x0),
 fBackgroundBranch(""),
 fNonStdFile(""),
@@ -70,240 +73,67 @@ fVtxZMax(10.),
 fEvtClassMin(0),
 fEvtClassMax(4),
 fFilterMask(0),
+fFilterMaskBestPt(0),
+fFilterType(0),
 fRadioFrac(0.2),
 fMinDist(0.1),
 fCentMin(0.),
 fCentMax(100.),
 fNInputTracksMin(0),
 fNInputTracksMax(-1),
+fRequireITSRefit(0),
+fApplySharedClusterCut(0),
 fAngStructCloseTracks(0),
+fCheckMethods(0),
+fDoEventMixing(0), 
+fFlagPhiBkg(0),
+fFlagEtaBkg(0),
+fFlagJetHadron(0),
+fFrac(0.8),
+fTTLowRef(11),
+fTTUpRef(13),
+fTTLowSig(15),
+fTTUpSig(19),
+fHardest(0),
+fFlagRandom(0),
+fFlagOnlyRecoil(0),
+fFlagOnlyHardest(1),
+fTrackTypeRec(kTrackUndef),
+fRPAngle(0),
+fNRPBins(50),
+fSemigoodCorrect(0),
+fHolePos(4.71),
+fHoleWidth(0.2),
 fJetEtaMin(-.5),
 fJetEtaMax(.5),
+fNevents(0),
+fTindex(0),
+fTrigBufferIndex(0),
+fCountAgain(0), 
 fJetPtMin(20.),
 fJetTriggerExcludeMask(AliAODJet::kHighTrackPtTriggered),
 fJetPtFractionMin(0.5),
 fNMatchJets(4),
 fMatchMaxDist(0.8),
 fKeepJets(kFALSE),
+fRunAnaAzimuthalCorrelation(kFALSE),
 fkNbranches(2),
 fkEvtClasses(12),
 fOutputList(0x0),
 fbEvent(kTRUE),
 fHistEvtSelection(0x0),
-fHistJetSelection(0x0),
-fh2JetSelection(0x0),
+fhnDeltaR(0x0),
+fhnMixedEvents(0x0),
 fh2JetCoreMethod1C10(0x0),
 fh2JetCoreMethod2C10(0x0),
-fh2JetCoreMethod3C10(0x0),
 fh2JetCoreMethod1C20(0x0),
 fh2JetCoreMethod2C20(0x0),
-fh2JetCoreMethod3C20(0x0),
 fh2JetCoreMethod1C30(0x0),
 fh2JetCoreMethod2C30(0x0),
-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),
-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),
-fh2DeltaRC60pt1(0x0),  
-fh2DeltaRC10pt2(0x0),
-fh2DeltaRC20pt2(0x0),
-fh2DeltaRC30pt2(0x0),
-fh2DeltaRC60pt2(0x0),  
-fh2DeltaRC10pt3(0x0),
-fh2DeltaRC20pt3(0x0),
-fh2DeltaRC30pt3(0x0),
-fh2DeltaRC60pt3(0x0),  
-fh2DeltaRC10pt4(0x0),
-fh2DeltaRC20pt4(0x0),
-fh2DeltaRC30pt4(0x0),
-fh2DeltaRC60pt4(0x0),
-fh2DeltaEtaC10pt1(0x0),
-fh2DeltaEtaC20pt1(0x0),
-fh2DeltaEtaC30pt1(0x0),
-fh2DeltaEtaC60pt1(0x0),  
-fh2DeltaEtaC10pt2(0x0),
-fh2DeltaEtaC20pt2(0x0),
-fh2DeltaEtaC30pt2(0x0),
-fh2DeltaEtaC60pt2(0x0),  
-fh2DeltaEtaC10pt3(0x0),
-fh2DeltaEtaC20pt3(0x0),
-fh2DeltaEtaC30pt3(0x0),
-fh2DeltaEtaC60pt3(0x0),  
-fh2DeltaEtaC10pt4(0x0),
-fh2DeltaEtaC20pt4(0x0),
-fh2DeltaEtaC30pt4(0x0),
-fh2DeltaEtaC60pt4(0x0),
-fh2DeltaPhiC10pt1(0x0),
-fh2DeltaPhiC20pt1(0x0),
-fh2DeltaPhiC30pt1(0x0),
-fh2DeltaPhiC60pt1(0x0),  
-fh2DeltaPhiC10pt2(0x0),
-fh2DeltaPhiC20pt2(0x0),
-fh2DeltaPhiC30pt2(0x0),
-fh2DeltaPhiC60pt2(0x0),  
-fh2DeltaPhiC10pt3(0x0),
-fh2DeltaPhiC20pt3(0x0),
-fh2DeltaPhiC30pt3(0x0),
-fh2DeltaPhiC60pt3(0x0),  
-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),
+fh3JetTrackC3060(0x0),
+fh3JetTrackC20(0x0),
 fh2AngStructpt1C10(0x0),
 fh2AngStructpt2C10(0x0),
 fh2AngStructpt3C10(0x0),
@@ -319,10 +149,38 @@ fh2AngStructpt4C30(0x0),
 fh2AngStructpt1C60(0x0),
 fh2AngStructpt2C60(0x0),
 fh2AngStructpt3C60(0x0),
-fh2AngStructpt4C60(0x0)
+fh2AngStructpt4C60(0x0),
+fh1TrigRef(0x0),
+fh1TrigSig(0x0),
+fh2Ntriggers(0x0),
+fh2Ntriggers2C10(0x0),
+fh2Ntriggers2C20(0x0), 
+fh3JetDensity(0x0),
+fh3JetDensityA4(0x0),
+fh2RPJetsC10(0x0),
+fh2RPJetsC20(0x0),
+fh2RPTC10(0x0),
+fh2RPTC20(0x0), 
+fHJetSpec(0x0),
+fhTTPt(0x0),
+fHJetPhiCorr(0x0),
+fRandom(0x0)
+ 
 {
    // default Constructor
 
+
+ // Trigger buffer.
+   for(Int_t i=0; i<10; i++) {
+               for(Int_t j=0; j<6; j++) {
+                       fTrigBuffer[i][j]=0;
+               }                               
+       }       
+
+
+
+
+
    fJetBranchName[0] = "";
    fJetBranchName[1] = "";
 
@@ -333,7 +191,8 @@ fh2AngStructpt4C60(0x0)
 AliAnalysisTaskJetCore::AliAnalysisTaskJetCore(const char *name) :
 AliAnalysisTaskSE(name),
 fESD(0x0),
-fAOD(0x0),
+fAODIn(0x0),
+fAODOut(0x0),
 fAODExtension(0x0),
 fBackgroundBranch(""),
 fNonStdFile(""),
@@ -345,240 +204,67 @@ fVtxZMax(10.),
 fEvtClassMin(0),
 fEvtClassMax(4),
 fFilterMask(0),
+fFilterMaskBestPt(0),
+fFilterType(0),
 fRadioFrac(0.2),
 fMinDist(0.1),
 fCentMin(0.),
 fCentMax(100.),
 fNInputTracksMin(0),
 fNInputTracksMax(-1),
+fRequireITSRefit(0),
+fApplySharedClusterCut(0),
 fAngStructCloseTracks(0),
+fCheckMethods(0),
+fDoEventMixing(0),
+fFlagPhiBkg(0),
+fFlagEtaBkg(0),
+fFlagJetHadron(0),
+fFrac(0.8),
+fTTLowRef(11),
+fTTUpRef(13),
+fTTLowSig(15),
+fTTUpSig(19),
+fHardest(0),
+fFlagRandom(0),
+fFlagOnlyRecoil(0),
+fFlagOnlyHardest(1),
+fTrackTypeRec(kTrackUndef),
+fRPAngle(0),
+fNRPBins(50),
+fSemigoodCorrect(0),
+fHolePos(4.71),
+fHoleWidth(0.2),
 fJetEtaMin(-.5),
 fJetEtaMax(.5),
+fNevents(0),
+fTindex(0),
+fTrigBufferIndex(0),
+fCountAgain(0),
 fJetPtMin(20.),
 fJetTriggerExcludeMask(AliAODJet::kHighTrackPtTriggered),
 fJetPtFractionMin(0.5),
 fNMatchJets(4),
 fMatchMaxDist(0.8),
 fKeepJets(kFALSE),
+fRunAnaAzimuthalCorrelation(kFALSE),
 fkNbranches(2),
 fkEvtClasses(12),
 fOutputList(0x0),
 fbEvent(kTRUE),
 fHistEvtSelection(0x0),
-fHistJetSelection(0x0),
-fh2JetSelection(0x0),
+fhnDeltaR(0x0),
+fhnMixedEvents(0x0),
 fh2JetCoreMethod1C10(0x0),
 fh2JetCoreMethod2C10(0x0),
-fh2JetCoreMethod3C10(0x0),
 fh2JetCoreMethod1C20(0x0),
 fh2JetCoreMethod2C20(0x0),
-fh2JetCoreMethod3C20(0x0),
 fh2JetCoreMethod1C30(0x0),
 fh2JetCoreMethod2C30(0x0),
-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),
-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),
-fh2DeltaRC60pt1(0x0),  
-fh2DeltaRC10pt2(0x0),
-fh2DeltaRC20pt2(0x0),
-fh2DeltaRC30pt2(0x0),
-fh2DeltaRC60pt2(0x0),  
-fh2DeltaRC10pt3(0x0),
-fh2DeltaRC20pt3(0x0),
-fh2DeltaRC30pt3(0x0),
-fh2DeltaRC60pt3(0x0),  
-fh2DeltaRC10pt4(0x0),
-fh2DeltaRC20pt4(0x0),
-fh2DeltaRC30pt4(0x0),
-fh2DeltaRC60pt4(0x0),
-fh2DeltaEtaC10pt1(0x0),
-fh2DeltaEtaC20pt1(0x0),
-fh2DeltaEtaC30pt1(0x0),
-fh2DeltaEtaC60pt1(0x0),  
-fh2DeltaEtaC10pt2(0x0),
-fh2DeltaEtaC20pt2(0x0),
-fh2DeltaEtaC30pt2(0x0),
-fh2DeltaEtaC60pt2(0x0),  
-fh2DeltaEtaC10pt3(0x0),
-fh2DeltaEtaC20pt3(0x0),
-fh2DeltaEtaC30pt3(0x0),
-fh2DeltaEtaC60pt3(0x0),  
-fh2DeltaEtaC10pt4(0x0),
-fh2DeltaEtaC20pt4(0x0),
-fh2DeltaEtaC30pt4(0x0),
-fh2DeltaEtaC60pt4(0x0),
-fh2DeltaPhiC10pt1(0x0),
-fh2DeltaPhiC20pt1(0x0),
-fh2DeltaPhiC30pt1(0x0),
-fh2DeltaPhiC60pt1(0x0),  
-fh2DeltaPhiC10pt2(0x0),
-fh2DeltaPhiC20pt2(0x0),
-fh2DeltaPhiC30pt2(0x0),
-fh2DeltaPhiC60pt2(0x0),  
-fh2DeltaPhiC10pt3(0x0),
-fh2DeltaPhiC20pt3(0x0),
-fh2DeltaPhiC30pt3(0x0),
-fh2DeltaPhiC60pt3(0x0),  
-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),
+fh3JetTrackC3060(0x0),
+fh3JetTrackC20(0x0),
 fh2AngStructpt1C10(0x0),
 fh2AngStructpt2C10(0x0),
 fh2AngStructpt3C10(0x0),
@@ -594,11 +280,35 @@ fh2AngStructpt4C30(0x0),
 fh2AngStructpt1C60(0x0),
 fh2AngStructpt2C60(0x0),
 fh2AngStructpt3C60(0x0),
-fh2AngStructpt4C60(0x0)    
+fh2AngStructpt4C60(0x0),    
+fh1TrigRef(0x0),
+fh1TrigSig(0x0),
+fh2Ntriggers(0x0),
+fh2Ntriggers2C10(0x0),
+fh2Ntriggers2C20(0x0),
+fh3JetDensity(0x0),
+fh3JetDensityA4(0x0),
+fh2RPJetsC10(0x0),
+fh2RPJetsC20(0x0),
+fh2RPTC10(0x0),
+fh2RPTC20(0x0), 
+fHJetSpec(0x0),
+fhTTPt(0x0),
+fHJetPhiCorr(0x0),
+fRandom(0x0)
 
  {
    // Constructor
 
+
+    for(Int_t i=0; i<10; i++) {
+       for(Int_t j=0; j<6; j++) {
+           fTrigBuffer[i][j]=0;
+               }                               
+    }  
+
+
+
    fJetBranchName[0] = "";
    fJetBranchName[1] = "";
 
@@ -642,6 +352,12 @@ void AliAnalysisTaskJetCore::UserCreateOutputObjects()
    TH1::AddDirectory(kFALSE);
 
 
+        // set seed
+   fRandom = new TRandom3(0);
+  
+
+
+
    fHistEvtSelection = new TH1I("fHistEvtSelection", "event selection", 6, -0.5, 5.5);
    fHistEvtSelection->GetXaxis()->SetBinLabel(1,"ACCEPTED");
    fHistEvtSelection->GetXaxis()->SetBinLabel(2,"events IN");
@@ -650,267 +366,56 @@ void AliAnalysisTaskJetCore::UserCreateOutputObjects()
    fHistEvtSelection->GetXaxis()->SetBinLabel(5,"centrality (rejected)");
    fHistEvtSelection->GetXaxis()->SetBinLabel(6,"multiplicity (rejected)");
 
-   fHistJetSelection = new TH1I("fHistJetSelection", "jet selection", 8, -0.5, 7.5);
-   fHistJetSelection->GetXaxis()->SetBinLabel(1,"ACCEPTED");
-   fHistJetSelection->GetXaxis()->SetBinLabel(2,"probes IN");
-   fHistJetSelection->GetXaxis()->SetBinLabel(3,"no matching jet");
-   fHistJetSelection->GetXaxis()->SetBinLabel(4,"not in list");
-   fHistJetSelection->GetXaxis()->SetBinLabel(5,"fraction cut");
-   fHistJetSelection->GetXaxis()->SetBinLabel(6,"acceptance cut");
-   fHistJetSelection->GetXaxis()->SetBinLabel(7,"p_{T} cut");
-   fHistJetSelection->GetXaxis()->SetBinLabel(8,"trigger exclude mask");
-
-   fh2JetSelection = new TH2F("fh2JetSelection", "jet selection", 8, -0.5, 7.5,100,0.,200.);
-   fh2JetSelection->GetXaxis()->SetBinLabel(1,"ACCEPTED");
-   fh2JetSelection->GetXaxis()->SetBinLabel(2,"probes IN");
-   fh2JetSelection->GetXaxis()->SetBinLabel(3,"no matching jet");
-   fh2JetSelection->GetXaxis()->SetBinLabel(4,"not in list");
-   fh2JetSelection->GetXaxis()->SetBinLabel(5,"fraction cut");
-   fh2JetSelection->GetXaxis()->SetBinLabel(6,"acceptance cut");
-   fh2JetSelection->GetXaxis()->SetBinLabel(7,"p_{T} cut");
-   fh2JetSelection->GetXaxis()->SetBinLabel(8,"trigger exclude mask");
+     UInt_t entries = 0; // bit coded, see GetDimParams() below 
+     entries = 1<<0 | 1<<1 | 1<<2 | 1<<3 | 1<<4 | 1<<5 | 1<<6 |1<<7; 
+     fhnDeltaR = NewTHnSparseF("fhnDeltaR", entries);
+
+     //change binning in pTtrack
+     Double_t *xPt3 = new Double_t[10];
+     xPt3[0] = 0.;
+     for(int i = 1; i<=9;i++){
+      if(xPt3[i-1]<2)xPt3[i] = xPt3[i-1] + 0.4; // 1 - 5
+      else if(xPt3[i-1]<11)xPt3[i] = xPt3[i-1] + 3; // 5 - 12
+      else xPt3[i] = xPt3[i-1] + 150.; // 18 
+     }
+    fhnDeltaR->SetBinEdges(2,xPt3);
+    delete [] xPt3;
+
+    //change binning in HTI
+     Double_t *xPt4 = new Double_t[14];
+     xPt4[0] = 0.;
+     for(int i = 1; i<=13;i++){
+      if(xPt4[i-1]<10)xPt4[i] = xPt4[i-1] + 1; // 1 - 10
+      else if(xPt4[i-1]<20)xPt4[i] = xPt4[i-1] + 5; // 10 - 12
+      else xPt4[i] = xPt4[i-1] + 30.; // 13
+     }
+    fhnDeltaR->SetBinEdges(6,xPt4);
+    delete [] xPt4;
 
+    
 
-   //UInt_t entries = 0; // bit coded, see GetDimParams() below
-   //UInt_t opt = 0;  // bit coded, default (0) or high resolution (1)
 
-   //  Int_t bins[5]={10,10,150,150,15};
-   //Double_t xmin[5]={0.,0.,0.,0.,0.};
-   //Double_t xmax[5]={100.,6.5,150.,1.5,1.5};     
-   //fhnDeltaRjets = new THnSparseF("fhnDeltaRjets", "fhnDeltaRjets",5,bins,xmin,xmax);
 
 
+   
+     if(fDoEventMixing){    
+     UInt_t cifras = 0; // bit coded, see GetDimParams() below 
+     cifras = 1<<0 | 1<<1 | 1<<2 | 1<<3 | 1<<4 | 1<<5 | 1<<7; 
+     fhnMixedEvents = NewTHnSparseF("fhnMixedEvents", cifras);}
 
+    if(fCheckMethods){
     fh2JetCoreMethod1C10 = new TH2F("JetCoreMethod1C10","",150, 0., 150.,100, 0., 1.5);
     fh2JetCoreMethod2C10 = new TH2F("JetCoreMethod2C10","",150, 0., 150.,100, 0., 1.5);
-    fh2JetCoreMethod3C10 = new TH2F("JetCoreMethod3C10","",150, 0., 150.,100, 0., 1.5); 
     fh2JetCoreMethod1C20 = new TH2F("JetCoreMethod1C20","",150, 0., 150.,100, 0., 1.5);
     fh2JetCoreMethod2C20 = new TH2F("JetCoreMethod2C20","",150, 0., 150.,100, 0., 1.5);
-    fh2JetCoreMethod3C20 = new TH2F("JetCoreMethod3C20","",150, 0., 150.,100, 0., 1.5); 
     fh2JetCoreMethod1C30 = new TH2F("JetCoreMethod1C30","",150, 0., 150.,100, 0., 1.5);
     fh2JetCoreMethod2C30 = new TH2F("JetCoreMethod2C30","",150, 0., 150.,100, 0., 1.5);
-    fh2JetCoreMethod3C30 = new TH2F("JetCoreMethod3C30","",150, 0., 150.,100, 0., 1.5); 
     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., 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.); 
-
-
-
-    fh2DeltaRC10pt1 = new TH2F("DeltaRC10pt1","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC20pt1 = new TH2F("DeltaRC20pt1","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC30pt1 = new TH2F("DeltaRC30pt1","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC60pt1 = new TH2F("DeltaRC60pt1","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC10pt2 = new TH2F("DeltaRC10pt2","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC20pt2 = new TH2F("DeltaRC20pt2","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC30pt2 = new TH2F("DeltaRC30pt2","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC60pt2 = new TH2F("DeltaRC60pt2","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC10pt3 = new TH2F("DeltaRC10pt3","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC20pt3 = new TH2F("DeltaRC20pt3","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC30pt3 = new TH2F("DeltaRC30pt3","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC60pt3 = new TH2F("DeltaRC60pt3","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC10pt4 = new TH2F("DeltaRC10pt4","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC20pt4 = new TH2F("DeltaRC20pt4","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC30pt4 = new TH2F("DeltaRC30pt4","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC60pt4 = new TH2F("DeltaRC60pt4","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaEtaC10pt1 = new TH2F("DeltaEtaC10pt1","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC20pt1 = new TH2F("DeltaEtaC20pt1","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC30pt1 = new TH2F("DeltaEtaC30pt1","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC60pt1 = new TH2F("DeltaEtaC60pt1","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC10pt2 = new TH2F("DeltaEtaC10pt2","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC20pt2 = new TH2F("DeltaEtaC20pt2","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC30pt2 = new TH2F("DeltaEtaC30pt2","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC60pt2 = new TH2F("DeltaEtaC60pt2","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC10pt3 = new TH2F("DeltaEtaC10pt3","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC20pt3 = new TH2F("DeltaEtaC20pt3","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC30pt3 = new TH2F("DeltaEtaC30pt3","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC60pt3 = new TH2F("DeltaEtaC60pt3","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC10pt4 = new TH2F("DeltaEtaC10pt4","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC20pt4 = new TH2F("DeltaEtaC20pt4","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC30pt4 = new TH2F("DeltaEtaC30pt4","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC60pt4 = new TH2F("DeltaEtaC60pt4","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaPhiC10pt1 = new TH2F("DeltaPhiC10pt1","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC20pt1 = new TH2F("DeltaPhiC20pt1","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC30pt1 = new TH2F("DeltaPhiC30pt1","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC60pt1 = new TH2F("DeltaPhiC60pt1","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC10pt2 = new TH2F("DeltaPhiC10pt2","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC20pt2 = new TH2F("DeltaPhiC20pt2","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC30pt2 = new TH2F("DeltaPhiC30pt2","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC60pt2 = new TH2F("DeltaPhiC60pt2","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC10pt3 = new TH2F("DeltaPhiC10pt3","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC20pt3 = new TH2F("DeltaPhiC20pt3","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC30pt3 = new TH2F("DeltaPhiC30pt3","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC60pt3 = new TH2F("DeltaPhiC60pt3","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC10pt4 = new TH2F("DeltaPhiC10pt4","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC20pt4 = new TH2F("DeltaPhiC20pt4","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC30pt4 = new TH2F("DeltaPhiC30pt4","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC60pt4 = new TH2F("DeltaPhiC60pt4","",150, 0., 15.,100,-6.5,6.5); 
-
-    fh2DeltaRC10pt1lead = new TH2F("DeltaRC10pt1lead","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC20pt1lead = new TH2F("DeltaRC20pt1lead","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC30pt1lead = new TH2F("DeltaRC30pt1lead","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC60pt1lead = new TH2F("DeltaRC60pt1lead","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC10pt2lead = new TH2F("DeltaRC10pt2lead","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC20pt2lead = new TH2F("DeltaRC20pt2lead","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC30pt2lead = new TH2F("DeltaRC30pt2lead","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC60pt2lead = new TH2F("DeltaRC60pt2lead","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC10pt3lead = new TH2F("DeltaRC10pt3lead","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC20pt3lead = new TH2F("DeltaRC20pt3lead","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC30pt3lead = new TH2F("DeltaRC30pt3lead","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC60pt3lead = new TH2F("DeltaRC60pt3lead","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC10pt4lead = new TH2F("DeltaRC10pt4lead","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC20pt4lead = new TH2F("DeltaRC20pt4lead","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC30pt4lead = new TH2F("DeltaRC30pt4lead","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC60pt4lead = new TH2F("DeltaRC60pt4lead","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaEtaC10pt1lead = new TH2F("DeltaEtaC10pt1lead","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC20pt1lead = new TH2F("DeltaEtaC20pt1lead","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC30pt1lead = new TH2F("DeltaEtaC30pt1lead","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC60pt1lead = new TH2F("DeltaEtaC60pt1lead","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC10pt2lead = new TH2F("DeltaEtaC10pt2lead","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC20pt2lead = new TH2F("DeltaEtaC20pt2lead","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC30pt2lead = new TH2F("DeltaEtaC30pt2lead","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC60pt2lead = new TH2F("DeltaEtaC60pt2lead","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC10pt3lead = new TH2F("DeltaEtaC10pt3lead","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC20pt3lead = new TH2F("DeltaEtaC20pt3lead","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC30pt3lead = new TH2F("DeltaEtaC30pt3lead","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC60pt3lead = new TH2F("DeltaEtaC60pt3lead","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC10pt4lead = new TH2F("DeltaEtaC10pt4lead","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC20pt4lead = new TH2F("DeltaEtaC20pt4lead","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC30pt4lead = new TH2F("DeltaEtaC30pt4lead","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC60pt4lead = new TH2F("DeltaEtaC60pt4lead","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaPhiC10pt1lead = new TH2F("DeltaPhiC10pt1lead","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC20pt1lead = new TH2F("DeltaPhiC20pt1lead","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC30pt1lead = new TH2F("DeltaPhiC30pt1lead","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC60pt1lead = new TH2F("DeltaPhiC60pt1lead","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC10pt2lead = new TH2F("DeltaPhiC10pt2lead","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC20pt2lead = new TH2F("DeltaPhiC20pt2lead","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC30pt2lead = new TH2F("DeltaPhiC30pt2lead","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC60pt2lead = new TH2F("DeltaPhiC60pt2lead","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC10pt3lead = new TH2F("DeltaPhiC10pt3lead","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC20pt3lead = new TH2F("DeltaPhiC20pt3lead","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC30pt3lead = new TH2F("DeltaPhiC30pt3lead","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC60pt3lead = new TH2F("DeltaPhiC60pt3lead","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC10pt4lead = new TH2F("DeltaPhiC10pt4lead","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC20pt4lead = new TH2F("DeltaPhiC20pt4lead","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC30pt4lead = new TH2F("DeltaPhiC30pt4lead","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC60pt4lead = new TH2F("DeltaPhiC60pt4lead","",150, 0., 15.,100,-6.5,6.5); 
-
-
-    fh2DeltaRC10pt1sublead = new TH2F("DeltaRC10pt1sublead","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC20pt1sublead = new TH2F("DeltaRC20pt1sublead","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC30pt1sublead = new TH2F("DeltaRC30pt1sublead","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC60pt1sublead = new TH2F("DeltaRC60pt1sublead","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC10pt2sublead = new TH2F("DeltaRC10pt2sublead","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC20pt2sublead = new TH2F("DeltaRC20pt2sublead","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC30pt2sublead = new TH2F("DeltaRC30pt2sublead","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC60pt2sublead = new TH2F("DeltaRC60pt2sublead","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC10pt3sublead = new TH2F("DeltaRC10pt3sublead","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC20pt3sublead = new TH2F("DeltaRC20pt3sublead","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC30pt3sublead = new TH2F("DeltaRC30pt3sublead","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC60pt3sublead = new TH2F("DeltaRC60pt3sublead","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC10pt4sublead = new TH2F("DeltaRC10pt4sublead","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC20pt4sublead = new TH2F("DeltaRC20pt4sublead","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC30pt4sublead = new TH2F("DeltaRC30pt4sublead","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaRC60pt4sublead = new TH2F("DeltaRC60pt4sublead","",150, 0., 15.,100,0.,1.5); 
-    fh2DeltaEtaC10pt1sublead = new TH2F("DeltaEtaC10pt1sublead","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC20pt1sublead = new TH2F("DeltaEtaC20pt1sublead","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC30pt1sublead = new TH2F("DeltaEtaC30pt1sublead","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC60pt1sublead = new TH2F("DeltaEtaC60pt1sublead","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC10pt2sublead = new TH2F("DeltaEtaC10pt2sublead","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC20pt2sublead = new TH2F("DeltaEtaC20pt2sublead","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC30pt2sublead = new TH2F("DeltaEtaC30pt2sublead","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC60pt2sublead = new TH2F("DeltaEtaC60pt2sublead","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC10pt3sublead = new TH2F("DeltaEtaC10pt3sublead","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC20pt3sublead = new TH2F("DeltaEtaC20pt3sublead","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC30pt3sublead = new TH2F("DeltaEtaC30pt3sublead","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC60pt3sublead = new TH2F("DeltaEtaC60pt3sublead","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC10pt4sublead = new TH2F("DeltaEtaC10pt4sublead","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC20pt4sublead = new TH2F("DeltaEtaC20pt4sublead","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC30pt4sublead = new TH2F("DeltaEtaC30pt4sublead","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaEtaC60pt4sublead = new TH2F("DeltaEtaC60pt4sublead","",150, 0., 15.,100,-1.5,1.5); 
-    fh2DeltaPhiC10pt1sublead = new TH2F("DeltaPhiC10pt1sublead","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC20pt1sublead = new TH2F("DeltaPhiC20pt1sublead","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC30pt1sublead = new TH2F("DeltaPhiC30pt1sublead","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC60pt1sublead = new TH2F("DeltaPhiC60pt1sublead","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC10pt2sublead = new TH2F("DeltaPhiC10pt2sublead","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC20pt2sublead = new TH2F("DeltaPhiC20pt2sublead","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC30pt2sublead = new TH2F("DeltaPhiC30pt2sublead","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC60pt2sublead = new TH2F("DeltaPhiC60pt2sublead","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC10pt3sublead = new TH2F("DeltaPhiC10pt3sublead","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC20pt3sublead = new TH2F("DeltaPhiC20pt3sublead","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC30pt3sublead = new TH2F("DeltaPhiC30pt3sublead","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC60pt3sublead = new TH2F("DeltaPhiC60pt3sublead","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC10pt4sublead = new TH2F("DeltaPhiC10pt4sublead","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC20pt4sublead = new TH2F("DeltaPhiC20pt4sublead","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC30pt4sublead = new TH2F("DeltaPhiC30pt4sublead","",150, 0., 15.,100,-6.5,6.5); 
-    fh2DeltaPhiC60pt4sublead = new TH2F("DeltaPhiC60pt4sublead","",150, 0., 15.,100,-6.5,6.5); 
-
-
-
-
-
-      if(fAngStructCloseTracks>0){
+    fh2JetCoreMethod2C60 = new TH2F("JetCoreMethod2C60","",150, 0., 150.,100, 0., 1.5);}
+
+    fh3JetTrackC3060=new TH3F("JetTrackC3060","",50,0,50,150,0.,150.,35,0.,3.5);
+    fh3JetTrackC20=new TH3F("JetTrackC20","",50,0,50,150,0.,150.,35,0.,3.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.);
@@ -928,235 +433,46 @@ void AliAnalysisTaskJetCore::UserCreateOutputObjects()
     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.); }
 
+   
+    fh1TrigRef=new TH1D("Trig Ref","",10,0.,10);
+    fh1TrigSig=new TH1D("Trig Sig","",10,0.,10);  
+    fh2Ntriggers=new TH2F("# of triggers","",100,0.,100.,50,0.,50.);
+    fh2Ntriggers2C10=new TH2F("# of triggers2C10","",50,0.,50.,50,0.,50.);
+    fh2Ntriggers2C20=new TH2F("# of triggers2C20","",50,0.,50.,50,0.,50.);
+    fh3JetDensity=new TH3F("Jet density vs mutliplicity A>0.07","",100,0.,4000.,100,0.,5.,25,0.,50.);
+    fh3JetDensityA4=new TH3F("Jet density vs multiplicity A>0.4","",100,0.,4000.,100,0.,5.,25,0.,50.);
+    fh2RPJetsC10=new TH2F("RPJetC10","",35,0.,3.5,100,0.,100.);
+    fh2RPJetsC20=new TH2F("RPJetC20","",35,0.,3.5,100,0.,100.); 
+    fh2RPTC10=new TH2F("RPTriggerC10","",35,0.,3.5,50,0.,50.); 
+    fh2RPTC20=new TH2F("RPTriggerC20","",35,0.,3.5,50,0.,50.);  
+
+
+    
+    
    fOutputList->Add(fHistEvtSelection);
-   fOutputList->Add(fHistJetSelection);
-   fOutputList->Add(fh2JetSelection);
+
+   fOutputList->Add(fhnDeltaR);
+   
+   fOutputList->Add(fhnMixedEvents);
+         
+     
+  
+      if(fCheckMethods){
 
       fOutputList->Add(fh2JetCoreMethod1C10);
       fOutputList->Add(fh2JetCoreMethod2C10);
-      fOutputList->Add(fh2JetCoreMethod3C10);
       fOutputList->Add(fh2JetCoreMethod1C20);
       fOutputList->Add(fh2JetCoreMethod2C20);
-      fOutputList->Add(fh2JetCoreMethod3C20);
       fOutputList->Add(fh2JetCoreMethod1C30);
       fOutputList->Add(fh2JetCoreMethod2C30);
-      fOutputList->Add(fh2JetCoreMethod3C30);
       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(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(fh2JetCoreMethod2C60);}
       
-      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);
-      fOutputList->Add(fh2DeltaRC60pt1); 
-      fOutputList->Add(fh2DeltaRC10pt2);
-      fOutputList->Add(fh2DeltaRC20pt2);
-      fOutputList->Add(fh2DeltaRC30pt2);
-      fOutputList->Add(fh2DeltaRC60pt2); 
-      fOutputList->Add(fh2DeltaRC10pt3);
-      fOutputList->Add(fh2DeltaRC20pt3);
-      fOutputList->Add(fh2DeltaRC30pt3);
-      fOutputList->Add(fh2DeltaRC60pt3); 
-      fOutputList->Add(fh2DeltaRC10pt4);
-      fOutputList->Add(fh2DeltaRC20pt4);
-      fOutputList->Add(fh2DeltaRC30pt4);
-      fOutputList->Add(fh2DeltaRC60pt4);
-      fOutputList->Add(fh2DeltaEtaC10pt1);
-      fOutputList->Add(fh2DeltaEtaC20pt1);
-      fOutputList->Add(fh2DeltaEtaC30pt1);
-      fOutputList->Add(fh2DeltaEtaC60pt1); 
-      fOutputList->Add(fh2DeltaEtaC10pt2);
-      fOutputList->Add(fh2DeltaEtaC20pt2);
-      fOutputList->Add(fh2DeltaEtaC30pt2);
-      fOutputList->Add(fh2DeltaEtaC60pt2); 
-      fOutputList->Add(fh2DeltaEtaC10pt3);
-      fOutputList->Add(fh2DeltaEtaC20pt3);
-      fOutputList->Add(fh2DeltaEtaC30pt3);
-      fOutputList->Add(fh2DeltaEtaC60pt3); 
-      fOutputList->Add(fh2DeltaEtaC10pt4);
-      fOutputList->Add(fh2DeltaEtaC20pt4);
-      fOutputList->Add(fh2DeltaEtaC30pt4);
-      fOutputList->Add(fh2DeltaEtaC60pt4); 
-      fOutputList->Add(fh2DeltaPhiC10pt1);
-      fOutputList->Add(fh2DeltaPhiC20pt1);
-      fOutputList->Add(fh2DeltaPhiC30pt1);
-      fOutputList->Add(fh2DeltaPhiC60pt1); 
-      fOutputList->Add(fh2DeltaPhiC10pt2);
-      fOutputList->Add(fh2DeltaPhiC20pt2);
-      fOutputList->Add(fh2DeltaPhiC30pt2);
-      fOutputList->Add(fh2DeltaPhiC60pt2); 
-      fOutputList->Add(fh2DeltaPhiC10pt3);
-      fOutputList->Add(fh2DeltaPhiC20pt3);
-      fOutputList->Add(fh2DeltaPhiC30pt3);
-      fOutputList->Add(fh2DeltaPhiC60pt3); 
-      fOutputList->Add(fh2DeltaPhiC10pt4);
-      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);      
+      fOutputList->Add(fh3JetTrackC3060);
+      fOutputList->Add(fh3JetTrackC20);
+     
+            
      
-       
-
 
 
         if(fAngStructCloseTracks>0){
@@ -1180,6 +496,59 @@ void AliAnalysisTaskJetCore::UserCreateOutputObjects()
 
 
 
+        fOutputList->Add(fh1TrigRef);
+        fOutputList->Add(fh1TrigSig); 
+       fOutputList->Add(fh2Ntriggers);
+        fOutputList->Add(fh2Ntriggers2C10);
+        fOutputList->Add(fh2Ntriggers2C20); 
+        fOutputList->Add(fh3JetDensity);
+        fOutputList->Add(fh3JetDensityA4);
+        fOutputList->Add(fh2RPJetsC10);
+        fOutputList->Add(fh2RPJetsC20);
+         fOutputList->Add(fh2RPTC10);
+        fOutputList->Add(fh2RPTC20);
+
+        const Int_t dimSpec = 5;
+       const Int_t nBinsSpec[dimSpec]     = {100,6, 140, 50, fNRPBins};
+       const Double_t lowBinSpec[dimSpec] = {0,0,-80, 0, 0};
+       const Double_t hiBinSpec[dimSpec]  = {100,1, 200, 50,  static_cast<Double_t>(fNRPBins)};
+       fHJetSpec = new THnSparseF("fHJetSpec","Recoil jet spectrum",dimSpec,nBinsSpec,lowBinSpec,hiBinSpec);
+
+             //change binning in jet area
+     Double_t *xPt6 = new Double_t[7];
+     xPt6[0] = 0.;
+     xPt6[1]=0.07;
+     xPt6[2]=0.2;
+     xPt6[3]=0.4;
+     xPt6[4]=0.6;
+     xPt6[5]=0.8; 
+     xPt6[6]=1;
+    fHJetSpec->SetBinEdges(1,xPt6);
+    delete [] xPt6;
+
+
+
+
+
+       fOutputList->Add(fHJetSpec);  
+
+
+       if(fRunAnaAzimuthalCorrelation)
+         {
+           fhTTPt = new TH2F("fhTTPt","Trigger track p_{T} vs centrality",10,0,100,100,0,100);
+           fOutputList->Add(fhTTPt);
+
+           const Int_t dimCor = 5;
+           const Int_t nBinsCor[dimCor]     = {50, 200, 100,              8,   10};
+           const Double_t lowBinCor[dimCor] = {0,  -50, -0.5*TMath::Pi(), 0,   0};
+           const Double_t hiBinCor[dimCor]  = {50, 150, 1.5*TMath::Pi(),  0.8, 100};
+           fHJetPhiCorr = new THnSparseF("fHJetPhiCorr","TT p_{T} vs jet p_{T} vs dPhi vs area vs centrality",dimCor,nBinsCor,lowBinCor,hiBinCor);
+           fOutputList->Add(fHJetPhiCorr);
+         }
+
+
+        
+     
    // =========== Switch on Sumw2 for all histos ===========
    for (Int_t i=0; i<fOutputList->GetEntries(); ++i) {
       TH1 *h1 = dynamic_cast<TH1*>(fOutputList->At(i));
@@ -1209,10 +578,17 @@ void AliAnalysisTaskJetCore::UserExec(Option_t *)
    fESD=dynamic_cast<AliESDEvent*>(InputEvent());
    if (!fESD) {
       AliError("ESD not available");
-      fAOD = dynamic_cast<AliAODEvent*>(InputEvent());
-   } else {
-      fAOD = dynamic_cast<AliAODEvent*>(AODEvent());
-   }
+      fAODIn = dynamic_cast<AliAODEvent*>(InputEvent());
+   } 
+      fAODOut = dynamic_cast<AliAODEvent*>(AODEvent());
+
+       static AliAODEvent* aod = 0;
+       // take all other information from the aod we take the tracks from
+       if(!aod){
+       if(!fESD)aod = fAODIn;
+       else aod = fAODOut;}
+
+   
  
     if(fNonStdFile.Length()!=0){
     // case that we have an AOD extension we need can fetch the jets from the extended output
@@ -1223,15 +599,23 @@ void AliAnalysisTaskJetCore::UserExec(Option_t *)
     }}
     
 
+   // -- event selection --
+   fHistEvtSelection->Fill(1); // number of events before event selection
 
 
+       Bool_t selected=kTRUE;
+       selected = AliAnalysisHelperJetTasks::Selected();
+       if(!selected){
+      // no selection by the service task, we continue
+       PostData(1,fOutputList);
+       return;}
+    
 
-   // -- event selection --
-   fHistEvtSelection->Fill(1); // number of events before event selection
 
-   // physics selection
+  // physics selection: this is now redundant, all should appear as accepted after service task selection
    AliInputEventHandler* inputHandler = (AliInputEventHandler*)
    ((AliAnalysisManager::GetAnalysisManager())->GetInputEventHandler());
+        std::cout<<inputHandler->IsEventSelected()<<" "<<fOfflineTrgMask<<std::endl;
    if(!(inputHandler->IsEventSelected() & fOfflineTrgMask)){
       if(fDebug) Printf(" Trigger Selection: event REJECTED ... ");
       fHistEvtSelection->Fill(2);
@@ -1239,8 +623,23 @@ void AliAnalysisTaskJetCore::UserExec(Option_t *)
       return;
    }
 
+
+      
    // vertex selection
-   AliAODVertex* primVtx = fAOD->GetPrimaryVertex();
+   if(!aod){
+     if(fDebug) Printf("%s:%d No AOD",(char*)__FILE__,__LINE__);
+     fHistEvtSelection->Fill(3);
+      PostData(1, fOutputList);
+   }
+   AliAODVertex* primVtx = aod->GetPrimaryVertex();
+
+   if(!primVtx){
+     if(fDebug) Printf("%s:%d No primVtx",(char*)__FILE__,__LINE__);
+     fHistEvtSelection->Fill(3);
+     PostData(1, fOutputList);
+     return;
+   }
+
    Int_t nTracksPrim = primVtx->GetNContributors();
    if ((nTracksPrim < fMinContribVtx) ||
          (primVtx->GetZ() < fVtxZMin) ||
@@ -1262,76 +661,95 @@ void AliAnalysisTaskJetCore::UserExec(Option_t *)
 
    // centrality selection
    AliCentrality *cent = 0x0;
-   Float_t centValue = 0.; 
-   if(fESD) cent = fESD->GetCentrality();
-   if(cent) centValue = cent->GetCentralityPercentile("V0M");
+   Double_t centValue = 0.; 
+   if(fIsPbPb){
+   if(fESD) {cent = fESD->GetCentrality();
+     if(cent) centValue = cent->GetCentralityPercentile("V0M");}
+   else     centValue=((AliVAODHeader*)aod->GetHeader())->GetCentrality();
+   
    if(fDebug) printf("centrality: %f\n", centValue);
-   if (centValue < fCentMin || centValue > fCentMax){
+      if (centValue < fCentMin || centValue > fCentMax){
       fHistEvtSelection->Fill(4);
       PostData(1, fOutputList);
       return;
-   }
+      }}
 
 
-   // multiplicity due to input tracks
-   //Int_t nInputTracks = GetNInputTracks();
-   //if (nInputTracks < fNInputTracksMin || (fNInputTracksMax > -1 && nInputTracks > fNInputTracksMax)){
-   //   fHistEvtSelection->Fill(5);
-   //   PostData(1, fOutputList);
-   //   return;
-   // }
-
-   
    fHistEvtSelection->Fill(0); 
    // accepted events  
    // -- end event selection --
   
    // get background
    AliAODJetEventBackground* externalBackground = 0;
-   if(fAOD&&!externalBackground&&fBackgroundBranch.Length()){
-      externalBackground =  (AliAODJetEventBackground*)(fAOD->FindListObject(fBackgroundBranch.Data()));
+   if(fAODOut&&!externalBackground&&fBackgroundBranch.Length()){
+      externalBackground =  (AliAODJetEventBackground*)(fAODOut->FindListObject(fBackgroundBranch.Data()));
       if(!externalBackground)Printf("%s:%d Background branch not found %s",(char*)__FILE__,__LINE__,fBackgroundBranch.Data());;
    }
    if(fAODExtension&&!externalBackground&&fBackgroundBranch.Length()){
      externalBackground =  (AliAODJetEventBackground*)(fAODExtension->GetAOD()->FindListObject(fBackgroundBranch.Data()));
       if(!externalBackground)Printf("%s:%d Background branch not found %s",(char*)__FILE__,__LINE__,fBackgroundBranch.Data());;
    }
+
+    if(fAODIn&&!externalBackground&&fBackgroundBranch.Length()){
+      externalBackground =  (AliAODJetEventBackground*)(fAODIn->FindListObject(fBackgroundBranch.Data()));
+      if(!externalBackground)Printf("%s:%d Background branch not found %s",(char*)__FILE__,__LINE__,fBackgroundBranch.Data());;
+    }
    
    Float_t rho = 0;
-   if(externalBackground)rho = externalBackground->GetBackground(0);
-
 
+   if(fFlagRandom==0){
+     if(externalBackground)rho = externalBackground->GetBackground(0);}
+   if(fFlagRandom==2){
+      if(externalBackground)rho = externalBackground->GetBackground(2);}
+   if(fFlagRandom==3){
+      if(externalBackground)rho = externalBackground->GetBackground(3);}
    // fetch jets
    TClonesArray *aodJets[2];
    aodJets[0]=0;
-   if(fAOD&&!aodJets[0]){
-   aodJets[0] = dynamic_cast<TClonesArray*>(fAOD->FindListObject(fJetBranchName[0].Data())); 
-   aodJets[1] = dynamic_cast<TClonesArray*>(fAOD->FindListObject(fJetBranchName[1].Data()));  }
+   if(fAODOut&&!aodJets[0]){
+   aodJets[0] = dynamic_cast<TClonesArray*>(fAODOut->FindListObject(fJetBranchName[0].Data())); 
+   aodJets[1] = dynamic_cast<TClonesArray*>(fAODOut->FindListObject(fJetBranchName[1].Data()));  }
    if(fAODExtension && !aodJets[0]){ 
-     aodJets[0] = dynamic_cast<TClonesArray*>(fAODExtension->GetAOD()->FindListObject(fJetBranchName[0].Data())); 
-     aodJets[1] = dynamic_cast<TClonesArray*>(fAODExtension->GetAOD()->FindListObject(fJetBranchName[1].Data()));  }
+   aodJets[0] = dynamic_cast<TClonesArray*>(fAODExtension->GetAOD()->FindListObject(fJetBranchName[0].Data())); 
+   aodJets[1] = dynamic_cast<TClonesArray*>(fAODExtension->GetAOD()->FindListObject(fJetBranchName[1].Data()));  }
+     if(fAODIn&&!aodJets[0]){
+   aodJets[0] = dynamic_cast<TClonesArray*>(fAODIn->FindListObject(fJetBranchName[0].Data())); 
+   aodJets[1] = dynamic_cast<TClonesArray*>(fAODIn->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;}   
 
+   Int_t nT=0;
    TList ParticleList;
-   Int_t nT = GetListOfTracks(&ParticleList);
+   Double_t minT=0;
+   Double_t maxT=0;
+   Int_t number=0;
+   Double_t dice=fRandom->Uniform(0,1);
+   if(dice>fFrac){ minT=fTTLowRef;
+                   maxT=fTTUpRef;}
+   if(dice<=fFrac){minT=fTTLowSig;
+                   maxT=fTTUpSig;} 
+
+
+
+   if(fHardest==1 || fHardest==2) nT = GetListOfTracks(&ParticleList);
+   if(fHardest==0) nT=SelectTrigger(&ParticleList,minT,maxT,number);
+   if(nT<0){  
+   PostData(1, fOutputList);
+   return;}   
+
+      if(dice>fFrac) fh1TrigRef->Fill(number);
+      if(dice<=fFrac)fh1TrigSig->Fill(number);
+
      for (Int_t iJetType = 0; iJetType < 2; iJetType++) {
       fListJets[iJetType]->Clear();
       if (!aodJets[iJetType]) continue;
-
       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);
-         if(iJetType==0){
-          ptsub[iJet]=jet->Pt()-rho*jet->EffectiveAreaCharged();
-        }}}
+        // if(iJetType==0){
+        // ptsub[iJet]=jet->Pt()-rho*jet->EffectiveAreaCharged();}
+      }}
    
    Double_t etabig=0;
    Double_t ptbig=0;
@@ -1339,39 +757,53 @@ void AliAnalysisTaskJetCore::UserExec(Option_t *)
    Double_t phibig=0.;
    Double_t etasmall=0;
    Double_t ptsmall=0;
-   Double_t areasmall=0;
-   Double_t distr=0.;
+   //   Double_t areasmall=0;
    Double_t phismall=0.;
-   Int_t indexlead=-1;
-   Int_t indexsublead=-1;
-   Int_t indexstop=-1;
-  
-   if(fListJets[0]->GetEntries()>0) 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;}}
          
+  
+
+   Int_t iCount=0; 
+   Int_t trigJet=-1;
+   Int_t trigBBTrack=-1;
+          // Int_t trigInTrack=-1;
+   fRPAngle = ((AliVAODHeader*)aod->GetHeader())->GetEventplane();     
+
+   if(fHardest==0 || fHardest==1){
+   AliVParticle *partback = (AliVParticle*)ParticleList.At(nT);     
+   if(!partback){  
+   PostData(1, fOutputList);
+   return;}
+
+    if(fSemigoodCorrect){
+   Double_t disthole=RelativePhi(partback->Phi(),fHolePos);
+   if(TMath::Abs(disthole)+fHoleWidth>TMath::Pi()-0.6){ 
+   PostData(1, fOutputList);
+   return;}} 
+
+
+    }
+
+
+   for(Int_t tt=0;tt<ParticleList.GetEntries();tt++){
+     if(fHardest==0||fHardest==1){if(tt!=nT) continue;}
+     AliVParticle *partback = (AliVParticle*)ParticleList.At(tt);     
+     if(!partback) continue;
+     if(partback->Pt()<8) continue;
 
-   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.};
-   Double_t up3[15]={0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
-   Double_t up4[15]={0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
-   Double_t down1[15]={0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
-   Double_t down2[15]={0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
-   Double_t down3[15]={0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
-   Double_t down4[15]={0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
+   Double_t accep=2.*TMath::Pi()*1.8;
+   Int_t injet4=0;
+   Int_t injet=0; 
+
+   if(fSemigoodCorrect){
+   Double_t disthole=RelativePhi(partback->Phi(),fHolePos);
+   if(TMath::Abs(disthole)+fHoleWidth>TMath::Pi()-0.6) continue;}
+
+   
+
+   fh2Ntriggers->Fill(centValue,partback->Pt());
+   Double_t phiBinT = RelativePhi(partback->Phi(),fRPAngle);
+   if(centValue<20.) fh2RPTC20->Fill(TMath::Abs(phiBinT),partback->Pt());
+   if(centValue<10.) fh2RPTC10->Fill(TMath::Abs(phiBinT),partback->Pt());
 
 
 
@@ -1381,57 +813,100 @@ void AliAnalysisTaskJetCore::UserExec(Option_t *)
            phibig  = jetbig->Phi();
            ptbig   = jetbig->Pt();
            if(ptbig==0) continue; 
+           Double_t phiBin = RelativePhi(phibig,fRPAngle);       
            areabig = jetbig->EffectiveAreaCharged();
            Double_t ptcorr=ptbig-rho*areabig;
-           if(ptcorr<=0) continue;
           if((etabig<fJetEtaMin)||(etabig>fJetEtaMax)) continue;
+           if(areabig>=0.07) injet=injet+1;
+           if(areabig>=0.4) injet4=injet4+1;   
+           Double_t dphi=RelativePhi(partback->Phi(),phibig); 
+
+           if(fFlagEtaBkg==1){
+          Double_t etadif= partback->Eta()-etabig;
+           if(TMath::Abs(etadif)<=0.5){             
+          
+           if(centValue<20.) fh3JetTrackC20->Fill(partback->Pt(),ptcorr,TMath::Abs(dphi));
+           if(centValue>30. && centValue<60.) fh3JetTrackC3060->Fill(partback->Pt(),ptcorr,TMath::Abs(dphi));}}
+           if(fFlagEtaBkg==0){
+           if(centValue<20.) fh3JetTrackC20->Fill(partback->Pt(),ptcorr,TMath::Abs(dphi));
+           if(centValue>30. && centValue<60.) fh3JetTrackC3060->Fill(partback->Pt(),ptcorr,TMath::Abs(dphi));}
+
+
+           if(fFlagJetHadron==0){
+           if(fFlagPhiBkg==1) if((TMath::Abs(dphi)<TMath::Pi()/2.-0.1)||(TMath::Abs(dphi)>TMath::Pi()/2.+0.1)) continue;
+           if(fFlagPhiBkg==0) if(TMath::Abs(dphi)<TMath::Pi()-0.6) continue;
+           if(fFlagPhiBkg==2) if(TMath::Abs(dphi)<TMath::Pi()-0.7) continue;
+           if(fFlagPhiBkg==3) if(TMath::Abs(dphi)<TMath::Pi()-0.5) continue;}
+ 
+           if(fFlagJetHadron!=0) if(TMath::Abs(dphi)>0.4) continue;
+
+
+          if(centValue<10.) fh2RPJetsC10->Fill(TMath::Abs(phiBin), ptcorr);
+          if(centValue<20.) fh2RPJetsC20->Fill(TMath::Abs(phiBin), ptcorr);
                    Double_t dismin=100.;
                    Double_t ptmax=-10.; 
                    Int_t index1=-1;
                    Int_t index2=-1;
-                   Double_t fracin=0.;
-                   Double_t sumPtIn=0.;
-                   Double_t sumPtOut=0.; 
-          //compute sum of the pt of the tracks in a concentric cone
-           TRefArray *genTrackList = jetbig->GetRefTracks();
-           Int_t nTracksGenJet = genTrackList->GetEntriesFast();
-           AliAODTrack* genTrack;
-             for(Int_t ir=0; ir<nTracksGenJet; ++ir){
-             genTrack = (AliAODTrack*)(genTrackList->At(ir));
-            Float_t etr=genTrack->Eta();
-             Float_t phir=genTrack->Phi();
-             distr=(etr-etabig)*(etr-etabig)+(phir-phibig)*(phir-phibig);
-             distr=TMath::Sqrt(distr);
-            if(distr<=fRadioFrac){ fracin=fracin+genTrack->Pt();}}
-    
-            if(centValue<10) {fh2JetCoreMethod3C10->Fill(ptcorr,fracin/ptbig);
-              if(i==indexlead) fh2JetCoreMethod3C10lead->Fill(ptcorr,fracin/ptbig);
-               if(i==indexsublead) fh2JetCoreMethod3C10sublead->Fill(ptcorr,fracin/ptbig);}
+         
+                  Float_t phitt=partback->Phi();
+                   if(phitt<0)phitt+=TMath::Pi()*2.; 
+                   Int_t phiBintt = GetPhiBin(phitt-fRPAngle);
 
-             
-             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);}
+                  Double_t fillspec[] = {centValue,jetbig->EffectiveAreaCharged(),ptcorr,partback->Pt(), static_cast<Double_t>(phiBintt)};
+                 fHJetSpec->Fill(fillspec);
+           
+          
+
+                   if(ptcorr<=0) continue;
+
+                       AliAODTrack* leadtrack=0; 
+                       Int_t ippt=0;
+                       Double_t ppt=-10;
+                       if(fFlagJetHadron==0){   
+                        TRefArray *genTrackList = jetbig->GetRefTracks();
+                        Int_t nTracksGenJet = genTrackList->GetEntriesFast();
+                        AliAODTrack* genTrack;
+                        for(Int_t ir=0; ir<nTracksGenJet; ++ir){
+                          genTrack = (AliAODTrack*)(genTrackList->At(ir));
+                          if(genTrack->Pt()>ppt){ppt=genTrack->Pt();
+                            ippt=ir;}}
+                        leadtrack=(AliAODTrack*)(genTrackList->At(ippt));
+                        if(!leadtrack) continue;
+                      }
+
+
+
+                      AliVParticle* leadtrackb=0;
+                       if(fFlagJetHadron!=0){
+                        Int_t nTb = GetHardestTrackBackToJet(jetbig);
+                         leadtrackb = (AliVParticle*)ParticleList.At(nTb);
+                         if(!leadtrackb) continue;  
+                      }
+
+
+
+
+                       
+                      //store one trigger info                   
+                      if(iCount==0){                        
+                        trigJet=i;
+                        trigBBTrack=nT;
+                        //                      trigInTrack=ippt;
+                        iCount=iCount+1;} 
+                      
    
-            
-
-                for(Int_t j=0; j<fListJets[1]->GetEntries(); ++j){
-                  AliAODJet* jetsmall = (AliAODJet*)(fListJets[1]->At(j));
-                  etasmall  = jetsmall->Eta();
-                  phismall = jetsmall->Phi();
-                  ptsmall   = jetsmall->Pt();
-                  areasmall = jetsmall->EffectiveAreaCharged();
-                  Double_t tmpDeltaR=(phismall-phibig)*(phismall-phibig)+(etasmall-etabig)*(etasmall-etabig);
-                 tmpDeltaR=TMath::Sqrt(tmpDeltaR);
-                    //Fraction in the jet core  
-                    if((ptsmall>ptmax)&&(tmpDeltaR<=fRadioFrac)){ptmax=ptsmall;  
-                   index2=j;}  
+                      if(fCheckMethods){
+                        for(Int_t j=0; j<fListJets[1]->GetEntries(); ++j){
+                          AliAODJet* jetsmall = (AliAODJet*)(fListJets[1]->At(j));
+                          etasmall  = jetsmall->Eta();
+                          phismall = jetsmall->Phi();
+                          ptsmall   = jetsmall->Pt();
+                          //                      areasmall = jetsmall->EffectiveAreaCharged();
+                          Double_t tmpDeltaR=(phismall-phibig)*(phismall-phibig)+(etasmall-etabig)*(etasmall-etabig);
+                          tmpDeltaR=TMath::Sqrt(tmpDeltaR);
+                          //Fraction in the jet core  
+                          if((ptsmall>ptmax)&&(tmpDeltaR<=fRadioFrac)){ptmax=ptsmall;  
+                            index2=j;}  
                     if(tmpDeltaR<=dismin){ dismin=tmpDeltaR;
                      index1=j;}} //en of loop over R=0.2 jets
                 //method1:most concentric jet=core 
@@ -1446,390 +921,372 @@ void AliAnalysisTaskJetCore::UserExec(Option_t *)
                   if(centValue<10) fh2JetCoreMethod2C10->Fill(ptcorr,jetmethod2->Pt()/ptbig);
                   if((centValue>20)&&(centValue<40)) fh2JetCoreMethod2C20->Fill(ptcorr,jetmethod2->Pt()/ptbig); 
                  if((centValue>30)&&(centValue<60)) fh2JetCoreMethod2C30->Fill(ptcorr,jetmethod2->Pt()/ptbig);
-                 if(centValue>60) fh2JetCoreMethod2C60->Fill(ptcorr,jetmethod2->Pt()/ptbig); }  
-
-     Double_t R=fRadioFrac*2.;      
-  
-     for(int it = 0;it<nT;++it){
-
+                 if(centValue>60) fh2JetCoreMethod2C60->Fill(ptcorr,jetmethod2->Pt()/ptbig); }}  
+                 if(centValue<10&&leadtrack) fh2Ntriggers2C10->Fill(leadtrack->Pt(),partback->Pt());                  
+                  if(centValue<20&&leadtrack) fh2Ntriggers2C20->Fill(leadtrack->Pt(),partback->Pt());  
+         if(fDoEventMixing==0 && fFlagOnlyRecoil==0){ 
+        for(int it = 0;it<ParticleList.GetEntries();++it){
          AliVParticle *part = (AliVParticle*)ParticleList.At(it);
-         Float_t deltaR = jetbig->DeltaR(part);
-          Float_t deltaEta = part->Eta()-etabig;
-          Float_t deltaPhi = part->Phi()-phibig;  
-         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);}
-              if((ptcorr>=85.)&&(ptcorr<100.)) {fh2DeltaRC10pt2->Fill(part->Pt(),deltaR);
-               if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC10pt2->Fill(part->Pt(),deltaEta);
-               if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC10pt2->Fill(part->Pt(),deltaPhi); }
-              if((ptcorr>=100.)&&(ptcorr<120.)) {fh2DeltaRC10pt3->Fill(part->Pt(),deltaR);
-              if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC10pt3->Fill(part->Pt(),deltaEta);
-               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(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); }} 
-
- 
-
-
-         }
+                 Double_t deltaR = jetbig->DeltaR(part);
+          Double_t deltaEta = etabig-part->Eta();
+         
+          Double_t deltaPhi=phibig-part->Phi();
+          if(deltaPhi<-0.5*TMath::Pi()) deltaPhi+=2.*TMath::Pi();
+          if(deltaPhi>3./2.*TMath::Pi()) deltaPhi-=2.*TMath::Pi();
+         Double_t pTcont=0;
+          if(fFlagJetHadron==0) pTcont=leadtrack->Pt();
+          if(fFlagJetHadron!=0) pTcont=leadtrackb->Pt(); 
+          Double_t jetEntries[8] = {centValue,ptcorr,part->Pt(),deltaR,deltaEta,deltaPhi,pTcont,partback->Pt()};  
+           fhnDeltaR->Fill(jetEntries);}
 
-          if((centValue>20.)&&(centValue<40.)){
-              if((ptcorr>=70.)&&(ptcorr<85.)) {fh2DeltaRC20pt1->Fill(part->Pt(),deltaR);
-             if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R))fh2DeltaEtaC20pt1->Fill(part->Pt(),deltaEta);
-             if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC20pt1->Fill(part->Pt(),deltaPhi);}
-                     if((ptcorr>=85.)&&(ptcorr<100.)) {fh2DeltaRC20pt2->Fill(part->Pt(),deltaR);
-              if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC20pt2->Fill(part->Pt(),deltaEta);
-              if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC20pt2->Fill(part->Pt(),deltaPhi); }
-             if((ptcorr>=100.)&&(ptcorr<120.)) {fh2DeltaRC20pt3->Fill(part->Pt(),deltaR);
-             if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC20pt3->Fill(part->Pt(),deltaEta);
-              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(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); }} 
 
          }
+        
+        
+          //end of track loop, we only do it if EM is switched off
+         
 
-         if((centValue>30.)&&(centValue<60.)){
-
-              if((ptcorr>=70.)&&(ptcorr<85.)) {fh2DeltaRC30pt1->Fill(part->Pt(),deltaR);
-             if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R))fh2DeltaEtaC30pt1->Fill(part->Pt(),deltaEta);
-             if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC30pt1->Fill(part->Pt(),deltaPhi);}
-             if((ptcorr>=85.)&&(ptcorr<100.)) {fh2DeltaRC30pt2->Fill(part->Pt(),deltaR);
-              if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC30pt2->Fill(part->Pt(),deltaEta);
-              if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC30pt2->Fill(part->Pt(),deltaPhi); }
-             if((ptcorr>=100.)&&(ptcorr<120.)) {fh2DeltaRC30pt3->Fill(part->Pt(),deltaR);
-             if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC30pt3->Fill(part->Pt(),deltaEta);
-              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(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.){
-              if((ptcorr>=70.)&&(ptcorr<85.)) {fh2DeltaRC60pt1->Fill(part->Pt(),deltaR);
-             if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R))fh2DeltaEtaC60pt1->Fill(part->Pt(),deltaEta);
-             if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC60pt1->Fill(part->Pt(),deltaPhi);}
-             if((ptcorr>=85.)&&(ptcorr<100.)) {fh2DeltaRC60pt2->Fill(part->Pt(),deltaR);
-              if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC60pt2->Fill(part->Pt(),deltaEta);
-              if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC60pt2->Fill(part->Pt(),deltaPhi); }
-             if((ptcorr>=100.)&&(ptcorr<120.)) {fh2DeltaRC60pt3->Fill(part->Pt(),deltaR);
-             if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC60pt3->Fill(part->Pt(),deltaEta);
-              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(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); }} 
 
+   }
+   if(injet>0) fh3JetDensity->Fill(ParticleList.GetEntries(),injet/accep,partback->Pt());
+   if(injet4>0)fh3JetDensityA4->Fill(ParticleList.GetEntries(),injet4/accep,partback->Pt());
+          //end of jet loop
+
+   //}
+
+
+          if(fDoEventMixing>0){
+            //check before if the trigger exists
+            // fTrigBuffer[i][0] = zvtx
+            // fTrigBuffer[i][1] = phi
+            // fTrigBuffer[i][2] = eta
+            // fTrigBuffer[i][3] = pt_jet
+            // fTrigBuffer[i][4] = pt_trig
+            // fTrigBuffer[i][5]= centrality
+            if(fTindex==10) fTindex=0;
+            if(fTrigBuffer[fTindex][3]>0){
+           if (TMath::Abs(fTrigBuffer[fTindex][0]-primVtx->GetZ()<2.)){
+           if (TMath::Abs(fTrigBuffer[fTindex][5]-centValue<5)){  
+             
+                        for(int it = 0;it<nT;++it){
+                       AliVParticle *part = (AliVParticle*)ParticleList.At(it);         
+                        Double_t DPhi = fTrigBuffer[fTindex][1] - part->Phi();
+                        Double_t DEta = fTrigBuffer[fTindex][2] - part->Eta();
+                        Double_t DR=TMath::Sqrt(DPhi*DPhi+DEta*DEta);
+                        if(DPhi<-0.5*TMath::Pi()) DPhi+=2.*TMath::Pi();
+                        if(DPhi>3./2.*TMath::Pi()) DPhi-=2.*TMath::Pi();
+                        Double_t triggerEntries[7] = {centValue,fTrigBuffer[fTindex][3],part->Pt(),DR,DEta,DPhi,fTrigBuffer[fTindex][4]};                      
+                        fhnMixedEvents->Fill(triggerEntries);
+                        }
+                        fNevents=fNevents+1;  
+                        if(fNevents==10) fTindex=fTindex+1; 
+           }}}
+
+              if(fTindex==10&&fNevents==10) fCountAgain=0;
+
+               // Copy the triggers from the current event into the buffer.
+               //again, only if the trigger exists:
+              if(fCountAgain==0){
+               if(trigJet>-1){
+                AliAODJet* jetT = (AliAODJet*)(fListJets[0]->At(trigJet));                      AliVParticle *partT = (AliVParticle*)ParticleList.At(trigBBTrack);         
+                fTrigBuffer[fTrigBufferIndex][0] = primVtx->GetZ();
+                fTrigBuffer[fTrigBufferIndex][1] = jetT->Phi();
+                fTrigBuffer[fTrigBufferIndex][2] = jetT->Eta();
+                fTrigBuffer[fTrigBufferIndex][3] = jetT->Pt()-rho*jetT->EffectiveAreaCharged();
+                fTrigBuffer[fTrigBufferIndex][4] = partT->Pt();
+                fTrigBuffer[fTrigBufferIndex][5] = centValue;
+                fTrigBufferIndex++;
+                if(fTrigBufferIndex==9) {fTrigBufferIndex=0; 
+                                        fCountAgain=1;}
+               }
+              }
+         
          }
-                 
-     } 
-     //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.){  
-     fh2SumPtInC10bkg->Fill(ptcorr,coronain/ptbig);       
-     fh2SumPtOutC10bkg->Fill(ptcorr,coronaout/ptbig);           
-     fh2SumPtInC10->Fill(ptcorr,sumPtIn/ptbig);
-     fh2SumPtOutC10->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); 
-      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);}
-
 
+         /////////////////////////////////////////////////////////////////////////////
+         ////////////////////// Rongrong's analysis //////////////////////////////////
+         if(fRunAnaAzimuthalCorrelation)
+           {
+             fhTTPt->Fill(centValue,partback->Pt());
+             for(Int_t ij=0; ij<fListJets[0]->GetEntries(); ij++)
+               {
+                 AliAODJet* jet = (AliAODJet*)(fListJets[0]->At(ij));
+                 Double_t jetPt   = jet->Pt();
+                 Double_t jetEta  = jet->Eta();
+                 Double_t jetPhi  = jet->Phi();
+                 if(jetPt==0) continue; 
+                 if((jetEta<fJetEtaMin)||(jetEta>fJetEtaMax)) continue;
+                 Double_t jetArea = jet->EffectiveAreaCharged();
+                 Double_t jetPtCorr=jetPt-rho*jetArea;
+                 Double_t dPhi=jetPhi-partback->Phi();
+                 if(dPhi>2*TMath::Pi()) dPhi -= 2*TMath::Pi();
+                 if(dPhi<-2*TMath::Pi()) dPhi += 2*TMath::Pi();
+                 if(dPhi<-0.5*TMath::Pi()) dPhi += 2*TMath::Pi();
+                 if(dPhi>1.5*TMath::Pi()) dPhi -= 2*TMath::Pi();
+
+                 Double_t fill[] = {partback->Pt(),jetPtCorr,dPhi,jetArea,centValue};
+                 fHJetPhiCorr->Fill(fill);
+               }
+           }
+         /////////////////////////////////////////////////////////////////////////////
+         /////////////////////////////////////////////////////////////////////////////
 
-     }
-     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); 
-
-     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);
-
-     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//////////////////////////////////////
 
      //tracks up to R=0.8 distant from the jet axis
-     if(fAngStructCloseTracks==1){
-      TList CloseTrackList;
-      Int_t nn=GetListOfTracksCloseToJet(&CloseTrackList,jetbig);
-      Double_t difR=0.04;
-      for(Int_t l=0;l<15;l++){
-      Double_t rr=l*0.1+0.1;
-       for(int it = 0;it<nn;++it){
-           AliVParticle *part1 = (AliVParticle*)CloseTrackList.At(it);
-           for(int itu=it+1;itu<CloseTrackList.GetEntries();itu++){      
-           AliVParticle *part2 = (AliVParticle*)CloseTrackList.At(itu);  
-           Double_t ptm=part1->Pt();
-           Double_t ptn=part2->Pt();   
-           Double_t Rnm = (part1->Eta()-part2->Eta())*(part1->Eta()-part2->Eta())+(part1->Phi()-part2->Phi())*(part1->Phi()-part2->Phi());
-                      Rnm=TMath::Sqrt(Rnm);
-                      Double_t deltag=(1./(TMath::Sqrt(2*TMath::Pi())*difR))*TMath::Exp(-1.*(rr-Rnm)*(rr-Rnm)/(2.*difR*difR));      
-                      Double_t stepf=0.5*(1.+TMath::Erf((rr-Rnm)/(TMath::Sqrt(2.)*difR)));
-                      if((ptcorr<85.) && (ptcorr>=70.)){up1[l]=up1[l]+ptm*ptn*Rnm*Rnm*deltag;
-                                                       down1[l]=down1[l]+ptm*ptn*Rnm*Rnm*stepf;}
-                      if((ptcorr<100.) && (ptcorr>=85.)){up2[l]=up2[l]+ptm*ptn*Rnm*Rnm*deltag;
-                                                       down2[l]=down2[l]+ptm*ptn*Rnm*Rnm*stepf;}  
-                      if((ptcorr<120.) && (ptcorr>=100.)){up3[l]=up3[l]+ptm*ptn*Rnm*Rnm*deltag;
-                                                        down3[l]=down3[l]+ptm*ptn*Rnm*Rnm*stepf;}
-                      if((ptcorr<140.) && (ptcorr>=120.)){up4[l]=up4[l]+ptm*ptn*Rnm*Rnm*deltag;
-                       down4[l]=down4[l]+ptm*ptn*Rnm*Rnm*stepf;}}}}
-     }
+   //   if(fAngStructCloseTracks==1){
+   //    TList CloseTrackList;
+   //    Int_t nn=GetListOfTracksCloseToJet(&CloseTrackList,jetbig);
+   //    Double_t difR=0.04;
+   //    for(Int_t l=0;l<15;l++){
+   //    Double_t rr=l*0.1+0.1;
+   //     for(int it = 0;it<nn;++it){
+   //         AliVParticle *part1 = (AliVParticle*)CloseTrackList.At(it);
+   //         for(int itu=it+1;itu<CloseTrackList.GetEntries();itu++){      
+   //         AliVParticle *part2 = (AliVParticle*)CloseTrackList.At(itu);  
+   //         Double_t ptm=part1->Pt();
+   //         Double_t ptn=part2->Pt();        
+   //         Double_t Rnm = (part1->Eta()-part2->Eta())*(part1->Eta()-part2->Eta())+(part1->Phi()-part2->Phi())*(part1->Phi()-part2->Phi());
+   //                    Rnm=TMath::Sqrt(Rnm);
+   //                    Double_t deltag=(1./(TMath::Sqrt(2*TMath::Pi())*difR))*TMath::Exp(-1.*(rr-Rnm)*(rr-Rnm)/(2.*difR*difR));      
+   //                    Double_t stepf=0.5*(1.+TMath::Erf((rr-Rnm)/(TMath::Sqrt(2.)*difR)));
+   //                    if((ptcorr<85.) && (ptcorr>=70.)){up1[l]=up1[l]+ptm*ptn*Rnm*Rnm*deltag;
+   //                                                  down1[l]=down1[l]+ptm*ptn*Rnm*Rnm*stepf;}
+   //                    if((ptcorr<100.) && (ptcorr>=85.)){up2[l]=up2[l]+ptm*ptn*Rnm*Rnm*deltag;
+   //                                                  down2[l]=down2[l]+ptm*ptn*Rnm*Rnm*stepf;}  
+   //                    if((ptcorr<120.) && (ptcorr>=100.)){up3[l]=up3[l]+ptm*ptn*Rnm*Rnm*deltag;
+   //                                                   down3[l]=down3[l]+ptm*ptn*Rnm*Rnm*stepf;}
+   //                    if((ptcorr<140.) && (ptcorr>=120.)){up4[l]=up4[l]+ptm*ptn*Rnm*Rnm*deltag;
+   //                  down4[l]=down4[l]+ptm*ptn*Rnm*Rnm*stepf;}}}}
+   //   }
     
-     //only jet constituents
-      if(fAngStructCloseTracks==2){
+   //   //only jet constituents
+   //    if(fAngStructCloseTracks==2){
 
-      Double_t difR=0.04;
-      for(Int_t l=0;l<15;l++){
-      Double_t rr=l*0.1+0.1;
+   //    Double_t difR=0.04;
+   //    for(Int_t l=0;l<15;l++){
+   //    Double_t rr=l*0.1+0.1;
 
     
-      AliAODTrack* part1;
-      AliAODTrack* part2;
-          for(Int_t it=0; it<nTracksGenJet; ++it){
-             part1 = (AliAODTrack*)(genTrackList->At(it));
-           for(Int_t itu=0; itu<nTracksGenJet; ++itu){
-             part2 = (AliAODTrack*)(genTrackList->At(itu));
-           Double_t ptm=part1->Pt();
-           Double_t ptn=part2->Pt();   
-           Double_t Rnm = (part1->Eta()-part2->Eta())*(part1->Eta()-part2->Eta())+(part1->Phi()-part2->Phi())*(part1->Phi()-part2->Phi());
-                      Rnm=TMath::Sqrt(Rnm);
-                      Double_t deltag=(1./(TMath::Sqrt(2*TMath::Pi())*difR))*TMath::Exp(-1.*(rr-Rnm)*(rr-Rnm)/(2.*difR*difR));
-                      Double_t stepf=0.5*(1.+TMath::Erf((rr-Rnm)/(TMath::Sqrt(2.)*difR))); 
-                      if((ptcorr<85.) && (ptcorr>=70.)){up1[l]=up1[l]+ptm*ptn*Rnm*Rnm*deltag;
-                                                       down1[l]=down1[l]+ptm*ptn*Rnm*Rnm*stepf;}
-                      if((ptcorr<100.) && (ptcorr>=85.)){up2[l]=up2[l]+ptm*ptn*Rnm*Rnm*deltag;
-                                                       down2[l]=down2[l]+ptm*ptn*Rnm*Rnm*stepf;}  
-                      if((ptcorr<120.) && (ptcorr>=100.)){up3[l]=up3[l]+ptm*ptn*Rnm*Rnm*deltag;
-                                                        down3[l]=down3[l]+ptm*ptn*Rnm*Rnm*stepf;}
-                      if((ptcorr<140.) && (ptcorr>=120.)){up4[l]=up4[l]+ptm*ptn*Rnm*Rnm*deltag;
-                       down4[l]=down4[l]+ptm*ptn*Rnm*Rnm*stepf;}}}}}
+   //    AliAODTrack* part1;
+   //    AliAODTrack* part2;
+
+   //        TRefArray *genTrackListb = jetbig->GetRefTracks();
+   //        Int_t nTracksGenJetb = genTrackListb->GetEntriesFast();
+          
+             
+
+   //        for(Int_t it=0; it<nTracksGenJetb; ++it){
+   //           part1 = (AliAODTrack*)(genTrackListb->At(it));
+   //         for(Int_t itu=0; itu<nTracksGenJetb; ++itu){
+   //           part2 = (AliAODTrack*)(genTrackListb->At(itu));
+   //         Double_t ptm=part1->Pt();
+   //         Double_t ptn=part2->Pt();        
+   //         Double_t Rnm = (part1->Eta()-part2->Eta())*(part1->Eta()-part2->Eta())+(part1->Phi()-part2->Phi())*(part1->Phi()-part2->Phi());
+   //                    Rnm=TMath::Sqrt(Rnm);
+   //                    Double_t deltag=(1./(TMath::Sqrt(2*TMath::Pi())*difR))*TMath::Exp(-1.*(rr-Rnm)*(rr-Rnm)/(2.*difR*difR));
+   //                    Double_t stepf=0.5*(1.+TMath::Erf((rr-Rnm)/(TMath::Sqrt(2.)*difR))); 
+   //                    if((ptcorr<85.) && (ptcorr>=70.)){up1[l]=up1[l]+ptm*ptn*Rnm*Rnm*deltag;
+   //                                                  down1[l]=down1[l]+ptm*ptn*Rnm*Rnm*stepf;}
+   //                    if((ptcorr<100.) && (ptcorr>=85.)){up2[l]=up2[l]+ptm*ptn*Rnm*Rnm*deltag;
+   //                                                  down2[l]=down2[l]+ptm*ptn*Rnm*Rnm*stepf;}  
+   //                    if((ptcorr<120.) && (ptcorr>=100.)){up3[l]=up3[l]+ptm*ptn*Rnm*Rnm*deltag;
+   //                                                   down3[l]=down3[l]+ptm*ptn*Rnm*Rnm*stepf;}
+   //                    if((ptcorr<140.) && (ptcorr>=120.)){up4[l]=up4[l]+ptm*ptn*Rnm*Rnm*deltag;
+   //                  down4[l]=down4[l]+ptm*ptn*Rnm*Rnm*stepf;}}}}}
+   // }
+     // //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){  
+     //        if(centValue<10.)fh2AngStructpt1C10->Fill(rr,rr*up1[l]/down1[l]);
+     //    if(centValue>20. && centValue<40.) fh2AngStructpt1C20->Fill(rr,rr*up1[l]/down1[l]);
+     //    if(centValue>30. && centValue<60.) fh2AngStructpt1C30->Fill(rr,rr*up1[l]/down1[l]);
+     //    if(centValue>60.) fh2AngStructpt1C60->Fill(rr,rr*up1[l]/down1[l]);}
+     //    if(down2[l]!=0){  
+     //        if(centValue<10.) fh2AngStructpt2C10->Fill(rr,rr*up2[l]/down2[l]);
+     //    if(centValue>20. && centValue<40.) fh2AngStructpt2C20->Fill(rr,rr*up2[l]/down2[l]);
+     //    if(centValue>30. && centValue<60.) fh2AngStructpt2C30->Fill(rr,rr*up2[l]/down2[l]);
+     //    if(centValue>60.) fh2AngStructpt2C60->Fill(rr,rr*up2[l]/down2[l]);}
+     //    if(down3[l]!=0){  
+     //        if(centValue<10.) fh2AngStructpt3C10->Fill(rr,rr*up3[l]/down3[l]);
+     //    if(centValue>20. && centValue<40.) fh2AngStructpt3C20->Fill(rr,rr*up3[l]/down3[l]);
+     //    if(centValue>30. && centValue<60.) fh2AngStructpt3C30->Fill(rr,rr*up3[l]/down3[l]);
+     //    if(centValue>60.) fh2AngStructpt3C60->Fill(rr,rr*up3[l]/down3[l]);}
+     //    if(down4[l]!=0){  
+     //        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]);}}}
 
+    
+
+   }
 
-   
 
 
+   PostData(1, fOutputList);
+}
 
+void AliAnalysisTaskJetCore::Terminate(const Option_t *)
+{
+   // Draw result to the screen
+   // Called once at the end of the query
 
+   if (!GetOutputData(1))
+   return;
+}
 
 
 
+  
 
 
+Int_t  AliAnalysisTaskJetCore::GetListOfTracks(TList *list){
 
+      Int_t iCount = 0;
+      AliAODEvent *aod = 0;
 
+     if(!fESD)aod = fAODIn;
+     else aod = fAODOut;   
 
-   }
+     if(!aod)return 0;
 
+     Int_t index=-1;
+     Double_t ptmax=-10;
 
-     //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){  
-       if(centValue<10.)fh2AngStructpt1C10->Fill(rr,rr*up1[l]/down1[l]);
-        if(centValue>20. && centValue<40.) fh2AngStructpt1C20->Fill(rr,rr*up1[l]/down1[l]);
-        if(centValue>30. && centValue<60.) fh2AngStructpt1C30->Fill(rr,rr*up1[l]/down1[l]);
-        if(centValue>60.) fh2AngStructpt1C60->Fill(rr,rr*up1[l]/down1[l]);}
-        if(down2[l]!=0){  
-       if(centValue<10.) fh2AngStructpt2C10->Fill(rr,rr*up2[l]/down2[l]);
-        if(centValue>20. && centValue<40.) fh2AngStructpt2C20->Fill(rr,rr*up2[l]/down2[l]);
-        if(centValue>30. && centValue<60.) fh2AngStructpt2C30->Fill(rr,rr*up2[l]/down2[l]);
-        if(centValue>60.) fh2AngStructpt2C60->Fill(rr,rr*up2[l]/down2[l]);}
-        if(down3[l]!=0){  
-       if(centValue<10.) fh2AngStructpt3C10->Fill(rr,rr*up3[l]/down3[l]);
-        if(centValue>20. && centValue<40.) fh2AngStructpt3C20->Fill(rr,rr*up3[l]/down3[l]);
-        if(centValue>30. && centValue<60.) fh2AngStructpt3C30->Fill(rr,rr*up3[l]/down3[l]);
-        if(centValue>60.) fh2AngStructpt3C60->Fill(rr,rr*up3[l]/down3[l]);}
-        if(down4[l]!=0){  
-       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]);}}}
 
     
+     for(int it = 0;it < aod->GetNumberOfTracks();++it){
+      AliAODTrack *tr = dynamic_cast<AliAODTrack*>(aod->GetTrack(it));
+      if(!tr) AliFatal("Not a standard AOD");
+      Bool_t bGood = false;
+      if(fFilterType == 0)bGood = true;
+      else if(fFilterType == 1)bGood = tr->IsHybridTPCConstrainedGlobal();
+      else if(fFilterType == 2)bGood = tr->IsHybridGlobalConstrainedGlobal();    
+      if((fFilterMask>0)&&!(tr->TestFilterBit(fFilterMask)))continue;
+      if(fRequireITSRefit==1){if((tr->GetStatus()&AliESDtrack::kITSrefit)==0)continue;}
+      if(bGood==false) continue;
+      if (fApplySharedClusterCut) {
+           Double_t frac = Double_t(tr->GetTPCnclsS()) /Double_t(tr->GetTPCncls());
+           if (frac > 0.4) continue;
+      }
+     if(TMath::Abs(tr->Eta())>0.9)continue;
+      if(tr->Pt()<0.15)continue;
+      list->Add(tr);
+      iCount++;
+      if(fFilterType==2 && fFilterMaskBestPt>0){// only set the trigger track index for good quality tracks
+       if(tr->TestFilterBit(fFilterMaskBestPt)){
+         if(tr->Pt()>ptmax){ 
+           ptmax=tr->Pt();     
+           index=iCount-1;
+         }
+       }
+      }
+      else{
+       if(tr->Pt()>ptmax){ 
+         ptmax=tr->Pt();       
+         index=iCount-1;
+       }
+      }
+     }
+  
+   
+    // else if (type == kTrackAODMCCharged) {
+    // TClonesArray *tca = dynamic_cast<TClonesArray*>(aod->FindListObject(AliAODMCParticle::StdBranchName()));
+    // if(!tca)return iCount;
+    // for(int it = 0;it < tca->GetEntriesFast();++it){
+    //   AliAODMCParticle *part = dynamic_cast<AliAODMCParticle*>(tca->At(it));
+    //   if(!part)continue;
+    //   if(part->Pt()<0.15)continue;
+    //   if(!part->IsPhysicalPrimary())continue;
+    //   if(part->Charge()==0)continue;
+    //   if(TMath::Abs(part->Eta())>0.9)continue;
+    //   list->Add(part);
+    //   iCount++;
+    //   if(part->Pt()>ptmax){ ptmax=part->Pt();
+    //         index=iCount-1;}}}
+      return index;
+ 
+}
 
 
 
+Int_t  AliAnalysisTaskJetCore::SelectTrigger(TList *list,Double_t minT,Double_t maxT,Int_t &number){
+     Int_t iCount = 0;
+     AliAODEvent *aod = 0;
+     if(!fESD)aod = fAODIn;
+     else aod = fAODOut;   
+     if(!aod)return 0;
+     Int_t index=-1;
+     Int_t triggers[100];
+     Int_t triggers2[100];
+     for(Int_t cr=0;cr<100;cr++){triggers[cr]=-1;
+       triggers2[cr]=-1;}
+     Int_t im=0;
+     Int_t im2=0;
+     for(int it = 0;it < aod->GetNumberOfTracks();++it){
+      AliAODTrack *tr = dynamic_cast<AliAODTrack*>(aod->GetTrack(it));
+      if(!tr) AliFatal("Not a standard AOD");
+      Bool_t bGood = false;
+      if(fFilterType == 0)bGood = true;
+      else if(fFilterType == 1)bGood = tr->IsHybridTPCConstrainedGlobal();
+      else if(fFilterType == 2)bGood = tr->IsHybridGlobalConstrainedGlobal();
+      if((fFilterMask>0)&&!(tr->TestFilterBit(fFilterMask)))continue;
+       if(fRequireITSRefit==1){if((tr->GetStatus()&AliESDtrack::kITSrefit)==0)continue;}
+      if(bGood==false) continue;
+      if (fApplySharedClusterCut) {
+           Double_t frac = Double_t(tr->GetTPCnclsS()) /Double_t(tr->GetTPCncls());
+           if (frac > 0.4) continue;
+      }
+      if(TMath::Abs(tr->Eta())>0.9)continue;
+      if(tr->Pt()<0.15)continue;
+      list->Add(tr);
+      iCount++;
+      
+      if(tr->Pt()>=minT && tr->Pt()<maxT){
+       triggers[im]=iCount-1;
+        im=im+1;
+       if(tr->Pt()>=20.){triggers2[im2]=iCount-1;
+         im2=im2+1;}
 
+      }}
+      
+      number=im;
+      Int_t rd=0;
+      if(im2==0) rd=0;
+      if(im2>0) rd=fRandom->Integer(im2);
+      index=triggers2[rd];
+      AliVParticle *tr1 = (AliVParticle*)list->At(index);     
+      
+    
 
-   PostData(1, fOutputList);
-     }
-
-void AliAnalysisTaskJetCore::Terminate(const Option_t *)
-{
-   // Draw result to the screen
-   // Called once at the end of the query
+      for(Int_t kk=0;kk<im;kk++){
+       //if(kk==rd) continue;
+        if(index==triggers[kk]) continue;
+       Int_t lab=triggers[kk];
+         AliVParticle *tr2 = (AliVParticle*)list->At(lab);     
+       
+       Double_t detat=tr1->Eta()-tr2->Eta();
+       Double_t dphit=RelativePhi(tr1->Phi(),tr2->Phi());
+       Double_t deltaRt=TMath::Sqrt(detat*detat+dphit*dphit);
+     
+        if(deltaRt>0.4) number=number-1;}      
 
-   if (!GetOutputData(1))
-   return;
+     
+  
+   
+  
+      return index;
+ 
 }
 
 
@@ -1842,33 +1299,57 @@ void AliAnalysisTaskJetCore::Terminate(const Option_t *)
 
 
 
-Int_t  AliAnalysisTaskJetCore::GetListOfTracks(TList *list){
 
-    Int_t iCount = 0;
+
+
+
+
+   Int_t  AliAnalysisTaskJetCore::GetHardestTrackBackToJet(AliAODJet *jetbig){
  
-    
-    for(int it = 0;it < fAOD->GetNumberOfTracks();++it){
-      AliAODTrack *tr = fAOD->GetTrack(it);
+    AliAODEvent *aod = 0;
+    if(!fESD)aod = fAODIn;
+    else aod = fAODOut;     
+    Int_t index=-1;
+    Double_t ptmax=-10;
+    Double_t dphi=0;
+    // Double_t dif=0;
+    Int_t iCount=0;
+    for(int it = 0;it < aod->GetNumberOfTracks();++it){
+      AliAODTrack *tr = dynamic_cast<AliAODTrack*>(aod->GetTrack(it));
+      if(!tr) AliFatal("Not a standard AOD");
       if((fFilterMask>0)&&!(tr->TestFilterBit(fFilterMask)))continue;
       if(TMath::Abs(tr->Eta())>0.9)continue;
       if(tr->Pt()<0.15)continue;
-      list->Add(tr);
-      //cout<<fAOD->GetNumberOfTracks()<<" "<<tr->Pt()<<endl;
-      iCount++;
-    }
+      iCount=iCount+1;
+      dphi=RelativePhi(tr->Phi(),jetbig->Phi());  
+      if(TMath::Abs(dphi)<TMath::Pi()-0.6) continue;
+      if(tr->Pt()>ptmax){ ptmax=tr->Pt();
+      index=iCount-1;
+      //  dif=dphi; 
+      }}
   
-   list->Sort();
-  return iCount;
- 
-}
+      return index;
+
+   }
+
+
+
+
+
+
+
+
 
  Int_t  AliAnalysisTaskJetCore::GetListOfTracksCloseToJet(TList *list,AliAODJet *jetbig){
 
     Int_t iCount = 0;
- 
+      AliAODEvent *aod = 0;
+     if(!fESD)aod = fAODIn;
+     else aod = fAODOut;   
   
-    for(int it = 0;it < fAOD->GetNumberOfTracks();++it){
-      AliAODTrack *tr = fAOD->GetTrack(it);
+      for(int it = 0;it < aod->GetNumberOfTracks();++it){
+      AliAODTrack *tr = dynamic_cast<AliAODTrack*>(aod->GetTrack(it));
+      if(!tr) AliFatal("Not a standard AOD");
       if((fFilterMask>0)&&!(tr->TestFilterBit(fFilterMask)))continue;
       if(TMath::Abs(tr->Eta())>0.9)continue;
       if(tr->Pt()<0.15)continue;
@@ -1898,7 +1379,9 @@ Int_t AliAnalysisTaskJetCore::GetNInputTracks()
 {
 
    Int_t nInputTracks = 0;
-
+     AliAODEvent *aod = 0;
+     if(!fESD)aod = fAODIn;
+     else aod = fAODOut;   
    TString jbname(fJetBranchName[1]);
    //needs complete event, use jets without background subtraction
    for(Int_t i=1; i<=3; ++i){
@@ -1909,7 +1392,7 @@ Int_t AliAnalysisTaskJetCore::GetNInputTracks()
    if(jbname.Contains("AODextra")) jbname.ReplaceAll("AODextra","AOD");
 
    if(fDebug) Printf("Multiplicity from jet branch %s", jbname.Data());
-   TClonesArray *tmpAODjets = dynamic_cast<TClonesArray*>(fAOD->FindListObject(jbname.Data()));
+   TClonesArray *tmpAODjets = dynamic_cast<TClonesArray*>(aod->FindListObject(jbname.Data()));
    if(!tmpAODjets){
       Printf("Jet branch %s not found", jbname.Data());
       Printf("AliAnalysisTaskJetCore::GetNInputTracks FAILED");
@@ -1943,5 +1426,137 @@ Double_t AliAnalysisTaskJetCore::RelativePhi(Double_t mphi,Double_t vphi){
   return dphi;//dphi in [-Pi, Pi]
 }
 
+Int_t AliAnalysisTaskJetCore::GetPhiBin(Double_t phi)
+{
+    Int_t phibin=-1;
+    if(!(TMath::Abs(phi)<=2*TMath::Pi())){AliError("phi w.r.t. RP out of defined range");return -1;}
+    Double_t phiwrtrp=TMath::ACos(TMath::Abs(TMath::Cos(phi)));
+    phibin=Int_t(fNRPBins*phiwrtrp/(0.5*TMath::Pi()));
+    if(phibin<0||phibin>=fNRPBins){AliError("Phi Bin not defined");}
+    return phibin;
+}
+
+
+
+
+THnSparse* AliAnalysisTaskJetCore::NewTHnSparseF(const char* name, UInt_t entries)
+{
+   // generate new THnSparseF, axes are defined in GetDimParams()
+
+   Int_t count = 0;
+   UInt_t tmp = entries;
+   while(tmp!=0){
+      count++;
+      tmp = tmp &~ -tmp;  // clear lowest bit
+   }
+
+   TString hnTitle(name);
+   const Int_t dim = count;
+   Int_t nbins[dim];
+   Double_t xmin[dim];
+   Double_t xmax[dim];
+
+   Int_t i=0;
+   Int_t c=0;
+   while(c<dim && i<32){
+      if(entries&(1<<i)){
+      
+         TString label("");
+         GetDimParams(i, label, nbins[c], xmin[c], xmax[c]);
+         hnTitle += Form(";%s",label.Data());
+         c++;
+      }
+      
+      i++;
+   }
+   hnTitle += ";";
+
+   return new THnSparseF(name, hnTitle.Data(), dim, nbins, xmin, xmax);
+}
+
+void AliAnalysisTaskJetCore::GetDimParams(Int_t iEntry, TString &label, Int_t &nbins, Double_t &xmin, Double_t &xmax)
+{
+   // stores label and binning of axis for THnSparse
+
+   const Double_t pi = TMath::Pi();
+   
+   switch(iEntry){
+      
+   case 0:
+      label = "V0 centrality (%)";
+     
+         nbins = 10;
+         xmin = 0.;
+         xmax = 100.;
+         break;
+      
+      
+   case 1:
+      label = "corrected jet pt";
+         nbins = 20;
+         xmin = 0.;
+         xmax = 200.;
+          break;
+      
+      
+   case 2:
+      label = "track pT";
+     
+         nbins = 9;
+         xmin = 0.;
+         xmax = 150;
+         break;
+      
+      
+    case 3:
+      label = "deltaR";
+      nbins = 15;
+      xmin = 0.;
+      xmax = 1.5;
+      break;
+
 
 
+   case 4:
+      label = "deltaEta";
+      nbins = 8;
+      xmin = -1.6;
+      xmax = 1.6;
+      break;
+
+
+  case 5:
+      label = "deltaPhi";
+      nbins = 90;
+      xmin = -0.5*pi;
+      xmax = 1.5*pi;
+      break;   
+   
+      
+        
+    case 6:
+      label = "leading track";
+      nbins = 13;
+      xmin = 0;
+      xmax = 50;
+      break;
+           
+     case 7:
+    
+      label = "trigger track";
+      nbins =10;
+      xmin = 0;
+      xmax = 50;
+      break;
+
+
+   
+  
+
+
+
+
+   }
+
+}
+