Add hypertriton kink analysis into task

Author: wang-yuanzhe

PWGLF/DataModel/LFHyperNucleiKinkTables.h

@@ -24,6 +24,7 @@ namespace o2::aod
 
 namespace hyperkink
 {
+DECLARE_SOA_COLUMN(MagPolarity, magPolarity, int8_t);                   //! Magnetic field polarity
 DECLARE_SOA_COLUMN(XPV, xPV, float);                                    //! Primary vertex of the candidate (x direction)
 DECLARE_SOA_COLUMN(YPV, yPV, float);                                    //! Primary vertex of the candidate (y direction)
 DECLARE_SOA_COLUMN(ZPV, zPV, float);                                    //! Primary vertex of the candidate (z direction)
@@ -82,6 +83,7 @@ DECLARE_SOA_COLUMN(UpdatePzMothPV, updatePzMothPV, float);      //! updated pz o
 
 DECLARE_SOA_TABLE(HypKinkCand, "AOD", "HYPKINKCANDS",
                   o2::soa::Index<>,
+                  hyperkink::MagPolarity,
                   hyperkink::XPV, hyperkink::YPV, hyperkink::ZPV,
                   hyperkink::XSV, hyperkink::YSV, hyperkink::ZSV,
                   hyperkink::IsMatter,
@@ -96,6 +98,7 @@ DECLARE_SOA_TABLE(HypKinkCand, "AOD", "HYPKINKCANDS",
 
 DECLARE_SOA_TABLE(MCHypKinkCand, "AOD", "MCHYPKINKCANDS",
                   o2::soa::Index<>,
+                  hyperkink::MagPolarity,
                   hyperkink::XPV, hyperkink::YPV, hyperkink::ZPV,
                   hyperkink::XSV, hyperkink::YSV, hyperkink::ZSV,
                   hyperkink::IsMatter,

PWGLF/TableProducer/Common/kinkBuilder.cxx

@@ -91,6 +91,7 @@ struct kinkCandidate {
 struct kinkBuilder {
 
   enum PartType { kSigmaMinus = 0,
+                  kHypertriton,
                   kHyperhelium4sigma };
 
   Produces<aod::KinkCands> outputDataTable;
@@ -161,6 +162,11 @@ struct kinkBuilder {
       mothMass = o2::constants::physics::MassSigmaMinus;
       chargedDauMass = o2::constants::physics::MassPiMinus;
       neutDauMass = o2::constants::physics::MassNeutron;
+    } else if (hypoMoth == kHypertriton) {
+      charge = 1;
+      mothMass = o2::constants::physics::MassHyperTriton;
+      chargedDauMass = o2::constants::physics::MassTriton;
+      neutDauMass = o2::constants::physics::MassPi0;
     } else if (hypoMoth == kHyperhelium4sigma) {
       charge = 2;
       mothMass = o2::constants::physics::MassHyperHelium4;
@@ -199,6 +205,8 @@ struct kinkBuilder {
     AxisSpec massAxis(100, 1.1, 1.4, "m (GeV/#it{c}^{2})");
     if (hypoMoth == kSigmaMinus) {
       massAxis = AxisSpec{100, 1.1, 1.4, "m (GeV/#it{c}^{2})"};
+    } else if (hypoMoth == kHypertriton) {
+      massAxis = AxisSpec{100, 2.94, 3.2, "m (GeV/#it{c}^{2})"};
     } else if (hypoMoth == kHyperhelium4sigma) {
       massAxis = AxisSpec{100, 3.85, 4.25, "m (GeV/#it{c}^{2})"};
     }

PWGLF/TableProducer/Nuspex/CMakeLists.txt

@@ -94,8 +94,8 @@ o2physics_add_dpl_workflow(nuclei-flow-trees
     PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore O2::DetectorsBase O2Physics::EventFilteringUtils
     COMPONENT_NAME Analysis)
 
-o2physics_add_dpl_workflow(hyperhelium4sigma-reco-task
-    SOURCES hyperhelium4sigmaRecoTask.cxx
+o2physics_add_dpl_workflow(hyperkink-reco-task
+    SOURCES hyperkinkRecoTask.cxx
     PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore
     COMPONENT_NAME Analysis)