@@ -46,6 +46,9 @@ enum Observable { kstar ///< kstar
4646enum EventType { same, // /< Pair from same event
4747 mixed // /< Pair from mixed event
4848};
49+
50+ static constexpr int nMinMomComponents{4 }; // /< Minimum number of momentum components to compute the femto observable
51+
4952}; // namespace femtoDreamContainer
5053
5154// / \class FemtoDreamContainer
@@ -520,7 +523,7 @@ class FemtoDreamContainer
520523 {
521524
522525 std::vector<double > k3d = FemtoDreamMath::newpairfunc (part1, mMassOne , part2, mMassTwo , IsSameSpecies);
523- if (k3d.size () < 4 ) {
526+ if (k3d.size () < femtoDreamContainer::nMinMomComponents ) {
524527 LOG (error) << " newpairfunc returned size=" << k3d.size ();
525528 return ;
526529 }
@@ -539,7 +542,7 @@ class FemtoDreamContainer
539542 if (part1.has_fdMCParticle () && part2.has_fdMCParticle ()) {
540543
541544 std::vector<double > k3dMC = FemtoDreamMath::newpairfuncMC (part1.fdMCParticle (), mMassOne , part2.fdMCParticle (), mMassTwo , IsSameSpecies);
542- if (k3dMC.size () < 4 ) {
545+ if (k3dMC.size () < femtoDreamContainer::nMinMomComponents ) {
543546 LOG (error) << " newpairfunc returned size=" << k3d.size ();
544547 return ;
545548 }
@@ -565,7 +568,7 @@ class FemtoDreamContainer
565568 {
566569
567570 std::vector<double > k3d = FemtoDreamMath::newpairfunc (part1, mMassOne , part2, mMassTwo , IsSameSpecies);
568- if (k3d.size () < 4 ) {
571+ if (k3d.size () < femtoDreamContainer::nMinMomComponents ) {
569572 LOG (error) << " newpairfunc returned size=" << k3d.size ();
570573 return ;
571574 }
@@ -584,7 +587,7 @@ class FemtoDreamContainer
584587 if (part1.has_fdMCParticle () && part2.has_fdMCParticle ()) {
585588
586589 std::vector<double > k3dMC = FemtoDreamMath::newpairfuncMC (part1.fdMCParticle (), mMassOne , part2.fdMCParticle (), mMassTwo , IsSameSpecies);
587- if (k3dMC.size () < 4 ) {
590+ if (k3dMC.size () < femtoDreamContainer::nMinMomComponents ) {
588591 LOG (error) << " newpairfunc returned size=" << k3d.size ();
589592 return ;
590593 }
0 commit comments