[PWGHF] tree creator for Sigmac correlated bkg. sources.

PWGHF/D2H/Tasks/taskSigmac.cxx

@@ -17,6 +17,7 @@
 
 #include "PWGHF/Core/DecayChannels.h"
 #include "PWGHF/Core/HfHelper.h"
+#include "PWGHF/D2H/Utils/utilsSigmac.h"
 #include "PWGHF/DataModel/CandidateReconstructionTables.h"
 #include "PWGHF/DataModel/CandidateSelectionTables.h"
 
@@ -340,27 +341,6 @@ struct HfTaskSigmac {
 
   }; /// end init
 
-  /// @brief Function to determine if the reconstructed candidate Σc0,++ decays into Λc+ → pK-π+, Λc+ → π+K-p or both
-  /// @tparam L template for Lc daughter of Sc candidate
-  /// @tparam S template for Sc candidate
-  /// @param candidateLc Lc daughter of Sc candidate
-  /// @param candSc Sc candidate
-  /// @return 0: none; 1: only Λc+ → pK-π+ possible; 2: Λc+ → π+K-p possible; 3: both possible
-  template <typename L, typename S>
-  int8_t isDecayToPKPiToPiKP(L& candidateLc, S& candSc)
-  {
-    int8_t channel = 0;
-    if ((candidateLc.isSelLcToPKPi() >= 1) && candSc.statusSpreadLcMinvPKPiFromPDG()) {
-      // Λc+ → pK-π+ and within the requested mass to build the Σc0,++
-      SETBIT(channel, o2::aod::hf_cand_sigmac::Decays::PKPi);
-    }
-    if ((candidateLc.isSelLcToPiKP() >= 1) && candSc.statusSpreadLcMinvPiKPFromPDG()) {
-      // Λc+ → π+K-p and within the requested mass to build the Σc0,++
-      SETBIT(channel, o2::aod::hf_cand_sigmac::Decays::PiKP);
-    }
-    return channel; /// 0: none; 1: pK-π+ only; 2: π+K-p only; 3: both possible
-  }
-
   /// @brief function to fill the histograms needed in analysis (data)
   /// @param candidatesSc are the reconstructed candidate Σc0,++
   /// @param
@@ -386,7 +366,7 @@ struct HfTaskSigmac {
       const auto& candidateLc = candSc.prongLc_as<CandsLc>();
       // const int iscandidateLcpKpi = (candidateLc.isSelLcToPKPi() >= 1) && candSc.statusSpreadLcMinvPKPiFromPDG(); // Λc+ → pK-π+ and within the requested mass to build the Σc0,++
       // const int iscandidateLcpiKp = (candidateLc.isSelLcToPiKP() >= 1) && candSc.statusSpreadLcMinvPiKPFromPDG(); // Λc+ → π+K-p and within the requested mass to build the Σc0,++
-      const int8_t isCandPKPiPiKP = isDecayToPKPiToPiKP(candidateLc, candSc);
+      const int8_t isCandPKPiPiKP = hf_sigmac_utils::isDecayToPKPiToPiKP(candidateLc, candSc);
       double massSc(-1.), massLc(-1.), deltaMass(-1.);
       double ptSc(candSc.pt()), ptLc(candidateLc.pt());
       double etaSc(candSc.eta()), etaLc(candidateLc.eta());
@@ -822,7 +802,7 @@ struct HfTaskSigmac {
       /// get the candidate Λc+ used to build the Σc0
       /// and understand which mass hypotheses are possible
       const auto& candidateLc = candSc.prongLc_as<CandsLc>();
-      const int8_t isCandPKPiPiKP = isDecayToPKPiToPiKP(candidateLc, candSc);
+      const int8_t isCandPKPiPiKP = hf_sigmac_utils::isDecayToPKPiToPiKP(candidateLc, candSc);
 
       // candidateLc.flagMcDecayChanRec();
 

PWGHF/D2H/Utils/utilsSigmac.h

@@ -0,0 +1,47 @@
+// Copyright 2019-2020 CERN and copyright holders of ALICE O2.
+// See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.
+// All rights not expressly granted are reserved.
+//
+// This software is distributed under the terms of the GNU General Public
+// License v3 (GPL Version 3), copied verbatim in the file "COPYING".
+//
+// In applying this license CERN does not waive the privileges and immunities
+// granted to it by virtue of its status as an Intergovernmental Organization
+// or submit itself to any jurisdiction.
+
+/// \file utilsSigmac.h
+/// \brief Utilities for Sigmac analysis
+/// \author Mattia Faggin <mattia.faggin@cern.ch>, INFN Padova
+
+#ifndef PWGHF_D2H_UTILS_UTILSSIGMAC_H_
+#define PWGHF_D2H_UTILS_UTILSSIGMAC_H_
+
+#include "PWGHF/DataModel/CandidateReconstructionTables.h"
+
+#include <cstdint>
+
+namespace o2::hf_sigmac_utils
+{
+/// @brief Function to determine if the reconstructed candidate Σc0,++ decays into Λc+ → pK-π+, Λc+ → π+K-p or both
+/// @tparam L template for Lc daughter of Sc candidate
+/// @tparam S template for Sc candidate
+/// @param candidateLc Lc daughter of Sc candidate
+/// @param candSc Sc candidate
+/// @return 0: none; 1: only Λc+ → pK-π+ possible; 2: Λc+ → π+K-p possible; 3: both possible
+template <typename L, typename S>
+int8_t isDecayToPKPiToPiKP(L& candidateLc, S& candSc)
+{
+  int8_t channel = 0;
+  if ((candidateLc.isSelLcToPKPi() >= 1) && candSc.statusSpreadLcMinvPKPiFromPDG()) {
+    // Λc+ → pK-π+ and within the requested mass to build the Σc0,++
+    SETBIT(channel, o2::aod::hf_cand_sigmac::Decays::PKPi);
+  }
+  if ((candidateLc.isSelLcToPiKP() >= 1) && candSc.statusSpreadLcMinvPiKPFromPDG()) {
+    // Λc+ → π+K-p and within the requested mass to build the Σc0,++
+    SETBIT(channel, o2::aod::hf_cand_sigmac::Decays::PiKP);
+  }
+  return channel; /// 0: none; 1: pK-π+ only; 2: π+K-p only; 3: both possible
+}
+} // namespace o2::hf_sigmac_utils
+
+#endif // PWGHF_D2H_UTILS_UTILSSIGMAC_H_

PWGHF/TableProducer/CMakeLists.txt

@@ -249,6 +249,11 @@ o2physics_add_dpl_workflow(tree-creator-lc-to-p-k-pi
                     PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore
                     COMPONENT_NAME Analysis)
 
+o2physics_add_dpl_workflow(tree-creator-sigmac-corr-bkg
+                    SOURCES treeCreatorSigmacCorrBkg.cxx
+                    PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore
+                    COMPONENT_NAME Analysis)
+
 o2physics_add_dpl_workflow(tree-creator-omegac-st
                     SOURCES treeCreatorOmegacSt.cxx
                     PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore O2::DCAFitter O2Physics::EventFilteringUtils