AliPhysics  251aa1e (251aa1e)
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
AliAnalysisTaskPhiFlow.h
Go to the documentation of this file.
1 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. */
2 /* See cxx source for full Copyright notice */
3 /* $Id$ */
4 
5 // AliAnalysisTaskPhiFlow:
6 // origin: Redmer Alexander Bertens (rbertens@nikhef.nl)
7 // analyis task for phi-meson reconstruction and estimation of v_n
8 
9 #ifndef ALIANALYSISTASKPHIFLOW_H
10 #define ALIANALYSISTASKPHIFLOW_H
11 
12 class TH1F;
13 class TH2F;
14 class TProfile;
15 class AliESDEvent;
16 class AliESDtrackCuts;
17 class AliFlowTrackCuts;
18 class AliFlowEvent;
20 class AliFlowBayesianPID;
21 class AliEventPoolManager;
22 class AliPIDCombined;
23 class AliVVertex;
24 
25 #include "AliAnalysisTaskSE.h"
26 
28 {
29  public:
30  AliPhiMesonHelperTrack(Float_t eta, Float_t phi, Float_t p, Float_t px, Float_t py, Float_t pz, Float_t pt, Int_t charge) : fEta(eta), fPhi(phi), fp(p), fpX(px), fpY(py), fpZ(pz), fpT(pt), fCharge(charge) { }
32  virtual Double_t P() const { return fp; }
33  virtual Double_t Px() const { return fpX; }
34  virtual Double_t Py() const { return fpY; }
35  virtual Double_t Pz() const { return fpZ; }
36  virtual Double_t Pt() const { return fpT; }
37  virtual Double_t Phi() const { return fPhi; }
38  virtual Double_t Eta() const { return fEta; }
39  virtual Int_t Charge() const { return fCharge; }
40  void InitializeHelperTrack(Float_t eta, Float_t phi, Float_t p, Float_t px, Float_t py, Float_t pz, Float_t pt, Int_t charge) { fEta = eta; fPhi = phi; fp = p; fpX = px; fpY = py; fpZ = pz; fpT = pt; fCharge = charge; }
41  private:
42  Float_t fEta; // eta
43  Float_t fPhi; // phi
44  Float_t fp; // p
45  Float_t fpX; // pX
46  Float_t fpY; // pY
47  Float_t fpZ; // pZ
48  Float_t fpT; // pT
49  Int_t fCharge; // charge
50  ClassDef(AliPhiMesonHelperTrack, 1); // lightweight helper track for phi reconstruction
51 };
52 
54 {
55  public:
57  AliAnalysisTaskPhiFlow(const char *name);
58  virtual ~AliAnalysisTaskPhiFlow();
59  Int_t SetDebugLevelPhiTask(Int_t debug) {fDebug = debug; return fDebug; }
60  Bool_t SetIsMC(Bool_t ismc) {fIsMC = ismc; return fIsMC; }
61  Bool_t UseEventMixing(Bool_t mix, Bool_t type) {fEventMixing = mix; fTypeMixing = type; return mix; }
62  Bool_t SetVZEROSubEvents(Bool_t v0) { fV0 = v0; return v0; }
63  TH1F* BookHistogram(const char * name);
64  TH2F* BookPIDHistogram(const char * name, Bool_t TPC);
65  TH1F* InitPtSpectraHistograms(Float_t nmin, Float_t nmax);
66  TH1F* BookPtHistogram(const char* name);
68  virtual void UserCreateOutputObjects();
69  AliEventPoolManager* InitializeEventMixing();
70  void SetPtBins(Float_t bin[19], Int_t n) { for(Int_t i = 0; i < n+1; i++) fPtBins[i] = bin[i]; fNPtBins = n; }
71  template <typename T> Double_t InvariantMass(const T* track1, const T* track2) const;
72  // template <typename T> Double_t DeltaDipAngle(const T* track1, const T* track2) const;
73  // template <typename T> Bool_t CheckDeltaDipAngle(const T* track1, const T* track2) const;
74  template <typename T> Bool_t CheckCandidateEtaPtCut(const T* track1, const T* track2) const;
76  void SetCentralityParameters(Double_t min, Double_t max, const char* a, const char* b, Bool_t c, Bool_t d) {
77  fCentralityMin = min;
78  fCentralityMax = max;
82  fCentralityCut2011 = d;}
84  Double_t GetCenMin() const {return fCentralityMin; }
85  Double_t GetCenMax() const {return fCentralityMax; }
86  const char* GetCentralityMethod() const {return fkCentralityMethodA; }
87  void SetVertexZ(Float_t z) { fVertexRange = z; }
88  Float_t GetVertexZ() const { return fVertexRange; }
89  // void SetMaxDeltaDipAngleAndPt(Float_t a, Float_t pt) { fDeltaDipAngle = a;
90  // fDeltaDipPt = pt;
91  // fApplyDeltaDipCut = kTRUE; };
92  // Float_t GetDeltaDipAngle() const {return fDeltaDipAngle; }
93  // Float_t GetDeltaDipPt() const {return fDeltaDipPt; }
94  template <typename T> Bool_t EventCut(T* event);
95  template <typename T> void PlotMultiplcities(const T* event) const;
96  template <typename T> Bool_t CheckVertex(const T* event);
97  template <typename T> Bool_t CheckCentrality(T* event);
99  template <typename T> Bool_t PassesTPCbayesianCut(T* track) const;
100  Bool_t PassesDCACut(AliAODTrack* track) const;
101  Bool_t IsKaon(AliAODTrack* track) const;
102  template <typename T> Double_t PhiPt(const T* track_1, const T* track_2) const;
103  template <typename T> void PtSelector(Int_t _track_type, const T* track_1, const T* track_2) const;
104  template <typename T> Bool_t PhiTrack(T* track) const;
105  template <typename T> void SetNullCuts(T* esd);
106  void PrepareFlowEvent(Int_t iMulti);
107  void VZEROSubEventAnalysis();
108  virtual void UserExec(Option_t *option);
109  virtual void Exec(Option_t *);
110  void ReconstructionWithEventMixing(TObjArray* MixingCandidates) const;
111  virtual void Terminate(Option_t *);
112  void SetPOICuts(AliFlowTrackCuts *cutsPOI) { fPOICuts = cutsPOI; }
113  void SetRPCuts(AliFlowTrackCuts *cutsRP) { fCutsRP = cutsRP; }
115  AliFlowTrackCuts* GetRPCuts() const {return fCutsRP;}
116  void SetPIDConfiguration(Double_t prob[7]) { for(Int_t i = 0; i < 7; i++) fPIDConfig[i] = prob[i]; }
117  void GetPIDConfiguration(Double_t prob[7]) const {for(Int_t i = 0; i < 7; i++) prob[i] = fPIDConfig[i]; }
118  void SetPOIDCAXYZ(Double_t dca[5]) { for(Int_t i = 0; i < 5; i++) fDCAConfig[i] = dca[i]; }
119  void GetPOIDCZXYZ(Double_t dca[5]) const { for(Int_t i = 0; i < 5; i++) dca[i] = fDCAConfig[i]; }
120  void SetMixingBins(Int_t c[20], Int_t v[20]) {for(Int_t i = 0; i < 20; i++) { fCentralityMixingBins[i] = c[i];
121  fVertexMixingBins[i] = v[i]; } }
122  void SetMixingParameters(Int_t p[3]) { for(Int_t i = 0; i < 3; i++) fMixingParameters[i] = p[i]; }
123  void GetMixingParameters(Int_t p[3]) const { for(Int_t i = 0; i < 3; i++) p[i] = fMixingParameters[i]; }
124  void SetCandidateEtaAndPt(Double_t mineta, Double_t maxeta, Double_t minpt, Double_t maxpt) { fCandidateMinEta = mineta;
125  fCandidateMaxEta = maxeta;
126  fCandidateMinPt = minpt;
127  fCandidateMaxPt = maxpt;
128  fCandidateEtaPtCut = kTRUE;}
130  fCandidateMinY = min;
131  fCandidateMaxY = max;
132  fCandidateYCut = kTRUE;
133  }
134  void GetCandidateEtaAndPt(Double_t etapt[4]) const { etapt[0] = fCandidateMinEta;
135  etapt[1] = fCandidateMaxEta;
136  etapt[2] = fCandidateMinPt;
137  etapt[3] = fCandidateMaxPt; }
139  fMinMass = minMass;
140  fMaxMass= maxMass; }
141  void IsMC();
142  Bool_t SetQA(Bool_t qa) {fQA = qa; return fQA;}
145  fUsePidResponse = kTRUE; // bayesian pid object will require some event info
146  fCentrality = 5.;} // should be set by user, skipping event selection will also forego centrality selection
148  Double_t GetWDist(const AliVVertex* v0, const AliVVertex* v1);
149  Bool_t plpMV(const AliAODEvent* aod);
151  Short_t FindMinNSigma(Double_t nSpi, Double_t nSk, Double_t nSp) const;
152  Bool_t GetDoubleCountingK(Double_t nSk, Short_t minNSigma) const;
153 
154  private:
155 
156  Int_t fDebug; // debug level (0 none, 1 fcn calls, 2 verbose)
157  Bool_t fIsMC; // use mc mode
158  Bool_t fEventMixing; // use event mixing
159  Bool_t fTypeMixing; // select type: kTRUE for unlike sign background, kFALSE for like sign background
160  Bool_t fQA; // make qa plots
161  Bool_t fV0; // use three subevents including vzero
162  Int_t fMassBins; // mass bins
163  Double_t fMinMass; // mass range
164  Double_t fMaxMass; // mass range
165  AliFlowTrackCuts *fCutsRP; // track cuts for reference particles
166  AliFlowTrackCuts *fNullCuts; // dummy cuts for flow event tracks
167  AliPIDResponse *fPIDResponse;
170  TObjArray *fCandidates; // candidate array
171  Bool_t fCandidateEtaPtCut; // set eta and pt cut for candidate tracks and combinatorial background
172  Double_t fCandidateMinEta; // minimum eta for candidates
173  Double_t fCandidateMaxEta; // maximum eta for candidates
174  Double_t fCandidateMinPt; // minimum pt for candidates
175  Double_t fCandidateMaxPt; // maximum pt for candidates
176  Bool_t fCandidateYCut; // y cut on candidates
177  Double_t fCandidateMinY; // min y value for candidate
178  Double_t fCandidateMaxY; // max y value for canddiates
179  Double_t fPIDConfig[7]; // configure pid routine
180  Double_t fDCAConfig[5]; // configure dca routine
181  Int_t fMixingParameters[3]; // mixing: poolsize, mixing tracks, pool buffer
182  Int_t fCentralityMixingBins[20]; // configure centrality bins for event mixing
183  Int_t fVertexMixingBins[20]; // configure vertex bins for event mixing
184  Float_t fPtBins[19]; // pt bin borders
185  Int_t fNPtBins; // no of pt bins + 1
186  Double_t fCentrality; // event centrality
187  Double_t fVertex; // event vertex z
189  AliEventPoolManager *fPoolManager;
190  TList *fOutputList; // ! Output list
191  TH1F *fEventStats; // ! Histogram for event statistics
192  TH1F *fCentralityPass; // ! QA histogram of events that pass centrality cut
198  TH1F *fInvMNP[18];
199  TH1F *fInvMNN[18];
200  TH1F *fInvMPP[18];
201  TH1F *fPtSpectra[18];
202  TH1F *fPtP;
203  TH1F *fPtN;
204  TH1F *fPtKP;
205  TH1F *fPtKN;
208  Double_t fCentralityMin; // lower bound of cenrality bin
209  Double_t fCentralityMax; // upper bound of centrality bin
210  const char *fkCentralityMethodA; // method used to determine centrality, default
211  const char *fkCentralityMethodB; // method used to determine centrality, fallback
212  Bool_t fCentralityCut2010; // cut away the multiplicity outliers 2010
213  Bool_t fCentralityCut2011; // cut away the multiplicity outliers 2011
214  AliFlowTrackCuts *fPOICuts; // cuts for particles of interest (flow package)
215  Float_t fVertexRange; // absolute value of maximum distance of vertex along the z-axis
216  TH1F *fPhi;
217  TH1F *fPt;
218  TH1F *fEta;
219  TH1F *fVZEROA;
220  TH1F *fVZEROC;
221  TH1F *fTPCM;
222  // Float_t fDeltaDipAngle; // absolute value of delta dip angle to be excluded
223  // Float_t fDeltaDipPt; // upper value of pt range in which delta dip angle must be applied
224  // Bool_t fApplyDeltaDipCut; // enforce delta dip cut
226  TH1F *fDCAXYQA;
227  TH1F *fDCAZQA;
231  TProfile *fSubEventDPhiv2;
232  TProfile *fV0Data[18][2];
233  Bool_t fSkipEventSelection;// skip event selection and set bayesian pid object to MC mode
234  Bool_t fUsePidResponse;//use pid response instead of aliflowbayesianpid object for pid
235  AliPIDCombined* fPIDCombined; // pid combined
236  Bool_t fPileUp; // cut pile up
237  TF1 *fLowCut;
238  TF1 *fHighCut;
243  TH1F *fVertexZ;
244 
245  AliAnalysisTaskPhiFlow(const AliAnalysisTaskPhiFlow&); // Not implemented
246  AliAnalysisTaskPhiFlow& operator=(const AliAnalysisTaskPhiFlow&); // Not implemented
247  void MakeTrack(Double_t, Double_t, Double_t, Double_t, Int_t , Int_t[], Double_t p = 0., Double_t pz = 0.) const;
248 
250 };
251 
252 #endif
253 
254 
255 
Int_t charge
virtual Double_t Pt() const
Bool_t plpMV(const AliAODEvent *aod)
Bool_t PassesTPCbayesianCut(T *track) const
TProfile * fV0Data[18][2]
subevent resolution info for v2
Double_t fCentralityMin
QA profile of centralty vs multiplicity.
double Double_t
Definition: External.C:58
Bool_t CheckVertex(const T *event)
TH2F * fDCASecondaryWeak
dca of primaries (mc) or kaons (data)
TH1F * fPtKP
QA histogram of p_t distribution of negative particles.
Definition: External.C:236
AliFlowEvent * fFlowEvent
pid response object
void ReconstructionWithEventMixing(TObjArray *MixingCandidates) const
Double_t PhiPt(const T *track_1, const T *track_2) const
TH2F * fPIDk
QA histogram of TPC response of all charged particles.
void SetCurrentCentralityBin(Double_t c)
TH1F * fPtKN
QA histogram of p_t distribution of positive kaons.
ClassDef(AliAnalysisTaskPhiFlow, 8)
TH2F * fDCAPrim
qa plot of dca z
void PlotMultiplcities(const T *event) const
TH2F * fDCAMaterial
dca of weak (mc)
TCanvas * c
Definition: TestFitELoss.C:172
TH1F * fPtSpectra[18]
like-sign kaon pairs
TH1F * fVZEROC
QA plot vzeroa multiplicity (all tracks in event)
Double_t InvariantMass(const T *track1, const T *track2) const
TH2F * fPIDTOF
QA histo of TOF repsonse charged particles.
AliFlowTrackCuts * GetRPCuts() const
void SetMixingParameters(Int_t p[3])
void SetPOICuts(AliFlowTrackCuts *cutsPOI)
Double_t GetWDist(const AliVVertex *v0, const AliVVertex *v1)
UShort_t T(UShort_t m, UShort_t t)
Definition: RingBits.C:60
Bool_t fSkipEventSelection
profiles for vzero vn(minv)
void SetPOIDCAXYZ(Double_t dca[5])
TH1F * fPt
QA plot of azimuthal distribution of tracks used for event plane estimation.
virtual void Exec(Option_t *)
void GetMixingParameters(Int_t p[3]) const
int Int_t
Definition: External.C:63
Bool_t PassesDCACut(AliAODTrack *track) const
void PrepareFlowEvent(Int_t iMulti)
float Float_t
Definition: External.C:68
virtual Double_t Pz() const
void GetPIDConfiguration(Double_t prob[7]) const
virtual Double_t P() const
TH1F * fVZEROA
QA plot of eta distribution of tracks used for event plane estimation.
const char * GetCentralityMethod() const
AliPhiMesonHelperTrack(Float_t eta, Float_t phi, Float_t p, Float_t px, Float_t py, Float_t pz, Float_t pt, Int_t charge)
Bool_t SetIsMC(Bool_t ismc)
virtual Double_t Eta() const
TH2F * fMultvsCentr
QA profile global and tpc multiplicity after outlier cut.
TH1F * fEta
QA plot of p_t sectrum of tracks used for event plane estimation.
TH1F * fInvMNN[18]
unlike sign kaon pairs
AliEventPoolManager * InitializeEventMixing()
void InitializeHelperTrack(Float_t eta, Float_t phi, Float_t p, Float_t px, Float_t py, Float_t pz, Float_t pt, Int_t charge)
Bool_t PhiTrack(T *track) const
virtual Double_t Px() const
TH2F * fMultCorAfterCuts
QA histogram of p_t distribution of negative kaons.
short Short_t
Definition: External.C:23
TH1F * fDCAZQA
qa plot of dca xz
Bool_t fUsePidResponse
profiles for vzero vn(minv)
TObjArray * fCandidates
PID response object.
AliFlowTrackCuts * fPOICuts
QA profile of centralty vs multiplicity.
TH1F * fInvMNP[18]
QA histo of TOF response kaons.
Bool_t CheckCandidateEtaPtCut(const T *track1, const T *track2) const
Bool_t SetVZEROSubEvents(Bool_t v0)
void SetRPCuts(AliFlowTrackCuts *cutsRP)
virtual Double_t Phi() const
void SetCentralityWeights(TH1 *hist)
void InitializeBayesianPID(AliAODEvent *event)
void SetCommonConstants(Int_t massBins, Double_t minMass, Double_t maxMass)
TH1F * fPtN
QA histogram of p_t distribution of positive particles.
TH1F * BookHistogram(const char *name)
AliFlowBayesianPID * fBayesianResponse
flow events (one for each inv mass band)
void GetCandidateEtaAndPt(Double_t etapt[4]) const
ClassDef(AliPhiMesonHelperTrack, 1)
Bool_t GetDoubleCountingK(Double_t nSk, Short_t minNSigma) const
Double_t minMass
TH1F * InitPtSpectraHistograms(Float_t nmin, Float_t nmax)
void SetCandidateEtaAndPt(Double_t mineta, Double_t maxeta, Double_t minpt, Double_t maxpt)
TH1F * fDCAXYQA
qa dca of all charged particles
virtual Int_t Charge() const
void SetCentralityParameters(Double_t min, Double_t max, const char *a, const char *b, Bool_t c, Bool_t d)
TH2F * fNOPIDTOF
QA histogram of TPC response of kaons.
const char Option_t
Definition: External.C:48
Bool_t IsKaon(AliAODTrack *track) const
TH2F * fDCAAll
QA plot TPC multiplicity (tracks used for event plane estimation)
Double_t maxMass
TList * fOutputList
event pool manager
bool Bool_t
Definition: External.C:53
virtual void Terminate(Option_t *)
void SetPIDConfiguration(Double_t prob[7])
TProfile * fSubEventDPhiv2
dca material (mc) all (data)
void GetPOIDCZXYZ(Double_t dca[5]) const
TH1F * BookPtHistogram(const char *name)
TH1F * fTPCM
QA plot vzeroc multiplicity (all tracks in event)
Short_t FindMinNSigma(Double_t nSpi, Double_t nSk, Double_t nSp) const
TH2F * BookPIDHistogram(const char *name, Bool_t TPC)
void SetMixingBins(Int_t c[20], Int_t v[20])
Bool_t UseEventMixing(Bool_t mix, Bool_t type)
virtual void UserExec(Option_t *option)
AliEventPoolManager * fPoolManager
AOD oject.
Definition: External.C:196
virtual Double_t Py() const
AliFlowTrackCuts * GetPOICuts() const
void SetCandidateMinMaxY(Double_t min, Double_t max)
TH1F * fInvMPP[18]
like-sign kaon pairs
AliAnalysisTaskPhiFlow & operator=(const AliAnalysisTaskPhiFlow &)
Int_t SetDebugLevelPhiTask(Int_t debug)
void PtSelector(Int_t _track_type, const T *track_1, const T *track_2) const
TH2F * fNOPID
QA histogram of events that do not pass centrality cut.
void MakeTrack(Double_t, Double_t, Double_t, Double_t, Int_t, Int_t[], Double_t p=0., Double_t pz=0.) const
void SetPtBins(Float_t bin[19], Int_t n)