]>
Commit | Line | Data |
---|---|---|
ddc7ef22 | 1 | // Analysis task for basic QA exploiting symmetries |
2 | //of global, TPC, and ITS tracks | |
3 | ||
8d49bc89 | 4 | #ifndef ALIANALYSISTASKQASYM_H |
5 | #define ALIANALYSISTASKQASYM_H | |
9f9631f3 | 6 | |
7 | ||
8 | class TH1F; | |
9 | class TH2F; | |
1d55ca60 | 10 | class TH3F; |
9f9631f3 | 11 | class TList; |
9f9631f3 | 12 | |
13 | class AliESDEvent; | |
14 | class AliESDtrack; | |
15 | class AliESDtrackCuts; | |
16 | ||
17 | ||
18 | #include "AliAnalysisTaskSE.h" | |
1d55ca60 | 19 | #include "TFile.h" |
9f9631f3 | 20 | |
21 | class AliAnalysisTaskQASym : public AliAnalysisTaskSE { | |
22 | public: | |
1d55ca60 | 23 | AliAnalysisTaskQASym(const char *name = "AliAnalysisTaskQASym"); |
24 | virtual ~AliAnalysisTaskQASym() {} | |
25 | ||
26 | virtual void UserCreateOutputObjects(); | |
27 | virtual void UserExec(Option_t *option); | |
28 | virtual void Terminate(Option_t *); | |
1c91c6a4 | 29 | virtual void SetTrackType(Int_t type) {fTrackType = type;} // set the track type: global, its and tpc |
b676e981 | 30 | virtual void SetStandAloneTrack(Bool_t standAlone = kFALSE) {fStandAlone = standAlone;} //needed for ITS tracks |
1c91c6a4 | 31 | |
32 | virtual void SetNChargedRange(Int_t low = 0, Int_t high=1.0*1e7)//set multiplicity region of analysis | |
33 | // in order to compare different trigger settings | |
34 | { | |
35 | fLow = low; | |
36 | fHigh = high; | |
37 | } | |
9f9631f3 | 38 | |
1d55ca60 | 39 | virtual void SetCuts(AliESDtrackCuts* cuts) |
40 | {fCuts = cuts;} | |
41 | ||
42 | virtual void SetFieldOn(Bool_t b = kTRUE){fFieldOn = b;} | |
43 | ||
9f9631f3 | 44 | |
45 | private: | |
46 | ||
8d49bc89 | 47 | Int_t fTrackType; // track type |
1c91c6a4 | 48 | Bool_t fStandAlone; // needed for ITS tracks |
49 | Int_t fLow; // low Ncharges cut | |
50 | Int_t fHigh; // high Ncharges cut | |
8d49bc89 | 51 | Bool_t fFieldOn; // field status |
9f9631f3 | 52 | |
53 | TList *fHists; // List of histos | |
54 | ||
55 | //old | |
56 | TH1F *fHistRECpt; // pt | |
57 | TH1F *fEta; // eta | |
1c91c6a4 | 58 | |
59 | TH1F *fEtaWidth; // eta of tracks in Nch range (fLow <= Ncharged <= fHigh) | |
60 | TH1F *fPhiWidth; // phi of tracks in Nch range (fLow <= Ncharged <= fHigh) | |
61 | TH1F *fDcaWidth; // dca of tracks in Nch range (fLow <= Ncharged <= fHigh) | |
62 | TH1F *fPtWidth; // pt of tracks in Nch range (fLow <= Ncharged <= fHigh) | |
63 | ||
e4ab7356 | 64 | TH2F *fEtaPhi; // eta-phi |
9f9631f3 | 65 | TH1F *fEtaPt; // eta over pt |
66 | TH1F *fQPt; // charge over pt | |
67 | TH1F *fDca; // distance of closest approach | |
9f9631f3 | 68 | TH1F *fqRec; // reconstrcuted charge |
ddc7ef22 | 69 | TH1F *fSigmaPtHist; // sigma_pT |
1d55ca60 | 70 | |
71 | //positive und negative tracks | |
72 | TH1F *fRecPtPos; // pt of pos tracks | |
73 | TH1F *fRecPtNeg; // pt of neg tracks | |
74 | TH1F *fRecPhiPos; // phi of pos. tracks | |
75 | TH1F *fRecPhiNeg; // phi of neg. tracks | |
76 | TH1F *fRecEtaPos; // eta of neg. tracks | |
77 | TH1F *fRecEtaNeg; // eta of neg. tracks | |
78 | TH1F *fRecEtaPtPos; // eta over pt of neg. tracks | |
79 | TH1F *fRecEtaPtNeg; // eta over pt of neg. tracks | |
80 | TH1F *fRecDcaPos; // distance of closest approach of neg. tracks | |
81 | TH1F *fRecDcaNeg; // distance of closest of neg. tracks | |
82 | TH1F *fRecDcaNegInv; // invers dca of neg. tracks | |
83 | TH1F *fRecDPos; // impact parameter of neg. tracks | |
84 | TH1F *fRecDNeg; // impact parameter of neg. tracks | |
9f9631f3 | 85 | |
9f9631f3 | 86 | // two sides of TPC -> Eta/Theta |
1d55ca60 | 87 | TH1F *fRecQPtPosEta; //charge/pT for pos. eta |
88 | TH1F *fRecQPtNegEta; //charge/pT for neg. eta | |
89 | TH1F *fRecPtPosEta; //pT for pos. eta | |
90 | TH1F *fRecPtNegEta; //pT for neg. eta | |
91 | TH1F *fRecPhiPosEta; //phi for pos. eta | |
92 | TH1F *fRecPhiNegEta; //phi for neg. eta | |
93 | TH1F *fRecDcaPosEta; //dca for pos. eta | |
94 | TH1F *fRecDcaNegEta; //dca for neg. eta | |
95 | TH1F *fRecDPosEta; //d for pos. eta | |
96 | TH1F *fRecDNegEta; //d for neg. eta | |
97 | ||
98 | // 2D: all measures as function of z of first trackpoint | |
99 | TH2F *fRecPtPosVz; //pt-zfirst of pos tracks | |
8d49bc89 | 100 | TH2F *fRecPtNegVz; //pt-zfirst of neg tracks |
1d55ca60 | 101 | TH2F *fRecEtaPosVz; //eta-zfirst of pos tracks |
102 | TH2F *fRecEtaNegVz; //eta-zfirst of neg tracks | |
103 | TH2F *fRecPhiPosVz; //phi-zfirst of pos tracks | |
104 | TH2F *fRecPhiNegVz; //phi-zfirst of neg tracks | |
105 | TH2F *fSignedDcaPosVz; //dca-zfirst of pos tracks | |
106 | TH2F *fSignedDcaNegVz; //dca-zfirst of neg tracks | |
107 | TH2F *fRecQPtPosEtaVz; //charge/pT-zfirst of pos tracks | |
108 | TH2F *fRecQPtNegEtaVz; //charge/pT-zfirst of neg tracks | |
109 | TH2F *fRecEtaPtPosVz; //eta/pT-zfirst of pos tracks | |
110 | TH2F *fRecEtaPtNegVz; //eta/pT-zfirst of neg tracks | |
9f9631f3 | 111 | |
9f9631f3 | 112 | |
113 | //high | |
1d55ca60 | 114 | TH1F * fDeltaPhiAll; // phiLeaingTracks-phiOthers |
115 | TH2F * fDeltaPhiLeading; // phiLeaingTracks-phiOthers vs. phiLeading | |
116 | TH1F * fDiffDcaD; // d-dca | |
9f9631f3 | 117 | |
118 | //sim | |
1d55ca60 | 119 | TH1F * fPhiRec; //phi |
120 | TH1F * fThetaRec; //theta | |
121 | TH1F * fNumber; //Number of tracks per event | |
ddc7ef22 | 122 | TH1F * fNumberAfterCut; //Number of tracks per event after cuts |
1d55ca60 | 123 | TH1F * fVx; // x of first track point |
124 | TH1F * fVy; // y of first track point | |
125 | TH1F * fVz; // z of first track point | |
cbac53a7 | 126 | TH1F * fVertexX; // x of vertex |
127 | TH1F * fVertexY; // y of vertex | |
128 | TH1F * fVertexZ; // z of vertex | |
ddc7ef22 | 129 | TH1F * fNVertexSPD; //number of vertices SPD |
130 | TH1F * fNVertexTracks; //number of vertices of Tracks | |
1d55ca60 | 131 | |
132 | //new | |
133 | TH2F *fRecDcaPosPhi; //dca-phi for pos. | |
134 | TH2F *fRecDcaNegPhi; //dca-phi for neg. | |
135 | TH2F *fRecPtPosPhi; //pt-phi for pos. | |
136 | TH2F *fRecPtNegPhi; //pt-phi for neg. | |
137 | TH2F *fRecEtaPosPhi; //eta-phi for pos. | |
138 | TH2F *fRecEtaNegPhi; //eta-phi for neg. | |
139 | TH2F *fRecQPtPhi; //charge/pt-phi | |
140 | TH2F *fRecEtaPtPosPhi; //eta/pt-phi for neg. | |
141 | TH2F *fRecEtaPtNegPhi; //eta/pt-phi for pos. | |
142 | ||
143 | TH1F *fRecPtPosEtaPos; //pt for pos tracks and pos eta | |
144 | TH1F *fRecPtNegEtaPos; //pt for neg tracks and pos eta | |
145 | TH1F *fRecPtPosEtaNeg; //pt for pos tracks and neg eta | |
146 | TH1F *fRecPtNegEtaNeg; //pt for neg tracks and neg eta | |
147 | ||
148 | TH1F *fRec1PtPosEtaPos; //1/pt for pos tracks and pos eta | |
149 | TH1F *fRec1PtNegEtaPos; //1/pt for neg tracks and pos eta | |
150 | TH1F *fRec1PtPosEtaNeg; //1/pt for pos tracks and neg eta | |
151 | TH1F *fRec1PtNegEtaNeg; //1/pt for neg tracks and neg eta | |
152 | ||
153 | TH1F *fRecPhiPosEtaPos; //phi for pos tracks and pos eta | |
154 | TH1F *fRecPhiNegEtaPos; //phi for neg tracks and pos eta | |
155 | TH1F *fRecPhiPosEtaNeg; //phi for pos tracks and neg eta | |
156 | TH1F *fRecPhiNegEtaNeg; //phi for neg tracks and neg eta | |
157 | ||
158 | TH2F *fRecDcaPosPhiEtaPos; //dca-phi for pos tracks and pos eta | |
159 | TH2F *fRecDcaNegPhiEtaPos; //dca-phi for neg tracks and pos eta | |
160 | TH2F *fRecDcaPosPhiEtaNeg; //dca-phi for pos tracks and neg eta | |
161 | TH2F *fRecDcaNegPhiEtaNeg; //dca-phi for neg tracks and neg eta | |
162 | ||
98a5f582 | 163 | TH2F *fRecDcaPosPtEtaPos; //dca-pt for pos tracks and pos eta |
164 | TH2F *fRecDcaNegPtEtaPos; //dca-pt for neg tracks and pos eta | |
165 | TH2F *fRecDcaPosPtEtaNeg; //dca-pt for pos tracks and neg eta | |
166 | TH2F *fRecDcaNegPtEtaNeg; //dca-pt for neg tracks and neg eta | |
167 | ||
1d55ca60 | 168 | TH2F *fRecPtPosPhiEtaPos; //pt-phi for pos tracks and pos eta |
169 | TH2F *fRecPtNegPhiEtaPos; //pt-phi for neg tracks and pos eta | |
170 | TH2F *fRecPtPosPhiEtaNeg; //pt-phi for pos tracks and neg eta | |
171 | TH2F *fRecPtNegPhiEtaNeg; //pt-phi for neg tracks and neg eta | |
172 | ||
8ae6f5d0 | 173 | // TH3F *fRecDcaPhiPtPosEtaPos; //dca-pt-phi for pos tracks and pos eta |
174 | // TH3F *fRecDcaPhiPtNegEtaPos; //dca-pt-phi for neg tracks and pos eta | |
175 | // TH3F *fRecDcaPhiPtPosEtaNeg; //dca-pt-phi for pos tracks and neg eta | |
176 | // TH3F *fRecDcaPhiPtNegEtaNeg; //dca-pt-phi for neg tracks and neg eta | |
1d55ca60 | 177 | |
8d49bc89 | 178 | TH2F *fEtavPt; // eta vs pt |
ddc7ef22 | 179 | TH2F *fPhivPt; // phi vs pt |
180 | ||
8d49bc89 | 181 | TH2F *fCompareTPCparam; // TPC param |
ddc7ef22 | 182 | |
8d49bc89 | 183 | TH1F *fITSlayer; // ITS layer |
ddc7ef22 | 184 | TH2F *fITSlayerEta; // ITS layer vs eta |
185 | TH2F *fITSlayerPhi; // ITS layer vs phi | |
1d55ca60 | 186 | |
8d49bc89 | 187 | AliESDtrackCuts* fCuts; // List of cuts |
1d55ca60 | 188 | |
189 | // sectors of TPC | |
190 | TH1F *fRecPtTpcSector[18]; //pt for TPC sectors | |
191 | TH1F *fRecEtaTpcSector[18]; //eta for TPC sectors | |
192 | TH1F *fRecQPtTpcSector[18]; //charge/pt for TPC sectors | |
193 | TH1F *fRecEtaPtTpcSector[18]; //eta/pt for TPC sectors | |
194 | TH1F *fSignedDcaTpcSector[18]; //dca for TPC sectors | |
195 | ||
196 | ||
197 | // 7 different case of hit in ITS ladders | |
198 | TH1F *fRecPtPosLadder[7]; //pt for pos tracks | |
199 | TH1F *fRecPtNegLadder[7]; //pt for neg tracks | |
200 | TH1F *fRecPhiPosLadder[7]; //phi for pos tracks | |
201 | TH1F *fRecPhiNegLadder[7]; //phi for neg tracks | |
202 | TH1F *fRecEtaPosLadder[7]; //eta for pos tracks | |
203 | TH1F *fRecEtaNegLadder[7]; //eta for neg tracks | |
204 | TH1F *fSignDcaPos[7]; //dca for pos tracks | |
205 | TH1F *fSignDcaNeg[7]; //dca for neg tracks | |
206 | TH1F *fSignDcaNegInv[7]; //-dca for neg tracks | |
207 | TH1F *fPtSigmaPos[7]; //sigma_pT for pos tracks | |
208 | TH1F *fPtSigmaNeg[7]; //sigma_pT for neg tracks | |
209 | TH1F *fqPtRec[7]; // charge/pt | |
210 | TH2F *fDcaSigmaPos[7]; // dca - sigma_pT for pos tracks | |
211 | TH2F *fDcaSigmaNeg[7]; // dca - sigma_pT for neg tracks | |
212 | ||
1c91c6a4 | 213 | TH1F *fEtaBinPt[3][2]; // eta histogram for different pt bins and different charges |
214 | TH1F *fPhiBinPt[3][2]; // phi histogram for different pt bins and different charges | |
215 | TH1F *fDcaBinPt[3][2]; // dca histogram for different pt bins and different charges | |
216 | TH2F *fEtaPhiBinPt[3][2]; // eta-phi histogram for different pt bins and different charges | |
217 | ||
1d55ca60 | 218 | |
219 | ||
220 | ||
221 | ||
9f9631f3 | 222 | AliAnalysisTaskQASym(const AliAnalysisTaskQASym&); // not implemented |
223 | AliAnalysisTaskQASym& operator=(const AliAnalysisTaskQASym&); // not implemented | |
224 | ||
225 | ClassDef(AliAnalysisTaskQASym, 1); // Basic QA exploiting symmetries | |
226 | }; | |
227 | ||
228 | #endif |