]>
Commit | Line | Data |
---|---|---|
45ee0bcc | 1 | |
2 | /************************************************************************** | |
3 | * Authors : Antonin Maire, Boris Hippolyte * | |
4 | * Contributors are mentioned in the code where appropriate. * | |
5 | * * | |
6 | * Permission to use, copy, modify and distribute this software and its * | |
7 | * documentation strictly for non-commercial purposes is hereby granted * | |
8 | * without fee, provided that the above copyright notice appears in all * | |
9 | * copies and that both the copyright notice and this permission notice * | |
10 | * appear in the supporting documentation. The authors make no claims * | |
11 | * about the suitability of this software for any purpose. It is * | |
12 | * provided "as is" without express or implied warranty. * | |
13 | **************************************************************************/ | |
14 | ||
15 | //----------------------------------------------------------------- | |
16 | // AliAnalysisTaskCheckCascade class | |
3873f6ee | 17 | // (AliAnalysisTaskCheckCascade) |
32599cd4 | 18 | // This task has four roles : |
3873f6ee | 19 | // 1. QAing the Cascades from ESD and AOD |
20 | // Origin: AliAnalysisTaskESDCheckV0 by B.H. Nov2007, hippolyt@in2p3.fr | |
32599cd4 | 21 | // 2. Prepare the plots which stand as raw material for yield extraction (wi/wo PID) |
22 | // 3. Supply an AliCFContainer meant to define the optimised topological selections | |
23 | // 4. Rough azimuthal correlation study (Eta, Phi) | |
3873f6ee | 24 | // Adapted to Cascade : A.Maire Mar2008, antonin.maire@ires.in2p3.fr |
c525f360 | 25 | // Modified : A.Maire Nov2010, antonin.maire@ires.in2p3.fr |
45ee0bcc | 26 | //----------------------------------------------------------------- |
27 | ||
28 | ||
29 | ||
45ee0bcc | 30 | class TTree; |
31 | class TParticle; | |
45ee0bcc | 32 | |
33 | //class AliMCEventHandler; | |
34 | //class AliMCEvent; | |
35 | //class AliStack; | |
36 | ||
37 | class AliESDVertex; | |
38 | class AliAODVertex; | |
39 | class AliESDv0; | |
40 | class AliAODv0; | |
41 | ||
f87cd3db | 42 | #include <Riostream.h> |
45ee0bcc | 43 | #include "TList.h" |
6dc8b253 | 44 | #include "TH1.h" |
45 | #include "TH2.h" | |
ff0753e1 | 46 | #include "TH3.h" |
1bf66ceb | 47 | //#include "THnSparse.h" |
3873f6ee | 48 | #include "TVector3.h" |
45ee0bcc | 49 | #include "TCanvas.h" |
50 | #include "TMath.h" | |
f3dc9369 | 51 | #include "TLegend.h" |
45ee0bcc | 52 | |
3873f6ee | 53 | |
45ee0bcc | 54 | #include "AliLog.h" |
55 | ||
56 | #include "AliESDEvent.h" | |
57 | #include "AliAODEvent.h" | |
c525f360 | 58 | // #include "AliV0vertexer.h" |
59 | // #include "AliCascadeVertexer.h" | |
10d100d4 | 60 | #include "AliESDpid.h" |
c525f360 | 61 | #include "AliESDtrackCuts.h" |
7e504402 | 62 | |
63 | #include "AliInputEventHandler.h" | |
64 | #include "AliAnalysisManager.h" | |
90926f02 | 65 | #include "AliMCEventHandler.h" |
7e504402 | 66 | |
32599cd4 | 67 | #include "AliCFContainer.h" |
68 | #include "AliMultiplicity.h" | |
45ee0bcc | 69 | |
70 | #include "AliESDcascade.h" | |
71 | #include "AliAODcascade.h" | |
72 | ||
73 | #include "AliAnalysisTaskCheckCascade.h" | |
74 | ||
75 | ClassImp(AliAnalysisTaskCheckCascade) | |
76 | ||
77 | ||
78 | ||
79 | //________________________________________________________________________ | |
80 | AliAnalysisTaskCheckCascade::AliAnalysisTaskCheckCascade() | |
c525f360 | 81 | : AliAnalysisTaskSE(), fAnalysisType("ESD"), fTriggerMaskType("kMB"), fCollidingSystems(0), fESDpid(0), fESDtrackCuts(0), /*fPaveTextBookKeeping(0),*/ |
82 | fkRerunV0CascVertexers (0), | |
83 | fkQualityCutZprimVtxPos (kTRUE), | |
84 | fkRejectEventPileUp (kTRUE), | |
85 | fkQualityCutNoTPConlyPrimVtx (kTRUE), | |
86 | fkQualityCutTPCrefit (kTRUE), | |
87 | fkQualityCut80TPCcls (kTRUE), | |
88 | fkIsDataRecoWith1PadTPCCluster (kTRUE), | |
89 | fkExtraSelections (0), | |
90 | fAngularCorrelationType ("TrigLeadingTrck-AssoCasc"), | |
7e504402 | 91 | |
45ee0bcc | 92 | // - Cascade part initialisation |
93 | fListHistCascade(0), | |
9d670198 | 94 | fHistCascadeMultiplicityBeforeTrigSel(0), |
a786bd57 | 95 | fHistCascadeMultiplicityForTrigEvt(0), fHistTrackMultiplicityForTrigEvt(0), fHistTPCrefitTrackMultiplicityForTrigEvt(0), fHistPrimaryTrackMultiplicityForTrigEvt(0), |
96 | fHistEstimateITSTPCMultiplicityForTrigEvt(0), | |
c525f360 | 97 | fHistCascadeMultiplicityForTrigEvtAndZprimVtx(0), fHistCascadeMultiplicityForTrigEvtNonPiledUpAndZprimVtx(0), |
9d670198 | 98 | fHistCascadeMultiplicityForSelEvt(0), |
99 | fHistPosBestPrimaryVtxXForSelEvt(0), fHistPosBestPrimaryVtxYForSelEvt(0), fHistPosBestPrimaryVtxZForSelEvt(0), | |
c525f360 | 100 | fHistTPCrefitTrackMultiplicityForCascadeEvt(0), fHistPrimaryTrackMultiplicityForCascadeEvt(0), |
7e504402 | 101 | fHistPosV0TPCClusters(0), fHistNegV0TPCClusters(0), fHistBachTPCClusters(0), |
45ee0bcc | 102 | fHistVtxStatus(0), |
103 | ||
9d670198 | 104 | fHistPosTrkgPrimaryVtxXForCascadeEvt(0), fHistPosTrkgPrimaryVtxYForCascadeEvt(0), fHistPosTrkgPrimaryVtxZForCascadeEvt(0), fHistTrkgPrimaryVtxRadius(0), |
105 | fHistPosBestPrimaryVtxXForCascadeEvt(0), fHistPosBestPrimaryVtxYForCascadeEvt(0), fHistPosBestPrimaryVtxZForCascadeEvt(0), fHistBestPrimaryVtxRadius(0), | |
45ee0bcc | 106 | f2dHistTrkgPrimVtxVsBestPrimVtx(0), |
107 | ||
108 | fHistEffMassXi(0), fHistChi2Xi(0), | |
109 | fHistDcaXiDaughters(0), fHistDcaBachToPrimVertex(0), fHistXiCosineOfPointingAngle(0), fHistXiRadius(0), | |
110 | ||
111 | fHistMassLambdaAsCascDghter(0), | |
112 | fHistV0Chi2Xi(0), | |
113 | fHistDcaV0DaughtersXi(0), | |
114 | fHistDcaV0ToPrimVertexXi(0), | |
115 | fHistV0CosineOfPointingAngleXi(0), | |
116 | fHistV0RadiusXi(0), | |
117 | fHistDcaPosToPrimVertexXi(0), fHistDcaNegToPrimVertexXi(0), | |
118 | ||
119 | fHistMassXiMinus(0), fHistMassXiPlus(0), | |
120 | fHistMassOmegaMinus(0), fHistMassOmegaPlus(0), | |
f87cd3db | 121 | fHistMassWithCombPIDXiMinus(0), fHistMassWithCombPIDXiPlus(0), |
122 | fHistMassWithCombPIDOmegaMinus(0), fHistMassWithCombPIDOmegaPlus(0), | |
45ee0bcc | 123 | |
124 | fHistXiTransvMom(0), fHistXiTotMom(0), | |
90926f02 | 125 | fHistBachTransvMomXi(0), fHistBachTotMomXi(0), |
45ee0bcc | 126 | |
127 | fHistChargeXi(0), | |
128 | fHistV0toXiCosineOfPointingAngle(0), | |
129 | ||
90926f02 | 130 | fHistRapXi(0), fHistRapOmega(0), fHistEtaXi(0), |
131 | fHistThetaXi(0), fHistPhiXi(0), | |
45ee0bcc | 132 | |
c525f360 | 133 | fHistcTauXiMinus(0), fHistcTauXiPlus(0), fHistcTauOmegaMinus(0), fHistcTauOmegaPlus(0), |
134 | ||
45ee0bcc | 135 | f2dHistArmenteros(0), |
136 | f2dHistEffMassLambdaVsEffMassXiMinus(0), f2dHistEffMassXiVsEffMassOmegaMinus(0), | |
137 | f2dHistEffMassLambdaVsEffMassXiPlus(0), f2dHistEffMassXiVsEffMassOmegaPlus(0), | |
e9c3f0b6 | 138 | f2dHistXiRadiusVsEffMassXiMinus(0), f2dHistXiRadiusVsEffMassXiPlus(0), |
139 | f2dHistXiRadiusVsEffMassOmegaMinus(0), f2dHistXiRadiusVsEffMassOmegaPlus(0), | |
140 | ||
e8fd19de | 141 | f2dHistTPCdEdxOfCascDghters(0), |
142 | ||
ff0753e1 | 143 | f3dHistXiPtVsEffMassVsYXiMinus(0), f3dHistXiPtVsEffMassVsYXiPlus(0), |
3873f6ee | 144 | f3dHistXiPtVsEffMassVsYOmegaMinus(0), f3dHistXiPtVsEffMassVsYOmegaPlus(0), |
145 | ||
32599cd4 | 146 | fCFContCascadePIDXiMinus(0), |
147 | fCFContCascadePIDXiPlus(0), | |
148 | fCFContCascadePIDOmegaMinus(0), | |
149 | fCFContCascadePIDOmegaPlus(0), | |
150 | fCFContCascadeCuts(0), | |
3aedd4a5 | 151 | |
3873f6ee | 152 | fHnSpAngularCorrXiMinus(0), fHnSpAngularCorrXiPlus(0), |
153 | fHnSpAngularCorrOmegaMinus(0), fHnSpAngularCorrOmegaPlus(0) | |
45ee0bcc | 154 | |
155 | { | |
156 | // Dummy Constructor | |
9d670198 | 157 | for(Int_t iAlephIdx = 0; iAlephIdx < 5; iAlephIdx++ ) { fAlephParameters [iAlephIdx] = -1.; } |
158 | for(Int_t iV0selIdx = 0; iV0selIdx < 7; iV0selIdx++ ) { fV0Sels [iV0selIdx ] = -1.; } | |
159 | for(Int_t iCascSelIdx = 0; iCascSelIdx < 8; iCascSelIdx++ ) { fCascSels [iCascSelIdx ] = -1.; } | |
45ee0bcc | 160 | } |
161 | ||
162 | ||
163 | ||
164 | ||
165 | ||
166 | ||
167 | ||
168 | ||
169 | //________________________________________________________________________ | |
170 | AliAnalysisTaskCheckCascade::AliAnalysisTaskCheckCascade(const char *name) | |
c525f360 | 171 | : AliAnalysisTaskSE(name), fAnalysisType("ESD"), fTriggerMaskType("kMB"), fCollidingSystems(0), fESDpid(0), fESDtrackCuts(0), /*fPaveTextBookKeeping(0),*/ |
172 | fkRerunV0CascVertexers (0), | |
173 | fkQualityCutZprimVtxPos (kTRUE), | |
174 | fkRejectEventPileUp (kTRUE), | |
175 | fkQualityCutNoTPConlyPrimVtx (kTRUE), | |
176 | fkQualityCutTPCrefit (kTRUE), | |
177 | fkQualityCut80TPCcls (kTRUE), | |
178 | fkIsDataRecoWith1PadTPCCluster (kTRUE), | |
179 | fkExtraSelections (0), | |
180 | fAngularCorrelationType ("TrigLeadingTrck-AssoCasc"), | |
7e504402 | 181 | |
45ee0bcc | 182 | // - Cascade part initialisation |
183 | fListHistCascade(0), | |
9d670198 | 184 | fHistCascadeMultiplicityBeforeTrigSel(0), |
c525f360 | 185 | fHistCascadeMultiplicityForTrigEvt(0), fHistTrackMultiplicityForTrigEvt(0), fHistTPCrefitTrackMultiplicityForTrigEvt(0), fHistPrimaryTrackMultiplicityForTrigEvt(0), |
a786bd57 | 186 | fHistEstimateITSTPCMultiplicityForTrigEvt(0), |
c525f360 | 187 | fHistCascadeMultiplicityForTrigEvtAndZprimVtx(0), fHistCascadeMultiplicityForTrigEvtNonPiledUpAndZprimVtx(0), |
9d670198 | 188 | fHistCascadeMultiplicityForSelEvt(0), |
189 | fHistPosBestPrimaryVtxXForSelEvt(0), fHistPosBestPrimaryVtxYForSelEvt(0), fHistPosBestPrimaryVtxZForSelEvt(0), | |
c525f360 | 190 | fHistTPCrefitTrackMultiplicityForCascadeEvt(0), fHistPrimaryTrackMultiplicityForCascadeEvt(0), |
7e504402 | 191 | fHistPosV0TPCClusters(0), fHistNegV0TPCClusters(0), fHistBachTPCClusters(0), |
45ee0bcc | 192 | fHistVtxStatus(0), |
193 | ||
9d670198 | 194 | fHistPosTrkgPrimaryVtxXForCascadeEvt(0), fHistPosTrkgPrimaryVtxYForCascadeEvt(0), fHistPosTrkgPrimaryVtxZForCascadeEvt(0), fHistTrkgPrimaryVtxRadius(0), |
195 | fHistPosBestPrimaryVtxXForCascadeEvt(0), fHistPosBestPrimaryVtxYForCascadeEvt(0), fHistPosBestPrimaryVtxZForCascadeEvt(0), fHistBestPrimaryVtxRadius(0), | |
45ee0bcc | 196 | f2dHistTrkgPrimVtxVsBestPrimVtx(0), |
197 | ||
198 | fHistEffMassXi(0), fHistChi2Xi(0), | |
199 | fHistDcaXiDaughters(0), fHistDcaBachToPrimVertex(0), fHistXiCosineOfPointingAngle(0), fHistXiRadius(0), | |
200 | ||
201 | fHistMassLambdaAsCascDghter(0), | |
202 | fHistV0Chi2Xi(0), | |
203 | fHistDcaV0DaughtersXi(0), | |
204 | fHistDcaV0ToPrimVertexXi(0), | |
205 | fHistV0CosineOfPointingAngleXi(0), | |
206 | fHistV0RadiusXi(0), | |
207 | fHistDcaPosToPrimVertexXi(0), fHistDcaNegToPrimVertexXi(0), | |
208 | ||
209 | fHistMassXiMinus(0), fHistMassXiPlus(0), | |
210 | fHistMassOmegaMinus(0), fHistMassOmegaPlus(0), | |
f87cd3db | 211 | fHistMassWithCombPIDXiMinus(0), fHistMassWithCombPIDXiPlus(0), |
212 | fHistMassWithCombPIDOmegaMinus(0), fHistMassWithCombPIDOmegaPlus(0), | |
45ee0bcc | 213 | |
214 | fHistXiTransvMom(0), fHistXiTotMom(0), | |
90926f02 | 215 | fHistBachTransvMomXi(0), fHistBachTotMomXi(0), |
45ee0bcc | 216 | |
217 | fHistChargeXi(0), | |
218 | fHistV0toXiCosineOfPointingAngle(0), | |
219 | ||
90926f02 | 220 | fHistRapXi(0), fHistRapOmega(0), fHistEtaXi(0), |
221 | fHistThetaXi(0), fHistPhiXi(0), | |
45ee0bcc | 222 | |
c525f360 | 223 | fHistcTauXiMinus(0), fHistcTauXiPlus(0), fHistcTauOmegaMinus(0), fHistcTauOmegaPlus(0), |
224 | ||
45ee0bcc | 225 | f2dHistArmenteros(0), |
226 | f2dHistEffMassLambdaVsEffMassXiMinus(0), f2dHistEffMassXiVsEffMassOmegaMinus(0), | |
227 | f2dHistEffMassLambdaVsEffMassXiPlus(0), f2dHistEffMassXiVsEffMassOmegaPlus(0), | |
e9c3f0b6 | 228 | f2dHistXiRadiusVsEffMassXiMinus(0), f2dHistXiRadiusVsEffMassXiPlus(0), |
229 | f2dHistXiRadiusVsEffMassOmegaMinus(0), f2dHistXiRadiusVsEffMassOmegaPlus(0), | |
230 | ||
e8fd19de | 231 | f2dHistTPCdEdxOfCascDghters(0), |
232 | ||
ff0753e1 | 233 | f3dHistXiPtVsEffMassVsYXiMinus(0), f3dHistXiPtVsEffMassVsYXiPlus(0), |
3873f6ee | 234 | f3dHistXiPtVsEffMassVsYOmegaMinus(0), f3dHistXiPtVsEffMassVsYOmegaPlus(0), |
235 | ||
32599cd4 | 236 | fCFContCascadePIDXiMinus(0), |
237 | fCFContCascadePIDXiPlus(0), | |
238 | fCFContCascadePIDOmegaMinus(0), | |
239 | fCFContCascadePIDOmegaPlus(0), | |
240 | fCFContCascadeCuts(0), | |
3aedd4a5 | 241 | |
3873f6ee | 242 | fHnSpAngularCorrXiMinus(0), fHnSpAngularCorrXiPlus(0), |
243 | fHnSpAngularCorrOmegaMinus(0), fHnSpAngularCorrOmegaPlus(0) | |
45ee0bcc | 244 | |
245 | { | |
246 | // Constructor | |
247 | ||
248 | // Define input and output slots here | |
249 | // Input slot #0 works with a TChain | |
9d670198 | 250 | // DefineInput(0, TChain::Class()); |
251 | // Output slot #1 writes into a TList container (cascade) | |
252 | ||
253 | for(Int_t iAlephIdx = 0; iAlephIdx < 5; iAlephIdx++ ) { fAlephParameters [iAlephIdx] = -1.; } | |
254 | ||
a786bd57 | 255 | // New Loose : 1st step for the 7 TeV pp analysis |
e8fd19de | 256 | |
a786bd57 | 257 | fV0Sels[0] = 33. ; // max allowed chi2 |
258 | fV0Sels[1] = 0.02; // min allowed impact parameter for the 1st daughter (LHC09a4 : 0.05) | |
259 | fV0Sels[2] = 0.02; // min allowed impact parameter for the 2nd daughter (LHC09a4 : 0.05) | |
260 | fV0Sels[3] = 2.0 ; // max allowed DCA between the daughter tracks (LHC09a4 : 0.5) | |
261 | fV0Sels[4] = 0.95; // min allowed cosine of V0's pointing angle (LHC09a4 : 0.99) | |
262 | fV0Sels[5] = 1.0 ; // min radius of the fiducial volume (LHC09a4 : 0.2) | |
263 | fV0Sels[6] = 100. ; // max radius of the fiducial volume (LHC09a4 : 100.0) | |
264 | ||
265 | fCascSels[0] = 33. ; // max allowed chi2 (same as PDC07) | |
266 | fCascSels[1] = 0.05 ; // min allowed V0 impact parameter (PDC07 : 0.05 / LHC09a4 : 0.025 ) | |
267 | fCascSels[2] = 0.010; // "window" around the Lambda mass (PDC07 : 0.008 / LHC09a4 : 0.010 ) | |
268 | fCascSels[3] = 0.03 ; // min allowed bachelor's impact parameter (PDC07 : 0.035 / LHC09a4 : 0.025 ) | |
269 | fCascSels[4] = 2.0 ; // max allowed DCA between the V0 and the bachelor (PDC07 : 0.1 / LHC09a4 : 0.2 ) | |
270 | fCascSels[5] = 0.95 ; // min allowed cosine of the cascade pointing angle (PDC07 : 0.9985 / LHC09a4 : 0.998 ) | |
271 | fCascSels[6] = 0.4 ; // min radius of the fiducial volume (PDC07 : 0.9 / LHC09a4 : 0.2 ) | |
272 | fCascSels[7] = 100. ; // max radius of the fiducial volume (PDC07 : 100 / LHC09a4 : 100 ) | |
9d670198 | 273 | |
e8fd19de | 274 | |
a786bd57 | 275 | |
276 | // Hyper Loose "Ã la 900 GeV 2009 data", with lower cosine of pointing angle for Xi (0.95 down to 0.82) = 900 GeV paper | |
277 | /* | |
e8fd19de | 278 | fV0Sels[0] = 33. ; // max allowed chi2 |
279 | fV0Sels[1] = 0.001; // min allowed impact parameter for the 1st daughter (LHC09a4 : 0.05) | |
280 | fV0Sels[2] = 0.001; // min allowed impact parameter for the 2nd daughter (LHC09a4 : 0.05) | |
281 | fV0Sels[3] = 5.0 ; // max allowed DCA between the daughter tracks (LHC09a4 : 0.5) | |
282 | fV0Sels[4] = 0.0 ; // min allowed cosine of V0's pointing angle (LHC09a4 : 0.99) | |
283 | fV0Sels[5] = 0.1 ; // min radius of the fiducial volume (LHC09a4 : 0.2) | |
284 | fV0Sels[6] = 100. ; // max radius of the fiducial volume (LHC09a4 : 100.0) | |
285 | ||
286 | fCascSels[0] = 33. ; // max allowed chi2 (same as PDC07) | |
287 | fCascSels[1] = 0.001; // min allowed V0 impact parameter (PDC07 : 0.05 / LHC09a4 : 0.025 ) | |
288 | fCascSels[2] = 0.008; // "window" around the Lambda mass (PDC07 : 0.008 / LHC09a4 : 0.010 ) | |
289 | fCascSels[3] = 0.001; // min allowed bachelor's impact parameter (PDC07 : 0.035 / LHC09a4 : 0.025 ) | |
290 | fCascSels[4] = 5.0 ; // max allowed DCA between the V0 and the bachelor (PDC07 : 0.1 / LHC09a4 : 0.2 ) | |
291 | fCascSels[5] = 0.82 ; //FIXME min allowed cosine of the cascade pointing angle (PDC07 : 0.9985 / LHC09a4 : 0.998 ) | |
292 | fCascSels[6] = 0.1 ; // min radius of the fiducial volume (PDC07 : 0.9 / LHC09a4 : 0.2 ) | |
293 | fCascSels[7] = 100. ; // max radius of the fiducial volume (PDC07 : 100 / LHC09a4 : 100 ) | |
a786bd57 | 294 | */ |
c525f360 | 295 | |
e8fd19de | 296 | //New default vtxR (http://alisoft.cern.ch/viewvc?view=rev&root=AliRoot&revision=40955, 5 May 2010) |
297 | /* | |
298 | fV0Sels[0] = 33. ; // max allowed chi2 | |
299 | fV0Sels[1] = 0.05; // min allowed impact parameter for the 1st daughter (LHC09a4 : 0.05) | |
300 | fV0Sels[2] = 0.05; // min allowed impact parameter for the 2nd daughter (LHC09a4 : 0.05) | |
301 | fV0Sels[3] = 1.5 ; // max allowed DCA between the daughter tracks (LHC09a4 : 0.5) | |
302 | fV0Sels[4] = 0.9 ; // min allowed cosine of V0's pointing angle (LHC09a4 : 0.99) | |
303 | fV0Sels[5] = 0.2 ; // min radius of the fiducial volume (LHC09a4 : 0.2) | |
304 | fV0Sels[6] = 100. ; // max radius of the fiducial volume (LHC09a4 : 100.0) | |
305 | ||
306 | fCascSels[0] = 33. ; // max allowed chi2 (same as PDC07) | |
307 | fCascSels[1] = 0.01 ; // min allowed V0 impact parameter (PDC07 : 0.05 / LHC09a4 : 0.025 ) | |
308 | fCascSels[2] = 0.008; // "window" around the Lambda mass (PDC07 : 0.008 / LHC09a4 : 0.010 ) | |
309 | fCascSels[3] = 0.01 ; // min allowed bachelor's impact parameter (PDC07 : 0.035 / LHC09a4 : 0.025 ) | |
310 | fCascSels[4] = 2.0 ; // max allowed DCA between the V0 and the bachelor (PDC07 : 0.1 / LHC09a4 : 0.2 ) | |
311 | fCascSels[5] = 0.98 ; // min allowed cosine of the cascade pointing angle (PDC07 : 0.9985 / LHC09a4 : 0.998 ) | |
312 | fCascSels[6] = 0.2 ; // min radius of the fiducial volume (PDC07 : 0.9 / LHC09a4 : 0.2 ) | |
313 | fCascSels[7] = 100. ; // max radius of the fiducial volume (PDC07 : 100 / LHC09a4 : 100 ) | |
314 | */ | |
315 | ||
316 | // Tight for Xi in p-p (http://alisoft.cern.ch/viewvc?view=rev&root=AliRoot&revision=40955, 5 May 2010) | |
317 | /* | |
318 | fV0Sels[0] = 33. ; // max allowed chi2 | |
319 | fV0Sels[1] = 0.05; // min allowed impact parameter for the 1st daughter (LHC09a4 : 0.05) | |
320 | fV0Sels[2] = 0.05; // min allowed impact parameter for the 2nd daughter (LHC09a4 : 0.05) | |
321 | fV0Sels[3] = 0.5 ; // max allowed DCA between the daughter tracks (LHC09a4 : 0.5) | |
322 | fV0Sels[4] = 0.99 ; // min allowed cosine of V0's pointing angle (LHC09a4 : 0.99) | |
323 | fV0Sels[5] = 0.2 ; // min radius of the fiducial volume (LHC09a4 : 0.2) | |
324 | fV0Sels[6] = 100. ; // max radius of the fiducial volume (LHC09a4 : 100.0) | |
325 | ||
326 | fCascSels[0] = 33. ; // max allowed chi2 (same as PDC07) | |
327 | fCascSels[1] = 0.05 ; // min allowed V0 impact parameter (PDC07 : 0.05 / LHC09a4 : 0.025 ) | |
328 | fCascSels[2] = 0.01; // "window" around the Lambda mass (PDC07 : 0.008 / LHC09a4 : 0.010 ) | |
329 | fCascSels[3] = 0.035 ; // min allowed bachelor's impact parameter (PDC07 : 0.035 / LHC09a4 : 0.025 ) | |
330 | fCascSels[4] = 0.10 ; // max allowed DCA between the V0 and the bachelor (PDC07 : 0.1 / LHC09a4 : 0.2 ) | |
331 | fCascSels[5] = 0.9985 ; // min allowed cosine of the cascade pointing angle (PDC07 : 0.9985 / LHC09a4 : 0.998 ) | |
332 | fCascSels[6] = 0.2 ; // min radius of the fiducial volume (PDC07 : 0.9 / LHC09a4 : 0.2 ) | |
333 | fCascSels[7] = 100. ; // max radius of the fiducial volume (PDC07 : 100 / LHC09a4 : 100 ) | |
334 | */ | |
335 | ||
336 | //NOTE | |
337 | // For PbPb coming data, take a look at Iouri's proposal | |
338 | // https://savannah.cern.ch/bugs/index.php?69877 | |
339 | ||
9d670198 | 340 | |
45ee0bcc | 341 | // Output slot #0 writes into a TList container (Cascade) |
342 | DefineOutput(1, TList::Class()); | |
9d670198 | 343 | /*DefineOutput(2, TPaveText::Class());*/ |
45ee0bcc | 344 | } |
345 | ||
346 | ||
9d670198 | 347 | AliAnalysisTaskCheckCascade::~AliAnalysisTaskCheckCascade() |
348 | { | |
349 | // | |
350 | // Destructor | |
351 | // | |
45ee0bcc | 352 | |
9d670198 | 353 | // For all TH1, 2, 3 HnSparse and CFContainer are in the fListCascade TList. |
354 | // They will be deleted when fListCascade is deleted by the TSelector dtor | |
355 | // Because of TList::SetOwner() ... | |
356 | ||
357 | if (fListHistCascade) { delete fListHistCascade; fListHistCascade = 0x0; } | |
358 | if (fESDpid) { delete fESDpid; fESDpid = 0x0;} // fESDpid is not stored into the TList | |
c525f360 | 359 | if (fESDtrackCuts) { delete fESDtrackCuts; fESDtrackCuts = 0x0; } |
9d670198 | 360 | //if (fPaveTextBookKeeping) { delete fPaveTextBookKeeping; fPaveTextBookKeeping = 0x0;} // fPaveTextBookKeeping is not strored in the TList |
361 | } | |
45ee0bcc | 362 | |
363 | ||
364 | ||
365 | //________________________________________________________________________ | |
366 | void AliAnalysisTaskCheckCascade::UserCreateOutputObjects() | |
367 | { | |
368 | // Create histograms | |
369 | // Called once | |
370 | ||
371 | ||
372 | ||
373 | fListHistCascade = new TList(); | |
9d670198 | 374 | fListHistCascade->SetOwner(); // See http://root.cern.ch/root/html/TCollection.html#TCollection:SetOwner |
7e504402 | 375 | |
376 | ||
9d670198 | 377 | if(! fESDpid){ |
378 | ||
379 | AliMCEventHandler *lmcEvtHandler = dynamic_cast<AliMCEventHandler*>( (AliAnalysisManager::GetAnalysisManager())->GetMCtruthEventHandler() ); | |
45ee0bcc | 380 | |
9d670198 | 381 | if( !lmcEvtHandler ){ // !0x0 = real data or !1 = there is an MC handler available (useMC = kTRUE in AnalysisTrainNew), so = data from MC |
382 | ||
c525f360 | 383 | if(fkIsDataRecoWith1PadTPCCluster){ |
384 | // Parameters extracted for LHC10d Pass2 - A. Kalweit | |
385 | fAlephParameters[0] = 1.28949/50.; // = 0,0257898 | |
386 | fAlephParameters[1] = 2.74095e+01; | |
387 | fAlephParameters[2] = TMath::Exp(-3.21763e+01); | |
388 | fAlephParameters[3] = 2.44026; | |
389 | fAlephParameters[4] = 6.58800; | |
390 | } | |
391 | else { | |
392 | // Reasonable parameters extracted for real p-p event (Dec 2009 - GSI Pass5) - A.Kalweit | |
393 | fAlephParameters[0] = 0.0283086; // No extra-division to apply in SetBlochParam | |
394 | fAlephParameters[1] = 2.63394e+01; | |
395 | fAlephParameters[2] = 5.04114e-11; | |
396 | fAlephParameters[3] = 2.12543e+00; | |
397 | fAlephParameters[4] = 4.88663e+00; | |
398 | } | |
399 | ||
400 | Printf("CheckCascade - Check Aleph Param in case of REAL Data (fAlephParameters[3] = %f) (To be compared with : 2.44026 for 1-pad-cluster prod. / 2.12543 otherwise)\n", fAlephParameters[3]); | |
401 | ||
9d670198 | 402 | } |
c525f360 | 403 | else { // MC reco |
404 | if(fkIsDataRecoWith1PadTPCCluster){ | |
405 | // Home made parameterization for LHC10f6a production = p+p 7 TeV | |
406 | fAlephParameters[0] = 0.04; | |
407 | fAlephParameters[1] = 17.5; | |
408 | fAlephParameters[2] = 3.4e-09; | |
409 | fAlephParameters[3] = 2.15; | |
410 | fAlephParameters[4] = 3.91720e+00; | |
411 | ||
412 | // Home made parameterization for LHC10e13 production = p+p 900 GeV/c | |
413 | ||
414 | } | |
415 | else { | |
416 | // Reasonable parameters extracted for p-p simulation (LHC09a4) - A.Kalweit | |
417 | // fAlephParameters[0] = 4.23232575531564326e+00;//50*0.76176e-1; // do not forget to divide this value by 50 in SetBlochParam ! | |
418 | // fAlephParameters[1] = 8.68482806165147636e+00;//10.632; | |
419 | // fAlephParameters[2] = 1.34000000000000005e-05;//0.13279e-4; | |
420 | // fAlephParameters[3] = 2.30445734159456084e+00;//1.8631; | |
421 | // fAlephParameters[4] = 2.25624744086878559e+00;//1.9479; | |
422 | ||
423 | // Reasonable parameters extracted for MC LHC09d10 event (Jan 2010) - A.Kalweit | |
424 | fAlephParameters[0] = 2.15898e+00/50.; | |
425 | fAlephParameters[1] = 1.75295e+01; | |
426 | fAlephParameters[2] = 3.40030e-09; | |
427 | fAlephParameters[3] = 1.96178e+00; | |
428 | fAlephParameters[4] = 3.91720e+00; | |
429 | } | |
430 | Printf("CheckCascade - Check Aleph Param in case of MC Data (fAlephParameters[3] = %f) (To be compared with : 2.15 for 1-pad-cluster prod. / 1.96178 otherwise)\n", fAlephParameters[3]); | |
9d670198 | 431 | } |
432 | ||
433 | fESDpid = new AliESDpid(); | |
434 | fESDpid->GetTPCResponse().SetBetheBlochParameters( fAlephParameters[0], | |
435 | fAlephParameters[1], | |
436 | fAlephParameters[2], | |
437 | fAlephParameters[3], | |
438 | fAlephParameters[4] ); | |
439 | } | |
440 | ||
c525f360 | 441 | if(! fESDtrackCuts ){ |
442 | fESDtrackCuts = AliESDtrackCuts::GetStandardITSTPCTrackCuts2010(kTRUE); // Std definition of primary (see kTRUE argument) tracks for 2010 | |
443 | fESDtrackCuts->SetEtaRange(-0.8,+0.8); | |
444 | fESDtrackCuts->SetPtRange(0.15, 1e10); | |
445 | Printf("CheckCascade - ESDtrackCuts set up to 2010 std ITS-TPC cuts..."); | |
446 | } | |
9d670198 | 447 | |
448 | /* | |
449 | if( !fPaveTextBookKeeping){ | |
e8fd19de | 450 | // FIXME : prepare a field with the AliRoot+root distrib ... |
451 | ||
9d670198 | 452 | fPaveTextBookKeeping = new TPaveText(0.1, 0.1, 0.9, 0.9,"NDC"); |
453 | fPaveTextBookKeeping->SetName("fPaveTextBookKeeping"); | |
454 | fPaveTextBookKeeping->SetBorderSize(0); | |
455 | fPaveTextBookKeeping->SetTextAlign(12); | |
456 | fPaveTextBookKeeping->SetFillColor(kWhite); | |
457 | fPaveTextBookKeeping->SetTextFont(42); // regular Arial or Helvetica, | |
458 | fPaveTextBookKeeping->SetTextColor(kBlue+4); | |
459 | ||
460 | fPaveTextBookKeeping->AddText( "Task CHECK CASCADE analysis" ); | |
461 | fPaveTextBookKeeping->AddText("- - - - - - - - - - - "); | |
462 | fPaveTextBookKeeping->AddText( Form("AnalysisType : %s ", fAnalysisType.Data() )); | |
463 | if(!fCollidingSystems) fPaveTextBookKeeping->AddText("Colliding system : p-p collisions "); | |
464 | else fPaveTextBookKeeping->AddText("Colliding system : A-A collisions "); | |
465 | ||
466 | fPaveTextBookKeeping->AddText("- - - - - - - - - - - "); | |
467 | ||
468 | if(fkRerunV0CascVertexers){ | |
469 | fPaveTextBookKeeping->AddText("A.1. With V0 vertexer : "); | |
470 | fPaveTextBookKeeping->AddText( Form(" - V0 #chi^{2} _________________ < %.3f ", fV0Sels[0] )); | |
471 | fPaveTextBookKeeping->AddText( Form(" - DCA(prim. Vtx/ 1^{st} daughter) ___ > %.3f cm ", fV0Sels[1] )); | |
472 | fPaveTextBookKeeping->AddText( Form(" - DCA(prim. Vtx/ 2^{nd} daughter) __ > %.3f cm", fV0Sels[2] )); | |
473 | fPaveTextBookKeeping->AddText( Form(" - DCA between V0 daughters ___ < %.3f cm", fV0Sels[3] )); | |
474 | fPaveTextBookKeeping->AddText( Form(" - cos(V0 pointing angle) _______ > %.3f ", fV0Sels[4] )); | |
475 | fPaveTextBookKeeping->AddText( Form(" - R_{transv}(V0 decay) ________ > %.3f cm", fV0Sels[5] )); | |
476 | fPaveTextBookKeeping->AddText( Form(" - R_{transv}(V0 decay) ________ < %.3f cm", fV0Sels[6] )); | |
477 | ||
478 | fPaveTextBookKeeping->AddText(" "); | |
479 | ||
480 | fPaveTextBookKeeping->AddText("A.2. With Casc. vertexer : "); | |
481 | fPaveTextBookKeeping->AddText( Form(" - Casc. #chi^{2} ______________ < %.3f ", fCascSels[0] )); | |
482 | fPaveTextBookKeeping->AddText( Form(" - DCA(prim. Vtx/ V0) _________ > %.3f cm", fCascSels[1] )); | |
483 | fPaveTextBookKeeping->AddText( Form(" - | M_{#Lambda}(reco) - M_{#Lambda}(pdg) | _______ < %.3f GeV/c^{2}", fCascSels[2] )); | |
484 | fPaveTextBookKeeping->AddText( Form(" - DCA(prim. Vtx/ Bach) _______ > %.3f cm", fCascSels[3] )); | |
485 | fPaveTextBookKeeping->AddText( Form(" - DCA between Bach/ #Lambda ______ < %.3f cm", fCascSels[4] )); | |
486 | fPaveTextBookKeeping->AddText( Form(" - cos(Casc. pointing angle) ____ > %.3f ", fCascSels[5] )); | |
487 | fPaveTextBookKeeping->AddText( Form(" - R_{transv}(Casc. decay) ______ > %.3f cm", fCascSels[6] )); | |
488 | fPaveTextBookKeeping->AddText( Form(" - R_{transv}(Casc. decay) ______ < %.3f cm", fCascSels[7] )); | |
489 | } | |
490 | else{ fPaveTextBookKeeping->AddText("A. No rerunning of the V0/Casc. vertexers ... See std cuts in (AliRoot+Rec.C) used for this prod. cycle");} | |
491 | ||
492 | fPaveTextBookKeeping->AddText("- - - - - - - - - - - "); | |
493 | ||
494 | if(fkQualityCutZprimVtxPos) fPaveTextBookKeeping->AddText("B. Quality Cut(prim. Vtx z-Pos) = ON "); | |
495 | else fPaveTextBookKeeping->AddText("B. Quality Cut(prim. Vtx z-Pos) = Off "); | |
496 | if(fkQualityCutNoTPConlyPrimVtx) fPaveTextBookKeeping->AddText("C. Quality Cut(No TPC prim. vtx) = ON "); | |
497 | else fPaveTextBookKeeping->AddText("C. Quality Cut(No TPC prim. vtx) = Off "); | |
498 | if(fkQualityCutTPCrefit) fPaveTextBookKeeping->AddText("D. Quality Cut(TPCrefit) = ON "); | |
499 | else fPaveTextBookKeeping->AddText("D. Quality Cut(TPCrefit) = Off "); | |
500 | if(fkQualityCut80TPCcls) fPaveTextBookKeeping->AddText("E. Quality Cut(80 TPC clusters) = ON "); | |
501 | else fPaveTextBookKeeping->AddText("E. Quality Cut(80 TPC clusters) = Off "); | |
502 | if(fkExtraSelections) fPaveTextBookKeeping->AddText("F. Extra Analysis Selections = ON "); | |
503 | else fPaveTextBookKeeping->AddText("F. Extra Analysis Selections = Off "); | |
504 | ||
505 | fPaveTextBookKeeping->AddText("- - - - - - - - - - - "); | |
506 | ||
507 | fPaveTextBookKeeping->AddText("G. TPC Aleph Param : "); | |
508 | fPaveTextBookKeeping->AddText( Form(" - fAlephParam [0] = %.5g", fAlephParameters[0] )); | |
509 | fPaveTextBookKeeping->AddText( Form(" - fAlephParam [1] = %.5g", fAlephParameters[1] )); | |
510 | fPaveTextBookKeeping->AddText( Form(" - fAlephParam [2] = %.5g", fAlephParameters[2] )); | |
511 | fPaveTextBookKeeping->AddText( Form(" - fAlephParam [3] = %.5g", fAlephParameters[3] )); | |
512 | fPaveTextBookKeeping->AddText( Form(" - fAlephParam [4] = %.5g", fAlephParameters[4] )); | |
513 | } | |
514 | */ | |
515 | ||
516 | // - General histos | |
517 | //-------------- | |
518 | ||
519 | if(! fHistCascadeMultiplicityBeforeTrigSel) { | |
520 | if(fCollidingSystems)// AA collisions | |
521 | fHistCascadeMultiplicityBeforeTrigSel = new TH1F("fHistCascadeMultiplicityBeforeTrigSel", | |
522 | "Cascades per event (before Trig. Sel.);Nbr of Cascades/Evt;Events", | |
523 | 100, 0, 100); | |
524 | else // pp collisions | |
525 | fHistCascadeMultiplicityBeforeTrigSel = new TH1F("fHistCascadeMultiplicityBeforeTrigSel", | |
526 | "Cascades per event (before Trig. Sel.);Nbr of Cascades/Evt;Events", | |
527 | 25, 0, 25); | |
528 | fListHistCascade->Add(fHistCascadeMultiplicityBeforeTrigSel); | |
529 | } | |
530 | ||
531 | // - Histos for events passing the trigger selection | |
532 | //-------------- | |
533 | ||
534 | if(! fHistCascadeMultiplicityForTrigEvt) { | |
535 | if(fCollidingSystems)// AA collisions | |
536 | fHistCascadeMultiplicityForTrigEvt = new TH1F("fHistCascadeMultiplicityForTrigEvt", | |
537 | "Cascades per event (for triggered evt);Nbr of Cascades/Evt;Events", | |
538 | 100, 0, 100); | |
539 | else // pp collisions | |
540 | fHistCascadeMultiplicityForTrigEvt = new TH1F("fHistCascadeMultiplicityForTrigEvt", | |
541 | "Cascades per event (for triggered evt);Nbr of Cascades/Evt;Events", | |
542 | 25, 0, 25); | |
543 | fListHistCascade->Add(fHistCascadeMultiplicityForTrigEvt); | |
544 | } | |
545 | ||
546 | ||
547 | if(! fHistTrackMultiplicityForTrigEvt) { | |
45ee0bcc | 548 | if(fCollidingSystems)// AA collisions |
9d670198 | 549 | fHistTrackMultiplicityForTrigEvt = new TH1F("fHistTrackMultiplicityForTrigEvt", |
550 | "Track Multiplicity (for triggered evt);Nbr of tracks/Evt;Events", | |
32599cd4 | 551 | 200, 0, 20000); |
45ee0bcc | 552 | else // pp collisions |
9d670198 | 553 | fHistTrackMultiplicityForTrigEvt = new TH1F("fHistTrackMultiplicityForTrigEvt", |
554 | "Track Multiplicity (for triggered evt);Nbr of tracks/Evt;Events", | |
e8fd19de | 555 | 300, 0, 300); |
9d670198 | 556 | fListHistCascade->Add(fHistTrackMultiplicityForTrigEvt); |
45ee0bcc | 557 | } |
558 | ||
9d670198 | 559 | if(! fHistTPCrefitTrackMultiplicityForTrigEvt) { |
7e504402 | 560 | if(fCollidingSystems)// AA collisions |
9d670198 | 561 | fHistTPCrefitTrackMultiplicityForTrigEvt = new TH1F("fHistTPCrefitTrackMultiplicityForTrigEvt", |
562 | "TPCrefit track Multiplicity (for triggered evt);Nbr of TPCrefit tracks/Evt;Events", | |
7e504402 | 563 | 200, 0, 20000); |
564 | else // pp collisions | |
9d670198 | 565 | fHistTPCrefitTrackMultiplicityForTrigEvt = new TH1F("fHistTPCrefitTrackMultiplicityForTrigEvt", |
566 | "TPCrefit track Multiplicity (for triggered evt);Nbr of TPCrefit tracks/Evt;Events", | |
e8fd19de | 567 | 300, 0, 300); |
9d670198 | 568 | fListHistCascade->Add(fHistTPCrefitTrackMultiplicityForTrigEvt); |
7e504402 | 569 | } |
570 | ||
9d670198 | 571 | |
c525f360 | 572 | if(! fHistPrimaryTrackMultiplicityForTrigEvt) { |
573 | if(fCollidingSystems)// AA collisions | |
574 | fHistPrimaryTrackMultiplicityForTrigEvt = new TH1F("fHistPrimaryTrackMultiplicityForTrigEvt", | |
575 | "Primary track Multiplicity (for triggered evt);Nbr of primary tracks/Evt;Events", | |
576 | 100, 0, 10000); | |
577 | else // pp collisions | |
578 | fHistPrimaryTrackMultiplicityForTrigEvt = new TH1F("fHistPrimaryTrackMultiplicityForTrigEvt", | |
579 | "Primary track Multiplicity (for triggered evt);Nbr of primary tracks/Evt;Events", | |
580 | 200, 0, 200); | |
581 | fListHistCascade->Add(fHistPrimaryTrackMultiplicityForTrigEvt); | |
582 | } | |
e8fd19de | 583 | |
584 | ||
a786bd57 | 585 | if(! fHistEstimateITSTPCMultiplicityForTrigEvt) { |
586 | if(fCollidingSystems)// AA collisions | |
587 | fHistEstimateITSTPCMultiplicityForTrigEvt = new TH1F("fHistEstimateITSTPCMultiplicityForTrigEvt", | |
588 | "(ITS+TPC tracks + SPD tracklets) multiplicity (for triggered evt);Nbr of (ITS+TPC tracks + SPD tracklets)/Evt;Events", | |
589 | 100, 0, 10000); | |
590 | else // pp collisions | |
591 | fHistEstimateITSTPCMultiplicityForTrigEvt = new TH1F("fHistEstimateITSTPCMultiplicityForTrigEvt", | |
592 | "(ITS+TPC tracks + SPD tracklets) multiplicity (for triggered evt);Nbr of (ITS+TPC tracks + SPD tracklets)/Evt;Events", | |
593 | 200, 0, 200); | |
594 | fListHistCascade->Add(fHistEstimateITSTPCMultiplicityForTrigEvt); | |
595 | } | |
596 | ||
597 | ||
598 | ||
e8fd19de | 599 | // - Histos for events passing the trigger selection + |z(prim. vertex)| < XX cm |
600 | //-------------- | |
601 | ||
602 | if(! fHistCascadeMultiplicityForTrigEvtAndZprimVtx) { | |
603 | if(fCollidingSystems)// AA collisions | |
604 | fHistCascadeMultiplicityForTrigEvtAndZprimVtx = new TH1F("fHistCascadeMultiplicityForTrigEvtAndZprimVtx", | |
605 | "Cascades per event;Nbr of Cascades/Evt;Events", | |
606 | 100, 0, 100); | |
607 | else // pp collisions | |
608 | fHistCascadeMultiplicityForTrigEvtAndZprimVtx = new TH1F("fHistCascadeMultiplicityForTrigEvtAndZprimVtx", | |
609 | "Cascades per event;Nbr of Cascades/Evt;Events", | |
610 | 25, 0, 25); | |
611 | fListHistCascade->Add(fHistCascadeMultiplicityForTrigEvtAndZprimVtx); | |
612 | } | |
613 | ||
614 | ||
c525f360 | 615 | // - Histos for events passing the trigger selection + |z(prim. vertex)| < XX cm + after pile-up rejection |
616 | //-------------- | |
617 | ||
618 | if(! fHistCascadeMultiplicityForTrigEvtNonPiledUpAndZprimVtx) { | |
619 | if(fCollidingSystems)// AA collisions | |
620 | fHistCascadeMultiplicityForTrigEvtNonPiledUpAndZprimVtx = new TH1F("fHistCascadeMultiplicityForTrigEvtNonPiledUpAndZprimVtx", | |
621 | "Cascades per event;Nbr of Cascades/Evt;Events", | |
622 | 100, 0, 100); | |
623 | else // pp collisions | |
624 | fHistCascadeMultiplicityForTrigEvtNonPiledUpAndZprimVtx = new TH1F("fHistCascadeMultiplicityForTrigEvtNonPiledUpAndZprimVtx", | |
625 | "Cascades per event;Nbr of Cascades/Evt;Events", | |
626 | 25, 0, 25); | |
627 | fListHistCascade->Add(fHistCascadeMultiplicityForTrigEvtNonPiledUpAndZprimVtx); | |
628 | } | |
629 | ||
630 | ||
9d670198 | 631 | // - Histos for events passing the event selection at the analysis level |
632 | //-------------- | |
633 | ||
634 | if(! fHistCascadeMultiplicityForSelEvt) { | |
45ee0bcc | 635 | if(fCollidingSystems)// AA collisions |
9d670198 | 636 | fHistCascadeMultiplicityForSelEvt = new TH1F("fHistCascadeMultiplicityForSelEvt", |
7e504402 | 637 | "Cascades per event;Nbr of Cascades/Evt;Events", |
f3dc9369 | 638 | 100, 0, 100); |
45ee0bcc | 639 | else // pp collisions |
9d670198 | 640 | fHistCascadeMultiplicityForSelEvt = new TH1F("fHistCascadeMultiplicityForSelEvt", |
45ee0bcc | 641 | "Cascades per event;Nbr of Cascades/Evt;Events", |
f3dc9369 | 642 | 25, 0, 25); |
9d670198 | 643 | fListHistCascade->Add(fHistCascadeMultiplicityForSelEvt); |
644 | } | |
645 | ||
c525f360 | 646 | |
9d670198 | 647 | if(! fHistPosBestPrimaryVtxXForSelEvt ){ |
e8fd19de | 648 | fHistPosBestPrimaryVtxXForSelEvt = new TH1F( "fHistPosBestPrimaryVtxXForSelEvt" , "Best Prim. Vertex Position in x; x (cm); Events" , 360, -0.9, 0.9 ); |
9d670198 | 649 | fListHistCascade->Add(fHistPosBestPrimaryVtxXForSelEvt); |
650 | } | |
651 | ||
652 | if(! fHistPosBestPrimaryVtxYForSelEvt){ | |
e8fd19de | 653 | fHistPosBestPrimaryVtxYForSelEvt = new TH1F( "fHistPosBestPrimaryVtxYForSelEvt" , "Best Prim. Vertex Position in y; y (cm); Events" , 360, -0.9, 0.9 ); |
9d670198 | 654 | fListHistCascade->Add(fHistPosBestPrimaryVtxYForSelEvt); |
45ee0bcc | 655 | } |
656 | ||
9d670198 | 657 | if(! fHistPosBestPrimaryVtxZForSelEvt ){ |
e8fd19de | 658 | fHistPosBestPrimaryVtxZForSelEvt = new TH1F( "fHistPosBestPrimaryVtxZForSelEvt" , "Best Prim. Vertex Position in z; z (cm); Events" , 300, -30.0, 30.0 ); |
9d670198 | 659 | fListHistCascade->Add(fHistPosBestPrimaryVtxZForSelEvt); |
660 | } | |
661 | ||
662 | ||
45ee0bcc | 663 | |
664 | ||
9d670198 | 665 | // - Histos for events containing at least ONE CASCADE |
666 | //-------------- | |
667 | ||
7e504402 | 668 | if(! fHistTPCrefitTrackMultiplicityForCascadeEvt) { |
669 | if(fCollidingSystems)// AA collisions | |
670 | fHistTPCrefitTrackMultiplicityForCascadeEvt = new TH1F("fHistTPCrefitTrackMultiplicityForCascadeEvt", | |
671 | "TPCrefit track Multiplicity (for evt with Casc.);Nbr of TPCrefit tracks/Evt with cascade(s);Events", | |
672 | 200, 0, 20000); | |
673 | else // pp collisions | |
674 | fHistTPCrefitTrackMultiplicityForCascadeEvt = new TH1F("fHistTPCrefitTrackMultiplicityForCascadeEvt", | |
675 | "TPCrefit track Multiplicity (for evt with Casc.);Nbr of TPCrefit tracks/Evt with cascade(s);Events", | |
e8fd19de | 676 | 300, 0, 300); |
7e504402 | 677 | fListHistCascade->Add(fHistTPCrefitTrackMultiplicityForCascadeEvt); |
45ee0bcc | 678 | } |
679 | ||
c525f360 | 680 | if(! fHistPrimaryTrackMultiplicityForCascadeEvt) { |
681 | if(fCollidingSystems)// AA collisions | |
682 | fHistPrimaryTrackMultiplicityForCascadeEvt = new TH1F("fHistPrimaryTrackMultiplicityForCascadeEvt", | |
683 | "Primary track Multiplicity (for evt with Casc.);Nbr of primary tracks/Evt;Events", | |
684 | 100, 0, 10000); | |
685 | else // pp collisions | |
686 | fHistPrimaryTrackMultiplicityForCascadeEvt = new TH1F("fHistPrimaryTrackMultiplicityForCascadeEvt", | |
687 | "Primary track Multiplicity (for evt with Casc.);Nbr of primary tracks/Evt;Events", | |
688 | 200, 0, 200); | |
689 | fListHistCascade->Add(fHistPrimaryTrackMultiplicityForCascadeEvt); | |
690 | } | |
691 | ||
7e504402 | 692 | if(! fHistPosV0TPCClusters ){ |
693 | 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); | |
694 | fListHistCascade->Add(fHistPosV0TPCClusters); | |
695 | } | |
696 | ||
697 | if(! fHistNegV0TPCClusters ){ | |
698 | 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); | |
699 | fListHistCascade->Add(fHistNegV0TPCClusters); | |
700 | } | |
701 | ||
702 | if(! fHistBachTPCClusters ){ | |
703 | fHistBachTPCClusters = new TH1F("fHistBachTPCClusters", "TPC clusters for Bachelor track; Nbr of TPC clusters (Bach); Track counts", 165, 0.0 ,165.0); | |
704 | fListHistCascade->Add(fHistBachTPCClusters); | |
705 | } | |
45ee0bcc | 706 | |
707 | ||
708 | ||
709 | ||
710 | ||
7e504402 | 711 | |
712 | ||
713 | if(! fHistVtxStatus ){ | |
714 | fHistVtxStatus = new TH1F( "fHistVtxStatus" , "Does a Trckg Prim.vtx exist ?; true=1 or false=0; Nb of Events" , 4, -1.0, 3.0 ); | |
715 | fListHistCascade->Add(fHistVtxStatus); | |
716 | } | |
717 | ||
718 | ||
45ee0bcc | 719 | // - Vertex Positions |
720 | ||
9d670198 | 721 | if(! fHistPosTrkgPrimaryVtxXForCascadeEvt ){ |
e8fd19de | 722 | fHistPosTrkgPrimaryVtxXForCascadeEvt = new TH1F( "fHistPosTrkgPrimaryVtxXForCascadeEvt" , "Trkg Prim. Vertex Position in x; x (cm); Events" , 360, -0.9, 0.9 ); |
9d670198 | 723 | fListHistCascade->Add(fHistPosTrkgPrimaryVtxXForCascadeEvt); |
45ee0bcc | 724 | } |
725 | ||
726 | ||
9d670198 | 727 | if(! fHistPosTrkgPrimaryVtxYForCascadeEvt){ |
e8fd19de | 728 | fHistPosTrkgPrimaryVtxYForCascadeEvt = new TH1F( "fHistPosTrkgPrimaryVtxYForCascadeEvt" , "Trkg Prim. Vertex Position in y; y (cm); Events" , 360, -0.9, 0.9 ); |
9d670198 | 729 | fListHistCascade->Add(fHistPosTrkgPrimaryVtxYForCascadeEvt); |
45ee0bcc | 730 | } |
731 | ||
9d670198 | 732 | if(! fHistPosTrkgPrimaryVtxZForCascadeEvt ){ |
e8fd19de | 733 | fHistPosTrkgPrimaryVtxZForCascadeEvt = new TH1F( "fHistPosTrkgPrimaryVtxZForCascadeEvt" , "Trkg Prim. Vertex Position in z; z (cm); Events" , 300, -30.0, 30.0 ); |
9d670198 | 734 | fListHistCascade->Add(fHistPosTrkgPrimaryVtxZForCascadeEvt); |
45ee0bcc | 735 | } |
736 | ||
737 | if(! fHistTrkgPrimaryVtxRadius ){ | |
738 | fHistTrkgPrimaryVtxRadius = new TH1F( "fHistTrkgPrimaryVtxRadius", "Trkg Prim. Vertex radius; r (cm); Events" , 150, 0., 15.0 ); | |
739 | fListHistCascade->Add(fHistTrkgPrimaryVtxRadius); | |
740 | } | |
741 | ||
742 | ||
743 | ||
744 | ||
9d670198 | 745 | if(! fHistPosBestPrimaryVtxXForCascadeEvt ){ |
e8fd19de | 746 | fHistPosBestPrimaryVtxXForCascadeEvt = new TH1F( "fHistPosBestPrimaryVtxXForCascadeEvt" , "Best Prim. Vertex Position in x; x (cm); Events" , 360, -0.9, 0.9 ); |
9d670198 | 747 | fListHistCascade->Add(fHistPosBestPrimaryVtxXForCascadeEvt); |
45ee0bcc | 748 | } |
749 | ||
9d670198 | 750 | if(! fHistPosBestPrimaryVtxYForCascadeEvt){ |
e8fd19de | 751 | fHistPosBestPrimaryVtxYForCascadeEvt = new TH1F( "fHistPosBestPrimaryVtxYForCascadeEvt" , "Best Prim. Vertex Position in y; y (cm); Events" , 360, -0.9, 0.9 ); |
9d670198 | 752 | fListHistCascade->Add(fHistPosBestPrimaryVtxYForCascadeEvt); |
45ee0bcc | 753 | } |
754 | ||
9d670198 | 755 | if(! fHistPosBestPrimaryVtxZForCascadeEvt ){ |
e8fd19de | 756 | fHistPosBestPrimaryVtxZForCascadeEvt = new TH1F( "fHistPosBestPrimaryVtxZForCascadeEvt" , "Best Prim. Vertex Position in z; z (cm); Events" , 300, -30.0, 30.0 ); |
9d670198 | 757 | fListHistCascade->Add(fHistPosBestPrimaryVtxZForCascadeEvt); |
45ee0bcc | 758 | } |
759 | ||
760 | if(! fHistBestPrimaryVtxRadius ){ | |
761 | fHistBestPrimaryVtxRadius = new TH1F( "fHistBestPrimaryVtxRadius", "Best Prim. vertex radius; r (cm); Events" , 150, 0., 15.0 ); | |
762 | fListHistCascade->Add(fHistBestPrimaryVtxRadius); | |
763 | } | |
764 | ||
765 | if(! f2dHistTrkgPrimVtxVsBestPrimVtx) { | |
766 | f2dHistTrkgPrimVtxVsBestPrimVtx = new TH2F( "f2dHistTrkgPrimVtxVsBestPrimVtx", "r_{Trck Prim. Vtx} Vs r_{Best Prim. Vtx}; r_{Track Vtx} (cm); r_{Best Vtx} (cm)", 300, 0., 15.0, 300, 0., 15.); | |
767 | fListHistCascade->Add(f2dHistTrkgPrimVtxVsBestPrimVtx); | |
768 | } | |
769 | ||
770 | ||
771 | ||
772 | ||
773 | // - Typical histos for cascades | |
774 | ||
775 | ||
776 | if(! fHistEffMassXi) { | |
90926f02 | 777 | fHistEffMassXi = new TH1F("fHistEffMassXi", "Cascade candidates ; Invariant Mass (GeV/c^{2}) ; Counts", 400, 1.2, 2.0); |
45ee0bcc | 778 | fListHistCascade->Add(fHistEffMassXi); |
779 | } | |
780 | ||
781 | if(! fHistChi2Xi ){ | |
782 | fHistChi2Xi = new TH1F("fHistChi2Xi", "Cascade #chi^{2}; #chi^{2}; Number of Cascades", 160, 0, 40); | |
783 | fListHistCascade->Add(fHistChi2Xi); | |
784 | } | |
785 | ||
786 | if(! fHistDcaXiDaughters ){ | |
787 | fHistDcaXiDaughters = new TH1F( "fHistDcaXiDaughters", "DCA between Xi Daughters; DCA (cm) ; Number of Cascades", 100, 0., 0.5); | |
788 | fListHistCascade->Add(fHistDcaXiDaughters); | |
789 | } | |
790 | ||
791 | if(! fHistDcaBachToPrimVertex) { | |
792 | fHistDcaBachToPrimVertex = new TH1F("fHistDcaBachToPrimVertex", "DCA of Bach. to Prim. Vertex;DCA (cm);Number of Cascades", 250, 0., 0.25); | |
793 | fListHistCascade->Add(fHistDcaBachToPrimVertex); | |
794 | } | |
795 | ||
796 | if(! fHistXiCosineOfPointingAngle) { | |
32599cd4 | 797 | fHistXiCosineOfPointingAngle = new TH1F("fHistXiCosineOfPointingAngle", "Cosine of Xi Pointing Angle; Cos (Xi Point.Angl);Number of Xis", 200, 0.99, 1.0); |
45ee0bcc | 798 | fListHistCascade->Add(fHistXiCosineOfPointingAngle); |
799 | } | |
800 | ||
801 | if(! fHistXiRadius ){ | |
32599cd4 | 802 | fHistXiRadius = new TH1F( "fHistXiRadius", "Casc. decay transv. radius; r (cm); Counts" , 1050, 0., 105.0 ); |
45ee0bcc | 803 | fListHistCascade->Add(fHistXiRadius); |
804 | } | |
805 | ||
806 | ||
807 | // - Histos about ~ the "V0 part" of the cascade, coming by inheritance from AliESDv0 | |
808 | ||
809 | ||
810 | ||
811 | if (! fHistMassLambdaAsCascDghter) { | |
3aedd4a5 | 812 | fHistMassLambdaAsCascDghter = new TH1F("fHistMassLambdaAsCascDghter","#Lambda associated to Casc. candidates;Eff. Mass (GeV/c^{2});Counts", 300,1.00,1.3); |
45ee0bcc | 813 | fListHistCascade->Add(fHistMassLambdaAsCascDghter); |
814 | } | |
815 | ||
816 | if (! fHistV0Chi2Xi) { | |
817 | fHistV0Chi2Xi = new TH1F("fHistV0Chi2Xi", "V0 #chi^{2}, in cascade; #chi^{2};Counts", 160, 0, 40); | |
818 | fListHistCascade->Add(fHistV0Chi2Xi); | |
819 | } | |
820 | ||
821 | if (! fHistDcaV0DaughtersXi) { | |
822 | fHistDcaV0DaughtersXi = new TH1F("fHistDcaV0DaughtersXi", "DCA between V0 daughters, in cascade;DCA (cm);Number of V0s", 120, 0., 0.6); | |
823 | fListHistCascade->Add(fHistDcaV0DaughtersXi); | |
824 | } | |
825 | ||
826 | if (! fHistDcaV0ToPrimVertexXi) { | |
827 | fHistDcaV0ToPrimVertexXi = new TH1F("fHistDcaV0ToPrimVertexXi", "DCA of V0 to Prim. Vertex, in cascade;DCA (cm);Number of Cascades", 200, 0., 1.); | |
828 | fListHistCascade->Add(fHistDcaV0ToPrimVertexXi); | |
829 | } | |
830 | ||
831 | if (! fHistV0CosineOfPointingAngleXi) { | |
832 | fHistV0CosineOfPointingAngleXi = new TH1F("fHistV0CosineOfPointingAngleXi", "Cosine of V0 Pointing Angle, in cascade;Cos(V0 Point. Angl); Counts", 200, 0.98, 1.0); | |
833 | fListHistCascade->Add(fHistV0CosineOfPointingAngleXi); | |
834 | } | |
835 | ||
836 | if (! fHistV0RadiusXi) { | |
32599cd4 | 837 | fHistV0RadiusXi = new TH1F("fHistV0RadiusXi", "V0 decay radius, in cascade; radius (cm); Counts", 1050, 0., 105.0); |
45ee0bcc | 838 | fListHistCascade->Add(fHistV0RadiusXi); |
839 | } | |
840 | ||
841 | if (! fHistDcaPosToPrimVertexXi) { | |
842 | fHistDcaPosToPrimVertexXi = new TH1F("fHistDcaPosToPrimVertexXi", "DCA of V0 pos daughter to Prim. Vertex;DCA (cm);Counts", 300, 0, 3); | |
843 | fListHistCascade->Add(fHistDcaPosToPrimVertexXi); | |
844 | } | |
845 | ||
846 | if (! fHistDcaNegToPrimVertexXi) { | |
847 | fHistDcaNegToPrimVertexXi = new TH1F("fHistDcaNegToPrimVertexXi", "DCA of V0 neg daughter to Prim. Vertex;DCA (cm);Counts", 300, 0, 3); | |
848 | fListHistCascade->Add(fHistDcaNegToPrimVertexXi); | |
849 | } | |
850 | ||
851 | ||
852 | ||
853 | ||
854 | // - Effective mass histos for cascades. | |
f87cd3db | 855 | // By cascade hyp |
45ee0bcc | 856 | if (! fHistMassXiMinus) { |
90926f02 | 857 | fHistMassXiMinus = new TH1F("fHistMassXiMinus","#Xi^{-} candidates;M( #Lambda , #pi^{-} ) (GeV/c^{2});Counts", 400,1.2,2.0); |
45ee0bcc | 858 | fListHistCascade->Add(fHistMassXiMinus); |
859 | } | |
860 | ||
861 | if (! fHistMassXiPlus) { | |
a786bd57 | 862 | fHistMassXiPlus = new TH1F("fHistMassXiPlus","#bar{#Xi}^{+} candidates;M( #bar{#Lambda}^{0} , #pi^{+} ) (GeV/c^{2});Counts",400,1.2,2.0); |
45ee0bcc | 863 | fListHistCascade->Add(fHistMassXiPlus); |
864 | } | |
865 | ||
866 | if (! fHistMassOmegaMinus) { | |
90926f02 | 867 | fHistMassOmegaMinus = new TH1F("fHistMassOmegaMinus","#Omega^{-} candidates;M( #Lambda , K^{-} ) (GeV/c^{2});Counts", 500,1.5,2.5); |
45ee0bcc | 868 | fListHistCascade->Add(fHistMassOmegaMinus); |
869 | } | |
870 | ||
871 | if (! fHistMassOmegaPlus) { | |
a786bd57 | 872 | fHistMassOmegaPlus = new TH1F("fHistMassOmegaPlus","#bar{#Omega}^{+} candidates;M( #bar{#Lambda}^{0} , K^{+} ) (GeV/c^{2});Counts", 500,1.5,2.5); |
45ee0bcc | 873 | fListHistCascade->Add(fHistMassOmegaPlus); |
874 | } | |
875 | ||
f87cd3db | 876 | // By cascade hyp + bachelor PID |
877 | if (! fHistMassWithCombPIDXiMinus) { | |
90926f02 | 878 | fHistMassWithCombPIDXiMinus = new TH1F("fHistMassWithCombPIDXiMinus","#Xi^{-} candidates, with Bach. comb. PID;M( #Lambda , #pi^{-} ) (GeV/c^{2});Counts", 400,1.2,2.0); |
f87cd3db | 879 | fListHistCascade->Add(fHistMassWithCombPIDXiMinus); |
880 | } | |
881 | ||
882 | if (! fHistMassWithCombPIDXiPlus) { | |
a786bd57 | 883 | fHistMassWithCombPIDXiPlus = new TH1F("fHistMassWithCombPIDXiPlus","#bar{#Xi}^{+} candidates, with Bach. comb. PID;M( #bar{#Lambda}^{0} , #pi^{+} ) (GeV/c^{2});Counts",400,1.2,2.0); |
f87cd3db | 884 | fListHistCascade->Add(fHistMassWithCombPIDXiPlus); |
885 | } | |
886 | ||
887 | if (! fHistMassWithCombPIDOmegaMinus) { | |
90926f02 | 888 | fHistMassWithCombPIDOmegaMinus = new TH1F("fHistMassWithCombPIDOmegaMinus","#Omega^{-} candidates, with Bach. comb. PID;M( #Lambda , K^{-} ) (GeV/c^{2});Counts", 500,1.5,2.5); |
f87cd3db | 889 | fListHistCascade->Add(fHistMassWithCombPIDOmegaMinus); |
890 | } | |
891 | ||
892 | if (! fHistMassWithCombPIDOmegaPlus) { | |
a786bd57 | 893 | fHistMassWithCombPIDOmegaPlus = new TH1F("fHistMassWithCombPIDOmegaPlus","#bar{#Omega}^{+} candidates, with Bach. comb. PID;M( #bar{#Lambda}^{0} , K^{+} ) (GeV/c^{2});Counts", 500,1.5,2.5); |
f87cd3db | 894 | fListHistCascade->Add(fHistMassWithCombPIDOmegaPlus); |
895 | } | |
896 | ||
45ee0bcc | 897 | |
898 | ||
899 | // - Complements for QA | |
900 | ||
901 | if(! fHistXiTransvMom ){ | |
90926f02 | 902 | fHistXiTransvMom = new TH1F( "fHistXiTransvMom" , "#Xi transverse momentum (cand. around the mass peak) ; p_{t}(#Xi) (GeV/c); Counts", 100, 0.0, 10.0); |
45ee0bcc | 903 | fListHistCascade->Add(fHistXiTransvMom); |
904 | } | |
905 | ||
906 | if(! fHistXiTotMom ){ | |
90926f02 | 907 | fHistXiTotMom = new TH1F( "fHistXiTotMom" , "#Xi momentum norm (cand. around the mass peak); p_{tot}(#Xi) (GeV/c); Counts", 150, 0.0, 15.0); |
45ee0bcc | 908 | fListHistCascade->Add(fHistXiTotMom); |
909 | } | |
910 | ||
911 | ||
90926f02 | 912 | if(! fHistBachTransvMomXi ){ |
913 | fHistBachTransvMomXi = new TH1F( "fHistBachTransvMomXi" , "#Xi Bach. transverse momentum (cand. around the mass peak) ; p_{t}(Bach.) (GeV/c); Counts", 100, 0.0, 5.0); | |
914 | fListHistCascade->Add(fHistBachTransvMomXi); | |
45ee0bcc | 915 | } |
916 | ||
90926f02 | 917 | if(! fHistBachTotMomXi ){ |
918 | fHistBachTotMomXi = new TH1F( "fHistBachTotMomXi" , "#Xi Bach. momentum norm (cand. around the mass peak); p_{tot}(Bach.) (GeV/c); Counts", 100, 0.0, 5.0); | |
919 | fListHistCascade->Add(fHistBachTotMomXi); | |
45ee0bcc | 920 | } |
921 | ||
922 | ||
923 | if(! fHistChargeXi ){ | |
924 | fHistChargeXi = new TH1F( "fHistChargeXi" , "Charge of casc. candidates ; Sign ; Counts", 5, -2.0, 3.0); | |
925 | fListHistCascade->Add(fHistChargeXi); | |
926 | } | |
927 | ||
928 | ||
929 | if (! fHistV0toXiCosineOfPointingAngle) { | |
930 | fHistV0toXiCosineOfPointingAngle = new TH1F("fHistV0toXiCosineOfPointingAngle", "Cos. of V0 Ptng Angl / Xi vtx ;Cos(V0 Point. Angl / Xi vtx); Counts", 100, 0.99, 1.0); | |
931 | fListHistCascade->Add(fHistV0toXiCosineOfPointingAngle); | |
932 | } | |
933 | ||
934 | ||
935 | if(! fHistRapXi ){ | |
90926f02 | 936 | fHistRapXi = new TH1F( "fHistRapXi" , "Rapidity of #Xi candidates (around the mass peak); y ; Counts", 200, -5.0, 5.0); |
45ee0bcc | 937 | fListHistCascade->Add(fHistRapXi); |
938 | } | |
939 | ||
940 | if(! fHistRapOmega ){ | |
90926f02 | 941 | fHistRapOmega = new TH1F( "fHistRapOmega" , "Rapidity of #Omega candidates (around the mass peak); y ; Counts", 200, -5.0, 5.0); |
45ee0bcc | 942 | fListHistCascade->Add(fHistRapOmega); |
943 | } | |
944 | ||
90926f02 | 945 | if(! fHistEtaXi ){ |
946 | fHistEtaXi = new TH1F( "fHistEtaXi" , "Pseudo-rap. of #Xi candidates (around the mass peak) ; #eta ; Counts", 120, -3.0, 3.0); | |
947 | fListHistCascade->Add(fHistEtaXi); | |
45ee0bcc | 948 | } |
949 | ||
90926f02 | 950 | if(! fHistThetaXi ){ |
951 | fHistThetaXi = new TH1F( "fHistThetaXi" , "#theta of #Xi candidates (around the mass peak); #theta (deg) ; Counts", 180, 0., 180.0); | |
952 | fListHistCascade->Add(fHistThetaXi); | |
45ee0bcc | 953 | } |
954 | ||
90926f02 | 955 | if(! fHistPhiXi ){ |
956 | fHistPhiXi = new TH1F( "fHistPhiXi" , "#phi of #Xi candidates (around the mass peak); #phi (deg) ; Counts", 360, 0., 360.); | |
957 | fListHistCascade->Add(fHistPhiXi); | |
45ee0bcc | 958 | } |
959 | ||
960 | ||
c525f360 | 961 | if(! fHistcTauXiMinus){ |
962 | fHistcTauXiMinus = new TH1F("fHistcTauXiMinus", "Lifetime c.#tau for #Xi^{-}; L_{3D}.m_{PDG}(#Xi^{-}) / p_{3D} (cm); Counts", 100, 0., 50.); | |
963 | fListHistCascade->Add(fHistcTauXiMinus); | |
964 | } | |
965 | ||
966 | if(! fHistcTauXiPlus){ | |
a786bd57 | 967 | fHistcTauXiPlus = new TH1F("fHistcTauXiPlus", "Lifetime c.#tau for #bar{#Xi}^{+}; L_{3D}.m_{PDG}(#bar{#Xi}^{+}) / p_{3D} (cm); Counts", 100, 0., 50.); |
c525f360 | 968 | fListHistCascade->Add(fHistcTauXiPlus); |
969 | } | |
970 | ||
971 | if(! fHistcTauOmegaMinus){ | |
972 | fHistcTauOmegaMinus = new TH1F("fHistcTauOmegaMinus", "Lifetime c.#tau for #Omega^{-}; L_{3D}.m_{PDG}(#Omega^{-}) / p_{3D} (cm); Counts", 100, 0., 50.); | |
973 | fListHistCascade->Add(fHistcTauOmegaMinus); | |
974 | } | |
975 | ||
976 | if(! fHistcTauOmegaPlus){ | |
a786bd57 | 977 | fHistcTauOmegaPlus = new TH1F("fHistcTauOmegaPlus", "Lifetime c.#tau for #bar{#Omega}^{+}; L_{3D}.m_{PDG}(#bar{#Omega}^{+}) / p_{3D} (cm); Counts", 100, 0., 50.); |
c525f360 | 978 | fListHistCascade->Add(fHistcTauOmegaPlus); |
979 | } | |
980 | ||
981 | ||
45ee0bcc | 982 | if(! f2dHistArmenteros) { |
983 | 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); | |
984 | fListHistCascade->Add(f2dHistArmenteros); | |
985 | } | |
986 | ||
e9c3f0b6 | 987 | //------- |
988 | ||
45ee0bcc | 989 | if(! f2dHistEffMassLambdaVsEffMassXiMinus) { |
90926f02 | 990 | 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); |
45ee0bcc | 991 | fListHistCascade->Add(f2dHistEffMassLambdaVsEffMassXiMinus); |
992 | } | |
993 | ||
994 | if(! f2dHistEffMassXiVsEffMassOmegaMinus) { | |
90926f02 | 995 | 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); |
45ee0bcc | 996 | fListHistCascade->Add(f2dHistEffMassXiVsEffMassOmegaMinus); |
997 | } | |
998 | ||
999 | if(! f2dHistEffMassLambdaVsEffMassXiPlus) { | |
a786bd57 | 1000 | f2dHistEffMassLambdaVsEffMassXiPlus = new TH2F( "f2dHistEffMassLambdaVsEffMassXiPlus", "M_{#Lambda} Vs M_{#bar{#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); |
45ee0bcc | 1001 | fListHistCascade->Add(f2dHistEffMassLambdaVsEffMassXiPlus); |
1002 | } | |
1003 | ||
1004 | if(! f2dHistEffMassXiVsEffMassOmegaPlus) { | |
a786bd57 | 1005 | f2dHistEffMassXiVsEffMassOmegaPlus = new TH2F( "f2dHistEffMassXiVsEffMassOmegaPlus", "M_{#bar{#Xi}^{+} candidates} Vs M_{#bar{#Omega}^{+} candidates} ; M( #Lambda , #pi^{+} ) (GeV/c^{2}) ; M( #Lambda , K^{+} ) (GeV/c^{2})", 400, 1.2, 2.0, 500, 1.5, 2.5); |
45ee0bcc | 1006 | fListHistCascade->Add(f2dHistEffMassXiVsEffMassOmegaPlus); |
1007 | } | |
1008 | ||
e9c3f0b6 | 1009 | //------- |
1010 | ||
45ee0bcc | 1011 | if(! f2dHistXiRadiusVsEffMassXiMinus) { |
90926f02 | 1012 | 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); |
45ee0bcc | 1013 | fListHistCascade->Add(f2dHistXiRadiusVsEffMassXiMinus); |
1014 | } | |
1015 | ||
1016 | if(! f2dHistXiRadiusVsEffMassXiPlus) { | |
a786bd57 | 1017 | f2dHistXiRadiusVsEffMassXiPlus = new TH2F( "f2dHistXiRadiusVsEffMassXiPlus", "Transv. R_{Xi Decay} Vs M_{#bar{#Xi}^{+} candidates}; r_{cascade} (cm); M( #Lambda , #pi^{+} ) (GeV/c^{2}) ", 450, 0., 45.0, 400, 1.2, 2.0); |
45ee0bcc | 1018 | fListHistCascade->Add(f2dHistXiRadiusVsEffMassXiPlus); |
1019 | } | |
1020 | ||
e9c3f0b6 | 1021 | if(! f2dHistXiRadiusVsEffMassOmegaMinus) { |
90926f02 | 1022 | 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); |
e9c3f0b6 | 1023 | fListHistCascade->Add(f2dHistXiRadiusVsEffMassOmegaMinus); |
1024 | } | |
1025 | ||
1026 | if(! f2dHistXiRadiusVsEffMassOmegaPlus) { | |
a786bd57 | 1027 | f2dHistXiRadiusVsEffMassOmegaPlus = new TH2F( "f2dHistXiRadiusVsEffMassOmegaPlus", "Transv. R_{Xi Decay} Vs M_{#bar{#Omega}^{+} candidates}; r_{cascade} (cm); M( #Lambda , K^{+} ) (GeV/c^{2}) ", 450, 0., 45.0, 500, 1.5, 2.5); |
e9c3f0b6 | 1028 | fListHistCascade->Add(f2dHistXiRadiusVsEffMassOmegaPlus); |
1029 | } | |
1030 | ||
e8fd19de | 1031 | //------ |
1032 | ||
1033 | if(! f2dHistTPCdEdxOfCascDghters){ | |
1034 | 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.); | |
1035 | fListHistCascade->Add(f2dHistTPCdEdxOfCascDghters); | |
1036 | } | |
1037 | ||
1038 | ||
e9c3f0b6 | 1039 | |
3aedd4a5 | 1040 | // Part 2 : Raw material for yield extraction ------- |
45ee0bcc | 1041 | |
ff0753e1 | 1042 | if(! f3dHistXiPtVsEffMassVsYXiMinus) { |
e8fd19de | 1043 | f3dHistXiPtVsEffMassVsYXiMinus = new TH3F( "f3dHistXiPtVsEffMassVsYXiMinus", "Pt_{cascade} Vs M_{#Xi^{-} candidates} Vs Y_{#Xi}; Pt_{cascade} (GeV/c); M( #Lambda , #pi^{-} ) (GeV/c^{2}) ;Y_{#Xi} ", 100, 0., 10.0, 400, 1.2, 2.0, 44, -1.1,1.1); |
ff0753e1 | 1044 | fListHistCascade->Add(f3dHistXiPtVsEffMassVsYXiMinus); |
e9c3f0b6 | 1045 | } |
1046 | ||
ff0753e1 | 1047 | if(! f3dHistXiPtVsEffMassVsYXiPlus) { |
a786bd57 | 1048 | f3dHistXiPtVsEffMassVsYXiPlus = new TH3F( "f3dHistXiPtVsEffMassVsYXiPlus", "Pt_{cascade} Vs M_{#bar{#Xi}^{+} candidates} Vs Y_{#Xi}; Pt_{cascade} (GeV/c); M( #Lambda , #pi^{+} ) (GeV/c^{2}); Y_{#Xi}", 100, 0., 10.0, 400, 1.2, 2.0, 44, -1.1,1.1); |
ff0753e1 | 1049 | fListHistCascade->Add(f3dHistXiPtVsEffMassVsYXiPlus); |
e9c3f0b6 | 1050 | } |
1051 | ||
ff0753e1 | 1052 | if(! f3dHistXiPtVsEffMassVsYOmegaMinus) { |
e8fd19de | 1053 | f3dHistXiPtVsEffMassVsYOmegaMinus = new TH3F( "f3dHistXiPtVsEffMassVsYOmegaMinus", "Pt_{cascade} Vs M_{#Omega^{-} candidates} Vs Y_{#Omega}; Pt_{cascade} (GeV/c); M( #Lambda , K^{-} ) (GeV/c^{2}); Y_{#Omega}", 100, 0., 10.0, 500, 1.5, 2.5, 44, -1.1,1.1); |
ff0753e1 | 1054 | fListHistCascade->Add(f3dHistXiPtVsEffMassVsYOmegaMinus); |
e9c3f0b6 | 1055 | } |
1056 | ||
ff0753e1 | 1057 | if(! f3dHistXiPtVsEffMassVsYOmegaPlus) { |
a786bd57 | 1058 | f3dHistXiPtVsEffMassVsYOmegaPlus = new TH3F( "f3dHistXiPtVsEffMassVsYOmegaPlus", "Pt_{cascade} Vs M_{#bar{#Omega}^{+} candidates} Vs Y_{#Omega}; Pt_{cascade} (GeV/c); M( #Lambda , K^{+} ) (GeV/c^{2}); Y_{#Omega}", 100, 0., 10.0, 500, 1.5, 2.5, 44, -1.1,1.1); |
ff0753e1 | 1059 | fListHistCascade->Add(f3dHistXiPtVsEffMassVsYOmegaPlus); |
e9c3f0b6 | 1060 | } |
45ee0bcc | 1061 | |
3aedd4a5 | 1062 | //-- |
32599cd4 | 1063 | if(! fCFContCascadePIDXiMinus) { |
1064 | const Int_t lNbSteps = 7 ; | |
1065 | const Int_t lNbVariables = 4 ; | |
1066 | ||
1067 | //array for the number of bins in each dimension : | |
e8fd19de | 1068 | Int_t lNbBinsPerVar[4] = {0}; |
32599cd4 | 1069 | lNbBinsPerVar[0] = 100; |
a786bd57 | 1070 | lNbBinsPerVar[1] = 75; |
e8fd19de | 1071 | lNbBinsPerVar[2] = 44; |
f3dc9369 | 1072 | lNbBinsPerVar[3] = 250; |
32599cd4 | 1073 | |
1074 | ||
1075 | fCFContCascadePIDXiMinus = new AliCFContainer("fCFContCascadePIDXiMinus","Pt_{cascade} Vs M_{#Xi^{-} candidates} Vs Y_{#Xi}", lNbSteps, lNbVariables, lNbBinsPerVar ); | |
1076 | ||
1077 | //setting the bin limits (valid for v4-18-10-AN) | |
a786bd57 | 1078 | fCFContCascadePIDXiMinus->SetBinLimits(0, 0.0 , 10.0 ); // Pt(Cascade) |
1079 | fCFContCascadePIDXiMinus->SetBinLimits(1, 1.25 , 1.40 ); // Xi Effective mass | |
1080 | fCFContCascadePIDXiMinus->SetBinLimits(2, -1.1 , 1.1 ); // Rapidity | |
32599cd4 | 1081 | if(fCollidingSystems) |
c525f360 | 1082 | fCFContCascadePIDXiMinus->SetBinLimits(3, 0.0, 20000.0 ); // Primary track Multiplicity |
32599cd4 | 1083 | else |
c525f360 | 1084 | fCFContCascadePIDXiMinus->SetBinLimits(3, 0.0, 250.0 ); // Primary track Multiplicity |
32599cd4 | 1085 | |
1086 | // Setting the step title : one per PID case | |
1087 | fCFContCascadePIDXiMinus->SetStepTitle(0, "No PID"); | |
e8fd19de | 1088 | fCFContCascadePIDXiMinus->SetStepTitle(1, "TPC PID / 4-#sigma cut on Bachelor track"); |
1089 | fCFContCascadePIDXiMinus->SetStepTitle(2, "TPC PID / 4-#sigma cut on Bachelor+Baryon tracks"); | |
1090 | fCFContCascadePIDXiMinus->SetStepTitle(3, "TPC PID / 4-#sigma cut on Bachelor+Baryon+Meson tracks"); | |
32599cd4 | 1091 | fCFContCascadePIDXiMinus->SetStepTitle(4, "Comb. PID / Bachelor"); |
1092 | fCFContCascadePIDXiMinus->SetStepTitle(5, "Comb. PID / Bachelor+Baryon"); | |
1093 | fCFContCascadePIDXiMinus->SetStepTitle(6, "Comb. PID / Bachelor+Baryon+Meson"); | |
1094 | ||
1095 | // Setting the variable title, per axis | |
1096 | fCFContCascadePIDXiMinus->SetVarTitle(0, "Pt_{cascade} (GeV/c)"); | |
1097 | fCFContCascadePIDXiMinus->SetVarTitle(1, "M( #Lambda , #pi^{-} ) (GeV/c^{2})"); | |
1098 | fCFContCascadePIDXiMinus->SetVarTitle(2, "Y_{#Xi}"); | |
c525f360 | 1099 | fCFContCascadePIDXiMinus->SetVarTitle(3, "Primary track Multiplicity"); |
32599cd4 | 1100 | |
1101 | fListHistCascade->Add(fCFContCascadePIDXiMinus); | |
1102 | ||
1103 | } | |
1104 | ||
1105 | if(! fCFContCascadePIDXiPlus) { | |
1106 | const Int_t lNbSteps = 7 ; | |
1107 | const Int_t lNbVariables = 4 ; | |
1108 | ||
1109 | //array for the number of bins in each dimension : | |
e8fd19de | 1110 | Int_t lNbBinsPerVar[4] = {0}; |
32599cd4 | 1111 | lNbBinsPerVar[0] = 100; |
a786bd57 | 1112 | lNbBinsPerVar[1] = 75; |
e8fd19de | 1113 | lNbBinsPerVar[2] = 44; |
f3dc9369 | 1114 | lNbBinsPerVar[3] = 250; |
32599cd4 | 1115 | |
1116 | ||
a786bd57 | 1117 | fCFContCascadePIDXiPlus = new AliCFContainer("fCFContCascadePIDXiPlus","Pt_{cascade} Vs M_{#bar{#Xi}^{+} candidates} Vs Y_{#Xi}", lNbSteps, lNbVariables, lNbBinsPerVar ); |
32599cd4 | 1118 | |
1119 | ||
1120 | //setting the bin limits (valid for v4-18-10-AN) | |
a786bd57 | 1121 | fCFContCascadePIDXiPlus->SetBinLimits(0, 0.0 , 10.0 ); // Pt(Cascade) |
1122 | fCFContCascadePIDXiPlus->SetBinLimits(1, 1.25 , 1.40 ); // Xi Effective mass | |
1123 | fCFContCascadePIDXiPlus->SetBinLimits(2, -1.1 , 1.1 ); // Rapidity | |
32599cd4 | 1124 | if(fCollidingSystems) |
c525f360 | 1125 | fCFContCascadePIDXiPlus->SetBinLimits(3, 0.0, 20000.0 ); // Primary track Multiplicity |
32599cd4 | 1126 | else |
c525f360 | 1127 | fCFContCascadePIDXiPlus->SetBinLimits(3, 0.0, 250.0 ); // Primary track Multiplicity |
32599cd4 | 1128 | |
1129 | // Setting the step title : one per PID case | |
1130 | fCFContCascadePIDXiPlus->SetStepTitle(0, "No PID"); | |
e8fd19de | 1131 | fCFContCascadePIDXiPlus->SetStepTitle(1, "TPC PID / 4-#sigma cut on Bachelor track"); |
1132 | fCFContCascadePIDXiPlus->SetStepTitle(2, "TPC PID / 4-#sigma cut on Bachelor+Baryon tracks"); | |
1133 | fCFContCascadePIDXiPlus->SetStepTitle(3, "TPC PID / 4-#sigma cut on Bachelor+Baryon+Meson tracks"); | |
32599cd4 | 1134 | fCFContCascadePIDXiPlus->SetStepTitle(4, "Comb. PID / Bachelor"); |
1135 | fCFContCascadePIDXiPlus->SetStepTitle(5, "Comb. PID / Bachelor+Baryon"); | |
1136 | fCFContCascadePIDXiPlus->SetStepTitle(6, "Comb. PID / Bachelor+Baryon+Meson"); | |
1137 | ||
1138 | // Setting the variable title, per axis | |
1139 | fCFContCascadePIDXiPlus->SetVarTitle(0, "Pt_{cascade} (GeV/c)"); | |
1140 | fCFContCascadePIDXiPlus->SetVarTitle(1, "M( #Lambda , #pi^{+} ) (GeV/c^{2})"); | |
1141 | fCFContCascadePIDXiPlus->SetVarTitle(2, "Y_{#Xi}"); | |
c525f360 | 1142 | fCFContCascadePIDXiPlus->SetVarTitle(3, "Primary track Multiplicity"); |
32599cd4 | 1143 | |
1144 | fListHistCascade->Add(fCFContCascadePIDXiPlus); | |
1145 | ||
1146 | } | |
1147 | ||
1148 | ||
1149 | if(! fCFContCascadePIDOmegaMinus) { | |
1150 | const Int_t lNbSteps = 7 ; | |
1151 | const Int_t lNbVariables = 4 ; | |
1152 | ||
1153 | //array for the number of bins in each dimension : | |
e8fd19de | 1154 | Int_t lNbBinsPerVar[4] = {0}; |
32599cd4 | 1155 | lNbBinsPerVar[0] = 100; |
a786bd57 | 1156 | lNbBinsPerVar[1] = 60; |
e8fd19de | 1157 | lNbBinsPerVar[2] = 44; |
f3dc9369 | 1158 | lNbBinsPerVar[3] = 250; |
32599cd4 | 1159 | |
1160 | ||
1161 | fCFContCascadePIDOmegaMinus = new AliCFContainer("fCFContCascadePIDOmegaMinus","Pt_{cascade} Vs M_{#Omega^{-} candidates} Vs Y_{#Omega}", lNbSteps, lNbVariables, lNbBinsPerVar ); | |
1162 | ||
1163 | ||
1164 | //setting the bin limits (valid for v4-18-10-AN) | |
a786bd57 | 1165 | fCFContCascadePIDOmegaMinus->SetBinLimits(0, 0.0 , 10.0 ); // Pt(Cascade) |
1166 | fCFContCascadePIDOmegaMinus->SetBinLimits(1, 1.62 , 1.74 ); // Omega Effective mass | |
1167 | fCFContCascadePIDOmegaMinus->SetBinLimits(2, -1.1 , 1.1 ); // Rapidity | |
32599cd4 | 1168 | if(fCollidingSystems) |
c525f360 | 1169 | fCFContCascadePIDOmegaMinus->SetBinLimits(3, 0.0, 20000.0 ); //Primary track Multiplicity |
32599cd4 | 1170 | else |
c525f360 | 1171 | fCFContCascadePIDOmegaMinus->SetBinLimits(3, 0.0, 250.0 ); // Primary track Multiplicity |
32599cd4 | 1172 | |
1173 | // Setting the step title : one per PID case | |
1174 | fCFContCascadePIDOmegaMinus->SetStepTitle(0, "No PID"); | |
e8fd19de | 1175 | fCFContCascadePIDOmegaMinus->SetStepTitle(1, "TPC PID / 4-#sigma cut on Bachelor track"); |
1176 | fCFContCascadePIDOmegaMinus->SetStepTitle(2, "TPC PID / 4-#sigma cut on Bachelor+Baryon tracks"); | |
1177 | fCFContCascadePIDOmegaMinus->SetStepTitle(3, "TPC PID / 4-#sigma cut on Bachelor+Baryon+Meson tracks"); | |
32599cd4 | 1178 | fCFContCascadePIDOmegaMinus->SetStepTitle(4, "Comb. PID / Bachelor"); |
1179 | fCFContCascadePIDOmegaMinus->SetStepTitle(5, "Comb. PID / Bachelor+Baryon"); | |
1180 | fCFContCascadePIDOmegaMinus->SetStepTitle(6, "Comb. PID / Bachelor+Baryon+Meson"); | |
1181 | ||
1182 | // Setting the variable title, per axis | |
1183 | fCFContCascadePIDOmegaMinus->SetVarTitle(0, "Pt_{cascade} (GeV/c)"); | |
1184 | fCFContCascadePIDOmegaMinus->SetVarTitle(1, "M( #Lambda , K^{-} ) (GeV/c^{2})"); | |
1185 | fCFContCascadePIDOmegaMinus->SetVarTitle(2, "Y_{#Omega}"); | |
c525f360 | 1186 | fCFContCascadePIDOmegaMinus->SetVarTitle(3, "Primary track Multiplicity"); |
32599cd4 | 1187 | |
1188 | fListHistCascade->Add(fCFContCascadePIDOmegaMinus); | |
1189 | ||
1190 | } | |
1191 | ||
1192 | if(! fCFContCascadePIDOmegaPlus) { | |
1193 | const Int_t lNbSteps = 7 ; | |
1194 | const Int_t lNbVariables = 4 ; | |
1195 | ||
1196 | //array for the number of bins in each dimension : | |
e8fd19de | 1197 | Int_t lNbBinsPerVar[4] = {0}; |
32599cd4 | 1198 | lNbBinsPerVar[0] = 100; |
a786bd57 | 1199 | lNbBinsPerVar[1] = 60; |
e8fd19de | 1200 | lNbBinsPerVar[2] = 44; |
f3dc9369 | 1201 | lNbBinsPerVar[3] = 250; |
32599cd4 | 1202 | |
1203 | ||
a786bd57 | 1204 | fCFContCascadePIDOmegaPlus = new AliCFContainer("fCFContCascadePIDOmegaPlus","Pt_{cascade} Vs M_{#bar{#Omega}^{+} candidates} Vs Y_{#Omega}", lNbSteps, lNbVariables, lNbBinsPerVar ); |
32599cd4 | 1205 | |
1206 | ||
1207 | //setting the bin limits (valid for v4-18-10-AN) | |
a786bd57 | 1208 | fCFContCascadePIDOmegaPlus->SetBinLimits(0, 0.0 , 10.0 ); // Pt(Cascade) |
1209 | fCFContCascadePIDOmegaPlus->SetBinLimits(1, 1.62 , 1.74 ); // Omega Effective mass | |
1210 | fCFContCascadePIDOmegaPlus->SetBinLimits(2, -1.1 , 1.1 ); // Rapidity | |
32599cd4 | 1211 | if(fCollidingSystems) |
c525f360 | 1212 | fCFContCascadePIDOmegaPlus->SetBinLimits(3, 0.0, 20000.0 ); // Primary track Multiplicity |
32599cd4 | 1213 | else |
c525f360 | 1214 | fCFContCascadePIDOmegaPlus->SetBinLimits(3, 0.0, 250.0 ); // Primary track Multiplicity |
32599cd4 | 1215 | |
1216 | // Setting the step title : one per PID case | |
1217 | fCFContCascadePIDOmegaPlus->SetStepTitle(0, "No PID"); | |
e8fd19de | 1218 | fCFContCascadePIDOmegaPlus->SetStepTitle(1, "TPC PID / 4-#sigma cut on Bachelor track"); |
1219 | fCFContCascadePIDOmegaPlus->SetStepTitle(2, "TPC PID / 4-#sigma cut on Bachelor+Baryon tracks"); | |
1220 | fCFContCascadePIDOmegaPlus->SetStepTitle(3, "TPC PID / 4-#sigma cut on Bachelor+Baryon+Meson tracks"); | |
32599cd4 | 1221 | fCFContCascadePIDOmegaPlus->SetStepTitle(4, "Comb. PID / Bachelor"); |
1222 | fCFContCascadePIDOmegaPlus->SetStepTitle(5, "Comb. PID / Bachelor+Baryon"); | |
1223 | fCFContCascadePIDOmegaPlus->SetStepTitle(6, "Comb. PID / Bachelor+Baryon+Meson"); | |
1224 | ||
1225 | // Setting the variable title, per axis | |
1226 | fCFContCascadePIDOmegaPlus->SetVarTitle(0, "Pt_{cascade} (GeV/c)"); | |
1227 | fCFContCascadePIDOmegaPlus->SetVarTitle(1, "M( #Lambda , K^{+} ) (GeV/c^{2})"); | |
1228 | fCFContCascadePIDOmegaPlus->SetVarTitle(2, "Y_{#Omega}"); | |
c525f360 | 1229 | fCFContCascadePIDOmegaPlus->SetVarTitle(3, "Primary track Multiplicity"); |
32599cd4 | 1230 | |
1231 | fListHistCascade->Add(fCFContCascadePIDOmegaPlus); | |
1232 | ||
1233 | } | |
1234 | ||
1235 | ||
1236 | ||
1237 | ||
1238 | ||
32599cd4 | 1239 | // Part 3 : Towards the optimisation of topological selections ------- |
1240 | if(! fCFContCascadeCuts){ | |
1241 | ||
1242 | // Container meant to store all the relevant distributions corresponding to the cut variables. | |
9d670198 | 1243 | // So far, 20 variables have been identified. |
1244 | // The following will be done in quite a brut force way ... | |
1245 | // FIXME Improvement expected later (before Pb-Pb data at least) | |
1246 | // - Define a user binning to have less bins in each dimension | |
1247 | // - boolean for enabling/disbaling this CFContainer | |
1248 | const Int_t lNbSteps = 4 ; | |
1249 | const Int_t lNbVariables = 20 ; | |
32599cd4 | 1250 | |
1251 | //array for the number of bins in each dimension : | |
e8fd19de | 1252 | Int_t lNbBinsPerVar[20] = {0}; |
32599cd4 | 1253 | lNbBinsPerVar[0] = 25; |
1254 | lNbBinsPerVar[1] = 25; | |
1255 | lNbBinsPerVar[2] = 20; | |
1256 | lNbBinsPerVar[3] = 40; | |
e8fd19de | 1257 | lNbBinsPerVar[4] = 30; |
1258 | lNbBinsPerVar[5] = 25; | |
32599cd4 | 1259 | |
1260 | lNbBinsPerVar[6] = 20; | |
1261 | lNbBinsPerVar[7] = 40; | |
1262 | lNbBinsPerVar[8] = 40; | |
1263 | lNbBinsPerVar[9] = 25; | |
1264 | lNbBinsPerVar[10] = 25; | |
1265 | ||
e8fd19de | 1266 | lNbBinsPerVar[11] = 75; // 2-MeV/c2 bins |
1267 | lNbBinsPerVar[12] = 60; // 2-MeV/c2 bins | |
32599cd4 | 1268 | |
e8fd19de | 1269 | lNbBinsPerVar[13] = 100; |
9d670198 | 1270 | |
1271 | lNbBinsPerVar[14] = 44; // 0.05 in rapidity units | |
1272 | lNbBinsPerVar[15] = 44; // 0.05 in rapidity units | |
1273 | ||
1274 | lNbBinsPerVar[16] = 20; | |
32599cd4 | 1275 | |
9d670198 | 1276 | lNbBinsPerVar[17] = 50; |
e8fd19de | 1277 | lNbBinsPerVar[18] = 100; |
1278 | lNbBinsPerVar[19] = 24; | |
32599cd4 | 1279 | |
1280 | fCFContCascadeCuts = new AliCFContainer("fCFContCascadeCuts","Container for Cascade cuts", lNbSteps, lNbVariables, lNbBinsPerVar ); | |
1281 | ||
1282 | ||
e8fd19de | 1283 | //setting the bin limits |
1284 | ||
1285 | //0 | |
1286 | Double_t *lBinLim0 = new Double_t[ lNbBinsPerVar[0]+1 ]; | |
1287 | for(Int_t i=0; i< lNbBinsPerVar[0];i++) lBinLim0[i] = (Double_t)0.0 + (4.8 - 0.0 )/(lNbBinsPerVar[0]-1) * (Double_t)i ; | |
1288 | lBinLim0[ lNbBinsPerVar[0] ] = 20.0; | |
1289 | fCFContCascadeCuts -> SetBinLimits(0, lBinLim0 ); // DcaXiDaughters : 0.0 to 5.0 | |
1290 | delete [] lBinLim0; | |
1291 | //1 | |
1292 | Double_t *lBinLim1 = new Double_t[ lNbBinsPerVar[1]+1 ]; | |
1293 | for(Int_t i=0; i< lNbBinsPerVar[1];i++) lBinLim1[i] = (Double_t)0.0 + (0.24 - 0.0 )/(lNbBinsPerVar[1]-1) * (Double_t)i ; | |
1294 | lBinLim1[ lNbBinsPerVar[1] ] = 100.0; | |
1295 | fCFContCascadeCuts -> SetBinLimits(1, lBinLim1 ); // DcaBachToPrimVertexXi : 0.0 to 0.25 | |
1296 | delete [] lBinLim1; | |
1297 | //2 | |
1298 | Double_t *lBinLim2 = new Double_t[ lNbBinsPerVar[2]+1 ]; | |
1299 | for(Int_t i=1; i< lNbBinsPerVar[2]+1;i++) lBinLim2[i] = (Double_t)0.81 + (1.0 - 0.81 )/(lNbBinsPerVar[2]-1) * (Double_t) (i-1) ; | |
1300 | lBinLim2[0] = 0.0; | |
1301 | fCFContCascadeCuts -> SetBinLimits(2, lBinLim2 ); // XiCosineOfPointingAngle : 0.80 to 1.0 | |
1302 | delete [] lBinLim2; | |
1303 | //3 | |
1304 | Double_t *lBinLim3 = new Double_t[ lNbBinsPerVar[3]+1 ]; | |
1305 | for(Int_t i=0; i< lNbBinsPerVar[3];i++) lBinLim3[i] = (Double_t)0.0 + (3.9 - 0.0 )/(lNbBinsPerVar[3]-1) * (Double_t)i ; | |
1306 | lBinLim3[ lNbBinsPerVar[3] ] = 110.0; | |
1307 | fCFContCascadeCuts -> SetBinLimits(3, lBinLim3 ); // XiRadius : 0.0 to 4.0 | |
1308 | delete [] lBinLim3; | |
1309 | //4 | |
1310 | fCFContCascadeCuts->SetBinLimits(4, 1.1 , 1.13 ); // InvMassLambdaAsCascDghter | |
1311 | //5 | |
1312 | Double_t *lBinLim5 = new Double_t[ lNbBinsPerVar[5]+1 ]; | |
1313 | for(Int_t i=0; i< lNbBinsPerVar[5];i++) lBinLim5[i] = (Double_t)0.0 + (4.8 - 0.0 )/(lNbBinsPerVar[5]-1) * (Double_t)i ; | |
1314 | lBinLim5[ lNbBinsPerVar[5] ] = 20.0; | |
1315 | fCFContCascadeCuts -> SetBinLimits(5, lBinLim5 ); // DcaV0DaughtersXi : 0.0 to 5.0 | |
1316 | delete [] lBinLim5; | |
1317 | ||
1318 | ||
1319 | //6 | |
1320 | Double_t *lBinLim6 = new Double_t[ lNbBinsPerVar[6]+1 ]; | |
1321 | for(Int_t i=1; i< lNbBinsPerVar[6]+1 ;i++) lBinLim6[i] = (Double_t)0.81 + (1.0 - 0.81 )/(lNbBinsPerVar[6]-1) * (Double_t) (i-1) ; | |
1322 | lBinLim6[0] = 0.0; | |
1323 | fCFContCascadeCuts -> SetBinLimits(6, lBinLim6 ); // V0CosineOfPointingAngleXi : 0.80 to 1.0 | |
1324 | delete [] lBinLim6; | |
1325 | //7 | |
1326 | Double_t *lBinLim7 = new Double_t[ lNbBinsPerVar[7]+1 ]; | |
1327 | for(Int_t i=0; i< lNbBinsPerVar[7];i++) lBinLim7[i] = (Double_t)0.0 + (7.8 - 0.0 )/(lNbBinsPerVar[7]-1) * (Double_t)i ; | |
1328 | lBinLim7[ lNbBinsPerVar[7] ] = 100.0; | |
1329 | fCFContCascadeCuts -> SetBinLimits(7, lBinLim7 ); // V0RadiusXi : 0.0 to 8.0 | |
1330 | delete [] lBinLim7; | |
1331 | //8 | |
1332 | Double_t *lBinLim8 = new Double_t[ lNbBinsPerVar[8]+1 ]; | |
1333 | for(Int_t i=0; i< lNbBinsPerVar[8];i++) lBinLim8[i] = (Double_t)0.0 + (0.39 - 0.0 )/(lNbBinsPerVar[8]-1) * (Double_t)i ; | |
1334 | lBinLim8[ lNbBinsPerVar[8] ] = 100.0; | |
1335 | fCFContCascadeCuts -> SetBinLimits(8, lBinLim8 ); // DcaV0ToPrimVertexXi : 0.0 to 0.40 | |
1336 | delete [] lBinLim8; | |
1337 | //9 | |
1338 | Double_t *lBinLim9 = new Double_t[ lNbBinsPerVar[9]+1 ]; | |
1339 | for(Int_t i=0; i< lNbBinsPerVar[9];i++) lBinLim9[i] = (Double_t)0.0 + (0.24 - 0.0 )/(lNbBinsPerVar[9]-1) * (Double_t)i ; | |
1340 | lBinLim9[ lNbBinsPerVar[9] ] = 100.0; | |
1341 | fCFContCascadeCuts -> SetBinLimits(9, lBinLim9 ); // DcaPosToPrimVertexXi : 0.0 to 0.25 | |
1342 | delete [] lBinLim9; | |
1343 | //10 | |
1344 | Double_t *lBinLim10 = new Double_t[ lNbBinsPerVar[10]+1 ]; | |
1345 | for(Int_t i=0; i< lNbBinsPerVar[10];i++) lBinLim10[i] = (Double_t)0.0 + (0.24 - 0.0 )/(lNbBinsPerVar[10]-1) * (Double_t)i ; | |
1346 | lBinLim10[ lNbBinsPerVar[10] ] = 100.0; | |
1347 | fCFContCascadeCuts -> SetBinLimits(10, lBinLim10 ); // DcaPosToPrimVertexXi : 0.0 to 0.25 | |
1348 | delete [] lBinLim10; | |
1349 | ||
1350 | //11 | |
1351 | fCFContCascadeCuts->SetBinLimits(11, 1.25 , 1.40 ); // InvMassXi | |
1352 | fCFContCascadeCuts->SetBinLimits(12, 1.62 , 1.74 ); // InvMassOmega | |
1353 | fCFContCascadeCuts->SetBinLimits(13, 0.0 , 10.0 ); // XiTransvMom | |
1354 | fCFContCascadeCuts->SetBinLimits(14, -1.1 , 1.1 ); // Y(Xi) | |
1355 | fCFContCascadeCuts->SetBinLimits(15, -1.1 , 1.1 ); // Y(Omega) | |
1356 | fCFContCascadeCuts->SetBinLimits(16, -10.0 , 10.0 ); // BestPrimaryVtxPosZ | |
1357 | if(fCollidingSystems){ | |
c525f360 | 1358 | fCFContCascadeCuts->SetBinLimits(17, 0.0, 10000.0 ); // nTrackPrimaryMultiplicity |
a786bd57 | 1359 | fCFContCascadeCuts->SetBinLimits(18, 0.0, 10000.0 ); // nITSandTPCtracksAndSPDtracklets |
e8fd19de | 1360 | } |
1361 | else{ | |
a786bd57 | 1362 | //17 |
1363 | //fCFContCascadeCuts->SetBinLimits(17, 0.0, 250.0 ); // nTrackPrimaryMultiplicity | |
1364 | Double_t *lBinLim17 = new Double_t[ lNbBinsPerVar[17]+1 ]; | |
1365 | lBinLim17[0] = 0; lBinLim17[10] = 10; lBinLim17[20] = 24; lBinLim17[30] = 45; lBinLim17[40] = 95; lBinLim17[50] = 250; | |
1366 | lBinLim17[1] = 1; lBinLim17[11] = 11; lBinLim17[21] = 25; lBinLim17[31] = 50; lBinLim17[41] = 100; | |
1367 | lBinLim17[2] = 2; lBinLim17[12] = 13; lBinLim17[22] = 27; lBinLim17[32] = 55; lBinLim17[42] = 105; | |
1368 | lBinLim17[3] = 3; lBinLim17[13] = 14; lBinLim17[23] = 30; lBinLim17[33] = 60; lBinLim17[43] = 110; | |
1369 | lBinLim17[4] = 4; lBinLim17[14] = 15; lBinLim17[24] = 31; lBinLim17[34] = 65; lBinLim17[44] = 115; | |
1370 | lBinLim17[5] = 5; lBinLim17[15] = 16; lBinLim17[25] = 32; lBinLim17[35] = 70; lBinLim17[45] = 120; | |
1371 | lBinLim17[6] = 6; lBinLim17[16] = 20; lBinLim17[26] = 33; lBinLim17[36] = 75; lBinLim17[46] = 125; | |
1372 | lBinLim17[7] = 7; lBinLim17[17] = 21; lBinLim17[27] = 34; lBinLim17[37] = 80; lBinLim17[47] = 130; | |
1373 | lBinLim17[8] = 8; lBinLim17[18] = 22; lBinLim17[28] = 35; lBinLim17[38] = 85; lBinLim17[48] = 135; | |
1374 | lBinLim17[9] = 9; lBinLim17[19] = 23; lBinLim17[29] = 40; lBinLim17[39] = 90; lBinLim17[49] = 140; | |
1375 | ||
1376 | fCFContCascadeCuts -> SetBinLimits(17, lBinLim17 ); // nTrackPrimaryMultiplicity : 0 to 250 | |
1377 | delete [] lBinLim17; | |
1378 | ||
1379 | fCFContCascadeCuts->SetBinLimits(18, 0.0, 200.0 ); // nITSandTPCtracksAndSPDtracklets | |
e8fd19de | 1380 | } |
1381 | fCFContCascadeCuts->SetBinLimits(19, 68.0 ,164.0 ); // BachTPCClusters | |
1382 | ||
1383 | ||
1384 | // Regular binning definition (valid for v4-18-10-AN on) | |
1385 | /* | |
1386 | fCFContCascadeCuts->SetBinLimits(0, 0.0 , 2.5 ); // DcaXiDaughters | |
32599cd4 | 1387 | fCFContCascadeCuts->SetBinLimits(1, 0.0 , 0.25 ); // DcaBachToPrimVertexXi |
e8fd19de | 1388 | fCFContCascadeCuts->SetBinLimits(2, 0.99 , 1.0 ); // XiCosineOfPointingAngle |
32599cd4 | 1389 | fCFContCascadeCuts->SetBinLimits(3, 0.0 , 4.0 ); // XiRadius |
9d670198 | 1390 | fCFContCascadeCuts->SetBinLimits(4, 1.1 , 1.15 ); // InvMassLambdaAsCascDghter |
e8fd19de | 1391 | fCFContCascadeCuts->SetBinLimits(5, 0.0 , 1.0 ); // DcaV0DaughtersXi |
32599cd4 | 1392 | fCFContCascadeCuts->SetBinLimits(6, 0.98 , 1.0 ); // V0CosineOfPointingAngleXi |
1393 | fCFContCascadeCuts->SetBinLimits(7, 0.0 , 20.0 ); // V0RadiusXi | |
1394 | fCFContCascadeCuts->SetBinLimits(8, 0.0 , 1.0 ); // DcaV0ToPrimVertexXi | |
e8fd19de | 1395 | fCFContCascadeCuts->SetBinLimits(9, 0.0 , 0.25 ); // DcaPosToPrimVertexXi |
1396 | fCFContCascadeCuts->SetBinLimits(10, 0.0 , 0.25 ); // DcaNegToPrimVertexXi | |
1397 | fCFContCascadeCuts->SetBinLimits(11, 1.25 , 1.40 ); // InvMassXi | |
1398 | fCFContCascadeCuts->SetBinLimits(12, 1.62 , 1.74 ); // InvMassOmega | |
9d670198 | 1399 | fCFContCascadeCuts->SetBinLimits(13, 0.0 , 10.0 ); // XiTransvMom |
1400 | fCFContCascadeCuts->SetBinLimits(14, -1.1 , 1.1 ); // Y(Xi) | |
1401 | fCFContCascadeCuts->SetBinLimits(15, -1.1 , 1.1 ); // Y(Omega) | |
1402 | fCFContCascadeCuts->SetBinLimits(16, -10.0 , 10.0 ); // BestPrimaryVtxPosZ | |
32599cd4 | 1403 | if(fCollidingSystems){ |
c525f360 | 1404 | fCFContCascadeCuts->SetBinLimits(17, 0.0, 10000.0 ); // nTrackPrimaryMultiplicity |
a786bd57 | 1405 | fCFContCascadeCuts->SetBinLimits(18, 0.0, 10000.0 ); // nITSandTPCtracksAndSPDtracklets |
32599cd4 | 1406 | } |
1407 | else{ | |
c525f360 | 1408 | fCFContCascadeCuts->SetBinLimits(17, 0.0, 250.0 ); // nTrackPrimaryMultiplicity |
a786bd57 | 1409 | fCFContCascadeCuts->SetBinLimits(18, 0.0, 200.0 ); // nITSandTPCtracksAndSPDtracklets |
32599cd4 | 1410 | } |
9d670198 | 1411 | fCFContCascadeCuts->SetBinLimits(19, 25.0 ,165.0 ); // BachTPCClusters |
e8fd19de | 1412 | */ |
32599cd4 | 1413 | |
1414 | ||
9d670198 | 1415 | // Setting the number of steps : one for each cascade species (Xi-, Xi+ and Omega-, Omega+) |
1416 | fCFContCascadeCuts->SetStepTitle(0, "#Xi^{-} candidates"); | |
1417 | fCFContCascadeCuts->SetStepTitle(1, "#bar{#Xi}^{+} candidates"); | |
1418 | fCFContCascadeCuts->SetStepTitle(2, "#Omega^{-} candidates"); | |
1419 | fCFContCascadeCuts->SetStepTitle(3, "#bar{#Omega}^{+} candidates"); | |
32599cd4 | 1420 | |
1421 | // Setting the variable title, per axis | |
9d670198 | 1422 | // fCFContCascadeCuts->SetVarTitle(40, "Chi2Xi"); |
1423 | fCFContCascadeCuts->SetVarTitle(0, "Dca(XiDaughters) (cm)"); | |
1424 | fCFContCascadeCuts->SetVarTitle(1, "Dca(Bach/PrimVertex) (cm)"); | |
1425 | fCFContCascadeCuts->SetVarTitle(2, "cos(Xi pointing angle)"); | |
1426 | fCFContCascadeCuts->SetVarTitle(3, "R_{2d}(Xi decay) (cm)"); | |
1427 | fCFContCascadeCuts->SetVarTitle(4, "M_{#Lambda}(As Casc Dghter) (GeV/c^{2})"); | |
1428 | // fCFContCascadeCuts->SetVarTitle(40, "V0Chi2Xi"); | |
e8fd19de | 1429 | fCFContCascadeCuts->SetVarTitle(5, "Dca(V0 Daughters) in Xi (cm)"); |
9d670198 | 1430 | |
1431 | fCFContCascadeCuts->SetVarTitle(6, "cos(V0 pointing Angle) in Casc"); | |
1432 | fCFContCascadeCuts->SetVarTitle(7, "R_{2d}(V0 decay) (cm)"); | |
1433 | fCFContCascadeCuts->SetVarTitle(8, "Dca(V0/PrimVertex) (cm)"); | |
1434 | fCFContCascadeCuts->SetVarTitle(9, "Dca(Pos/PrimVertex) (cm)"); | |
1435 | fCFContCascadeCuts->SetVarTitle(10, "Dca(Neg/PrimVertex) (cm)"); | |
1436 | ||
1437 | fCFContCascadeCuts->SetVarTitle(11, "Inv. Mass(Xi) (GeV/c^{2})"); | |
1438 | fCFContCascadeCuts->SetVarTitle(12, "Inv. Mass(Omega) (GeV/c^{2})"); | |
1439 | ||
1440 | fCFContCascadeCuts->SetVarTitle(13, "pt(Casc.) (GeV/c)"); | |
1441 | //fCFContCascadeCuts->SetVarTitle(40, "V0toXiCosineOfPointingAngle"); | |
1442 | ||
1443 | fCFContCascadeCuts->SetVarTitle(14, "Y(Xi)"); | |
1444 | fCFContCascadeCuts->SetVarTitle(15, "Y(Omega)"); | |
1445 | ||
1446 | fCFContCascadeCuts->SetVarTitle(16, "Z-position(BestPrimVtx) (cm)"); | |
1447 | ||
c525f360 | 1448 | fCFContCascadeCuts->SetVarTitle(17, "Primary Track Multiplicity"); |
a786bd57 | 1449 | fCFContCascadeCuts->SetVarTitle(18, "(ITS+TPC tracks + SPD tracklets) Multiplicity"); |
9d670198 | 1450 | fCFContCascadeCuts->SetVarTitle(19, "Bach.TPC Clusters"); |
32599cd4 | 1451 | |
1452 | fListHistCascade->Add(fCFContCascadeCuts); | |
1453 | } | |
1454 | ||
3873f6ee | 1455 | |
3aedd4a5 | 1456 | |
32599cd4 | 1457 | // Part 4 : Angular correlation study ------- |
3873f6ee | 1458 | |
1459 | if(! fHnSpAngularCorrXiMinus){ | |
1460 | // Delta Phi(Casc,any trck) Vs Delta Eta(Casc,any trck) Vs Casc Pt Vs Pt of the tracks | |
1461 | // Delta Phi = 360 bins de -180., 180. | |
1462 | // Delta Eta = 120 bins de -3.0, 3.0 | |
1463 | // Pt Cascade = 100 bins de 0., 10.0, | |
1464 | // Pt track = 150 bins de 0., 15.0 | |
1465 | ||
1466 | Int_t bins[5] = { 360, 120, 100, 150, 40}; | |
3aedd4a5 | 1467 | Double_t xmin[5] = {-50., -3., 0., 0., 1.30}; |
1468 | Double_t xmax[5] = { 310., 3., 10., 15., 1.34}; | |
c525f360 | 1469 | |
1470 | TString strHnSparseTitle(""); | |
1471 | TString strAxisTitle[5]; | |
1472 | if(fAngularCorrelationType == "TrigAnyCasc-AssoAnyPrim" ){ | |
1473 | strHnSparseTitle = "Angular Correlation for #Xi^{-}: Trig = Casc. / Asso = all prim. tracks"; | |
1474 | strAxisTitle[0] = " #Delta#phi(Casc,Track) (deg)"; | |
1475 | strAxisTitle[1] = " #Delta#eta(Casc,Track)"; | |
1476 | strAxisTitle[2] = " Pt_{Casc} (GeV/c)"; | |
1477 | strAxisTitle[3] = " Pt_{asso. track} (GeV/c)"; | |
1478 | } | |
1479 | else if(fAngularCorrelationType == "TrigCascLeading-AssoAnyPrim"){ | |
1480 | strHnSparseTitle = "Angular Correlation for #Xi^{-}: Trig = Casc. (leading part.) / Asso = all prim. tracks"; | |
1481 | strAxisTitle[0] = " #Delta#phi(Casc_{LEADING},Track) (deg)"; | |
1482 | strAxisTitle[1] = " #Delta#eta(Casc_{LEADING},Track)"; | |
1483 | strAxisTitle[2] = " Pt(Casc_{LEADING}) (GeV/c)"; | |
1484 | strAxisTitle[3] = " Pt_{asso. track} (GeV/c)"; | |
1485 | } | |
1486 | else if(fAngularCorrelationType == "TrigLeadingTrck-AssoCasc"){ | |
1487 | strHnSparseTitle = "Angular Correlation for #Xi^{-}: Trig = leading track / Asso = any cascade"; | |
1488 | strAxisTitle[0] = " #Delta#phi(Leading Track,Casc) (deg)"; | |
1489 | strAxisTitle[1] = " #Delta#eta(Leading Track,Casc)"; | |
1490 | strAxisTitle[2] = " Pt(asso. Casc) (GeV/c)"; | |
1491 | strAxisTitle[3] = " Pt_{Leading track} (GeV/c)"; | |
1492 | ||
1493 | } | |
1494 | strAxisTitle[4] = " Eff. Inv Mass (GeV/c^{2})"; | |
1495 | ||
1496 | fHnSpAngularCorrXiMinus = new THnSparseF("fHnSpAngularCorrXiMinus", strHnSparseTitle.Data(), 5, bins, xmin, xmax); | |
1497 | fHnSpAngularCorrXiMinus->GetAxis(0)->SetTitle( strAxisTitle[0].Data() ); | |
1498 | fHnSpAngularCorrXiMinus->GetAxis(1)->SetTitle( strAxisTitle[1].Data() ); | |
1499 | fHnSpAngularCorrXiMinus->GetAxis(2)->SetTitle( strAxisTitle[2].Data() ); | |
1500 | fHnSpAngularCorrXiMinus->GetAxis(3)->SetTitle( strAxisTitle[3].Data() ); | |
1501 | fHnSpAngularCorrXiMinus->GetAxis(4)->SetTitle( strAxisTitle[4].Data() ); | |
1502 | fHnSpAngularCorrXiMinus->Sumw2(); | |
3873f6ee | 1503 | fListHistCascade->Add(fHnSpAngularCorrXiMinus); |
1504 | } | |
1505 | ||
1506 | if(! fHnSpAngularCorrXiPlus){ | |
1507 | // Delta Phi(Casc,any trck) Vs Delta Eta(Casc,any trck) Vs Casc Pt Vs Pt of the tracks | |
1508 | // Delta Phi = 360 bins de -180., 180. | |
1509 | // Delta Eta = 120 bins de -3.0, 3.0 | |
1510 | // Pt Cascade = 100 bins de 0., 10.0, | |
1511 | // Pt track = 150 bins de 0., 15.0 | |
1512 | Int_t bins[5] = { 360, 120, 100, 150, 40}; | |
3aedd4a5 | 1513 | Double_t xmin[5] = {-50., -3., 0., 0., 1.30}; |
1514 | Double_t xmax[5] = { 310., 3., 10., 15., 1.34}; | |
c525f360 | 1515 | |
1516 | TString strHnSparseTitle(""); | |
1517 | TString strAxisTitle[5]; | |
1518 | if(fAngularCorrelationType == "TrigAnyCasc-AssoAnyPrim" ){ | |
1519 | strHnSparseTitle = "Angular Correlation for #bar{#Xi}^{+}: Trig = Casc. / Asso = all prim. tracks"; | |
1520 | strAxisTitle[0] = " #Delta#phi(Casc,Track) (deg)"; | |
1521 | strAxisTitle[1] = " #Delta#eta(Casc,Track)"; | |
1522 | strAxisTitle[2] = " Pt_{Casc} (GeV/c)"; | |
1523 | strAxisTitle[3] = " Pt_{asso. track} (GeV/c)"; | |
1524 | } | |
1525 | else if(fAngularCorrelationType == "TrigCascLeading-AssoAnyPrim"){ | |
1526 | strHnSparseTitle = "Angular Correlation for #bar{#Xi}^{+}: Trig = Casc. (leading part.) / Asso = all prim. tracks"; | |
1527 | strAxisTitle[0] = " #Delta#phi(Casc_{LEADING},Track) (deg)"; | |
1528 | strAxisTitle[1] = " #Delta#eta(Casc_{LEADING},Track)"; | |
1529 | strAxisTitle[2] = " Pt(Casc_{LEADING}) (GeV/c)"; | |
1530 | strAxisTitle[3] = " Pt_{asso. track} (GeV/c)"; | |
1531 | } | |
1532 | else if(fAngularCorrelationType == "TrigLeadingTrck-AssoCasc"){ | |
1533 | strHnSparseTitle = "Angular Correlation for #bar{#Xi}^{+}: Trig = leading track / Asso = any cascade"; | |
1534 | strAxisTitle[0] = " #Delta#phi(Leading Track,Casc) (deg)"; | |
1535 | strAxisTitle[1] = " #Delta#eta(Leading Track,Casc)"; | |
1536 | strAxisTitle[2] = " Pt(asso. Casc) (GeV/c)"; | |
1537 | strAxisTitle[3] = " Pt_{Leading track} (GeV/c)"; | |
1538 | ||
1539 | } | |
1540 | strAxisTitle[4] = " Eff. Inv Mass (GeV/c^{2})"; | |
1541 | ||
1542 | fHnSpAngularCorrXiPlus = new THnSparseF("fHnSpAngularCorrXiPlus", strHnSparseTitle.Data(), 5, bins, xmin, xmax); | |
1543 | fHnSpAngularCorrXiPlus->GetAxis(0)->SetTitle( strAxisTitle[0].Data() ); | |
1544 | fHnSpAngularCorrXiPlus->GetAxis(1)->SetTitle( strAxisTitle[1].Data() ); | |
1545 | fHnSpAngularCorrXiPlus->GetAxis(2)->SetTitle( strAxisTitle[2].Data() ); | |
1546 | fHnSpAngularCorrXiPlus->GetAxis(3)->SetTitle( strAxisTitle[3].Data() ); | |
1547 | fHnSpAngularCorrXiPlus->GetAxis(4)->SetTitle( strAxisTitle[4].Data() ); | |
3873f6ee | 1548 | fHnSpAngularCorrXiPlus->Sumw2(); |
1549 | fListHistCascade->Add(fHnSpAngularCorrXiPlus); | |
1550 | } | |
1551 | ||
1552 | if(! fHnSpAngularCorrOmegaMinus){ | |
1553 | // Delta Phi(Casc,any trck) Vs Delta Eta(Casc,any trck) Vs Casc Pt Vs Pt of the tracks | |
1554 | // Delta Phi = 360 bins de -180., 180. | |
1555 | // Delta Eta = 120 bins de -3.0, 3.0 | |
1556 | // Pt Cascade = 100 bins de 0., 10.0, | |
1557 | // Pt track = 150 bins de 0., 15.0 | |
1558 | ||
1559 | Int_t bins[5] = { 360, 120, 100, 150, 40}; | |
3aedd4a5 | 1560 | Double_t xmin[5] = {-50., -3., 0., 0., 1.65}; |
1561 | Double_t xmax[5] = { 310., 3., 10., 15., 1.69}; | |
c525f360 | 1562 | |
1563 | TString strHnSparseTitle(""); | |
1564 | TString strAxisTitle[5]; | |
1565 | if(fAngularCorrelationType == "TrigAnyCasc-AssoAnyPrim" ){ | |
1566 | strHnSparseTitle = "Angular Correlation for #Omega^{-}: Trig = Casc. / Asso = all prim. tracks"; | |
1567 | strAxisTitle[0] = " #Delta#phi(Casc,Track) (deg)"; | |
1568 | strAxisTitle[1] = " #Delta#eta(Casc,Track)"; | |
1569 | strAxisTitle[2] = " Pt_{Casc} (GeV/c)"; | |
1570 | strAxisTitle[3] = " Pt_{asso. track} (GeV/c)"; | |
1571 | } | |
1572 | else if(fAngularCorrelationType == "TrigCascLeading-AssoAnyPrim"){ | |
1573 | strHnSparseTitle = "Angular Correlation for #Omega^{-}: Trig = Casc. (leading part.) / Asso = all prim. tracks"; | |
1574 | strAxisTitle[0] = " #Delta#phi(Casc_{LEADING},Track) (deg)"; | |
1575 | strAxisTitle[1] = " #Delta#eta(Casc_{LEADING},Track)"; | |
1576 | strAxisTitle[2] = " Pt(Casc_{LEADING}) (GeV/c)"; | |
1577 | strAxisTitle[3] = " Pt_{asso. track} (GeV/c)"; | |
1578 | } | |
1579 | else if(fAngularCorrelationType == "TrigLeadingTrck-AssoCasc"){ | |
1580 | strHnSparseTitle = "Angular Correlation for #Omega^{-}: Trig = leading track / Asso = any cascade"; | |
1581 | strAxisTitle[0] = " #Delta#phi(Leading Track,Casc) (deg)"; | |
1582 | strAxisTitle[1] = " #Delta#eta(Leading Track,Casc)"; | |
1583 | strAxisTitle[2] = " Pt(asso. Casc) (GeV/c)"; | |
1584 | strAxisTitle[3] = " Pt_{Leading track} (GeV/c)"; | |
1585 | ||
1586 | } | |
1587 | strAxisTitle[4] = " Eff. Inv Mass (GeV/c^{2})"; | |
1588 | ||
1589 | fHnSpAngularCorrOmegaMinus = new THnSparseF("fHnSpAngularCorrOmegaMinus", strHnSparseTitle.Data(), 5, bins, xmin, xmax); | |
1590 | fHnSpAngularCorrOmegaMinus->GetAxis(0)->SetTitle( strAxisTitle[0].Data() ); | |
1591 | fHnSpAngularCorrOmegaMinus->GetAxis(1)->SetTitle( strAxisTitle[1].Data() ); | |
1592 | fHnSpAngularCorrOmegaMinus->GetAxis(2)->SetTitle( strAxisTitle[2].Data() ); | |
1593 | fHnSpAngularCorrOmegaMinus->GetAxis(3)->SetTitle( strAxisTitle[3].Data() ); | |
1594 | fHnSpAngularCorrOmegaMinus->GetAxis(4)->SetTitle( strAxisTitle[4].Data() ); | |
3873f6ee | 1595 | fHnSpAngularCorrOmegaMinus->Sumw2(); |
1596 | fListHistCascade->Add(fHnSpAngularCorrOmegaMinus); | |
1597 | } | |
1598 | ||
1599 | if(! fHnSpAngularCorrOmegaPlus){ | |
1600 | // Delta Phi(Casc,any trck) Vs Delta Eta(Casc,any trck) Vs Casc Pt Vs Pt of the tracks | |
1601 | // Delta Phi = 360 bins de -180., 180. | |
1602 | // Delta Eta = 120 bins de -3.0, 3.0 | |
1603 | // Pt Cascade = 100 bins de 0., 10.0, | |
1604 | // Pt track = 150 bins de 0., 15.0 | |
1605 | Int_t bins[5] = { 360, 120, 100, 150, 40}; | |
3aedd4a5 | 1606 | Double_t xmin[5] = {-50., -3., 0., 0., 1.65}; |
1607 | Double_t xmax[5] = { 310., 3., 10., 15., 1.69}; | |
c525f360 | 1608 | |
1609 | TString strHnSparseTitle(""); | |
1610 | TString strAxisTitle[5]; | |
1611 | if(fAngularCorrelationType == "TrigAnyCasc-AssoAnyPrim" ){ | |
1612 | strHnSparseTitle = "Angular Correlation for #bar{#Omega}^{+}: Trig = Casc. / Asso = all prim. tracks"; | |
1613 | strAxisTitle[0] = " #Delta#phi(Casc,Track) (deg)"; | |
1614 | strAxisTitle[1] = " #Delta#eta(Casc,Track)"; | |
1615 | strAxisTitle[2] = " Pt_{Casc} (GeV/c)"; | |
1616 | strAxisTitle[3] = " Pt_{asso. track} (GeV/c)"; | |
1617 | } | |
1618 | else if(fAngularCorrelationType == "TrigCascLeading-AssoAnyPrim"){ | |
1619 | strHnSparseTitle = "Angular Correlation for #bar{#Omega}^{+}: Trig = Casc. (leading part.) / Asso = all prim. tracks"; | |
1620 | strAxisTitle[0] = " #Delta#phi(Casc_{LEADING},Track) (deg)"; | |
1621 | strAxisTitle[1] = " #Delta#eta(Casc_{LEADING},Track)"; | |
1622 | strAxisTitle[2] = " Pt(Casc_{LEADING}) (GeV/c)"; | |
1623 | strAxisTitle[3] = " Pt_{asso. track} (GeV/c)"; | |
1624 | } | |
1625 | else if(fAngularCorrelationType == "TrigLeadingTrck-AssoCasc"){ | |
1626 | strHnSparseTitle = "Angular Correlation for #bar{#Omega}^{+}: Trig = leading track / Asso = any cascade"; | |
1627 | strAxisTitle[0] = " #Delta#phi(Leading Track,Casc) (deg)"; | |
1628 | strAxisTitle[1] = " #Delta#eta(Leading Track,Casc)"; | |
1629 | strAxisTitle[2] = " Pt(asso. Casc) (GeV/c)"; | |
1630 | strAxisTitle[3] = " Pt_{Leading track} (GeV/c)"; | |
1631 | ||
1632 | } | |
1633 | strAxisTitle[4] = " Eff. Inv Mass (GeV/c^{2})"; | |
1634 | ||
1635 | fHnSpAngularCorrOmegaPlus = new THnSparseF("fHnSpAngularCorrOmegaPlus", strHnSparseTitle.Data(), 5, bins, xmin, xmax); | |
1636 | fHnSpAngularCorrOmegaPlus->GetAxis(0)->SetTitle( strAxisTitle[0].Data() ); | |
1637 | fHnSpAngularCorrOmegaPlus->GetAxis(1)->SetTitle( strAxisTitle[1].Data() ); | |
1638 | fHnSpAngularCorrOmegaPlus->GetAxis(2)->SetTitle( strAxisTitle[2].Data() ); | |
1639 | fHnSpAngularCorrOmegaPlus->GetAxis(3)->SetTitle( strAxisTitle[3].Data() ); | |
1640 | fHnSpAngularCorrOmegaPlus->GetAxis(4)->SetTitle( strAxisTitle[4].Data() ); | |
3873f6ee | 1641 | fHnSpAngularCorrOmegaPlus->Sumw2(); |
1642 | fListHistCascade->Add(fHnSpAngularCorrOmegaPlus); | |
1643 | } | |
1644 | ||
1645 | ||
9d670198 | 1646 | PostData(1, fListHistCascade); |
1647 | /* PostData(2, fPaveTextBookKeeping);*/ | |
45ee0bcc | 1648 | }// end UserCreateOutputObjects |
1649 | ||
1650 | ||
1651 | ||
1652 | ||
1653 | ||
1654 | ||
1655 | //________________________________________________________________________ | |
1656 | void AliAnalysisTaskCheckCascade::UserExec(Option_t *) | |
1657 | { | |
1658 | // Main loop | |
1659 | // Called for each event | |
7e504402 | 1660 | |
1661 | AliESDEvent *lESDevent = 0x0; | |
1662 | AliAODEvent *lAODevent = 0x0; | |
9d670198 | 1663 | Int_t ncascades = -1; |
1664 | Int_t nTrackMultiplicity = -1; | |
c525f360 | 1665 | Int_t nTrackWithTPCrefitMultiplicity = -1; |
1666 | Int_t nTrackPrimaryMultiplicity = -1; | |
a786bd57 | 1667 | Int_t nSPDTracklets = 0; // AliESDEvent::EstimateMultiplicity will re-initialise the value to 0 |
1668 | Int_t nITSandTPCtracksAndSPDtracklets = 0; // AliESDEvent::EstimateMultiplicity will re-initialise the value to 0 | |
1669 | Int_t nTracksITSSApure = 0; // AliESDEvent::EstimateMultiplicity will re-initialise the value to 0 | |
9d670198 | 1670 | |
1671 | Short_t lStatusTrackingPrimVtx = -2; | |
1672 | Double_t lTrkgPrimaryVtxPos[3] = {-100.0, -100.0, -100.0}; | |
1673 | Double_t lBestPrimaryVtxPos[3] = {-100.0, -100.0, -100.0}; | |
1674 | Double_t lMagneticField = -10.; | |
1675 | ||
45ee0bcc | 1676 | |
1677 | ||
1678 | // Connect to the InputEvent | |
1679 | // After these lines, we should have an ESD/AOD event + the number of cascades in it. | |
1680 | ||
1681 | if(fAnalysisType == "ESD"){ | |
1682 | lESDevent = dynamic_cast<AliESDEvent*>( InputEvent() ); | |
1683 | if (!lESDevent) { | |
7e504402 | 1684 | AliWarning("ERROR: lESDevent not available \n"); |
45ee0bcc | 1685 | return; |
1686 | } | |
9d670198 | 1687 | |
1688 | fHistCascadeMultiplicityBeforeTrigSel->Fill ( lESDevent->GetNumberOfCascades() ); | |
1689 | ||
7e504402 | 1690 | //------------------------------------------------- |
1691 | // 0 - Trigger managment | |
9d670198 | 1692 | // NOTE : Check the availability of the proper trigger |
7e504402 | 1693 | |
1694 | // 1st option | |
e8fd19de | 1695 | //AliMCEventHandler *lmcEvtHandler = dynamic_cast<AliMCEventHandler*>( (AliAnalysisManager::GetAnalysisManager())->GetMCtruthEventHandler() ); |
1696 | //if( !lmcEvtHandler ){ // !0x0 = real data or !1 = there is an MC handler available (useMC = kTRUE in AnalysisTrainNew), so = data from MC | |
7e504402 | 1697 | // if ( !( lESDevent->IsTriggerClassFired("CINT1B-ABCE-NOPF-ALL")) ) return; |
1698 | //} | |
1699 | ||
e8fd19de | 1700 | // 2nd option - Presuppose the presence of AliPhysicsSelectionTask FIXME |
9d670198 | 1701 | |
e8fd19de | 1702 | UInt_t maskIsSelected = ((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected(); |
1703 | Bool_t isSelected = 0; | |
c525f360 | 1704 | if( fTriggerMaskType == "kMB") isSelected = (maskIsSelected & AliVEvent::kMB) == AliVEvent::kMB; |
1705 | else if(fTriggerMaskType == "kHighMult") isSelected = (maskIsSelected & AliVEvent::kHighMult) == AliVEvent::kHighMult; | |
1706 | else isSelected = 1; // default = select anyway (use case = run without Phys Selection task) | |
1707 | ||
9d670198 | 1708 | if ( ! isSelected ) { |
1709 | PostData(1, fListHistCascade); | |
1710 | return; | |
1711 | } | |
1712 | ||
7e504402 | 1713 | //else Printf("Event selected ... \n"); |
1714 | ||
1715 | ||
1716 | ||
1717 | //------------------------------------------------- | |
1718 | // 1 - Cascade vertexer (ESD) | |
9d670198 | 1719 | |
1720 | if(fkRerunV0CascVertexers){ // FIXME : relaunch V0 and Cascade vertexers | |
c525f360 | 1721 | if(fAnalysisType == "ESD" ){ |
9d670198 | 1722 | // lESDevent->ResetCascades(); |
1723 | // lESDevent->ResetV0s(); | |
1724 | // | |
1725 | // AliV0vertexer lV0vtxer; | |
1726 | // AliCascadeVertexer lCascVtxer; | |
1727 | // | |
1728 | // lV0vtxer.SetDefaultCuts(fV0Sels); | |
1729 | // lCascVtxer.SetDefaultCuts(fCascSels); | |
1730 | // | |
1731 | // lV0vtxer.Tracks2V0vertices(lESDevent); | |
1732 | // lCascVtxer.V0sTracks2CascadeVertices(lESDevent); | |
c525f360 | 1733 | } |
9d670198 | 1734 | }// end if(RelaunchV0CascVertexers) |
1735 | ||
7e504402 | 1736 | //------------------------------------------------- |
1737 | ncascades = lESDevent->GetNumberOfCascades(); | |
1738 | nTrackWithTPCrefitMultiplicity = DoESDTrackWithTPCrefitMultiplicity(lESDevent); | |
c525f360 | 1739 | nTrackPrimaryMultiplicity = fESDtrackCuts->CountAcceptedTracks(lESDevent); |
a786bd57 | 1740 | //EstimateMultiplicity(Int_t &tracklets, Int_t &trITSTPC, Int_t &trITSSApure, Double_t eta, Bool_t useDCAFlag,Bool_t useV0Flag) |
1741 | lESDevent->EstimateMultiplicity( nSPDTracklets, nITSandTPCtracksAndSPDtracklets, nTracksITSSApure, 1.0, kTRUE, kTRUE); | |
7e504402 | 1742 | |
7e504402 | 1743 | }//if (fAnalysisType == "ESD") |
45ee0bcc | 1744 | |
9d670198 | 1745 | |
45ee0bcc | 1746 | if(fAnalysisType == "AOD"){ |
1747 | lAODevent = dynamic_cast<AliAODEvent*>( InputEvent() ); | |
9d670198 | 1748 | if (!lAODevent) { |
7e504402 | 1749 | AliWarning("ERROR: lAODevent not available \n"); |
45ee0bcc | 1750 | return; |
1751 | } | |
9d670198 | 1752 | ncascades = lAODevent->GetNumberOfCascades(); |
1753 | nTrackWithTPCrefitMultiplicity = -1; | |
c525f360 | 1754 | nTrackPrimaryMultiplicity = -1; |
9d670198 | 1755 | |
1756 | fHistCascadeMultiplicityBeforeTrigSel->Fill ( ncascades ); | |
45ee0bcc | 1757 | } |
9d670198 | 1758 | |
1759 | // For AOD or ESD ... | |
32599cd4 | 1760 | nTrackMultiplicity = (InputEvent())->GetNumberOfTracks(); |
9d670198 | 1761 | |
1762 | ||
1763 | //------------------------------------------------- | |
1764 | fHistTrackMultiplicityForTrigEvt ->Fill( nTrackMultiplicity ); | |
1765 | fHistTPCrefitTrackMultiplicityForTrigEvt ->Fill( nTrackWithTPCrefitMultiplicity ); | |
c525f360 | 1766 | fHistPrimaryTrackMultiplicityForTrigEvt ->Fill( nTrackPrimaryMultiplicity ); |
a786bd57 | 1767 | fHistEstimateITSTPCMultiplicityForTrigEvt->Fill( nITSandTPCtracksAndSPDtracklets); |
9d670198 | 1768 | fHistCascadeMultiplicityForTrigEvt ->Fill( ncascades ); |
1769 | ||
1770 | ||
1771 | ||
1772 | ||
45ee0bcc | 1773 | // --------------------------------------------------------------- |
9d670198 | 1774 | // I - Global characteristics of the events + general histos (filled for any selected events and/or for the analysed events) |
1775 | ||
1776 | // - I.Step 1 : Characteristics of the event : prim. Vtx + magnetic field (ESD) | |
1777 | //------------- | |
1778 | ||
1779 | if(fAnalysisType == "ESD"){ | |
1780 | const AliESDVertex *lPrimaryTrackingESDVtx = lESDevent->GetPrimaryVertexTracks(); | |
1781 | // get the vtx stored in ESD found with tracks | |
1782 | lPrimaryTrackingESDVtx->GetXYZ( lTrkgPrimaryVtxPos ); | |
1783 | ||
1784 | ||
1785 | const AliESDVertex *lPrimaryBestESDVtx = lESDevent->GetPrimaryVertex(); | |
1786 | // get the best primary vertex available for the event | |
1787 | // As done in AliCascadeVertexer, we keep the one which is the best one available. | |
1788 | // between : Tracking vertex > SPD vertex > TPC vertex > default SPD vertex | |
1789 | // This one will be used for next calculations (DCA essentially) | |
1790 | lPrimaryBestESDVtx->GetXYZ( lBestPrimaryVtxPos ); | |
1791 | lStatusTrackingPrimVtx = lPrimaryTrackingESDVtx->GetStatus(); | |
1792 | ||
1793 | // FIXME : quality cut on the z-position of the prim vertex. | |
1794 | if(fkQualityCutZprimVtxPos) { | |
1795 | if(TMath::Abs(lBestPrimaryVtxPos[2]) > 10.0 ) { | |
1796 | AliWarning("Pb / | Z position of Best Prim Vtx | > 10.0 cm ... return !"); | |
1797 | PostData(1, fListHistCascade); | |
1798 | return; | |
1799 | } | |
1800 | } | |
e8fd19de | 1801 | |
1802 | fHistCascadeMultiplicityForTrigEvtAndZprimVtx->Fill( ncascades ); | |
1803 | ||
c525f360 | 1804 | // FIXME : quality selection regarding pile-up rejection |
1805 | if(fkRejectEventPileUp) { | |
a786bd57 | 1806 | if(lESDevent->IsPileupFromSPDInMultBins() ){// minContributors=3, minZdist=0.8, nSigmaZdist=3., nSigmaDiamXY=2., nSigmaDiamZ=5. -> see http://alisoft.cern.ch/viewvc/trunk/STEER/AliESDEvent.h?root=AliRoot&r1=41914&r2=42199&pathrev=42199 |
c525f360 | 1807 | AliWarning("Pb / Event tagged as pile-up by SPD... return !"); |
1808 | PostData(1, fListHistCascade); | |
1809 | return; | |
1810 | } | |
1811 | } | |
1812 | ||
1813 | fHistCascadeMultiplicityForTrigEvtNonPiledUpAndZprimVtx->Fill( ncascades ); | |
1814 | ||
9d670198 | 1815 | // FIXME : remove TPC-only primary vertex : retain only events with tracking + SPD vertex |
1816 | if(fkQualityCutNoTPConlyPrimVtx) { | |
1817 | const AliESDVertex *lPrimarySPDVtx = lESDevent->GetPrimaryVertexSPD(); | |
1818 | if (!lPrimarySPDVtx->GetStatus() && !lPrimaryTrackingESDVtx->GetStatus() ){ | |
1819 | AliWarning("Pb / No SPD prim. vertex nor prim. Tracking vertex ... return !"); | |
1820 | PostData(1, fListHistCascade); | |
1821 | return; | |
1822 | } | |
1823 | } | |
1824 | ||
1825 | // NOTE : For older evts | |
1826 | ||
1827 | // As previously done in AliCascadeVertexer, we keep, between both retrieved vertices (SPD or Tracking) | |
1828 | // the one which is the best one available. | |
1829 | // This one will be used for next calculations (DCA essentially) | |
1830 | // At that time, the TPC-only primary vertex was not considered | |
1831 | ||
1832 | ||
1833 | lMagneticField = lESDevent->GetMagneticField( ); | |
1834 | // FIXME if(TMath::Abs(lMagneticField ) < 10e-6) continue; | |
1835 | ||
1836 | }// end if(ESD) | |
1837 | ||
1838 | if(fAnalysisType == "AOD"){ | |
1839 | // To be developed | |
1840 | const AliAODVertex *lPrimaryBestAODVtx = lAODevent->GetPrimaryVertex(); | |
1841 | // get the best primary vertex available for the event | |
1842 | // We may keep the one which is the best one available = GetVertex(0) | |
1843 | // Pb with pile-up to expect | |
1844 | // This one will be used for next calculations (DCA essentially) | |
1845 | lPrimaryBestAODVtx->GetXYZ( lBestPrimaryVtxPos ); | |
1846 | ||
1847 | lStatusTrackingPrimVtx = -1; | |
1848 | lTrkgPrimaryVtxPos[0] = -100.0; | |
1849 | lTrkgPrimaryVtxPos[1] = -100.0; | |
1850 | lTrkgPrimaryVtxPos[2] = -100.0; | |
1851 | lMagneticField = 0.; | |
1852 | } | |
45ee0bcc | 1853 | |
9d670198 | 1854 | |
1855 | // - I.Step 2 : Filling histos that characterize the selected event : x,y,z prim. Vtx distrib. (ESD) | |
1856 | //------------- | |
1857 | ||
1858 | fHistCascadeMultiplicityForSelEvt ->Fill( ncascades ); | |
1859 | fHistPosBestPrimaryVtxXForSelEvt ->Fill( lBestPrimaryVtxPos[0] ); | |
1860 | fHistPosBestPrimaryVtxYForSelEvt ->Fill( lBestPrimaryVtxPos[1] ); | |
1861 | fHistPosBestPrimaryVtxZForSelEvt ->Fill( lBestPrimaryVtxPos[2] ); | |
1862 | ||
7e504402 | 1863 | |
45ee0bcc | 1864 | |
32599cd4 | 1865 | // --------------------------------------------------------------- |
1866 | // II - Calcultaion Part dedicated to Xi vertices | |
1867 | ||
45ee0bcc | 1868 | for (Int_t iXi = 0; iXi < ncascades; iXi++) |
1869 | {// This is the begining of the Cascade loop (ESD or AOD) | |
1870 | ||
32599cd4 | 1871 | // ------------------------------------- |
1872 | // II.Init - Initialisation of the local variables that will be needed for ESD/AOD | |
1873 | ||
1874 | ||
7e504402 | 1875 | // - 0th part of initialisation : around primary vertex ... |
32599cd4 | 1876 | |
9d670198 | 1877 | Double_t lTrkgPrimaryVtxRadius3D = -500.0; |
32599cd4 | 1878 | Double_t lBestPrimaryVtxRadius3D = -500.0; |
1879 | ||
7e504402 | 1880 | // - 1st part of initialisation : variables needed to store AliESDCascade data members |
32599cd4 | 1881 | Double_t lEffMassXi = 0. ; |
9d670198 | 1882 | Double_t lChi2Xi = -1. ; |
1883 | Double_t lDcaXiDaughters = -1. ; | |
1884 | Double_t lXiCosineOfPointingAngle = -1. ; | |
32599cd4 | 1885 | Double_t lPosXi[3] = { -1000.0, -1000.0, -1000.0 }; |
c525f360 | 1886 | Double_t lXiRadius2D = -1000. ; |
1887 | Double_t lXiRadius3D = -1000. ; | |
7e504402 | 1888 | |
1889 | // - 2nd part of initialisation : Nbr of clusters within TPC for the 3 daughter cascade tracks | |
1890 | Int_t lPosTPCClusters = -1; // For ESD only ...//FIXME : wait for availability in AOD | |
1891 | Int_t lNegTPCClusters = -1; // For ESD only ... | |
1892 | Int_t lBachTPCClusters = -1; // For ESD only ... | |
e8fd19de | 1893 | |
1894 | Double_t lInnerWallMomCascDghters[3] = {-100., -100., -100.}; | |
1895 | Double_t lTPCSignalCascDghters [3] = {-100., -100., -100.}; | |
1896 | ||
1897 | ||
7e504402 | 1898 | // - 3rd part of initialisation : about V0 part in cascades |
32599cd4 | 1899 | Double_t lInvMassLambdaAsCascDghter = 0.; |
9d670198 | 1900 | Double_t lV0Chi2Xi = -1. ; |
1901 | Double_t lDcaV0DaughtersXi = -1.; | |
32599cd4 | 1902 | |
9d670198 | 1903 | Double_t lDcaBachToPrimVertexXi = -1., lDcaV0ToPrimVertexXi = -1.; |
1904 | Double_t lDcaPosToPrimVertexXi = -1.; | |
1905 | Double_t lDcaNegToPrimVertexXi = -1.; | |
1906 | Double_t lV0CosineOfPointingAngleXi = -1. ; | |
32599cd4 | 1907 | Double_t lPosV0Xi[3] = { -1000. , -1000., -1000. }; // Position of VO coming from cascade |
1908 | Double_t lV0RadiusXi = -1000.0; | |
1909 | Double_t lV0quality = 0.; | |
1910 | ||
1911 | ||
7e504402 | 1912 | // - 4th part of initialisation : Effective masses |
32599cd4 | 1913 | Double_t lInvMassXiMinus = 0.; |
1914 | Double_t lInvMassXiPlus = 0.; | |
1915 | Double_t lInvMassOmegaMinus = 0.; | |
1916 | Double_t lInvMassOmegaPlus = 0.; | |
1917 | ||
7e504402 | 1918 | // - 5th part of initialisation : PID treatment |
32599cd4 | 1919 | Bool_t lIsPosInXiProton = kFALSE; |
1920 | Bool_t lIsPosInXiPion = kFALSE; | |
1921 | Bool_t lIsPosInOmegaProton = kFALSE; | |
1922 | Bool_t lIsPosInOmegaPion = kFALSE; | |
1923 | ||
1924 | Bool_t lIsNegInXiProton = kFALSE; | |
1925 | Bool_t lIsNegInXiPion = kFALSE; | |
1926 | Bool_t lIsNegInOmegaProton = kFALSE; | |
1927 | Bool_t lIsNegInOmegaPion = kFALSE; | |
1928 | ||
1929 | Bool_t lIsBachelorKaon = kFALSE; | |
1930 | Bool_t lIsBachelorPion = kFALSE; | |
1931 | ||
1932 | Bool_t lIsBachelorKaonForTPC = kFALSE; // For ESD only ...//FIXME : wait for availability in AOD | |
1933 | Bool_t lIsBachelorPionForTPC = kFALSE; // For ESD only ... | |
1934 | Bool_t lIsNegPionForTPC = kFALSE; // For ESD only ... | |
1935 | Bool_t lIsPosPionForTPC = kFALSE; // For ESD only ... | |
1936 | Bool_t lIsNegProtonForTPC = kFALSE; // For ESD only ... | |
1937 | Bool_t lIsPosProtonForTPC = kFALSE; // For ESD only ... | |
1938 | ||
7e504402 | 1939 | // - 6th part of initialisation : extra info for QA |
1940 | Double_t lXiMomX = 0. , lXiMomY = 0., lXiMomZ = 0.; | |
32599cd4 | 1941 | Double_t lXiTransvMom = 0. ; |
1942 | Double_t lXiTotMom = 0. ; | |
1943 | ||
1944 | Double_t lBachMomX = 0., lBachMomY = 0., lBachMomZ = 0.; | |
1945 | Double_t lBachTransvMom = 0.; | |
1946 | Double_t lBachTotMom = 0.; | |
1947 | ||
1948 | Short_t lChargeXi = -2; | |
1949 | Double_t lV0toXiCosineOfPointingAngle = 0. ; | |
1950 | ||
1951 | Double_t lRapXi = -20.0, lRapOmega = -20.0, lEta = -20.0, lTheta = 360., lPhi = 720. ; | |
1952 | Double_t lAlphaXi = -200., lPtArmXi = -200.0; | |
1953 | ||
7e504402 | 1954 | // - 7th part of initialisation : variables for the AliCFContainer dedicated to cascade cut optmisiation |
a786bd57 | 1955 | // Int_t lSPDTrackletsMultiplicity = -1; // now provided by AliESDEvent::EstimateMultiplicity ... |
32599cd4 | 1956 | |
7e504402 | 1957 | // - 8th part of initialisation : variables needed for Angular correlations |
32599cd4 | 1958 | TVector3 lTVect3MomXi(0.,0.,0.); |
1959 | Int_t lArrTrackID[3] = {-1, -1, -1}; | |
1960 | ||
45ee0bcc | 1961 | |
1962 | if(fAnalysisType == "ESD"){ | |
1963 | ||
1964 | // ------------------------------------- | |
32599cd4 | 1965 | // II.ESD - Calcultaion Part dedicated to Xi vertices (ESD) |
45ee0bcc | 1966 | |
1967 | AliESDcascade *xi = lESDevent->GetCascade(iXi); | |
1968 | if (!xi) continue; | |
e8fd19de | 1969 | |
9d670198 | 1970 | |
1971 | // - II.Step 1 : around primary vertex | |
1972 | //------------- | |
1973 | lTrkgPrimaryVtxRadius3D = TMath::Sqrt( lTrkgPrimaryVtxPos[0] * lTrkgPrimaryVtxPos[0] + | |
1974 | lTrkgPrimaryVtxPos[1] * lTrkgPrimaryVtxPos[1] + | |
1975 | lTrkgPrimaryVtxPos[2] * lTrkgPrimaryVtxPos[2] ); | |
45ee0bcc | 1976 | |
9d670198 | 1977 | lBestPrimaryVtxRadius3D = TMath::Sqrt( lBestPrimaryVtxPos[0] * lBestPrimaryVtxPos[0] + |
1978 | lBestPrimaryVtxPos[1] * lBestPrimaryVtxPos[1] + | |
1979 | lBestPrimaryVtxPos[2] * lBestPrimaryVtxPos[2] ); | |
45ee0bcc | 1980 | |
1981 | ||
45ee0bcc | 1982 | |
1983 | // - II.Step 2 : Assigning the necessary variables for specific AliESDcascade data members (ESD) | |
1984 | //------------- | |
1985 | lV0quality = 0.; | |
1986 | xi->ChangeMassHypothesis(lV0quality , 3312); // default working hypothesis : cascade = Xi- decay | |
1987 | ||
1988 | lEffMassXi = xi->GetEffMassXi(); | |
1989 | lChi2Xi = xi->GetChi2Xi(); | |
1990 | lDcaXiDaughters = xi->GetDcaXiDaughters(); | |
1991 | lXiCosineOfPointingAngle = xi->GetCascadeCosineOfPointingAngle( lBestPrimaryVtxPos[0], | |
9d670198 | 1992 | lBestPrimaryVtxPos[1], |
1993 | lBestPrimaryVtxPos[2] ); | |
45ee0bcc | 1994 | // Take care : the best available vertex should be used (like in AliCascadeVertexer) |
1995 | ||
1996 | xi->GetXYZcascade( lPosXi[0], lPosXi[1], lPosXi[2] ); | |
c525f360 | 1997 | lXiRadius2D = TMath::Sqrt( lPosXi[0]*lPosXi[0] + lPosXi[1]*lPosXi[1] ); |
1998 | lXiRadius3D = TMath::Sqrt( lPosXi[0]*lPosXi[0] + lPosXi[1]*lPosXi[1] + lPosXi[2]*lPosXi[2]); | |
45ee0bcc | 1999 | |
2000 | ||
2001 | ||
2002 | // - II.Step 3 : around the tracks : Bach + V0 (ESD) | |
2003 | // ~ Necessary variables for ESDcascade data members coming from the ESDv0 part (inheritance) | |
2004 | //------------- | |
2005 | ||
7e504402 | 2006 | UInt_t lIdxPosXi = (UInt_t) TMath::Abs( xi->GetPindex() ); |
2007 | UInt_t lIdxNegXi = (UInt_t) TMath::Abs( xi->GetNindex() ); | |
2008 | UInt_t lBachIdx = (UInt_t) TMath::Abs( xi->GetBindex() ); | |
2009 | // Care track label can be negative in MC production (linked with the track quality) | |
2010 | // However = normally, not the case for track index ... | |
2011 | ||
2012 | // FIXME : rejection of a double use of a daughter track (nothing but just a crosscheck of what is done in the cascade vertexer) | |
2013 | if(lBachIdx == lIdxNegXi) { | |
2014 | AliWarning("Pb / Idx(Bach. track) = Idx(Neg. track) ... continue!"); continue; | |
2015 | } | |
2016 | if(lBachIdx == lIdxPosXi) { | |
2017 | AliWarning("Pb / Idx(Bach. track) = Idx(Pos. track) ... continue!"); continue; | |
2018 | } | |
2019 | ||
45ee0bcc | 2020 | AliESDtrack *pTrackXi = lESDevent->GetTrack( lIdxPosXi ); |
2021 | AliESDtrack *nTrackXi = lESDevent->GetTrack( lIdxNegXi ); | |
2022 | AliESDtrack *bachTrackXi = lESDevent->GetTrack( lBachIdx ); | |
7e504402 | 2023 | if (!pTrackXi || !nTrackXi || !bachTrackXi ) { |
2024 | AliWarning("ERROR: Could not retrieve one of the 3 ESD daughter tracks of the cascade ..."); | |
2025 | continue; | |
2026 | } | |
e8fd19de | 2027 | |
2028 | ||
7e504402 | 2029 | lPosTPCClusters = pTrackXi->GetTPCNcls(); |
2030 | lNegTPCClusters = nTrackXi->GetTPCNcls(); | |
e8fd19de | 2031 | lBachTPCClusters = bachTrackXi->GetTPCNcls(); |
2032 | ||
7e504402 | 2033 | // FIXME : rejection of a poor quality tracks |
9d670198 | 2034 | if(fkQualityCutTPCrefit){ |
7e504402 | 2035 | // 1 - Poor quality related to TPCrefit |
2036 | ULong_t pStatus = pTrackXi->GetStatus(); | |
2037 | ULong_t nStatus = nTrackXi->GetStatus(); | |
2038 | ULong_t bachStatus = bachTrackXi->GetStatus(); | |
2039 | if ((pStatus&AliESDtrack::kTPCrefit) == 0) { AliWarning("Pb / V0 Pos. track has no TPCrefit ... continue!"); continue; } | |
2040 | if ((nStatus&AliESDtrack::kTPCrefit) == 0) { AliWarning("Pb / V0 Neg. track has no TPCrefit ... continue!"); continue; } | |
2041 | if ((bachStatus&AliESDtrack::kTPCrefit) == 0) { AliWarning("Pb / Bach. track has no TPCrefit ... continue!"); continue; } | |
2042 | } | |
9d670198 | 2043 | if(fkQualityCut80TPCcls){ |
7e504402 | 2044 | // 2 - Poor quality related to TPC clusters |
2045 | if(lPosTPCClusters < 80) { AliWarning("Pb / V0 Pos. track has less than 80 TPC clusters ... continue!"); continue; } | |
2046 | if(lNegTPCClusters < 80) { AliWarning("Pb / V0 Neg. track has less than 80 TPC clusters ... continue!"); continue; } | |
2047 | if(lBachTPCClusters < 80) { AliWarning("Pb / Bach. track has less than 80 TPC clusters ... continue!"); continue; } | |
2048 | } | |
e8fd19de | 2049 | |
2050 | const AliExternalTrackParam *pExtTrack = pTrackXi ->GetInnerParam(); | |
2051 | const AliExternalTrackParam *nExtTrack = nTrackXi ->GetInnerParam(); | |
2052 | const AliExternalTrackParam *bachExtTrack = bachTrackXi ->GetInnerParam(); | |
2053 | ||
2054 | if (pExtTrack) { | |
2055 | lInnerWallMomCascDghters[0] = pExtTrack ->GetP() * pExtTrack ->Charge(); | |
2056 | lTPCSignalCascDghters [0] = pTrackXi ->GetTPCsignal(); | |
2057 | } | |
2058 | if (nExtTrack) { | |
2059 | lInnerWallMomCascDghters[1] = nExtTrack ->GetP() * nExtTrack ->Charge(); | |
2060 | lTPCSignalCascDghters [1] = nTrackXi ->GetTPCsignal(); | |
2061 | } | |
2062 | if (bachExtTrack) { | |
2063 | lInnerWallMomCascDghters[2] = bachExtTrack ->GetP() * bachExtTrack ->Charge(); | |
2064 | lTPCSignalCascDghters [2] = bachTrackXi ->GetTPCsignal(); | |
2065 | } | |
2066 | ||
7e504402 | 2067 | |
45ee0bcc | 2068 | lInvMassLambdaAsCascDghter = xi->GetEffMass(); |
2069 | // This value shouldn't change, whatever the working hyp. is : Xi-, Xi+, Omega-, Omega+ | |
2070 | lDcaV0DaughtersXi = xi->GetDcaV0Daughters(); | |
2071 | lV0Chi2Xi = xi->GetChi2V0(); | |
2072 | ||
2073 | lV0CosineOfPointingAngleXi = xi->GetV0CosineOfPointingAngle( lBestPrimaryVtxPos[0], | |
2074 | lBestPrimaryVtxPos[1], | |
2075 | lBestPrimaryVtxPos[2] ); | |
2076 | ||
2077 | lDcaV0ToPrimVertexXi = xi->GetD( lBestPrimaryVtxPos[0], | |
2078 | lBestPrimaryVtxPos[1], | |
2079 | lBestPrimaryVtxPos[2] ); | |
2080 | ||
2081 | lDcaBachToPrimVertexXi = TMath::Abs( bachTrackXi->GetD( lBestPrimaryVtxPos[0], | |
2082 | lBestPrimaryVtxPos[1], | |
2083 | lMagneticField ) ); | |
2084 | // Note : AliExternalTrackParam::GetD returns an algebraic value ... | |
2085 | ||
2086 | xi->GetXYZ( lPosV0Xi[0], lPosV0Xi[1], lPosV0Xi[2] ); | |
2087 | lV0RadiusXi = TMath::Sqrt( lPosV0Xi[0]*lPosV0Xi[0] + lPosV0Xi[1]*lPosV0Xi[1] ); | |
2088 | ||
2089 | lDcaPosToPrimVertexXi = TMath::Abs( pTrackXi ->GetD( lBestPrimaryVtxPos[0], | |
2090 | lBestPrimaryVtxPos[1], | |
2091 | lMagneticField ) ); | |
2092 | ||
2093 | lDcaNegToPrimVertexXi = TMath::Abs( nTrackXi ->GetD( lBestPrimaryVtxPos[0], | |
2094 | lBestPrimaryVtxPos[1], | |
2095 | lMagneticField ) ); | |
7e504402 | 2096 | |
32599cd4 | 2097 | // - II.Step 3' : extra-selection for cascade candidates |
2098 | // Towards optimisation of AA selection | |
7e504402 | 2099 | // FIXME |
9d670198 | 2100 | if(fkExtraSelections){ |
32599cd4 | 2101 | // if(lChi2Xi > 2000) continue; |
2102 | // if(lV0Chi2Xi > 2000) continue; | |
2103 | ||
2104 | if(lDcaXiDaughters > 0.05) continue; // > 0.1 by default | |
2105 | //if(lXiCosineOfPointingAngle < 0.999 ) continue; | |
c525f360 | 2106 | if(lXiRadius2D < 1.0) continue; |
2107 | if(lXiRadius2D > 100) continue; | |
2108 | if(TMath::Abs(lInvMassLambdaAsCascDghter-1.11568) > 0.008) continue; | |
32599cd4 | 2109 | if(lDcaV0DaughtersXi > 0.3) continue; |
2110 | ||
2111 | if(lV0CosineOfPointingAngleXi > 0.9999) continue; | |
2112 | //if(lDcaV0ToPrimVertexXi < 0.09) continue; | |
2113 | if(lDcaBachToPrimVertexXi < 0.04) continue; | |
2114 | ||
2115 | if(lV0RadiusXi < 1.0) continue; | |
2116 | if(lV0RadiusXi > 100) continue; | |
2117 | //if(lDcaPosToPrimVertexXi < 0.6) continue; | |
2118 | //if(lDcaNegToPrimVertexXi < 0.6) continue; | |
2119 | } | |
2120 | ||
2121 | ||
2122 | ||
45ee0bcc | 2123 | // - II.Step 4 : around effective masses (ESD) |
2124 | // ~ change mass hypotheses to cover all the possibilities : Xi-/+, Omega -/+ | |
2125 | //------------- | |
2126 | ||
2dd4214d | 2127 | |
2dd4214d | 2128 | if( bachTrackXi->Charge() < 0 ) { |
2129 | lV0quality = 0.; | |
2130 | xi->ChangeMassHypothesis(lV0quality , 3312); | |
2131 | // Calculate the effective mass of the Xi- candidate. | |
2132 | // pdg code 3312 = Xi- | |
2133 | lInvMassXiMinus = xi->GetEffMassXi(); | |
45ee0bcc | 2134 | |
2dd4214d | 2135 | lV0quality = 0.; |
2136 | xi->ChangeMassHypothesis(lV0quality , 3334); | |
2137 | // Calculate the effective mass of the Xi- candidate. | |
2138 | // pdg code 3334 = Omega- | |
2139 | lInvMassOmegaMinus = xi->GetEffMassXi(); | |
2140 | ||
2141 | lV0quality = 0.; | |
2142 | xi->ChangeMassHypothesis(lV0quality , 3312); // Back to default hyp. | |
f87cd3db | 2143 | }// end if negative bachelor |
2144 | ||
2dd4214d | 2145 | |
2146 | if( bachTrackXi->Charge() > 0 ){ | |
2147 | lV0quality = 0.; | |
2148 | xi->ChangeMassHypothesis(lV0quality , -3312); | |
2149 | // Calculate the effective mass of the Xi+ candidate. | |
2150 | // pdg code -3312 = Xi+ | |
2151 | lInvMassXiPlus = xi->GetEffMassXi(); | |
2152 | ||
2153 | lV0quality = 0.; | |
2154 | xi->ChangeMassHypothesis(lV0quality , -3334); | |
2155 | // Calculate the effective mass of the Xi+ candidate. | |
2156 | // pdg code -3334 = Omega+ | |
2157 | lInvMassOmegaPlus = xi->GetEffMassXi(); | |
2158 | ||
2159 | lV0quality = 0.; | |
2160 | xi->ChangeMassHypothesis(lV0quality , -3312); // Back to "default" hyp. | |
f87cd3db | 2161 | }// end if positive bachelor |
2162 | ||
2163 | ||
2164 | ||
32599cd4 | 2165 | // - II.Step 5 : PID on the daughter tracks |
f87cd3db | 2166 | //------------- |
2167 | ||
32599cd4 | 2168 | // A - Combined PID |
a786bd57 | 2169 | // Reasonable guess for the priors for the cascade track sample (e-, mu, pi, K, p) |
7475515f | 2170 | Double_t lPriorsGuessXi [AliPID::kSPECIESCN] = {0, 0, 2, 0, 1, 0,0,0,0}; |
2171 | Double_t lPriorsGuessOmega[AliPID::kSPECIESCN] = {0, 0, 1, 1, 1, 0,0,0,0}; | |
2172 | // Trick Jan, 15th, 2013 : Coverity fix - defect 11329 | |
2173 | // If the boolean for charged option is activated in AliPID::SetPriors, only the first 9 double_t in the array are read anyway. | |
2174 | // But Coverity sniffs an issue (a priori safe but the logic is a bit twisted; Coverity is not totally wrong). | |
2175 | // Use AliPID::kSPECIESCN (=14) is size of the array + initialise what I am interested in (charged particle = 9 first particles) | |
2176 | // The rest of the table will be initialised to 0 anyway. | |
2177 | ||
3aedd4a5 | 2178 | |
32599cd4 | 2179 | // Combined VO-positive-daughter PID |
a786bd57 | 2180 | AliPID pPidXi; pPidXi.SetPriors( lPriorsGuessXi , kTRUE); // kTRUE = for charged particle PID |
2181 | AliPID pPidOmega; pPidOmega.SetPriors( lPriorsGuessOmega , kTRUE); // kTRUE = for charged particle PID | |
3aedd4a5 | 2182 | |
2183 | if( pTrackXi->IsOn(AliESDtrack::kESDpid) ){ // Combined PID exists | |
2184 | Double_t r[10] = {0.}; pTrackXi->GetESDpid(r); | |
2185 | pPidXi.SetProbabilities(r); | |
2186 | pPidOmega.SetProbabilities(r); | |
32599cd4 | 2187 | |
3aedd4a5 | 2188 | // Check if the V0 positive track is a proton (case for Xi-) |
2189 | Double_t pproton = pPidXi.GetProbability(AliPID::kProton); | |
2190 | if (pproton > pPidXi.GetProbability(AliPID::kElectron) && | |
2191 | pproton > pPidXi.GetProbability(AliPID::kMuon) && | |
2192 | pproton > pPidXi.GetProbability(AliPID::kPion) && | |
2193 | pproton > pPidXi.GetProbability(AliPID::kKaon) ) lIsPosInXiProton = kTRUE; | |
32599cd4 | 2194 | |
2195 | // Check if the V0 positive track is a pi+ (case for Xi+) | |
2196 | Double_t ppion = pPidXi.GetProbability(AliPID::kPion); | |
2197 | if (ppion > pPidXi.GetProbability(AliPID::kElectron) && | |
2198 | ppion > pPidXi.GetProbability(AliPID::kMuon) && | |
2199 | ppion > pPidXi.GetProbability(AliPID::kKaon) && | |
2200 | ppion > pPidXi.GetProbability(AliPID::kProton) ) lIsPosInXiPion = kTRUE; | |
2201 | ||
2202 | ||
3aedd4a5 | 2203 | // Check if the V0 positive track is a proton (case for Omega-) |
2204 | pproton = 0.; | |
32599cd4 | 2205 | pproton = pPidOmega.GetProbability(AliPID::kProton); |
3aedd4a5 | 2206 | if (pproton > pPidOmega.GetProbability(AliPID::kElectron) && |
2207 | pproton > pPidOmega.GetProbability(AliPID::kMuon) && | |
2208 | pproton > pPidOmega.GetProbability(AliPID::kPion) && | |
2209 | pproton > pPidOmega.GetProbability(AliPID::kKaon) ) lIsPosInOmegaProton = kTRUE; | |
32599cd4 | 2210 | |
2211 | // Check if the V0 positive track is a pi+ (case for Omega+) | |
2212 | ppion = 0.; | |
2213 | ppion = pPidOmega.GetProbability(AliPID::kPion); | |
2214 | if (ppion > pPidOmega.GetProbability(AliPID::kElectron) && | |
2215 | ppion > pPidOmega.GetProbability(AliPID::kMuon) && | |
2216 | ppion > pPidOmega.GetProbability(AliPID::kKaon) && | |
2217 | ppion > pPidOmega.GetProbability(AliPID::kProton) ) lIsPosInOmegaPion = kTRUE; | |
2218 | ||
3aedd4a5 | 2219 | }// end if V0 positive track with existing combined PID |
2220 | ||
32599cd4 | 2221 | |
2222 | // Combined VO-negative-daughter PID | |
a786bd57 | 2223 | AliPID nPidXi; nPidXi.SetPriors( lPriorsGuessXi , kTRUE); // kTRUE = for charged particle PID |
2224 | AliPID nPidOmega; nPidOmega.SetPriors( lPriorsGuessOmega , kTRUE); // kTRUE = for charged particle PID | |
3aedd4a5 | 2225 | |
2226 | if( nTrackXi->IsOn(AliESDtrack::kESDpid) ){ // Combined PID exists | |
2227 | Double_t r[10] = {0.}; nTrackXi->GetESDpid(r); | |
2228 | nPidXi.SetProbabilities(r); | |
2229 | nPidOmega.SetProbabilities(r); | |
32599cd4 | 2230 | |
2231 | // Check if the V0 negative track is a pi- (case for Xi-) | |
2232 | Double_t ppion = nPidXi.GetProbability(AliPID::kPion); | |
2233 | if (ppion > nPidXi.GetProbability(AliPID::kElectron) && | |
2234 | ppion > nPidXi.GetProbability(AliPID::kMuon) && | |
2235 | ppion > nPidXi.GetProbability(AliPID::kKaon) && | |
2236 | ppion > nPidXi.GetProbability(AliPID::kProton) ) lIsNegInXiPion = kTRUE; | |
2237 | ||
3aedd4a5 | 2238 | // Check if the V0 negative track is an anti-proton (case for Xi+) |
2239 | Double_t pproton = nPidXi.GetProbability(AliPID::kProton); | |
2240 | if (pproton > nPidXi.GetProbability(AliPID::kElectron) && | |
2241 | pproton > nPidXi.GetProbability(AliPID::kMuon) && | |
2242 | pproton > nPidXi.GetProbability(AliPID::kPion) && | |
2243 | pproton > nPidXi.GetProbability(AliPID::kKaon) ) lIsNegInXiProton = kTRUE; | |
32599cd4 | 2244 | |
2245 | // Check if the V0 negative track is a pi- (case for Omega-) | |
2246 | ppion = 0.; | |
2247 | ppion = nPidOmega.GetProbability(AliPID::kPion); | |
2248 | if (ppion > nPidOmega.GetProbability(AliPID::kElectron) && | |
2249 | ppion > nPidOmega.GetProbability(AliPID::kMuon) && | |
2250 | ppion > nPidOmega.GetProbability(AliPID::kKaon) && | |
2251 | ppion > nPidOmega.GetProbability(AliPID::kProton) ) lIsNegInOmegaPion = kTRUE; | |
2252 | ||
3aedd4a5 | 2253 | // Check if the V0 negative track is an anti-proton (case for Omega+) |
2254 | pproton = 0.; | |
2255 | pproton = nPidOmega.GetProbability(AliPID::kProton); | |
2256 | if (pproton > nPidOmega.GetProbability(AliPID::kElectron) && | |
2257 | pproton > nPidOmega.GetProbability(AliPID::kMuon) && | |
2258 | pproton > nPidOmega.GetProbability(AliPID::kPion) && | |
2259 | pproton > nPidOmega.GetProbability(AliPID::kKaon) ) lIsNegInOmegaProton = kTRUE; | |
32599cd4 | 2260 | |
3aedd4a5 | 2261 | }// end if V0 negative track with existing combined PID |
2262 | ||
2263 | ||
32599cd4 | 2264 | // Combined bachelor PID |
a786bd57 | 2265 | AliPID bachPidXi; bachPidXi.SetPriors( lPriorsGuessXi , kTRUE ); // kTRUE = for charged particle PID |
2266 | AliPID bachPidOmega; bachPidOmega.SetPriors( lPriorsGuessOmega, kTRUE ); // kTRUE = for charged particle PID | |
f87cd3db | 2267 | |
2268 | if( bachTrackXi->IsOn(AliESDtrack::kESDpid) ){ // Combined PID exists | |
3aedd4a5 | 2269 | Double_t r[10] = {0.}; bachTrackXi->GetESDpid(r); |
2270 | bachPidXi.SetProbabilities(r); | |
2271 | bachPidOmega.SetProbabilities(r); | |
f87cd3db | 2272 | // Check if the bachelor track is a pion |
3aedd4a5 | 2273 | Double_t ppion = bachPidXi.GetProbability(AliPID::kPion); |
2274 | if (ppion > bachPidXi.GetProbability(AliPID::kElectron) && | |
2275 | ppion > bachPidXi.GetProbability(AliPID::kMuon) && | |
2276 | ppion > bachPidXi.GetProbability(AliPID::kKaon) && | |
2277 | ppion > bachPidXi.GetProbability(AliPID::kProton) ) lIsBachelorPion = kTRUE; | |
f87cd3db | 2278 | // Check if the bachelor track is a kaon |
3aedd4a5 | 2279 | Double_t pkaon = bachPidOmega.GetProbability(AliPID::kKaon); |
2280 | if (pkaon > bachPidOmega.GetProbability(AliPID::kElectron) && | |
2281 | pkaon > bachPidOmega.GetProbability(AliPID::kMuon) && | |
2282 | pkaon > bachPidOmega.GetProbability(AliPID::kPion) && | |
2283 | pkaon > bachPidOmega.GetProbability(AliPID::kProton) ) lIsBachelorKaon = kTRUE; | |
f87cd3db | 2284 | }// end if bachelor track with existing combined PID |
2dd4214d | 2285 | |
3aedd4a5 | 2286 | |
32599cd4 | 2287 | // B - TPC PID : 3-sigma bands on Bethe-Bloch curve |
32599cd4 | 2288 | |
e8fd19de | 2289 | // Bachelor |
2290 | if (TMath::Abs(fESDpid->NumberOfSigmasTPC( bachTrackXi,AliPID::kKaon)) < 4) lIsBachelorKaonForTPC = kTRUE; | |
2291 | if (TMath::Abs(fESDpid->NumberOfSigmasTPC( bachTrackXi,AliPID::kPion)) < 4) lIsBachelorPionForTPC = kTRUE; | |
2292 | ||
2293 | // Negative V0 daughter | |
2294 | if (TMath::Abs(fESDpid->NumberOfSigmasTPC( nTrackXi,AliPID::kPion )) < 4) lIsNegPionForTPC = kTRUE; | |
2295 | if (TMath::Abs(fESDpid->NumberOfSigmasTPC( nTrackXi,AliPID::kProton )) < 4) lIsNegProtonForTPC = kTRUE; | |
2296 | ||
2297 | // Positive V0 daughter | |
2298 | if (TMath::Abs(fESDpid->NumberOfSigmasTPC( pTrackXi,AliPID::kPion )) < 4) lIsPosPionForTPC = kTRUE; | |
2299 | if (TMath::Abs(fESDpid->NumberOfSigmasTPC( pTrackXi,AliPID::kProton )) < 4) lIsPosProtonForTPC = kTRUE; | |
2300 | ||
2301 | /* | |
2302 | const AliExternalTrackParam *pInnerWallTrackXi = pTrackXi ->GetInnerParam(); | |
2303 | const AliExternalTrackParam *nInnerWallTrackXi = nTrackXi ->GetInnerParam(); | |
2304 | const AliExternalTrackParam *bachInnerWallTrackXi = bachTrackXi ->GetInnerParam(); | |
2305 | if(pInnerWallTrackXi && nInnerWallTrackXi && bachInnerWallTrackXi ){ | |
2306 | ||
2307 | Double_t pMomInnerWall = pInnerWallTrackXi ->GetP(); | |
2308 | Double_t nMomInnerWall = nInnerWallTrackXi ->GetP(); | |
2309 | Double_t bachMomInnerWall = bachInnerWallTrackXi->GetP(); | |
2310 | ||
2311 | // Bachelor | |
2312 | if (TMath::Abs(fESDpid->NumberOfSigmasTPC( bachTrackXi,AliPID::kPion)) < 3) lIsBachelorPionForTPC = kTRUE; | |
2313 | if (bachMomInnerWall < 0.350 && TMath::Abs(fESDpid->NumberOfSigmasTPC( bachTrackXi,AliPID::kKaon)) < 5) lIsBachelorKaonForTPC = kTRUE; | |
2314 | if (bachMomInnerWall > 0.350 && TMath::Abs(fESDpid->NumberOfSigmasTPC( bachTrackXi,AliPID::kKaon)) < 3) lIsBachelorKaonForTPC = kTRUE; | |
2315 | ||
2316 | // Negative V0 daughter | |
2317 | if (TMath::Abs(fESDpid->NumberOfSigmasTPC( nTrackXi,AliPID::kPion )) < 3 ) lIsNegPionForTPC = kTRUE; | |
2318 | if (nMomInnerWall < 0.6 && TMath::Abs(fESDpid->NumberOfSigmasTPC( nTrackXi,AliPID::kProton ) ) < 5 ) lIsNegProtonForTPC = kTRUE; | |
2319 | if (nMomInnerWall > 0.6 && TMath::Abs(fESDpid->NumberOfSigmasTPC( nTrackXi,AliPID::kProton ) ) < 3 ) lIsNegProtonForTPC = kTRUE; | |
2320 | ||
2321 | // Positive V0 daughter | |
2322 | if (TMath::Abs(fESDpid->NumberOfSigmasTPC( pTrackXi,AliPID::kPion )) < 3 ) lIsPosPionForTPC = kTRUE; | |
2323 | if (pMomInnerWall < 0.6 && TMath::Abs(fESDpid->NumberOfSigmasTPC( pTrackXi,AliPID::kProton )) < 5) lIsPosProtonForTPC = kTRUE; | |
2324 | if (pMomInnerWall > 0.6 && TMath::Abs(fESDpid->NumberOfSigmasTPC( pTrackXi,AliPID::kProton )) < 3) lIsPosProtonForTPC = kTRUE; | |
2325 | } | |
2326 | */ | |
2327 | ||
2328 | ||
2329 | ||
f87cd3db | 2330 | |
2331 | // - II.Step 6 : extra info for QA (ESD) | |
3873f6ee | 2332 | // miscellaneous pieces of info that may help regarding data quality assessment. |
45ee0bcc | 2333 | //------------- |
2334 | ||
2335 | xi->GetPxPyPz( lXiMomX, lXiMomY, lXiMomZ ); | |
2336 | lXiTransvMom = TMath::Sqrt( lXiMomX*lXiMomX + lXiMomY*lXiMomY ); | |
2337 | lXiTotMom = TMath::Sqrt( lXiMomX*lXiMomX + lXiMomY*lXiMomY + lXiMomZ*lXiMomZ ); | |
2338 | ||
2339 | xi->GetBPxPyPz( lBachMomX, lBachMomY, lBachMomZ ); | |
2340 | lBachTransvMom = TMath::Sqrt( lBachMomX*lBachMomX + lBachMomY*lBachMomY ); | |
2341 | lBachTotMom = TMath::Sqrt( lBachMomX*lBachMomX + lBachMomY*lBachMomY + lBachMomZ*lBachMomZ ); | |
2342 | ||
2343 | lChargeXi = xi->Charge(); | |
2344 | ||
2345 | lV0toXiCosineOfPointingAngle = xi->GetV0CosineOfPointingAngle( lPosXi[0], lPosXi[1], lPosXi[2] ); | |
2346 | ||
2347 | lRapXi = xi->RapXi(); | |
2348 | lRapOmega = xi->RapOmega(); | |
2349 | lEta = xi->Eta(); | |
2350 | lTheta = xi->Theta() *180.0/TMath::Pi(); | |
2351 | lPhi = xi->Phi() *180.0/TMath::Pi(); | |
2352 | lAlphaXi = xi->AlphaXi(); | |
2353 | lPtArmXi = xi->PtArmXi(); | |
2354 | ||
2355 | ||
f3dc9369 | 2356 | //FIXME : Extra-cut = Anti-splitting cut for lambda daughters |
7e504402 | 2357 | Bool_t kAntiSplittingLambda = kFALSE; |
f3dc9369 | 2358 | |
2359 | if(kAntiSplittingLambda){ | |
2360 | Double_t lNMomX = 0., lNMomY = 0., lNMomZ = 0.; | |
2361 | Double_t lPMomX = 0., lPMomY = 0., lPMomZ = 0.; | |
2362 | ||
2363 | xi->GetPPxPyPz(lPMomX, lPMomY, lPMomZ); | |
2364 | xi->GetNPxPyPz(lNMomX, lNMomY, lNMomZ); | |
2365 | ||
2366 | if( xi->Charge() < 0){// Xi- or Omega- | |
7e504402 | 2367 | if(TMath::Abs(lBachTransvMom - TMath::Sqrt( lNMomX*lNMomX + lNMomY*lNMomY ) ) < 0.075) continue; |
f3dc9369 | 2368 | } |
2369 | else { //Xi+ or Omega+ | |
7e504402 | 2370 | if(TMath::Abs(lBachTransvMom - TMath::Sqrt( lPMomX*lPMomX + lPMomY*lPMomY ) ) < 0.075) continue; |
f3dc9369 | 2371 | } |
2372 | } | |
2373 | ||
e8fd19de | 2374 | //FIXME : Just to know which file is currently open : locate the file containing Xi |
2375 | // cout << "Name of the file containing Xi candidate(s) :" | |
2376 | // << CurrentFileName() | |
2377 | // << " / entry: " << Entry() | |
2378 | // << " / in file: " << lESDevent->GetEventNumberInFile() // <- Cvetan / From Mihaela: AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()->GetTree()->GetReadEntry(); | |
2379 | // << " : mass(Xi) = " << xi->GetEffMassXi() | |
2380 | // << " / charge = " << lChargeXi | |
2381 | // << " / pt(Casc) = " << lXiTransvMom | |
c525f360 | 2382 | // << " / Decay 2d R(Xi) = " << lXiRadius2D |
e8fd19de | 2383 | // << " / Track Index(Pos) = " << lIdxPosXi << "/ Nb(TPC clusters) = " << lPosTPCClusters |
2384 | // << " / Track Index(Neg) = " << lIdxNegXi << "/ Nb(TPC clusters) = " << lNegTPCClusters | |
2385 | // << " / Track Index(Bach) = " << lBachIdx << "/ Nb(TPC clusters) = " << lBachTPCClusters | |
2386 | // << endl; | |
90926f02 | 2387 | |
32599cd4 | 2388 | |
2389 | // II.Step 7 - Complementary info for monitoring the cascade cut variables | |
2390 | ||
a786bd57 | 2391 | //const AliMultiplicity *lAliMult = lESDevent->GetMultiplicity(); |
2392 | //lSPDTrackletsMultiplicity = lAliMult->GetNumberOfTracklets(); | |
2393 | // --> now provided by AliESDEvent::EstimateMultiplicity ... | |
32599cd4 | 2394 | |
2395 | // II.Step 8 - Azimuthal correlation study | |
3873f6ee | 2396 | //------------- |
2397 | ||
2398 | lTVect3MomXi.SetXYZ( lXiMomX, lXiMomY, lXiMomZ ); | |
2399 | lArrTrackID[0] = pTrackXi ->GetID(); | |
2400 | lArrTrackID[1] = nTrackXi ->GetID(); | |
2401 | lArrTrackID[2] = bachTrackXi->GetID(); | |
2402 | ||
2403 | ||
45ee0bcc | 2404 | }// end of ESD treatment |
2405 | ||
2406 | ||
2407 | if(fAnalysisType == "AOD"){ | |
2408 | ||
2409 | // ------------------------------------- | |
32599cd4 | 2410 | // II.AOD - Calcultaion Part dedicated to Xi vertices (ESD) |
45ee0bcc | 2411 | |
2412 | const AliAODcascade *xi = lAODevent->GetCascade(iXi); | |
2413 | if (!xi) continue; | |
2414 | ||
2415 | // Just to know which file is currently open : locate the file containing Xi | |
2416 | // cout << "Name of the file containing Xi candidate(s) :" << fesdH->GetTree()->GetCurrentFile()->GetName() << endl; | |
2417 | ||
2418 | ||
2419 | // - II.Step 1 : Characteristics of the event : prim. Vtx + magnetic field (AOD) | |
2420 | //------------- | |
2421 | ||
45ee0bcc | 2422 | |
9d670198 | 2423 | lTrkgPrimaryVtxRadius3D = -500. ; |
2424 | // FIXME : We don't have the different prim. vertex at the AOD level -> nothing to do. | |
45ee0bcc | 2425 | |
45ee0bcc | 2426 | lBestPrimaryVtxRadius3D = TMath::Sqrt( lBestPrimaryVtxPos[0] * lBestPrimaryVtxPos[0] + |
2427 | lBestPrimaryVtxPos[1] * lBestPrimaryVtxPos[1] + | |
2428 | lBestPrimaryVtxPos[2] * lBestPrimaryVtxPos[2] ); | |
2429 | ||
2430 | ||
2431 | // - II.Step 2 : Assigning the necessary variables for specific AliAODcascade data members (AOD) | |
2432 | //------------- | |
2433 | ||
2434 | lEffMassXi = xi->MassXi(); // default working hypothesis : cascade = Xi- decay | |
2435 | lChi2Xi = xi->Chi2Xi(); | |
2436 | lDcaXiDaughters = xi->DcaXiDaughters(); | |
2437 | lXiCosineOfPointingAngle = xi->CosPointingAngleXi( lBestPrimaryVtxPos[0], | |
2438 | lBestPrimaryVtxPos[1], | |
2439 | lBestPrimaryVtxPos[2] ); | |
2440 | // Take care : | |
2441 | // the best available vertex should be used (like in AliCascadeVertexer) | |
2442 | ||
2443 | lPosXi[0] = xi->DecayVertexXiX(); | |
2444 | lPosXi[1] = xi->DecayVertexXiY(); | |
2445 | lPosXi[2] = xi->DecayVertexXiZ(); | |
c525f360 | 2446 | lXiRadius2D = TMath::Sqrt( lPosXi[0]*lPosXi[0] + lPosXi[1]*lPosXi[1] ); |
2447 | lXiRadius3D = TMath::Sqrt( lPosXi[0]*lPosXi[0] + lPosXi[1]*lPosXi[1] + lPosXi[2]*lPosXi[2] ); | |
45ee0bcc | 2448 | |
2449 | ||
2450 | // - II.Step 3 : around the tracks : Bach + V0 (AOD) | |
2451 | // ~ Necessary variables for AODcascade data members coming from the AODv0 part (inheritance) | |
2452 | //------------- | |
2453 | ||
2454 | lChargeXi = xi->ChargeXi(); | |
2455 | ||
2456 | if( lChargeXi < 0) | |
2457 | lInvMassLambdaAsCascDghter = xi->MassLambda(); | |
2458 | else | |
2459 | lInvMassLambdaAsCascDghter = xi->MassAntiLambda(); | |
2460 | ||
2461 | lDcaV0DaughtersXi = xi->DcaV0Daughters(); | |
2462 | lV0Chi2Xi = xi->Chi2V0(); | |
2463 | ||
2464 | lV0CosineOfPointingAngleXi = xi->CosPointingAngle( lBestPrimaryVtxPos ); | |
2465 | lDcaV0ToPrimVertexXi = xi->DcaV0ToPrimVertex(); | |
2466 | ||
2467 | lDcaBachToPrimVertexXi = xi->DcaBachToPrimVertex(); | |
2468 | ||
2469 | ||
2470 | lPosV0Xi[0] = xi->DecayVertexV0X(); | |
2471 | lPosV0Xi[1] = xi->DecayVertexV0Y(); | |
2472 | lPosV0Xi[2] = xi->DecayVertexV0Z(); | |
2473 | lV0RadiusXi = TMath::Sqrt( lPosV0Xi[0]*lPosV0Xi[0] + lPosV0Xi[1]*lPosV0Xi[1] ); | |
2474 | ||
2475 | lDcaPosToPrimVertexXi = xi->DcaPosToPrimVertex(); | |
2476 | lDcaNegToPrimVertexXi = xi->DcaNegToPrimVertex(); | |
2477 | ||
2478 | ||
2479 | // - II.Step 4 : around effective masses (AOD) | |
2480 | // ~ change mass hypotheses to cover all the possibilities : Xi-/+, Omega -/+ | |
2481 | //------------- | |
2482 | ||
2483 | if( lChargeXi < 0 ) lInvMassXiMinus = xi->MassXi(); | |
2484 | if( lChargeXi > 0 ) lInvMassXiPlus = xi->MassXi(); | |
2485 | if( lChargeXi < 0 ) lInvMassOmegaMinus = xi->MassOmega(); | |
2486 | if( lChargeXi > 0 ) lInvMassOmegaPlus = xi->MassOmega(); | |
2487 | ||
2488 | ||
32599cd4 | 2489 | // - II.Step 5 : PID on the daughters (To be developed ...) |
f87cd3db | 2490 | //------------- |
2491 | ||
32599cd4 | 2492 | // Combined PID |
2493 | ||
3873f6ee | 2494 | /* |
ff0753e1 | 2495 | // Reasonable guess for the priors for the cascade track sample |
a786bd57 | 2496 | Double_t lPriorsGuessXi[5] = {0.0, 0.0, 2, 0, 1, 0,0,0,0,0}; |
2497 | Double_t lPriorsGuessOmega[5] = {0.0, 0.0, 1, 1, 1, 0,0,0,0,0}; | |
2498 | AliPID bachPidXi; bachPidXi.SetPriors( lPriorsGuessXi , kTRUE); // kTRUE = for charged particle PID | |
2499 | AliPID bachPidOmega; bachPidOmega.SetPriors( lPriorsGuessOmega , kTRUE); // kTRUE = for charged particle PID | |
ff0753e1 | 2500 | |
3873f6ee | 2501 | const AliAODTrack *bachTrackXi = lAODevent->GetTrack( xi->GetBachID() ); // FIXME : GetBachID not implemented ? |
ff0753e1 | 2502 | |
2503 | if( bachTrackXi->IsOn(AliESDtrack::kESDpid) ){ // Combined PID exists, the AOD flags = a copy of the ESD ones | |
2504 | Double_t r[10]; bachTrackXi->GetPID(r); | |
3aedd4a5 | 2505 | bachPidXi.SetProbabilities(r); |
2506 | bachPidOmega.SetProbabilities(r); | |
ff0753e1 | 2507 | // Check if the bachelor track is a pion |
3aedd4a5 | 2508 | Double_t ppion = bachPidXi.GetProbability(AliPID::kPion); |
2509 | if (ppion > bachPidXi.GetProbability(AliPID::kElectron) && | |
2510 | ppion > bachPidXi.GetProbability(AliPID::kMuon) && | |
2511 | ppion > bachPidXi.GetProbability(AliPID::kKaon) && | |
2512 | ppion > bachPidXi.GetProbability(AliPID::kProton) ) lIsBachelorPion = kTRUE; | |
ff0753e1 | 2513 | // Check if the bachelor track is a kaon |
3aedd4a5 | 2514 | Double_t pkaon = bachPidOmega.GetProbability(AliPID::kKaon); |
2515 | if (pkaon > bachPidOmega.GetProbability(AliPID::kElectron) && | |
2516 | pkaon > bachPidOmega.GetProbability(AliPID::kMuon) && | |
2517 | pkaon > bachPidOmega.GetProbability(AliPID::kPion) && | |
2518 | pkaon > bachPidOmega.GetProbability(AliPID::kProton) ) lIsBachelorKaon = kTRUE; | |
ff0753e1 | 2519 | |
2520 | }// end if bachelor track with existing combined PID | |
2521 | */ | |
f87cd3db | 2522 | |
32599cd4 | 2523 | // TPC PID |
f87cd3db | 2524 | |
2525 | // - II.Step 6 : extra info for QA (AOD) | |
45ee0bcc | 2526 | // miscellaneous pieces onf info that may help regarding data quality assessment. |
2527 | //------------- | |
2528 | ||
2529 | lXiMomX = xi->MomXiX(); | |
2530 | lXiMomY = xi->MomXiY(); | |
2531 | lXiMomZ = xi->MomXiZ(); | |
2532 | lXiTransvMom = TMath::Sqrt( lXiMomX*lXiMomX + lXiMomY*lXiMomY ); | |
2533 | lXiTotMom = TMath::Sqrt( lXiMomX*lXiMomX + lXiMomY*lXiMomY + lXiMomZ*lXiMomZ ); | |
2534 | ||
2535 | lBachMomX = xi->MomBachX(); | |
2536 | lBachMomY = xi->MomBachY(); | |
2537 | lBachMomZ = xi->MomBachZ(); | |
2538 | lBachTransvMom = TMath::Sqrt( lBachMomX*lBachMomX + lBachMomY*lBachMomY ); | |
2539 | lBachTotMom = TMath::Sqrt( lBachMomX*lBachMomX + lBachMomY*lBachMomY + lBachMomZ*lBachMomZ ); | |
2540 | ||
2541 | ||
2542 | lV0toXiCosineOfPointingAngle = xi->CosPointingAngle( xi->GetDecayVertexXi() ); | |
2543 | ||
2544 | lRapXi = xi->RapXi(); | |
2545 | lRapOmega = xi->RapOmega(); | |
2546 | lEta = xi->Eta(); // Will not work ! need a method Pz(), Py() Px() | |
2547 | lTheta = xi->Theta() *180.0/TMath::Pi(); // in AODcascade. | |
2548 | lPhi = xi->Phi() *180.0/TMath::Pi(); // Here, we will get eta, theta, phi for the V0 ... | |
2549 | lAlphaXi = xi->AlphaXi(); | |
2550 | lPtArmXi = xi->PtArmXi(); | |
2551 | ||
32599cd4 | 2552 | // II.Step 7 - Complementary info for monitoring the cascade cut variables |
f3dc9369 | 2553 | //FIXME : missing for AOD : Tacklet Multiplicity + TPCCluster |
3873f6ee | 2554 | |
32599cd4 | 2555 | // II.Step 8 - Azimuthal correlation study |
3873f6ee | 2556 | //------------- |
2557 | ||
2558 | lTVect3MomXi.SetXYZ( lXiMomX, lXiMomY, lXiMomZ ); | |
2559 | ||
2560 | AliAODTrack *pTrackXi = dynamic_cast<AliAODTrack*>( xi->GetDaughter(0) ); | |
2561 | AliAODTrack *nTrackXi = dynamic_cast<AliAODTrack*>( xi->GetDaughter(1) ); | |
2562 | AliAODTrack *bachTrackXi = dynamic_cast<AliAODTrack*>( xi->GetDecayVertexXi()->GetDaughter(0) ); | |
2563 | if (!pTrackXi || !nTrackXi || !bachTrackXi ) { | |
7e504402 | 2564 | AliWarning("ERROR: Could not retrieve one of the 3 AOD daughter tracks of the cascade ..."); |
3873f6ee | 2565 | continue; |
2566 | } | |
2567 | ||
2568 | lArrTrackID[0] = pTrackXi ->GetID(); | |
2569 | lArrTrackID[1] = nTrackXi ->GetID(); | |
2570 | lArrTrackID[2] = bachTrackXi->GetID(); | |
2571 | ||
45ee0bcc | 2572 | }// end of AOD treatment |
2573 | ||
2574 | ||
32599cd4 | 2575 | // ------------------------------------- |
9d670198 | 2576 | // II.Fill - Filling the TH1,2,3Fs, HnSparses, CFContainers, FOR events with CASCADES ! |
e8fd19de | 2577 | |
2578 | // if( lIsBachelorKaonForTPC ) | |
2579 | // // FIXME : Just to know which file is currently open : locate the file containing Xi | |
2580 | // cout << "Name of the file containing Omega candidate(s) :" | |
2581 | // << CurrentFileName() | |
2582 | // << " / entry: " << Entry() | |
2583 | // << " / in file: " << lESDevent->GetEventNumberInFile() // <- Cvetan / From Mihaela: AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()->GetTree()->GetReadEntry(); | |
2584 | // << " : mass(Omega+) = " << lInvMassOmegaPlus | |
2585 | // << " : mass(Omega-) = " << lInvMassOmegaMinus | |
2586 | // << " / charge = " << lChargeXi | |
2587 | // << " / pt(Casc) = " << lXiTransvMom | |
c525f360 | 2588 | // << " / Decay 2d R(Xi) = " << lXiRadius2D |
e8fd19de | 2589 | // << endl; |
45ee0bcc | 2590 | |
2591 | ||
32599cd4 | 2592 | // - II.Fill.Step 1 : primary vertex |
7e504402 | 2593 | |
2594 | fHistTPCrefitTrackMultiplicityForCascadeEvt->Fill( nTrackWithTPCrefitMultiplicity ); | |
c525f360 | 2595 | fHistPrimaryTrackMultiplicityForCascadeEvt ->Fill( nTrackPrimaryMultiplicity ); |
7e504402 | 2596 | |
2597 | fHistPosV0TPCClusters ->Fill( lPosTPCClusters ); | |
2598 | fHistNegV0TPCClusters ->Fill( lNegTPCClusters ); | |
2599 | fHistBachTPCClusters ->Fill( lBachTPCClusters ); | |
e8fd19de | 2600 | |
2601 | f2dHistTPCdEdxOfCascDghters ->Fill( lInnerWallMomCascDghters[0] , lTPCSignalCascDghters[0] ); | |
2602 | f2dHistTPCdEdxOfCascDghters ->Fill( lInnerWallMomCascDghters[1] , lTPCSignalCascDghters[1] ); | |
2603 | f2dHistTPCdEdxOfCascDghters ->Fill( lInnerWallMomCascDghters[2] , lTPCSignalCascDghters[2] ); | |
2604 | ||
45ee0bcc | 2605 | fHistVtxStatus ->Fill( lStatusTrackingPrimVtx ); // 1 if tracking vtx = ok |
2606 | ||
2607 | if( lStatusTrackingPrimVtx ){ | |
9d670198 | 2608 | fHistPosTrkgPrimaryVtxXForCascadeEvt ->Fill( lTrkgPrimaryVtxPos[0] ); |
2609 | fHistPosTrkgPrimaryVtxYForCascadeEvt ->Fill( lTrkgPrimaryVtxPos[1] ); | |
2610 | fHistPosTrkgPrimaryVtxZForCascadeEvt ->Fill( lTrkgPrimaryVtxPos[2] ); | |
2611 | fHistTrkgPrimaryVtxRadius ->Fill( lTrkgPrimaryVtxRadius3D ); | |
45ee0bcc | 2612 | } |
2613 | ||
9d670198 | 2614 | fHistPosBestPrimaryVtxXForCascadeEvt ->Fill( lBestPrimaryVtxPos[0] ); |
2615 | fHistPosBestPrimaryVtxYForCascadeEvt ->Fill( lBestPrimaryVtxPos[1] ); | |
2616 | fHistPosBestPrimaryVtxZForCascadeEvt ->Fill( lBestPrimaryVtxPos[2] ); | |
2617 | fHistBestPrimaryVtxRadius ->Fill( lBestPrimaryVtxRadius3D ); | |
45ee0bcc | 2618 | |
2619 | f2dHistTrkgPrimVtxVsBestPrimVtx->Fill( lTrkgPrimaryVtxRadius3D, lBestPrimaryVtxRadius3D ); | |
2620 | ||
e8fd19de | 2621 | // **************************** With PID on ? ... for the signal region ? ************FIXME************************************** |
2622 | if( ( (lChargeXi<0) && lIsBachelorPionForTPC && lIsPosProtonForTPC && lIsNegPionForTPC ) || | |
2623 | ( (lChargeXi>0) && lIsBachelorPionForTPC && lIsNegProtonForTPC && lIsPosPionForTPC ) ) | |
2624 | // NOTE : | |
c525f360 | 2625 | // with this PID condition, it could happen that a cascade candidate satisfies the wrong requirement, |
e8fd19de | 2626 | // e.g. one looks at a Xi- candidate for which lIsBachelorPionForTPC && lIsPosProtonForTPC && lIsNegPionForTPC = kFALSE |
2627 | // Expectation: it should be excluded. | |
2628 | // but lIsBachelorPionForTPC && lIsNegProtonForTPC && lIsPosPionForTPC = kTRUE | |
2629 | // then this bad Xi-candidate will contribute anyway (OR condition). | |
2630 | // Hence : the extra condition on the sign of the Cascade | |
2631 | { | |
2632 | // if( TMath::Abs( lInvMassXiMinus-1.3217 ) < 0.010 || TMath::Abs( lInvMassXiPlus-1.3217 ) < 0.010){} } | |
2633 | ||
2634 | // II.Fill.Step 2 | |
2635 | fHistEffMassXi ->Fill( lEffMassXi ); | |
2636 | fHistChi2Xi ->Fill( lChi2Xi ); // Flag CascadeVtxer: Cut Variable a | |
2637 | fHistDcaXiDaughters ->Fill( lDcaXiDaughters ); // Flag CascadeVtxer: Cut Variable e | |
2638 | fHistDcaBachToPrimVertex ->Fill( lDcaBachToPrimVertexXi ); // Flag CascadeVtxer: Cut Variable d | |
2639 | fHistXiCosineOfPointingAngle ->Fill( lXiCosineOfPointingAngle ); // Flag CascadeVtxer: Cut Variable f | |
c525f360 | 2640 | fHistXiRadius ->Fill( lXiRadius2D ); // Flag CascadeVtxer: Cut Variable g+h |
e8fd19de | 2641 | |
2642 | ||
2643 | // II.Fill.Step 3 | |
2644 | fHistMassLambdaAsCascDghter ->Fill( lInvMassLambdaAsCascDghter ); // Flag CascadeVtxer: Cut Variable c | |
2645 | fHistV0Chi2Xi ->Fill( lV0Chi2Xi ); | |
2646 | fHistDcaV0DaughtersXi ->Fill( lDcaV0DaughtersXi ); | |
2647 | fHistV0CosineOfPointingAngleXi ->Fill( lV0CosineOfPointingAngleXi ); | |
2648 | fHistV0RadiusXi ->Fill( lV0RadiusXi ); | |
2649 | ||
2650 | fHistDcaV0ToPrimVertexXi ->Fill( lDcaV0ToPrimVertexXi ); // Flag CascadeVtxer: Cut Variable b | |
2651 | fHistDcaPosToPrimVertexXi ->Fill( lDcaPosToPrimVertexXi ); | |
2652 | fHistDcaNegToPrimVertexXi ->Fill( lDcaNegToPrimVertexXi ); | |
2653 | ||
90926f02 | 2654 | |
e8fd19de | 2655 | // II.Fill.Step 4 : extra QA info |
2656 | ||
2657 | fHistChargeXi ->Fill( lChargeXi ); | |
2658 | fHistV0toXiCosineOfPointingAngle->Fill( lV0toXiCosineOfPointingAngle ); | |
2659 | ||
c525f360 | 2660 | if( TMath::Abs( lInvMassXiMinus-1.3217 ) < 0.010 || TMath::Abs( lInvMassXiPlus-1.3217 ) < 0.010){// One InvMass should be different from 0 |
e8fd19de | 2661 | fHistXiTransvMom ->Fill( lXiTransvMom ); |
2662 | fHistXiTotMom ->Fill( lXiTotMom ); | |
2663 | ||
2664 | fHistBachTransvMomXi ->Fill( lBachTransvMom ); | |
2665 | fHistBachTotMomXi ->Fill( lBachTotMom ); | |
2666 | ||
2667 | fHistRapXi ->Fill( lRapXi ); | |
2668 | fHistEtaXi ->Fill( lEta ); | |
2669 | fHistThetaXi ->Fill( lTheta ); | |
2670 | fHistPhiXi ->Fill( lPhi ); | |
c525f360 | 2671 | |
e8fd19de | 2672 | } |
c525f360 | 2673 | |
2674 | if( (lChargeXi < 0) && (TMath::Abs( lInvMassXiMinus-1.3217 ) < 0.010) ) fHistcTauXiMinus ->Fill( lXiRadius3D * 1.3271/lXiTotMom ); | |
2675 | if( (lChargeXi > 0) && (TMath::Abs( lInvMassXiPlus -1.3217 ) < 0.010) ) fHistcTauXiPlus ->Fill( lXiRadius3D * 1.3271/lXiTotMom ); | |
45ee0bcc | 2676 | |
c525f360 | 2677 | if( TMath::Abs( lInvMassOmegaMinus-1.672 ) < 0.010 || TMath::Abs( lInvMassOmegaPlus-1.672 ) < 0.010 ){// One InvMass should be different from 0 |
e8fd19de | 2678 | fHistRapOmega ->Fill( lRapOmega ); |
c525f360 | 2679 | |
e8fd19de | 2680 | } |
c525f360 | 2681 | |
2682 | if( (lChargeXi < 0) && (TMath::Abs( lInvMassOmegaMinus-1.672 ) < 0.010) ) fHistcTauOmegaMinus ->Fill( lXiRadius3D * 1.67245/lXiTotMom ); | |
2683 | if( (lChargeXi > 0) && (TMath::Abs( lInvMassOmegaPlus- 1.672 ) < 0.010) ) fHistcTauOmegaPlus ->Fill( lXiRadius3D * 1.67245/lXiTotMom ); | |
e8fd19de | 2684 | |
c525f360 | 2685 | f2dHistArmenteros ->Fill( lChargeXi*lAlphaXi, lPtArmXi ); |
e8fd19de | 2686 | }// end with PID ... |
32599cd4 | 2687 | |
2688 | // II.Fill.Step 5 : inv mass plots 1D | |
2689 | if( lChargeXi < 0 ){ | |
2690 | fHistMassXiMinus ->Fill( lInvMassXiMinus ); | |
2691 | fHistMassOmegaMinus ->Fill( lInvMassOmegaMinus ); | |
2692 | if(lIsBachelorPion) fHistMassWithCombPIDXiMinus ->Fill( lInvMassXiMinus ); | |
2693 | if(lIsBachelorKaon) fHistMassWithCombPIDOmegaMinus ->Fill( lInvMassOmegaMinus ); | |
2694 | } | |
2695 | ||
2696 | if( lChargeXi > 0 ){ | |
2697 | fHistMassXiPlus ->Fill( lInvMassXiPlus ); | |
2698 | fHistMassOmegaPlus ->Fill( lInvMassOmegaPlus ); | |
2699 | if(lIsBachelorPion) fHistMassWithCombPIDXiPlus ->Fill( lInvMassXiPlus ); | |
2700 | if(lIsBachelorKaon) fHistMassWithCombPIDOmegaPlus ->Fill( lInvMassOmegaPlus ); | |
2701 | } | |
2702 | ||
2703 | ||
2704 | // II.Fill.Step 6 : inv mass plots 2D, 3D | |
45ee0bcc | 2705 | if( lChargeXi < 0 ) { |
2706 | f2dHistEffMassLambdaVsEffMassXiMinus->Fill( lInvMassLambdaAsCascDghter, lInvMassXiMinus ); | |
2707 | f2dHistEffMassXiVsEffMassOmegaMinus ->Fill( lInvMassXiMinus, lInvMassOmegaMinus ); | |
c525f360 | 2708 | f2dHistXiRadiusVsEffMassXiMinus ->Fill( lXiRadius2D, lInvMassXiMinus ); |
2709 | f2dHistXiRadiusVsEffMassOmegaMinus ->Fill( lXiRadius2D, lInvMassOmegaMinus ); | |
ff0753e1 | 2710 | f3dHistXiPtVsEffMassVsYXiMinus ->Fill( lXiTransvMom, lInvMassXiMinus, lRapXi ); |
2711 | f3dHistXiPtVsEffMassVsYOmegaMinus ->Fill( lXiTransvMom, lInvMassOmegaMinus, lRapOmega ); | |
45ee0bcc | 2712 | } |
2713 | else{ | |
2714 | f2dHistEffMassLambdaVsEffMassXiPlus ->Fill( lInvMassLambdaAsCascDghter, lInvMassXiPlus ); | |
2715 | f2dHistEffMassXiVsEffMassOmegaPlus ->Fill( lInvMassXiPlus, lInvMassOmegaPlus ); | |
c525f360 | 2716 | f2dHistXiRadiusVsEffMassXiPlus ->Fill( lXiRadius2D, lInvMassXiPlus); |
2717 | f2dHistXiRadiusVsEffMassOmegaPlus ->Fill( lXiRadius2D, lInvMassOmegaPlus ); | |
ff0753e1 | 2718 | f3dHistXiPtVsEffMassVsYXiPlus ->Fill( lXiTransvMom, lInvMassXiPlus, lRapXi ); |
2719 | f3dHistXiPtVsEffMassVsYOmegaPlus ->Fill( lXiTransvMom, lInvMassOmegaPlus, lRapOmega ); | |
45ee0bcc | 2720 | } |
3873f6ee | 2721 | |
32599cd4 | 2722 | // - Filling the AliCFContainers related to PID |
2723 | ||
2724 | Double_t lContainerPIDVars[4] = {0.0}; | |
2725 | ||
2726 | // Xi Minus | |
2727 | if( lChargeXi < 0 ) { | |
2728 | lContainerPIDVars[0] = lXiTransvMom ; | |
2729 | lContainerPIDVars[1] = lInvMassXiMinus ; | |
2730 | lContainerPIDVars[2] = lRapXi ; | |
c525f360 | 2731 | lContainerPIDVars[3] = nTrackPrimaryMultiplicity; // FIXME : nTrackPrimaryMultiplicity not set for AOD ... = -1 |
32599cd4 | 2732 | |
2733 | // No PID | |
2734 | fCFContCascadePIDXiMinus->Fill(lContainerPIDVars, 0); // No PID | |
2735 | // TPC PID | |
2736 | if( lIsBachelorPionForTPC ) | |
e8fd19de | 2737 | fCFContCascadePIDXiMinus->Fill(lContainerPIDVars, 1); // TPC PID / 4-#sigma cut on Bachelor track |
32599cd4 | 2738 | |
2739 | if( lIsBachelorPionForTPC && | |
2740 | lIsPosProtonForTPC ) | |
e8fd19de | 2741 | fCFContCascadePIDXiMinus->Fill(lContainerPIDVars, 2); // TPC PID / 4-#sigma cut on Bachelor+Baryon tracks |
32599cd4 | 2742 | |
2743 | if( lIsBachelorPionForTPC && | |
2744 | lIsPosProtonForTPC && | |
2745 | lIsNegPionForTPC ) | |
e8fd19de | 2746 | fCFContCascadePIDXiMinus->Fill(lContainerPIDVars, 3); // TPC PID / 4-#sigma cut on Bachelor+Baryon+Meson tracks |
32599cd4 | 2747 | |
2748 | // Combined PID | |
2749 | if( lIsBachelorPion ) | |
2750 | fCFContCascadePIDXiMinus->Fill(lContainerPIDVars, 4); // Comb. PID / Bachelor | |
2751 | ||
2752 | if( lIsBachelorPion && | |
2753 | lIsPosInXiProton ) | |
2754 | fCFContCascadePIDXiMinus->Fill(lContainerPIDVars, 5); // Comb. PID / Bachelor+Baryon | |
2755 | ||
2756 | if(lIsBachelorPion && | |
2757 | lIsPosInXiProton && | |
2758 | lIsNegInXiPion ) | |
2759 | fCFContCascadePIDXiMinus->Fill(lContainerPIDVars, 6); // Comb. PID / Bachelor+Baryon+Meson | |
2760 | } | |
2761 | ||
2762 | lContainerPIDVars[0] = 0.; lContainerPIDVars[1] = 0.; lContainerPIDVars[2] = 0.; lContainerPIDVars[3] = 0.; | |
2763 | ||
2764 | // Xi Plus | |
2765 | if( lChargeXi > 0 ) { | |
2766 | lContainerPIDVars[0] = lXiTransvMom ; | |
2767 | lContainerPIDVars[1] = lInvMassXiPlus ; | |
2768 | lContainerPIDVars[2] = lRapXi ; | |
c525f360 | 2769 | lContainerPIDVars[3] = nTrackPrimaryMultiplicity; // FIXME : nTrackPrimaryMultiplicity not set for AOD ... = -1 |
32599cd4 | 2770 | |
2771 | // No PID | |
2772 | fCFContCascadePIDXiPlus->Fill(lContainerPIDVars, 0); // No PID | |
2773 | // TPC PID | |
2774 | if( lIsBachelorPionForTPC ) | |
e8fd19de | 2775 | fCFContCascadePIDXiPlus->Fill(lContainerPIDVars, 1); // TPC PID / 4-#sigma cut on Bachelor track |
32599cd4 | 2776 | |
2777 | if( lIsBachelorPionForTPC && | |
2778 | lIsNegProtonForTPC ) | |
e8fd19de | 2779 | fCFContCascadePIDXiPlus->Fill(lContainerPIDVars, 2); // TPC PID / 4-#sigma cut on Bachelor+Baryon tracks |
32599cd4 | 2780 | |
2781 | if( lIsBachelorPionForTPC && | |
2782 | lIsNegProtonForTPC && | |
2783 | lIsPosPionForTPC ) | |
e8fd19de | 2784 | fCFContCascadePIDXiPlus->Fill(lContainerPIDVars, 3); // TPC PID / 4-#sigma cut on Bachelor+Baryon+Meson tracks |
32599cd4 | 2785 | |
2786 | // Combined PID | |
2787 | if( lIsBachelorPion ) | |
2788 | fCFContCascadePIDXiPlus->Fill(lContainerPIDVars, 4); // Comb. PID / Bachelor | |
2789 | ||
2790 | if( lIsBachelorPion && | |
2791 | lIsNegInXiProton ) | |
2792 | fCFContCascadePIDXiPlus->Fill(lContainerPIDVars, 5); // Comb. PID / Bachelor+Baryon | |
2793 | ||
2794 | if(lIsBachelorPion && | |
2795 | lIsNegInXiProton && | |
2796 | lIsPosInXiPion ) | |
2797 | fCFContCascadePIDXiPlus->Fill(lContainerPIDVars, 6); // Comb. PID / Bachelor+Baryon+Meson | |
2798 | } | |
2799 | ||
2800 | lContainerPIDVars[0] = 0.; lContainerPIDVars[1] = 0.; lContainerPIDVars[2] = 0.; lContainerPIDVars[3] = 0.; | |
2801 | ||
2802 | // Omega Minus | |
a786bd57 | 2803 | if( lChargeXi < 0 && (TMath::Abs( lInvMassXiMinus-1.3217 ) > 0.008) ) { |
2804 | // In order to reduce the background contribution from Xi to Omega, | |
2805 | // --> trial to reject the neg. candidate whose mass under Xi hypothesis | |
2806 | // falls into the Xi window (Mpdg +- 8 MeV/c2) | |
32599cd4 | 2807 | lContainerPIDVars[0] = lXiTransvMom ; |
2808 | lContainerPIDVars[1] = lInvMassOmegaMinus ; | |
2809 | lContainerPIDVars[2] = lRapOmega ; | |
c525f360 | 2810 | lContainerPIDVars[3] = nTrackPrimaryMultiplicity; // FIXME : nTrackPrimaryMultiplicity not set for AOD ... = -1 |
32599cd4 | 2811 | |
2812 | // No PID | |
2813 | fCFContCascadePIDOmegaMinus->Fill(lContainerPIDVars, 0); // No PID | |
2814 | // TPC PID | |
2815 | if( lIsBachelorKaonForTPC ) | |
e8fd19de | 2816 | fCFContCascadePIDOmegaMinus->Fill(lContainerPIDVars, 1); // TPC PID / 4-#sigma cut on Bachelor track |
32599cd4 | 2817 | |
2818 | if( lIsBachelorKaonForTPC && | |
2819 | lIsPosProtonForTPC ) | |
e8fd19de | 2820 | fCFContCascadePIDOmegaMinus->Fill(lContainerPIDVars, 2); // TPC PID / 4-#sigma cut on Bachelor+Baryon tracks |
32599cd4 | 2821 | |
2822 | if( lIsBachelorKaonForTPC && | |
2823 | lIsPosProtonForTPC && | |
2824 | lIsNegPionForTPC ) | |
e8fd19de | 2825 | fCFContCascadePIDOmegaMinus->Fill(lContainerPIDVars, 3); // TPC PID / 4-#sigma cut on Bachelor+Baryon+Meson tracks |
32599cd4 | 2826 | |
2827 | // Combined PID | |
2828 | if( lIsBachelorKaon ) | |
2829 | fCFContCascadePIDOmegaMinus->Fill(lContainerPIDVars, 4); // Comb. PID / Bachelor | |
2830 | ||
2831 | if( lIsBachelorKaon && | |
2832 | lIsPosInOmegaProton ) | |
2833 | fCFContCascadePIDOmegaMinus->Fill(lContainerPIDVars, 5); // Comb. PID / Bachelor+Baryon | |
2834 | ||
2835 | if(lIsBachelorKaon && | |
2836 | lIsPosInOmegaProton && | |
2837 | lIsNegInOmegaPion ) | |
2838 | fCFContCascadePIDOmegaMinus->Fill(lContainerPIDVars, 6); // Comb. PID / Bachelor+Baryon+Meson | |
2839 | } | |
2840 | ||
2841 | lContainerPIDVars[0] = 0.; lContainerPIDVars[1] = 0.; lContainerPIDVars[2] = 0.; lContainerPIDVars[3] = 0.; | |
2842 | ||
2843 | // Omega Plus | |
a786bd57 | 2844 | if( lChargeXi > 0 && (TMath::Abs( lInvMassXiPlus-1.3217 ) > 0.008) ) { |
2845 | // In order to reduce the background contribution from Xi to Omega, | |
2846 | // --> trial to reject the neg. candidate whose mass under Xi hypothesis | |
2847 | // falls into the Xi window (Mpdg +- 8 MeV/c2) | |
32599cd4 | 2848 | lContainerPIDVars[0] = lXiTransvMom ; |
2849 | lContainerPIDVars[1] = lInvMassOmegaPlus ; | |
2850 | lContainerPIDVars[2] = lRapOmega ; | |
c525f360 | 2851 | lContainerPIDVars[3] = nTrackPrimaryMultiplicity; // FIXME : nTrackPrimaryMultiplicity not set for AOD ... = -1 |
32599cd4 | 2852 | |
2853 | // No PID | |
2854 | fCFContCascadePIDOmegaPlus->Fill(lContainerPIDVars, 0); // No PID | |
2855 | // TPC PID | |
2856 | if( lIsBachelorKaonForTPC ) | |
e8fd19de | 2857 | fCFContCascadePIDOmegaPlus->Fill(lContainerPIDVars, 1); // TPC PID / 4-#sigma cut on Bachelor track |
32599cd4 | 2858 | |
2859 | if( lIsBachelorKaonForTPC && | |
2860 | lIsNegProtonForTPC ) | |
e8fd19de | 2861 | fCFContCascadePIDOmegaPlus->Fill(lContainerPIDVars, 2); // TPC PID / 4-#sigma cut on Bachelor+Baryon tracks |
32599cd4 | 2862 | |
2863 | if( lIsBachelorKaonForTPC && | |
2864 | lIsNegProtonForTPC && | |
2865 | lIsPosPionForTPC ) | |
e8fd19de | 2866 | fCFContCascadePIDOmegaPlus->Fill(lContainerPIDVars, 3); // TPC PID / 4-#sigma cut on Bachelor+Baryon+Meson tracks |
32599cd4 | 2867 | |
2868 | // Combined PID | |
2869 | if( lIsBachelorKaon ) | |
2870 | fCFContCascadePIDOmegaPlus->Fill(lContainerPIDVars, 4); // Comb. PID / Bachelor | |
2871 | ||
2872 | if( lIsBachelorKaon && | |
2873 | lIsNegInOmegaProton ) | |
2874 | fCFContCascadePIDOmegaPlus->Fill(lContainerPIDVars, 5); // Comb. PID / Bachelor+Baryon | |
2875 | ||
2876 | if(lIsBachelorKaon && | |
2877 | lIsNegInOmegaProton && | |
2878 | lIsPosInOmegaPion ) | |
2879 | fCFContCascadePIDOmegaPlus->Fill(lContainerPIDVars, 6); // Comb. PID / Bachelor+Baryon+Meson | |
2880 | } | |
2881 | ||
e8fd19de | 2882 | |
2883 | // II.Fill.Step 7 : filling the AliCFContainer (optimisation of topological selections) | |
2884 | Double_t lContainerCutVars[20] = {0.0}; | |
2885 | ||
2886 | lContainerCutVars[0] = lDcaXiDaughters; | |
2887 | lContainerCutVars[1] = lDcaBachToPrimVertexXi; | |
2888 | lContainerCutVars[2] = lXiCosineOfPointingAngle; | |
c525f360 | 2889 | lContainerCutVars[3] = lXiRadius2D; |
e8fd19de | 2890 | lContainerCutVars[4] = lInvMassLambdaAsCascDghter; |
2891 | lContainerCutVars[5] = lDcaV0DaughtersXi; | |
2892 | lContainerCutVars[6] = lV0CosineOfPointingAngleXi; | |
2893 | lContainerCutVars[7] = lV0RadiusXi; | |
2894 | lContainerCutVars[8] = lDcaV0ToPrimVertexXi; | |
2895 | lContainerCutVars[9] = lDcaPosToPrimVertexXi; | |
2896 | lContainerCutVars[10] = lDcaNegToPrimVertexXi; | |
2897 | ||
2898 | lContainerCutVars[13] = lXiTransvMom; | |
2899 | ||
2900 | lContainerCutVars[16] = lBestPrimaryVtxPos[2]; | |
a786bd57 | 2901 | lContainerCutVars[17] = nTrackPrimaryMultiplicity; // FIXME : nTrackPrimaryMultiplicity not checked for AOD ... = -1 |
2902 | lContainerCutVars[18] = nITSandTPCtracksAndSPDtracklets; // FIXME : nITSandTPCtracksAndSPDtracklets is not available for AOD ... = -1 | |
e8fd19de | 2903 | lContainerCutVars[19] = lBachTPCClusters; // FIXME : BachTPCClusters is not available for AOD ... = -1 |
2904 | ||
2905 | if( lChargeXi < 0 ) { | |
9d670198 | 2906 | lContainerCutVars[11] = lInvMassXiMinus; |
e8fd19de | 2907 | lContainerCutVars[12] = 1.63; |
9d670198 | 2908 | lContainerCutVars[14] = lRapXi; |
2909 | lContainerCutVars[15] = -1.; | |
e8fd19de | 2910 | if( lIsBachelorPionForTPC && lIsPosProtonForTPC && lIsNegPionForTPC ) fCFContCascadeCuts->Fill(lContainerCutVars,0); // for Xi- |
9d670198 | 2911 | |
e8fd19de | 2912 | lContainerCutVars[11] = 1.26; |
9d670198 | 2913 | lContainerCutVars[12] = lInvMassOmegaMinus; |
2914 | lContainerCutVars[14] = -1.; | |
2915 | lContainerCutVars[15] = lRapOmega; | |
a786bd57 | 2916 | if( lIsBachelorKaonForTPC && lIsPosProtonForTPC && lIsNegPionForTPC |
2917 | && (TMath::Abs( lInvMassXiMinus-1.3217 ) > 0.008) ) fCFContCascadeCuts->Fill(lContainerCutVars,2); // for Omega- | |
9d670198 | 2918 | } |
e8fd19de | 2919 | else{ |
9d670198 | 2920 | lContainerCutVars[11] = lInvMassXiPlus; |
e8fd19de | 2921 | lContainerCutVars[12] = 1.63; |
9d670198 | 2922 | lContainerCutVars[14] = lRapXi; |
2923 | lContainerCutVars[15] = -1.; | |
e8fd19de | 2924 | if( lIsBachelorPionForTPC && lIsNegProtonForTPC && lIsPosPionForTPC ) fCFContCascadeCuts->Fill(lContainerCutVars,1); // for Xi+ |
9d670198 | 2925 | |
e8fd19de | 2926 | lContainerCutVars[11] = 1.26; |
9d670198 | 2927 | lContainerCutVars[12] = lInvMassOmegaPlus; |
2928 | lContainerCutVars[14] = -1.; | |
2929 | lContainerCutVars[15] = lRapOmega; | |
a786bd57 | 2930 | if( lIsBachelorKaonForTPC && lIsNegProtonForTPC && lIsPosPionForTPC |
2931 | && (TMath::Abs( lInvMassXiPlus-1.3217 ) > 0.008) ) fCFContCascadeCuts->Fill(lContainerCutVars,3); // for Omega+ | |
9d670198 | 2932 | } |
e8fd19de | 2933 | |
2934 | ||
2935 | // II.Fill.Step 8 : angular correlations | |
2936 | ||
2937 | if( lChargeXi < 0 ){ | |
c525f360 | 2938 | if( lIsBachelorPionForTPC && lIsPosProtonForTPC && lIsNegPionForTPC ) DoAngularCorrelation("Xi-", lInvMassXiMinus, lArrTrackID, lTVect3MomXi, lEta, lRapXi); |
a786bd57 | 2939 | if( lIsBachelorKaonForTPC && lIsPosProtonForTPC && lIsNegPionForTPC |
2940 | && (TMath::Abs( lInvMassXiMinus-1.3217 ) > 0.008) ) DoAngularCorrelation("Omega-", lInvMassOmegaMinus, lArrTrackID, lTVect3MomXi, lEta, lRapOmega); | |
e8fd19de | 2941 | } |
2942 | else{ | |
c525f360 | 2943 | if( lIsBachelorPionForTPC && lIsNegProtonForTPC && lIsPosPionForTPC ) DoAngularCorrelation("Xi+", lInvMassXiPlus, lArrTrackID, lTVect3MomXi, lEta, lRapXi); |
a786bd57 | 2944 | if( lIsBachelorKaonForTPC && lIsNegProtonForTPC && lIsPosPionForTPC |
2945 | && (TMath::Abs( lInvMassXiPlus-1.3217 ) > 0.008) ) DoAngularCorrelation("Omega+", lInvMassOmegaPlus, lArrTrackID, lTVect3MomXi, lEta, lRapOmega); | |
e8fd19de | 2946 | } |
2947 | ||
3873f6ee | 2948 | |
45ee0bcc | 2949 | }// end of the Cascade loop (ESD or AOD) |
2950 | ||
2951 | ||
2952 | // Post output data. | |
2953 | PostData(1, fListHistCascade); | |
2954 | } | |
2955 | ||
2956 | ||
3873f6ee | 2957 | void AliAnalysisTaskCheckCascade::DoAngularCorrelation( const Char_t *lCascType, |
2958 | Double_t lInvMassCascade, | |
2959 | const Int_t *lArrTrackID, | |
2960 | TVector3 &lTVect3MomXi, | |
c525f360 | 2961 | Double_t lEtaXi, |
2962 | Double_t lRapCascade){ | |
3aedd4a5 | 2963 | // Perform the Delta(Phi)Delta(Eta) analysis |
2964 | // by properly filling the THnSparseF | |
c525f360 | 2965 | |
2966 | if( fAnalysisType == "AOD") return; // FIXME : AOD development lost, because of AliESDtrack needed by AliESDtrackCuts | |
2967 | ||
2968 | TString lStrCascType( lCascType ); | |
2969 | ||
2970 | // Check the Xi- candidate is within the proper mass window m0 +- 10 MeV | |
2971 | Double_t lCascPdgMass = 0.0; | |
2972 | if( lStrCascType.Contains("Xi") ) lCascPdgMass = 1.3217; | |
2973 | if( lStrCascType.Contains("Omega") ) lCascPdgMass = 1.6724; | |
2974 | ||
2975 | if( lInvMassCascade > lCascPdgMass + 0.010) return; | |
2976 | if( lInvMassCascade < lCascPdgMass - 0.010) return; | |
2977 | ||
2978 | // Check the Xi- candidate is within the proper rapidity window (flat efficiency) | |
2979 | if( lRapCascade > 0.5 ) return; | |
2980 | if( lRapCascade < -0.5 ) return; | |
2981 | ||
2982 | ||
2983 | // 3 options to follow for the correlations: | |
2984 | // 1.1 - "TrigAnyCasc-AssoAnyPrim" | |
2985 | // 1.2 - "TrigCascLeading-AssoAnyPrim" | |
2986 | // 2. - "TrigLeadingTrck-AssoCasc" | |
2987 | ||
2988 | if(fAngularCorrelationType.Contains("AssoAnyPrim") ){ | |
2989 | //----------------------- Option 1 --------------------------------------------------------------------------------------------------- | |
a786bd57 | 2990 | // Trigger = cascade, |
c525f360 | 2991 | // Associated = all the primary tracks in the event. |
2992 | // 1.1 Cascade = trigger above a certain pt but nothing more complicated than that | |
2993 | // 1.2 Cascade = leading particle (Ã la heavy-ion -> test for coalescence) | |
2994 | ||
2995 | Bool_t kRejectLowPtCascades = kTRUE; | |
2996 | if(kRejectLowPtCascades && (lTVect3MomXi.Pt() < 1.7) ) return; | |
2997 | // Do not even consider the cascade of low pt for the correlation ... | |
2998 | ||
a786bd57 | 2999 | if(fAngularCorrelationType == "TrigCascLeading-AssoAnyPrim"){// Require the Cascade to be the Leading Part. in the event |
c525f360 | 3000 | |
3001 | // 1st loop: check there is no primary track with a higher pt ... | |
3002 | // = The cascade is meant to be a leading particle : Pt(Casc) > any primary track in the event | |
3003 | for(Int_t TrckIdx = 0; TrckIdx < (InputEvent())->GetNumberOfTracks() ; TrckIdx++ ) | |
3004 | {// Loop over all the tracks of the event | |
3005 | ||
3006 | AliESDtrack *lCurrentTrck = dynamic_cast<AliESDtrack*>( (InputEvent())->GetTrack( TrckIdx ) ); | |
3007 | if (!lCurrentTrck ) { | |
3008 | AliWarning("ERROR Correl. Study : Could not retrieve a track while looping over the event tracks ..."); | |
3009 | continue; | |
3010 | } | |
3011 | ||
3012 | if( !fESDtrackCuts->AcceptTrack(lCurrentTrck) ) continue; | |
3013 | // Just consider primary tracks (= reject track that are not primary ones) | |
3014 | if(lTVect3MomXi.Pt() < lCurrentTrck->Pt() ) return; | |
3015 | // Room for improvement: //FIXME | |
3016 | // 1. There is a given resolution on pt : maybe release the cut Pt(casc) < Pt(track)*90% ? | |
3017 | // 2. Apply this cut only when DeltaPhi(casc, track) > 90 deg = when track is in the near-side ? | |
3018 | // 3. Anti-splitting cut (like in Femto analysis) ? = now done via ESDtrackCuts ... | |
3019 | ||
3020 | }// end control loop | |
3021 | }// end of prelim. check : Cascade = leading part in the event | |
3022 | ||
3023 | // 2nd loop: filling loop | |
3024 | for(Int_t TrckIdx = 0; TrckIdx < (InputEvent())->GetNumberOfTracks() ; TrckIdx++ ) | |
3025 | {// Loop over all the tracks of the event | |
3026 | ||
3027 | AliESDtrack *lCurrentTrck = dynamic_cast<AliESDtrack*>( (InputEvent())->GetTrack( TrckIdx ) ); | |
3028 | if (!lCurrentTrck ) { | |
3029 | AliWarning("ERROR Correl. Study : Could not retrieve a track while looping over the event tracks ..."); | |
3030 | continue; | |
3031 | } | |
3032 | // Just consider primary tracks (= reject track that are not primary ones) | |
3033 | if( !fESDtrackCuts->AcceptTrack(lCurrentTrck) ) continue; | |
3034 | ||
3035 | // Room for improvement: //FIXME | |
3036 | // 1. | |
3037 | // 2. Exclude the tracks that build the condisdered cascade = the bachelor + the V0 dghters | |
3038 | // This may bias the outcome, especially for low multplicity events. | |
3039 | // Note : For ESD event, track ID == track index. | |
3040 | if(lCurrentTrck->GetID() == lArrTrackID[0]) continue; | |
3041 | if(lCurrentTrck->GetID() == lArrTrackID[1]) continue; | |
3042 | if(lCurrentTrck->GetID() == lArrTrackID[2]) continue; | |
3043 | ||
3044 | TVector3 lTVect3MomTrck(lCurrentTrck->Px(), lCurrentTrck->Py(), lCurrentTrck->Pz() ); | |
3045 | ||
3046 | // 2 hypotheses made here : | |
3047 | // - The Xi trajectory is a straight line, | |
3048 | // - The Xi doesn't loose any energy by crossing the first layer(s) of ITS, if ever; | |
3049 | // So, meaning hyp: vect p(Xi) at the emission = vect p(Xi) at the decay vertex | |
3050 | // By doing this, we introduce a systematic error on the cascade Phi ... | |
3051 | // Room for improvement: take into account the curvature of the Xi trajectory ? | |
3052 | // or rather, the resolution in space of the decay vertex ... | |
3053 | //FIXME | |
3054 | ||
3055 | Double_t lHnSpFillVar[5] = {0.}; | |
3056 | lHnSpFillVar[0] = lTVect3MomXi.DeltaPhi(lTVect3MomTrck) * 180.0/TMath::Pi(); // Delta phi(Casc,Track) (deg) | |
3057 | if(lHnSpFillVar[0] < -50.0) lHnSpFillVar[0] += 360.0; | |
3058 | lHnSpFillVar[1] = lEtaXi - lCurrentTrck->Eta(); // Delta eta(Casc,Track) | |
3059 | lHnSpFillVar[2] = lTVect3MomXi.Pt(); // Pt_{Casc} | |
3060 | lHnSpFillVar[3] = lCurrentTrck->Pt(); // Pt_{any track} | |
3061 | lHnSpFillVar[4] = lInvMassCascade; // Eff. Inv Mass (control var) | |
3062 | ||
3063 | if( lStrCascType.Contains("Xi-") ) fHnSpAngularCorrXiMinus ->Fill( lHnSpFillVar ); | |
3064 | else if( lStrCascType.Contains("Xi+") ) fHnSpAngularCorrXiPlus ->Fill( lHnSpFillVar ); | |
3065 | else if( lStrCascType.Contains("Omega-") ) fHnSpAngularCorrOmegaMinus ->Fill( lHnSpFillVar ); | |
3066 | else if( lStrCascType.Contains("Omega+") ) fHnSpAngularCorrOmegaPlus ->Fill( lHnSpFillVar ); | |
3067 | ||
3068 | }// end - Loop over all the tracks in the event | |
3069 | ||
3070 | }// end of correlation type : "Trig = Casc - Asso = AnyPrim", (cases 1.1 and 1.2) | |
3071 | ||
3072 | ||
3073 | ||
3074 | else if(fAngularCorrelationType == "TrigLeadingTrck-AssoCasc"){ | |
3075 | ||
3076 | //----------------------- Option 2 --------------------------------------------------------------------------------------------------- | |
a786bd57 | 3077 | // Trigger = leading track, |
c525f360 | 3078 | // Associated = the cascade'S' in the event |
3079 | // NOTE : several good cascades could be present in the event (e.g. one leading track as trigger, 2 associated Xi) ... | |
3080 | // The present function will then be called several times. | |
3081 | // = issue for the normalisation ... | |
3082 | ||
3083 | // 1st loop: | |
3084 | // find the index of the (1) primary track (2) which is the leading particle in pt | |
3085 | // NOTE : we do not take into account the Cascade pt, i.e. pt(Casc) could be greater or lower than pt(Leading) ... | |
3086 | Int_t lLeadingPrimTrackIdx = -1; | |
3087 | Double_t lPtMax = 0.1; | |
3088 | ||
3089 | for(Int_t TrckIdx = 0; TrckIdx < (InputEvent())->GetNumberOfTracks() ; TrckIdx++ ) | |
3090 | {// Loop over all the tracks of the event | |
3091 | ||
3092 | AliESDtrack *lCurrentTrck = dynamic_cast<AliESDtrack*>( (InputEvent())->GetTrack( TrckIdx ) ); | |
3093 | if (!lCurrentTrck ) { | |
3094 | AliWarning("ERROR Correl. Study : Could not retrieve a track while looping over the event tracks ..."); | |
3095 | continue; | |
3096 | } | |
3097 | ||
3098 | // Primary track selection | |
3099 | if( !fESDtrackCuts->AcceptTrack(lCurrentTrck) ) continue; | |
3100 | ||
3101 | // Exclude the tracks that build the condisdered cascade = the bachelor + the V0 dghters | |
3102 | // This may bias the outcome, especially for low multplicity events. | |
3103 | // Note : For ESD event, track ID == track index. | |
3104 | if(lCurrentTrck->GetID() == lArrTrackID[0]) continue; | |
3105 | if(lCurrentTrck->GetID() == lArrTrackID[1]) continue; | |
3106 | if(lCurrentTrck->GetID() == lArrTrackID[2]) continue; | |
3107 | ||
3108 | // Towards the leading track | |
3109 | if( lPtMax < lCurrentTrck->Pt() ){ | |
3110 | lLeadingPrimTrackIdx = TMath::Abs( lCurrentTrck->GetID() ); | |
3111 | lPtMax = lCurrentTrck->Pt(); | |
3112 | } | |
3113 | }// end leading track finding loop | |
3114 | ||
3115 | if( lLeadingPrimTrackIdx < 0 ) return; | |
3116 | if( lPtMax < 0.101 ) return; | |
3117 | ||
3118 | ||
3119 | // 2nd step: filling ONCE the THnSparse | |
14f5660a | 3120 | AliESDtrack *lLeadingTrck = 0x0; |
3121 | lLeadingTrck =dynamic_cast<AliESDtrack*>( (InputEvent())->GetTrack( lLeadingPrimTrackIdx ) ); | |
3122 | if( !lLeadingTrck ) return; | |
c525f360 | 3123 | |
3124 | TVector3 lTVect3MomLeadingTrck( lLeadingTrck->Px(), lLeadingTrck->Py(), lLeadingTrck->Pz() ); | |
3125 | ||
3126 | // 2 hypotheses made here : | |
3127 | // - The Xi trajectory is a straight line, | |
3128 | // - The Xi doesn't loose any energy by crossing the first layer(s) of ITS, if ever; | |
3129 | // So, meaning hyp: vect p(Xi) at the emission = vect p(Xi) at the decay vertex | |
3130 | // By doing this, we introduce a systematic error on the cascade Phi ... | |
3131 | // Room for improvement: take into account the curvature of the Xi trajectory ? | |
3132 | // or rather, the resolution in space of the decay vertex ... | |
3133 | //FIXME | |
3134 | ||
3135 | Double_t lHnSpFillVar[5] = {0.}; | |
3136 | lHnSpFillVar[0] = lTVect3MomLeadingTrck.DeltaPhi(lTVect3MomXi) * 180.0/TMath::Pi(); // Delta phi(leading Track, Casc) (deg) | |
3137 | if(lHnSpFillVar[0] < -50.0) lHnSpFillVar[0] += 360.0; | |
3138 | lHnSpFillVar[1] = lLeadingTrck->Eta() - lEtaXi; // Delta eta(leading Track, Casc) | |
3139 | lHnSpFillVar[2] = lTVect3MomXi.Pt(); // Pt_{Casc} | |
3140 | lHnSpFillVar[3] = lLeadingTrck->Pt(); // Pt_{leading track} | |
3141 | lHnSpFillVar[4] = lInvMassCascade; // Eff. Inv Mass (control var) | |
3142 | ||
3143 | if( lStrCascType.Contains("Xi-") ) fHnSpAngularCorrXiMinus ->Fill( lHnSpFillVar ); | |
3144 | else if( lStrCascType.Contains("Xi+") ) fHnSpAngularCorrXiPlus ->Fill( lHnSpFillVar ); | |
3145 | else if( lStrCascType.Contains("Omega-") ) fHnSpAngularCorrOmegaMinus ->Fill( lHnSpFillVar ); | |
3146 | else if( lStrCascType.Contains("Omega+") ) fHnSpAngularCorrOmegaPlus ->Fill( lHnSpFillVar ); | |
3147 | }// end of correlation type : "Trig = LeadingTrck -Asso = Casc" | |
3148 | else | |
3149 | return; | |
45ee0bcc | 3150 | |
3873f6ee | 3151 | } |
45ee0bcc | 3152 | |
7e504402 | 3153 | Int_t AliAnalysisTaskCheckCascade::DoESDTrackWithTPCrefitMultiplicity(const AliESDEvent *lESDevent) |
3154 | { | |
3155 | // Checking the number of tracks with TPCrefit for each event | |
3156 | // Needed for a rough assessment of the event multiplicity | |
3157 | ||
3158 | Int_t nTrackWithTPCrefitMultiplicity = 0; | |
3159 | for(Int_t iTrackIdx = 0; iTrackIdx < (InputEvent())->GetNumberOfTracks(); iTrackIdx++){ | |
3160 | AliESDtrack *esdTrack = 0x0; | |
3161 | esdTrack = lESDevent->GetTrack( iTrackIdx ); | |
3162 | if (!esdTrack) { AliWarning("Pb / Could not retrieve one track within the track loop for TPCrefit check ..."); continue; } | |
3163 | ||
3164 | ULong_t lTrackStatus = esdTrack->GetStatus(); | |
3165 | if ((lTrackStatus&AliESDtrack::kTPCrefit) == 0) continue; | |
3166 | else nTrackWithTPCrefitMultiplicity++; | |
3167 | // FIXME : | |
3168 | // The goal here is to get a better assessment of the event multiplicity. | |
3169 | // (InputEvent())->GetNumberOfTracks() takes into account ITS std alone tracks + global tracks | |
3170 | // This may introduce a bias. Hence the number of TPC refit tracks. | |
3171 | // Note : the event multiplicity = analysis on its own... See Jacek's or Jan Fiete's analysis on dN/d(eta) | |
45ee0bcc | 3172 | |
7e504402 | 3173 | }// end loop over all event tracks |
3174 | return nTrackWithTPCrefitMultiplicity; | |
3175 | } | |
45ee0bcc | 3176 | |
3177 | ||
3178 | ||
3179 | ||
3180 | ||
3181 | //________________________________________________________________________ | |
3182 | void AliAnalysisTaskCheckCascade::Terminate(Option_t *) | |
3183 | { | |
3184 | // Draw result to the screen | |
3185 | // Called once at the end of the query | |
3186 | ||
cd9188c9 | 3187 | TList *cRetrievedList = 0x0; |
3188 | cRetrievedList = (TList*)GetOutputData(1); | |
3189 | if(!cRetrievedList){ | |
7e504402 | 3190 | AliWarning("ERROR - AliAnalysisTaskCheckCascade: ouput data container list not available\n"); return; |
cd9188c9 | 3191 | } |
3192 | ||
9d670198 | 3193 | fHistTrackMultiplicityForTrigEvt = dynamic_cast<TH1F*> ( cRetrievedList->FindObject("fHistTrackMultiplicityForTrigEvt") ); |
3194 | if (!fHistTrackMultiplicityForTrigEvt) { | |
3195 | AliWarning("ERROR - AliAnalysisTaskCheckCascade: fHistTrackMultiplicityForTrigEvt not available\n"); return; | |
cd9188c9 | 3196 | } |
45ee0bcc | 3197 | |
9d670198 | 3198 | fHistCascadeMultiplicityForTrigEvt = dynamic_cast<TH1F*> ( cRetrievedList->FindObject("fHistCascadeMultiplicityForTrigEvt") ); |
3199 | if (!fHistCascadeMultiplicityForTrigEvt) { | |
3200 | AliWarning("ERROR - AliAnalysisTaskCheckCascade: fHistCascadeMultiplicityForTrigEvt not available\n"); return; | |
cd9188c9 | 3201 | } |
32599cd4 | 3202 | |
3203 | fHistMassXiMinus = dynamic_cast<TH1F*> ( cRetrievedList->FindObject("fHistMassXiMinus") ); | |
3204 | if (!fHistMassXiMinus) { | |
7e504402 | 3205 | AliWarning("ERROR - AliAnalysisTaskCheckCascade: fHistMassXiMinus not available\n"); return; |
32599cd4 | 3206 | } |
3207 | fHistMassXiPlus = dynamic_cast<TH1F*> ( cRetrievedList->FindObject("fHistMassXiPlus") ); | |
3208 | if (!fHistMassXiPlus) { | |
7e504402 | 3209 | AliWarning("ERROR - AliAnalysisTaskCheckCascade: fHistMassXiPlus not available\n"); return; |
32599cd4 | 3210 | } |
3211 | fHistMassOmegaMinus = dynamic_cast<TH1F*> ( cRetrievedList->FindObject("fHistMassOmegaMinus") ); | |
3212 | if (!fHistMassOmegaMinus) { | |
7e504402 | 3213 | AliWarning("ERROR - AliAnalysisTaskCheckCascade: fHistMassOmegaMinus not available\n"); return; |
32599cd4 | 3214 | } |
3215 | fHistMassOmegaPlus = dynamic_cast<TH1F*> ( cRetrievedList->FindObject("fHistMassOmegaPlus") ); | |
3216 | if (!fHistMassOmegaPlus) { | |
7e504402 | 3217 | AliWarning("ERROR - AliAnalysisTaskCheckCascade: fHistMassOmegaPlus not available\n"); return; |
32599cd4 | 3218 | } |
3219 | ||
f3dc9369 | 3220 | TCanvas *canCheckCascade = new TCanvas("AliAnalysisTaskCheckCascade","CheckCascade overview",10,10,1010,660); |
32599cd4 | 3221 | canCheckCascade->Divide(2,2); |
3222 | ||
3223 | canCheckCascade->cd(1); | |
cd9188c9 | 3224 | canCheckCascade->cd(1)->SetLogy(); |
9d670198 | 3225 | fHistTrackMultiplicityForTrigEvt->SetMarkerStyle(kFullStar); |
3226 | fHistTrackMultiplicityForTrigEvt->GetXaxis()->SetLabelFont(42); | |
3227 | fHistTrackMultiplicityForTrigEvt->GetYaxis()->SetLabelFont(42); | |
3228 | fHistTrackMultiplicityForTrigEvt->SetTitleFont(42, "xy"); | |
3229 | fHistTrackMultiplicityForTrigEvt->GetXaxis()->SetTitleOffset(1.1); | |
3230 | fHistTrackMultiplicityForTrigEvt->DrawCopy("H"); | |
45ee0bcc | 3231 | |
32599cd4 | 3232 | canCheckCascade->cd(2); |
3233 | canCheckCascade->cd(2)->SetLogy(); | |
9d670198 | 3234 | fHistCascadeMultiplicityForTrigEvt->SetMarkerStyle(kOpenSquare); |
3235 | fHistCascadeMultiplicityForTrigEvt->GetXaxis()->SetLabelFont(42); | |
3236 | fHistCascadeMultiplicityForTrigEvt->GetYaxis()->SetLabelFont(42); | |
3237 | fHistCascadeMultiplicityForTrigEvt->SetTitleFont(42, "xy"); | |
3238 | fHistCascadeMultiplicityForTrigEvt->GetXaxis()->SetTitleOffset(1.1); | |
3239 | fHistCascadeMultiplicityForTrigEvt->DrawCopy("E"); | |
45ee0bcc | 3240 | |
32599cd4 | 3241 | canCheckCascade->cd(3); |
f3dc9369 | 3242 | fHistMassXiMinus ->SetMarkerStyle(kFullCircle); |
32599cd4 | 3243 | fHistMassXiMinus ->SetMarkerSize(0.5); |
f3dc9369 | 3244 | fHistMassXiMinus ->GetXaxis()->SetLabelFont(42); |
3245 | fHistMassXiMinus ->GetYaxis()->SetLabelFont(42); | |
3246 | fHistMassXiMinus ->SetTitleFont(42, "xy"); | |
3247 | fHistMassXiMinus ->GetXaxis()->SetTitleOffset(1.1); | |
3248 | fHistMassXiMinus ->GetYaxis()->SetTitleOffset(1.3); | |
3249 | // fHistMassXiMinus->Rebin(2); | |
3250 | fHistMassXiMinus ->GetXaxis()->SetRangeUser(1.24, 1.42); | |
3251 | fHistMassXiMinus ->DrawCopy("E"); | |
3252 | ||
32599cd4 | 3253 | fHistMassXiPlus ->SetMarkerStyle(kOpenCircle); |
f3dc9369 | 3254 | fHistMassXiPlus ->SetMarkerColor(kRed+2); |
3255 | fHistMassXiPlus ->SetLineColor(kRed+2); | |
32599cd4 | 3256 | fHistMassXiPlus ->SetMarkerSize(0.5); |
f3dc9369 | 3257 | // fHistMassXiPlus ->Rebin(2); |
32599cd4 | 3258 | fHistMassXiPlus ->DrawCopy("ESAME"); |
3259 | ||
f3dc9369 | 3260 | |
3261 | TLegend *legendeXi =new TLegend(0.67,0.34,0.97,0.54); | |
3262 | legendeXi->SetTextFont(42); | |
3263 | legendeXi->SetTextSize(0.05); | |
3264 | legendeXi->SetFillColor(kWhite); | |
3265 | legendeXi->AddEntry( fHistMassXiMinus,"#Xi^{-} candidates","lp"); | |
a786bd57 | 3266 | legendeXi->AddEntry( fHistMassXiPlus,"#bar{#Xi}^{+} candidates","lp"); |
f3dc9369 | 3267 | legendeXi->Draw(); |
3268 | ||
3269 | ||
32599cd4 | 3270 | canCheckCascade->cd(4); |
32599cd4 | 3271 | fHistMassOmegaPlus ->SetMarkerStyle(kOpenCircle); |
f3dc9369 | 3272 | fHistMassOmegaPlus ->SetMarkerColor(kRed+2); |
3273 | fHistMassOmegaPlus ->SetLineColor(kRed+2); | |
32599cd4 | 3274 | fHistMassOmegaPlus ->SetMarkerSize(0.5); |
f3dc9369 | 3275 | fHistMassOmegaPlus ->GetXaxis()->SetLabelFont(42); |
3276 | fHistMassOmegaPlus ->GetYaxis()->SetLabelFont(42); | |
3277 | fHistMassOmegaPlus ->SetTitleFont(42, "xy"); | |
3278 | fHistMassOmegaPlus ->GetXaxis()->SetTitleOffset(1.1); | |
3279 | fHistMassOmegaPlus ->GetYaxis()->SetTitleOffset(1.25); | |
3280 | // fHistMassOmegaPlus ->Rebin(2); | |
3281 | fHistMassOmegaPlus ->GetXaxis()->SetRangeUser(1.6, 1.84); | |
3282 | fHistMassOmegaPlus ->DrawCopy("E"); | |
3283 | ||
3284 | fHistMassOmegaMinus->SetMarkerStyle(kFullCircle); | |
3285 | fHistMassOmegaMinus->SetMarkerSize(0.5); | |
3286 | // fHistMassOmegaMinus->Rebin(2); | |
3287 | fHistMassOmegaMinus->DrawCopy("ESAME"); | |
3288 | ||
3289 | ||
3290 | TLegend *legendeOmega = new TLegend(0.67,0.34,0.97,0.54); | |
3291 | legendeOmega->SetTextFont(42); | |
3292 | legendeOmega->SetTextSize(0.05); | |
3293 | legendeOmega->SetFillColor(kWhite); | |
3294 | legendeOmega->AddEntry( fHistMassOmegaMinus,"#Omega^{-} candidates","lp"); | |
a786bd57 | 3295 | legendeOmega->AddEntry( fHistMassOmegaPlus,"#bar{#Omega}^{+} candidates","lp"); |
f3dc9369 | 3296 | legendeOmega->Draw(); |
32599cd4 | 3297 | |
45ee0bcc | 3298 | } |