Skip to content

Commit 686c3db

Browse files
omvazqueomvazque
authored
[PWGLF] Loosen the gamma conversions criteria + dEdx dist. (#12830)
1 parent cd5a1c4 commit 686c3db

File tree

1 file changed

+95
-81
lines changed

1 file changed

+95
-81
lines changed

‎PWGLF/Tasks/Nuspex/piKpRAA.cxx‎

Lines changed: 95 additions & 81 deletions
Original file line numberDiff line numberDiff line change
@@ -77,6 +77,8 @@ std::array<std::shared_ptr<TH3>, kNEtaHists> dEdxPiV0{};
7777
std::array<std::shared_ptr<TH3>, kNEtaHists> dEdxPrV0{};
7878
std::array<std::shared_ptr<TH3>, kNEtaHists> dEdxElV0{};
7979
std::array<std::shared_ptr<TH3>, kNEtaHists> dEdxPiTOF{};
80+
std::array<std::shared_ptr<TH3>, kNEtaHists> dEdx{};
81+
std::array<std::shared_ptr<TH2>, kNEtaHists> pTVsP{};
8082

8183
struct PiKpRAA {
8284

@@ -126,12 +128,13 @@ struct PiKpRAA {
126128
// original equation: lArmPt*5>TMath::Abs(lArmAlpha)
127129
Configurable<float> armPodCut{"armPodCut", 5.0f, "pT * (cut) > |alpha|, AP cut. Negative: no cut"};
128130
Configurable<float> armAlphaSel{"armAlphaSel", 0.45f, "Armenteros alpha selection (Gammas)"};
129-
Configurable<float> qTSel{"qTSel", 0.1f, "Armenteros qT select (Gammas)"};
131+
Configurable<float> qTSel{"qTSel", 0.01f, "Armenteros qT select (Gammas)"};
130132

131133
// Selection
132134
Configurable<bool> applyInvMassSel{"applyInvMassSel", false, "Select V0s close to the Inv. mass value"};
133135
Configurable<float> dMassSel{"dMassSel", 0.01f, "Invariant mass selection"};
134136
Configurable<float> dMassSelG{"dMassSelG", 0.1f, "Inv mass selection gammas"};
137+
Configurable<float> dMassGcut{"dMassGcut", 0.01f, "Inv mass selection gammas"};
135138

136139
// PID (TPC/TOF)
137140
Configurable<float> tpcPidNsigmaCut{"tpcPidNsigmaCut", 5, "tpcPidNsigmaCut"};
@@ -173,7 +176,7 @@ struct PiKpRAA {
173176
ConfigurableAxis binsPt{"binsPt", {VARIABLE_WIDTH, 0.0, 0.1, 0.12}, "pT binning"};
174177
ConfigurableAxis binsCent{"binsCent", {VARIABLE_WIDTH, 0., 5., 10., 20., 30., 40., 50., 60., 70., 80., 90., 100.}, "T0C binning"};
175178
ConfigurableAxis axisArmAlpha{"axisArmAlpha", {200, -1.0, 1.0}, "Armenteros alpha"};
176-
ConfigurableAxis axisArmqT{"axisArmqT", {200, 0.0f, 0.3f}, "Armenteros qT"};
179+
ConfigurableAxis axisArmqT{"axisArmqT", {600, 0.0f, 0.3f}, "Armenteros qT"};
177180
ConfigurableAxis axisK0Mass{"axisK0Mass", {200, 0.4f, 0.6f}, "Mass K0Short"};
178181
ConfigurableAxis axisLambdaMass{"axisLambdaMass", {200, 1.101f, 1.131f}, "Mass Lambda"};
179182
ConfigurableAxis axisGammaMass{"axisGammaMass", {200, 0.0f, 0.5f}, "Mass Gamma"};
@@ -289,6 +292,7 @@ struct PiKpRAA {
289292
registry.add("ExcludedEvtVsNch", ";Nch;Entries;", kTH1F, {{nBinsNch, minNch, maxNch}});
290293
registry.add("ExcludedEvtVsNPV", ";NPV;Entries;", kTH1F, {{nBinsNPV, minNpv, maxNpv}});
291294

295+
registry.add("dcaVsPt", "Primary particles;#it{p}_{T} (GeV/#it{c});DCA_{xy} (cm);", kTH2F, {{{axisPt}, {40, -1.0, 1.0}}});
292296
registry.add("dcaDauVsPt", ";V0 #it{p}_{T} (GeV/#it{c});DCA_{xy} (cm) daughters;", kTH2F, {{{axisPtV0s}, {200, -10., 10.}}});
293297
registry.add("nSigPiFromK0s", ";#it{n#sigma};;", kTH2F, {axisPtV0s, axisNsigmaTPC});
294298
registry.add("nSigPiFromL", ";#it{n#sigma};;", kTH2F, {axisPtV0s, axisNsigmaTPC});
@@ -314,7 +318,8 @@ struct PiKpRAA {
314318
registry.add("NclFoundp", ";;Found #LTNcl#GT TPC", kTProfile, {axisXNcl});
315319
registry.add("NclFoundVsPhipBeforeCut", Form("Found #LTNcl#GT TPC;%s (GeV/#it{c});#varphi", titlePorPt.data()), kTProfile2D, {{{axisXPhiCut}, {350, 0.0, 0.35}}});
316320
registry.add("NclFoundVsPhipAfterCut", Form("Found #LTNcl#GT TPC;%s (GeV/#it{c});#varphi", titlePorPt.data()), kTProfile2D, {{{axisXPhiCut}, {350, 0.0, 0.35}}});
317-
registry.add("NclVsEta", ";#eta; Ncl TPC", kTH2F, {{{axisEta}, {161, -0.5, 160.5}}});
321+
registry.add("NclVsEta", ";#eta;Found Ncl TPC", kTH2F, {{{axisEta}, {161, -0.5, 160.5}}});
322+
registry.add("NclVsEtap", ";#eta;Found #LTNcl#GT TPC", kTProfile, {axisEta});
318323

319324
registry.add("NclVsEtaPiMIP", "MIP #pi^{+} + #pi^{-} (0.4 < #it{p} < 0.6 GeV/#it{c}, 40 < dE/dx < 60);#eta; Found Ncl TPC", kTH2F, {{{axisEta}, {161, -0.5, 160.5}}});
320325
registry.add("NclVsEtaPiMIPp", "MIP #pi^{+} + #pi^{-} (0.4 < #it{p} < 0.6 GeV/#it{c}, 40 < dE/dx < 60);#eta; Found #LTNcl#GT TPC", kTProfile, {axisEta});
@@ -335,7 +340,6 @@ struct PiKpRAA {
335340
registry.add("TOFExpEl2TOF", ";Momentum (GeV/#it{c});t^{#pi}_{Exp}/t_{TOF}", kTH2F, {{{axisPtV0s}, {100, 0.2, 1.2}}});
336341

337342
registry.add("betaVsMomentum", ";Momentum (GeV/#it{c}); #beta", kTH2F, {{{axisPtV0s}, {500, 0, 1.2}}});
338-
registry.add("dEdxVsMomentum", ";Momentum (GeV/#it{c}); dE/dx", kTH2F, {axisPtV0s, axisdEdx});
339343
registry.add("dEdxVsEtaPiMIP", "MIP #pi^{+} + #pi^{-} (0.4 < #it{p} < 0.6 GeV/#it{c});#eta; dE/dx;", kTH2F, {{{axisEta}, {100, 0, 100}}});
340344
registry.add("dEdxVsEtaPiMIPp", "MIP #pi^{+} + #pi^{-} (0.4 < #it{p} < 0.6 GeV/#it{c});#eta; #LTdE/dx#GT", kTProfile, {axisEta});
341345
registry.add("dEdxVsEtaElMIP", "MIP e^{+} + e^{-} (0.3 < #it{p} < 0.4 GeV/#it{c});#eta; dE/dx;", kTH2F, {{{axisEta}, {100, 0, 100}}});
@@ -346,6 +350,8 @@ struct PiKpRAA {
346350
registry.add("dEdxVsEtaElMIPV0p", "e^{+} + e^{-} (0.4 < #it{p} < 0.6 GeV/#it{c});#eta; #LTdE/dx#GT", kTProfile, {axisEta});
347351

348352
for (int i = 0; i < kNEtaHists; ++i) {
353+
dEdx[i] = registry.add<TH3>(Form("dEdx_%s", endingEta[i]), Form("%s;Momentum (GeV/#it{c});dE/dx;", latexEta[i]), kTH3F, {axisPt, axisdEdx, axisCent});
354+
pTVsP[i] = registry.add<TH2>(Form("pTVsP_%s", endingEta[i]), Form("%s;Momentum (GeV/#it{c});#it{p}_{T} (GeV/#it{c});", latexEta[i]), kTH2F, {axisPt, axisPt});
349355
dEdxPiV0[i] = registry.add<TH3>(Form("dEdxPiV0_%s", endingEta[i]), Form("#pi^{+} + #pi^{-}, %s;Momentum (GeV/#it{c});dE/dx;", latexEta[i]), kTH3F, {axisPtV0s, axisdEdx, axisCent});
350356
dEdxPrV0[i] = registry.add<TH3>(Form("dEdxPrV0_%s", endingEta[i]), Form("p + #bar{p}, %s;Momentum (GeV/#it{c});dE/dx;", latexEta[i]), kTH3F, {axisPtV0s, axisdEdx, axisCent});
351357
dEdxElV0[i] = registry.add<TH3>(Form("dEdxElV0_%s", endingEta[i]), Form("e^{+} + e^{-}, %s;Momentum (GeV/#it{c});dE/dx;", latexEta[i]), kTH3F, {axisPtV0s, axisdEdx, axisCent});
@@ -358,6 +364,8 @@ struct PiKpRAA {
358364
LOG(info) << "\tdetector4Calibration=" << detector4Calibration.value;
359365
LOG(info) << "\tminPt=" << v0Selections.minPt;
360366
LOG(info) << "\tmaxPt=" << v0Selections.maxPt;
367+
LOG(info) << "\tqTSel=" << v0Selections.qTSel;
368+
LOG(info) << "\tarmAlphaSel=" << v0Selections.armAlphaSel;
361369
LOG(info) << "\tapplyTPCTOFCombinedCut=" << v0Selections.applyTPCTOFCombinedCut;
362370
LOG(info) << "\tapplyInvMassSel=" << v0Selections.applyInvMassSel;
363371
LOG(info) << "\tapplyNclSel=" << v0Selections.applyNclSel;
@@ -515,6 +523,14 @@ struct PiKpRAA {
515523
continue;
516524
}
517525

526+
int indexEta{0};
527+
for (int i = 1; i < kNEtaHists; ++i) {
528+
if (std::abs(eta) >= kLowEta[i] && std::abs(eta) < kHighEta[i]) {
529+
indexEta = i;
530+
break;
531+
}
532+
}
533+
518534
if (momentum > kMinPMIP && momentum < kMaxPMIP && dedx > kMindEdxMIP && dedx < kMaxdEdxMIP) {
519535
registry.fill(HIST("dEdxVsEtaPiMIP"), eta, dedx);
520536
registry.fill(HIST("dEdxVsEtaPiMIPp"), eta, dedx);
@@ -527,22 +543,17 @@ struct PiKpRAA {
527543
registry.fill(HIST("dEdxVsEtaElMIPp"), eta, dedx);
528544
}
529545

530-
registry.fill(HIST("dEdxVsMomentum"), momentum, dedx);
546+
dEdx[indexEta]->Fill(momentum, dedx, centrality);
547+
pTVsP[indexEta]->Fill(momentum, pt);
548+
registry.fill(HIST("dcaVsPt"), pt, track.dcaXY());
531549
registry.fill(HIST("NclVsEta"), eta, track.tpcNClsFound());
550+
registry.fill(HIST("NclVsEtap"), eta, track.tpcNClsFound());
532551
registry.fill(HIST("NclFound"), pOrPt, track.tpcNClsFound());
533552
registry.fill(HIST("NclFoundp"), pOrPt, track.tpcNClsFound());
534553
registry.fill(HIST("NclFindable"), pOrPt, track.tpcNClsFindable());
535554
registry.fill(HIST("NclFindablep"), pOrPt, track.tpcNClsFindable());
536555
registry.fill(HIST("NclFoundVsPhipAfterCut"), pOrPt, phiPrime, track.tpcNClsFound());
537556

538-
int indexEta{0};
539-
for (int i = 1; i < kNEtaHists; ++i) {
540-
if (std::abs(eta) >= kLowEta[i] && std::abs(eta) < kHighEta[i]) {
541-
indexEta = i;
542-
break;
543-
}
544-
}
545-
546557
if (track.hasTOF() && track.goodTOFMatch()) {
547558
const float tTOF{track.tofSignal()};
548559
const float trkLength{track.length()};
@@ -631,19 +642,16 @@ struct PiKpRAA {
631642
getArmeterosVariables(ppos, pneg, alpha, qT);
632643
registry.fill(HIST("ArmAll"), alpha, qT);
633644

645+
bool passesTopoSel{false};
634646
// Passes V0 topological cuts?
635-
if (!passesV0TopologicalSelection(v0))
636-
continue;
647+
if (passesV0TopologicalSelection(v0))
648+
passesTopoSel = true;
637649

638650
const double dMassK0s{std::abs(v0.mK0Short() - o2::constants::physics::MassK0Short)};
639651
const double dMassL{std::abs(v0.mLambda() - o2::constants::physics::MassLambda0)};
640652
const double dMassAL{std::abs(v0.mAntiLambda() - o2::constants::physics::MassLambda0)};
641653
const double dMassG{std::abs(v0.mGamma() - o2::constants::physics::MassGamma)};
642654

643-
registry.fill(HIST("ArmAfterTopoSel"), alpha, qT);
644-
registry.fill(HIST("dcaDauVsPt"), v0.pt(), v0.dcapostopv());
645-
registry.fill(HIST("dcaDauVsPt"), v0.pt(), v0.dcanegtopv());
646-
647655
int posIndexEta{0};
648656
int negIndexEta{0};
649657
for (int i = 1; i < kNEtaHists; ++i) {
@@ -660,81 +668,87 @@ struct PiKpRAA {
660668
}
661669
}
662670

663-
if (v0Selections.applyInvMassSel) { // apply Inv. Mass selection?
664-
if (dMassK0s < v0Selections.dMassSel && dMassL > v0Selections.dMassSel && dMassAL > v0Selections.dMassSel && dMassG > v0Selections.dMassSelG) { // Mass cut
665-
if (passesK0Selection(collision, v0)) { // nSigma TPC and y cuts
666-
registry.fill(HIST("ArmK0NOSel"), alpha, qT);
667-
if (v0Selections.armPodCut * qT > std::abs(alpha)) { // Armenters selection
668-
registry.fill(HIST("ArmK0"), alpha, qT);
669-
registry.fill(HIST("MassK0sVsPt"), v0.pt(), v0.mK0Short());
670-
registry.fill(HIST("nSigPiFromK0s"), posTrkPt, posTrack.tpcNSigmaPi());
671-
registry.fill(HIST("nSigPiFromK0s"), negTrkPt, negTrack.tpcNSigmaPi());
672-
registry.fill(HIST("NclVsEtaPiV0"), posTrkEta, posTrack.tpcNClsFound());
673-
registry.fill(HIST("NclVsEtaPiV0p"), posTrkEta, posTrack.tpcNClsFound());
671+
if (passesTopoSel) {
672+
registry.fill(HIST("ArmAfterTopoSel"), alpha, qT);
673+
registry.fill(HIST("dcaDauVsPt"), v0.pt(), v0.dcapostopv());
674+
registry.fill(HIST("dcaDauVsPt"), v0.pt(), v0.dcanegtopv());
675+
676+
if (v0Selections.applyInvMassSel) { // apply Inv. Mass selection?
677+
if (dMassK0s < v0Selections.dMassSel && dMassL > v0Selections.dMassSel && dMassAL > v0Selections.dMassSel && dMassG > v0Selections.dMassSelG) { // Mass cut
678+
if (passesK0Selection(collision, v0)) { // nSigma TPC and y cuts
679+
registry.fill(HIST("ArmK0NOSel"), alpha, qT);
680+
if (v0Selections.armPodCut * qT > std::abs(alpha)) { // Armenters selection
681+
registry.fill(HIST("ArmK0"), alpha, qT);
682+
registry.fill(HIST("MassK0sVsPt"), v0.pt(), v0.mK0Short());
683+
registry.fill(HIST("nSigPiFromK0s"), posTrkPt, posTrack.tpcNSigmaPi());
684+
registry.fill(HIST("nSigPiFromK0s"), negTrkPt, negTrack.tpcNSigmaPi());
685+
registry.fill(HIST("NclVsEtaPiV0"), posTrkEta, posTrack.tpcNClsFound());
686+
registry.fill(HIST("NclVsEtaPiV0p"), posTrkEta, posTrack.tpcNClsFound());
687+
registry.fill(HIST("NclVsEtaPiV0"), negTrkEta, negTrack.tpcNClsFound());
688+
registry.fill(HIST("NclVsEtaPiV0p"), negTrkEta, negTrack.tpcNClsFound());
689+
registry.fill(HIST("NclPiV0"), posPorPt, posTrack.tpcNClsFound());
690+
registry.fill(HIST("NclPiV0p"), posPorPt, posTrack.tpcNClsFound());
691+
registry.fill(HIST("NclPiV0"), negPorPt, negTrack.tpcNClsFound());
692+
registry.fill(HIST("NclPiV0p"), negPorPt, negTrack.tpcNClsFound());
693+
694+
dEdxPiV0[posIndexEta]->Fill(posTrkP, posTrkdEdx, centrality);
695+
dEdxPiV0[negIndexEta]->Fill(negTrkP, negTrkdEdx, centrality);
696+
697+
if (posTrkP > kMinPMIP && posTrkP < kMaxPMIP && posTrkdEdx > kMindEdxMIP && posTrkdEdx < kMaxdEdxMIP) {
698+
registry.fill(HIST("dEdxVsEtaPiMIPV0"), posTrkEta, posTrkdEdx);
699+
registry.fill(HIST("dEdxVsEtaPiMIPV0p"), posTrkEta, posTrkdEdx);
700+
}
701+
if (negTrkP > kMinPMIP && negTrkP < kMaxPMIP && negTrkdEdx > kMindEdxMIP && negTrkdEdx < kMaxdEdxMIP) {
702+
registry.fill(HIST("dEdxVsEtaPiMIPV0"), negTrkEta, negTrkdEdx);
703+
registry.fill(HIST("dEdxVsEtaPiMIPV0p"), negTrkEta, negTrkdEdx);
704+
}
705+
}
706+
}
707+
}
708+
}
709+
710+
if (v0Selections.applyInvMassSel) {
711+
if (dMassL < v0Selections.dMassSel && dMassK0s > v0Selections.dMassSel && dMassG > v0Selections.dMassSelG) {
712+
if (passesLambdaSelection(collision, v0)) {
713+
registry.fill(HIST("ArmL"), alpha, qT);
714+
registry.fill(HIST("MassLVsPt"), v0.pt(), v0.mLambda());
715+
registry.fill(HIST("nSigPrFromL"), posTrkPt, posTrack.tpcNSigmaPr());
716+
registry.fill(HIST("nSigPiFromL"), negTrkPt, negTrack.tpcNSigmaPi());
717+
registry.fill(HIST("NclVsEtaPrV0"), posTrkEta, posTrack.tpcNClsFound());
718+
registry.fill(HIST("NclVsEtaPrV0p"), posTrkEta, posTrack.tpcNClsFound());
674719
registry.fill(HIST("NclVsEtaPiV0"), negTrkEta, negTrack.tpcNClsFound());
675720
registry.fill(HIST("NclVsEtaPiV0p"), negTrkEta, negTrack.tpcNClsFound());
676-
registry.fill(HIST("NclPiV0"), posPorPt, posTrack.tpcNClsFound());
677-
registry.fill(HIST("NclPiV0p"), posPorPt, posTrack.tpcNClsFound());
721+
registry.fill(HIST("NclPrV0"), posPorPt, posTrack.tpcNClsFound());
722+
registry.fill(HIST("NclPrV0p"), posPorPt, posTrack.tpcNClsFound());
678723
registry.fill(HIST("NclPiV0"), negPorPt, negTrack.tpcNClsFound());
679724
registry.fill(HIST("NclPiV0p"), negPorPt, negTrack.tpcNClsFound());
680-
681-
dEdxPiV0[posIndexEta]->Fill(posTrkP, posTrkdEdx, centrality);
725+
dEdxPrV0[posIndexEta]->Fill(posTrkP, posTrkdEdx, centrality);
682726
dEdxPiV0[negIndexEta]->Fill(negTrkP, negTrkdEdx, centrality);
683-
684-
if (posTrkP > kMinPMIP && posTrkP < kMaxPMIP && posTrkdEdx > kMindEdxMIP && posTrkdEdx < kMaxdEdxMIP) {
685-
registry.fill(HIST("dEdxVsEtaPiMIPV0"), posTrkEta, posTrkdEdx);
686-
registry.fill(HIST("dEdxVsEtaPiMIPV0p"), posTrkEta, posTrkdEdx);
687-
}
688-
if (negTrkP > kMinPMIP && negTrkP < kMaxPMIP && negTrkdEdx > kMindEdxMIP && negTrkdEdx < kMaxdEdxMIP) {
689-
registry.fill(HIST("dEdxVsEtaPiMIPV0"), negTrkEta, negTrkdEdx);
690-
registry.fill(HIST("dEdxVsEtaPiMIPV0p"), negTrkEta, negTrkdEdx);
691-
}
692727
}
693728
}
694729
}
695-
}
696730

697-
if (v0Selections.applyInvMassSel) {
698-
if (dMassL < v0Selections.dMassSel && dMassK0s > v0Selections.dMassSel && dMassG > v0Selections.dMassSelG) {
699-
if (passesLambdaSelection(collision, v0)) {
700-
registry.fill(HIST("ArmL"), alpha, qT);
701-
registry.fill(HIST("MassLVsPt"), v0.pt(), v0.mLambda());
702-
registry.fill(HIST("nSigPrFromL"), posTrkPt, posTrack.tpcNSigmaPr());
703-
registry.fill(HIST("nSigPiFromL"), negTrkPt, negTrack.tpcNSigmaPi());
704-
registry.fill(HIST("NclVsEtaPrV0"), posTrkEta, posTrack.tpcNClsFound());
705-
registry.fill(HIST("NclVsEtaPrV0p"), posTrkEta, posTrack.tpcNClsFound());
706-
registry.fill(HIST("NclVsEtaPiV0"), negTrkEta, negTrack.tpcNClsFound());
707-
registry.fill(HIST("NclVsEtaPiV0p"), negTrkEta, negTrack.tpcNClsFound());
708-
registry.fill(HIST("NclPrV0"), posPorPt, posTrack.tpcNClsFound());
709-
registry.fill(HIST("NclPrV0p"), posPorPt, posTrack.tpcNClsFound());
710-
registry.fill(HIST("NclPiV0"), negPorPt, negTrack.tpcNClsFound());
711-
registry.fill(HIST("NclPiV0p"), negPorPt, negTrack.tpcNClsFound());
712-
dEdxPrV0[posIndexEta]->Fill(posTrkP, posTrkdEdx, centrality);
713-
dEdxPiV0[negIndexEta]->Fill(negTrkP, negTrkdEdx, centrality);
731+
if (v0Selections.applyInvMassSel && dMassAL < v0Selections.dMassSel && dMassK0s > v0Selections.dMassSel && dMassG > v0Selections.dMassSelG) {
732+
if (passesAntiLambdaSelection(collision, v0)) {
733+
registry.fill(HIST("ArmAL"), alpha, qT);
734+
registry.fill(HIST("MassALVsPt"), v0.pt(), v0.mAntiLambda());
735+
registry.fill(HIST("nSigPrFromAL"), negTrkPt, negTrack.tpcNSigmaPr());
736+
registry.fill(HIST("nSigPiFromAL"), posTrkPt, posTrack.tpcNSigmaPi());
737+
registry.fill(HIST("NclVsEtaPiV0"), posTrkEta, posTrack.tpcNClsFound());
738+
registry.fill(HIST("NclVsEtaPiV0p"), posTrkEta, posTrack.tpcNClsFound());
739+
registry.fill(HIST("NclVsEtaPrV0"), negTrkEta, negTrack.tpcNClsFound());
740+
registry.fill(HIST("NclVsEtaPrV0p"), negTrkEta, negTrack.tpcNClsFound());
741+
registry.fill(HIST("NclPiV0"), posPorPt, posTrack.tpcNClsFound());
742+
registry.fill(HIST("NclPiV0p"), posPorPt, posTrack.tpcNClsFound());
743+
registry.fill(HIST("NclPrV0"), negPorPt, negTrack.tpcNClsFound());
744+
registry.fill(HIST("NclPrV0p"), negPorPt, negTrack.tpcNClsFound());
745+
dEdxPrV0[negIndexEta]->Fill(negTrkP, negTrkdEdx, centrality);
746+
dEdxPiV0[posIndexEta]->Fill(posTrkP, posTrkdEdx, centrality);
714747
}
715748
}
716749
}
717750

718-
if (v0Selections.applyInvMassSel && dMassAL < v0Selections.dMassSel && dMassK0s > v0Selections.dMassSel && dMassG > v0Selections.dMassSelG) {
719-
if (passesAntiLambdaSelection(collision, v0)) {
720-
registry.fill(HIST("ArmAL"), alpha, qT);
721-
registry.fill(HIST("MassALVsPt"), v0.pt(), v0.mAntiLambda());
722-
registry.fill(HIST("nSigPrFromAL"), negTrkPt, negTrack.tpcNSigmaPr());
723-
registry.fill(HIST("nSigPiFromAL"), posTrkPt, posTrack.tpcNSigmaPi());
724-
registry.fill(HIST("NclVsEtaPiV0"), posTrkEta, posTrack.tpcNClsFound());
725-
registry.fill(HIST("NclVsEtaPiV0p"), posTrkEta, posTrack.tpcNClsFound());
726-
registry.fill(HIST("NclVsEtaPrV0"), negTrkEta, negTrack.tpcNClsFound());
727-
registry.fill(HIST("NclVsEtaPrV0p"), negTrkEta, negTrack.tpcNClsFound());
728-
registry.fill(HIST("NclPiV0"), posPorPt, posTrack.tpcNClsFound());
729-
registry.fill(HIST("NclPiV0p"), posPorPt, posTrack.tpcNClsFound());
730-
registry.fill(HIST("NclPrV0"), negPorPt, negTrack.tpcNClsFound());
731-
registry.fill(HIST("NclPrV0p"), negPorPt, negTrack.tpcNClsFound());
732-
dEdxPrV0[negIndexEta]->Fill(negTrkP, negTrkdEdx, centrality);
733-
dEdxPiV0[posIndexEta]->Fill(posTrkP, posTrkdEdx, centrality);
734-
}
735-
}
736-
737-
if (v0Selections.applyInvMassSel && dMassK0s > v0Selections.dMassSel && dMassL > v0Selections.dMassSel && dMassAL > v0Selections.dMassSel && dMassG < v0Selections.dMassSel) {
751+
if (v0Selections.applyInvMassSel && dMassK0s > v0Selections.dMassSel && dMassL > v0Selections.dMassSel && dMassAL > v0Selections.dMassSel && dMassG < v0Selections.dMassGcut) {
738752
if (passesGammaSelection(collision, v0)) {
739753
if (std::abs(alpha) < v0Selections.armAlphaSel && qT < v0Selections.qTSel) {
740754
registry.fill(HIST("ArmG"), alpha, qT);

0 commit comments

Comments
 (0)