2020#ifndef PWGLF_UTILS_MCPARTICLE_H_
2121#define PWGLF_UTILS_MCPARTICLE_H_
2222
23+ #include " ReconstructionDataFormats/PID.h"
24+
2325#include < string>
2426#include < vector>
2527
26- #include " ReconstructionDataFormats/PID.h"
27-
2828namespace o2
2929{
3030namespace pwglf
@@ -55,6 +55,7 @@ class PIDExtended
5555 static constexpr ID OmegaMinus = 16 ;
5656 static constexpr ID HyperHelium4 = 17 ;
5757 static constexpr ID HyperHelium5 = 18 ;
58+ static constexpr ID Lambda1520 = 19 ;
5859
5960 static_assert (Electron == o2::track::PID::Electron, " PID::Electron mismatch" );
6061 static_assert (Muon == o2::track::PID::Muon, " PID::Muon mismatch" );
@@ -79,11 +80,11 @@ class PIDExtended
7980 static constexpr ID PIDCountsUntilAl = 9 ; // Number of indices defined in PID.h equivalent to o2::track::PID::NIDs
8081 // static_assert(PIDCountsUntilAl == o2::track::PID::NIDs, "PID::NIDs mismatch");
8182
82- static constexpr ID PIDCounts = 19 ; // Number of indices defined in PID.h
83+ static constexpr ID PIDCounts = 20 ; // Number of indices defined in PID.h
8384 // static_assert(PIDCounts == o2::track::PID::NIDsTot, "PID::NIDsTot mismatch");
8485 // Define an array of IDs
8586 static constexpr std::array<ID, PIDCountsUntilAl> mIDsUntilAl = {Electron, Muon, Pion, Kaon, Proton, Deuteron, Triton, Helium3, Alpha};
86- static constexpr std::array<ID, PIDCounts> mIDs = {Electron, Muon, Pion, Kaon, Proton, Deuteron, Triton, Helium3, Alpha, PI0, Photon, K0, Lambda, HyperTriton, Hyperhydrog4, XiMinus, OmegaMinus, HyperHelium4, HyperHelium5};
87+ static constexpr std::array<ID, PIDCounts> mIDs = {Electron, Muon, Pion, Kaon, Proton, Deuteron, Triton, Helium3, Alpha, PI0, Photon, K0, Lambda, HyperTriton, Hyperhydrog4, XiMinus, OmegaMinus, HyperHelium4, HyperHelium5, Lambda1520 };
8788
8889 // Define the antiparticles
8990 static constexpr ID Positron = PIDCounts;
@@ -126,7 +127,8 @@ class PIDExtended
126127 static constexpr ID XiCCPlusPlus = PIDCounts + 36 ;
127128 static constexpr ID XiCPlus = PIDCounts + 37 ;
128129 static constexpr ID XiC0 = PIDCounts + 38 ;
129- static constexpr ID NIDsTot = PIDCounts + 39 ;
130+ static constexpr ID AntiLambda1520 = PIDCounts + 39 ;
131+ static constexpr ID NIDsTot = PIDCounts + 40 ;
130132
131133 static constexpr const char * sNames [NIDsTot + 1 ] = {
132134 o2::track::pid_constants::sNames [Electron], // Electron
@@ -148,6 +150,7 @@ class PIDExtended
148150 o2::track::pid_constants::sNames [OmegaMinus], // OmegaMinus
149151 " HyperHelium4" , // HyperHelium4
150152 " HyperHelium5" , // HyperHelium5
153+ " Lambda1520" , // Lambda1520
151154 " Positron" , // Positron
152155 " MuonPlus" , // MuonPlus
153156 " PionMinus" , // PionMinus
@@ -187,6 +190,7 @@ class PIDExtended
187190 " XiCCPlusPlus" , // XiCCPlusPlus
188191 " XiCPlus" , // XiCPlus
189192 " XiC0" , // XiC0
193+ " AntiLambda1520" , // AntiLambda1520
190194 nullptr };
191195
192196 static std::vector<std::string> arrayNames ()
@@ -321,6 +325,10 @@ class PIDExtended
321325 return XiCPlus;
322326 case o2::constants::physics::Pdg::kXiC0 :
323327 return XiC0;
328+ case 3124 :
329+ return Lambda1520;
330+ case -3124 :
331+ return AntiLambda1520;
324332 default :
325333 LOG (debug) << " Cannot identify particle with PDG code " << particle.pdgCode ();
326334 break ;
0 commit comments