Polish Xic femto comments for HF guidelines

Author: FDUEnrich

PWGCF/DataModel/FemtoDerived.h

@@ -260,8 +260,8 @@ DECLARE_SOA_COLUMN(BDTBkg, bdtBkg, float);                       //! Background
 DECLARE_SOA_COLUMN(BDTPrompt, bdtPrompt, float);                 //! Prompt signal score using Boosted Decision Tree for charm hadron
 DECLARE_SOA_COLUMN(BDTFD, bdtFD, float);                         //! Feed-down score using Boosted Decision Tree for charm hadron
 DECLARE_SOA_COLUMN(CascBachelorTrackId, cascBachelorTrackId, int); //! Bachelor track ID from Xi cascade (Xic -> Xi pi pi)
-DECLARE_SOA_COLUMN(CascPosTrackId, cascPosTrackId, int);        //! Positive track ID from Λ in Ξ cascade (Ξc± → Ξππ)
-DECLARE_SOA_COLUMN(CascNegTrackId, cascNegTrackId, int);       //! Negative track ID from Λ in Ξ cascade (Ξc± → Ξππ)
+DECLARE_SOA_COLUMN(CascPosTrackId, cascPosTrackId, int);        //! Positive track ID from Lambda in Xi cascade (Xic -> Xi pi pi)
+DECLARE_SOA_COLUMN(CascNegTrackId, cascNegTrackId, int);       //! Negative track ID from Lambda in Xi cascade (Xic -> Xi pi pi)
 DECLARE_SOA_COLUMN(CascBachelorPt, cascBachelorPt, float);       //! pT of the bachelor track from the Xi cascade
 DECLARE_SOA_COLUMN(CascBachelorPhi, cascBachelorPhi, float);     //! phi of the bachelor track from the Xi cascade
 DECLARE_SOA_COLUMN(CascBachelorEta, cascBachelorEta, float);     //! eta of the bachelor track from the Xi cascade
@@ -271,7 +271,7 @@ DECLARE_SOA_COLUMN(CascPosEta, cascPosEta, float);               //! eta of the
 DECLARE_SOA_COLUMN(CascNegPt, cascNegPt, float);                 //! pT of the negative Lambda daughter track from the Xi cascade
 DECLARE_SOA_COLUMN(CascNegPhi, cascNegPhi, float);               //! phi of the negative Lambda daughter track from the Xi cascade
 DECLARE_SOA_COLUMN(CascNegEta, cascNegEta, float);               //! eta of the negative Lambda daughter track from the Xi cascade
-DECLARE_SOA_COLUMN(FlagMc, flagMc, int);                         //! To select MC particle among charm hadrons, { DplusToPiKPi = 1, LcToPKPi = 17, DsToKKPi = 6, XicToPKPi = 21, XicToXiPiPi = 1, N3ProngD = 2ecays };
+DECLARE_SOA_COLUMN(FlagMc, flagMc, int);                         //! MC matching flag for the selected charm hadron decay channel
 DECLARE_SOA_COLUMN(OriginMcRec, originMcRec, int);               //! flag for reconstruction level matching (1 for prompt, 2 for non-prompt)
 DECLARE_SOA_COLUMN(OriginMcGen, originMcGen, int);               //! flag for generator level matching (1 for prompt, 2 for non-prompt)
 DECLARE_SOA_COLUMN(IsCandidateSwapped, isCandidateSwapped, int); //! swapping of the prongs order (0 for Lc -> pkpi, 1 for Lc -> pikp)
@@ -283,12 +283,12 @@ DECLARE_SOA_COLUMN(KT, kT, float);                               //! kT distribu
 DECLARE_SOA_COLUMN(MT, mT, float);                               //! Transverse mass distribution
 DECLARE_SOA_COLUMN(CharmM, charmM, float);                       //! Charm hadron mass
 DECLARE_SOA_COLUMN(CharmDaughM, charmDaughM, float);             //! Charm hadron daughter mass
-DECLARE_SOA_COLUMN(CharmTrkM, charmtrkM, float);                 //! Charm hadron track mass
+DECLARE_SOA_COLUMN(CharmTrkM, charmtrkM, float);                 //! Invariant-mass difference of the charm-track pair
 DECLARE_SOA_COLUMN(CharmPt, charmPt, float);                     //! Transverse momentum of charm hadron for result task
 DECLARE_SOA_COLUMN(CharmEta, charmEta, float);                   //! Eta of charm hadron for result task
 DECLARE_SOA_COLUMN(CharmPhi, charmPhi, float);                   //! Phi of charm hadron for result task
 DECLARE_SOA_COLUMN(Mult, mult, int);                             //! Charge particle multiplicity
-DECLARE_SOA_COLUMN(MultPercentile, multPercentile, float);       //! Multiplicity precentile
+DECLARE_SOA_COLUMN(MultPercentile, multPercentile, float);       //! Multiplicity percentile
 DECLARE_SOA_COLUMN(PairSign, pairSign, int8_t);                  //! Selection between like sign (1) and unlike sign pair (2)
 DECLARE_SOA_COLUMN(ProcessType, processType, int64_t);           //! Selection between same-event (1), and mixed-event (2)
 DECLARE_SOA_DYNAMIC_COLUMN(M, m,                                 //!

PWGHF/HFC/TableProducer/producerCharmHadronsTrackFemtoDream.cxx

@@ -146,8 +146,6 @@ struct HfProducerCharmHadronsTrackFemtoDream {
     Configurable<bool> loadModelsFromCCDB{"loadModelsFromCCDB", false, "Flag to enable or disable the loading of models from CCDB"};
   } ccdbCfg;
 
-  // Configurable<bool> isForceGRP{"isForceGRP", false, "Set true if the magnetic field configuration is not available in the usual CCDB directory (e.g. for Run 2 converted data or unanchorad Monte Carlo)"};
-
   // -------------------------
   // Kaon PID cut parameters
   // -------------------------
@@ -169,7 +167,7 @@ struct HfProducerCharmHadronsTrackFemtoDream {
   struct : ConfigurableGroup {
     Configurable<bool> isDebug{"isDebug", true, "Enable Debug tables"};
     Configurable<bool> isRun3{"isRun3", true, "Running on Run3 or pilot"};
-    Configurable<int> selectionFlagHadron{"selectionFlagHadron", 1, "Selection Flag for Charm Hadron: 1 for Lc, 7 for Dplus (Topologic and PID cuts)"};
+    Configurable<int> selectionFlagHadron{"selectionFlagHadron", 1, "Selection flag for charm hadrons; applied to D0, Dplus, Lc and Xic selector decisions"};
     Configurable<bool> useCent{"useCent", false, "Enable centrality for Charm Hadron"};
   } generalCfg;
 
@@ -221,7 +219,6 @@ struct HfProducerCharmHadronsTrackFemtoDream {
   o2::hf_evsel::HfEventSelection hfEvSel;
   Service<o2::ccdb::BasicCCDBManager> ccdb{}; /// Accessing the CCDB
   o2::base::MatLayerCylSet* lut{};
-  // if (doPvRefit){ lut = o2::base::MatLayerCylSet::rectifyPtrFromFile(ccdb->get<o2::base::MatLayerCylSet>(ccdbCfg.ccdbPathLut));} //! may be it useful, will check later
 
   float magField{};
   int runNumber{};

PWGHF/HFC/Tasks/taskCharmHadronsTrackFemtoDream.cxx

@@ -108,7 +108,7 @@ struct HfTaskCharmHadronsTrackFemtoDream {
   struct : ConfigurableGroup {
     Configurable<float> charmHadBkgBDTmax{"charmHadBkgBDTmax", 1., "Maximum background bdt score for Charm Hadron (particle 2)"};
     Configurable<int> charmHadCandSel{"charmHadCandSel", 1, "candidate selection for charm hadron"};
-    Configurable<int> charmHadMcSel{"charmHadMcSel", DecayChannelMain::LcToPKPi, "charm hadron selection for mc, DplusToPiKPi = 1, XicToXiPiPi = 1, LcToPKPi = 17"};
+    Configurable<int> charmHadMcSel{"charmHadMcSel", DecayChannelMain::LcToPKPi, "MC matching flag for the selected charm hadron decay channel"};
     Configurable<float> charmHadFdBDTmin{"charmHadFdBDTmin", 0., "Minimum feed-down bdt score Charm Hadron (particle 2)"};
     Configurable<float> charmHadFdBDTmax{"charmHadFdBDTmax", 1., "Maximum feed-down bdt score Charm Hadron (particle 2)"};
     Configurable<float> charmHadMaxInvMass{"charmHadMaxInvMass", 2.45, "Maximum invariant mass of Charm Hadron (particle 2)"};
@@ -375,18 +375,18 @@ struct HfTaskCharmHadronsTrackFemtoDream {
       }
       invMass = cand.m(std::array{MassPiPlus, MassKPlus, MassProton});
       return invMass;
-    } else if constexpr (Channel == DecayChannel::DplusToPiKPi) { // D+ → π K π (PDG: 411)
+    } else if constexpr (Channel == DecayChannel::DplusToPiKPi) { // D+ -> pi K pi
       invMass = cand.m(std::array{MassPiPlus, MassKPlus, MassPiPlus});
       return invMass;
-    } else if constexpr (Channel == DecayChannel::D0ToPiK) { // D0 → π K  (PDG: 421)
+    } else if constexpr (Channel == DecayChannel::D0ToPiK) { // D0 -> pi K
       if (cand.candidateSelFlag() == 1) {
         invMass = cand.m(std::array{MassPiPlus, MassKPlus});
         return invMass;
       } else {
         invMass = cand.m(std::array{MassKPlus, MassPiPlus});
         return invMass;
       }
-    } else if constexpr (Channel == DecayChannel::DstarToD0Pi) { // D* → D0π (PDG: 413)
+    } else if constexpr (Channel == DecayChannel::DstarToD0Pi) { // D* -> D0 pi
       float mDstar = 0.f;
       float mD0 = 0.f;
       if (cand.charge() > 0.f) {
@@ -436,7 +436,7 @@ struct HfTaskCharmHadronsTrackFemtoDream {
         LOG(fatal) << "Invalid PDG code for track mass hypothesis: " << trkPDGCode;
     }
 
-    // D0 → K π + track (2-prong)
+      // D0 -> K pi + track (2-prong)
     if constexpr (Channel == DecayChannel::D0ToPiK) {
       const auto pVecCharmTrk = std::array{pVecProng0, pVecProng1, pVecTrack};
       std::array<double, 3> massCharmTrk{};
@@ -465,17 +465,17 @@ struct HfTaskCharmHadronsTrackFemtoDream {
       std::array<double, 4> massCharmTrk{};
 
       if constexpr (Channel == DecayChannel::LcToPKPi) {
-        // Λc⁺ → p K π
+        // Lc+ -> p K pi
         if (cand.candidateSelFlag() == 1) {
           massCharmTrk = {MassProton, MassKPlus, MassPiPlus, trackMassHyp};
         } else {
           massCharmTrk = {MassPiPlus, MassKPlus, MassProton, trackMassHyp};
         }
       } else if constexpr (Channel == DecayChannel::DplusToPiKPi) {
-        // D⁺ → π K π
+        // D+ -> pi K pi
         massCharmTrk = {MassPiPlus, MassKPlus, MassPiPlus, trackMassHyp};
       } else if constexpr (Channel == DecayChannel::DstarToD0Pi) {
-        // D* → D0π
+        // D* -> D0 pi
         if (cand.candidateSelFlag() == 1) {
           massCharmTrk = {MassPiPlus, MassKPlus, MassPiPlus, trackMassHyp};
         } else {
@@ -571,7 +571,7 @@ struct HfTaskCharmHadronsTrackFemtoDream {
 
       float deltaInvMassPair = getCharmHadronTrackMass<Channel>(p2, p1, trackSel.pdgCodeTrack1.value) - invMass;
 
-      // proton track charge
+      // associated track charge
       float chargeTrack = 0.;
       if ((p1.cut() & CutBitChargePositive) == CutBitChargePositive) {
         chargeTrack = PositiveCharge;
@@ -708,7 +708,7 @@ struct HfTaskCharmHadronsTrackFemtoDream {
 
         float deltaInvMassPair = getCharmHadronTrackMass<Channel>(p2, p1, trackSel.pdgCodeTrack1.value) - invMass;
 
-        // proton track charge
+        // associated track charge
         float chargeTrack = 0.;
         if ((p1.cut() & CutBitChargePositive) == CutBitChargePositive) {
           chargeTrack = PositiveCharge;