#include "AliAnalysisTaskFastEmbedding.h"
#include "AliAODEvent.h"
+#include "AliAODHandler.h"
#include "AliAODJet.h"
#include "AliAnalysisTaskJetCore.h"
AliAnalysisTaskSE(),
fESD(0x0),
fAOD(0x0),
+fAODExtension(0x0),
fBackgroundBranch(""),
+fNonStdFile(""),
fIsPbPb(kTRUE),
fOfflineTrgMask(AliVEvent::kAny),
fMinContribVtx(1),
fVtxZMax(8.),
fEvtClassMin(0),
fEvtClassMax(4),
+fRadioFrac(0.2),
+fMinDist(0.1),
fCentMin(0.),
fCentMax(100.),
fNInputTracksMin(0),
fNInputTracksMax(-1),
+fAngStructCloseTracks(0),
fJetEtaMin(-.5),
fJetEtaMax(.5),
fJetPtMin(20.),
fNMatchJets(4),
fMatchMaxDist(0.8),
fKeepJets(kFALSE),
-fRadioFrac(0.2),
-fMinDist(0.1),
fkNbranches(2),
fkEvtClasses(12),
fOutputList(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)
-
-
+fh2JetCoreMethod3C60(0x0),
+fh2SumPtInC10(0x0),
+fh2SumPtInC20(0x0),
+fh2SumPtInC30(0x0),
+fh2SumPtInC60(0x0),
+fh2SumPtOutC10(0x0),
+fh2SumPtOutC10b(0x0),
+fh2SumPtOutC20(0x0),
+fh2SumPtOutC30(0x0),
+fh2SumPtOutC60(0x0),
+fh2SumPtInC10bkg(0x0),
+fh2SumPtInC20bkg(0x0),
+fh2SumPtInC30bkg(0x0),
+fh2SumPtInC60bkg(0x0),
+fh2SumPtOutC10bkg(0x0),
+fh2SumPtOutC20bkg(0x0),
+fh2SumPtOutC30bkg(0x0),
+fh2SumPtOutC60bkg(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),
+fh2AngStructpt1C10(0x0),
+fh2AngStructpt2C10(0x0),
+fh2AngStructpt3C10(0x0),
+fh2AngStructpt4C10(0x0),
+fh2AngStructpt1C20(0x0),
+fh2AngStructpt2C20(0x0),
+fh2AngStructpt3C20(0x0),
+fh2AngStructpt4C20(0x0),
+fh2AngStructpt1C30(0x0),
+fh2AngStructpt2C30(0x0),
+fh2AngStructpt3C30(0x0),
+fh2AngStructpt4C30(0x0),
+fh2AngStructpt1C60(0x0),
+fh2AngStructpt2C60(0x0),
+fh2AngStructpt3C60(0x0),
+fh2AngStructpt4C60(0x0)
{
// default Constructor
AliAnalysisTaskSE(name),
fESD(0x0),
fAOD(0x0),
+fAODExtension(0x0),
fBackgroundBranch(""),
+fNonStdFile(""),
fIsPbPb(kTRUE),
fOfflineTrgMask(AliVEvent::kAny),
fMinContribVtx(1),
fVtxZMax(8.),
fEvtClassMin(0),
fEvtClassMax(4),
+fRadioFrac(0.2),
+fMinDist(0.1),
fCentMin(0.),
fCentMax(100.),
fNInputTracksMin(0),
fNInputTracksMax(-1),
+fAngStructCloseTracks(0),
fJetEtaMin(-.5),
fJetEtaMax(.5),
fJetPtMin(20.),
fNMatchJets(4),
fMatchMaxDist(0.8),
fKeepJets(kFALSE),
-fRadioFrac(0.2),
-fMinDist(0.1),
fkNbranches(2),
fkEvtClasses(12),
fOutputList(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)
-
+fh2JetCoreMethod3C60(0x0),
+fh2SumPtInC10(0x0),
+fh2SumPtInC20(0x0),
+fh2SumPtInC30(0x0),
+fh2SumPtInC60(0x0),
+fh2SumPtOutC10(0x0),
+fh2SumPtOutC10b(0x0),
+fh2SumPtOutC20(0x0),
+fh2SumPtOutC30(0x0),
+fh2SumPtOutC60(0x0),
+fh2SumPtInC10bkg(0x0),
+fh2SumPtInC20bkg(0x0),
+fh2SumPtInC30bkg(0x0),
+fh2SumPtInC60bkg(0x0),
+fh2SumPtOutC10bkg(0x0),
+fh2SumPtOutC20bkg(0x0),
+fh2SumPtOutC30bkg(0x0),
+fh2SumPtOutC60bkg(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),
+fh2AngStructpt1C10(0x0),
+fh2AngStructpt2C10(0x0),
+fh2AngStructpt3C10(0x0),
+fh2AngStructpt4C10(0x0),
+fh2AngStructpt1C20(0x0),
+fh2AngStructpt2C20(0x0),
+fh2AngStructpt3C20(0x0),
+fh2AngStructpt4C20(0x0),
+fh2AngStructpt1C30(0x0),
+fh2AngStructpt2C30(0x0),
+fh2AngStructpt3C30(0x0),
+fh2AngStructpt4C30(0x0),
+fh2AngStructpt1C60(0x0),
+fh2AngStructpt2C60(0x0),
+fh2AngStructpt3C60(0x0),
+fh2AngStructpt4C60(0x0)
{
// Constructor
fh2JetSelection->GetXaxis()->SetBinLabel(8,"trigger exclude mask");
- // UInt_t entries = 0; // bit coded, see GetDimParams() below
- // UInt_t opt = 0; // bit coded, default (0) or high resolution (1)
+ //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);
+
+
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);
-
+
+ fh2SumPtInC10 = new TH2F("PtSumInC10","",150, 0., 150.,5000, 0., 50.);
+ fh2SumPtInC20 = new TH2F("PtSumInC20","",150, 0., 150.,5000, 0., 50.);
+ fh2SumPtInC30 = new TH2F("PtSumInC30","",150, 0., 150.,5000, 0., 50.);
+ fh2SumPtInC60 = new TH2F("PtSumInC60","",150, 0., 150.,5000, 0., 50.);
+ fh2SumPtOutC10 = new TH2F("PtSumOutC10","",150, 0., 150.,5000, 0., 50.);
+ fh2SumPtOutC10b= new TH2F("PtSumOutC10b","",150, 0., 150.,500, 0., 50.);
+ fh2SumPtOutC20 = new TH2F("PtSumOutC20","",150, 0., 150.,5000, 0., 50.);
+ fh2SumPtOutC30 = new TH2F("PtSumOutC30","",150, 0., 150.,5000, 0., 50.);
+ fh2SumPtOutC60 = new TH2F("PtSumOutC60","",150, 0., 150.,5000, 0., 50.);
+
+ fh2SumPtInC10bkg = new TH2F("PtSumInC10bkg","",150, 0., 150.,5000, 0., 50.);
+ fh2SumPtInC20bkg = new TH2F("PtSumInC20bkg","",150, 0., 150.,5000, 0., 50.);
+ fh2SumPtInC30bkg = new TH2F("PtSumInC30bkg","",150, 0., 150.,5000, 0., 50.);
+ fh2SumPtInC60bkg = new TH2F("PtSumInC60bkg","",150, 0., 150.,5000, 0., 50.);
+ fh2SumPtOutC10bkg = new TH2F("PtSumOutC10bkg","",150, 0., 150.,5000, 0., 50.);
+ fh2SumPtOutC20bkg = new TH2F("PtSumOutC20bkg","",150, 0., 150.,5000, 0., 50.);
+ fh2SumPtOutC30bkg = new TH2F("PtSumOutC30bkg","",150, 0., 150.,5000, 0., 50.);
+ fh2SumPtOutC60bkg = new TH2F("PtSumOutC60bkg","",150, 0., 150.,5000, 0., 50.);
+
+
+
+ fh2DeltaRC10pt1 = new TH2F("DeltaRC10pt1","",150, 0., 150.,100,0.,1.5);
+ fh2DeltaRC20pt1 = new TH2F("DeltaRC20pt1","",150, 0., 150.,100,0.,1.5);
+ fh2DeltaRC30pt1 = new TH2F("DeltaRC30pt1","",150, 0., 150.,100,0.,1.5);
+ fh2DeltaRC60pt1 = new TH2F("DeltaRC60pt1","",150, 0., 150.,100,0.,1.5);
+ fh2DeltaRC10pt2 = new TH2F("DeltaRC10pt2","",150, 0., 150.,100,0.,1.5);
+ fh2DeltaRC20pt2 = new TH2F("DeltaRC20pt2","",150, 0., 150.,100,0.,1.5);
+ fh2DeltaRC30pt2 = new TH2F("DeltaRC30pt2","",150, 0., 150.,100,0.,1.5);
+ fh2DeltaRC60pt2 = new TH2F("DeltaRC60pt2","",150, 0., 150.,100,0.,1.5);
+ fh2DeltaRC10pt3 = new TH2F("DeltaRC10pt3","",150, 0., 150.,100,0.,1.5);
+ fh2DeltaRC20pt3 = new TH2F("DeltaRC20pt3","",150, 0., 150.,100,0.,1.5);
+ fh2DeltaRC30pt3 = new TH2F("DeltaRC30pt3","",150, 0., 150.,100,0.,1.5);
+ fh2DeltaRC60pt3 = new TH2F("DeltaRC60pt3","",150, 0., 150.,100,0.,1.5);
+ fh2DeltaRC10pt4 = new TH2F("DeltaRC10pt4","",150, 0., 150.,100,0.,1.5);
+ fh2DeltaRC20pt4 = new TH2F("DeltaRC20pt4","",150, 0., 150.,100,0.,1.5);
+ fh2DeltaRC30pt4 = new TH2F("DeltaRC30pt4","",150, 0., 150.,100,0.,1.5);
+ fh2DeltaRC60pt4 = new TH2F("DeltaRC60pt4","",150, 0., 150.,100,0.,1.5);
+
+ fh2DeltaEtaC10pt1 = new TH2F("DeltaEtaC10pt1","",150, 0., 150.,100,-1.5,1.5);
+ fh2DeltaEtaC20pt1 = new TH2F("DeltaEtaC20pt1","",150, 0., 150.,100,-1.5,1.5);
+ fh2DeltaEtaC30pt1 = new TH2F("DeltaEtaC30pt1","",150, 0., 150.,100,-1.5,1.5);
+ fh2DeltaEtaC60pt1 = new TH2F("DeltaEtaC60pt1","",150, 0., 150.,100,-1.5,1.5);
+ fh2DeltaEtaC10pt2 = new TH2F("DeltaEtaC10pt2","",150, 0., 150.,100,-1.5,1.5);
+ fh2DeltaEtaC20pt2 = new TH2F("DeltaEtaC20pt2","",150, 0., 150.,100,-1.5,1.5);
+ fh2DeltaEtaC30pt2 = new TH2F("DeltaEtaC30pt2","",150, 0., 150.,100,-1.5,1.5);
+ fh2DeltaEtaC60pt2 = new TH2F("DeltaEtaC60pt2","",150, 0., 150.,100,-1.5,1.5);
+ fh2DeltaEtaC10pt3 = new TH2F("DeltaEtaC10pt3","",150, 0., 150.,100,-1.5,1.5);
+ fh2DeltaEtaC20pt3 = new TH2F("DeltaEtaC20pt3","",150, 0., 150.,100,-1.5,1.5);
+ fh2DeltaEtaC30pt3 = new TH2F("DeltaEtaC30pt3","",150, 0., 150.,100,-1.5,1.5);
+ fh2DeltaEtaC60pt3 = new TH2F("DeltaEtaC60pt3","",150, 0., 150.,100,-1.5,1.5);
+ fh2DeltaEtaC10pt4 = new TH2F("DeltaEtaC10pt4","",150, 0., 150.,100,-1.5,1.5);
+ fh2DeltaEtaC20pt4 = new TH2F("DeltaEtaC20pt4","",150, 0., 150.,100,-1.5,1.5);
+ fh2DeltaEtaC30pt4 = new TH2F("DeltaEtaC30pt4","",150, 0., 150.,100,-1.5,1.5);
+ fh2DeltaEtaC60pt4 = new TH2F("DeltaEtaC60pt4","",150, 0., 150.,100,-1.5,1.5);
+ fh2DeltaPhiC10pt1 = new TH2F("DeltaPhiC10pt1","",150, 0., 150.,100,-6.5,6.5);
+ fh2DeltaPhiC20pt1 = new TH2F("DeltaPhiC20pt1","",150, 0., 150.,100,-6.5,6.5);
+ fh2DeltaPhiC30pt1 = new TH2F("DeltaPhiC30pt1","",150, 0., 150.,100,-6.5,6.5);
+ fh2DeltaPhiC60pt1 = new TH2F("DeltaPhiC60pt1","",150, 0., 150.,100,-6.5,6.5);
+ fh2DeltaPhiC10pt2 = new TH2F("DeltaPhiC10pt2","",150, 0., 150.,100,-6.5,6.5);
+ fh2DeltaPhiC20pt2 = new TH2F("DeltaPhiC20pt2","",150, 0., 150.,100,-6.5,6.5);
+ fh2DeltaPhiC30pt2 = new TH2F("DeltaPhiC30pt2","",150, 0., 150.,100,-6.5,6.5);
+ fh2DeltaPhiC60pt2 = new TH2F("DeltaPhiC60pt2","",150, 0., 150.,100,-6.5,6.5);
+ fh2DeltaPhiC10pt3 = new TH2F("DeltaPhiC10pt3","",150, 0., 150.,100,-6.5,6.5);
+ fh2DeltaPhiC20pt3 = new TH2F("DeltaPhiC20pt3","",150, 0., 150.,100,-6.5,6.5);
+ fh2DeltaPhiC30pt3 = new TH2F("DeltaPhiC30pt3","",150, 0., 150.,100,-6.5,6.5);
+ fh2DeltaPhiC60pt3 = new TH2F("DeltaPhiC60pt3","",150, 0., 150.,100,-6.5,6.5);
+ fh2DeltaPhiC10pt4 = new TH2F("DeltaPhiC10pt4","",150, 0., 150.,100,-6.5,6.5);
+ fh2DeltaPhiC20pt4 = new TH2F("DeltaPhiC20pt4","",150, 0., 150.,100,-6.5,6.5);
+ fh2DeltaPhiC30pt4 = new TH2F("DeltaPhiC30pt4","",150, 0., 150.,100,-6.5,6.5);
+ fh2DeltaPhiC60pt4 = new TH2F("DeltaPhiC60pt4","",150, 0., 150.,100,-6.5,6.5);
+
+ fh2AngStructpt1C10 = new TH2F("Ang struct pt1 C10","",15,0.,1.5,150,0.,5.);
+ fh2AngStructpt2C10 = new TH2F("Ang struct pt2 C10","",15,0.,1.5,150,0.,5.);
+ fh2AngStructpt3C10 = new TH2F("Ang struct pt3 C10","",15,0.,1.5,150,0.,5.);
+ fh2AngStructpt4C10 = new TH2F("Ang struct pt4 C10","",15,0.,1.5,150,0.,5.);
+ fh2AngStructpt1C20 = new TH2F("Ang struct pt1 C20","",15,0.,1.5,150,0.,5.);
+ fh2AngStructpt2C20 = new TH2F("Ang struct pt2 C20","",15,0.,1.5,150,0.,5.);
+ fh2AngStructpt3C20 = new TH2F("Ang struct pt3 C20","",15,0.,1.5,150,0.,5.);
+ fh2AngStructpt4C20 = new TH2F("Ang struct pt4 C20","",15,0.,1.5,150,0.,5.);
+ fh2AngStructpt1C30 = new TH2F("Ang struct pt1 C30","",15,0.,1.5,150,0.,5.);
+ fh2AngStructpt2C30 = new TH2F("Ang struct pt2 C30","",15,0.,1.5,150,0.,5.);
+ fh2AngStructpt3C30 = new TH2F("Ang struct pt3 C30","",15,0.,1.5,150,0.,5.);
+ fh2AngStructpt4C30 = new TH2F("Ang struct pt4 C30","",15,0.,1.5,150,0.,5.);
+ fh2AngStructpt1C60 = new TH2F("Ang struct pt1 C60","",15,0.,1.5,150,0.,5.);
+ fh2AngStructpt2C60 = new TH2F("Ang struct pt2 C60","",15,0.,1.5,150,0.,5.);
+ fh2AngStructpt3C60 = new TH2F("Ang struct pt3 C60","",15,0.,1.5,150,0.,5.);
+ fh2AngStructpt4C60 = new TH2F("Ang struct pt4 C60","",15,0.,1.5,150,0.,5.);
fOutputList->Add(fHistEvtSelection);
fOutputList->Add(fHistJetSelection);
fOutputList->Add(fh2JetSelection);
-
-
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(fh2SumPtInC10);
+ fOutputList->Add(fh2SumPtInC20);
+ fOutputList->Add(fh2SumPtInC30);
+ fOutputList->Add(fh2SumPtInC60);
+ fOutputList->Add(fh2SumPtOutC10);
+ fOutputList->Add(fh2SumPtOutC10b);
+ fOutputList->Add(fh2SumPtOutC20);
+ fOutputList->Add(fh2SumPtOutC30);
+ fOutputList->Add(fh2SumPtOutC60);
+
+ fOutputList->Add(fh2SumPtInC10bkg);
+ fOutputList->Add(fh2SumPtInC20bkg);
+ fOutputList->Add(fh2SumPtInC30bkg);
+ fOutputList->Add(fh2SumPtInC60bkg);
+ fOutputList->Add(fh2SumPtOutC10bkg);
+ fOutputList->Add(fh2SumPtOutC20bkg);
+ fOutputList->Add(fh2SumPtOutC30bkg);
+ fOutputList->Add(fh2SumPtOutC60bkg);
+ 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(fh2AngStructpt1C10);
+ fOutputList->Add(fh2AngStructpt2C10);
+ fOutputList->Add(fh2AngStructpt3C10);
+ fOutputList->Add(fh2AngStructpt4C10);
+ fOutputList->Add(fh2AngStructpt1C20);
+ fOutputList->Add(fh2AngStructpt2C20);
+ fOutputList->Add(fh2AngStructpt3C20);
+ fOutputList->Add(fh2AngStructpt4C20);
+ fOutputList->Add(fh2AngStructpt1C30);
+ fOutputList->Add(fh2AngStructpt2C30);
+ fOutputList->Add(fh2AngStructpt3C30);
+ fOutputList->Add(fh2AngStructpt4C30);
+ fOutputList->Add(fh2AngStructpt1C60);
+ fOutputList->Add(fh2AngStructpt2C60);
+ fOutputList->Add(fh2AngStructpt3C60);
+ fOutputList->Add(fh2AngStructpt4C60);
+
+
+
// =========== Switch on Sumw2 for all histos ===========
for (Int_t i=0; i<fOutputList->GetEntries(); ++i) {
h1->Sumw2();
continue;
}
-
+ THnSparse *hn = dynamic_cast<THnSparse*>(fOutputList->At(i));
+ if (hn){
+ hn->Sumw2();
+ }
}
TH1::AddDirectory(oldStatus);
} else {
fAOD = dynamic_cast<AliAODEvent*>(AODEvent());
}
- if (!fAOD) {
- AliError("AOD not available");
- return;
- }
+
+ if(fNonStdFile.Length()!=0){
+ // case that we have an AOD extension we need can fetch the jets from the extended output
+ AliAODHandler *aodH = dynamic_cast<AliAODHandler*>(AliAnalysisManager::GetAnalysisManager()->GetOutputEventHandler());
+ fAODExtension = (aodH?aodH->GetExtension(fNonStdFile.Data()):0);
+ if(!fAODExtension){
+ if(fDebug>1)Printf("AODExtension found for %s",fNonStdFile.Data());
+ }}
+
+
+
+
// -- event selection --
fHistEvtSelection->Fill(1); // number of events before event selection
// multiplicity due to input tracks
- Int_t nInputTracks = GetNInputTracks();
-
- if (nInputTracks < fNInputTracksMin || (fNInputTracksMax > -1 && nInputTracks > fNInputTracksMax)){
- fHistEvtSelection->Fill(5);
- PostData(1, fOutputList);
- return;
- }
+ //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(!externalBackground&&fBackgroundBranch.Length()){
+ if(fAOD&&!externalBackground&&fBackgroundBranch.Length()){
externalBackground = (AliAODJetEventBackground*)(fAOD->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());;
+ }
+
Float_t rho = 0;
if(externalBackground)rho = externalBackground->GetBackground(0);
// 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())); //
-
-
-
- for (Int_t iJetType = 0; iJetType < 2; iJetType++) {
+ aodJets[1] = dynamic_cast<TClonesArray*>(fAOD->FindListObject(fJetBranchName[1].Data())); }
+ if(fAODExtension && !aodJets[0]){
+ aodJets[0] = dynamic_cast<TClonesArray*>(fAODExtension->GetAOD()->FindListObject(fJetBranchName[0].Data())); //
+ aodJets[1] = dynamic_cast<TClonesArray*>(fAODExtension->GetAOD()->FindListObject(fJetBranchName[1].Data())); }
+
+ TList ParticleList;
+ Int_t nT = GetListOfTracks(&ParticleList);
+ for (Int_t iJetType = 0; iJetType < 2; iJetType++) {
fListJets[iJetType]->Clear();
if (!aodJets[iJetType]) continue;
Double_t areasmall=0;
Double_t distr=0.;
Double_t phismall=0.;
+
+ 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.};
+
+
+
for(Int_t i=0; i<fListJets[0]->GetEntries(); ++i){
AliAODJet* jetbig = (AliAODJet*)(fListJets[0]->At(i));
etabig = jetbig->Eta();
ptbig = jetbig->Pt();
if(ptbig==0) continue;
areabig = jetbig->EffectiveAreaCharged();
+ Double_t ptcorr=ptbig-rho*areabig;
+ if(ptcorr<=0) continue;
if((etabig<fJetEtaMin)||(etabig>fJetEtaMax)) continue;
Double_t dismin=100.;
Double_t ptmax=-10.;
Int_t index1=-1;
Int_t index2=-1;
- Double_t fracin=0.;
+ 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();
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(ptbig-rho*areabig,fracin/ptbig);
- if((centValue>30)&&(centValue<60)) fh2JetCoreMethod3C30->Fill(ptbig-rho*areabig,fracin/ptbig);
- if(centValue>60) fh2JetCoreMethod3C60->Fill(ptbig-rho*areabig,fracin/ptbig);
- //cout<<"method3 "<<fracin/ptbig<<" "<<ptbig-rho*areabig<<endl;
- /////////////////////////////////////////////////////////////
+ if(distr<=fRadioFrac){ fracin=fracin+genTrack->Pt();}}
+
+ if(centValue<10) fh2JetCoreMethod3C10->Fill(ptcorr,fracin/ptbig);
+
+ if((centValue>20)&&(centValue<40)) fh2JetCoreMethod3C20->Fill(ptcorr,fracin/ptbig);
+
+ if((centValue>30)&&(centValue<60)) fh2JetCoreMethod3C30->Fill(ptcorr,fracin/ptbig);
+
+ if(centValue>60) fh2JetCoreMethod3C60->Fill(ptcorr,fracin/ptbig);
+
+
+
for(Int_t j=0; j<fListJets[1]->GetEntries(); ++j){
AliAODJet* jetsmall = (AliAODJet*)(fListJets[1]->At(j));
phismall = jetsmall->Phi();
ptsmall = jetsmall->Pt();
areasmall = jetsmall->EffectiveAreaCharged();
-
- if((etasmall<fJetEtaMin)||(etasmall>fJetEtaMax)) continue;
-
- Float_t tmpDeltaR=(phismall-phibig)*(phismall-phibig)+(etasmall-etabig)*(etasmall-etabig);
- tmpDeltaR=TMath::Sqrt(tmpDeltaR);
-
- if((ptsmall>ptmax)&&(tmpDeltaR<=fRadioFrac)){ptmax=ptsmall;
+ 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;}}
-
+ if(tmpDeltaR<=dismin){ dismin=tmpDeltaR;
+ index1=j;}} //en of loop over R=0.2 jets
//method1:most concentric jet=core
if(dismin<fMinDist){ AliAODJet* jetmethod1 = (AliAODJet*)(fListJets[1]->At(index1));
-
- if(centValue<10) fh2JetCoreMethod1C10->Fill(ptbig-rho*areabig,jetmethod1->Pt()/ptbig);
- if((centValue>30)&&(centValue<60)) fh2JetCoreMethod1C30->Fill(ptbig-rho*areabig,jetmethod1->Pt()/ptbig);
- if(centValue>60) fh2JetCoreMethod1C60->Fill(ptbig-rho*areabig,jetmethod1->Pt()/ptbig);
-
- //cout<<"method1 "<<jetmethod1->Pt()/ptbig<<" "<<ptbig<<endl;
- }
+ if(centValue<10) fh2JetCoreMethod1C10->Fill(ptcorr,jetmethod1->Pt()/ptbig);
+ if((centValue>20)&&(centValue<40)) fh2JetCoreMethod1C20->Fill(ptcorr,jetmethod1->Pt()/ptbig);
+ if((centValue>30)&&(centValue<60)) fh2JetCoreMethod1C30->Fill(ptcorr,jetmethod1->Pt()/ptbig);
+ if(centValue>60) fh2JetCoreMethod1C60->Fill(ptcorr,jetmethod1->Pt()/ptbig); }
//method2:hardest contained jet=core
if(index2!=-1){
AliAODJet* jetmethod2 = (AliAODJet*)(fListJets[1]->At(index2));
- if(centValue<10) fh2JetCoreMethod2C10->Fill(ptbig-rho*areabig,jetmethod2->Pt()/ptbig);
- if((centValue>30)&&(centValue<60)) fh2JetCoreMethod2C30->Fill(ptbig-rho*areabig,jetmethod2->Pt()/ptbig);
- if(centValue>60) fh2JetCoreMethod2C60->Fill(ptbig-rho*areabig,jetmethod2->Pt()/ptbig);
+ 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){
+
+ 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((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((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(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);}}
+
+ } //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);
+ fh2SumPtOutC10b->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((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(centValue>60.){
+ fh2SumPtInC60bkg->Fill(ptcorr,coronain/ptbig);
+ fh2SumPtOutC60bkg->Fill(ptcorr,coronaout/ptbig);
+ fh2SumPtInC60->Fill(ptcorr,sumPtIn/ptbig);
+ fh2SumPtOutC60->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;}}}}
+ }
+
+ //only jet constituents
+ if(fAngStructCloseTracks==0){
+
+ 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=it+1; 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;}}}}}
- //cout<<"method2 "<<jetmethod2->Pt()/ptbig<<" "<<ptbig<<endl;
- }
-
-
+
+
+
+
+
+
+
+
+
+
+
+
}
-
-
+
+
+ //end loop over R=0.4 jets
+
+ 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);
return;
}
+
+
+
+
+
+
+
+
+
+
+Int_t AliAnalysisTaskJetCore::GetListOfTracks(TList *list){
+
+ Int_t iCount = 0;
+
+ Int_t fFilterMask=272;
+ for(int it = 0;it < fAOD->GetNumberOfTracks();++it){
+ AliAODTrack *tr = fAOD->GetTrack(it);
+ 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++;
+ }
+
+ list->Sort();
+ return iCount;
+
+}
+
+ Int_t AliAnalysisTaskJetCore::GetListOfTracksCloseToJet(TList *list,AliAODJet *jetbig){
+
+ Int_t iCount = 0;
+
+ Int_t fFilterMask=272;
+ for(int it = 0;it < fAOD->GetNumberOfTracks();++it){
+ AliAODTrack *tr = fAOD->GetTrack(it);
+ if((fFilterMask>0)&&!(tr->TestFilterBit(fFilterMask)))continue;
+ if(TMath::Abs(tr->Eta())>0.9)continue;
+ if(tr->Pt()<0.15)continue;
+ Double_t disR=jetbig->DeltaR(tr);
+ if(disR>0.8) continue;
+ list->Add(tr);
+ //cout<<fAOD->GetNumberOfTracks()<<" "<<tr->Pt()<<endl;
+ iCount++;
+ }
+
+ list->Sort();
+ return iCount;
+
+}
+
+
+
+
+
+
+
+
+
+
+
Int_t AliAnalysisTaskJetCore::GetNInputTracks()
{
+Double_t AliAnalysisTaskJetCore::RelativePhi(Double_t mphi,Double_t vphi){
+
+ if (vphi < -1*TMath::Pi()) vphi += (2*TMath::Pi());
+ else if (vphi > TMath::Pi()) vphi -= (2*TMath::Pi());
+ if (mphi < -1*TMath::Pi()) mphi += (2*TMath::Pi());
+ else if (mphi > TMath::Pi()) mphi -= (2*TMath::Pi());
+ double dphi = mphi-vphi;
+ if (dphi < -1*TMath::Pi()) dphi += (2*TMath::Pi());
+ else if (dphi > TMath::Pi()) dphi -= (2*TMath::Pi());
+ return dphi;//dphi in [-Pi, Pi]
+}
+