AliPhysics  2c6b7ad (2c6b7ad)
AliAnalysisTaskCmeEse.h
Go to the documentation of this file.
1 #ifndef ALIANALYSISTASKCMEESE_H
2 #define ALIANALYSISTASKCMEESE_H
3 
4 // ROOT includes
5 #include <TList.h>
6 #include <TH1.h>
7 #include <TH2.h>
8 #include <TProfile.h>
9 #include <TF1.h>
10 #include <TSpline.h>
11 #include <THnSparse.h>
12 
13 // AliRoot includes
14 #include "AliAnalysisTaskSE.h"
15 #include "AliAODEvent.h"
16 #include "AliVHeader.h"
17 #include "AliVVertex.h"
18 #include "AliVEvent.h"
19 #include "AliVTrack.h"
20 
21 
23  public:
25  AliAnalysisTaskCmeEse(const char *name);
26 
27  virtual ~AliAnalysisTaskCmeEse();
28 
29  virtual void UserCreateOutputObjects();
30  virtual void UserExec(Option_t *option);
31  virtual void Terminate(Option_t *);
32 
33  virtual void SetVtxCut(Double_t vtxCut){fVtxCut = vtxCut;}
34  virtual void SetFilterbit(UInt_t filterbit){fFilterbit = filterbit;}
35  virtual void SetNoClus(Int_t noclus){fNoClus = noclus;}
36  virtual void SetEtaCut(Double_t etaCut){fEtaCut = etaCut;}
37  virtual void SetMinPt(Double_t minPt){fMinPt = minPt;}
38  virtual void SetMaxPt(Double_t maxPt){fMaxPt = maxPt;}
39  virtual void SetFlagLHC10h(Bool_t IsLHC10h){fLHC10h = IsLHC10h;}
40  virtual void SetFlagPileUp(Bool_t IsPileUP){fPileUp = IsPileUP;}
41  virtual void SetCentFlag(Short_t nCent){fCent = nCent;}
42  virtual void SetFlagRecEff(Bool_t IsRec){fRecEff = IsRec;}
43  virtual void SetFlagQAV0(Bool_t isQAV0){fQAV0 = isQAV0;}
44  virtual void SetFlagQATrk(Bool_t isTrkQA){fTrkQA = isTrkQA;}
45 
46 
47  private:
48  virtual Float_t GetVertex(AliAODEvent* aod) const;
49  virtual void Analyze(AliAODEvent* aod, Float_t vtxZ);
50  void OpenInfoCalbration(Int_t run);
51  Short_t GetPercCode(Double_t perc) const;
52  virtual Int_t GetTPCMult(AliVEvent* ev) const;
53  virtual Int_t GetGlobalMult(AliVEvent* ev) const;
54  Double_t GetWDist(const AliVVertex* v0, const AliVVertex* v1);
55  Bool_t plpMV(const AliVEvent *event);
56 
57 
59 
60  Int_t fRun; // run number - for calibration
61  TH1D* fMultV0; // profile from V0 multiplicity
62 
63  TH1D* fQx2mV0A; // <Qx2> V0A
64  TH1D* fQy2mV0A; // <Qy2> V0A
65  TH1D* fQx2sV0A; // sigma Qx2 V0A
66  TH1D* fQy2sV0A; // sigma Qy2 V0A
67  TH1D* fQx2mV0C; // <Qx2> V0C
68  TH1D* fQy2mV0C; // <Qy2> V0C
69  TH1D* fQx2sV0C; // sigma Qx2 V0C
70  TH1D* fQy2sV0C; // sigma Qy2 V0C
71 
72  // Cuts and options
73  TF1* fFitRecLow; // rec eff fit 1
74  TF1* fFitRecHigh; // rec eff fit 2
75  Bool_t fRecEff; // rec eff flag
76  Double_t fVtxCut; // Vtx cut on z position in cm
77  UInt_t fFilterbit; // filter bit
78  Double_t fEtaCut; // Eta cut used to select particles
79  Int_t fNoClus; // No of TPC clusters
80  Double_t fMinPt; // Min pt - for histogram limits
81  Double_t fMaxPt; // Max pt - for histogram limits
82  Bool_t fLHC10h; // flag to LHC10h data
83  Bool_t fPileUp; // flag for pileup
84  Short_t fCent; // centrality flag
85  Bool_t fQAV0; // qa V0 flag
86  Bool_t fTrkQA; // qa for tracks
87  TSpline3* fSplQ2c[90]; // splines for q2 V0C
88 
89  // Output objects
93  TH1F* fVtxAfterCuts;
96 
98 
99  THnSparseF* fAllQA;
100 
105 
110 
111 
116 
121 
122 
123  TProfile* fV2A;
124  TProfile* fV0AV0Cv2;
125  TProfile* fV0ATPCv2;
126  TProfile* fV0CTPCv2;
127 
128  TProfile* fCMESameQPos;
129  TProfile* fCMESameQNeg;
130  TProfile* fCMEOppQ;
131 
132  TProfile* fCMESameQPosCos;
133  TProfile* fCMESameQPosSin;
134  TProfile* fCMESameQPosv1;
135 
136  TProfile* fCMESameQNegCos;
137  TProfile* fCMESameQNegSin;
138  TProfile* fCMESameQNegv1;
139 
140  TProfile* fCMEOppQCos;
141  TProfile* fCMEOppQSin;
142  TProfile* fCMEOppQv1;
143 
144 
145 
146  TProfile* fV2Aqc2[10];
147  TProfile* fV0AV0Cv2qc2[10];
148  TProfile* fV0ATPCv2qc2[10];
149  TProfile* fV0CTPCv2qc2[10];
150 
151  TProfile* fCMESameQPosqc2[10];
152  TProfile* fCMESameQNegqc2[10];
153  TProfile* fCMEOppQqc2[10];
154 
155  TProfile* fCMESameQPosCosqc2[10];
156  TProfile* fCMESameQPosSinqc2[10];
157  TProfile* fCMESameQPosv1qc2[10];
158 
159  TProfile* fCMESameQNegCosqc2[10];
160  TProfile* fCMESameQNegSinqc2[10];
161  TProfile* fCMESameQNegv1qc2[10];
162 
163  TProfile* fCMEOppQCosqc2[10];
164  TProfile* fCMEOppQSinqc2[10];
165  TProfile* fCMEOppQv1qc2[10];
166 
167 
168  TProfile* fV2Pt[9];
169 
170  TProfile* fCMESameQPosEta[9];
171  TProfile* fCMESameQPosPtDif[9];
172  TProfile* fCMESameQPosPtSum[9];
173 
174  TProfile* fCMESameQNegEta[9];
175  TProfile* fCMESameQNegPtDif[9];
176  TProfile* fCMESameQNegPtSum[9];
177 
178  TProfile* fCMEOppQEta[9];
179  TProfile* fCMEOppQPtDif[9];
180  TProfile* fCMEOppQPtSum[9];
181 
182 
183  TH1D* fPsiA[10];
184  TH1D* fPsiC[10];
186 
187 
188  TProfile* fV2Ptqc2[9][10];
189 
190  TProfile* fCMESameQPosEtaqc2[9][10];
191  TProfile* fCMESameQPosPtDifqc2[9][10];
192  TProfile* fCMESameQPosPtSumqc2[9][10];
193 
194  TProfile* fCMESameQNegEtaqc2[9][10];
195  TProfile* fCMESameQNegPtDifqc2[9][10];
196  TProfile* fCMESameQNegPtSumqc2[9][10];
197 
198  TProfile* fCMEOppQEtaqc2[9][10];
199  TProfile* fCMEOppQPtDifqc2[9][10];
200  TProfile* fCMEOppQPtSumqc2[9][10];
201 
202 
205 
206  ClassDef(AliAnalysisTaskCmeEse, 1);
207 };
208 
209 #endif
TProfile * fCMESameQNegEta[9]
cme same charge pos pt sum
TProfile * fCMESameQPosEta[9]
v2 pt unbiased
virtual void SetMaxPt(Double_t maxPt)
double Double_t
Definition: External.C:58
TH1F * fVtxAfterCuts
Vertex z dist before cuts.
Definition: External.C:236
TProfile * fCMESameQNeg
cme same charge pos
TProfile * fV2Pt[9]
cme opp charge v1
TH2D * fQxavsV0Aft
corrected Qx V0C vtx
virtual Int_t GetGlobalMult(AliVEvent *ev) const
TH1D * fPsiA[10]
cme opp charge pt sum
TProfile * fV2A
corrected Qx V0C vtx
TProfile * fCMESameQNegPtDif[9]
cme same charge neg eta
TProfile * fCMEOppQSin
cme opp charge <cos*cos>
TProfile * fCMESameQNegPtDifqc2[9][10]
cme same charge neg ESE
virtual Float_t GetVertex(AliAODEvent *aod) const
virtual void SetFilterbit(UInt_t filterbit)
TProfile * fCMESameQPosSin
cme same charge pos <cos*cos>
TProfile * fV0CTPCv2
correlation V0A-TPC for resolution
TProfile * fCMEOppQCos
cme same charge neg v1
TProfile * fCMEOppQv1qc2[10]
cme opp charge <sin*sin>
TProfile * fCMESameQNegCosqc2[10]
cme same charge pos v1
TProfile * fCMESameQNegSin
cme same charge neg <cos*cos>
TProfile * fCMEOppQPtSumqc2[9][10]
cme opp charge ESE
TProfile * fCMEOppQPtDifqc2[9][10]
cme opp charge ESE
TProfile * fCMESameQPosPtSumqc2[9][10]
cme same charge pos ESE
virtual void UserCreateOutputObjects()
TH2D * fQycvsVtxZBef
corrected Qx V0C vtx
Bool_t plpMV(const AliVEvent *event)
TH2D * fQycvsV0Bef
corrected Qx V0C cent
TProfile * fCMEOppQEtaqc2[9][10]
cme same charge neg ESE
TProfile * fCMESameQPosqc2[10]
correlation V0C-TPC for resolution ESE
TProfile * fCMESameQPos
correlation V0C-TPC for resolution
TProfile * fCMEOppQPtSum[9]
cme opp charge pt dif
TH1F * fVtxBeforeCuts
Event vertex info.
TProfile * fCMEOppQ
cme same charge neg
TH2D * fQxcvsVtxZAft
corrected Qy V0A vtx
THnSparseF * fAllQA
percentile distribution
TProfile * fCMEOppQCosqc2[10]
cme same charge neg v1
TH2D * fQycvsVtxZAft
corrected Qx V0C vtx
int Int_t
Definition: External.C:63
TProfile * fV0ATPCv2qc2[10]
correlation V0A-V0C for resolution ESE
TProfile * fCMESameQNegCos
cme same charge pos v1
TProfile * fCMEOppQPtDif[9]
cme opp charge eta
Definition: External.C:204
virtual Int_t GetTPCMult(AliVEvent *ev) const
TH2D * fQyavsVtxZAft
corrected Qx V0A vtx
unsigned int UInt_t
Definition: External.C:33
TH2D * fPsiAvsPsiC[10]
Psi V0C.
TH2D * fQycvsV0Aft
corrected Qx V0C cent
virtual void SetFlagRecEff(Bool_t IsRec)
float Float_t
Definition: External.C:68
TH1D * fPsiC[10]
Psi V0A.
TProfile * fV0CTPCv2qc2[10]
correlation V0A-TPC for resolution ESE
TH2D * fQxcvsV0Aft
corrected Qy V0A cent
TProfile * fCMESameQPosPtDifqc2[9][10]
cme same charge pos ESE
TProfile * fCMESameQPosSinqc2[10]
cme same charge pos <cos*cos>
virtual void SetVtxCut(Double_t vtxCut)
TProfile * fV0AV0Cv2qc2[10]
v2 V0A ESE
Definition: External.C:228
Definition: External.C:212
TH2D * fQxavsVtxZBef
corrected Qx V0C cent
TProfile * fCMESameQPosEtaqc2[9][10]
v2 pt ESE
virtual void Analyze(AliAODEvent *aod, Float_t vtxZ)
virtual void SetEtaCut(Double_t etaCut)
TH2D * fPercqc2
correlation between TPC and global track multiplicity after cuts
TProfile * fCMESameQPosPtSum[9]
cme same charge pos pt dif
virtual void UserExec(Option_t *option)
TProfile * fV0ATPCv2
correlation V0A-V0C for resolution
virtual void SetFlagPileUp(Bool_t IsPileUP)
TProfile * fCMESameQPosPtDif[9]
cme same charge pos eta
TProfile * fCMESameQPosv1
cme same charge pos <sin*sin>
TH2D * fQxavsVtxZAft
corrected Qx V0C cent
short Short_t
Definition: External.C:23
TH2D * fQxavsV0Bef
qa histo for tracks
TProfile * fCMEOppQv1
cme opp charge <sin*sin>
TProfile * fCMESameQPosCosqc2[10]
cme opp charge ESE
TH1I * fVtx
Output list of objects.
TProfile * fCMESameQNegv1
cme same charge neg <sin*sin>
TProfile * fCMEOppQSinqc2[10]
cme opp charge <cos*cos>
TProfile * fCMESameQNegPtSum[9]
cme same charge neg pt dif
AliAnalysisTaskCmeEse & operator=(const AliAnalysisTaskCmeEse &)
TH2D * fQxcvsV0Bef
corrected Qy V0A cent
TProfile * fCMESameQNegqc2[10]
cme same charge pos ESE
virtual void SetCentFlag(Short_t nCent)
TProfile * fV2Aqc2[10]
cme opp charge v1
Double_t GetWDist(const AliVVertex *v0, const AliVVertex *v1)
TH2F * fMultCorAfterCuts
correlation between TPC and global track multiplicity before cuts
TProfile * fCMESameQNegSinqc2[10]
cme same charge neg <cos*cos>
const char Option_t
Definition: External.C:48
TProfile * fCMEOppQEta[9]
cme same charge neg pt sum
TH2D * fQxcvsVtxZBef
corrected Qy V0A vtx
Short_t GetPercCode(Double_t perc) const
TProfile * fCMESameQNegv1qc2[10]
cme same charge neg <sin*sin>
bool Bool_t
Definition: External.C:53
virtual void SetNoClus(Int_t noclus)
TH2D * fQyavsV0Bef
corrected Qx V0A cent
TProfile * fCMESameQNegPtSumqc2[9][10]
cme same charge neg ESE
TProfile * fCMESameQPosv1qc2[10]
cme same charge pos <sin*sin>
virtual void SetFlagQATrk(Bool_t isTrkQA)
TProfile * fCMEOppQqc2[10]
cme same charge neg ESE
TProfile * fV2Ptqc2[9][10]
correlation PsiA-PsiC
virtual void Terminate(Option_t *)
TH2D * fQyavsV0Aft
corrected Qx V0A cent
TProfile * fCMESameQPosCos
cme opp charge
virtual void SetFlagLHC10h(Bool_t IsLHC10h)
virtual void SetMinPt(Double_t minPt)
TProfile * fCMESameQNegEtaqc2[9][10]
cme same charge pos ESE
TH2D * fQyavsVtxZBef
corrected Qx V0A vtx
TH2F * fMultCorBeforeCuts
Vertex z dist after cuts.
virtual void SetFlagQAV0(Bool_t isQAV0)