Skip to content

Commit dff9639

Browse files
ddobrigkddobrigk
committed
Minor adjustments
1 parent 3d665fe commit dff9639

File tree

1 file changed

+34
-2
lines changed

1 file changed

+34
-2
lines changed

‎PWGLF/Tasks/QC/strangepidqa.cxx‎

Lines changed: 34 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -220,6 +220,13 @@ struct strangepidqa {
220220
histos.add("h2dTimeVsLengthProtonProng", "h2dTimeVsLengthProtonProng", {HistType::kTH2F, {axisLength, axisTime}});
221221
histos.add("h2dTimeVsLengthPionProng", "h2dTimeVsLengthPionProng", {HistType::kTH2F, {axisLength, axisTime}});
222222

223+
histos.add("h1dMassK0Short", "h1dMassK0Short", {HistType::kTH1F, {axisK0ShortMass}});
224+
histos.add("h1dMassLambda", "h1dMassLambda", {HistType::kTH1F, {axisLambdaMass}});
225+
histos.add("h1dMassAntiLambda", "h1dMassAntiLambda", {HistType::kTH1F, {axisLambdaMass}});
226+
histos.add("h1dMassCompatibleK0Short", "h1dMassCompatibleK0Short", {HistType::kTH1F, {axisK0ShortMass}});
227+
histos.add("h1dMassCompatibleLambda", "h1dMassCompatibleLambda", {HistType::kTH1F, {axisLambdaMass}});
228+
histos.add("h1dMassCompatibleAntiLambda", "h1dMassCompatibleAntiLambda", {HistType::kTH1F, {axisLambdaMass}});
229+
223230
histos.add("h3dMassK0Short", "h3dMassK0Short", {HistType::kTH3F, {centAxis, axisPt, axisK0ShortMass}});
224231
histos.add("h3dMassLambda", "h3dMassLambda", {HistType::kTH3F, {centAxis, axisPt, axisLambdaMass}});
225232
histos.add("h3dMassAntiLambda", "h3dMassAntiLambda", {HistType::kTH3F, {centAxis, axisPt, axisLambdaMass}});
@@ -256,6 +263,15 @@ struct strangepidqa {
256263
}
257264

258265
if (doprocessCascades) {
266+
histos.add("h1dMassXiMinus", "h1dMassXiMinus", {HistType::kTH1F, {massAxisXi}});
267+
histos.add("h1dMassXiPlus", "h1dMassXiPlus", {HistType::kTH1F, {massAxisXi}});
268+
histos.add("h1dMassOmegaMinus", "h1dMassOmegaMinus", {HistType::kTH1F, {massAxisXi}});
269+
histos.add("h1dMassOmegaPlus", "h1dMassOmegaPlus", {HistType::kTH1F, {massAxisXi}});
270+
histos.add("h1dMassCompatibleXiMinus", "h1dMassCompatibleXiMinus", {HistType::kTH1F, {massAxisXi}});
271+
histos.add("h1dMassCompatibleXiPlus", "h1dMassCompatibleXiPlus", {HistType::kTH1F, {massAxisXi}});
272+
histos.add("h1dMassCompatibleOmegaMinus", "h1dMassCompatibleOmegaMinus", {HistType::kTH1F, {massAxisXi}});
273+
histos.add("h1dMassCompatibleOmegaPlus", "h1dMassCompatibleOmegaPlus", {HistType::kTH1F, {massAxisXi}});
274+
259275
histos.add("h3dMassXiMinus", "h3dMassXiMinus", {HistType::kTH3F, {centAxis, axisPt, massAxisXi}});
260276
histos.add("h3dMassXiPlus", "h3dMassXiPlus", {HistType::kTH3F, {centAxis, axisPt, massAxisXi}});
261277
histos.add("h3dMassOmegaMinus", "h3dMassOmegaMinus", {HistType::kTH3F, {centAxis, axisPt, massAxisOmega}});
@@ -284,24 +300,30 @@ struct strangepidqa {
284300
if (TMath::Abs(posExtra.tpcNSigmaPr()) < tpcNsigmaProton && TMath::Abs(negExtra.tpcNSigmaPi()) < tpcNsigmaPion) {
285301
// lambda case
286302
histos.fill(HIST("h3dMassLambda"), coll.centFT0C(), lambda.pt(), lambda.mLambda());
303+
histos.fill(HIST("h1dMassLambda"), lambda.mLambda());
287304
if(lambda.tofLambdaCompatibility(tofNsigmaCompatibility.value)){
288305
histos.fill(HIST("h3dMassCompatibleLambda"), coll.centFT0C(), lambda.pt(), lambda.mLambda());
306+
histos.fill(HIST("h1dMassCompatibleLambda"), lambda.mLambda());
289307
}
290308
}
291309

292310
if (TMath::Abs(posExtra.tpcNSigmaPi()) < tpcNsigmaProton && TMath::Abs(negExtra.tpcNSigmaPr()) < tpcNsigmaPion) {
293311
// lambda case
294-
histos.fill(HIST("h3dMassAntiLambda"), coll.centFT0C(), lambda.pt(), lambda.mLambda());
312+
histos.fill(HIST("h3dMassAntiLambda"), coll.centFT0C(), lambda.pt(), lambda.mAntiLambda());
313+
histos.fill(HIST("h1dMassAntiLambda"), lambda.mAntiLambda());
295314
if(lambda.tofAntiLambdaCompatibility(tofNsigmaCompatibility.value)){
296-
histos.fill(HIST("h3dMassCompatibleAntiLambda"), coll.centFT0C(), lambda.pt(), lambda.mLambda());
315+
histos.fill(HIST("h3dMassCompatibleAntiLambda"), coll.centFT0C(), lambda.pt(), lambda.mAntiLambda());
316+
histos.fill(HIST("h1dMassCompatibleAntiLambda"), lambda.mAntiLambda());
297317
}
298318
}
299319

300320
if (TMath::Abs(posExtra.tpcNSigmaPi()) < tpcNsigmaPion && TMath::Abs(negExtra.tpcNSigmaPr()) < tpcNsigmaPion) {
301321
// lambda case
302322
histos.fill(HIST("h3dMassK0Short"), coll.centFT0C(), lambda.pt(), lambda.mK0Short());
323+
histos.fill(HIST("h1dMassK0Short"), lambda.mK0Short());
303324
if(lambda.tofAntiLambdaCompatibility(tofNsigmaCompatibility.value)){
304325
histos.fill(HIST("h3dMassCompatibleK0Short"), coll.centFT0C(), lambda.pt(), lambda.mK0Short());
326+
histos.fill(HIST("h1dMassCompatibleK0Short"), lambda.mK0Short());
305327
}
306328
}
307329

@@ -421,6 +443,8 @@ struct strangepidqa {
421443
void processCascades(soa::Join<aod::StraCollisions, aod::StraCents>::iterator const& col, soa::Filtered<soa::Join<aod::CascCores, aod::CascCollRefs, aod::CascExtras, aod::CascTOFNSigmas>> const& Cascades, soa::Join<aod::DauTrackExtras, aod::DauTrackTPCPIDs> const&)
422444
{
423445
histos.fill(HIST("hEventCentrality"), col.centFT0C());
446+
if (coll.centFT0C() > maxCentrality || coll.centFT0C() < minCentrality)
447+
return;
424448

425449
for (auto& casc : Cascades) {
426450
// major selections here
@@ -441,28 +465,36 @@ struct strangepidqa {
441465
if (casc.sign() < 0) {
442466
if (TMath::Abs(posExtra.tpcNSigmaPr()) < tpcNsigmaProton && TMath::Abs(negExtra.tpcNSigmaPi()) < tpcNsigmaPion && TMath::Abs(bachExtra.tpcNSigmaPi()) < tpcNsigmaBachelor) {
443467
histos.fill(HIST("h3dMassXiMinus"), col.centFT0C(), casc.pt(), casc.mXi());
468+
histos.fill(HIST("h1dMassXiMinus"), casc.mXi());
444469
if(casc.tofXiCompatibility(tofNsigmaCompatibilityCascades.value)){
445470
histos.fill(HIST("h3dMassCompatibleXiMinus"), col.centFT0C(), casc.pt(), casc.mXi());
471+
histos.fill(HIST("h1dMassCompatibleXiMinus"), casc.mXi());
446472
}
447473
}
448474
if (TMath::Abs(posExtra.tpcNSigmaPr()) < tpcNsigmaProton && TMath::Abs(negExtra.tpcNSigmaPi()) < tpcNsigmaPion && TMath::Abs(bachExtra.tpcNSigmaKa()) < tpcNsigmaBachelor) {
449475
histos.fill(HIST("h3dMassOmegaMinus"), col.centFT0C(), casc.pt(), casc.mOmega());
476+
histos.fill(HIST("h1dMassOmegaMinus"), casc.mOmega());
450477
if(casc.tofOmCompatibility(tofNsigmaCompatibilityCascades.value)){
451478
histos.fill(HIST("h3dMassCompatibleOmegaMinus"), col.centFT0C(), casc.pt(), casc.mOmega());
479+
histos.fill(HIST("h1dMassCompatibleOmegaMinus"), casc.mOmega());
452480
}
453481
}
454482
} else {
455483
if (TMath::Abs(posExtra.tpcNSigmaPi()) < tpcNsigmaPion && TMath::Abs(negExtra.tpcNSigmaPr()) < tpcNsigmaProton && TMath::Abs(bachExtra.tpcNSigmaPi()) < tpcNsigmaBachelor) {
456484
histos.fill(HIST("h3dMassXiPlus"), col.centFT0C(), casc.pt(), casc.mXi());
485+
histos.fill(HIST("h1dMassXiPlus"), casc.mXi());
457486
if(casc.tofXiCompatibility(tofNsigmaCompatibilityCascades.value)){
458487
histos.fill(HIST("h3dMassCompatibleXiPlus"), col.centFT0C(), casc.pt(), casc.mXi());
488+
histos.fill(HIST("h1dCompatibleMassXiPlus"), casc.mXi());
459489
}
460490
}
461491

462492
if (TMath::Abs(posExtra.tpcNSigmaPi()) < tpcNsigmaPion && TMath::Abs(negExtra.tpcNSigmaPr()) < tpcNsigmaProton && TMath::Abs(bachExtra.tpcNSigmaKa()) < tpcNsigmaBachelor) {
463493
histos.fill(HIST("h3dMassOmegaPlus"), col.centFT0C(), casc.pt(), casc.mOmega());
494+
histos.fill(HIST("h1dMassOmegaPlus"), casc.mOmega());
464495
if(casc.tofOmCompatibility(tofNsigmaCompatibilityCascades.value)){
465496
histos.fill(HIST("h3dMassCompatibleOmegaPlus"), col.centFT0C(), casc.pt(), casc.mOmega());
497+
histos.fill(HIST("h1dMassCompatibleOmegaPlus"), casc.mOmega());
466498
}
467499
}
468500
}

0 commit comments

Comments
 (0)