@@ -45,7 +45,7 @@ struct sigmaminustask {
4545 // Configurable for event selection
4646 Configurable<float > cutzvertex{" cutzvertex" , 10 .0f , " Accepted z-vertex range (cm)" };
4747 Configurable<float > cutNSigmaPi{" cutNSigmaPi" , 4 , " NSigmaTPCPion" };
48- Configurable<float > cutEtaMotherMC{ " cutEtaMotherMC " , 1 .0f , " Eta cut for mother Sigma in MC" };
48+ Configurable<float > cutRapMotherMC{ " cutRapMotherMC " , 1 .0f , " Rapidity cut for mother Sigma in MC" };
4949
5050 Configurable<bool > fillOutputTree{" fillOutputTree" , true , " If true, fill the output tree with Kink candidates" };
5151
@@ -153,7 +153,7 @@ struct sigmaminustask {
153153 if (piMother.globalIndex () != mcTrackSigma.globalIndex ()) {
154154 continue ;
155155 }
156- if (std::abs (mcTrackSigma.pdgCode ()) != 3112 ) {
156+ if (std::abs (mcTrackSigma.pdgCode ()) != 3112 && std::abs (mcTrackSigma. pdgCode ()) != 3222 ) {
157157 continue ;
158158 }
159159 if (std::abs (mcTrackPiDau.pdgCode ()) != 211 && std::abs (mcTrackPiDau.pdgCode ()) != 2212 ) {
@@ -200,7 +200,7 @@ struct sigmaminustask {
200200
201201 // Loop over all generated particles to fill MC histograms
202202 for (const auto & mcPart : particlesMC) {
203- if (std::abs (mcPart.pdgCode ()) != 3112 || std::abs (mcPart.y ()) > cutEtaMotherMC ) { // only sigma mothers and rapidity cut
203+ if (( std::abs (mcPart.pdgCode ()) != 3112 && std::abs (mcPart. pdgCode ()) != 3222 ) || std::abs (mcPart.y ()) > cutRapMotherMC ) { // only sigma mothers and rapidity cut
204204 continue ;
205205 }
206206 if (!mcPart.has_daughters ()) {
0 commit comments