[PWGHF] Added trigger for Sc-proton correlation #14786

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

Open

singhra1994 wants to merge 10 commits into AliceO2Group:master from singhra1994:master

+234 −53

Collaborator

singhra1994 commented Feb 3, 2026

I have modified the block trigger SigmaCPPK in such a way that both SigmaCPPK and SigmaCP can be filled in the same block

singhra1994 added 3 commits

February 3, 2026 17:00


          added SigmaC-proton trigger

21fcb3a

Added parameters for Sigma_C trigger selection and updated event tagging logic for SigmaC candidates.


          added SigmaC-proton trigger

4ac7d66

Added support for SigmaC particle selection in the HfFilterHelper class.


          Add HfSigmaCP column to filterTables.h

3b5f08b

singhra1994 requested review from alibuild, ariedel-cern, fgrosa, lietava, mpuccio and zhangbiao-phy as code owners

February 3, 2026 16:19

github-actions bot added pwghf trigger labels


          Please consider the following formatting changes

21cc8df

alibuild mentioned this pull request

[Trigger,PWGHF] Please consider the following formatting changes to #14786 singhra1994/O2Physics#43

Merged


          Merge pull request #43 from alibuild/alibot-cleanup-14786

568051e

[Trigger,PWGHF] Please consider the following formatting changes to AliceO2Group#14786

fgrosa requested changes

View reviewed changes

Collaborator

fgrosa left a comment

Hi @singhra1994 I still didn't look into the helper functions, but I already spotted several issues in the logic of the filter itself. Please implement the suggestions, and then I will take a look at the helper functions as well.

EventFiltering/filterTables.h Outdated

               DECLARE_SOA_COLUMN(HfPrCharm2P, hasHfPrCharm2P, bool);                           //! Charm baryon to 2-prong + bachelors
               DECLARE_SOA_COLUMN(HfSigmaCPPK, hasHfSigmaCPPK, bool);                           //! SigmaC(2455)++K- and SigmaC(2520)++K- + c.c.
               DECLARE_SOA_COLUMN(HfSigmaC0K0, hasHfSigmaC0K0, bool);                           //! SigmaC(2455)0KS0 and SigmaC(2520)0KS0
+              DECLARE_SOA_COLUMN(HfSigmaCP, hasHfSigmaCP, bool);                               //! SigmaC(2455)P

Collaborator

fgrosa Feb 3, 2026

Please use Pr and not P to avoid confusion since in the other SigmaC triggers P is used to indicate the charge state

Suggested change

      
            DECLARE_SOA_COLUMN(HfSigmaCP, hasHfSigmaCP, bool);                               //! SigmaC(2455)P
          
            DECLARE_SOA_COLUMN(HfSigmaCPr, hasHfSigmaCPr, bool);                               //! SigmaC(2455)-Proton pairs

EventFiltering/filterTables.h Outdated

                                 filtering::HfCharmBarToXiBach,
                                 filtering::HfSigmaCPPK,
                                 filtering::HfSigmaC0K0,
+                                filtering::HfSigmaCP,

Collaborator

fgrosa Feb 3, 2026

Suggested change

      
                              filtering::HfSigmaCP,
          
                              filtering::HfSigmaCPr,

EventFiltering/PWGHF/HFFilter.cxx Outdated

+                                auto tagBDT = helper.isBDTSelected(scores[2], thresholdBDTScores[5]);
+                                if (helper.selectionSigmaCForScPCorr(pVecFirst, pVecThird, pVecSecond, pVecSoftPi, ptSigmaC, is3Prong[2], hMassVsPtC[kNCharmParticles + 23], activateQA, configSigmaC.minMassSigmaCCorr, configSigmaC.maxMassSigmaCCorr, configSigmaC.minPtSigmaC, configSigmaC.maxPtSigmaC) && TESTBIT(tagBDT, RecoDecay::OriginType::Prompt)) {
+                                  keepEvent[kSigmaCP] = true;

Collaborator

fgrosa Feb 3, 2026

Suggested change

      
                                keepEvent[kSigmaCP] = true;
          
                                keepEvent[kSigmaCPr] = true;

EventFiltering/PWGHF/HFFilter.cxx Outdated

                     }
-                    tags(keepEvent[kHighPt2P], keepEvent[kHighPt3P], keepEvent[kBeauty3P], keepEvent[kBeauty4P], keepEvent[kFemto2P], keepEvent[kFemto3P], keepEvent[kDoubleCharm2P], keepEvent[kDoubleCharm3P], keepEvent[kDoubleCharmMix], keepEvent[kV0Charm2P], keepEvent[kV0Charm3P], keepEvent[kCharmBarToXiBach], keepEvent[kSigmaCPPK], keepEvent[kSigmaC0K0], keepEvent[kPhotonCharm2P], keepEvent[kPhotonCharm3P], keepEvent[kSingleCharm2P], keepEvent[kSingleCharm3P], keepEvent[kSingleNonPromptCharm2P], keepEvent[kSingleNonPromptCharm3P], keepEvent[kCharmBarToXi2Bach], keepEvent[kPrCharm2P], keepEvent[kBtoJPsiKa], keepEvent[kBtoJPsiKstar], keepEvent[kBtoJPsiPhi], keepEvent[kBtoJPsiPrKa], keepEvent[kBtoJPsiPi]);
+                    tags(keepEvent[kHighPt2P], keepEvent[kHighPt3P], keepEvent[kBeauty3P], keepEvent[kBeauty4P], keepEvent[kFemto2P], keepEvent[kFemto3P], keepEvent[kDoubleCharm2P], keepEvent[kDoubleCharm3P], keepEvent[kDoubleCharmMix], keepEvent[kV0Charm2P], keepEvent[kV0Charm3P], keepEvent[kCharmBarToXiBach], keepEvent[kSigmaCPPK], keepEvent[kSigmaC0K0], keepEvent[kSigmaCP], keepEvent[kPhotonCharm2P], keepEvent[kPhotonCharm3P], keepEvent[kSingleCharm2P], keepEvent[kSingleCharm3P], keepEvent[kSingleNonPromptCharm2P], keepEvent[kSingleNonPromptCharm3P], keepEvent[kCharmBarToXi2Bach], keepEvent[kPrCharm2P], keepEvent[kBtoJPsiKa], keepEvent[kBtoJPsiKstar], keepEvent[kBtoJPsiPhi], keepEvent[kBtoJPsiPrKa], keepEvent[kBtoJPsiPi]);

Collaborator

fgrosa Feb 3, 2026

Suggested change

      
                  tags(keepEvent[kHighPt2P], keepEvent[kHighPt3P], keepEvent[kBeauty3P], keepEvent[kBeauty4P], keepEvent[kFemto2P], keepEvent[kFemto3P], keepEvent[kDoubleCharm2P], keepEvent[kDoubleCharm3P], keepEvent[kDoubleCharmMix], keepEvent[kV0Charm2P], keepEvent[kV0Charm3P], keepEvent[kCharmBarToXiBach], keepEvent[kSigmaCPPK], keepEvent[kSigmaC0K0], keepEvent[kSigmaCP], keepEvent[kPhotonCharm2P], keepEvent[kPhotonCharm3P], keepEvent[kSingleCharm2P], keepEvent[kSingleCharm3P], keepEvent[kSingleNonPromptCharm2P], keepEvent[kSingleNonPromptCharm3P], keepEvent[kCharmBarToXi2Bach], keepEvent[kPrCharm2P], keepEvent[kBtoJPsiKa], keepEvent[kBtoJPsiKstar], keepEvent[kBtoJPsiPhi], keepEvent[kBtoJPsiPrKa], keepEvent[kBtoJPsiPi]);
          
                  tags(keepEvent[kHighPt2P], keepEvent[kHighPt3P], keepEvent[kBeauty3P], keepEvent[kBeauty4P], keepEvent[kFemto2P], keepEvent[kFemto3P], keepEvent[kDoubleCharm2P], keepEvent[kDoubleCharm3P], keepEvent[kDoubleCharmMix], keepEvent[kV0Charm2P], keepEvent[kV0Charm3P], keepEvent[kCharmBarToXiBach], keepEvent[kSigmaCPPK], keepEvent[kSigmaC0K0], keepEvent[kSigmaCPr], keepEvent[kPhotonCharm2P], keepEvent[kPhotonCharm3P], keepEvent[kSingleCharm2P], keepEvent[kSingleCharm3P], keepEvent[kSingleNonPromptCharm2P], keepEvent[kSingleNonPromptCharm3P], keepEvent[kCharmBarToXi2Bach], keepEvent[kPrCharm2P], keepEvent[kBtoJPsiKa], keepEvent[kBtoJPsiKstar], keepEvent[kBtoJPsiPhi], keepEvent[kBtoJPsiPrKa], keepEvent[kBtoJPsiPi]);

EventFiltering/PWGHF/HFFilter.cxx Outdated

                     bool isSelectedPvZ = (std::fabs(collision.posZ()) < evSel.maxPvZ);
                     if (!isSelectedTVX || !isSelectedTFBorder || !isSelectedITSROFBorder || !isSelectedPvZ) {
-                      tags(keepEvent[kHighPt2P], keepEvent[kHighPt3P], keepEvent[kBeauty3P], keepEvent[kBeauty4P], keepEvent[kFemto2P], keepEvent[kFemto3P], keepEvent[kDoubleCharm2P], keepEvent[kDoubleCharm3P], keepEvent[kDoubleCharmMix], keepEvent[kV0Charm2P], keepEvent[kV0Charm3P], keepEvent[kCharmBarToXiBach], keepEvent[kSigmaCPPK], keepEvent[kSigmaC0K0], keepEvent[kPhotonCharm2P], keepEvent[kPhotonCharm3P], keepEvent[kSingleCharm2P], keepEvent[kSingleCharm3P], keepEvent[kSingleNonPromptCharm2P], keepEvent[kSingleNonPromptCharm3P], keepEvent[kCharmBarToXi2Bach], keepEvent[kPrCharm2P], keepEvent[kBtoJPsiKa], keepEvent[kBtoJPsiKstar], keepEvent[kBtoJPsiPhi], keepEvent[kBtoJPsiPrKa], keepEvent[kBtoJPsiPi]);
+                      tags(keepEvent[kHighPt2P], keepEvent[kHighPt3P], keepEvent[kBeauty3P], keepEvent[kBeauty4P], keepEvent[kFemto2P], keepEvent[kFemto3P], keepEvent[kDoubleCharm2P], keepEvent[kDoubleCharm3P], keepEvent[kDoubleCharmMix], keepEvent[kV0Charm2P], keepEvent[kV0Charm3P], keepEvent[kCharmBarToXiBach], keepEvent[kSigmaCPPK], keepEvent[kSigmaC0K0], keepEvent[kSigmaCP], keepEvent[kPhotonCharm2P], keepEvent[kPhotonCharm3P], keepEvent[kSingleCharm2P], keepEvent[kSingleCharm3P], keepEvent[kSingleNonPromptCharm2P], keepEvent[kSingleNonPromptCharm3P], keepEvent[kCharmBarToXi2Bach], keepEvent[kPrCharm2P], keepEvent[kBtoJPsiKa], keepEvent[kBtoJPsiKstar], keepEvent[kBtoJPsiPhi], keepEvent[kBtoJPsiPrKa], keepEvent[kBtoJPsiPi]);

Collaborator

fgrosa Feb 3, 2026

Suggested change

      
                    tags(keepEvent[kHighPt2P], keepEvent[kHighPt3P], keepEvent[kBeauty3P], keepEvent[kBeauty4P], keepEvent[kFemto2P], keepEvent[kFemto3P], keepEvent[kDoubleCharm2P], keepEvent[kDoubleCharm3P], keepEvent[kDoubleCharmMix], keepEvent[kV0Charm2P], keepEvent[kV0Charm3P], keepEvent[kCharmBarToXiBach], keepEvent[kSigmaCPPK], keepEvent[kSigmaC0K0], keepEvent[kSigmaCP], keepEvent[kPhotonCharm2P], keepEvent[kPhotonCharm3P], keepEvent[kSingleCharm2P], keepEvent[kSingleCharm3P], keepEvent[kSingleNonPromptCharm2P], keepEvent[kSingleNonPromptCharm3P], keepEvent[kCharmBarToXi2Bach], keepEvent[kPrCharm2P], keepEvent[kBtoJPsiKa], keepEvent[kBtoJPsiKstar], keepEvent[kBtoJPsiPhi], keepEvent[kBtoJPsiPrKa], keepEvent[kBtoJPsiPi]);
          
                    tags(keepEvent[kHighPt2P], keepEvent[kHighPt3P], keepEvent[kBeauty3P], keepEvent[kBeauty4P], keepEvent[kFemto2P], keepEvent[kFemto3P], keepEvent[kDoubleCharm2P], keepEvent[kDoubleCharm3P], keepEvent[kDoubleCharmMix], keepEvent[kV0Charm2P], keepEvent[kV0Charm3P], keepEvent[kCharmBarToXiBach], keepEvent[kSigmaCPPK], keepEvent[kSigmaC0K0], keepEvent[kSigmaCPr], keepEvent[kPhotonCharm2P], keepEvent[kPhotonCharm3P], keepEvent[kSingleCharm2P], keepEvent[kSingleCharm3P], keepEvent[kSingleNonPromptCharm2P], keepEvent[kSingleNonPromptCharm3P], keepEvent[kCharmBarToXi2Bach], keepEvent[kPrCharm2P], keepEvent[kBtoJPsiKa], keepEvent[kBtoJPsiKstar], keepEvent[kBtoJPsiPhi], keepEvent[kBtoJPsiPrKa], keepEvent[kBtoJPsiPi]);

EventFiltering/PWGHF/HFFilter.cxx

+                              float deltaEta = std::abs(RecoDecay::eta(pVecSigmaC) - track.eta() && pt3Prong < 3.0);
+                              if (configSigmaC.isTrigSigmaCP && isTrackProton && deltaEta < 1.0) {
+                                auto tagBDT = helper.isBDTSelected(scores[2], thresholdBDTScores[5]);

Collaborator

fgrosa Feb 3, 2026

be careful that this only works if thresholdBDTScores[5] is tighter than thresholdBDTScores[4]

Collaborator Author

singhra1994 Feb 3, 2026

I know, I assumed that BDT cut on 3 prong creators should be looser, I only need stronger BDT on sigmaCpr trigger if there is the problem in reaching required selectivity

EventFiltering/PWGHF/HFFilter.cxx Outdated

                               int8_t whichSigmaC = helper.isSelectedSigmaCInDeltaMassRange<2>(pVecFirst, pVecThird, pVecSecond, pVecSoftPi, ptSigmaC, is3Prong[2], hMassVsPtC[kNCharmParticles + 9], activateQA);
                               if (whichSigmaC > 0) {
-                                /// let's build a candidate SigmaC++K- pair
+                                /// let's build a candidate SigmaC++K- pair and SigmaC++,0K - proton

Collaborator

fgrosa Feb 3, 2026

This if statement selects only SigmaC in the signal mass range, while for SigmaC++,0K - proton you should keep also the sidebands, so it should not be included here

Collaborator Author

singhra1994 Feb 3, 2026

I forgot to remove this comment, it was useless

EventFiltering/PWGHF/HFFilter.cxx Outdated

-                                // check the invariant mass
-                                float massSigmaCPKPi{-999.}, massSigmaCPiKP{-999.}, deltaMassXicResoPKPi{-999.}, deltaMassXicResoPiKP{-999.};
-                                float ptSigmaCKaon = RecoDecay::pt(pVecSigmaC, pVecFourth);
+                                if (std::abs(chargeSc + track.sign()) == 1 && std::abs(chargeSc) == 2 && configSigmaC.isTrigSigmaCPPK) {

Collaborator

fgrosa Feb 3, 2026

Here you could add !keepEvent[kSigmaCPPK] && in the if statement (no need to further test if we found already a SigmaCK pair in this event)

EventFiltering/PWGHF/HFFilter.cxx Outdated

                                 }
                               }
+                              float deltaEta = std::abs(RecoDecay::eta(pVecSigmaC) - track.eta() && pt3Prong < 3.0);
+                              if (configSigmaC.isTrigSigmaCP && isTrackProton && deltaEta < 1.0) {

Collaborator

fgrosa Feb 3, 2026

Here you could add !keepEvent[kSigmaCPr] && in the if statement (no need to further test if we found already a SigmaCPr pair in this event)

EventFiltering/PWGHF/HFFilter.cxx


		auto tagBDT = helper.isBDTSelected(scores[2], thresholdBDTScores[5]);

		if (helper.selectionSigmaCForScPCorr(pVecFirst, pVecThird, pVecSecond, pVecSoftPi, ptSigmaC, is3Prong[2], hMassVsPtC[kNCharmParticles + 23], activateQA, configSigmaC.minMassSigmaCCorr, configSigmaC.maxMassSigmaCCorr, configSigmaC.minPtSigmaC, configSigmaC.maxPtSigmaC) && TESTBIT(tagBDT, RecoDecay::OriginType::Prompt)) {

Collaborator

fgrosa Feb 3, 2026

For this version it's ok, but please then open a second PR when this is merged and follow the convention as for the other triggers that the cuts are set as data members of the helper class

Collaborator Author

singhra1994 Feb 3, 2026

Sure thanks


          Refactor SigmaC trigger parameters and conditions

98059a1

alibuild mentioned this pull request

[Trigger,PWGHF] Please consider the following formatting changes to #14786 singhra1994/O2Physics#44

Merged

singhra1994 and others added 4 commits

February 3, 2026 19:57


          Refactor trigger names and update SigmaC checks

933661f


          Rename HfSigmaCP to HfSigmaCPr in filterTables.h

7f3053a


          Please consider the following formatting changes

ae91eee


          Merge pull request #44 from alibuild/alibot-cleanup-14786

20a1b16

[Trigger,PWGHF] Please consider the following formatting changes to AliceO2Group#14786

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Reviewers

fgrosa fgrosa requested changes

alibuild Awaiting requested review from alibuild alibuild is a code owner

zhangbiao-phy Awaiting requested review from zhangbiao-phy zhangbiao-phy is a code owner

mpuccio Awaiting requested review from mpuccio mpuccio is a code owner

lietava Awaiting requested review from lietava lietava is a code owner

ariedel-cern Awaiting requested review from ariedel-cern ariedel-cern is a code owner

Requested changes must be addressed to merge this pull request.

Labels