-#ifndef ALIANALYSISTASKMINIJET_H\r
-#define ALIANALYSISTASKMINIJET_H\r
-\r
-// Two-particle correlations using all particles over pt threshold\r
-// Extract mini-jet yield and fragmentation properties via Delta-Phi histograms\r
-// Can use ESD or AOD, reconstructed and Monte Carlo data as input\r
-// Author: eva.sicking@cern.ch\r
-\r
-class TList;\r
-class TH1F;\r
-class TH2F;\r
-class TProfile;\r
-class THnSparse;\r
-class AliESDtrackCuts;\r
-\r
-#include "AliAnalysisTaskSE.h"\r
-#include <vector>\r
-\r
-class AliAnalysisTaskMinijet : public AliAnalysisTaskSE {\r
- public:\r
- AliAnalysisTaskMinijet(const char *name="<default name>");\r
- virtual ~AliAnalysisTaskMinijet();\r
- \r
- virtual void UserCreateOutputObjects();\r
- virtual void UserExec(Option_t* option);\r
- virtual void Terminate(Option_t *);\r
- virtual void SetCuts(AliESDtrackCuts* cuts){fCuts = cuts;}\r
- \r
- void SetUseMC(Bool_t useMC=kTRUE, Bool_t mcOnly=kFALSE) {fUseMC = useMC; fMcOnly=mcOnly;}\r
- void SetAnalyseOnlyPrimaries(Bool_t analysePrimOnly) {fAnalysePrimOnly = analysePrimOnly;} // not used anymore\r
- void SetPtRange(Float_t ptMin, Float_t ptMax) {fPtMin = ptMin; fPtMax = ptMax; }\r
- void SetTriggerPtCut(Float_t triggerPtCut) {fTriggerPtCut = triggerPtCut;} \r
- void SetAssociatePtCut(Float_t associatePtCut) {fAssociatePtCut = associatePtCut;} \r
- void SetModeEsdAod(Int_t mode) {fMode = mode;}\r
- void SetTriggerMask(Int_t triggerType) {fTriggerType = triggerType;}\r
- void SetFilterBit(Int_t filterBit) {fFilterBit = filterBit;}\r
- void SetMaxVertexZ(Float_t vertexZCut) {fVertexZCut = vertexZCut;}\r
- void SetMaxEta(Float_t etaCut) {fEtaCut = etaCut;}\r
- void SetMaxEtaSeed(Float_t etaCutSeed) {fEtaCutSeed = etaCutSeed;}\r
- void SetSelectParticles(Int_t selectParticles) {fSelectParticles = selectParticles;}\r
- void SetSelectParticlesAssoc(Int_t selectParticlesAssoc) {fSelectParticlesAssoc = selectParticlesAssoc;}\r
- void SetCheckSDD(Bool_t checkSDD, Int_t selOption) {fCheckSDD = checkSDD; fSelOption = selOption;}\r
- void SetCorrStrangeness(Bool_t corrStrangeness) {fCorrStrangeness = corrStrangeness;}\r
- void SetThreeParticleCorrelation(Bool_t threeParticleCorr) {fThreeParticleCorr = threeParticleCorr;}\r
- void SetRejectCorrupted(Bool_t rejectChunks, Int_t nTPC) {fRejectChunks = rejectChunks; fNTPC = nTPC;}\r
-\r
- void SetCentralityMethod(TString centralityMethod) {fCentralityMethod = centralityMethod;}\r
- \r
- \r
- private:\r
-\r
- Double_t ReadEventESD (std::vector<Float_t> &pt, std::vector<Float_t> &eta,\r
- std::vector<Float_t> &phi, std::vector<Short_t> &charge,\r
- std::vector<Float_t> &strangnessWeight,\r
- std::vector<Double_t> &nTracksTracklets, const Int_t step);\r
- Double_t ReadEventESDRecMcProp(std::vector<Float_t> &pt, std::vector<Float_t> &eta,\r
- std::vector<Float_t> &phi, std::vector<Short_t> &charge,\r
- std::vector<Float_t> &strangnessWeight,\r
- std::vector<Double_t> &nTracksTracklets, const Int_t step);\r
- Double_t ReadEventESDMC (std::vector<Float_t> &pt, std::vector<Float_t> &eta,\r
- std::vector<Float_t> &phi, std::vector<Short_t> &charge,\r
- std::vector<Float_t> &strangnessWeight,\r
- std::vector<Double_t> &nTracksTracklets, const Int_t step);\r
- \r
- Double_t ReadEventAOD (std::vector<Float_t> &pt, std::vector<Float_t> &eta,\r
- std::vector<Float_t> &phi, std::vector<Short_t> &charge,\r
- std::vector<Float_t> &strangnessWeight,\r
- std::vector<Double_t> &nTracksTracklets, const Int_t step);\r
- Double_t ReadEventAODRecMcProp(std::vector<Float_t> &pt, std::vector<Float_t> &eta,\r
- std::vector<Float_t> &phi, std::vector<Short_t> &charge,\r
- std::vector<Float_t> &strangnessWeight,\r
- std::vector<Double_t> &nTracksTracklets, const Int_t step);\r
- Double_t ReadEventAODMC (std::vector<Float_t> &pt, std::vector<Float_t> &eta,\r
- std::vector<Float_t> &phi, std::vector<Short_t> &charge,\r
- std::vector<Float_t> &strangnessWeight,\r
- std::vector<Double_t> &nTracksTracklets, const Int_t step);\r
- \r
- void Analyse (const std::vector<Float_t> &pt, \r
- const std::vector<Float_t> &eta, \r
- const std::vector<Float_t> &phi, \r
- const std::vector<Short_t> &charge, \r
- const std::vector<Float_t> &strangnessWeight,\r
- const Double_t ntacks, const Int_t ntacklets=0,\r
- const Int_t nAll=0, const Int_t step=0);\r
- \r
- Bool_t SelectParticlePlusCharged(const Short_t charge, const Int_t pdg, const Bool_t prim);\r
- Bool_t SelectParticle(const Short_t charge, const Int_t pdg, const Bool_t prim);\r
- Bool_t CheckEvent(const Bool_t recVertex);\r
- const Double_t* CreateLogAxis(const Int_t nbins, const Double_t xmin, const Double_t xmax); \r
- Bool_t CheckLikeSign(const Short_t chargeEventAxis, const Short_t chargeOthers);\r
-\r
-\r
- Bool_t fUseMC; // flag for Monte Carlo usages\r
- Bool_t fMcOnly; // flag defines, if only MC data is used in analysis or also reconstructed data\r
- Double_t fBSign; // magnetic field\r
- Bool_t fAnalysePrimOnly; // flag for analysis of primaries only (also in reconstructed data)\r
- Float_t fPtMin; // set lower limit for pt acceptance for mutliplicity defintion\r
- Float_t fPtMax; // set upper limit for pt acceptance for mutliplicity defintion\r
- AliESDtrackCuts* fCuts; // List of cuts for ESDs\r
- Float_t fTriggerPtCut; // cut on particle pt used as event axis\r
- Float_t fAssociatePtCut; // cut on particle pt used for correlations\r
- Int_t fMode; // ESD(=0) of AOD(=1) reading \r
- Int_t fTriggerType; // sets trigger -> AliVEvent::kMB, AliVEvent::kHighMult\r
- Int_t fFilterBit; // Filter bit written in ESD filter, select track type\r
- Float_t fVertexZCut; // vertex cut\r
- Float_t fEtaCut; // eta acceptance cut\r
- Float_t fEtaCutSeed; // eta acceptance cut for seed\r
- Int_t fSelectParticles; // only in cas of MC: use also neutral particles or not \r
- Int_t fSelectParticlesAssoc; // only in cas of MC: use also neutral particles or not \r
- Bool_t fCheckSDD; // check if SDD was in read out partition (needed for LHC11a)\r
- Int_t fSelOption; // 0 = use hit in SDD for event selection, 1 = use trigger for event selection\r
- Bool_t fCorrStrangeness; // for data correction -> Pythia simulations underestimate contamination from strangness\r
- Bool_t fThreeParticleCorr; // perform three particle correlation\r
- Bool_t fRejectChunks; // rejection of chunks in which no ITS tracks are reconstructed\r
- Int_t fNTPC; // track number limit for rejection decision.\r
-\r
- AliESDEvent *fESDEvent; //! esd event\r
- AliAODEvent *fAODEvent; //! aod event\r
- Double_t fNMcPrimAccept; // global variable for mc multiplucity\r
- Double_t fNRecAccept; // global variable for rec multiplucity\r
- Float_t fNRecAcceptStrangeCorr; // global variable for rec multiplucity\r
- Double_t fNMcPrimAcceptTracklet; // global variable for mc multiplucity\r
- Double_t fNRecAcceptTracklet; // global variable for rec multiplucity\r
- Float_t fVzEvent; // global variable for rec vertex position\r
- Double_t fMeanPtRec; // global variable for rec mean pt\r
- Double_t fLeadingPtRec; // global variable for rec mean pt\r
-\r
- TList *fHists; // output list\r
- TH1F *fStep; // how many events have passed which correction step\r
- TH1F *fEventStat; // how many events are accepted by trigger, vertex selection, 1 track in acceptance (for real data)\r
- TH1F *fHistPt; // Pt spectrum ESD\r
- TH1F *fHistPtMC; // Pt spectrum MC\r
- TH2F *fNContrNtracklets; // control histogram for vertex->nContributers and number of tracklets\r
- TH2F *fNContrNtracks; // control histogram for vertex->nContributers and number of tracks\r
- TH2F *fCorruptedChunks; // control histogram: TPC tracks versus ITS-TPC-tracks\r
- TH2F *fCorruptedChunksAfter; // control histogram: TPC tracks versus ITS-TPC-tracks\r
-\r
- TH2F *fNmcNch; // N mc - N ch rec\r
- TProfile *fPNmcNch; // N mc - N ch rec\r
- TH2F *fNmcNchVtx; // N mc - N ch rec for events with reconstructed vertex\r
- TH2F *fNmcNchVtxStrangeCorr; // N mc - N ch rec for events with reconstructed vertex + strangeness correction\r
- TProfile *fPNmcNchVtx; // N mc - N ch rec for events with reconstructed vertex\r
- TH2F *fNmcNchTracklet; // N mc - N ch rec\r
- TProfile *fPNmcNchTracklet; // N mc - N ch rec\r
- TH2F *fNmcNchVtxTracklet; // N mc - N ch rec for events with reconstructed vertex\r
- TProfile *fPNmcNchVtxTracklet; // N mc - N ch rec for events with reconstructed vertex\r
- TH2F *fChargedPi0; // charged versus charged+Pi0\r
- TH1F *fVertexCheck; // check which fraction of events has vtx_rec but no good vtx_mc\r
- TH1F *fPropagateDca; // check of AliAODtrack::PropagateToDca\r
-\r
- THnSparse *fMapSingleTrig[8]; //! multi-dim histo for trigger track properties\r
- THnSparse *fMapPair[8]; //! multi-dim histo for pair properties\r
- THnSparse *fMapEvent[8]; //! multi-dim histo for event properties\r
- THnSparse *fMapAll[8]; //! multi-dim histo for properties of all analysed tracks\r
- THnSparse *fMapThree[8]; //! multi-dim histo for properties of three particle correlations\r
- \r
- TH1F * fVertexZ[8]; // z of vertex\r
- TH1F * fNcharge[8]; // pt\r
- TH1F * fPt[8]; // pt\r
- TH1F * fEta[8]; // eta\r
- TH1F * fPhi[8]; // phi\r
- TH1F * fDcaXY[8]; // dca xy direction\r
- TH1F * fDcaZ[8]; // dca z direction\r
-\r
- TH1F * fPtSeed[8]; // pt of seed (event axis)\r
- TH1F * fEtaSeed[8]; // eta of seed \r
- TH1F * fPhiSeed[8]; // phi of seed\r
-\r
- TH1F * fPtOthers[8]; // pt of all other particels used in dEtadPhi\r
- TH1F * fEtaOthers[8]; // eta of all other particels used in dEtadPhi\r
- TH1F * fPhiOthers[8]; // phi of all other particels used in dEtadPhi\r
- TH2F * fPtEtaOthers[8]; // pt-eta of all other particels used in dEtadPhi\r
-\r
-\r
- TH2F * fPhiEta[8]; // eta - phi\r
- TH2F * fDPhiDEtaEventAxis[8]; // correlation dEta-dPhi towards event axis\r
- TH2F * fDPhiDEtaEventAxisSeeds[8]; // correlation dEta-dPhi towards event axis of trigger particles\r
- TH1F * fTriggerNch[8]; // number of triggers with accepted-track number\r
- TH2F * fTriggerNchSeeds[8]; // number of triggers with accepted-track number\r
- TH1F * fTriggerTracklet[8]; // number of triggers with accepted-tracklet number\r
- TH2F * fNch07Nch[8]; // nCharged with pT>fTriggerPtCut vs nCharged\r
- TProfile * fPNch07Nch[8]; // nCharged with pT>fTriggerPtCut vs nCharged\r
- \r
- TH2F * fNch07Tracklet[8]; // nCharged with pT>fTriggerPtCut vs nTracklet\r
- TH2F * fNchTracklet[8]; // nCharged vs nTracklet\r
- TProfile * fPNch07Tracklet[8]; // nCharged with pT>fTriggerPtCut vs nTracklet\r
-\r
- TH1F * fDPhiEventAxis[8]; // delta phi of associate tracks to event axis\r
- TH2F * fDPhi1DPhi2[8]; // dPhi1 versus dPhi2: three particle correlation test\r
- \r
- TString fCentralityMethod; //centrality pA\r
- \r
- AliAnalysisTaskMinijet(const AliAnalysisTaskMinijet&); // not implemented\r
- AliAnalysisTaskMinijet& operator=(const AliAnalysisTaskMinijet&); // not implemented\r
- \r
- ClassDef(AliAnalysisTaskMinijet, 2); // mini jet analysis with two particle correlations\r
-};\r
-\r
-#endif\r
+#ifndef ALIANALYSISTASKMINIJET_H
+#define ALIANALYSISTASKMINIJET_H
+
+// Two-particle correlations using all particles over pt threshold
+// Extract mini-jet yield and fragmentation properties via Delta-Phi histograms
+// Can use ESD or AOD, reconstructed and Monte Carlo data as input
+// Author: eva.sicking@cern.ch
+
+class TList;
+class TH1F;
+class TH2F;
+class TProfile;
+class THnSparse;
+class AliESDtrackCuts;
+
+#include "AliAnalysisTaskSE.h"
+#include <vector>
+
+class AliAnalysisTaskMinijet : public AliAnalysisTaskSE {
+ public:
+ AliAnalysisTaskMinijet(const char *name="<default name>");
+ virtual ~AliAnalysisTaskMinijet();
+
+ virtual void UserCreateOutputObjects();
+ virtual void UserExec(Option_t* option);
+ virtual void Terminate(Option_t *);
+ virtual void SetCuts(AliESDtrackCuts* cuts){fCuts = cuts;}
+
+ void SetUseMC(Bool_t useMC=kTRUE, Bool_t mcOnly=kFALSE) {fUseMC = useMC; fMcOnly=mcOnly;}
+ void SetAnalyseOnlyPrimaries(Bool_t analysePrimOnly) {fAnalysePrimOnly = analysePrimOnly;} // not used anymore
+ void SetPtRange(Float_t ptMin, Float_t ptMax) {fPtMin = ptMin; fPtMax = ptMax; }
+ void SetTriggerPtCut(Float_t triggerPtCut) {fTriggerPtCut = triggerPtCut;}
+ void SetAssociatePtCut(Float_t associatePtCut) {fAssociatePtCut = associatePtCut;}
+ void SetModeEsdAod(Int_t mode) {fMode = mode;}
+ void SetTriggerMask(Int_t triggerType) {fTriggerType = triggerType;}
+ void SetFilterBit(Int_t filterBit) {fFilterBit = filterBit;}
+ void SetMaxVertexZ(Float_t vertexZCut) {fVertexZCut = vertexZCut;}
+ void SetMaxEta(Float_t etaCut) {fEtaCut = etaCut;}
+ void SetMaxEtaSeed(Float_t etaCutSeed) {fEtaCutSeed = etaCutSeed;}
+ void SetSelectParticles(Int_t selectParticles) {fSelectParticles = selectParticles;}
+ void SetSelectParticlesAssoc(Int_t selectParticlesAssoc) {fSelectParticlesAssoc = selectParticlesAssoc;}
+ void SetCheckSDD(Bool_t checkSDD, Int_t selOption) {fCheckSDD = checkSDD; fSelOption = selOption;}
+ void SetCorrStrangeness(Bool_t corrStrangeness) {fCorrStrangeness = corrStrangeness;}
+ void SetThreeParticleCorrelation(Bool_t threeParticleCorr) {fThreeParticleCorr = threeParticleCorr;}
+ void SetRejectCorrupted(Bool_t rejectChunks, Int_t nTPC) {fRejectChunks = rejectChunks; fNTPC = nTPC;}
+
+ void SetCentralityMethod(TString centralityMethod) {fCentralityMethod = centralityMethod;}
+
+
+ private:
+
+ Double_t ReadEventESD (std::vector<Float_t> &pt, std::vector<Float_t> &eta,
+ std::vector<Float_t> &phi, std::vector<Short_t> &charge,
+ std::vector<Float_t> &strangnessWeight,
+ std::vector<Double_t> &nTracksTracklets, const Int_t step);
+ Double_t ReadEventESDRecMcProp(std::vector<Float_t> &pt, std::vector<Float_t> &eta,
+ std::vector<Float_t> &phi, std::vector<Short_t> &charge,
+ std::vector<Float_t> &strangnessWeight,
+ std::vector<Double_t> &nTracksTracklets, const Int_t step);
+ Double_t ReadEventESDMC (std::vector<Float_t> &pt, std::vector<Float_t> &eta,
+ std::vector<Float_t> &phi, std::vector<Short_t> &charge,
+ std::vector<Float_t> &strangnessWeight,
+ std::vector<Double_t> &nTracksTracklets, const Int_t step);
+
+ Double_t ReadEventAOD (std::vector<Float_t> &pt, std::vector<Float_t> &eta,
+ std::vector<Float_t> &phi, std::vector<Short_t> &charge,
+ std::vector<Float_t> &strangnessWeight,
+ std::vector<Double_t> &nTracksTracklets, const Int_t step);
+ Double_t ReadEventAODRecMcProp(std::vector<Float_t> &pt, std::vector<Float_t> &eta,
+ std::vector<Float_t> &phi, std::vector<Short_t> &charge,
+ std::vector<Float_t> &strangnessWeight,
+ std::vector<Double_t> &nTracksTracklets, const Int_t step);
+ Double_t ReadEventAODMC (std::vector<Float_t> &pt, std::vector<Float_t> &eta,
+ std::vector<Float_t> &phi, std::vector<Short_t> &charge,
+ std::vector<Float_t> &strangnessWeight,
+ std::vector<Double_t> &nTracksTracklets, const Int_t step);
+
+ void Analyse (const std::vector<Float_t> &pt,
+ const std::vector<Float_t> &eta,
+ const std::vector<Float_t> &phi,
+ const std::vector<Short_t> &charge,
+ const std::vector<Float_t> &strangnessWeight,
+ const Double_t ntacks, const Int_t ntacklets=0,
+ const Int_t nAll=0, const Int_t step=0);
+
+ Bool_t SelectParticlePlusCharged(const Short_t charge, const Int_t pdg, const Bool_t prim);
+ Bool_t SelectParticle(const Short_t charge, const Int_t pdg, const Bool_t prim);
+ Bool_t CheckEvent(const Bool_t recVertex);
+ const Double_t* CreateLogAxis(const Int_t nbins, const Double_t xmin, const Double_t xmax);
+ Bool_t CheckLikeSign(const Short_t chargeEventAxis, const Short_t chargeOthers);
+
+
+ Bool_t fUseMC; // flag for Monte Carlo usages
+ Bool_t fMcOnly; // flag defines, if only MC data is used in analysis or also reconstructed data
+ Double_t fBSign; // magnetic field
+ Bool_t fAnalysePrimOnly; // flag for analysis of primaries only (also in reconstructed data)
+ Float_t fPtMin; // set lower limit for pt acceptance for mutliplicity defintion
+ Float_t fPtMax; // set upper limit for pt acceptance for mutliplicity defintion
+ AliESDtrackCuts* fCuts; // List of cuts for ESDs
+ Float_t fTriggerPtCut; // cut on particle pt used as event axis
+ Float_t fAssociatePtCut; // cut on particle pt used for correlations
+ Int_t fMode; // ESD(=0) of AOD(=1) reading
+ Int_t fTriggerType; // sets trigger -> AliVEvent::kMB, AliVEvent::kHighMult
+ Int_t fFilterBit; // Filter bit written in ESD filter, select track type
+ Float_t fVertexZCut; // vertex cut
+ Float_t fEtaCut; // eta acceptance cut
+ Float_t fEtaCutSeed; // eta acceptance cut for seed
+ Int_t fSelectParticles; // only in cas of MC: use also neutral particles or not
+ Int_t fSelectParticlesAssoc; // only in cas of MC: use also neutral particles or not
+ Bool_t fCheckSDD; // check if SDD was in read out partition (needed for LHC11a)
+ Int_t fSelOption; // 0 = use hit in SDD for event selection, 1 = use trigger for event selection
+ Bool_t fCorrStrangeness; // for data correction -> Pythia simulations underestimate contamination from strangness
+ Bool_t fThreeParticleCorr; // perform three particle correlation
+ Bool_t fRejectChunks; // rejection of chunks in which no ITS tracks are reconstructed
+ Int_t fNTPC; // track number limit for rejection decision.
+
+ AliESDEvent *fESDEvent; //! esd event
+ AliAODEvent *fAODEvent; //! aod event
+ Double_t fNMcPrimAccept; // global variable for mc multiplucity
+ Double_t fNRecAccept; // global variable for rec multiplucity
+ Float_t fNRecAcceptStrangeCorr; // global variable for rec multiplucity
+ Double_t fNMcPrimAcceptTracklet; // global variable for mc multiplucity
+ Double_t fNRecAcceptTracklet; // global variable for rec multiplucity
+ Float_t fVzEvent; // global variable for rec vertex position
+ Double_t fMeanPtRec; // global variable for rec mean pt
+ Double_t fLeadingPtRec; // global variable for rec mean pt
+
+ TList *fHists; // output list
+ TH1F *fStep; // how many events have passed which correction step
+ TH1F *fEventStat; // how many events are accepted by trigger, vertex selection, 1 track in acceptance (for real data)
+ TH1F *fHistPt; // Pt spectrum ESD
+ TH1F *fHistPtMC; // Pt spectrum MC
+ TH2F *fNContrNtracklets; // control histogram for vertex->nContributers and number of tracklets
+ TH2F *fNContrNtracks; // control histogram for vertex->nContributers and number of tracks
+ TH2F *fCorruptedChunks; // control histogram: TPC tracks versus ITS-TPC-tracks
+ TH2F *fCorruptedChunksAfter; // control histogram: TPC tracks versus ITS-TPC-tracks
+
+ TH2F *fNmcNch; // N mc - N ch rec
+ TProfile *fPNmcNch; // N mc - N ch rec
+ TH2F *fNmcNchVtx; // N mc - N ch rec for events with reconstructed vertex
+ TH2F *fNmcNchVtxStrangeCorr; // N mc - N ch rec for events with reconstructed vertex + strangeness correction
+ TProfile *fPNmcNchVtx; // N mc - N ch rec for events with reconstructed vertex
+ TH2F *fNmcNchTracklet; // N mc - N ch rec
+ TProfile *fPNmcNchTracklet; // N mc - N ch rec
+ TH2F *fNmcNchVtxTracklet; // N mc - N ch rec for events with reconstructed vertex
+ TProfile *fPNmcNchVtxTracklet; // N mc - N ch rec for events with reconstructed vertex
+ TH2F *fChargedPi0; // charged versus charged+Pi0
+ TH1F *fVertexCheck; // check which fraction of events has vtx_rec but no good vtx_mc
+ TH1F *fPropagateDca; // check of AliAODtrack::PropagateToDca
+
+ THnSparse *fMapSingleTrig[8]; //! multi-dim histo for trigger track properties
+ THnSparse *fMapPair[8]; //! multi-dim histo for pair properties
+ THnSparse *fMapEvent[8]; //! multi-dim histo for event properties
+ THnSparse *fMapAll[8]; //! multi-dim histo for properties of all analysed tracks
+ THnSparse *fMapThree[8]; //! multi-dim histo for properties of three particle correlations
+
+ TH1F * fVertexZ[8]; // z of vertex
+ TH1F * fNcharge[8]; // pt
+ TH1F * fPt[8]; // pt
+ TH1F * fEta[8]; // eta
+ TH1F * fPhi[8]; // phi
+ TH1F * fDcaXY[8]; // dca xy direction
+ TH1F * fDcaZ[8]; // dca z direction
+
+ TH1F * fPtSeed[8]; // pt of seed (event axis)
+ TH1F * fEtaSeed[8]; // eta of seed
+ TH1F * fPhiSeed[8]; // phi of seed
+
+ TH1F * fPtOthers[8]; // pt of all other particels used in dEtadPhi
+ TH1F * fEtaOthers[8]; // eta of all other particels used in dEtadPhi
+ TH1F * fPhiOthers[8]; // phi of all other particels used in dEtadPhi
+ TH2F * fPtEtaOthers[8]; // pt-eta of all other particels used in dEtadPhi
+
+
+ TH2F * fPhiEta[8]; // eta - phi
+ TH2F * fDPhiDEtaEventAxis[8]; // correlation dEta-dPhi towards event axis
+ TH2F * fDPhiDEtaEventAxisSeeds[8]; // correlation dEta-dPhi towards event axis of trigger particles
+ TH1F * fTriggerNch[8]; // number of triggers with accepted-track number
+ TH2F * fTriggerNchSeeds[8]; // number of triggers with accepted-track number
+ TH1F * fTriggerTracklet[8]; // number of triggers with accepted-tracklet number
+ TH2F * fNch07Nch[8]; // nCharged with pT>fTriggerPtCut vs nCharged
+ TProfile * fPNch07Nch[8]; // nCharged with pT>fTriggerPtCut vs nCharged
+
+ TH2F * fNch07Tracklet[8]; // nCharged with pT>fTriggerPtCut vs nTracklet
+ TH2F * fNchTracklet[8]; // nCharged vs nTracklet
+ TProfile * fPNch07Tracklet[8]; // nCharged with pT>fTriggerPtCut vs nTracklet
+
+ TH1F * fDPhiEventAxis[8]; // delta phi of associate tracks to event axis
+ TH2F * fDPhi1DPhi2[8]; // dPhi1 versus dPhi2: three particle correlation test
+
+ TString fCentralityMethod; //centrality pA
+
+ AliAnalysisTaskMinijet(const AliAnalysisTaskMinijet&); // not implemented
+ AliAnalysisTaskMinijet& operator=(const AliAnalysisTaskMinijet&); // not implemented
+
+ ClassDef(AliAnalysisTaskMinijet, 2); // mini jet analysis with two particle correlations
+};
+
+#endif