]>
Commit | Line | Data |
---|---|---|
f9a6cab5 | 1 | /************************************************************************** |
2 | * Authors : Antonin Maire, Boris Hippolyte * | |
3 | * Contributors are mentioned in the code where appropriate. * | |
4 | * * | |
5 | * Permission to use, copy, modify and distribute this software and its * | |
6 | * documentation strictly for non-commercial purposes is hereby granted * | |
7 | * without fee, provided that the above copyright notice appears in all * | |
8 | * copies and that both the copyright notice and this permission notice * | |
9 | * appear in the supporting documentation. The authors make no claims * | |
10 | * about the suitability of this software for any purpose. It is * | |
11 | * provided "as is" without express or implied warranty. * | |
12 | **************************************************************************/ | |
13 | ||
14 | //----------------------------------------------------------------- | |
15 | // AliAnalysisTaskCheckCascadePbPb class | |
16 | // | |
17 | // Origin AliAnalysisTaskCheckCascade which has four roles : | |
18 | // 1. QAing the Cascades from ESD and AOD | |
19 | // Origin: AliAnalysisTaskESDCheckV0 by Boris Hippolyte Nov2007, hippolyt@in2p3.fr | |
20 | // 2. Prepare the plots which stand as raw material for yield extraction (wi/wo PID) | |
21 | // 3. Supply an AliCFContainer meant to define the optimised topological selections | |
22 | // 4. Rough azimuthal correlation study (Eta, Phi) | |
23 | // Adapted to Cascade : A.Maire Mar2008, antonin.maire@ires.in2p3.fr | |
24 | // Modified : A.Maire Mar2010 | |
25 | // | |
39b18694 | 26 | // Adapted to PbPb analysis: M. Nicassio, maria.nicassio@ba.infn.it |
27 | // Feb-August2011 | |
28 | // - Physics selection moved to the run.C macro | |
29 | // - Centrality selection added (+ setters) and histos | |
30 | // - flag and setters added (CF container usage, vertex range) | |
31 | // - histo added and histo/container binning changed | |
32 | // - protection in the destructor for CAF usage | |
33 | // - AliWarning disabled | |
34 | // - number of tracklets from AOD also | |
35 | // - automatic settings for PID | |
36 | // September2011 | |
37 | // - proper time histos/container added (V0 and Cascades) | |
38 | // - cosine PA V0 wrt Xi vertex in the container | |
39 | // November2011 | |
ff9a705e | 40 | // - re-run V0's and cascade's vertexers (SetCuts instead of SetDefaultCuts!!) |
39b18694 | 41 | // - AOD analysis part completed |
f9a6cab5 | 42 | //----------------------------------------------------------------- |
43 | ||
44 | class TTree; | |
45 | class TParticle; | |
46 | class TVector3; | |
47 | ||
48 | class AliESDVertex; | |
49 | class AliAODVertex; | |
50 | class AliESDv0; | |
51 | class AliAODv0; | |
52 | ||
53 | #include <Riostream.h> | |
54 | #include "TList.h" | |
55 | #include "TH1.h" | |
56 | #include "TH2.h" | |
57 | #include "TH3.h" | |
58 | #include "THnSparse.h" | |
59 | #include "TVector3.h" | |
60 | #include "TCanvas.h" | |
61 | #include "TMath.h" | |
62 | #include "TLegend.h" | |
63 | ||
64 | ||
65 | #include "AliLog.h" | |
66 | #include "AliCentrality.h" | |
67 | #include "AliESDEvent.h" | |
68 | #include "AliAODEvent.h" | |
39b18694 | 69 | #include "AliV0vertexer.h" |
70 | #include "AliCascadeVertexer.h" | |
11bcd1e4 | 71 | #include "AliESDtrackCuts.h" |
72 | #include "AliPIDResponse.h" | |
f9a6cab5 | 73 | |
74 | #include "AliESDVZERO.h" | |
75 | ||
76 | #include "AliInputEventHandler.h" | |
77 | #include "AliAnalysisManager.h" | |
78 | #include "AliMCEventHandler.h" | |
39b18694 | 79 | #include "AliESDInputHandler.h" |
80 | #include "AliAODInputHandler.h" | |
f9a6cab5 | 81 | #include "AliCFContainer.h" |
82 | #include "AliMultiplicity.h" | |
83 | ||
84 | #include "AliESDcascade.h" | |
85 | #include "AliAODcascade.h" | |
39b18694 | 86 | #include "AliAODTrack.h" |
f9a6cab5 | 87 | |
88 | #include "AliAnalysisTaskCheckCascadePbPb.h" | |
89 | ||
90 | ClassImp(AliAnalysisTaskCheckCascadePbPb) | |
91 | ||
92 | ||
93 | ||
94 | //________________________________________________________________________ | |
95 | AliAnalysisTaskCheckCascadePbPb::AliAnalysisTaskCheckCascadePbPb() | |
39b18694 | 96 | : AliAnalysisTaskSE(), fAnalysisType("ESD"), fESDtrackCuts(0),/*fPaveTextBookKeeping(0),*/ |
11bcd1e4 | 97 | fPIDResponse (0), |
f9a6cab5 | 98 | fkRerunV0CascVertexers (0), |
99 | fkQualityCutZprimVtxPos (kTRUE), | |
100 | fkQualityCutNoTPConlyPrimVtx(kTRUE), | |
101 | fkQualityCutTPCrefit (kTRUE), | |
102 | fkQualityCut80TPCcls (kTRUE), | |
103 | fkExtraSelections (0), | |
104 | fCentrLowLim (0), | |
105 | fCentrUpLim (0), | |
106 | fCentrEstimator (0), | |
107 | fVtxRange (0), | |
f9a6cab5 | 108 | |
109 | ||
110 | // - Cascade part initialisation | |
111 | fListHistCascade(0), | |
11bcd1e4 | 112 | fHistEvtsInCentralityBinsvsNtracks(0), |
f9a6cab5 | 113 | fHistCascadeMultiplicityBeforeEvSel(0), |
114 | fHistCascadeMultiplicityForCentrEvt(0), fHistTrackMultiplicityForCentrEvt(0), fHistTPCrefitTrackMultiplicityForCentrEvt(0), | |
f9a6cab5 | 115 | fHistCascadeMultiplicityForSelEvt(0), |
116 | fHistPosBestPrimaryVtxXForSelEvt(0), fHistPosBestPrimaryVtxYForSelEvt(0), fHistPosBestPrimaryVtxZForSelEvt(0), | |
117 | fHistTPCrefitTrackMultiplicityForCascadeEvt(0), | |
118 | fHistPosV0TPCClusters(0), fHistNegV0TPCClusters(0), fHistBachTPCClusters(0), | |
119 | fHistVtxStatus(0), | |
120 | ||
121 | fHistPosTrkgPrimaryVtxXForCascadeEvt(0), fHistPosTrkgPrimaryVtxYForCascadeEvt(0), fHistPosTrkgPrimaryVtxZForCascadeEvt(0), fHistTrkgPrimaryVtxRadius(0), | |
122 | fHistPosBestPrimaryVtxXForCascadeEvt(0), fHistPosBestPrimaryVtxYForCascadeEvt(0), fHistPosBestPrimaryVtxZForCascadeEvt(0), fHistBestPrimaryVtxRadius(0), | |
123 | f2dHistTrkgPrimVtxVsBestPrimVtx(0), | |
124 | ||
125 | fHistEffMassXi(0), fHistChi2Xi(0), | |
126 | ||
127 | fHistDcaXiDaughters(0), fHistDcaBachToPrimVertex(0), fHistXiCosineOfPointingAngle(0), fHistXiRadius(0), | |
128 | ||
129 | fHistMassLambdaAsCascDghter(0), | |
130 | fHistV0Chi2Xi(0), | |
131 | fHistDcaV0DaughtersXi(0), | |
132 | fHistDcaV0ToPrimVertexXi(0), | |
133 | fHistV0CosineOfPointingAngleXi(0), | |
134 | fHistV0RadiusXi(0), | |
135 | fHistDcaPosToPrimVertexXi(0), fHistDcaNegToPrimVertexXi(0), | |
136 | ||
137 | fHistMassXiMinus(0), fHistMassXiPlus(0), | |
138 | fHistMassOmegaMinus(0), fHistMassOmegaPlus(0), | |
139 | fHistMassWithCombPIDXiMinus(0), fHistMassWithCombPIDXiPlus(0), | |
140 | fHistMassWithCombPIDOmegaMinus(0), fHistMassWithCombPIDOmegaPlus(0), | |
141 | ||
142 | fHistXiTransvMom(0), fHistXiTotMom(0), | |
143 | fHistBachTransvMomXi(0), fHistBachTotMomXi(0), | |
144 | ||
145 | fHistChargeXi(0), | |
146 | fHistV0toXiCosineOfPointingAngle(0), | |
147 | ||
148 | fHistRapXi(0), fHistRapOmega(0), fHistEtaXi(0), | |
149 | fHistThetaXi(0), fHistPhiXi(0), | |
150 | ||
151 | f2dHistArmenteros(0), | |
152 | f2dHistEffMassLambdaVsEffMassXiMinus(0), f2dHistEffMassXiVsEffMassOmegaMinus(0), | |
153 | f2dHistEffMassLambdaVsEffMassXiPlus(0), f2dHistEffMassXiVsEffMassOmegaPlus(0), | |
154 | f2dHistXiRadiusVsEffMassXiMinus(0), f2dHistXiRadiusVsEffMassXiPlus(0), | |
155 | f2dHistXiRadiusVsEffMassOmegaMinus(0), f2dHistXiRadiusVsEffMassOmegaPlus(0), | |
156 | ||
157 | f2dHistTPCdEdxOfCascDghters(0), | |
158 | ||
159 | // f3dHistXiPtVsEffMassVsYXiMinus(0), f3dHistXiPtVsEffMassVsYXiPlus(0), | |
160 | // f3dHistXiPtVsEffMassVsYOmegaMinus(0), f3dHistXiPtVsEffMassVsYOmegaPlus(0), | |
161 | ||
162 | fCFContCascadePIDXiMinus(0), | |
163 | fCFContCascadePIDXiPlus(0), | |
164 | fCFContCascadePIDOmegaMinus(0), | |
165 | fCFContCascadePIDOmegaPlus(0), | |
166 | fCFContCascadeCuts(0), | |
167 | ||
168 | // fHnSpAngularCorrXiMinus(0), fHnSpAngularCorrXiPlus(0), | |
169 | // fHnSpAngularCorrOmegaMinus(0), fHnSpAngularCorrOmegaPlus(0), | |
170 | fV0Ampl(0), | |
171 | ||
172 | fHistDcaXiDaughtersvsInvMass(0), fHistDcaBachToPrimVertexvsInvMass(0), fHistXiCosineOfPointingAnglevsInvMass(0), | |
173 | fHistMassLambdaAsCascDghtervsInvMass(0),fHistDcaV0DaughtersXivsInvMass(0),fHistDcaV0ToPrimVertexXivsInvMass(0) | |
174 | ||
175 | ||
176 | { | |
177 | // Dummy Constructor | |
f9a6cab5 | 178 | for(Int_t iV0selIdx = 0; iV0selIdx < 7; iV0selIdx++ ) { fV0Sels [iV0selIdx ] = -1.; } |
179 | for(Int_t iCascSelIdx = 0; iCascSelIdx < 8; iCascSelIdx++ ) { fCascSels [iCascSelIdx ] = -1.; } | |
180 | } | |
181 | ||
182 | ||
183 | //________________________________________________________________________ | |
184 | AliAnalysisTaskCheckCascadePbPb::AliAnalysisTaskCheckCascadePbPb(const char *name) | |
39b18694 | 185 | : AliAnalysisTaskSE(name), fAnalysisType("ESD"), fESDtrackCuts(0), /*fPaveTextBookKeeping(0),*/ |
11bcd1e4 | 186 | fPIDResponse (0), |
f9a6cab5 | 187 | fkRerunV0CascVertexers (0), |
188 | fkQualityCutZprimVtxPos (kTRUE), | |
189 | fkQualityCutNoTPConlyPrimVtx(kTRUE), | |
190 | fkQualityCutTPCrefit (kTRUE), | |
191 | fkQualityCut80TPCcls (kTRUE), | |
192 | fkExtraSelections (0), | |
193 | fCentrLowLim (0), | |
194 | fCentrUpLim (0), | |
195 | fCentrEstimator (0), | |
196 | fVtxRange (0), | |
f9a6cab5 | 197 | |
198 | // - Cascade part initialisation | |
199 | fListHistCascade(0), | |
11bcd1e4 | 200 | fHistEvtsInCentralityBinsvsNtracks(0), |
f9a6cab5 | 201 | fHistCascadeMultiplicityBeforeEvSel(0), |
202 | fHistCascadeMultiplicityForCentrEvt(0), fHistTrackMultiplicityForCentrEvt(0), fHistTPCrefitTrackMultiplicityForCentrEvt(0), | |
f9a6cab5 | 203 | fHistCascadeMultiplicityForSelEvt(0), |
204 | fHistPosBestPrimaryVtxXForSelEvt(0), fHistPosBestPrimaryVtxYForSelEvt(0), fHistPosBestPrimaryVtxZForSelEvt(0), | |
205 | fHistTPCrefitTrackMultiplicityForCascadeEvt(0), | |
206 | fHistPosV0TPCClusters(0), fHistNegV0TPCClusters(0), fHistBachTPCClusters(0), | |
207 | fHistVtxStatus(0), | |
208 | ||
209 | fHistPosTrkgPrimaryVtxXForCascadeEvt(0), fHistPosTrkgPrimaryVtxYForCascadeEvt(0), fHistPosTrkgPrimaryVtxZForCascadeEvt(0), fHistTrkgPrimaryVtxRadius(0), | |
210 | fHistPosBestPrimaryVtxXForCascadeEvt(0), fHistPosBestPrimaryVtxYForCascadeEvt(0), fHistPosBestPrimaryVtxZForCascadeEvt(0), fHistBestPrimaryVtxRadius(0), | |
211 | f2dHistTrkgPrimVtxVsBestPrimVtx(0), | |
212 | ||
213 | fHistEffMassXi(0), fHistChi2Xi(0), | |
214 | fHistDcaXiDaughters(0), fHistDcaBachToPrimVertex(0), fHistXiCosineOfPointingAngle(0), fHistXiRadius(0), | |
215 | ||
216 | fHistMassLambdaAsCascDghter(0), | |
217 | fHistV0Chi2Xi(0), | |
218 | fHistDcaV0DaughtersXi(0), | |
219 | fHistDcaV0ToPrimVertexXi(0), | |
220 | fHistV0CosineOfPointingAngleXi(0), | |
221 | fHistV0RadiusXi(0), | |
222 | fHistDcaPosToPrimVertexXi(0), fHistDcaNegToPrimVertexXi(0), | |
223 | ||
224 | fHistMassXiMinus(0), fHistMassXiPlus(0), | |
225 | fHistMassOmegaMinus(0), fHistMassOmegaPlus(0), | |
226 | fHistMassWithCombPIDXiMinus(0), fHistMassWithCombPIDXiPlus(0), | |
227 | fHistMassWithCombPIDOmegaMinus(0), fHistMassWithCombPIDOmegaPlus(0), | |
228 | ||
229 | fHistXiTransvMom(0), fHistXiTotMom(0), | |
230 | fHistBachTransvMomXi(0), fHistBachTotMomXi(0), | |
231 | ||
232 | fHistChargeXi(0), | |
233 | fHistV0toXiCosineOfPointingAngle(0), | |
234 | ||
235 | fHistRapXi(0), fHistRapOmega(0), fHistEtaXi(0), | |
236 | fHistThetaXi(0), fHistPhiXi(0), | |
237 | ||
238 | f2dHistArmenteros(0), | |
239 | f2dHistEffMassLambdaVsEffMassXiMinus(0), f2dHistEffMassXiVsEffMassOmegaMinus(0), | |
240 | f2dHistEffMassLambdaVsEffMassXiPlus(0), f2dHistEffMassXiVsEffMassOmegaPlus(0), | |
241 | f2dHistXiRadiusVsEffMassXiMinus(0), f2dHistXiRadiusVsEffMassXiPlus(0), | |
242 | f2dHistXiRadiusVsEffMassOmegaMinus(0), f2dHistXiRadiusVsEffMassOmegaPlus(0), | |
243 | ||
244 | f2dHistTPCdEdxOfCascDghters(0), | |
245 | ||
246 | // f3dHistXiPtVsEffMassVsYXiMinus(0), f3dHistXiPtVsEffMassVsYXiPlus(0), | |
247 | // f3dHistXiPtVsEffMassVsYOmegaMinus(0), f3dHistXiPtVsEffMassVsYOmegaPlus(0), | |
248 | ||
249 | fCFContCascadePIDXiMinus(0), | |
250 | fCFContCascadePIDXiPlus(0), | |
251 | fCFContCascadePIDOmegaMinus(0), | |
252 | fCFContCascadePIDOmegaPlus(0), | |
253 | fCFContCascadeCuts(0), | |
254 | ||
255 | // fHnSpAngularCorrXiMinus(0), fHnSpAngularCorrXiPlus(0), | |
256 | // fHnSpAngularCorrOmegaMinus(0), fHnSpAngularCorrOmegaPlus(0), | |
257 | fV0Ampl(0), | |
258 | ||
259 | fHistDcaXiDaughtersvsInvMass(0), fHistDcaBachToPrimVertexvsInvMass(0), fHistXiCosineOfPointingAnglevsInvMass(0), | |
260 | fHistMassLambdaAsCascDghtervsInvMass(0),fHistDcaV0DaughtersXivsInvMass(0),fHistDcaV0ToPrimVertexXivsInvMass(0) | |
261 | ||
262 | ||
263 | { | |
264 | // Constructor | |
265 | ||
266 | // Define input and output slots here | |
267 | // Input slot #0 works with a TChain | |
268 | // DefineInput(0, TChain::Class()); | |
269 | // Output slot #1 writes into a TList container (cascade) | |
39b18694 | 270 | // default Pb-Pb values |
f9a6cab5 | 271 | fV0Sels[0] = 33. ; // max allowed chi2 |
39b18694 | 272 | fV0Sels[1] = 0.1; // min allowed impact parameter for the 1st daughter |
273 | fV0Sels[2] = 0.1; // min allowed impact parameter for the 2nd daughter | |
274 | fV0Sels[3] = 1.; // max allowed DCA between the daughter tracks | |
275 | fV0Sels[4] = .998 ;// min allowed cosine of V0's pointing angle | |
276 | fV0Sels[5] = 0.9 ; // min radius of the fiducial volume | |
277 | fV0Sels[6] = 100. ; // max radius of the fiducial volume | |
f9a6cab5 | 278 | |
279 | fCascSels[0] = 33. ; // max allowed chi2 (same as PDC07) | |
39b18694 | 280 | fCascSels[1] = 0.05; // min allowed V0 impact parameter |
281 | fCascSels[2] = 0.008; // "window" around the Lambda mass | |
282 | fCascSels[3] = 0.03; // min allowed bachelor's impact parameter //check cuts | |
283 | fCascSels[4] = 0.3 ; // max allowed DCA between the V0 and the bachelor | |
284 | fCascSels[5] = 0.999 ;// min allowed cosine of the cascade pointing angle | |
285 | fCascSels[6] = 0.9 ; // min radius of the fiducial volume | |
286 | fCascSels[7] = 100. ; // max radius of the fiducial volume | |
f9a6cab5 | 287 | |
288 | // Output slot #0 writes into a TList container (Cascade) | |
289 | DefineOutput(1, TList::Class()); | |
290 | /*DefineOutput(2, TPaveText::Class());*/ | |
291 | AliLog::SetClassDebugLevel("AliAnalysisTaskCheckCascadePbPb",1); // MN this should (?) enable only AliFatal | |
292 | } | |
293 | ||
294 | ||
39b18694 | 295 | AliAnalysisTaskCheckCascadePbPb::~AliAnalysisTaskCheckCascadePbPb() { |
f9a6cab5 | 296 | // |
297 | // Destructor | |
298 | // | |
299 | ||
300 | // For all TH1, 2, 3 HnSparse and CFContainer are in the fListCascade TList. | |
301 | // They will be deleted when fListCascade is deleted by the TSelector dtor | |
302 | // Because of TList::SetOwner() ... | |
303 | ||
304 | if (fListHistCascade && !AliAnalysisManager::GetAnalysisManager()->IsProofMode()) { delete fListHistCascade; fListHistCascade = 0x0; } | |
39b18694 | 305 | if (fESDtrackCuts) { delete fESDtrackCuts; fESDtrackCuts = 0x0; } |
f9a6cab5 | 306 | //if (fPaveTextBookKeeping) { delete fPaveTextBookKeeping; fPaveTextBookKeeping = 0x0;} // fPaveTextBookKeeping is not strored in the TList |
307 | } | |
308 | ||
309 | ||
310 | ||
311 | //________________________________________________________________________ | |
39b18694 | 312 | void AliAnalysisTaskCheckCascadePbPb::UserCreateOutputObjects() { |
f9a6cab5 | 313 | // Create histograms |
314 | // Called once | |
315 | ||
316 | ||
317 | ||
318 | fListHistCascade = new TList(); | |
319 | fListHistCascade->SetOwner(); // See http://root.cern.ch/root/html/TCollection.html#TCollection:SetOwner | |
39b18694 | 320 | |
11bcd1e4 | 321 | // New PID object |
39b18694 | 322 | AliAnalysisManager *man=AliAnalysisManager::GetAnalysisManager(); |
323 | AliInputEventHandler* inputHandler = (AliInputEventHandler*) (man->GetInputEventHandler()); | |
324 | fPIDResponse = inputHandler->GetPIDResponse(); | |
325 | ||
326 | // Only used to get the number of primary reconstructed tracks | |
327 | if (fAnalysisType == "ESD"&& (! fESDtrackCuts )){ | |
328 | fESDtrackCuts = AliESDtrackCuts::GetStandardITSTPCTrackCuts2010(kTRUE); // Std definition of primary (see kTRUE argument) tracks for 2010 | |
329 | // fESDtrackCuts->SetEtaRange(-0.8,+0.8); | |
330 | // fESDtrackCuts->SetPtRange(0.15, 1e10); | |
331 | Printf("CheckCascade - ESDtrackCuts set up to 2010 std ITS-TPC cuts..."); | |
332 | } | |
333 | ||
334 | // Initialize cuts to re-run V0 and cascade vertexers | |
335 | fV0Sels[0] = 33. ; // max allowed chi2 | |
336 | fV0Sels[1] = 0.1 ; // min allowed impact parameter for the 1st daughter | |
337 | fV0Sels[2] = 0.1 ; // min allowed impact parameter for the 2nd daughter | |
338 | fV0Sels[3] = 1. ; // max allowed DCA between the daughter tracks | |
339 | fV0Sels[4] = .998; // min allowed cosine of V0's pointing angle | |
340 | fV0Sels[5] = 0.9 ; // min radius of the fiducial volume | |
341 | fV0Sels[6] = 100. ; // max radius of the fiducial volume | |
342 | ||
343 | fCascSels[0] = 33. ; // max allowed chi2 (not used) | |
344 | fCascSels[1] = 0.05 ; // min allowed V0 impact parameter | |
345 | fCascSels[2] = 0.008 ; // "window" around the Lambda mass | |
346 | fCascSels[3] = 0.05 ; // min allowed bachelor's impact parameter | |
347 | fCascSels[4] = 0.3 ; // max allowed DCA between the V0 and the bachelor | |
348 | fCascSels[5] = 0.9989; // min allowed cosine of the cascade pointing angle | |
349 | fCascSels[6] = 0.9 ; // min radius of the fiducial volume | |
350 | fCascSels[7] = 100. ; | |
351 | ||
f9a6cab5 | 352 | |
f9a6cab5 | 353 | /* |
354 | if( !fPaveTextBookKeeping){ | |
355 | // FIXME : prepare a field with the AliRoot+root distrib ... | |
356 | ||
357 | fPaveTextBookKeeping = new TPaveText(0.1, 0.1, 0.9, 0.9,"NDC"); | |
358 | fPaveTextBookKeeping->SetName("fPaveTextBookKeeping"); | |
359 | fPaveTextBookKeeping->SetBorderSize(0); | |
360 | fPaveTextBookKeeping->SetTextAlign(12); | |
361 | fPaveTextBookKeeping->SetFillColor(kWhite); | |
362 | fPaveTextBookKeeping->SetTextFont(42); // regular Arial or Helvetica, | |
363 | fPaveTextBookKeeping->SetTextColor(kBlue+4); | |
364 | ||
365 | fPaveTextBookKeeping->AddText( "Task CHECK CASCADE analysis" ); | |
366 | fPaveTextBookKeeping->AddText("- - - - - - - - - - - "); | |
367 | fPaveTextBookKeeping->AddText( Form("AnalysisType : %s ", fAnalysisType.Data() )); | |
39b18694 | 368 | |
f9a6cab5 | 369 | else fPaveTextBookKeeping->AddText("Colliding system : A-A collisions "); |
370 | ||
371 | fPaveTextBookKeeping->AddText("- - - - - - - - - - - "); | |
372 | ||
373 | if(fkRerunV0CascVertexers){ | |
374 | fPaveTextBookKeeping->AddText("A.1. With V0 vertexer : "); | |
375 | fPaveTextBookKeeping->AddText( Form(" - V0 #chi^{2} _________________ < %.3f ", fV0Sels[0] )); | |
376 | fPaveTextBookKeeping->AddText( Form(" - DCA(prim. Vtx/ 1^{st} daughter) ___ > %.3f cm ", fV0Sels[1] )); | |
377 | fPaveTextBookKeeping->AddText( Form(" - DCA(prim. Vtx/ 2^{nd} daughter) __ > %.3f cm", fV0Sels[2] )); | |
378 | fPaveTextBookKeeping->AddText( Form(" - DCA between V0 daughters ___ < %.3f cm", fV0Sels[3] )); | |
379 | fPaveTextBookKeeping->AddText( Form(" - cos(V0 pointing angle) _______ > %.3f ", fV0Sels[4] )); | |
380 | fPaveTextBookKeeping->AddText( Form(" - R_{transv}(V0 decay) ________ > %.3f cm", fV0Sels[5] )); | |
381 | fPaveTextBookKeeping->AddText( Form(" - R_{transv}(V0 decay) ________ < %.3f cm", fV0Sels[6] )); | |
382 | ||
383 | fPaveTextBookKeeping->AddText(" "); | |
384 | ||
385 | fPaveTextBookKeeping->AddText("A.2. With Casc. vertexer : "); | |
386 | fPaveTextBookKeeping->AddText( Form(" - Casc. #chi^{2} ______________ < %.3f ", fCascSels[0] )); | |
387 | fPaveTextBookKeeping->AddText( Form(" - DCA(prim. Vtx/ V0) _________ > %.3f cm", fCascSels[1] )); | |
388 | fPaveTextBookKeeping->AddText( Form(" - | M_{#Lambda}(reco) - M_{#Lambda}(pdg) | _______ < %.3f GeV/c^{2}", fCascSels[2] )); | |
389 | fPaveTextBookKeeping->AddText( Form(" - DCA(prim. Vtx/ Bach) _______ > %.3f cm", fCascSels[3] )); | |
390 | fPaveTextBookKeeping->AddText( Form(" - DCA between Bach/ #Lambda ______ < %.3f cm", fCascSels[4] )); | |
391 | fPaveTextBookKeeping->AddText( Form(" - cos(Casc. pointing angle) ____ > %.3f ", fCascSels[5] )); | |
392 | fPaveTextBookKeeping->AddText( Form(" - R_{transv}(Casc. decay) ______ > %.3f cm", fCascSels[6] )); | |
393 | fPaveTextBookKeeping->AddText( Form(" - R_{transv}(Casc. decay) ______ < %.3f cm", fCascSels[7] )); | |
394 | } | |
395 | else{ fPaveTextBookKeeping->AddText("A. No rerunning of the V0/Casc. vertexers ... See std cuts in (AliRoot+Rec.C) used for this prod. cycle");} | |
396 | ||
397 | fPaveTextBookKeeping->AddText("- - - - - - - - - - - "); | |
398 | ||
399 | if(fkQualityCutZprimVtxPos) fPaveTextBookKeeping->AddText("B. Quality Cut(prim. Vtx z-Pos) = ON "); | |
400 | else fPaveTextBookKeeping->AddText("B. Quality Cut(prim. Vtx z-Pos) = Off "); | |
401 | if(fkQualityCutNoTPConlyPrimVtx) fPaveTextBookKeeping->AddText("C. Quality Cut(No TPC prim. vtx) = ON "); | |
402 | else fPaveTextBookKeeping->AddText("C. Quality Cut(No TPC prim. vtx) = Off "); | |
403 | if(fkQualityCutTPCrefit) fPaveTextBookKeeping->AddText("D. Quality Cut(TPCrefit) = ON "); | |
404 | else fPaveTextBookKeeping->AddText("D. Quality Cut(TPCrefit) = Off "); | |
405 | if(fkQualityCut80TPCcls) fPaveTextBookKeeping->AddText("E. Quality Cut(80 TPC clusters) = ON "); | |
406 | else fPaveTextBookKeeping->AddText("E. Quality Cut(80 TPC clusters) = Off "); | |
407 | if(fkExtraSelections) fPaveTextBookKeeping->AddText("F. Extra Analysis Selections = ON "); | |
408 | else fPaveTextBookKeeping->AddText("F. Extra Analysis Selections = Off "); | |
409 | ||
410 | fPaveTextBookKeeping->AddText("- - - - - - - - - - - "); | |
411 | ||
f9a6cab5 | 412 | } |
413 | */ | |
414 | ||
415 | // - General histos | |
416 | //-------------- | |
11bcd1e4 | 417 | // Events in centraity bins |
418 | Double_t centBinLimits[12] = {0.,5.,10.,20.,30.,40.,50.,60.,70.,80.,90.,100.}; | |
419 | fHistEvtsInCentralityBinsvsNtracks = new TH2F("fHistEvtsInCentralityBinsvsNtracks","",11,centBinLimits,100,0.,6000.); | |
420 | fListHistCascade->Add(fHistEvtsInCentralityBinsvsNtracks); | |
f9a6cab5 | 421 | |
422 | if(! fHistCascadeMultiplicityBeforeEvSel) { | |
39b18694 | 423 | fHistCascadeMultiplicityBeforeEvSel = new TH1F("fHistCascadeMultiplicityBeforeEvSel", |
f9a6cab5 | 424 | "Cascades per event (before vertex and centrality selections);Nbr of Cascades/Evt;Events", |
425 | 100, 0, 200000); | |
f9a6cab5 | 426 | fListHistCascade->Add(fHistCascadeMultiplicityBeforeEvSel); |
427 | } | |
428 | ||
429 | // - Histos for events passing the trigger selection | |
430 | //-------------- | |
431 | ||
432 | if(! fHistCascadeMultiplicityForCentrEvt) { | |
39b18694 | 433 | fHistCascadeMultiplicityForCentrEvt = new TH1F("fHistCascadeMultiplicityForCentrEvt", |
f9a6cab5 | 434 | "Cascades per event (for triggered evt);Nbr of Cascades/Evt;Events", |
435 | 100, 0, 200000); | |
f9a6cab5 | 436 | fListHistCascade->Add(fHistCascadeMultiplicityForCentrEvt); |
437 | } | |
438 | ||
439 | ||
440 | if(! fHistTrackMultiplicityForCentrEvt) { | |
39b18694 | 441 | fHistTrackMultiplicityForCentrEvt = new TH1F("fHistTrackMultiplicityForCentrEvt", |
f9a6cab5 | 442 | "Track Multiplicity (for triggered evt);Nbr of tracks/Evt;Events", |
443 | 200, 0, 12000); | |
f9a6cab5 | 444 | fListHistCascade->Add(fHistTrackMultiplicityForCentrEvt); |
445 | } | |
446 | ||
447 | if(! fHistTPCrefitTrackMultiplicityForCentrEvt) { | |
39b18694 | 448 | fHistTPCrefitTrackMultiplicityForCentrEvt = new TH1F("fHistTPCrefitTrackMultiplicityForCentrEvt", |
f9a6cab5 | 449 | "TPCrefit track Multiplicity (for triggered evt);Nbr of TPCrefit tracks/Evt;Events", |
450 | 200, 0, 12000); | |
f9a6cab5 | 451 | fListHistCascade->Add(fHistTPCrefitTrackMultiplicityForCentrEvt); |
452 | } | |
453 | ||
f9a6cab5 | 454 | // - Histos for events passing the event selection at the analysis level |
455 | //-------------- | |
456 | ||
457 | if(! fHistCascadeMultiplicityForSelEvt) { | |
39b18694 | 458 | fHistCascadeMultiplicityForSelEvt = new TH1F("fHistCascadeMultiplicityForSelEvt", |
f9a6cab5 | 459 | "Cascades per event;Nbr of Cascades/Evt;Events", |
460 | 100, 0, 200000); | |
f9a6cab5 | 461 | fListHistCascade->Add(fHistCascadeMultiplicityForSelEvt); |
462 | } | |
463 | ||
464 | if(! fHistPosBestPrimaryVtxXForSelEvt ){ | |
465 | fHistPosBestPrimaryVtxXForSelEvt = new TH1F( "fHistPosBestPrimaryVtxXForSelEvt" , "Best Prim. Vertex Position in x; x (cm); Events" , 360, -0.9, 0.9 ); | |
466 | fListHistCascade->Add(fHistPosBestPrimaryVtxXForSelEvt); | |
467 | } | |
468 | ||
469 | if(! fHistPosBestPrimaryVtxYForSelEvt){ | |
470 | fHistPosBestPrimaryVtxYForSelEvt = new TH1F( "fHistPosBestPrimaryVtxYForSelEvt" , "Best Prim. Vertex Position in y; y (cm); Events" , 360, -0.9, 0.9 ); | |
471 | fListHistCascade->Add(fHistPosBestPrimaryVtxYForSelEvt); | |
472 | } | |
473 | ||
474 | if(! fHistPosBestPrimaryVtxZForSelEvt ){ | |
475 | fHistPosBestPrimaryVtxZForSelEvt = new TH1F( "fHistPosBestPrimaryVtxZForSelEvt" , "Best Prim. Vertex Position in z; z (cm); Events" , 300, -30.0, 30.0 ); | |
476 | fListHistCascade->Add(fHistPosBestPrimaryVtxZForSelEvt); | |
477 | } | |
478 | ||
479 | ||
480 | ||
481 | ||
482 | // - Histos for events containing at least ONE CASCADE | |
483 | //-------------- | |
484 | ||
485 | if(! fHistTPCrefitTrackMultiplicityForCascadeEvt) { | |
39b18694 | 486 | fHistTPCrefitTrackMultiplicityForCascadeEvt = new TH1F("fHistTPCrefitTrackMultiplicityForCascadeEvt", |
f9a6cab5 | 487 | "TPCrefit track Multiplicity (for evt with Casc.);Nbr of TPCrefit tracks/Evt with cascade(s);Events", |
488 | 200, 0, 20000); | |
f9a6cab5 | 489 | fListHistCascade->Add(fHistTPCrefitTrackMultiplicityForCascadeEvt); |
490 | } | |
491 | ||
492 | if(! fHistPosV0TPCClusters ){ | |
493 | fHistPosV0TPCClusters = new TH1F("fHistPosV0TPCClusters", "TPC clusters for Pos. V0 daughter track, in Casc; Nbr of TPC clusters (V0 Pos.); Track counts", 165, 0.0 ,165.0); | |
494 | fListHistCascade->Add(fHistPosV0TPCClusters); | |
495 | } | |
496 | ||
497 | if(! fHistNegV0TPCClusters ){ | |
498 | fHistNegV0TPCClusters = new TH1F("fHistNegV0TPCClusters", "TPC clusters for Neg. V0 daughter track, in Casc; Nbr of TPC clusters (V0 Neg.); Track counts", 165, 0.0 ,165.0); | |
499 | fListHistCascade->Add(fHistNegV0TPCClusters); | |
500 | } | |
501 | ||
502 | if(! fHistBachTPCClusters ){ | |
503 | fHistBachTPCClusters = new TH1F("fHistBachTPCClusters", "TPC clusters for Bachelor track; Nbr of TPC clusters (Bach); Track counts", 165, 0.0 ,165.0); | |
504 | fListHistCascade->Add(fHistBachTPCClusters); | |
505 | } | |
506 | ||
507 | ||
508 | ||
509 | if(! fHistVtxStatus ){ | |
510 | fHistVtxStatus = new TH1F( "fHistVtxStatus" , "Does a Trckg Prim.vtx exist ?; true=1 or false=0; Nb of Events" , 4, -1.0, 3.0 ); | |
511 | fListHistCascade->Add(fHistVtxStatus); | |
512 | } | |
513 | ||
514 | ||
515 | // - Vertex Positions | |
516 | ||
517 | if(! fHistPosTrkgPrimaryVtxXForCascadeEvt ){ | |
518 | fHistPosTrkgPrimaryVtxXForCascadeEvt = new TH1F( "fHistPosTrkgPrimaryVtxXForCascadeEvt" , "Trkg Prim. Vertex Position in x; x (cm); Events" , 120, -0.6, 0.6 ); | |
519 | fListHistCascade->Add(fHistPosTrkgPrimaryVtxXForCascadeEvt); | |
520 | } | |
521 | ||
522 | ||
523 | if(! fHistPosTrkgPrimaryVtxYForCascadeEvt){ | |
524 | fHistPosTrkgPrimaryVtxYForCascadeEvt = new TH1F( "fHistPosTrkgPrimaryVtxYForCascadeEvt" , "Trkg Prim. Vertex Position in y; y (cm); Events" , 120, -0.6, 0.6 ); | |
525 | fListHistCascade->Add(fHistPosTrkgPrimaryVtxYForCascadeEvt); | |
526 | } | |
527 | ||
528 | if(! fHistPosTrkgPrimaryVtxZForCascadeEvt ){ | |
529 | fHistPosTrkgPrimaryVtxZForCascadeEvt = new TH1F( "fHistPosTrkgPrimaryVtxZForCascadeEvt" , "Trkg Prim. Vertex Position in z; z (cm); Events" , 200, -20.0, 20.0 ); | |
530 | fListHistCascade->Add(fHistPosTrkgPrimaryVtxZForCascadeEvt); | |
531 | } | |
532 | ||
533 | if(! fHistTrkgPrimaryVtxRadius ){ | |
534 | fHistTrkgPrimaryVtxRadius = new TH1F( "fHistTrkgPrimaryVtxRadius", "Trkg Prim. Vertex radius; r (cm); Events" , 150, 0., 15.0 ); | |
535 | fListHistCascade->Add(fHistTrkgPrimaryVtxRadius); | |
536 | } | |
537 | ||
538 | ||
539 | ||
540 | ||
541 | if(! fHistPosBestPrimaryVtxXForCascadeEvt ){ | |
542 | fHistPosBestPrimaryVtxXForCascadeEvt = new TH1F( "fHistPosBestPrimaryVtxXForCascadeEvt" , "Best Prim. Vertex Position in x; x (cm); Events" , 120, -0.6, 0.6 ); | |
543 | fListHistCascade->Add(fHistPosBestPrimaryVtxXForCascadeEvt); | |
544 | } | |
545 | ||
546 | if(! fHistPosBestPrimaryVtxYForCascadeEvt){ | |
547 | fHistPosBestPrimaryVtxYForCascadeEvt = new TH1F( "fHistPosBestPrimaryVtxYForCascadeEvt" , "Best Prim. Vertex Position in y; y (cm); Events" , 120, -0.6, 0.6 ); | |
548 | fListHistCascade->Add(fHistPosBestPrimaryVtxYForCascadeEvt); | |
549 | } | |
550 | ||
551 | if(! fHistPosBestPrimaryVtxZForCascadeEvt ){ | |
552 | fHistPosBestPrimaryVtxZForCascadeEvt = new TH1F( "fHistPosBestPrimaryVtxZForCascadeEvt" , "Best Prim. Vertex Position in z; z (cm); Events" , 200, -20.0, 20.0 ); | |
553 | fListHistCascade->Add(fHistPosBestPrimaryVtxZForCascadeEvt); | |
554 | } | |
555 | ||
556 | if(! fHistBestPrimaryVtxRadius ){ | |
557 | fHistBestPrimaryVtxRadius = new TH1F( "fHistBestPrimaryVtxRadius", "Best Prim. vertex radius; r (cm); Events" , 150, 0., 15.0 ); | |
558 | fListHistCascade->Add(fHistBestPrimaryVtxRadius); | |
559 | } | |
560 | ||
561 | if(! f2dHistTrkgPrimVtxVsBestPrimVtx) { | |
562 | f2dHistTrkgPrimVtxVsBestPrimVtx = new TH2F( "f2dHistTrkgPrimVtxVsBestPrimVtx", "r_{Trck Prim. Vtx} Vs r_{Best Prim. Vtx}; r_{Track Vtx} (cm); r_{Best Vtx} (cm)", 150, 0., 15.0, 150, 0., 15.); | |
563 | fListHistCascade->Add(f2dHistTrkgPrimVtxVsBestPrimVtx); | |
564 | } | |
565 | ||
566 | ||
567 | ||
568 | ||
569 | // - Typical histos for cascades | |
570 | ||
571 | ||
572 | if(! fHistEffMassXi) { | |
573 | fHistEffMassXi = new TH1F("fHistEffMassXi", "Cascade candidates ; Invariant Mass (GeV/c^{2}) ; Counts", 400, 1.2, 2.0); | |
574 | fListHistCascade->Add(fHistEffMassXi); | |
575 | } | |
576 | ||
577 | if(! fHistChi2Xi ){ | |
578 | fHistChi2Xi = new TH1F("fHistChi2Xi", "Cascade #chi^{2}; #chi^{2}; Number of Cascades", 160, 0, 40); | |
579 | fListHistCascade->Add(fHistChi2Xi); | |
580 | } | |
581 | ||
582 | if(! fHistDcaXiDaughters ){ | |
583 | fHistDcaXiDaughters = new TH1F( "fHistDcaXiDaughters", "DCA between Xi Daughters; DCA (cm) ; Number of Cascades", 100, 0., 0.5); | |
584 | fListHistCascade->Add(fHistDcaXiDaughters); | |
585 | } | |
586 | ||
587 | if(! fHistDcaBachToPrimVertex) { | |
588 | fHistDcaBachToPrimVertex = new TH1F("fHistDcaBachToPrimVertex", "DCA of Bach. to Prim. Vertex;DCA (cm);Number of Cascades", 250, 0., 0.25); | |
589 | fListHistCascade->Add(fHistDcaBachToPrimVertex); | |
590 | } | |
591 | ||
592 | if(! fHistXiCosineOfPointingAngle) { | |
39b18694 | 593 | fHistXiCosineOfPointingAngle = new TH1F("fHistXiCosineOfPointingAngle", "Cosine of Xi Pointing Angle; Cos (Xi Point.Angl);Number of Xis", 101, 0.99, 1.0001); |
f9a6cab5 | 594 | fListHistCascade->Add(fHistXiCosineOfPointingAngle); |
595 | } | |
596 | ||
597 | if(! fHistXiRadius ){ | |
598 | fHistXiRadius = new TH1F( "fHistXiRadius", "Casc. decay transv. radius; r (cm); Counts" , 1050, 0., 105.0 ); | |
599 | fListHistCascade->Add(fHistXiRadius); | |
600 | } | |
601 | ||
602 | ||
603 | // - Histos about ~ the "V0 part" of the cascade, coming by inheritance from AliESDv0 | |
604 | ||
605 | ||
606 | ||
607 | if (! fHistMassLambdaAsCascDghter) { | |
608 | fHistMassLambdaAsCascDghter = new TH1F("fHistMassLambdaAsCascDghter","#Lambda associated to Casc. candidates;Eff. Mass (GeV/c^{2});Counts", 300,1.00,1.3); | |
609 | fListHistCascade->Add(fHistMassLambdaAsCascDghter); | |
610 | } | |
611 | ||
612 | if (! fHistV0Chi2Xi) { | |
613 | fHistV0Chi2Xi = new TH1F("fHistV0Chi2Xi", "V0 #chi^{2}, in cascade; #chi^{2};Counts", 160, 0, 40); | |
614 | fListHistCascade->Add(fHistV0Chi2Xi); | |
615 | } | |
616 | ||
617 | if (! fHistDcaV0DaughtersXi) { | |
39b18694 | 618 | fHistDcaV0DaughtersXi = new TH1F("fHistDcaV0DaughtersXi", "DCA between V0 daughters, in cascade;DCA (cm);Number of V0s", 240, 0., 1.2); |
f9a6cab5 | 619 | fListHistCascade->Add(fHistDcaV0DaughtersXi); |
620 | } | |
621 | ||
622 | if (! fHistDcaV0ToPrimVertexXi) { | |
623 | fHistDcaV0ToPrimVertexXi = new TH1F("fHistDcaV0ToPrimVertexXi", "DCA of V0 to Prim. Vertex, in cascade;DCA (cm);Number of Cascades", 200, 0., 1.); | |
624 | fListHistCascade->Add(fHistDcaV0ToPrimVertexXi); | |
625 | } | |
626 | ||
627 | if (! fHistV0CosineOfPointingAngleXi) { | |
39b18694 | 628 | fHistV0CosineOfPointingAngleXi = new TH1F("fHistV0CosineOfPointingAngleXi", "Cosine of V0 Pointing Angle, in cascade;Cos(V0 Point. Angl); Counts", 200, 0.95, 1.0001); |
f9a6cab5 | 629 | fListHistCascade->Add(fHistV0CosineOfPointingAngleXi); |
630 | } | |
631 | ||
632 | if (! fHistV0RadiusXi) { | |
633 | fHistV0RadiusXi = new TH1F("fHistV0RadiusXi", "V0 decay radius, in cascade; radius (cm); Counts", 1050, 0., 105.0); | |
634 | fListHistCascade->Add(fHistV0RadiusXi); | |
635 | } | |
636 | ||
637 | if (! fHistDcaPosToPrimVertexXi) { | |
638 | fHistDcaPosToPrimVertexXi = new TH1F("fHistDcaPosToPrimVertexXi", "DCA of V0 pos daughter to Prim. Vertex;DCA (cm);Counts", 300, 0, 3); | |
639 | fListHistCascade->Add(fHistDcaPosToPrimVertexXi); | |
640 | } | |
641 | ||
642 | if (! fHistDcaNegToPrimVertexXi) { | |
643 | fHistDcaNegToPrimVertexXi = new TH1F("fHistDcaNegToPrimVertexXi", "DCA of V0 neg daughter to Prim. Vertex;DCA (cm);Counts", 300, 0, 3); | |
644 | fListHistCascade->Add(fHistDcaNegToPrimVertexXi); | |
645 | } | |
646 | ||
647 | ||
648 | ||
649 | ||
650 | // - Effective mass histos for cascades. | |
651 | // By cascade hyp | |
652 | if (! fHistMassXiMinus) { | |
653 | fHistMassXiMinus = new TH1F("fHistMassXiMinus","#Xi^{-} candidates;M( #Lambda , #pi^{-} ) (GeV/c^{2});Counts", 400,1.2,2.0); | |
654 | fListHistCascade->Add(fHistMassXiMinus); | |
655 | } | |
656 | ||
657 | if (! fHistMassXiPlus) { | |
658 | fHistMassXiPlus = new TH1F("fHistMassXiPlus","#Xi^{+} candidates;M( #bar{#Lambda}^{0} , #pi^{+} ) (GeV/c^{2});Counts",400,1.2,2.0); | |
659 | fListHistCascade->Add(fHistMassXiPlus); | |
660 | } | |
661 | ||
662 | if (! fHistMassOmegaMinus) { | |
663 | fHistMassOmegaMinus = new TH1F("fHistMassOmegaMinus","#Omega^{-} candidates;M( #Lambda , K^{-} ) (GeV/c^{2});Counts", 500,1.5,2.5); | |
664 | fListHistCascade->Add(fHistMassOmegaMinus); | |
665 | } | |
666 | ||
667 | if (! fHistMassOmegaPlus) { | |
668 | fHistMassOmegaPlus = new TH1F("fHistMassOmegaPlus","#Omega^{+} candidates;M( #bar{#Lambda}^{0} , K^{+} ) (GeV/c^{2});Counts", 500,1.5,2.5); | |
669 | fListHistCascade->Add(fHistMassOmegaPlus); | |
670 | } | |
671 | ||
672 | // By cascade hyp + bachelor PID | |
673 | if (! fHistMassWithCombPIDXiMinus) { | |
674 | fHistMassWithCombPIDXiMinus = new TH1F("fHistMassWithCombPIDXiMinus","#Xi^{-} candidates, with Bach. comb. PID;M( #Lambda , #pi^{-} ) (GeV/c^{2});Counts", 400,1.2,2.0); | |
675 | fListHistCascade->Add(fHistMassWithCombPIDXiMinus); | |
676 | } | |
677 | ||
678 | if (! fHistMassWithCombPIDXiPlus) { | |
679 | fHistMassWithCombPIDXiPlus = new TH1F("fHistMassWithCombPIDXiPlus","#Xi^{+} candidates, with Bach. comb. PID;M( #bar{#Lambda}^{0} , #pi^{+} ) (GeV/c^{2});Counts",400,1.2,2.0); | |
680 | fListHistCascade->Add(fHistMassWithCombPIDXiPlus); | |
681 | } | |
682 | ||
683 | if (! fHistMassWithCombPIDOmegaMinus) { | |
684 | fHistMassWithCombPIDOmegaMinus = new TH1F("fHistMassWithCombPIDOmegaMinus","#Omega^{-} candidates, with Bach. comb. PID;M( #Lambda , K^{-} ) (GeV/c^{2});Counts", 500,1.5,2.5); | |
685 | fListHistCascade->Add(fHistMassWithCombPIDOmegaMinus); | |
686 | } | |
687 | ||
688 | if (! fHistMassWithCombPIDOmegaPlus) { | |
689 | fHistMassWithCombPIDOmegaPlus = new TH1F("fHistMassWithCombPIDOmegaPlus","#Omega^{+} candidates, with Bach. comb. PID;M( #bar{#Lambda}^{0} , K^{+} ) (GeV/c^{2});Counts", 500,1.5,2.5); | |
690 | fListHistCascade->Add(fHistMassWithCombPIDOmegaPlus); | |
691 | } | |
692 | ||
693 | ||
694 | ||
695 | // - Complements for QA | |
696 | ||
697 | if(! fHistXiTransvMom ){ | |
698 | fHistXiTransvMom = new TH1F( "fHistXiTransvMom" , "#Xi transverse momentum (cand. around the mass peak) ; p_{t}(#Xi) (GeV/c); Counts", 100, 0.0, 10.0); | |
699 | fListHistCascade->Add(fHistXiTransvMom); | |
700 | } | |
701 | ||
702 | if(! fHistXiTotMom ){ | |
703 | fHistXiTotMom = new TH1F( "fHistXiTotMom" , "#Xi momentum norm (cand. around the mass peak); p_{tot}(#Xi) (GeV/c); Counts", 150, 0.0, 15.0); | |
704 | fListHistCascade->Add(fHistXiTotMom); | |
705 | } | |
706 | ||
707 | ||
708 | if(! fHistBachTransvMomXi ){ | |
709 | fHistBachTransvMomXi = new TH1F( "fHistBachTransvMomXi" , "#Xi Bach. transverse momentum (cand. around the mass peak) ; p_{t}(Bach.) (GeV/c); Counts", 100, 0.0, 5.0); | |
710 | fListHistCascade->Add(fHistBachTransvMomXi); | |
711 | } | |
712 | ||
713 | if(! fHistBachTotMomXi ){ | |
714 | fHistBachTotMomXi = new TH1F( "fHistBachTotMomXi" , "#Xi Bach. momentum norm (cand. around the mass peak); p_{tot}(Bach.) (GeV/c); Counts", 100, 0.0, 5.0); | |
715 | fListHistCascade->Add(fHistBachTotMomXi); | |
716 | } | |
717 | ||
718 | ||
719 | if(! fHistChargeXi ){ | |
720 | fHistChargeXi = new TH1F( "fHistChargeXi" , "Charge of casc. candidates ; Sign ; Counts", 5, -2.0, 3.0); | |
721 | fListHistCascade->Add(fHistChargeXi); | |
722 | } | |
723 | ||
724 | ||
725 | if (! fHistV0toXiCosineOfPointingAngle) { | |
39b18694 | 726 | fHistV0toXiCosineOfPointingAngle = new TH1F("fHistV0toXiCosineOfPointingAngle", "Cos. of V0 Ptng Angl / Xi vtx ;Cos(V0 Point. Angl / Xi vtx); Counts", 200, 0.95, 1.0001); |
f9a6cab5 | 727 | fListHistCascade->Add(fHistV0toXiCosineOfPointingAngle); |
728 | } | |
729 | ||
730 | ||
731 | if(! fHistRapXi ){ | |
732 | fHistRapXi = new TH1F( "fHistRapXi" , "Rapidity of #Xi candidates (around the mass peak); y ; Counts", 20, -1.0, 1.0); | |
733 | fListHistCascade->Add(fHistRapXi); | |
734 | } | |
735 | ||
736 | if(! fHistRapOmega ){ | |
737 | fHistRapOmega = new TH1F( "fHistRapOmega" , "Rapidity of #Omega candidates (around the mass peak); y ; Counts", 20, -1.0, 1.0); | |
738 | fListHistCascade->Add(fHistRapOmega); | |
739 | } | |
740 | ||
741 | if(! fHistEtaXi ){ | |
742 | fHistEtaXi = new TH1F( "fHistEtaXi" , "Pseudo-rap. of #Xi candidates (around the mass peak) ; #eta ; Counts", 20, -1.0, 1.0); | |
743 | fListHistCascade->Add(fHistEtaXi); | |
744 | } | |
745 | ||
746 | if(! fHistThetaXi ){ | |
747 | fHistThetaXi = new TH1F( "fHistThetaXi" , "#theta of #Xi candidates (around the mass peak); #theta (deg) ; Counts", 180, 0., 180.0); | |
748 | fListHistCascade->Add(fHistThetaXi); | |
749 | } | |
750 | ||
751 | if(! fHistPhiXi ){ | |
752 | fHistPhiXi = new TH1F( "fHistPhiXi" , "#phi of #Xi candidates (around the mass peak); #phi (deg) ; Counts", 360, 0., 360.); | |
753 | fListHistCascade->Add(fHistPhiXi); | |
754 | } | |
755 | ||
756 | ||
757 | if(! f2dHistArmenteros) { | |
758 | f2dHistArmenteros = new TH2F( "f2dHistArmenteros", "#alpha_{Arm}(casc. cand.) Vs Pt_{Arm}(casc. cand.); #alpha_{Arm} ; Pt_{Arm} (GeV/c)", 140, -1.2, 1.2, 300, 0., 0.3); | |
759 | fListHistCascade->Add(f2dHistArmenteros); | |
760 | } | |
761 | ||
762 | //------- | |
763 | ||
764 | if(! f2dHistEffMassLambdaVsEffMassXiMinus) { | |
765 | f2dHistEffMassLambdaVsEffMassXiMinus = new TH2F( "f2dHistEffMassLambdaVsEffMassXiMinus", "M_{#Lambda} Vs M_{#Xi^{-} candidates} ; Inv. M_{#Lambda^{0}} (GeV/c^{2}) ; M( #Lambda , #pi^{-} ) (GeV/c^{2})", 300, 1.1,1.13, 400, 1.2, 2.0); | |
766 | fListHistCascade->Add(f2dHistEffMassLambdaVsEffMassXiMinus); | |
767 | } | |
768 | ||
769 | if(! f2dHistEffMassXiVsEffMassOmegaMinus) { | |
770 | f2dHistEffMassXiVsEffMassOmegaMinus = new TH2F( "f2dHistEffMassXiVsEffMassOmegaMinus", "M_{#Xi^{-} candidates} Vs M_{#Omega^{-} candidates} ; M( #Lambda , #pi^{-} ) (GeV/c^{2}) ; M( #Lambda , K^{-} ) (GeV/c^{2})", 400, 1.2, 2.0, 500, 1.5, 2.5); | |
771 | fListHistCascade->Add(f2dHistEffMassXiVsEffMassOmegaMinus); | |
772 | } | |
773 | ||
774 | if(! f2dHistEffMassLambdaVsEffMassXiPlus) { | |
775 | f2dHistEffMassLambdaVsEffMassXiPlus = new TH2F( "f2dHistEffMassLambdaVsEffMassXiPlus", "M_{#Lambda} Vs M_{#Xi^{+} candidates} ; Inv. M_{#Lambda^{0}} (GeV/c^{2}) ; M( #Lambda , #pi^{+} ) (GeV/c^{2})", 300, 1.1,1.13, 400, 1.2, 2.0); | |
776 | fListHistCascade->Add(f2dHistEffMassLambdaVsEffMassXiPlus); | |
777 | } | |
778 | ||
779 | if(! f2dHistEffMassXiVsEffMassOmegaPlus) { | |
780 | f2dHistEffMassXiVsEffMassOmegaPlus = new TH2F( "f2dHistEffMassXiVsEffMassOmegaPlus", "M_{#Xi^{+} candidates} Vs M_{#Omega^{+} candidates} ; M( #Lambda , #pi^{+} ) (GeV/c^{2}) ; M( #Lambda , K^{+} ) (GeV/c^{2})", 400, 1.2, 2.0, 500, 1.5, 2.5); | |
781 | fListHistCascade->Add(f2dHistEffMassXiVsEffMassOmegaPlus); | |
782 | } | |
783 | ||
784 | //------- | |
785 | ||
786 | if(! f2dHistXiRadiusVsEffMassXiMinus) { | |
787 | f2dHistXiRadiusVsEffMassXiMinus = new TH2F( "f2dHistXiRadiusVsEffMassXiMinus", "Transv. R_{Xi Decay} Vs M_{#Xi^{-} candidates}; r_{cascade} (cm); M( #Lambda , #pi^{-} ) (GeV/c^{2}) ", 450, 0., 45.0, 400, 1.2, 2.0); | |
788 | fListHistCascade->Add(f2dHistXiRadiusVsEffMassXiMinus); | |
789 | } | |
790 | ||
791 | if(! f2dHistXiRadiusVsEffMassXiPlus) { | |
792 | f2dHistXiRadiusVsEffMassXiPlus = new TH2F( "f2dHistXiRadiusVsEffMassXiPlus", "Transv. R_{Xi Decay} Vs M_{#Xi^{+} candidates}; r_{cascade} (cm); M( #Lambda , #pi^{+} ) (GeV/c^{2}) ", 450, 0., 45.0, 400, 1.2, 2.0); | |
793 | fListHistCascade->Add(f2dHistXiRadiusVsEffMassXiPlus); | |
794 | } | |
795 | ||
796 | if(! f2dHistXiRadiusVsEffMassOmegaMinus) { | |
797 | f2dHistXiRadiusVsEffMassOmegaMinus = new TH2F( "f2dHistXiRadiusVsEffMassOmegaMinus", "Transv. R_{Xi Decay} Vs M_{#Omega^{-} candidates}; r_{cascade} (cm); M( #Lambda , K^{-} ) (GeV/c^{2}) ", 450, 0., 45.0, 500, 1.5, 2.5); | |
798 | fListHistCascade->Add(f2dHistXiRadiusVsEffMassOmegaMinus); | |
799 | } | |
800 | ||
801 | if(! f2dHistXiRadiusVsEffMassOmegaPlus) { | |
802 | f2dHistXiRadiusVsEffMassOmegaPlus = new TH2F( "f2dHistXiRadiusVsEffMassOmegaPlus", "Transv. R_{Xi Decay} Vs M_{#Omega^{+} candidates}; r_{cascade} (cm); M( #Lambda , K^{+} ) (GeV/c^{2}) ", 450, 0., 45.0, 500, 1.5, 2.5); | |
803 | fListHistCascade->Add(f2dHistXiRadiusVsEffMassOmegaPlus); | |
804 | } | |
805 | ||
f9a6cab5 | 806 | if(! f2dHistTPCdEdxOfCascDghters){ |
807 | f2dHistTPCdEdxOfCascDghters = new TH2F( "f2dHistTPCdEdxOfCascDghters", "TPC dE/dx of the cascade daughters; charge x || #vec{p}_{TPC inner wall}(Casc. daughter) || (GeV/c); TPC signal (ADC) ", 2000, -10.0, 10.0, 450, 0., 900.); | |
808 | fListHistCascade->Add(f2dHistTPCdEdxOfCascDghters); | |
809 | } | |
810 | ||
811 | ||
f9a6cab5 | 812 | if(!fCFContCascadePIDXiMinus) { |
813 | const Int_t lNbSteps = 7 ; | |
814 | const Int_t lNbVariables = 4 ; | |
815 | ||
816 | //array for the number of bins in each dimension : | |
817 | Int_t lNbBinsPerVar[4] = {0}; | |
818 | lNbBinsPerVar[0] = 100; | |
11bcd1e4 | 819 | lNbBinsPerVar[1] = 800; |
820 | lNbBinsPerVar[2] = 22; | |
39b18694 | 821 | lNbBinsPerVar[3] = 11; |
11bcd1e4 | 822 | |
f9a6cab5 | 823 | |
824 | fCFContCascadePIDXiMinus = new AliCFContainer("fCFContCascadePIDXiMinus","Pt_{cascade} Vs M_{#Xi^{-} candidates} Vs Y_{#Xi}", lNbSteps, lNbVariables, lNbBinsPerVar ); | |
825 | ||
11bcd1e4 | 826 | //setting the bin limits |
f9a6cab5 | 827 | fCFContCascadePIDXiMinus->SetBinLimits(0, 0.0 , 10.0 ); // Pt(Cascade) |
828 | fCFContCascadePIDXiMinus->SetBinLimits(1, 1.2 , 2.0 ); // Xi Effective mass | |
829 | fCFContCascadePIDXiMinus->SetBinLimits(2, -1.1 , 1.1 ); // Rapidity | |
39b18694 | 830 | |
831 | Double_t *lBinLim3 = new Double_t[ lNbBinsPerVar[3]+1 ]; | |
832 | for(Int_t i=3; i< lNbBinsPerVar[3]+1;i++) lBinLim3[i] = (Double_t)(i-1)*10.; | |
833 | lBinLim3[0] = 0.0; | |
834 | lBinLim3[1] = 5.0; | |
835 | lBinLim3[2] = 10.0; | |
836 | fCFContCascadePIDXiMinus->SetBinLimits(3, lBinLim3 ); // Centrality | |
f9a6cab5 | 837 | |
838 | // Setting the step title : one per PID case | |
839 | fCFContCascadePIDXiMinus->SetStepTitle(0, "No PID"); | |
840 | fCFContCascadePIDXiMinus->SetStepTitle(1, "TPC PID / 4-#sigma cut on Bachelor track"); | |
841 | fCFContCascadePIDXiMinus->SetStepTitle(2, "TPC PID / 4-#sigma cut on Bachelor+Baryon tracks"); | |
842 | fCFContCascadePIDXiMinus->SetStepTitle(3, "TPC PID / 4-#sigma cut on Bachelor+Baryon+Meson tracks"); | |
843 | fCFContCascadePIDXiMinus->SetStepTitle(4, "Comb. PID / Bachelor"); | |
844 | fCFContCascadePIDXiMinus->SetStepTitle(5, "Comb. PID / Bachelor+Baryon"); | |
845 | fCFContCascadePIDXiMinus->SetStepTitle(6, "Comb. PID / Bachelor+Baryon+Meson"); | |
846 | ||
847 | // Setting the variable title, per axis | |
848 | fCFContCascadePIDXiMinus->SetVarTitle(0, "Pt_{cascade} (GeV/c)"); | |
849 | fCFContCascadePIDXiMinus->SetVarTitle(1, "M( #Lambda , #pi^{-} ) (GeV/c^{2})"); | |
850 | fCFContCascadePIDXiMinus->SetVarTitle(2, "Y_{#Xi}"); | |
39b18694 | 851 | fCFContCascadePIDXiMinus->SetVarTitle(3, "Centrality"); |
f9a6cab5 | 852 | |
853 | fListHistCascade->Add(fCFContCascadePIDXiMinus); | |
854 | ||
855 | } | |
856 | ||
857 | if (!fCFContCascadePIDXiPlus) { | |
858 | const Int_t lNbSteps = 7 ; | |
859 | const Int_t lNbVariables = 4 ; | |
860 | ||
861 | //array for the number of bins in each dimension : | |
862 | Int_t lNbBinsPerVar[4] = {0}; | |
863 | lNbBinsPerVar[0] = 100; | |
11bcd1e4 | 864 | lNbBinsPerVar[1] = 800; |
865 | lNbBinsPerVar[2] = 22; | |
39b18694 | 866 | lNbBinsPerVar[3] = 11; |
f9a6cab5 | 867 | |
868 | fCFContCascadePIDXiPlus = new AliCFContainer("fCFContCascadePIDXiPlus","Pt_{cascade} Vs M_{#Xi^{+} candidates} Vs Y_{#Xi}", lNbSteps, lNbVariables, lNbBinsPerVar ); | |
869 | ||
870 | ||
11bcd1e4 | 871 | //setting the bin limits |
f9a6cab5 | 872 | fCFContCascadePIDXiPlus->SetBinLimits(0, 0.0 , 10.0 ); // Pt(Cascade) |
873 | fCFContCascadePIDXiPlus->SetBinLimits(1, 1.2 , 2.0 ); // Xi Effective mass | |
874 | fCFContCascadePIDXiPlus->SetBinLimits(2, -1.1 , 1.1 ); // Rapidity | |
39b18694 | 875 | Double_t *lBinLim3 = new Double_t[ lNbBinsPerVar[3]+1 ]; |
876 | for(Int_t i=3; i< lNbBinsPerVar[3]+1;i++) lBinLim3[i] = (Double_t)(i-1)*10.; | |
877 | lBinLim3[0] = 0.0; | |
878 | lBinLim3[1] = 5.0; | |
879 | lBinLim3[2] = 10.0; | |
880 | fCFContCascadePIDXiPlus->SetBinLimits(3,lBinLim3); // Centrality | |
f9a6cab5 | 881 | |
882 | // Setting the step title : one per PID case | |
883 | fCFContCascadePIDXiPlus->SetStepTitle(0, "No PID"); | |
884 | fCFContCascadePIDXiPlus->SetStepTitle(1, "TPC PID / 4-#sigma cut on Bachelor track"); | |
885 | fCFContCascadePIDXiPlus->SetStepTitle(2, "TPC PID / 4-#sigma cut on Bachelor+Baryon tracks"); | |
886 | fCFContCascadePIDXiPlus->SetStepTitle(3, "TPC PID / 4-#sigma cut on Bachelor+Baryon+Meson tracks"); | |
887 | fCFContCascadePIDXiPlus->SetStepTitle(4, "Comb. PID / Bachelor"); | |
888 | fCFContCascadePIDXiPlus->SetStepTitle(5, "Comb. PID / Bachelor+Baryon"); | |
889 | fCFContCascadePIDXiPlus->SetStepTitle(6, "Comb. PID / Bachelor+Baryon+Meson"); | |
890 | ||
891 | // Setting the variable title, per axis | |
892 | fCFContCascadePIDXiPlus->SetVarTitle(0, "Pt_{cascade} (GeV/c)"); | |
893 | fCFContCascadePIDXiPlus->SetVarTitle(1, "M( #Lambda , #pi^{+} ) (GeV/c^{2})"); | |
894 | fCFContCascadePIDXiPlus->SetVarTitle(2, "Y_{#Xi}"); | |
39b18694 | 895 | fCFContCascadePIDXiPlus->SetVarTitle(3, "Centrality"); |
f9a6cab5 | 896 | |
897 | fListHistCascade->Add(fCFContCascadePIDXiPlus); | |
898 | ||
899 | } | |
900 | ||
901 | ||
902 | if(!fCFContCascadePIDOmegaMinus) { | |
903 | const Int_t lNbSteps = 7 ; | |
904 | const Int_t lNbVariables = 4 ; | |
905 | ||
906 | //array for the number of bins in each dimension : | |
907 | Int_t lNbBinsPerVar[4] = {0}; | |
908 | lNbBinsPerVar[0] = 100; | |
11bcd1e4 | 909 | lNbBinsPerVar[1] = 1000; |
910 | lNbBinsPerVar[2] = 22; | |
39b18694 | 911 | lNbBinsPerVar[3] = 11; |
11bcd1e4 | 912 | |
f9a6cab5 | 913 | |
914 | fCFContCascadePIDOmegaMinus = new AliCFContainer("fCFContCascadePIDOmegaMinus","Pt_{cascade} Vs M_{#Omega^{-} candidates} Vs Y_{#Omega}", lNbSteps, lNbVariables, lNbBinsPerVar ); | |
915 | ||
916 | ||
11bcd1e4 | 917 | //setting the bin limits |
f9a6cab5 | 918 | fCFContCascadePIDOmegaMinus->SetBinLimits(0, 0.0 , 10.0 ); // Pt(Cascade) |
919 | fCFContCascadePIDOmegaMinus->SetBinLimits(1, 1.5 , 2.5 ); // Omega Effective mass | |
920 | fCFContCascadePIDOmegaMinus->SetBinLimits(2, -1.1 , 1.1 ); // Rapidity | |
39b18694 | 921 | Double_t *lBinLim3 = new Double_t[ lNbBinsPerVar[3]+1 ]; |
922 | for(Int_t i=3; i< lNbBinsPerVar[3]+1;i++) lBinLim3[i] = (Double_t)(i-1)*10.; | |
923 | lBinLim3[0] = 0.0; | |
924 | lBinLim3[1] = 5.0; | |
925 | lBinLim3[2] = 10.0; | |
926 | fCFContCascadePIDOmegaMinus->SetBinLimits(3,lBinLim3); // Centrality | |
f9a6cab5 | 927 | // Setting the step title : one per PID case |
928 | fCFContCascadePIDOmegaMinus->SetStepTitle(0, "No PID"); | |
929 | fCFContCascadePIDOmegaMinus->SetStepTitle(1, "TPC PID / 4-#sigma cut on Bachelor track"); | |
930 | fCFContCascadePIDOmegaMinus->SetStepTitle(2, "TPC PID / 4-#sigma cut on Bachelor+Baryon tracks"); | |
931 | fCFContCascadePIDOmegaMinus->SetStepTitle(3, "TPC PID / 4-#sigma cut on Bachelor+Baryon+Meson tracks"); | |
932 | fCFContCascadePIDOmegaMinus->SetStepTitle(4, "Comb. PID / Bachelor"); | |
933 | fCFContCascadePIDOmegaMinus->SetStepTitle(5, "Comb. PID / Bachelor+Baryon"); | |
934 | fCFContCascadePIDOmegaMinus->SetStepTitle(6, "Comb. PID / Bachelor+Baryon+Meson"); | |
935 | ||
936 | // Setting the variable title, per axis | |
937 | fCFContCascadePIDOmegaMinus->SetVarTitle(0, "Pt_{cascade} (GeV/c)"); | |
938 | fCFContCascadePIDOmegaMinus->SetVarTitle(1, "M( #Lambda , K^{-} ) (GeV/c^{2})"); | |
939 | fCFContCascadePIDOmegaMinus->SetVarTitle(2, "Y_{#Omega}"); | |
39b18694 | 940 | fCFContCascadePIDOmegaMinus->SetVarTitle(3, "Centrality"); |
f9a6cab5 | 941 | |
942 | fListHistCascade->Add(fCFContCascadePIDOmegaMinus); | |
943 | ||
944 | } | |
945 | ||
946 | if(!fCFContCascadePIDOmegaPlus) { | |
947 | const Int_t lNbSteps = 7 ; | |
948 | const Int_t lNbVariables = 4 ; | |
949 | ||
950 | //array for the number of bins in each dimension : | |
951 | Int_t lNbBinsPerVar[4] = {0}; | |
952 | lNbBinsPerVar[0] = 100; | |
11bcd1e4 | 953 | lNbBinsPerVar[1] = 1000; |
954 | lNbBinsPerVar[2] = 22; | |
955 | lNbBinsPerVar[3] = 100; | |
f9a6cab5 | 956 | |
957 | ||
958 | fCFContCascadePIDOmegaPlus = new AliCFContainer("fCFContCascadePIDOmegaPlus","Pt_{cascade} Vs M_{#Omega^{+} candidates} Vs Y_{#Omega}", lNbSteps, lNbVariables, lNbBinsPerVar ); | |
959 | ||
960 | ||
11bcd1e4 | 961 | //setting the bin limits |
f9a6cab5 | 962 | fCFContCascadePIDOmegaPlus->SetBinLimits(0, 0.0 , 10.0 ); // Pt(Cascade) |
963 | fCFContCascadePIDOmegaPlus->SetBinLimits(1, 1.5 , 2.5 ); // Omega Effective mass | |
964 | fCFContCascadePIDOmegaPlus->SetBinLimits(2, -1.1 , 1.1 ); // Rapidity | |
39b18694 | 965 | Double_t *lBinLim3 = new Double_t[ lNbBinsPerVar[3]+1 ]; |
966 | for(Int_t i=3; i< lNbBinsPerVar[3]+1;i++) lBinLim3[i] = (Double_t)(i-1)*10.; | |
967 | lBinLim3[0] = 0.0; | |
968 | lBinLim3[1] = 5.0; | |
969 | lBinLim3[2] = 10.0; | |
970 | fCFContCascadePIDOmegaPlus->SetBinLimits(3,lBinLim3); // Centrality | |
f9a6cab5 | 971 | |
972 | // Setting the step title : one per PID case | |
973 | fCFContCascadePIDOmegaPlus->SetStepTitle(0, "No PID"); | |
974 | fCFContCascadePIDOmegaPlus->SetStepTitle(1, "TPC PID / 4-#sigma cut on Bachelor track"); | |
975 | fCFContCascadePIDOmegaPlus->SetStepTitle(2, "TPC PID / 4-#sigma cut on Bachelor+Baryon tracks"); | |
976 | fCFContCascadePIDOmegaPlus->SetStepTitle(3, "TPC PID / 4-#sigma cut on Bachelor+Baryon+Meson tracks"); | |
977 | fCFContCascadePIDOmegaPlus->SetStepTitle(4, "Comb. PID / Bachelor"); | |
978 | fCFContCascadePIDOmegaPlus->SetStepTitle(5, "Comb. PID / Bachelor+Baryon"); | |
979 | fCFContCascadePIDOmegaPlus->SetStepTitle(6, "Comb. PID / Bachelor+Baryon+Meson"); | |
980 | ||
981 | // Setting the variable title, per axis | |
982 | fCFContCascadePIDOmegaPlus->SetVarTitle(0, "Pt_{cascade} (GeV/c)"); | |
983 | fCFContCascadePIDOmegaPlus->SetVarTitle(1, "M( #Lambda , K^{+} ) (GeV/c^{2})"); | |
984 | fCFContCascadePIDOmegaPlus->SetVarTitle(2, "Y_{#Omega}"); | |
39b18694 | 985 | fCFContCascadePIDOmegaPlus->SetVarTitle(3, "SPD tracklets Multiplicity"); |
f9a6cab5 | 986 | |
987 | fListHistCascade->Add(fCFContCascadePIDOmegaPlus); | |
988 | ||
989 | } | |
990 | ||
991 | ||
992 | // Part 3 : Towards the optimisation of topological selections ------- | |
39b18694 | 993 | if(! fCFContCascadeCuts) { |
f9a6cab5 | 994 | |
995 | // Container meant to store all the relevant distributions corresponding to the cut variables. | |
11bcd1e4 | 996 | // - NB overflow/underflow of variables on which we want to cut later should be 0!!! |
f9a6cab5 | 997 | const Int_t lNbSteps = 4 ; |
39b18694 | 998 | const Int_t lNbVariables = 22 ; |
f9a6cab5 | 999 | //array for the number of bins in each dimension : |
11bcd1e4 | 1000 | Int_t lNbBinsPerVar[lNbVariables] = {0}; |
f9a6cab5 | 1001 | lNbBinsPerVar[0] = 100; |
11bcd1e4 | 1002 | lNbBinsPerVar[1] = 126; |
39b18694 | 1003 | lNbBinsPerVar[2] = 24; |
1004 | lNbBinsPerVar[3] = 220; | |
f9a6cab5 | 1005 | lNbBinsPerVar[4] = 30; |
1006 | lNbBinsPerVar[5] = 50; | |
1007 | ||
39b18694 | 1008 | lNbBinsPerVar[6] = 101; |
1009 | ||
1010 | lNbBinsPerVar[7] = 102; | |
11bcd1e4 | 1011 | lNbBinsPerVar[8] = 101; |
1012 | lNbBinsPerVar[9] = 26; | |
1013 | lNbBinsPerVar[10] = 26; | |
f9a6cab5 | 1014 | |
11bcd1e4 | 1015 | lNbBinsPerVar[11] = 150; // 75 2-MeV/c2 bins |
1016 | lNbBinsPerVar[12] = 120; // 60 2-MeV/c2 bins | |
f9a6cab5 | 1017 | |
1018 | lNbBinsPerVar[13] = 100; | |
1019 | ||
1020 | lNbBinsPerVar[14] = 44; // 0.05 in rapidity units | |
1021 | lNbBinsPerVar[15] = 44; // 0.05 in rapidity units | |
1022 | ||
1023 | lNbBinsPerVar[16] = 20; | |
1024 | ||
39b18694 | 1025 | lNbBinsPerVar[17] = 11; |
11bcd1e4 | 1026 | lNbBinsPerVar[18] = 100; |
39b18694 | 1027 | lNbBinsPerVar[19] = 112; // Proper time of cascade |
1028 | lNbBinsPerVar[20] = 112; // Proper time of V0 | |
1029 | lNbBinsPerVar[21] = 112; // Distance V0-Xi in transverse plane | |
1030 | ||
11bcd1e4 | 1031 | fCFContCascadeCuts = new AliCFContainer("fCFContCascadeCuts","Container for Cascade cuts", lNbSteps, lNbVariables, lNbBinsPerVar ); |
f9a6cab5 | 1032 | |
1033 | ||
1034 | //setting the bin limits | |
1035 | ||
1036 | //0 | |
1037 | fCFContCascadeCuts->SetBinLimits(0, 0., 2.); // DcaXiDaughters : 0.0 to 2.0 | |
1038 | //1 | |
11bcd1e4 | 1039 | Double_t *lBinLim1 = new Double_t[ lNbBinsPerVar[1]+1 ]; |
39b18694 | 1040 | lBinLim1[0] = 0.0; |
1041 | lBinLim1[1] = 0.03; | |
1042 | for(Int_t i=2; i< lNbBinsPerVar[1];i++) lBinLim1[i] = (Double_t)0.03 + (5. - 0.03 )/(lNbBinsPerVar[1]-2) * (Double_t)(i-1) ; | |
11bcd1e4 | 1043 | lBinLim1[ lNbBinsPerVar[1] ] = 100.0; |
1044 | fCFContCascadeCuts -> SetBinLimits(1, lBinLim1 ); | |
1045 | delete [] lBinLim1; // DcaBachToPrimVertexXi : 0.0 to 0.5 | |
f9a6cab5 | 1046 | //2 |
39b18694 | 1047 | fCFContCascadeCuts->SetBinLimits(2, .9988, 1.); // XiCosineOfPointingAngle : 0.99 to 1.0 |
f9a6cab5 | 1048 | //3 |
39b18694 | 1049 | fCFContCascadeCuts -> SetBinLimits(3, 0., 110.); // XiRadius : 0.0 to 110.0 |
f9a6cab5 | 1050 | //4 |
1051 | fCFContCascadeCuts->SetBinLimits(4, 1.1, 1.13 ); // InvMassLambdaAsCascDghter | |
1052 | //5 | |
1053 | fCFContCascadeCuts->SetBinLimits(5, 0., 2.); // DcaV0DaughtersXi : 0.0 to 2.0 | |
1054 | //6 | |
39b18694 | 1055 | fCFContCascadeCuts->SetBinLimits(6, .95, 1.0005); // V0CosineOfPointingAngleXi : 0.99 to 1.0 |
f9a6cab5 | 1056 | //7 |
11bcd1e4 | 1057 | Double_t *lBinLim7 = new Double_t[ lNbBinsPerVar[7]+1 ]; |
39b18694 | 1058 | for(Int_t i=0; i< lNbBinsPerVar[7]-1;i++) lBinLim7[i] = (Double_t)0.0 + (100. - 0.0 )/(lNbBinsPerVar[7]-2) * (Double_t)i ; |
11bcd1e4 | 1059 | lBinLim7[ lNbBinsPerVar[7]-1] = 200.0; |
1060 | lBinLim7[ lNbBinsPerVar[7]] = 1000.0; | |
1061 | fCFContCascadeCuts -> SetBinLimits(7, lBinLim7 ); | |
39b18694 | 1062 | delete [] lBinLim7; // V0RadiusXi : 0.0 to 100.0 |
f9a6cab5 | 1063 | //8 |
11bcd1e4 | 1064 | Double_t *lBinLim8 = new Double_t[ lNbBinsPerVar[8]+1 ]; |
39b18694 | 1065 | for(Int_t i=0; i< lNbBinsPerVar[8];i++) lBinLim8[i] = (Double_t)0.0 + (0.4 - 0.0 )/(lNbBinsPerVar[8]-1) * (Double_t)i ; |
1066 | lBinLim8[ lNbBinsPerVar[8] ] = 100.0; | |
11bcd1e4 | 1067 | fCFContCascadeCuts -> SetBinLimits(8, lBinLim8 ); |
1068 | delete [] lBinLim8; // DcaV0ToPrimVertexXi : 0.0 to 0.4 | |
f9a6cab5 | 1069 | //9 |
11bcd1e4 | 1070 | Double_t *lBinLim9 = new Double_t[ lNbBinsPerVar[9]+1 ]; |
1071 | for(Int_t i=0; i< lNbBinsPerVar[9];i++) lBinLim9[i] = (Double_t)0.0 + (0.25 - 0.0 )/(lNbBinsPerVar[9]-1) * (Double_t)i ; | |
1072 | lBinLim9[ lNbBinsPerVar[9] ] = 100.0; | |
1073 | fCFContCascadeCuts -> SetBinLimits(9, lBinLim9 ); | |
1074 | delete [] lBinLim9; // DcaPosToPrimVertexXi : 0.0 to 0.25 | |
f9a6cab5 | 1075 | //10 |
11bcd1e4 | 1076 | Double_t *lBinLim10 = new Double_t[ lNbBinsPerVar[10]+1 ]; |
1077 | for(Int_t i=0; i< lNbBinsPerVar[10];i++) lBinLim10[i] = (Double_t)0.0 + (0.25 - 0.0 )/(lNbBinsPerVar[10]-1) * (Double_t)i ; | |
1078 | lBinLim10[ lNbBinsPerVar[10] ] = 100.0; | |
1079 | fCFContCascadeCuts -> SetBinLimits(10, lBinLim10 ); | |
39b18694 | 1080 | delete [] lBinLim10; // DcaNegToPrimVertexXi : 0.0 to 0.25 |
f9a6cab5 | 1081 | //11 |
1082 | fCFContCascadeCuts->SetBinLimits(11, 1.25, 1.40); // InvMassXi | |
1083 | fCFContCascadeCuts->SetBinLimits(12, 1.62, 1.74); // InvMassOmega | |
1084 | fCFContCascadeCuts->SetBinLimits(13, 0.0, 10.0); // XiTransvMom | |
1085 | fCFContCascadeCuts->SetBinLimits(14, -1.1, 1.1); // Y(Xi) | |
1086 | fCFContCascadeCuts->SetBinLimits(15, -1.1, 1.1); // Y(Omega) | |
1087 | fCFContCascadeCuts->SetBinLimits(16, -10.0, 10.0); // BestPrimaryVtxPosZ | |
39b18694 | 1088 | Double_t *lBinLim17 = new Double_t[ lNbBinsPerVar[17]+1 ]; |
1089 | for(Int_t i=3; i< lNbBinsPerVar[17]+1;i++) lBinLim17[i] = (Double_t)(i-1)*10.; | |
1090 | lBinLim17[0] = 0.0; | |
1091 | lBinLim17[1] = 5.0; | |
1092 | lBinLim17[2] = 10.0; | |
1093 | fCFContCascadeCuts -> SetBinLimits(17, lBinLim17 ); // Centrality | |
1094 | delete [] lBinLim17; | |
1095 | fCFContCascadeCuts->SetBinLimits(18, 0.0, 6000.0); // ESD track multiplicity | |
1096 | ||
1097 | Double_t *lBinLim19 = new Double_t[ lNbBinsPerVar[19]+1 ]; | |
1098 | for(Int_t i=0; i< lNbBinsPerVar[19];i++) lBinLim19[i] = (Double_t)-1. + (110. + 1.0 )/(lNbBinsPerVar[19]-1) * (Double_t)i ; | |
1099 | lBinLim19[ lNbBinsPerVar[19] ] = 2000.0; | |
1100 | fCFContCascadeCuts->SetBinLimits(19, lBinLim19); // Proper time cascade | |
1101 | ||
1102 | fCFContCascadeCuts->SetBinLimits(20, lBinLim19); // Proper time V0 | |
1103 | ||
1104 | fCFContCascadeCuts->SetBinLimits(21, lBinLim19); // Distance V0-Xi in tansverse plane | |
1105 | delete [] lBinLim19; | |
f9a6cab5 | 1106 | // Setting the number of steps : one for each cascade species (Xi-, Xi+ and Omega-, Omega+) |
1107 | fCFContCascadeCuts->SetStepTitle(0, "#Xi^{-} candidates"); | |
1108 | fCFContCascadeCuts->SetStepTitle(1, "#bar{#Xi}^{+} candidates"); | |
1109 | fCFContCascadeCuts->SetStepTitle(2, "#Omega^{-} candidates"); | |
1110 | fCFContCascadeCuts->SetStepTitle(3, "#bar{#Omega}^{+} candidates"); | |
1111 | ||
1112 | // Setting the variable title, per axis | |
f9a6cab5 | 1113 | fCFContCascadeCuts->SetVarTitle(0, "Dca(XiDaughters) (cm)"); |
1114 | fCFContCascadeCuts->SetVarTitle(1, "Dca(Bach/PrimVertex) (cm)"); | |
1115 | fCFContCascadeCuts->SetVarTitle(2, "cos(Xi pointing angle)"); | |
1116 | fCFContCascadeCuts->SetVarTitle(3, "R_{2d}(Xi decay) (cm)"); | |
1117 | fCFContCascadeCuts->SetVarTitle(4, "M_{#Lambda}(As Casc Dghter) (GeV/c^{2})"); | |
1118 | fCFContCascadeCuts->SetVarTitle(5, "Dca(V0 Daughters) in Xi (cm)"); | |
1119 | ||
1120 | fCFContCascadeCuts->SetVarTitle(6, "cos(V0 pointing Angle) in Casc"); | |
1121 | fCFContCascadeCuts->SetVarTitle(7, "R_{2d}(V0 decay) (cm)"); | |
1122 | fCFContCascadeCuts->SetVarTitle(8, "Dca(V0/PrimVertex) (cm)"); | |
1123 | fCFContCascadeCuts->SetVarTitle(9, "Dca(Pos/PrimVertex) (cm)"); | |
1124 | fCFContCascadeCuts->SetVarTitle(10, "Dca(Neg/PrimVertex) (cm)"); | |
1125 | ||
1126 | fCFContCascadeCuts->SetVarTitle(11, "Inv. Mass(Xi) (GeV/c^{2})"); | |
1127 | fCFContCascadeCuts->SetVarTitle(12, "Inv. Mass(Omega) (GeV/c^{2})"); | |
1128 | ||
1129 | fCFContCascadeCuts->SetVarTitle(13, "pt(Casc.) (GeV/c)"); | |
1130 | ||
1131 | fCFContCascadeCuts->SetVarTitle(14, "Y(Xi)"); | |
1132 | fCFContCascadeCuts->SetVarTitle(15, "Y(Omega)"); | |
1133 | ||
1134 | fCFContCascadeCuts->SetVarTitle(16, "Z-position(BestPrimVtx) (cm)"); | |
1135 | ||
39b18694 | 1136 | fCFContCascadeCuts->SetVarTitle(17, "Centrality"); |
11bcd1e4 | 1137 | |
1138 | fCFContCascadeCuts->SetVarTitle(18, "ESD track multiplicity"); | |
39b18694 | 1139 | |
1140 | fCFContCascadeCuts->SetVarTitle(19, "Proper time cascade"); | |
1141 | fCFContCascadeCuts->SetVarTitle(20, "Proper time V0"); | |
1142 | fCFContCascadeCuts->SetVarTitle(21, "Distance V0-Xi in transverse plane"); | |
f9a6cab5 | 1143 | |
1144 | fListHistCascade->Add(fCFContCascadeCuts); | |
1145 | } | |
1146 | ||
1147 | fV0Ampl = new TH1F("fV0Ampl","",500,0.,30000); | |
1148 | fListHistCascade->Add(fV0Ampl); | |
1149 | ||
1150 | ||
1151 | fHistDcaXiDaughtersvsInvMass = new TH2F( "fHistDcaXiDaughtersvsInvMass", "DCA between Xi Daughters; DCA (cm) ; Number of Cascades", 100, 0., 0.5,400,1.2,2.0); | |
1152 | fListHistCascade->Add(fHistDcaXiDaughtersvsInvMass); | |
1153 | fHistDcaBachToPrimVertexvsInvMass = new TH2F("fHistDcaBachToPrimVertexvsInvMass", "DCA of Bach. to Prim. Vertex;DCA (cm);Number of Cascades", 250, 0., 0.25,400,1.2,2.0); | |
1154 | fListHistCascade->Add(fHistDcaBachToPrimVertexvsInvMass); | |
1155 | fHistXiCosineOfPointingAnglevsInvMass= new TH2F("fHistXiCosineOfPointingAnglevsInvMass", "Cosine of Xi Pointing Angle; Cos (Xi Point.Angl);Number of Xis", 200, 0.99, 1.0,400,1.2,2.0); | |
1156 | fListHistCascade->Add(fHistXiCosineOfPointingAnglevsInvMass); | |
1157 | fHistMassLambdaAsCascDghtervsInvMass= new TH2F("fHistMassLambdaAsCascDghtervsInvMass","#Lambda associated to Casc. candidates;Eff. Mass (GeV/c^{2});Counts", 300,1.00,1.3,400,1.2,2.0); | |
1158 | fListHistCascade->Add(fHistMassLambdaAsCascDghtervsInvMass); | |
1159 | fHistDcaV0DaughtersXivsInvMass = new TH2F("fHistDcaV0DaughtersXivsInvMass", "DCA between V0 daughters, in cascade;DCA (cm);Number of V0s", 120, 0., 0.6,400,1.2,2.0); | |
1160 | fListHistCascade->Add(fHistDcaV0DaughtersXivsInvMass); | |
1161 | fHistDcaV0ToPrimVertexXivsInvMass = new TH2F("fHistDcaV0ToPrimVertexXivsInvMass", "DCA of V0 to Prim. Vertex, in cascade;DCA (cm);Number of Cascades", 200, 0., 1.,400,1.2,2.0); | |
1162 | fListHistCascade->Add(fHistDcaV0ToPrimVertexXivsInvMass); | |
1163 | ||
f9a6cab5 | 1164 | PostData(1, fListHistCascade); |
1165 | /* PostData(2, fPaveTextBookKeeping);*/ | |
1166 | }// end UserCreateOutputObjects | |
1167 | ||
1168 | ||
f9a6cab5 | 1169 | //________________________________________________________________________ |
1170 | void AliAnalysisTaskCheckCascadePbPb::UserExec(Option_t *) { | |
1171 | // Main loop | |
1172 | // Called for each event | |
1173 | ||
1174 | AliESDEvent *lESDevent = 0x0; | |
1175 | AliAODEvent *lAODevent = 0x0; | |
39b18694 | 1176 | |
1177 | if(!fPIDResponse) { | |
1178 | AliError("Cannot get pid response"); | |
1179 | return; | |
1180 | } | |
1181 | ||
1182 | ||
f9a6cab5 | 1183 | Int_t ncascades = -1; |
1184 | Int_t nTrackMultiplicity = -1; | |
1185 | Int_t nTrackWithTPCrefitMultiplicity = 0; | |
ff9a705e | 1186 | // Int_t lSPDTrackletsMultiplicity = -1; |
1187 | ||
39b18694 | 1188 | // Primary tracks from ESD/AOD |
1189 | Float_t lPrimaryTrackMultiplicity = -1.; | |
1190 | ||
f9a6cab5 | 1191 | |
1192 | Short_t lStatusTrackingPrimVtx = -2; | |
1193 | Double_t lTrkgPrimaryVtxPos[3] = {-100.0, -100.0, -100.0}; | |
1194 | Double_t lBestPrimaryVtxPos[3] = {-100.0, -100.0, -100.0}; | |
1195 | Double_t lMagneticField = -10.; | |
1196 | ||
39b18694 | 1197 | AliCentrality* centrality; |
1198 | AliESDVZERO* esdV0; | |
1199 | AliAODVZERO* aodV0; | |
1200 | Float_t multV0A = 0.; | |
1201 | Float_t multV0C = 0.; | |
1202 | ||
f9a6cab5 | 1203 | // Connect to the InputEvent |
1204 | // After these lines, we should have an ESD/AOD event + the number of cascades in it. | |
ff9a705e | 1205 | |
1206 | // For AOD or ESD ... | |
1207 | nTrackMultiplicity = (InputEvent())->GetNumberOfTracks(); | |
1208 | ||
f9a6cab5 | 1209 | if (fAnalysisType == "ESD") { |
1210 | lESDevent = dynamic_cast<AliESDEvent*>( InputEvent() ); | |
1211 | if (!lESDevent) { | |
1212 | AliWarning("ERROR: lESDevent not available \n"); | |
1213 | return; | |
1214 | } | |
1215 | ||
f9a6cab5 | 1216 | //------------------------------------------------- |
1217 | ||
39b18694 | 1218 | // 0 - Trigger management |
f9a6cab5 | 1219 | // NOTE : Check the availability of the proper trigger --> MN Physics selection moved to runProof macro |
1220 | ||
1221 | //------------------------------------------------- | |
1222 | ||
1223 | // 1 - Cascade vertexer (ESD) | |
39b18694 | 1224 | if (fkRerunV0CascVertexers) { // Relaunch V0 and Cascade vertexers |
1225 | ||
1226 | lESDevent->ResetCascades(); | |
1227 | lESDevent->ResetV0s(); | |
1228 | ||
1229 | AliV0vertexer *lV0vtxer = new AliV0vertexer(); | |
1230 | AliCascadeVertexer *lCascVtxer = new AliCascadeVertexer(); | |
1231 | ||
1232 | //lV0vtxer->GetCuts(fV0Sels); | |
1233 | //lCascVtxer->GetCuts(fCascSels); | |
1234 | ||
1235 | lV0vtxer->SetCuts(fV0Sels); // NB don't use SetDefaultCuts!! because it acts on static variables | |
1236 | lCascVtxer->SetCuts(fCascSels); | |
f9a6cab5 | 1237 | |
39b18694 | 1238 | |
1239 | lV0vtxer->Tracks2V0vertices(lESDevent); | |
1240 | lCascVtxer->V0sTracks2CascadeVertices(lESDevent); | |
1241 | ||
1242 | delete lV0vtxer; | |
1243 | delete lCascVtxer; | |
f9a6cab5 | 1244 | }// end if(RelaunchV0CascVertexers) |
1245 | ||
1246 | //------------------------------------------------- | |
1247 | ncascades = lESDevent->GetNumberOfCascades(); | |
ff9a705e | 1248 | nTrackWithTPCrefitMultiplicity = DoESDTrackWithTPCrefitMultiplicity(lESDevent); |
1249 | ||
1250 | // const AliMultiplicity *lAliMult = lESDevent->GetMultiplicity(); | |
1251 | // lSPDTrackletsMultiplicity = lAliMult->GetNumberOfTracklets(); | |
1252 | ||
39b18694 | 1253 | centrality = lESDevent->GetCentrality(); |
1254 | esdV0 = lESDevent->GetVZEROData(); | |
1255 | multV0A=esdV0->GetMTotV0A(); | |
1256 | multV0C=esdV0->GetMTotV0C(); | |
f9a6cab5 | 1257 | |
39b18694 | 1258 | } else if (fAnalysisType == "AOD") { |
1259 | lAODevent = dynamic_cast<AliAODEvent*>( InputEvent() ); | |
1260 | if (!lAODevent) { | |
1261 | AliWarning("ERROR: lAODevent not available \n"); | |
1262 | return; | |
1263 | } | |
1264 | ncascades = lAODevent->GetNumberOfCascades(); | |
ff9a705e | 1265 | nTrackWithTPCrefitMultiplicity = 0; |
1266 | lPrimaryTrackMultiplicity = 0; | |
1267 | for (Int_t itrack = 0; itrack<nTrackMultiplicity; itrack++) { | |
1268 | AliAODTrack* track = lAODevent->GetTrack(itrack); | |
1269 | if (track->TestFilterBit(AliAODTrack::kTrkGlobalNoDCA)) lPrimaryTrackMultiplicity++; // kTrkGlobal tight DCA cut --> mult is much lower than the one selectied with standard cuts in ESDs | |
1270 | if (track->IsOn(AliAODTrack::kTPCrefit)) nTrackWithTPCrefitMultiplicity++; | |
1271 | } | |
1272 | ||
1273 | // lSPDTrackletsMultiplicity = lAODevent->GetTracklets()->GetNumberOfTracklets(); | |
1274 | ||
39b18694 | 1275 | centrality = lAODevent->GetCentrality(); |
1276 | aodV0 = lAODevent->GetVZEROData(); | |
1277 | multV0A=aodV0->GetMTotV0A(); | |
1278 | multV0C=aodV0->GetMTotV0C(); | |
f9a6cab5 | 1279 | |
39b18694 | 1280 | } else { |
1281 | ||
1282 | Printf("Analysis type (ESD or AOD) not specified \n"); | |
1283 | return; | |
1284 | ||
1285 | } | |
1286 | ||
1287 | fHistCascadeMultiplicityBeforeEvSel->Fill (ncascades); | |
1288 | ||
11bcd1e4 | 1289 | // Printf("Centrality percentile V0M for this event %f)\n", centrality->GetCentralityPercentile("V0M")); |
11bcd1e4 | 1290 | Float_t lcentrality = centrality->GetCentralityPercentile(fCentrEstimator.Data()); |
39b18694 | 1291 | if (lcentrality<fCentrLowLim||lcentrality>=fCentrUpLim) { |
11bcd1e4 | 1292 | PostData(1, fListHistCascade); |
1293 | return; | |
f9a6cab5 | 1294 | } |
1295 | ||
1296 | fV0Ampl->Fill(multV0A+multV0C); | |
1297 | ||
f9a6cab5 | 1298 | //------------------------------------------------- |
1299 | ||
1300 | fHistTrackMultiplicityForCentrEvt ->Fill( nTrackMultiplicity ); | |
1301 | fHistTPCrefitTrackMultiplicityForCentrEvt ->Fill( nTrackWithTPCrefitMultiplicity ); | |
1302 | fHistCascadeMultiplicityForCentrEvt ->Fill( ncascades ); | |
1303 | ||
1304 | // --------------------------------------------------------------- | |
1305 | ||
1306 | // I - Global characteristics of the events + general histos (filled for any selected events and/or for the analysed events) | |
1307 | ||
1308 | // - I.Step 1 : Characteristics of the event : prim. Vtx + magnetic field (ESD) | |
1309 | //------------- | |
1310 | ||
1311 | if (fAnalysisType == "ESD") { | |
1312 | const AliESDVertex *lPrimaryTrackingESDVtx = lESDevent->GetPrimaryVertexTracks(); | |
1313 | // get the vtx stored in ESD found with tracks | |
1314 | lPrimaryTrackingESDVtx->GetXYZ( lTrkgPrimaryVtxPos ); | |
1315 | ||
1316 | const AliESDVertex *lPrimaryBestESDVtx = lESDevent->GetPrimaryVertex(); | |
1317 | // get the best primary vertex available for the event | |
1318 | // As done in AliCascadeVertexer, we keep the one which is the best one available. | |
1319 | // between : Tracking vertex > SPD vertex > TPC vertex > default SPD vertex | |
1320 | // This one will be used for next calculations (DCA essentially) | |
1321 | lPrimaryBestESDVtx->GetXYZ( lBestPrimaryVtxPos ); | |
1322 | lStatusTrackingPrimVtx = lPrimaryTrackingESDVtx->GetStatus(); | |
1323 | ||
f9a6cab5 | 1324 | // FIXME : remove TPC-only primary vertex : retain only events with tracking + SPD vertex |
1325 | if (fkQualityCutNoTPConlyPrimVtx) { | |
1326 | const AliESDVertex *lPrimarySPDVtx = lESDevent->GetPrimaryVertexSPD(); | |
1327 | if (!lPrimarySPDVtx->GetStatus() && !lPrimaryTrackingESDVtx->GetStatus() ){ | |
1328 | AliWarning("Pb / No SPD prim. vertex nor prim. Tracking vertex ... return !"); | |
1329 | PostData(1, fListHistCascade); | |
1330 | return; | |
1331 | } | |
1332 | } | |
1333 | ||
1334 | // NOTE : For older evts | |
1335 | ||
1336 | // As previously done in AliCascadeVertexer, we keep, between both retrieved vertices (SPD or Tracking) | |
1337 | // the one which is the best one available. | |
1338 | // This one will be used for next calculations (DCA essentially) | |
1339 | // At that time, the TPC-only primary vertex was not considered | |
1340 | ||
1341 | ||
1342 | lMagneticField = lESDevent->GetMagneticField( ); | |
1343 | // FIXME if(TMath::Abs(lMagneticField ) < 10e-6) continue; | |
39b18694 | 1344 | lPrimaryTrackMultiplicity = fESDtrackCuts->CountAcceptedTracks(lESDevent); |
ff9a705e | 1345 | } else if (fAnalysisType == "AOD") { |
39b18694 | 1346 | |
1347 | const AliAODVertex *lPrimaryBestAODVtx = lAODevent->GetPrimaryVertex(); | |
1348 | if (!lPrimaryBestAODVtx){ | |
1349 | AliWarning("No prim. vertex in AOD... return!"); | |
1350 | PostData(1, fListHistCascade); | |
1351 | return; | |
1352 | } | |
1353 | ||
1354 | // get the best primary vertex available for the event GetVertex(0) | |
f9a6cab5 | 1355 | // This one will be used for next calculations (DCA essentially) |
1356 | lPrimaryBestAODVtx->GetXYZ( lBestPrimaryVtxPos ); | |
1357 | ||
1358 | lStatusTrackingPrimVtx = -1; | |
1359 | lTrkgPrimaryVtxPos[0] = -100.0; | |
1360 | lTrkgPrimaryVtxPos[1] = -100.0; | |
39b18694 | 1361 | lTrkgPrimaryVtxPos[2] = -100.0; |
1362 | lMagneticField = lAODevent->GetMagneticField(); | |
1363 | ||
39b18694 | 1364 | } |
1365 | ||
1366 | // Quality cut on the z-position of the prim vertex. | |
1367 | if (fkQualityCutZprimVtxPos) { | |
1368 | if (TMath::Abs(lBestPrimaryVtxPos[2]) > fVtxRange ) { | |
1369 | AliWarning("Pb / | Z position of Best Prim Vtx | > 10.0 cm ... return !"); | |
1370 | PostData(1, fListHistCascade); | |
1371 | return; | |
1372 | } | |
f9a6cab5 | 1373 | } |
1374 | ||
1375 | // - I.Step 2 : Filling histos that characterize the selected event : x,y,z prim. Vtx distrib. (ESD) | |
1376 | //------------- | |
1377 | ||
1378 | fHistCascadeMultiplicityForSelEvt ->Fill( ncascades ); | |
1379 | fHistPosBestPrimaryVtxXForSelEvt ->Fill( lBestPrimaryVtxPos[0] ); | |
1380 | fHistPosBestPrimaryVtxYForSelEvt ->Fill( lBestPrimaryVtxPos[1] ); | |
1381 | fHistPosBestPrimaryVtxZForSelEvt ->Fill( lBestPrimaryVtxPos[2] ); | |
1382 | ||
1383 | // --------------------------------------------------------------- | |
11bcd1e4 | 1384 | |
39b18694 | 1385 | fHistEvtsInCentralityBinsvsNtracks->Fill(lcentrality,lPrimaryTrackMultiplicity); |
11bcd1e4 | 1386 | |
11bcd1e4 | 1387 | |
11bcd1e4 | 1388 | |
1389 | // II - Calculation Part dedicated to Xi vertices | |
f9a6cab5 | 1390 | |
1391 | for (Int_t iXi = 0; iXi < ncascades; iXi++) {// This is the begining of the Cascade loop (ESD or AOD) | |
1392 | ||
1393 | // ------------------------------------- | |
1394 | // II.Init - Initialisation of the local variables that will be needed for ESD/AOD | |
1395 | ||
1396 | ||
1397 | // - 0th part of initialisation : around primary vertex ... | |
1398 | ||
1399 | Double_t lTrkgPrimaryVtxRadius3D = -500.0; | |
1400 | Double_t lBestPrimaryVtxRadius3D = -500.0; | |
1401 | ||
1402 | // - 1st part of initialisation : variables needed to store AliESDCascade data members | |
1403 | Double_t lEffMassXi = 0. ; | |
1404 | Double_t lChi2Xi = -1. ; | |
1405 | Double_t lDcaXiDaughters = -1. ; | |
1406 | Double_t lXiCosineOfPointingAngle = -1. ; | |
1407 | Double_t lPosXi[3] = { -1000.0, -1000.0, -1000.0 }; | |
1408 | Double_t lXiRadius = -1000. ; | |
1409 | ||
1410 | // - 2nd part of initialisation : Nbr of clusters within TPC for the 3 daughter cascade tracks | |
39b18694 | 1411 | UShort_t lPosTPCClusters = -1; // For ESD only ... |
1412 | UShort_t lNegTPCClusters = -1; // For ESD only ... | |
1413 | UShort_t lBachTPCClusters = -1; // For ESD only ... | |
f9a6cab5 | 1414 | |
1415 | Double_t lInnerWallMomCascDghters[3] = {-100., -100., -100.}; | |
1416 | Double_t lTPCSignalCascDghters [3] = {-100., -100., -100.}; | |
1417 | ||
1418 | ||
1419 | // - 3rd part of initialisation : about V0 part in cascades | |
1420 | Double_t lInvMassLambdaAsCascDghter = 0.; | |
1421 | Double_t lV0Chi2Xi = -1. ; | |
1422 | Double_t lDcaV0DaughtersXi = -1.; | |
1423 | ||
1424 | Double_t lDcaBachToPrimVertexXi = -1., lDcaV0ToPrimVertexXi = -1.; | |
1425 | Double_t lDcaPosToPrimVertexXi = -1.; | |
1426 | Double_t lDcaNegToPrimVertexXi = -1.; | |
1427 | Double_t lV0CosineOfPointingAngleXi = -1. ; | |
1428 | Double_t lPosV0Xi[3] = { -1000. , -1000., -1000. }; // Position of VO coming from cascade | |
1429 | Double_t lV0RadiusXi = -1000.0; | |
1430 | Double_t lV0quality = 0.; | |
1431 | ||
1432 | ||
1433 | // - 4th part of initialisation : Effective masses | |
1434 | Double_t lInvMassXiMinus = 0.; | |
1435 | Double_t lInvMassXiPlus = 0.; | |
1436 | Double_t lInvMassOmegaMinus = 0.; | |
1437 | Double_t lInvMassOmegaPlus = 0.; | |
1438 | ||
1439 | // - 5th part of initialisation : PID treatment | |
1440 | Bool_t lIsPosInXiProton = kFALSE; | |
1441 | Bool_t lIsPosInXiPion = kFALSE; | |
1442 | Bool_t lIsPosInOmegaProton = kFALSE; | |
1443 | Bool_t lIsPosInOmegaPion = kFALSE; | |
1444 | ||
1445 | Bool_t lIsNegInXiProton = kFALSE; | |
1446 | Bool_t lIsNegInXiPion = kFALSE; | |
1447 | Bool_t lIsNegInOmegaProton = kFALSE; | |
1448 | Bool_t lIsNegInOmegaPion = kFALSE; | |
1449 | ||
1450 | Bool_t lIsBachelorKaon = kFALSE; | |
1451 | Bool_t lIsBachelorPion = kFALSE; | |
1452 | ||
39b18694 | 1453 | Bool_t lIsBachelorKaonForTPC = kFALSE; |
1454 | Bool_t lIsBachelorPionForTPC = kFALSE; | |
1455 | Bool_t lIsNegPionForTPC = kFALSE; | |
1456 | Bool_t lIsPosPionForTPC = kFALSE; | |
1457 | Bool_t lIsNegProtonForTPC = kFALSE; | |
1458 | Bool_t lIsPosProtonForTPC = kFALSE; | |
f9a6cab5 | 1459 | |
1460 | // - 6th part of initialisation : extra info for QA | |
1461 | Double_t lXiMomX = 0., lXiMomY = 0., lXiMomZ = 0.; | |
1462 | Double_t lXiTransvMom = 0. ; | |
1463 | Double_t lXiTotMom = 0. ; | |
39b18694 | 1464 | |
1465 | Double_t lV0PMomX = 0., lV0PMomY = 0., lV0PMomZ = 0.; | |
1466 | Double_t lV0NMomX = 0., lV0NMomY = 0., lV0NMomZ = 0.; | |
1467 | Double_t lV0TotMom = 0. ; | |
1468 | ||
f9a6cab5 | 1469 | |
1470 | Double_t lBachMomX = 0., lBachMomY = 0., lBachMomZ = 0.; | |
1471 | Double_t lBachTransvMom = 0.; | |
1472 | Double_t lBachTotMom = 0.; | |
1473 | ||
1474 | Short_t lChargeXi = -2; | |
1475 | Double_t lV0toXiCosineOfPointingAngle = 0. ; | |
1476 | ||
1477 | Double_t lRapXi = -20.0, lRapOmega = -20.0, lEta = -20.0, lTheta = 360., lPhi = 720. ; | |
1478 | Double_t lAlphaXi = -200., lPtArmXi = -200.0; | |
1479 | ||
1480 | // - 7th part of initialisation : variables for the AliCFContainer dedicated to cascade cut optmisiation | |
f9a6cab5 | 1481 | |
f9a6cab5 | 1482 | if (fAnalysisType == "ESD") { |
1483 | ||
1484 | // ------------------------------------- | |
1485 | // II.ESD - Calculation Part dedicated to Xi vertices (ESD) | |
1486 | ||
1487 | AliESDcascade *xi = lESDevent->GetCascade(iXi); | |
1488 | if (!xi) continue; | |
1489 | ||
1490 | ||
1491 | // - II.Step 1 : around primary vertex | |
1492 | //------------- | |
1493 | lTrkgPrimaryVtxRadius3D = TMath::Sqrt( lTrkgPrimaryVtxPos[0] * lTrkgPrimaryVtxPos[0] + | |
1494 | lTrkgPrimaryVtxPos[1] * lTrkgPrimaryVtxPos[1] + | |
1495 | lTrkgPrimaryVtxPos[2] * lTrkgPrimaryVtxPos[2] ); | |
1496 | ||
1497 | lBestPrimaryVtxRadius3D = TMath::Sqrt( lBestPrimaryVtxPos[0] * lBestPrimaryVtxPos[0] + | |
1498 | lBestPrimaryVtxPos[1] * lBestPrimaryVtxPos[1] + | |
1499 | lBestPrimaryVtxPos[2] * lBestPrimaryVtxPos[2] ); | |
1500 | ||
1501 | ||
1502 | ||
1503 | // - II.Step 2 : Assigning the necessary variables for specific AliESDcascade data members (ESD) | |
1504 | //------------- | |
1505 | lV0quality = 0.; | |
1506 | xi->ChangeMassHypothesis(lV0quality , 3312); // default working hypothesis : cascade = Xi- decay | |
1507 | ||
1508 | lEffMassXi = xi->GetEffMassXi(); | |
1509 | lChi2Xi = xi->GetChi2Xi(); | |
1510 | lDcaXiDaughters = xi->GetDcaXiDaughters(); | |
1511 | lXiCosineOfPointingAngle = xi->GetCascadeCosineOfPointingAngle( lBestPrimaryVtxPos[0], | |
1512 | lBestPrimaryVtxPos[1], | |
1513 | lBestPrimaryVtxPos[2] ); | |
1514 | // Take care : the best available vertex should be used (like in AliCascadeVertexer) | |
1515 | ||
1516 | xi->GetXYZcascade( lPosXi[0], lPosXi[1], lPosXi[2] ); | |
1517 | lXiRadius = TMath::Sqrt( lPosXi[0]*lPosXi[0] + lPosXi[1]*lPosXi[1] ); | |
1518 | ||
1519 | ||
1520 | ||
1521 | // - II.Step 3 : around the tracks : Bach + V0 (ESD) | |
1522 | // ~ Necessary variables for ESDcascade data members coming from the ESDv0 part (inheritance) | |
1523 | //------------- | |
1524 | ||
1525 | UInt_t lIdxPosXi = (UInt_t) TMath::Abs( xi->GetPindex() ); | |
1526 | UInt_t lIdxNegXi = (UInt_t) TMath::Abs( xi->GetNindex() ); | |
1527 | UInt_t lBachIdx = (UInt_t) TMath::Abs( xi->GetBindex() ); | |
1528 | // Care track label can be negative in MC production (linked with the track quality) | |
1529 | // However = normally, not the case for track index ... | |
1530 | ||
1531 | // FIXME : rejection of a double use of a daughter track (nothing but just a crosscheck of what is done in the cascade vertexer) | |
1532 | if (lBachIdx == lIdxNegXi) { | |
1533 | AliWarning("Pb / Idx(Bach. track) = Idx(Neg. track) ... continue!"); continue; | |
1534 | } | |
1535 | if (lBachIdx == lIdxPosXi) { | |
1536 | AliWarning("Pb / Idx(Bach. track) = Idx(Pos. track) ... continue!"); continue; | |
1537 | } | |
1538 | ||
1539 | AliESDtrack *pTrackXi = lESDevent->GetTrack( lIdxPosXi ); | |
1540 | AliESDtrack *nTrackXi = lESDevent->GetTrack( lIdxNegXi ); | |
1541 | AliESDtrack *bachTrackXi = lESDevent->GetTrack( lBachIdx ); | |
1542 | if (!pTrackXi || !nTrackXi || !bachTrackXi ) { | |
1543 | AliWarning("ERROR: Could not retrieve one of the 3 ESD daughter tracks of the cascade ..."); | |
1544 | continue; | |
1545 | } | |
1546 | ||
f9a6cab5 | 1547 | lPosTPCClusters = pTrackXi->GetTPCNcls(); |
1548 | lNegTPCClusters = nTrackXi->GetTPCNcls(); | |
1549 | lBachTPCClusters = bachTrackXi->GetTPCNcls(); | |
1550 | ||
1551 | // FIXME : rejection of a poor quality tracks | |
1552 | if (fkQualityCutTPCrefit) { | |
1553 | // 1 - Poor quality related to TPCrefit | |
1554 | ULong_t pStatus = pTrackXi->GetStatus(); | |
1555 | ULong_t nStatus = nTrackXi->GetStatus(); | |
1556 | ULong_t bachStatus = bachTrackXi->GetStatus(); | |
1557 | if ((pStatus&AliESDtrack::kTPCrefit) == 0) { AliWarning("Pb / V0 Pos. track has no TPCrefit ... continue!"); continue; } | |
1558 | if ((nStatus&AliESDtrack::kTPCrefit) == 0) { AliWarning("Pb / V0 Neg. track has no TPCrefit ... continue!"); continue; } | |
1559 | if ((bachStatus&AliESDtrack::kTPCrefit) == 0) { AliWarning("Pb / Bach. track has no TPCrefit ... continue!"); continue; } | |
1560 | } | |
1561 | if (fkQualityCut80TPCcls) { | |
1562 | // 2 - Poor quality related to TPC clusters | |
1563 | if (lPosTPCClusters < 80) { AliWarning("Pb / V0 Pos. track has less than 80 TPC clusters ... continue!"); continue; } | |
1564 | if (lNegTPCClusters < 80) { AliWarning("Pb / V0 Neg. track has less than 80 TPC clusters ... continue!"); continue; } | |
39b18694 | 1565 | if (lBachTPCClusters < 80) { AliWarning("Pb / Bach. track has less than 80 TPC clusters ... continue!"); continue; } |
f9a6cab5 | 1566 | } |
1567 | ||
1568 | const AliExternalTrackParam *pExtTrack = pTrackXi ->GetInnerParam(); | |
1569 | const AliExternalTrackParam *nExtTrack = nTrackXi ->GetInnerParam(); | |
1570 | const AliExternalTrackParam *bachExtTrack = bachTrackXi ->GetInnerParam(); | |
1571 | ||
1572 | if (pExtTrack) { | |
1573 | lInnerWallMomCascDghters[0] = pExtTrack ->GetP() * pExtTrack ->Charge(); | |
1574 | lTPCSignalCascDghters [0] = pTrackXi ->GetTPCsignal(); | |
1575 | } | |
1576 | if (nExtTrack) { | |
1577 | lInnerWallMomCascDghters[1] = nExtTrack ->GetP() * nExtTrack ->Charge(); | |
1578 | lTPCSignalCascDghters [1] = nTrackXi ->GetTPCsignal(); | |
1579 | } | |
1580 | if (bachExtTrack) { | |
1581 | lInnerWallMomCascDghters[2] = bachExtTrack ->GetP() * bachExtTrack ->Charge(); | |
1582 | lTPCSignalCascDghters [2] = bachTrackXi ->GetTPCsignal(); | |
1583 | } | |
1584 | ||
1585 | ||
1586 | lInvMassLambdaAsCascDghter = xi->GetEffMass(); | |
1587 | // This value shouldn't change, whatever the working hyp. is : Xi-, Xi+, Omega-, Omega+ | |
1588 | lDcaV0DaughtersXi = xi->GetDcaV0Daughters(); | |
1589 | lV0Chi2Xi = xi->GetChi2V0(); | |
1590 | ||
39b18694 | 1591 | lV0CosineOfPointingAngleXi = xi->GetV0CosineOfPointingAngle(lBestPrimaryVtxPos[0], |
1592 | lBestPrimaryVtxPos[1], | |
1593 | lBestPrimaryVtxPos[2] ); | |
1594 | //if (lV0CosineOfPointingAngleXi==1.) cout << "Cosine V0 PA wrt Xi pos ==1!" <<endl; | |
1595 | ||
f9a6cab5 | 1596 | lDcaV0ToPrimVertexXi = xi->GetD( lBestPrimaryVtxPos[0], |
1597 | lBestPrimaryVtxPos[1], | |
1598 | lBestPrimaryVtxPos[2] ); | |
1599 | ||
1600 | lDcaBachToPrimVertexXi = TMath::Abs( bachTrackXi->GetD( lBestPrimaryVtxPos[0], | |
1601 | lBestPrimaryVtxPos[1], | |
1602 | lMagneticField ) ); | |
1603 | // Note : AliExternalTrackParam::GetD returns an algebraic value ... | |
1604 | ||
1605 | xi->GetXYZ( lPosV0Xi[0], lPosV0Xi[1], lPosV0Xi[2] ); | |
1606 | lV0RadiusXi = TMath::Sqrt( lPosV0Xi[0]*lPosV0Xi[0] + lPosV0Xi[1]*lPosV0Xi[1] ); | |
1607 | ||
1608 | lDcaPosToPrimVertexXi = TMath::Abs( pTrackXi ->GetD( lBestPrimaryVtxPos[0], | |
1609 | lBestPrimaryVtxPos[1], | |
1610 | lMagneticField ) ); | |
1611 | ||
1612 | lDcaNegToPrimVertexXi = TMath::Abs( nTrackXi ->GetD( lBestPrimaryVtxPos[0], | |
1613 | lBestPrimaryVtxPos[1], | |
1614 | lMagneticField ) ); | |
1615 | ||
1616 | // - II.Step 3' : extra-selection for cascade candidates | |
f9a6cab5 | 1617 | if (fkExtraSelections) { |
39b18694 | 1618 | // if(lChi2Xi > 2000) continue; // in AliCascadeVertexer |
1619 | // if(lV0Chi2Xi > 2000) continue; // in AliV0vertexer | |
f9a6cab5 | 1620 | |
11bcd1e4 | 1621 | if (lDcaXiDaughters > 0.3) continue; // in AliCascadeVertexer |
f9a6cab5 | 1622 | if (lXiCosineOfPointingAngle < 0.999 ) continue; // in AliCascadeVertexer |
1623 | if (lDcaV0ToPrimVertexXi < 0.05) continue; // in AliCascadeVertexer | |
1624 | if (lDcaBachToPrimVertexXi < 0.03) continue; // in AliCascadeVertexer | |
39b18694 | 1625 | // if (TMath::Abs(lInvMassLambdaAsCascDghter-1.11568) > 0.006 ) continue; // in AliCascadeVertexer |
f9a6cab5 | 1626 | |
1627 | if (lDcaV0DaughtersXi > 1.) continue; // in AliV0vertexer | |
1628 | if (lV0CosineOfPointingAngleXi < 0.998) continue; // in AliV0vertexer | |
1629 | if (lDcaPosToPrimVertexXi < 0.1) continue; // in AliV0vertexer | |
1630 | if (lDcaNegToPrimVertexXi < 0.1) continue; // in AliV0vertexer | |
1631 | ||
1632 | ||
39b18694 | 1633 | if(lXiRadius < .9) continue; // in AliCascadeVertexer |
1634 | // if(lXiRadius > 100) continue; // in AliCascadeVertexer | |
1635 | if(lV0RadiusXi < 0.9) continue; // in AliV0vertexer | |
1636 | // if(lV0RadiusXi > 100) continue; // in AliV0vertexer | |
f9a6cab5 | 1637 | |
1638 | } | |
1639 | ||
1640 | ||
1641 | ||
1642 | // - II.Step 4 : around effective masses (ESD) | |
1643 | // ~ change mass hypotheses to cover all the possibilities : Xi-/+, Omega -/+ | |
1644 | //------------- | |
1645 | ||
1646 | ||
1647 | if ( bachTrackXi->Charge() < 0 ) { | |
1648 | lV0quality = 0.; | |
1649 | xi->ChangeMassHypothesis(lV0quality , 3312); | |
1650 | // Calculate the effective mass of the Xi- candidate. | |
1651 | // pdg code 3312 = Xi- | |
1652 | lInvMassXiMinus = xi->GetEffMassXi(); | |
1653 | ||
1654 | lV0quality = 0.; | |
1655 | xi->ChangeMassHypothesis(lV0quality , 3334); | |
1656 | // Calculate the effective mass of the Xi- candidate. | |
1657 | // pdg code 3334 = Omega- | |
1658 | lInvMassOmegaMinus = xi->GetEffMassXi(); | |
1659 | ||
1660 | lV0quality = 0.; | |
1661 | xi->ChangeMassHypothesis(lV0quality , 3312); // Back to default hyp. | |
1662 | }// end if negative bachelor | |
1663 | ||
1664 | ||
1665 | if ( bachTrackXi->Charge() > 0 ) { | |
1666 | lV0quality = 0.; | |
1667 | xi->ChangeMassHypothesis(lV0quality , -3312); | |
1668 | // Calculate the effective mass of the Xi+ candidate. | |
1669 | // pdg code -3312 = Xi+ | |
1670 | lInvMassXiPlus = xi->GetEffMassXi(); | |
1671 | ||
1672 | lV0quality = 0.; | |
1673 | xi->ChangeMassHypothesis(lV0quality , -3334); | |
1674 | // Calculate the effective mass of the Xi+ candidate. | |
1675 | // pdg code -3334 = Omega+ | |
1676 | lInvMassOmegaPlus = xi->GetEffMassXi(); | |
1677 | ||
1678 | lV0quality = 0.; | |
1679 | xi->ChangeMassHypothesis(lV0quality , -3312); // Back to "default" hyp. | |
1680 | }// end if positive bachelor | |
39b18694 | 1681 | |
f9a6cab5 | 1682 | // - II.Step 5 : PID on the daughter tracks |
1683 | //------------- | |
39b18694 | 1684 | /* |
1685 | // A - Combined PID */ // removed, add when will be used | |
f9a6cab5 | 1686 | |
1687 | // B - TPC PID : 3-sigma bands on Bethe-Bloch curve | |
1688 | ||
1689 | // Bachelor | |
11bcd1e4 | 1690 | if (TMath::Abs(fPIDResponse->NumberOfSigmasTPC( bachTrackXi,AliPID::kKaon)) < 4) lIsBachelorKaonForTPC = kTRUE; |
1691 | if (TMath::Abs(fPIDResponse->NumberOfSigmasTPC( bachTrackXi,AliPID::kPion)) < 4) lIsBachelorPionForTPC = kTRUE; | |
f9a6cab5 | 1692 | |
1693 | // Negative V0 daughter | |
11bcd1e4 | 1694 | if (TMath::Abs(fPIDResponse->NumberOfSigmasTPC( nTrackXi,AliPID::kPion )) < 4) lIsNegPionForTPC = kTRUE; |
1695 | if (TMath::Abs(fPIDResponse->NumberOfSigmasTPC( nTrackXi,AliPID::kProton )) < 4) lIsNegProtonForTPC = kTRUE; | |
f9a6cab5 | 1696 | |
1697 | // Positive V0 daughter | |
11bcd1e4 | 1698 | if (TMath::Abs(fPIDResponse->NumberOfSigmasTPC( pTrackXi,AliPID::kPion )) < 4) lIsPosPionForTPC = kTRUE; |
1699 | if (TMath::Abs(fPIDResponse->NumberOfSigmasTPC( pTrackXi,AliPID::kProton )) < 4) lIsPosProtonForTPC = kTRUE; | |
f9a6cab5 | 1700 | |
1701 | /* | |
1702 | const AliExternalTrackParam *pInnerWallTrackXi = pTrackXi ->GetInnerParam(); | |
1703 | const AliExternalTrackParam *nInnerWallTrackXi = nTrackXi ->GetInnerParam(); | |
1704 | const AliExternalTrackParam *bachInnerWallTrackXi = bachTrackXi ->GetInnerParam(); | |
1705 | if (pInnerWallTrackXi && nInnerWallTrackXi && bachInnerWallTrackXi ) { | |
1706 | ||
1707 | Double_t pMomInnerWall = pInnerWallTrackXi ->GetP(); | |
1708 | Double_t nMomInnerWall = nInnerWallTrackXi ->GetP(); | |
1709 | Double_t bachMomInnerWall = bachInnerWallTrackXi->GetP(); | |
1710 | ||
1711 | // Bachelor | |
1712 | if (TMath::Abs(fESDpid->NumberOfSigmasTPC( bachTrackXi,AliPID::kPion)) < 3) lIsBachelorPionForTPC = kTRUE; | |
1713 | if (bachMomInnerWall < 0.350 && TMath::Abs(fESDpid->NumberOfSigmasTPC( bachTrackXi,AliPID::kKaon)) < 5) lIsBachelorKaonForTPC = kTRUE; | |
1714 | if (bachMomInnerWall > 0.350 && TMath::Abs(fESDpid->NumberOfSigmasTPC( bachTrackXi,AliPID::kKaon)) < 3) lIsBachelorKaonForTPC = kTRUE; | |
1715 | ||
1716 | // Negative V0 daughter | |
1717 | if (TMath::Abs(fESDpid->NumberOfSigmasTPC( nTrackXi,AliPID::kPion )) < 3 ) lIsNegPionForTPC = kTRUE; | |
1718 | if (nMomInnerWall < 0.6 && TMath::Abs(fESDpid->NumberOfSigmasTPC( nTrackXi,AliPID::kProton ) ) < 5 ) lIsNegProtonForTPC = kTRUE; | |
1719 | if (nMomInnerWall > 0.6 && TMath::Abs(fESDpid->NumberOfSigmasTPC( nTrackXi,AliPID::kProton ) ) < 3 ) lIsNegProtonForTPC = kTRUE; | |
1720 | ||
1721 | // Positive V0 daughter | |
1722 | if (TMath::Abs(fESDpid->NumberOfSigmasTPC( pTrackXi,AliPID::kPion )) < 3 ) lIsPosPionForTPC = kTRUE; | |
1723 | if (pMomInnerWall < 0.6 && TMath::Abs(fESDpid->NumberOfSigmasTPC( pTrackXi,AliPID::kProton )) < 5) lIsPosProtonForTPC = kTRUE; | |
1724 | if (pMomInnerWall > 0.6 && TMath::Abs(fESDpid->NumberOfSigmasTPC( pTrackXi,AliPID::kProton )) < 3) lIsPosProtonForTPC = kTRUE; | |
1725 | } | |
1726 | */ | |
1727 | ||
1728 | // - II.Step 6 : extra info for QA (ESD) | |
1729 | // miscellaneous pieces of info that may help regarding data quality assessment. | |
1730 | //------------- | |
1731 | ||
1732 | xi->GetPxPyPz( lXiMomX, lXiMomY, lXiMomZ ); | |
1733 | lXiTransvMom = TMath::Sqrt( lXiMomX*lXiMomX + lXiMomY*lXiMomY ); | |
1734 | lXiTotMom = TMath::Sqrt( lXiMomX*lXiMomX + lXiMomY*lXiMomY + lXiMomZ*lXiMomZ ); | |
39b18694 | 1735 | |
1736 | xi->GetNPxPyPz(lV0NMomX,lV0NMomY,lV0NMomZ); | |
1737 | xi->GetPPxPyPz(lV0PMomX,lV0PMomY,lV0PMomZ); | |
1738 | lV0TotMom = TMath::Sqrt(TMath::Power(lV0NMomX+lV0PMomX,2)+TMath::Power(lV0NMomY+lV0PMomY,2)+TMath::Power(lV0NMomZ+lV0PMomZ,2)); | |
f9a6cab5 | 1739 | |
1740 | xi->GetBPxPyPz( lBachMomX, lBachMomY, lBachMomZ ); | |
1741 | lBachTransvMom = TMath::Sqrt( lBachMomX*lBachMomX + lBachMomY*lBachMomY ); | |
1742 | lBachTotMom = TMath::Sqrt( lBachMomX*lBachMomX + lBachMomY*lBachMomY + lBachMomZ*lBachMomZ ); | |
1743 | ||
1744 | lChargeXi = xi->Charge(); | |
1745 | ||
1746 | lV0toXiCosineOfPointingAngle = xi->GetV0CosineOfPointingAngle( lPosXi[0], lPosXi[1], lPosXi[2] ); | |
1747 | ||
1748 | lRapXi = xi->RapXi(); | |
1749 | lRapOmega = xi->RapOmega(); | |
1750 | lEta = xi->Eta(); | |
1751 | lTheta = xi->Theta() *180.0/TMath::Pi(); | |
1752 | lPhi = xi->Phi() *180.0/TMath::Pi(); | |
1753 | lAlphaXi = xi->AlphaXi(); | |
1754 | lPtArmXi = xi->PtArmXi(); | |
1755 | ||
1756 | ||
1757 | //FIXME : Extra-cut = Anti-splitting cut for lambda daughters | |
1758 | Bool_t kAntiSplittingLambda = kFALSE; | |
1759 | ||
1760 | if (kAntiSplittingLambda) { // not used | |
1761 | Double_t lNMomX = 0., lNMomY = 0., lNMomZ = 0.; | |
1762 | Double_t lPMomX = 0., lPMomY = 0., lPMomZ = 0.; | |
1763 | ||
1764 | xi->GetPPxPyPz(lPMomX, lPMomY, lPMomZ); | |
1765 | xi->GetNPxPyPz(lNMomX, lNMomY, lNMomZ); | |
1766 | ||
1767 | if ( xi->Charge() < 0) {// Xi- or Omega- | |
1768 | if (TMath::Abs(lBachTransvMom - TMath::Sqrt( lNMomX*lNMomX + lNMomY*lNMomY ) ) < 0.075) continue; | |
1769 | } else { //Xi+ or Omega+ | |
1770 | if(TMath::Abs(lBachTransvMom - TMath::Sqrt( lPMomX*lPMomX + lPMomY*lPMomY ) ) < 0.075) continue; | |
1771 | } | |
1772 | } | |
1773 | ||
1774 | //FIXME : Just to know which file is currently open : locate the file containing Xi | |
1775 | // cout << "Name of the file containing Xi candidate(s) :" | |
1776 | // << CurrentFileName() | |
1777 | // << " / entry: " << Entry() | |
1778 | // << " / in file: " << lESDevent->GetEventNumberInFile() // <- Cvetan / From Mihaela: AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()->GetTree()->GetReadEntry(); | |
1779 | // << " : mass(Xi) = " << xi->GetEffMassXi() | |
1780 | // << " / charge = " << lChargeXi | |
1781 | // << " / pt(Casc) = " << lXiTransvMom | |
1782 | // << " / Decay 2d R(Xi) = " << lXiRadius | |
1783 | // << " / Track Index(Pos) = " << lIdxPosXi << "/ Nb(TPC clusters) = " << lPosTPCClusters | |
1784 | // << " / Track Index(Neg) = " << lIdxNegXi << "/ Nb(TPC clusters) = " << lNegTPCClusters | |
1785 | // << " / Track Index(Bach) = " << lBachIdx << "/ Nb(TPC clusters) = " << lBachTPCClusters | |
1786 | // << endl; | |
1787 | ||
1788 | ||
ff9a705e | 1789 | } else if (fAnalysisType == "AOD") { |
39b18694 | 1790 | |
1791 | // II.AOD - Calculation Part dedicated to Xi vertices | |
1792 | ||
f9a6cab5 | 1793 | const AliAODcascade *xi = lAODevent->GetCascade(iXi); |
1794 | if (!xi) continue; | |
1795 | ||
f9a6cab5 | 1796 | // - II.Step 1 : Characteristics of the event : prim. Vtx + magnetic field (AOD) |
1797 | //------------- | |
1798 | ||
1799 | ||
1800 | lTrkgPrimaryVtxRadius3D = -500. ; | |
39b18694 | 1801 | // We don't have the different prim. vertex at the AOD level -> nothing to do. |
f9a6cab5 | 1802 | |
1803 | lBestPrimaryVtxRadius3D = TMath::Sqrt( lBestPrimaryVtxPos[0] * lBestPrimaryVtxPos[0] + | |
1804 | lBestPrimaryVtxPos[1] * lBestPrimaryVtxPos[1] + | |
1805 | lBestPrimaryVtxPos[2] * lBestPrimaryVtxPos[2] ); | |
1806 | ||
1807 | ||
1808 | // - II.Step 2 : Assigning the necessary variables for specific AliAODcascade data members (AOD) | |
1809 | //------------- | |
1810 | ||
39b18694 | 1811 | lEffMassXi = xi->MassXi(); // default working hypothesis : cascade = Xi- decay |
f9a6cab5 | 1812 | lChi2Xi = xi->Chi2Xi(); |
39b18694 | 1813 | lDcaXiDaughters = xi->DcaXiDaughters(); |
f9a6cab5 | 1814 | lXiCosineOfPointingAngle = xi->CosPointingAngleXi( lBestPrimaryVtxPos[0], |
39b18694 | 1815 | lBestPrimaryVtxPos[1], |
1816 | lBestPrimaryVtxPos[2] ); | |
f9a6cab5 | 1817 | // Take care : |
1818 | // the best available vertex should be used (like in AliCascadeVertexer) | |
1819 | ||
1820 | lPosXi[0] = xi->DecayVertexXiX(); | |
1821 | lPosXi[1] = xi->DecayVertexXiY(); | |
1822 | lPosXi[2] = xi->DecayVertexXiZ(); | |
1823 | lXiRadius = TMath::Sqrt( lPosXi[0]*lPosXi[0] + lPosXi[1]*lPosXi[1] ); | |
1824 | ||
39b18694 | 1825 | AliAODTrack *pTrackXi = dynamic_cast<AliAODTrack*>( xi->GetDaughter(0) ); |
1826 | AliAODTrack *nTrackXi = dynamic_cast<AliAODTrack*>( xi->GetDaughter(1) ); | |
1827 | AliAODTrack *bachTrackXi = dynamic_cast<AliAODTrack*>( xi->GetDecayVertexXi()->GetDaughter(0) ); | |
1828 | if (!pTrackXi || !nTrackXi || !bachTrackXi ) { | |
1829 | AliWarning("ERROR: Could not retrieve one of the 3 AOD daughter tracks of the cascade ..."); | |
1830 | continue; | |
1831 | } | |
1832 | ||
1833 | UInt_t lIdxPosXi = (UInt_t) TMath::Abs( pTrackXi->GetID() ); | |
1834 | UInt_t lIdxNegXi = (UInt_t) TMath::Abs( nTrackXi->GetID() ); | |
1835 | UInt_t lBachIdx = (UInt_t) TMath::Abs( bachTrackXi->GetID() ); | |
1836 | ||
1837 | // Care track label can be negative in MC production (linked with the track quality) | |
1838 | // However = normally, not the case for track index ... | |
1839 | ||
1840 | // FIXME : rejection of a double use of a daughter track (nothing but just a crosscheck of what is done in the cascade vertexer) | |
1841 | if (lBachIdx == lIdxNegXi) { | |
1842 | AliWarning("Pb / Idx(Bach. track) = Idx(Neg. track) ... continue!"); continue; | |
1843 | } | |
1844 | if (lBachIdx == lIdxPosXi) { | |
1845 | AliWarning("Pb / Idx(Bach. track) = Idx(Pos. track) ... continue!"); continue; | |
1846 | } | |
1847 | ||
1848 | lPosTPCClusters = pTrackXi->GetTPCNcls(); // FIXME: Is this ok? or something like in LambdaK0PbPb task AOD? | |
1849 | lNegTPCClusters = nTrackXi->GetTPCNcls(); | |
1850 | lBachTPCClusters = bachTrackXi->GetTPCNcls(); | |
1851 | ||
1852 | // Rejection of a poor quality tracks | |
1853 | if (fkQualityCutTPCrefit) { | |
1854 | // 1 - Poor quality related to TPCrefit | |
1855 | if (!(pTrackXi->IsOn(AliAODTrack::kTPCrefit))) { AliWarning("Pb / V0 Pos. track has no TPCrefit ... continue!"); continue; } | |
1856 | if (!(nTrackXi->IsOn(AliAODTrack::kTPCrefit))) { AliWarning("Pb / V0 Neg. track has no TPCrefit ... continue!"); continue; } | |
1857 | if (!(bachTrackXi->IsOn(AliAODTrack::kTPCrefit))) { AliWarning("Pb / Bach. track has no TPCrefit ... continue!"); continue; } | |
1858 | } | |
1859 | if (fkQualityCut80TPCcls) { | |
1860 | // 2 - Poor quality related to TPC clusters | |
1861 | if (lPosTPCClusters < 80) { //AliWarning("Pb / V0 Pos. track has less than 80 TPC clusters ... continue!"); | |
1862 | continue; } | |
1863 | if (lNegTPCClusters < 80) { //AliWarning("Pb / V0 Neg. track has less than 80 TPC clusters ... continue!"); | |
1864 | continue; } | |
1865 | if (lBachTPCClusters < 80) { //AliWarning("Pb / Bach. track has less than 80 TPC clusters ... continue!"); | |
1866 | continue; } | |
1867 | } | |
1868 | ||
f9a6cab5 | 1869 | |
1870 | // - II.Step 3 : around the tracks : Bach + V0 (AOD) | |
1871 | // ~ Necessary variables for AODcascade data members coming from the AODv0 part (inheritance) | |
1872 | //------------- | |
1873 | ||
1874 | lChargeXi = xi->ChargeXi(); | |
1875 | ||
1876 | if ( lChargeXi < 0) | |
1877 | lInvMassLambdaAsCascDghter = xi->MassLambda(); | |
1878 | else | |
1879 | lInvMassLambdaAsCascDghter = xi->MassAntiLambda(); | |
1880 | ||
1881 | lDcaV0DaughtersXi = xi->DcaV0Daughters(); | |
1882 | lV0Chi2Xi = xi->Chi2V0(); | |
1883 | ||
f9a6cab5 | 1884 | lDcaV0ToPrimVertexXi = xi->DcaV0ToPrimVertex(); |
1885 | ||
1886 | lDcaBachToPrimVertexXi = xi->DcaBachToPrimVertex(); | |
1887 | ||
1888 | ||
1889 | lPosV0Xi[0] = xi->DecayVertexV0X(); | |
1890 | lPosV0Xi[1] = xi->DecayVertexV0Y(); | |
1891 | lPosV0Xi[2] = xi->DecayVertexV0Z(); | |
1892 | lV0RadiusXi = TMath::Sqrt( lPosV0Xi[0]*lPosV0Xi[0] + lPosV0Xi[1]*lPosV0Xi[1] ); | |
1893 | ||
39b18694 | 1894 | lV0CosineOfPointingAngleXi = xi->CosPointingAngle( lBestPrimaryVtxPos ); |
1895 | ||
f9a6cab5 | 1896 | lDcaPosToPrimVertexXi = xi->DcaPosToPrimVertex(); |
1897 | lDcaNegToPrimVertexXi = xi->DcaNegToPrimVertex(); | |
1898 | ||
1899 | // - II.Step 4 : around effective masses (AOD) | |
1900 | // ~ change mass hypotheses to cover all the possibilities : Xi-/+, Omega -/+ | |
1901 | //------------- | |
1902 | ||
1903 | if ( lChargeXi < 0 ) lInvMassXiMinus = xi->MassXi(); | |
1904 | if ( lChargeXi > 0 ) lInvMassXiPlus = xi->MassXi(); | |
1905 | if ( lChargeXi < 0 ) lInvMassOmegaMinus = xi->MassOmega(); | |
1906 | if ( lChargeXi > 0 ) lInvMassOmegaPlus = xi->MassOmega(); | |
1907 | ||
39b18694 | 1908 | |
1909 | // - II.Step 5 : PID on the daughters | |
f9a6cab5 | 1910 | //------------- |
1911 | ||
1912 | // Combined PID | |
1913 | ||
f9a6cab5 | 1914 | |
39b18694 | 1915 | // TPC PID: 4-sigma bands on Bethe-Bloch curve |
1916 | ||
1917 | // Bachelor | |
1918 | if (TMath::Abs(fPIDResponse->NumberOfSigmasTPC( bachTrackXi,AliPID::kKaon)) < 4) lIsBachelorKaonForTPC = kTRUE; | |
1919 | if (TMath::Abs(fPIDResponse->NumberOfSigmasTPC( bachTrackXi,AliPID::kPion)) < 4) lIsBachelorPionForTPC = kTRUE; | |
1920 | ||
1921 | // Negative V0 daughter | |
1922 | if (TMath::Abs(fPIDResponse->NumberOfSigmasTPC( nTrackXi,AliPID::kPion )) < 4) lIsNegPionForTPC = kTRUE; | |
1923 | if (TMath::Abs(fPIDResponse->NumberOfSigmasTPC( nTrackXi,AliPID::kProton )) < 4) lIsNegProtonForTPC = kTRUE; | |
1924 | ||
1925 | // Positive V0 daughter | |
1926 | if (TMath::Abs(fPIDResponse->NumberOfSigmasTPC( pTrackXi,AliPID::kPion )) < 4) lIsPosPionForTPC = kTRUE; | |
1927 | if (TMath::Abs(fPIDResponse->NumberOfSigmasTPC( pTrackXi,AliPID::kProton )) < 4) lIsPosProtonForTPC = kTRUE; | |
1928 | ||
f9a6cab5 | 1929 | |
1930 | // - II.Step 6 : extra info for QA (AOD) | |
1931 | // miscellaneous pieces onf info that may help regarding data quality assessment. | |
1932 | //------------- | |
1933 | lXiMomX = xi->MomXiX(); | |
1934 | lXiMomY = xi->MomXiY(); | |
1935 | lXiMomZ = xi->MomXiZ(); | |
1936 | lXiTransvMom = TMath::Sqrt( lXiMomX*lXiMomX + lXiMomY*lXiMomY ); | |
1937 | lXiTotMom = TMath::Sqrt( lXiMomX*lXiMomX + lXiMomY*lXiMomY + lXiMomZ*lXiMomZ ); | |
39b18694 | 1938 | |
1939 | Double_t lV0MomX = xi->MomV0X(); | |
1940 | Double_t lV0MomY = xi->MomV0Y(); | |
1941 | Double_t lV0MomZ = xi->MomV0Z(); | |
1942 | lV0TotMom = TMath::Sqrt(TMath::Power(lV0MomX,2)+TMath::Power(lV0MomY,2)+TMath::Power(lV0MomZ,2)); | |
1943 | ||
f9a6cab5 | 1944 | lBachMomX = xi->MomBachX(); |
1945 | lBachMomY = xi->MomBachY(); | |
1946 | lBachMomZ = xi->MomBachZ(); | |
1947 | lBachTransvMom = TMath::Sqrt( lBachMomX*lBachMomX + lBachMomY*lBachMomY ); | |
1948 | lBachTotMom = TMath::Sqrt( lBachMomX*lBachMomX + lBachMomY*lBachMomY + lBachMomZ*lBachMomZ ); | |
1949 | ||
1950 | ||
1951 | lV0toXiCosineOfPointingAngle = xi->CosPointingAngle( xi->GetDecayVertexXi() ); | |
1952 | ||
1953 | lRapXi = xi->RapXi(); | |
1954 | lRapOmega = xi->RapOmega(); | |
39b18694 | 1955 | lEta = xi->Eta(); // Will not work ! need a method Pz(), Py() Px() |
1956 | lTheta = xi->Theta() *180.0/TMath::Pi(); // in AODcascade. | |
1957 | lPhi = xi->Phi() *180.0/TMath::Pi(); // Here, we will get eta, theta, phi for the V0 ... | |
f9a6cab5 | 1958 | lAlphaXi = xi->AlphaXi(); |
1959 | lPtArmXi = xi->PtArmXi(); | |
1960 | ||
1961 | // II.Step 7 - Complementary info for monitoring the cascade cut variables | |
39b18694 | 1962 | // TPC clusters |
f9a6cab5 | 1963 | |
1964 | }// end of AOD treatment | |
1965 | ||
1966 | ||
39b18694 | 1967 | // Calculate proper time for cascade |
1968 | ||
1969 | Double_t cascadeMass = 0.; | |
1970 | ||
1971 | if ( ( (lChargeXi<0) && lIsBachelorPionForTPC && lIsPosProtonForTPC && lIsNegPionForTPC ) || | |
1972 | ( (lChargeXi>0) && lIsBachelorPionForTPC && lIsNegProtonForTPC && lIsPosPionForTPC ) ) cascadeMass = 1.321; | |
1973 | ||
1974 | if ( ( (lChargeXi<0) && lIsBachelorKaonForTPC && lIsPosProtonForTPC && lIsNegPionForTPC ) || | |
1975 | ( (lChargeXi>0) && lIsBachelorKaonForTPC && lIsNegProtonForTPC && lIsPosPionForTPC ) ) cascadeMass = 1.672; | |
1976 | ||
1977 | Double_t lctau = TMath::Sqrt(TMath::Power((lPosXi[0]-lBestPrimaryVtxPos[0]),2)+TMath::Power((lPosXi[1]-lBestPrimaryVtxPos[1]),2)+TMath::Power(( lPosXi[2]-lBestPrimaryVtxPos[2]),2)); | |
1978 | if (lXiTotMom!=0) lctau = lctau*cascadeMass/lXiTotMom; | |
1979 | else lctau = -1.; | |
1980 | ||
1981 | ||
1982 | // Calculate proper time for Lambda (reconstructed) | |
1983 | Float_t lambdaMass = 1.115683; // PDG mass | |
1984 | Float_t distV0Xi = TMath::Sqrt(TMath::Power((lPosV0Xi[0]-lPosXi[0]),2)+TMath::Power((lPosV0Xi[1]-lPosXi[1]),2)+TMath::Power((lPosV0Xi[2]-lPosXi[2]),2)); | |
1985 | Float_t lctauV0 = -1.; | |
1986 | if (lV0TotMom!=0) lctauV0 = distV0Xi*lambdaMass/lV0TotMom; | |
1987 | ||
1988 | Float_t distTV0Xi = TMath::Sqrt(TMath::Power((lPosV0Xi[0]-lPosXi[0]),2)+TMath::Power((lPosV0Xi[1]-lPosXi[1]),2)); | |
1989 | ||
1990 | ||
f9a6cab5 | 1991 | // ------------------------------------- |
1992 | // II.Fill - Filling the TH1,2,3Fs, HnSparses, CFContainers, FOR events with CASCADES ! | |
39b18694 | 1993 | // For AliEVE |
1994 | /* if(lChargeXi < 0&& lIsBachelorPionForTPC && lIsPosProtonForTPC && lIsNegPionForTPC ) { | |
1995 | if (lXiTransvMom>2.&&lXiTransvMom<4.&&(lInvMassXiMinus<1.322&&lInvMassXiMinus>1.320)&&(lXiRadius<8.&&lXiRadius>3.)) { | |
f9a6cab5 | 1996 | // // FIXME : Just to know which file is currently open : locate the file containing Xi |
39b18694 | 1997 | cout << "Name of the file containing Xi candidate(s) :" |
1998 | << CurrentFileName() | |
1999 | << " / entry: " << Entry() | |
2000 | << " / in file: " << lESDevent->GetEventNumberInFile() // <- Cvetan / From Mihaela: AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()->GetTree()->GetReadEntry(); | |
2001 | << " AliESDcascade number " << iXi | |
2002 | << " : mass(Xi-) = " << lInvMassXiMinus | |
2003 | << " / charge = " << lChargeXi | |
2004 | << " / pt(Casc) = " << lXiTransvMom | |
2005 | << " / Decay 2d R(Xi) = " << lXiRadius | |
2006 | << endl; | |
2007 | } | |
2008 | } | |
2009 | if(lChargeXi < 0&& lIsBachelorKaonForTPC && lIsPosProtonForTPC && lIsNegPionForTPC ) { | |
2010 | if (lXiTransvMom>2&&lXiTransvMom<4&&(lInvMassOmegaMinus<1.674&&lInvMassOmegaMinus>1.670)&&(lXiRadius<8.&&lXiRadius>3.)) { | |
2011 | cout << "Name of the file containing Omega candidate(s) :" | |
2012 | << CurrentFileName() | |
2013 | << " / entry: " << Entry() | |
2014 | << " / in file: " << lESDevent->GetEventNumberInFile() // <- Cvetan / From Mihaela: AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()->GetTree()->GetReadEntry(); | |
2015 | << " AliESDcascade number " << iXi | |
2016 | << " : mass(Omega-) = " << lInvMassOmegaMinus | |
2017 | << " / charge = " << lChargeXi | |
2018 | << " / pt(Casc) = " << lXiTransvMom | |
2019 | << " / Decay 2d R(Xi) = " << lXiRadius | |
2020 | << endl; | |
2021 | ||
2022 | } | |
2023 | } | |
2024 | */ | |
2025 | ||
f9a6cab5 | 2026 | |
2027 | // - II.Fill.Step 1 : primary vertex | |
2028 | ||
2029 | fHistTPCrefitTrackMultiplicityForCascadeEvt->Fill( nTrackWithTPCrefitMultiplicity ); | |
2030 | ||
2031 | fHistPosV0TPCClusters ->Fill( lPosTPCClusters ); | |
2032 | fHistNegV0TPCClusters ->Fill( lNegTPCClusters ); | |
2033 | fHistBachTPCClusters ->Fill( lBachTPCClusters ); | |
2034 | ||
2035 | f2dHistTPCdEdxOfCascDghters ->Fill( lInnerWallMomCascDghters[0] , lTPCSignalCascDghters[0] ); | |
2036 | f2dHistTPCdEdxOfCascDghters ->Fill( lInnerWallMomCascDghters[1] , lTPCSignalCascDghters[1] ); | |
2037 | f2dHistTPCdEdxOfCascDghters ->Fill( lInnerWallMomCascDghters[2] , lTPCSignalCascDghters[2] ); | |
2038 | ||
2039 | fHistVtxStatus ->Fill( lStatusTrackingPrimVtx ); // 1 if tracking vtx = ok | |
2040 | ||
2041 | if ( lStatusTrackingPrimVtx ) { | |
2042 | fHistPosTrkgPrimaryVtxXForCascadeEvt ->Fill( lTrkgPrimaryVtxPos[0] ); | |
2043 | fHistPosTrkgPrimaryVtxYForCascadeEvt ->Fill( lTrkgPrimaryVtxPos[1] ); | |
2044 | fHistPosTrkgPrimaryVtxZForCascadeEvt ->Fill( lTrkgPrimaryVtxPos[2] ); | |
2045 | fHistTrkgPrimaryVtxRadius ->Fill( lTrkgPrimaryVtxRadius3D ); | |
2046 | } | |
2047 | ||
2048 | fHistPosBestPrimaryVtxXForCascadeEvt ->Fill( lBestPrimaryVtxPos[0] ); | |
2049 | fHistPosBestPrimaryVtxYForCascadeEvt ->Fill( lBestPrimaryVtxPos[1] ); | |
2050 | fHistPosBestPrimaryVtxZForCascadeEvt ->Fill( lBestPrimaryVtxPos[2] ); | |
2051 | fHistBestPrimaryVtxRadius ->Fill( lBestPrimaryVtxRadius3D ); | |
2052 | ||
2053 | f2dHistTrkgPrimVtxVsBestPrimVtx->Fill( lTrkgPrimaryVtxRadius3D, lBestPrimaryVtxRadius3D ); | |
2054 | ||
2055 | // **************************** With PID on ? ... for the signal region ? ************FIXME************************************** | |
2056 | if ( ( (lChargeXi<0) && lIsBachelorPionForTPC && lIsPosProtonForTPC && lIsNegPionForTPC ) || | |
2057 | ( (lChargeXi>0) && lIsBachelorPionForTPC && lIsNegProtonForTPC && lIsPosPionForTPC ) ) | |
2058 | // NOTE : | |
2059 | // with this condition, it could happen that a cascade candidate satisfies the wrong requirement, | |
2060 | // e.g. one looks at a Xi- candidate for which lIsBachelorPionForTPC && lIsPosProtonForTPC && lIsNegPionForTPC = kFALSE | |
2061 | // Expectation: it should be excluded. | |
2062 | // but lIsBachelorPionForTPC && lIsNegProtonForTPC && lIsPosPionForTPC = kTRUE | |
2063 | // then this bad Xi-candidate will contribute anyway (OR condition). | |
2064 | // Hence : the extra condition on the sign of the Cascade | |
2065 | { | |
2066 | // if( TMath::Abs( lInvMassXiMinus-1.3217 ) < 0.010 || TMath::Abs( lInvMassXiPlus-1.3217 ) < 0.010) | |
2067 | ||
2068 | // II.Fill.Step 2 | |
2069 | fHistEffMassXi ->Fill( lEffMassXi ); | |
2070 | fHistChi2Xi ->Fill( lChi2Xi ); // Flag CascadeVtxer: Cut Variable a | |
2071 | fHistDcaXiDaughters ->Fill( lDcaXiDaughters ); // Flag CascadeVtxer: Cut Variable e | |
2072 | fHistDcaBachToPrimVertex ->Fill( lDcaBachToPrimVertexXi ); // Flag CascadeVtxer: Cut Variable d | |
2073 | fHistXiCosineOfPointingAngle ->Fill( lXiCosineOfPointingAngle ); // Flag CascadeVtxer: Cut Variable f | |
2074 | fHistXiRadius ->Fill( lXiRadius ); // Flag CascadeVtxer: Cut Variable g+h | |
2075 | ||
2076 | ||
2077 | // II.Fill.Step 3 | |
2078 | fHistMassLambdaAsCascDghter ->Fill( lInvMassLambdaAsCascDghter ); // Flag CascadeVtxer: Cut Variable c | |
2079 | fHistV0Chi2Xi ->Fill( lV0Chi2Xi ); | |
2080 | fHistDcaV0DaughtersXi ->Fill( lDcaV0DaughtersXi ); | |
2081 | fHistV0CosineOfPointingAngleXi ->Fill( lV0CosineOfPointingAngleXi ); | |
2082 | fHistV0RadiusXi ->Fill( lV0RadiusXi ); | |
2083 | ||
2084 | fHistDcaV0ToPrimVertexXi ->Fill( lDcaV0ToPrimVertexXi ); // Flag CascadeVtxer: Cut Variable b | |
2085 | fHistDcaPosToPrimVertexXi ->Fill( lDcaPosToPrimVertexXi ); | |
2086 | fHistDcaNegToPrimVertexXi ->Fill( lDcaNegToPrimVertexXi ); | |
2087 | ||
2088 | ||
2089 | // II.Fill.Step 4 : extra QA info | |
2090 | ||
2091 | fHistChargeXi ->Fill( lChargeXi ); | |
2092 | fHistV0toXiCosineOfPointingAngle->Fill( lV0toXiCosineOfPointingAngle ); | |
2093 | ||
2094 | if ( TMath::Abs( lInvMassXiMinus-1.3217 ) < 0.012 || TMath::Abs( lInvMassXiPlus-1.3217 ) < 0.012) {// One InvMass should be different from 0 | |
2095 | fHistXiTransvMom ->Fill( lXiTransvMom ); | |
2096 | fHistXiTotMom ->Fill( lXiTotMom ); | |
2097 | ||
2098 | fHistBachTransvMomXi ->Fill( lBachTransvMom ); | |
2099 | fHistBachTotMomXi ->Fill( lBachTotMom ); | |
2100 | ||
2101 | fHistRapXi ->Fill( lRapXi ); | |
2102 | fHistEtaXi ->Fill( lEta ); | |
2103 | fHistThetaXi ->Fill( lTheta ); | |
2104 | fHistPhiXi ->Fill( lPhi ); | |
2105 | } | |
2106 | ||
2107 | if ( TMath::Abs( lInvMassOmegaMinus-1.672 ) < 0.012 || TMath::Abs( lInvMassOmegaPlus-1.672 ) < 0.012 ) {// One InvMass should be different from 0 | |
2108 | fHistRapOmega ->Fill( lRapOmega ); | |
2109 | } | |
2110 | ||
2111 | f2dHistArmenteros ->Fill( lAlphaXi, lPtArmXi ); | |
2112 | }// end with PID ... | |
2113 | ||
2114 | // II.Fill.Step 5 : inv mass plots 1D | |
2115 | if ( lChargeXi < 0 ) { | |
2116 | fHistMassXiMinus ->Fill( lInvMassXiMinus ); | |
11bcd1e4 | 2117 | fHistMassOmegaMinus ->Fill( lInvMassOmegaMinus ); |
f9a6cab5 | 2118 | if(lIsBachelorPion) fHistMassWithCombPIDXiMinus ->Fill( lInvMassXiMinus ); |
2119 | if(lIsBachelorKaon) fHistMassWithCombPIDOmegaMinus ->Fill( lInvMassOmegaMinus ); | |
2120 | ||
2121 | fHistDcaXiDaughtersvsInvMass->Fill(lDcaXiDaughters,lInvMassXiMinus); | |
2122 | fHistDcaBachToPrimVertexvsInvMass->Fill(lDcaBachToPrimVertexXi,lInvMassXiMinus); | |
2123 | fHistXiCosineOfPointingAnglevsInvMass->Fill(lXiCosineOfPointingAngle,lInvMassXiMinus); | |
2124 | fHistMassLambdaAsCascDghtervsInvMass->Fill(lInvMassLambdaAsCascDghter,lInvMassXiMinus); | |
2125 | fHistDcaV0DaughtersXivsInvMass->Fill(lDcaV0DaughtersXi,lInvMassXiMinus); | |
2126 | fHistDcaV0ToPrimVertexXivsInvMass->Fill(lDcaV0ToPrimVertexXi,lInvMassXiMinus); | |
2127 | } | |
2128 | ||
2129 | if ( lChargeXi > 0 ) { | |
2130 | fHistMassXiPlus ->Fill( lInvMassXiPlus ); | |
11bcd1e4 | 2131 | fHistMassOmegaPlus ->Fill( lInvMassOmegaPlus ); |
f9a6cab5 | 2132 | if(lIsBachelorPion) fHistMassWithCombPIDXiPlus ->Fill( lInvMassXiPlus ); |
2133 | if(lIsBachelorKaon) fHistMassWithCombPIDOmegaPlus ->Fill( lInvMassOmegaPlus ); | |
2134 | } | |
2135 | ||
2136 | ||
2137 | // II.Fill.Step 6 : inv mass plots 2D, 3D | |
2138 | if ( lChargeXi < 0 ) { | |
2139 | f2dHistEffMassLambdaVsEffMassXiMinus->Fill( lInvMassLambdaAsCascDghter, lInvMassXiMinus ); | |
2140 | f2dHistEffMassXiVsEffMassOmegaMinus ->Fill( lInvMassXiMinus, lInvMassOmegaMinus ); | |
2141 | f2dHistXiRadiusVsEffMassXiMinus ->Fill( lXiRadius, lInvMassXiMinus ); | |
2142 | f2dHistXiRadiusVsEffMassOmegaMinus ->Fill( lXiRadius, lInvMassOmegaMinus ); | |
2143 | // f3dHistXiPtVsEffMassVsYXiMinus ->Fill( lXiTransvMom, lInvMassXiMinus, lRapXi ); | |
2144 | // f3dHistXiPtVsEffMassVsYOmegaMinus ->Fill( lXiTransvMom, lInvMassOmegaMinus, lRapOmega ); | |
2145 | } else { | |
2146 | f2dHistEffMassLambdaVsEffMassXiPlus ->Fill( lInvMassLambdaAsCascDghter, lInvMassXiPlus ); | |
2147 | f2dHistEffMassXiVsEffMassOmegaPlus ->Fill( lInvMassXiPlus, lInvMassOmegaPlus ); | |
2148 | f2dHistXiRadiusVsEffMassXiPlus ->Fill( lXiRadius, lInvMassXiPlus); | |
2149 | f2dHistXiRadiusVsEffMassOmegaPlus ->Fill( lXiRadius, lInvMassOmegaPlus ); | |
2150 | // f3dHistXiPtVsEffMassVsYXiPlus ->Fill( lXiTransvMom, lInvMassXiPlus, lRapXi ); | |
2151 | // f3dHistXiPtVsEffMassVsYOmegaPlus ->Fill( lXiTransvMom, lInvMassOmegaPlus, lRapOmega ); | |
2152 | } | |
2153 | ||
2154 | // - Filling the AliCFContainers related to PID | |
2155 | ||
2156 | Double_t lContainerPIDVars[4] = {0.0}; | |
2157 | ||
2158 | // Xi Minus | |
2159 | if ( lChargeXi < 0 ) { | |
2160 | lContainerPIDVars[0] = lXiTransvMom ; | |
2161 | lContainerPIDVars[1] = lInvMassXiMinus ; | |
2162 | lContainerPIDVars[2] = lRapXi ; | |
39b18694 | 2163 | lContainerPIDVars[3] = lcentrality; |
f9a6cab5 | 2164 | |
2165 | // No PID | |
2166 | fCFContCascadePIDXiMinus->Fill(lContainerPIDVars, 0); // No PID | |
2167 | // TPC PID | |
2168 | if ( lIsBachelorPionForTPC ) | |
2169 | fCFContCascadePIDXiMinus->Fill(lContainerPIDVars, 1); // TPC PID / 4-#sigma cut on Bachelor track | |
2170 | ||
2171 | if ( lIsBachelorPionForTPC && | |
2172 | lIsPosProtonForTPC ) | |
2173 | fCFContCascadePIDXiMinus->Fill(lContainerPIDVars, 2); // TPC PID / 4-#sigma cut on Bachelor+Baryon tracks | |
2174 | ||
2175 | if ( lIsBachelorPionForTPC && | |
2176 | lIsPosProtonForTPC && | |
2177 | lIsNegPionForTPC ) | |
2178 | fCFContCascadePIDXiMinus->Fill(lContainerPIDVars, 3); // TPC PID / 4-#sigma cut on Bachelor+Baryon+Meson tracks | |
2179 | ||
2180 | // Combined PID | |
2181 | if ( lIsBachelorPion ) | |
2182 | fCFContCascadePIDXiMinus->Fill(lContainerPIDVars, 4); // Comb. PID / Bachelor | |
2183 | ||
2184 | if ( lIsBachelorPion && | |
2185 | lIsPosInXiProton ) | |
2186 | fCFContCascadePIDXiMinus->Fill(lContainerPIDVars, 5); // Comb. PID / Bachelor+Baryon | |
2187 | ||
2188 | if (lIsBachelorPion && | |
2189 | lIsPosInXiProton && | |
2190 | lIsNegInXiPion ) | |
2191 | fCFContCascadePIDXiMinus->Fill(lContainerPIDVars, 6); // Comb. PID / Bachelor+Baryon+Meson | |
2192 | } | |
2193 | ||
2194 | lContainerPIDVars[0] = 0.; lContainerPIDVars[1] = 0.; lContainerPIDVars[2] = 0.; lContainerPIDVars[3] = 0.; | |
2195 | ||
2196 | // Xi Plus | |
2197 | if ( lChargeXi > 0 ) { | |
2198 | lContainerPIDVars[0] = lXiTransvMom ; | |
2199 | lContainerPIDVars[1] = lInvMassXiPlus ; | |
2200 | lContainerPIDVars[2] = lRapXi ; | |
39b18694 | 2201 | lContainerPIDVars[3] = lcentrality; |
f9a6cab5 | 2202 | |
2203 | // No PID | |
2204 | fCFContCascadePIDXiPlus->Fill(lContainerPIDVars, 0); // No PID | |
2205 | // TPC PID | |
2206 | if ( lIsBachelorPionForTPC ) | |
2207 | fCFContCascadePIDXiPlus->Fill(lContainerPIDVars, 1); // TPC PID / 4-#sigma cut on Bachelor track | |
2208 | ||
2209 | if ( lIsBachelorPionForTPC && | |
2210 | lIsNegProtonForTPC ) | |
2211 | fCFContCascadePIDXiPlus->Fill(lContainerPIDVars, 2); // TPC PID / 4-#sigma cut on Bachelor+Baryon tracks | |
2212 | ||
2213 | if ( lIsBachelorPionForTPC && | |
2214 | lIsNegProtonForTPC && | |
2215 | lIsPosPionForTPC ) | |
2216 | fCFContCascadePIDXiPlus->Fill(lContainerPIDVars, 3); // TPC PID / 4-#sigma cut on Bachelor+Baryon+Meson tracks | |
2217 | ||
2218 | // Combined PID | |
2219 | if ( lIsBachelorPion ) | |
2220 | fCFContCascadePIDXiPlus->Fill(lContainerPIDVars, 4); // Comb. PID / Bachelor | |
2221 | ||
2222 | if ( lIsBachelorPion && | |
2223 | lIsNegInXiProton ) | |
2224 | fCFContCascadePIDXiPlus->Fill(lContainerPIDVars, 5); // Comb. PID / Bachelor+Baryon | |
2225 | ||
2226 | if (lIsBachelorPion && | |
2227 | lIsNegInXiProton && | |
2228 | lIsPosInXiPion ) | |
2229 | fCFContCascadePIDXiPlus->Fill(lContainerPIDVars, 6); // Comb. PID / Bachelor+Baryon+Meson | |
2230 | } | |
2231 | ||
2232 | lContainerPIDVars[0] = 0.; lContainerPIDVars[1] = 0.; lContainerPIDVars[2] = 0.; lContainerPIDVars[3] = 0.; | |
2233 | ||
2234 | // Omega Minus | |
2235 | if ( lChargeXi < 0 ) { | |
2236 | lContainerPIDVars[0] = lXiTransvMom ; | |
2237 | lContainerPIDVars[1] = lInvMassOmegaMinus ; | |
2238 | lContainerPIDVars[2] = lRapOmega ; | |
39b18694 | 2239 | lContainerPIDVars[3] = lcentrality; |
11bcd1e4 | 2240 | |
f9a6cab5 | 2241 | // No PID |
2242 | fCFContCascadePIDOmegaMinus->Fill(lContainerPIDVars, 0); // No PID | |
2243 | // TPC PID | |
2244 | if ( lIsBachelorKaonForTPC ) | |
2245 | fCFContCascadePIDOmegaMinus->Fill(lContainerPIDVars, 1); // TPC PID / 4-#sigma cut on Bachelor track | |
2246 | ||
2247 | if ( lIsBachelorKaonForTPC && | |
2248 | lIsPosProtonForTPC ) | |
2249 | fCFContCascadePIDOmegaMinus->Fill(lContainerPIDVars, 2); // TPC PID / 4-#sigma cut on Bachelor+Baryon tracks | |
2250 | ||
2251 | if ( lIsBachelorKaonForTPC && | |
2252 | lIsPosProtonForTPC && | |
2253 | lIsNegPionForTPC ) | |
2254 | fCFContCascadePIDOmegaMinus->Fill(lContainerPIDVars, 3); // TPC PID / 4-#sigma cut on Bachelor+Baryon+Meson tracks | |
2255 | ||
2256 | // Combined PID | |
2257 | if ( lIsBachelorKaon ) | |
2258 | fCFContCascadePIDOmegaMinus->Fill(lContainerPIDVars, 4); // Comb. PID / Bachelor | |
2259 | ||
2260 | if ( lIsBachelorKaon && | |
2261 | lIsPosInOmegaProton ) | |
2262 | fCFContCascadePIDOmegaMinus->Fill(lContainerPIDVars, 5); // Comb. PID / Bachelor+Baryon | |
2263 | ||
2264 | if (lIsBachelorKaon && | |
2265 | lIsPosInOmegaProton && | |
2266 | lIsNegInOmegaPion ) | |
2267 | fCFContCascadePIDOmegaMinus->Fill(lContainerPIDVars, 6); // Comb. PID / Bachelor+Baryon+Meson | |
2268 | } | |
2269 | ||
2270 | lContainerPIDVars[0] = 0.; lContainerPIDVars[1] = 0.; lContainerPIDVars[2] = 0.; lContainerPIDVars[3] = 0.; | |
2271 | ||
2272 | // Omega Plus | |
2273 | if ( lChargeXi > 0 ) { | |
2274 | lContainerPIDVars[0] = lXiTransvMom ; | |
2275 | lContainerPIDVars[1] = lInvMassOmegaPlus ; | |
2276 | lContainerPIDVars[2] = lRapOmega ; | |
39b18694 | 2277 | lContainerPIDVars[3] = lcentrality; |
f9a6cab5 | 2278 | |
2279 | // No PID | |
2280 | fCFContCascadePIDOmegaPlus->Fill(lContainerPIDVars, 0); // No PID | |
2281 | // TPC PID | |
2282 | if ( lIsBachelorKaonForTPC ) | |
2283 | fCFContCascadePIDOmegaPlus->Fill(lContainerPIDVars, 1); // TPC PID / 4-#sigma cut on Bachelor track | |
2284 | ||
2285 | if( lIsBachelorKaonForTPC && | |
2286 | lIsNegProtonForTPC ) | |
2287 | fCFContCascadePIDOmegaPlus->Fill(lContainerPIDVars, 2); // TPC PID / 4-#sigma cut on Bachelor+Baryon tracks | |
2288 | ||
2289 | if ( lIsBachelorKaonForTPC && | |
2290 | lIsNegProtonForTPC && | |
2291 | lIsPosPionForTPC ) | |
2292 | fCFContCascadePIDOmegaPlus->Fill(lContainerPIDVars, 3); // TPC PID / 4-#sigma cut on Bachelor+Baryon+Meson tracks | |
2293 | ||
2294 | // Combined PID | |
2295 | if ( lIsBachelorKaon ) | |
2296 | fCFContCascadePIDOmegaPlus->Fill(lContainerPIDVars, 4); // Comb. PID / Bachelor | |
2297 | ||
2298 | if ( lIsBachelorKaon && | |
2299 | lIsNegInOmegaProton ) | |
2300 | fCFContCascadePIDOmegaPlus->Fill(lContainerPIDVars, 5); // Comb. PID / Bachelor+Baryon | |
2301 | ||
2302 | if (lIsBachelorKaon && | |
2303 | lIsNegInOmegaProton && | |
2304 | lIsPosInOmegaPion ) | |
2305 | fCFContCascadePIDOmegaPlus->Fill(lContainerPIDVars, 6); // Comb. PID / Bachelor+Baryon+Meson | |
2306 | } | |
2307 | ||
2308 | ||
2309 | // II.Fill.Step 7 : filling the AliCFContainer (optimisation of topological selections) | |
39b18694 | 2310 | Double_t lContainerCutVars[22] = {0.0}; |
f9a6cab5 | 2311 | |
2312 | lContainerCutVars[0] = lDcaXiDaughters; | |
2313 | lContainerCutVars[1] = lDcaBachToPrimVertexXi; | |
2314 | lContainerCutVars[2] = lXiCosineOfPointingAngle; | |
2315 | lContainerCutVars[3] = lXiRadius; | |
2316 | lContainerCutVars[4] = lInvMassLambdaAsCascDghter; | |
2317 | lContainerCutVars[5] = lDcaV0DaughtersXi; | |
39b18694 | 2318 | lContainerCutVars[6] = lV0toXiCosineOfPointingAngle; |
f9a6cab5 | 2319 | lContainerCutVars[7] = lV0RadiusXi; |
2320 | lContainerCutVars[8] = lDcaV0ToPrimVertexXi; | |
2321 | lContainerCutVars[9] = lDcaPosToPrimVertexXi; | |
2322 | lContainerCutVars[10] = lDcaNegToPrimVertexXi; | |
2323 | ||
2324 | lContainerCutVars[13] = lXiTransvMom; | |
2325 | ||
2326 | lContainerCutVars[16] = lBestPrimaryVtxPos[2]; | |
39b18694 | 2327 | lContainerCutVars[17] = lcentrality; |
2328 | lContainerCutVars[18] = lPrimaryTrackMultiplicity; | |
2329 | // lContainerCutVars[19] = lBachTPCClusters; | |
2330 | ||
2331 | lContainerCutVars[19] = lctau; | |
2332 | lContainerCutVars[20] = lctauV0; | |
2333 | lContainerCutVars[21] = distTV0Xi; | |
2334 | ||
f9a6cab5 | 2335 | if ( lChargeXi < 0 ) { |
2336 | lContainerCutVars[11] = lInvMassXiMinus; | |
11bcd1e4 | 2337 | lContainerCutVars[12] = lInvMassOmegaMinus;//1.63; |
f9a6cab5 | 2338 | lContainerCutVars[14] = lRapXi; |
2339 | lContainerCutVars[15] = -1.; | |
11bcd1e4 | 2340 | if ( lIsBachelorPionForTPC && lIsPosProtonForTPC && lIsNegPionForTPC ) |
2341 | fCFContCascadeCuts->Fill(lContainerCutVars,0); // for Xi- | |
f9a6cab5 | 2342 | |
11bcd1e4 | 2343 | lContainerCutVars[11] = lInvMassXiMinus; |
f9a6cab5 | 2344 | lContainerCutVars[12] = lInvMassOmegaMinus; |
2345 | lContainerCutVars[14] = -1.; | |
2346 | lContainerCutVars[15] = lRapOmega; | |
11bcd1e4 | 2347 | if ( lIsBachelorKaonForTPC && lIsPosProtonForTPC && lIsNegPionForTPC ) |
2348 | fCFContCascadeCuts->Fill(lContainerCutVars,2); // for Omega- | |
f9a6cab5 | 2349 | } else { |
2350 | lContainerCutVars[11] = lInvMassXiPlus; | |
11bcd1e4 | 2351 | lContainerCutVars[12] = lInvMassOmegaPlus;//1.63; |
f9a6cab5 | 2352 | lContainerCutVars[14] = lRapXi; |
2353 | lContainerCutVars[15] = -1.; | |
11bcd1e4 | 2354 | if ( lIsBachelorPionForTPC && lIsNegProtonForTPC && lIsPosPionForTPC ) |
2355 | fCFContCascadeCuts->Fill(lContainerCutVars,1); // for Xi+ | |
f9a6cab5 | 2356 | |
11bcd1e4 | 2357 | lContainerCutVars[11] = lInvMassXiPlus;//1.26; |
f9a6cab5 | 2358 | lContainerCutVars[12] = lInvMassOmegaPlus; |
2359 | lContainerCutVars[14] = -1.; | |
2360 | lContainerCutVars[15] = lRapOmega; | |
11bcd1e4 | 2361 | if ( lIsBachelorKaonForTPC && lIsNegProtonForTPC && lIsPosPionForTPC ) |
2362 | fCFContCascadeCuts->Fill(lContainerCutVars,3); // for Omega+ | |
f9a6cab5 | 2363 | } |
f9a6cab5 | 2364 | |
f9a6cab5 | 2365 | }// end of the Cascade loop (ESD or AOD) |
2366 | ||
2367 | ||
2368 | // Post output data. | |
2369 | PostData(1, fListHistCascade); | |
2370 | } | |
2371 | ||
39b18694 | 2372 | //________________________________________________________________________ |
f9a6cab5 | 2373 | Int_t AliAnalysisTaskCheckCascadePbPb::DoESDTrackWithTPCrefitMultiplicity(const AliESDEvent *lESDevent) { |
2374 | // Checking the number of tracks with TPCrefit for each event | |
2375 | // Needed for a rough assessment of the event multiplicity | |
2376 | ||
2377 | Int_t nTrackWithTPCrefitMultiplicity = 0; | |
2378 | for (Int_t iTrackIdx = 0; iTrackIdx < (InputEvent())->GetNumberOfTracks(); iTrackIdx++) { | |
2379 | AliESDtrack *esdTrack = 0x0; | |
2380 | esdTrack = lESDevent->GetTrack( iTrackIdx ); | |
2381 | if (!esdTrack) { AliWarning("Pb / Could not retrieve one track within the track loop for TPCrefit check ..."); continue; } | |
2382 | ||
2383 | ULong_t lTrackStatus = esdTrack->GetStatus(); | |
2384 | if ((lTrackStatus&AliESDtrack::kTPCrefit) == 0) continue; | |
2385 | else nTrackWithTPCrefitMultiplicity++; | |
2386 | // FIXME : | |
2387 | // The goal here is to get a better assessment of the event multiplicity. | |
2388 | // (InputEvent())->GetNumberOfTracks() takes into account ITS std alone tracks + global tracks | |
2389 | // This may introduce a bias. Hence the number of TPC refit tracks. | |
2390 | // Note : the event multiplicity = analysis on its own... See Jacek's or Jan Fiete's analysis on dN/d(eta) | |
2391 | ||
2392 | }// end loop over all event tracks | |
2393 | return nTrackWithTPCrefitMultiplicity; | |
2394 | } | |
2395 | ||
f9a6cab5 | 2396 | //________________________________________________________________________ |
2397 | void AliAnalysisTaskCheckCascadePbPb::Terminate(Option_t *) | |
2398 | { | |
2399 | // Draw result to the screen | |
2400 | // Called once at the end of the query | |
2401 | ||
2402 | TList *cRetrievedList = 0x0; | |
2403 | cRetrievedList = (TList*)GetOutputData(1); | |
2404 | if(!cRetrievedList){ | |
2405 | AliWarning("ERROR - AliAnalysisTaskCheckCascadePbPb: ouput data container list not available\n"); return; | |
2406 | } | |
2407 | ||
2408 | fHistTrackMultiplicityForCentrEvt = dynamic_cast<TH1F*> ( cRetrievedList->FindObject("fHistTrackMultiplicityForCentrEvt") ); | |
2409 | if (!fHistTrackMultiplicityForCentrEvt) { | |
2410 | AliWarning("ERROR - AliAnalysisTaskCheckCascadePbPb: fHistTrackMultiplicityForCentrEvt not available\n"); return; | |
2411 | } | |
2412 | ||
2413 | fHistCascadeMultiplicityForCentrEvt = dynamic_cast<TH1F*> ( cRetrievedList->FindObject("fHistCascadeMultiplicityForCentrEvt") ); | |
2414 | if (!fHistCascadeMultiplicityForCentrEvt) { | |
2415 | AliWarning("ERROR - AliAnalysisTaskCheckCascadePbPb: fHistCascadeMultiplicityForCentrEvt not available\n"); return; | |
2416 | } | |
2417 | ||
2418 | fHistMassXiMinus = dynamic_cast<TH1F*> ( cRetrievedList->FindObject("fHistMassXiMinus") ); | |
2419 | if (!fHistMassXiMinus) { | |
2420 | AliWarning("ERROR - AliAnalysisTaskCheckCascadePbPb: fHistMassXiMinus not available\n"); return; | |
2421 | } | |
2422 | fHistMassXiPlus = dynamic_cast<TH1F*> ( cRetrievedList->FindObject("fHistMassXiPlus") ); | |
2423 | if (!fHistMassXiPlus) { | |
2424 | AliWarning("ERROR - AliAnalysisTaskCheckCascadePbPb: fHistMassXiPlus not available\n"); return; | |
2425 | } | |
2426 | fHistMassOmegaMinus = dynamic_cast<TH1F*> ( cRetrievedList->FindObject("fHistMassOmegaMinus") ); | |
2427 | if (!fHistMassOmegaMinus) { | |
2428 | AliWarning("ERROR - AliAnalysisTaskCheckCascadePbPb: fHistMassOmegaMinus not available\n"); return; | |
2429 | } | |
2430 | fHistMassOmegaPlus = dynamic_cast<TH1F*> ( cRetrievedList->FindObject("fHistMassOmegaPlus") ); | |
2431 | if (!fHistMassOmegaPlus) { | |
2432 | AliWarning("ERROR - AliAnalysisTaskCheckCascadePbPb: fHistMassOmegaPlus not available\n"); return; | |
2433 | } | |
2434 | ||
2435 | TCanvas *canCheckCascade = new TCanvas("AliAnalysisTaskCheckCascadePbPb","CheckCascade overview",10,10,1010,660); | |
2436 | canCheckCascade->Divide(2,2); | |
2437 | ||
2438 | canCheckCascade->cd(1); | |
2439 | canCheckCascade->cd(1)->SetLogy(); | |
2440 | fHistTrackMultiplicityForCentrEvt->SetMarkerStyle(kFullStar); | |
2441 | fHistTrackMultiplicityForCentrEvt->GetXaxis()->SetLabelFont(42); | |
2442 | fHistTrackMultiplicityForCentrEvt->GetYaxis()->SetLabelFont(42); | |
2443 | fHistTrackMultiplicityForCentrEvt->SetTitleFont(42, "xy"); | |
2444 | fHistTrackMultiplicityForCentrEvt->GetXaxis()->SetTitleOffset(1.1); | |
2445 | fHistTrackMultiplicityForCentrEvt->DrawCopy("H"); | |
2446 | ||
2447 | canCheckCascade->cd(2); | |
2448 | canCheckCascade->cd(2)->SetLogy(); | |
2449 | fHistCascadeMultiplicityForCentrEvt->SetMarkerStyle(kOpenSquare); | |
2450 | fHistCascadeMultiplicityForCentrEvt->GetXaxis()->SetLabelFont(42); | |
2451 | fHistCascadeMultiplicityForCentrEvt->GetYaxis()->SetLabelFont(42); | |
2452 | fHistCascadeMultiplicityForCentrEvt->SetTitleFont(42, "xy"); | |
2453 | fHistCascadeMultiplicityForCentrEvt->GetXaxis()->SetTitleOffset(1.1); | |
2454 | fHistCascadeMultiplicityForCentrEvt->DrawCopy("E"); | |
2455 | ||
2456 | canCheckCascade->cd(3); | |
2457 | fHistMassXiMinus ->SetMarkerStyle(kFullCircle); | |
2458 | fHistMassXiMinus ->SetMarkerSize(0.5); | |
2459 | fHistMassXiMinus ->GetXaxis()->SetLabelFont(42); | |
2460 | fHistMassXiMinus ->GetYaxis()->SetLabelFont(42); | |
2461 | fHistMassXiMinus ->SetTitleFont(42, "xy"); | |
2462 | fHistMassXiMinus ->GetXaxis()->SetTitleOffset(1.1); | |
2463 | fHistMassXiMinus ->GetYaxis()->SetTitleOffset(1.3); | |
2464 | // fHistMassXiMinus->Rebin(2); | |
2465 | fHistMassXiMinus ->GetXaxis()->SetRangeUser(1.24, 1.42); | |
2466 | fHistMassXiMinus ->DrawCopy("E"); | |
2467 | ||
2468 | fHistMassXiPlus ->SetMarkerStyle(kOpenCircle); | |
2469 | fHistMassXiPlus ->SetMarkerColor(kRed+2); | |
2470 | fHistMassXiPlus ->SetLineColor(kRed+2); | |
2471 | fHistMassXiPlus ->SetMarkerSize(0.5); | |
2472 | // fHistMassXiPlus ->Rebin(2); | |
2473 | fHistMassXiPlus ->DrawCopy("ESAME"); | |
2474 | ||
2475 | ||
39b18694 | 2476 | TLegend *legendXi =new TLegend(0.67,0.34,0.97,0.54); |
2477 | legendXi->SetTextFont(42); | |
2478 | legendXi->SetTextSize(0.05); | |
2479 | legendXi->SetFillColor(kWhite); | |
2480 | legendXi->AddEntry( fHistMassXiMinus,"#Xi^{-} candidates","lp"); | |
2481 | legendXi->AddEntry( fHistMassXiPlus,"#Xi^{+} candidates","lp"); | |
2482 | legendXi->Draw(); | |
f9a6cab5 | 2483 | |
2484 | ||
2485 | canCheckCascade->cd(4); | |
2486 | fHistMassOmegaPlus ->SetMarkerStyle(kOpenCircle); | |
2487 | fHistMassOmegaPlus ->SetMarkerColor(kRed+2); | |
2488 | fHistMassOmegaPlus ->SetLineColor(kRed+2); | |
2489 | fHistMassOmegaPlus ->SetMarkerSize(0.5); | |
2490 | fHistMassOmegaPlus ->GetXaxis()->SetLabelFont(42); | |
2491 | fHistMassOmegaPlus ->GetYaxis()->SetLabelFont(42); | |
2492 | fHistMassOmegaPlus ->SetTitleFont(42, "xy"); | |
2493 | fHistMassOmegaPlus ->GetXaxis()->SetTitleOffset(1.1); | |
2494 | fHistMassOmegaPlus ->GetYaxis()->SetTitleOffset(1.25); | |
2495 | // fHistMassOmegaPlus ->Rebin(2); | |
2496 | fHistMassOmegaPlus ->GetXaxis()->SetRangeUser(1.6, 1.84); | |
2497 | fHistMassOmegaPlus ->DrawCopy("E"); | |
2498 | ||
2499 | fHistMassOmegaMinus->SetMarkerStyle(kFullCircle); | |
2500 | fHistMassOmegaMinus->SetMarkerSize(0.5); | |
2501 | // fHistMassOmegaMinus->Rebin(2); | |
2502 | fHistMassOmegaMinus->DrawCopy("ESAME"); | |
2503 | ||
2504 | ||
39b18694 | 2505 | TLegend *legendOmega = new TLegend(0.67,0.34,0.97,0.54); |
2506 | legendOmega->SetTextFont(42); | |
2507 | legendOmega->SetTextSize(0.05); | |
2508 | legendOmega->SetFillColor(kWhite); | |
2509 | legendOmega->AddEntry( fHistMassOmegaMinus,"#Omega^{-} candidates","lp"); | |
2510 | legendOmega->AddEntry( fHistMassOmegaPlus,"#Omega^{+} candidates","lp"); | |
2511 | legendOmega->Draw(); | |
f9a6cab5 | 2512 | |
2513 | } |