PWGHF: first checked implementation of Xic->pKpi filtering (#5436)

* PWGHF: first checked code of Xic->pKpi filtering

* Update HFSecondaryVertex.h

* Update HFXicCandidateSelector.cxx

* Fix code with camelCase syntax conventions

* Fix code syntax conventions for Xic methods

* PWGHF: further fixes for Xic to pKpi filtering and analysis

 - (HFCandidateCreator3Prong.cxx) remove refuse std::move, forgot in the middle of the for loop
 - (HFCandidateCreator3Prong.cxx) fix if statement condition for isMatchedMCGen
 - (taskXic.cxx) fix MC histogram titles

* PWGHF: avoid storing return values of RecoDecay::getMatchedMCRec and isMatchedMCGen for Xic->pKpi

* PWGHF: pair Xic code with adaptAnalysisTask and charge() changes

* PWGHF: remove check for charge to be different from 0

* PWGHF: change label_as into mcParticle_as in taskXic

* PWGHF: fix particle label for ct distribution in taskXic

* PWGHF: remove selection flag for Xicbar (obsolete) from taskXic

Co-authored-by: Mattia Faggin <mfaggin@cern.ch>

Author: mfaggin

Analysis/Core/include/AnalysisCore/HFConfigurables.h

@@ -54,6 +54,12 @@ class HFTrackIndexSkimsCreatorConfigs
   double mInvMassDsToPiKKMax = 2.15; //original value 2.2
   double mCPADsToPiKKMin = 0.5;
   double mDecLenDsToPiKKMin = 0.;
+  // 3-prong cuts - XicToPKPi
+  double mPtXicToPKPiMin = 1.;        //
+  double mInvMassXicToPKPiMin = 2.25; //
+  double mInvMassXicToPKPiMax = 2.70; //
+  double mCPAXicToPKPiMin = 0.5;
+  double mDecLenXicToPKPiMin = 0.;
 
  private:
   ClassDef(HFTrackIndexSkimsCreatorConfigs, 1);

Analysis/DataModel/include/AnalysisDataModel/HFCandidateSelectionTables.h

@@ -39,4 +39,14 @@ DECLARE_SOA_COLUMN(IsSelJpsiToEE, isSelJpsiToEE, int);
 } // namespace hf_selcandidate_jpsi
 DECLARE_SOA_TABLE(HFSelJpsiToEECandidate, "AOD", "HFSELJPSICAND", hf_selcandidate_jpsi::IsSelJpsiToEE);
 } // namespace o2::aod
+
+namespace o2::aod
+{
+namespace hf_selcandidate_xic
+{
+DECLARE_SOA_COLUMN(IsSelXicToPKPi, isSelXicToPKPi, int);
+DECLARE_SOA_COLUMN(IsSelXicToPiKP, isSelXicToPiKP, int);
+} // namespace hf_selcandidate_xic
+DECLARE_SOA_TABLE(HFSelXicToPKPiCandidate, "AOD", "HFSELXICCAND", hf_selcandidate_xic::IsSelXicToPKPi, hf_selcandidate_xic::IsSelXicToPiKP);
+} // namespace o2::aod
 #endif // O2_ANALYSIS_HFCANDIDATESELECTIONTABLES_H_

Analysis/DataModel/include/AnalysisDataModel/HFSecondaryVertex.h

@@ -55,6 +55,7 @@ DECLARE_SOA_COLUMN(JpsiToEEFlag, jpsiToEEFlag, uint8_t);
 DECLARE_SOA_COLUMN(DPlusPiKPiFlag, dPlusPiKPiFlag, uint8_t);
 DECLARE_SOA_COLUMN(LcPKPiFlag, lcPKPiFlag, uint8_t);
 DECLARE_SOA_COLUMN(DsKKPiFlag, dsKKPiFlag, uint8_t);
+DECLARE_SOA_COLUMN(XicToPKPiFlag, xicToPKPiFlag, uint8_t);
 } // namespace hf_track_index
 
 DECLARE_SOA_TABLE(HfTrackIndexProng2, "AOD", "HFTRACKIDXP2",
@@ -75,7 +76,8 @@ DECLARE_SOA_TABLE(HfTrackIndexProng3, "AOD", "HFTRACKIDXP3",
 DECLARE_SOA_TABLE(HfCutStatusProng3, "AOD", "HFCUTSTATUSP3",
                   hf_track_index::DPlusPiKPiFlag,
                   hf_track_index::LcPKPiFlag,
-                  hf_track_index::DsKKPiFlag);
+                  hf_track_index::DsKKPiFlag,
+                  hf_track_index::XicToPKPiFlag);
 
 // general decay properties
 namespace hf_cand
@@ -312,6 +314,7 @@ DECLARE_SOA_DYNAMIC_COLUMN(MaxNormalisedDeltaIP, maxNormalisedDeltaIP, [](float
 // MC matching result:
 // - ±DPlusToPiKPi: D± → π± K∓ π±
 // - ±LcToPKPi: Λc± → p± K∓ π±
+// - ±XicToPKPi: Ξc± → p± K∓ π±
 DECLARE_SOA_COLUMN(FlagMCMatchRec, flagMCMatchRec, int8_t);         // reconstruction level
 DECLARE_SOA_COLUMN(FlagMCMatchGen, flagMCMatchGen, int8_t);         // generator level
 DECLARE_SOA_COLUMN(FlagMCDecayChanRec, flagMCDecayChanRec, int8_t); // Resonant decay channel flag, reconstruction level
@@ -321,6 +324,7 @@ DECLARE_SOA_COLUMN(FlagMCDecayChanGen, flagMCDecayChanGen, int8_t); // Resonant
 enum DecayType { DPlusToPiKPi = 0,
                  LcToPKPi,
                  DsToPiKK,
+                 XicToPKPi,
                  N3ProngDecays }; //always keep N3ProngDecays at the end
 
 // functions for specific particles
@@ -382,6 +386,38 @@ auto InvMassLcpiKp(const T& candidate)
 {
   return candidate.m(array{RecoDecay::getMassPDG(kPiPlus), RecoDecay::getMassPDG(kKPlus), RecoDecay::getMassPDG(kProton)});
 }
+
+// Ξc± → p± K∓ π±
+
+template <typename T>
+auto CtXic(const T& candidate)
+{
+  return candidate.ct(RecoDecay::getMassPDG(4232));
+}
+
+template <typename T>
+auto YXic(const T& candidate)
+{
+  return candidate.y(RecoDecay::getMassPDG(4232));
+}
+
+template <typename T>
+auto EXic(const T& candidate)
+{
+  return candidate.e(RecoDecay::getMassPDG(4232));
+}
+
+template <typename T>
+auto InvMassXicToPKPi(const T& candidate)
+{
+  return InvMassLcpKpi(candidate);
+}
+
+template <typename T>
+auto InvMassXicToPiKP(const T& candidate)
+{
+  return InvMassLcpiKp(candidate);
+}
 } // namespace hf_cand_prong3
 
 // 3-prong decay candidate table

Analysis/Tasks/PWGHF/CMakeLists.txt

@@ -68,6 +68,11 @@ o2_add_dpl_workflow(hf-jpsi-toee-candidate-selector
                     PUBLIC_LINK_LIBRARIES O2::Framework O2::AnalysisDataModel O2::AnalysisCore O2::DetectorsVertexing
                     COMPONENT_NAME Analysis)
 
+o2_add_dpl_workflow(hf-xic-topkpi-candidate-selector
+                    SOURCES HFXicToPKPiCandidateSelector.cxx
+                    PUBLIC_LINK_LIBRARIES O2::Framework O2::AnalysisDataModel O2::AnalysisCore O2::DetectorsVertexing
+                    COMPONENT_NAME Analysis)
+
 o2_add_dpl_workflow(hf-task-d0
                     SOURCES taskD0.cxx
                     PUBLIC_LINK_LIBRARIES O2::Framework O2::AnalysisDataModel O2::AnalysisCore O2::DetectorsVertexing
@@ -93,6 +98,11 @@ o2_add_dpl_workflow(hf-task-bplus
                     PUBLIC_LINK_LIBRARIES O2::Framework O2::AnalysisDataModel O2::AnalysisCore O2::DetectorsVertexing
                     COMPONENT_NAME Analysis)
 
+o2_add_dpl_workflow(hf-task-xic
+                    SOURCES taskXic.cxx
+                    PUBLIC_LINK_LIBRARIES O2::Framework O2::AnalysisDataModel O2::AnalysisCore O2::DetectorsVertexing
+                    COMPONENT_NAME Analysis)
+
 o2_add_dpl_workflow(hf-mc-validation
                     SOURCES HFMCValidation.cxx
                     PUBLIC_LINK_LIBRARIES O2::Framework O2::AnalysisDataModel O2::AnalysisCore O2::DetectorsVertexing

Analysis/Tasks/PWGHF/HFCandidateCreator3Prong.cxx

@@ -179,7 +179,7 @@ struct HFCandidateCreator3ProngMC {
       // Λc± → p± K∓ π±
       if (flag == 0) {
         //Printf("Checking Λc± → p± K∓ π±");
-        auto indexRecLc = RecoDecay::getMatchedMCRec(particlesMC, std::move(arrayDaughters), 4122, array{+kProton, -kKPlus, +kPiPlus}, true, &sign, 2);
+        auto indexRecLc = RecoDecay::getMatchedMCRec(particlesMC, arrayDaughters, 4122, array{+kProton, -kKPlus, +kPiPlus}, true, &sign, 2);
         if (indexRecLc > -1) {
           flag = sign * (1 << LcToPKPi);
 
@@ -201,6 +201,14 @@ struct HFCandidateCreator3ProngMC {
         }
       }
 
+      // Ξc± → p± K∓ π±
+      if (flag == 0) {
+        //Printf("Checking Ξc± → p± K∓ π±");
+        if (RecoDecay::getMatchedMCRec(particlesMC, std::move(arrayDaughters), 4232, array{+kProton, -kKPlus, +kPiPlus}, true, &sign) > -1) {
+          flag = sign * (1 << XicToPKPi);
+        }
+      }
+
       rowMCMatchRec(flag, DecayChannel);
     }
 
@@ -242,6 +250,14 @@ struct HFCandidateCreator3ProngMC {
         }
       }
 
+      // Ξc± → p± K∓ π±
+      if (flag == 0) {
+        //Printf("Checking Ξc± → p± K∓ π±");
+        if (RecoDecay::isMatchedMCGen(particlesMC, particle, 4232, array{+kProton, -kKPlus, +kPiPlus}, true, &sign)) {
+          flag = sign * (1 << XicToPKPi);
+        }
+      }
+
       rowMCMatchGen(flag, DecayChannel);
     }
   }

Analysis/Tasks/PWGHF/HFTrackIndexSkimsCreator.cxx

@@ -183,7 +183,8 @@ struct HFTrackIndexSkimsCreator {
      {"hNCand3ProngVsNTracks", "3-prong candidates preselected;# of selected tracks;# of candidates;entries", {HistType::kTH2F, {{2500, 0., 25000.}, {5000, 0., 500000.}}}},
      {"hmassDPlusToPiKPi", "D+ candidates;inv. mass (#pi K #pi) (GeV/#it{c}^{2});entries", {HistType::kTH1F, {{500, 0., 5.}}}},
      {"hmassLcToPKPi", "Lc candidates;inv. mass (p K #pi) (GeV/#it{c}^{2});entries", {HistType::kTH1F, {{500, 0., 5.}}}},
-     {"hmassDsToPiKK", "Ds candidates;inv. mass (K K #pi) (GeV/#it{c}^{2});entries", {HistType::kTH1F, {{500, 0., 5.}}}}}};
+     {"hmassDsToPiKK", "Ds candidates;inv. mass (K K #pi) (GeV/#it{c}^{2});entries", {HistType::kTH1F, {{500, 0., 5.}}}},
+     {"hmassXicToPKPi", "Xic candidates;inv. mass (p K #pi) (GeV/#it{c}^{2});entries", {HistType::kTH1F, {{500, 0., 5.}}}}}};
 
   Filter filterSelectTracks = (aod::hf_seltrack::isSelProng > 0);
   using SelectedTracks = soa::Filtered<soa::Join<aod::Tracks, aod::TracksCov, aod::TracksExtra, aod::HFSelTrack>>;
@@ -278,6 +279,12 @@ struct HFTrackIndexSkimsCreator {
     cut3ProngCPACandMin[DsToPiKK] = configs->mCPADsToPiKKMin;
     cut3ProngDecLenCandMin[DsToPiKK] = configs->mDecLenDsToPiKKMin;
 
+    cut3ProngPtCandMin[XicToPKPi] = configs->mPtXicToPKPiMin;
+    cut3ProngInvMassCandMin[XicToPKPi] = configs->mInvMassXicToPKPiMin;
+    cut3ProngInvMassCandMax[XicToPKPi] = configs->mInvMassXicToPKPiMax;
+    cut3ProngCPACandMin[XicToPKPi] = configs->mCPAXicToPKPiMin;
+    cut3ProngDecLenCandMin[XicToPKPi] = configs->mDecLenXicToPKPiMin;
+
     bool cutStatus2Prong[n2ProngDecays][nCuts2Prong];
     bool cutStatus3Prong[n3ProngDecays][nCuts3Prong];
     int nCutStatus2ProngBit = TMath::Power(2, nCuts2Prong) - 1; //bit value for selection status for each 2 prongs candidate where each selection is one bit and they are all set it 1
@@ -295,11 +302,13 @@ struct HFTrackIndexSkimsCreator {
     arr3Mass1[DPlusToPiKPi] = array{massPi, massK, massPi};
     arr3Mass1[LcToPKPi] = array{massProton, massK, massPi};
     arr3Mass1[DsToPiKK] = array{massK, massK, massPi};
+    arr3Mass1[XicToPKPi] = array{massProton, massK, massPi};
 
     array<array<double, 3>, n3ProngDecays> arr3Mass2;
     arr3Mass2[DPlusToPiKPi] = array{massPi, massK, massPi};
     arr3Mass2[LcToPKPi] = array{massPi, massK, massProton};
     arr3Mass2[DsToPiKK] = array{massPi, massK, massK};
+    arr3Mass2[XicToPKPi] = array{massPi, massK, massProton};
 
     double mass2ProngHypo1[n2ProngDecays];
     double mass2ProngHypo2[n2ProngDecays];
@@ -631,7 +640,7 @@ struct HFTrackIndexSkimsCreator {
                   }
                 }
               }
-              rowProng3CutStatus(Prong3CutStatus[0], Prong3CutStatus[1], Prong3CutStatus[2]); //FIXME when we can do this by looping over n3ProngDecays
+              rowProng3CutStatus(Prong3CutStatus[0], Prong3CutStatus[1], Prong3CutStatus[2], Prong3CutStatus[3]); //FIXME when we can do this by looping over n3ProngDecays
             }
 
             // fill histograms
@@ -654,6 +663,9 @@ struct HFTrackIndexSkimsCreator {
                     if (n3 == DsToPiKK) {
                       registry.get<TH1>(HIST("hmassDsToPiKK"))->Fill(mass3ProngHypo1[n3]);
                     }
+                    if (n3 == XicToPKPi) {
+                      registry.get<TH1>(HIST("hmassXicToPKPi"))->Fill(mass3ProngHypo1[n3]);
+                    }
                   }
                   if ((cut3ProngInvMassCandMin[n3] < 0. && cut3ProngInvMassCandMax[n3] <= 0.) || (mass3ProngHypo2[n3] >= cut3ProngInvMassCandMin[n3] && mass3ProngHypo2[n3] < cut3ProngInvMassCandMax[n3])) {
                     mass3ProngHypo2[n3] = RecoDecay::M(arr3Mom, arr3Mass2[n3]);
@@ -663,6 +675,9 @@ struct HFTrackIndexSkimsCreator {
                     if (n3 == DsToPiKK) {
                       registry.get<TH1>(HIST("hmassDsToPiKK"))->Fill(mass3ProngHypo2[n3]);
                     }
+                    if (n3 == XicToPKPi) {
+                      registry.get<TH1>(HIST("hmassXicToPKPi"))->Fill(mass3ProngHypo2[n3]);
+                    }
                   }
                 }
               }
@@ -798,7 +813,7 @@ struct HFTrackIndexSkimsCreator {
                   }
                 }
               }
-              rowProng3CutStatus(Prong3CutStatus[0], Prong3CutStatus[1], Prong3CutStatus[2]); //FIXME when we can do this by looping over n3ProngDecays
+              rowProng3CutStatus(Prong3CutStatus[0], Prong3CutStatus[1], Prong3CutStatus[2], Prong3CutStatus[3]); //FIXME when we can do this by looping over n3ProngDecays
             }
 
             // fill histograms
@@ -821,6 +836,9 @@ struct HFTrackIndexSkimsCreator {
                     if (n3 == DsToPiKK) {
                       registry.get<TH1>(HIST("hmassDsToPiKK"))->Fill(mass3ProngHypo1[n3]);
                     }
+                    if (n3 == XicToPKPi) {
+                      registry.get<TH1>(HIST("hmassXicToPKPi"))->Fill(mass3ProngHypo1[n3]);
+                    }
                   }
                   if ((cut3ProngInvMassCandMin[n3] < 0. && cut3ProngInvMassCandMax[n3] <= 0.) || (mass3ProngHypo2[n3] >= cut3ProngInvMassCandMin[n3] && mass3ProngHypo2[n3] < cut3ProngInvMassCandMax[n3])) {
                     mass3ProngHypo2[n3] = RecoDecay::M(arr3Mom, arr3Mass2[n3]);
@@ -830,6 +848,9 @@ struct HFTrackIndexSkimsCreator {
                     if (n3 == DsToPiKK) {
                       registry.get<TH1>(HIST("hmassDsToPiKK"))->Fill(mass3ProngHypo2[n3]);
                     }
+                    if (n3 == XicToPKPi) {
+                      registry.get<TH1>(HIST("hmassXicToPKPi"))->Fill(mass3ProngHypo2[n3]);
+                    }
                   }
                 }
               }