Skip to content

Commit f9d63f8

Browse files
author
Lucia Anna Tarasovicova
committed
fix of selection, numbering and new QA histograms
1 parent 8c2a0ac commit f9d63f8

File tree

1 file changed

+22
-19
lines changed

1 file changed

+22
-19
lines changed

‎ALICE3/Tasks/alice3-strangeness.cxx‎

Lines changed: 22 additions & 19 deletions
Original file line numberDiff line numberDiff line change
@@ -77,7 +77,7 @@ struct alice3strangeness {
7777
ConfigurableAxis axisDCAV0Daughters{"axisDCAV0Daughters", {20, 0, 5}, "DCA V0 daughters"};
7878
ConfigurableAxis axisPointingAngle{"axisPointingAngle", {40, 0.0f, 0.4f}, "pointing angle "};
7979
ConfigurableAxis axisProperLifeTime{"axisProperLifeTime", {100, 0.0f, 100.0f}, "proper lifetime (cm)"};
80-
80+
ConfigurableAxis axisEta{"axisEta", {100, -5.0f, 5.0f}, "eta"};
8181
} histAxes;
8282

8383
struct : ConfigurableGroup {
@@ -124,8 +124,9 @@ struct alice3strangeness {
124124
histos.add("K0/hDCAPosDaughter", "", kTH1D, {{200, -5, 5}});
125125
histos.add("hPVz", "hPVz", kTH1F, {histAxes.axisVertexZ});
126126
histos.add("hV0CandidateCounter", "hV0CandidateCounter", kTH1F, {{11, 0, 11}});
127-
histos.add("reconstructedCandidates/K0/hMass", "hMass", kTH2D, {histAxes.axisK0Mass, histAxes.axisPt});
128-
histos.add("reconstructedCandidates/Lambda/hMass", "hMass", kTH2D, {histAxes.axisLambdaMass, histAxes.axisPt});
127+
histos.add("reconstructedCandidates/hEtaDaughters", "hEtaDaughters", kTH1F, {histAxes.axisEta});
128+
histos.add("reconstructedCandidates/K0/hMass", "hMass", kTH3D, {histAxes.axisK0Mass, histAxes.axisPt, histAxes.axisEta});
129+
histos.add("reconstructedCandidates/Lambda/hMass", "hMass", kTH3D, {histAxes.axisLambdaMass, histAxes.axisPt, histAxes.axisEta});
129130
histos.add("reconstructedCandidates/hArmeterosBeforeAllSelections", "hArmeterosBeforeAllSelections", kTH2D, {{100, -1.0f, 1.0f}, {100, 0.0f, 0.25f}});
130131
histos.add("reconstructedCandidates/hArmeterosAfterAllSelections", "hArmeterosAfterAllSelections", kTH2D, {{100, -1.0f, 1.0f}, {100, 0.0f, 0.25f}});
131132
if (selectionFlags.doQAforSelectionVariables) {
@@ -217,28 +218,28 @@ struct alice3strangeness {
217218
bool isAntiLambda = (v0.mAntiLambda() - o2::constants::physics::MassLambda0) < 0.2;
218219

219220
histos.fill(HIST("reconstructedCandidates/hArmeterosBeforeAllSelections"), v0.alpha(), v0.qtarm());
220-
histos.fill(HIST("hV0CandidateCounter"), 1);
221+
histos.fill(HIST("hV0CandidateCounter"), 0.5);
221222
if (selectionFlags.applyRapiditySelection) {
222-
if (isK0 && std::abs(v0.yK0Short()) < selectionValues.yK0Selection)
223+
if (isK0 && std::abs(v0.yK0Short()) > selectionValues.yK0Selection)
223224
continue;
224225
if ((isLambda || isAntiLambda) && std::abs(v0.yLambda()) < selectionValues.yLambdaSelection)
225226
continue;
226227
}
227-
histos.fill(HIST("hV0CandidateCounter"), 2);
228+
histos.fill(HIST("hV0CandidateCounter"), 1.5);
228229
if (selectionFlags.applyDCAdaughterSelection) {
229230
if (std::abs(v0.dcaV0Daughters()) > selectionValues.dcaDaughterSelection)
230231
continue;
231232
} else {
232233
selectionCheck = v0.dcaV0Daughters();
233234
}
234-
histos.fill(HIST("hV0CandidateCounter"), 3);
235+
histos.fill(HIST("hV0CandidateCounter"), 2.5);
235236
if (selectionFlags.applyCosOfPAngleSelection) {
236237
if (v0.cosPA() < selectionValues.cosPAngleSelection)
237238
continue;
238239
} else {
239240
selectionCheck = std::acos(v0.cosPA());
240241
}
241-
histos.fill(HIST("hV0CandidateCounter"), 4);
242+
histos.fill(HIST("hV0CandidateCounter"), 3.5);
242243
if (selectionFlags.applyDCAdaughtersToPVSelection) {
243244
if ((std::abs(v0.dcaNegToPV()) < selectionValues.dcaDaughtersToPVSelection) ||
244245
(std::abs(v0.dcaPosToPV()) < selectionValues.dcaDaughtersToPVSelection))
@@ -247,21 +248,21 @@ struct alice3strangeness {
247248
selectionCheckPos = std::abs(v0.dcaPosToPV());
248249
selectionCheck = std::abs(v0.dcaNegToPV());
249250
}
250-
histos.fill(HIST("hV0CandidateCounter"), 5);
251+
histos.fill(HIST("hV0CandidateCounter"), 4.5);
251252
if (selectionFlags.applyV0RadiusSelection) {
252253
if (v0.v0radius() < selectionValues.v0RadiusSelection)
253254
continue;
254255
} else {
255256
selectionCheck = v0.v0radius();
256257
}
257-
histos.fill(HIST("hV0CandidateCounter"), 6);
258+
histos.fill(HIST("hV0CandidateCounter"), 5.5);
258259
if (isK0) {
259260
if (selectionFlags.applyArmenterosSelection) {
260261
if (v0.qtarm() < selectionValues.armenterosSelection * std::abs(v0.alpha()))
261262
continue;
262263
}
263264
}
264-
histos.fill(HIST("hV0CandidateCounter"), 7);
265+
histos.fill(HIST("hV0CandidateCounter"), 6.5);
265266
if (isK0 && selectionFlags.applyCompetingMassRejection) {
266267
if (std::abs(v0.mLambda() - o2::constants::physics::MassLambda0) < selectionValues.competingMassRejectionK0)
267268
continue;
@@ -272,7 +273,7 @@ struct alice3strangeness {
272273
if (std::abs(v0.mK0Short() - o2::constants::physics::MassK0Short) < selectionValues.competingMassRejectionLambda)
273274
continue;
274275
}
275-
histos.fill(HIST("hV0CandidateCounter"), 8);
276+
histos.fill(HIST("hV0CandidateCounter"), 7.5);
276277
if (selectionFlags.applyLifetimeSelection) {
277278
if (isK0 && v0.distovertotmom(collision.posX(), collision.posY(), collision.posZ()) * o2::constants::physics::MassK0Short < selectionValues.lifetimecutak0)
278279
continue;
@@ -284,14 +285,16 @@ struct alice3strangeness {
284285
else
285286
selectionCheck = v0.distovertotmom(collision.posX(), collision.posY(), collision.posZ()) * o2::constants::physics::MassLambda0;
286287
}
287-
histos.fill(HIST("hV0CandidateCounter"), 9);
288+
histos.fill(HIST("hV0CandidateCounter"), 8.5);
289+
auto posTrack = v0.template posTrack_as<alice3tracks>();
290+
auto negTrack = v0.template negTrack_as<alice3tracks>();
288291
if (selectionFlags.applyEtaDaughterSelection) {
289-
auto posTrack = v0.template posTrack_as<alice3tracks>();
290-
auto negTrack = v0.template negTrack_as<alice3tracks>();
291292
if (std::abs(posTrack.eta()) > selectionValues.etaDaughterSelection || std::abs(negTrack.eta()) > selectionValues.etaDaughterSelection)
292293
continue;
293294
}
294-
histos.fill(HIST("hV0CandidateCounter"), 10);
295+
histos.fill(HIST("reconstructedCandidates/hEtaDaughters"), posTrack.eta());
296+
histos.fill(HIST("reconstructedCandidates/hEtaDaughters"), negTrack.eta());
297+
histos.fill(HIST("hV0CandidateCounter"), 9.5);
295298

296299
histos.fill(HIST("reconstructedCandidates/hArmeterosAfterAllSelections"), v0.alpha(), v0.qtarm());
297300
if (selectionFlags.doQAforSelectionVariables) {
@@ -310,11 +313,11 @@ struct alice3strangeness {
310313
});
311314
}
312315
if (isK0)
313-
histos.fill(HIST("reconstructedCandidates/K0/hMass"), v0.mK0Short(), v0.pt());
316+
histos.fill(HIST("reconstructedCandidates/K0/hMass"), v0.mK0Short(), v0.pt(), v0.eta());
314317
if (isLambda)
315-
histos.fill(HIST("reconstructedCandidates/Lambda/hMass"), v0.mLambda(), v0.pt());
318+
histos.fill(HIST("reconstructedCandidates/Lambda/hMass"), v0.mLambda(), v0.pt(), v0.eta());
316319
if (isAntiLambda)
317-
histos.fill(HIST("reconstructedCandidates/AntiLambda/hMass"), v0.mAntiLambda(), v0.pt());
320+
histos.fill(HIST("reconstructedCandidates/AntiLambda/hMass"), v0.mAntiLambda(), v0.pt(), v0.eta());
318321
}
319322
}
320323
PROCESS_SWITCH(alice3strangeness, processAllFindableCandidates, "", false);

0 commit comments

Comments
 (0)