PWGCF:FemtoUniverse - D0/D0bar analysis update (#5541)

* PWGCF:FemtoUniverse - D0/D0bar analysis update

* PWGCF:FemtoUniverse - fixing clang

Author: kgwizdzi

PWGCF/FemtoUniverse/TableProducer/femtoUniverseProducerTask.cxx

@@ -246,6 +246,9 @@ struct femtoUniverseProducerTask {
   struct : o2::framework::ConfigurableGroup {
     Configurable<float> ConfD0D0barCandMaxY{"ConfD0D0barCandMaxY", -1., "max. cand. rapidity"};
     Configurable<float> ConfD0D0barCandEtaCut{"ConfD0D0barCandEtaCut", 0.8, "max. cand. pseudorapidity"};
+    Configurable<bool> ConfStoreD0D0barWithinTheMassRange{"ConfStoreD0D0barWithinTheMassRange", false, "Switch to save D0/D0bar within declared inv. mass range"};
+    Configurable<float> ConfStoreD0D0barInvMassLowLimit{"ConfStoreD0D0barInvMassLowLimit", 1.810, "Lower inv. mass limit of D0/D0bar candidate"};
+    Configurable<float> ConfStoreD0D0barInvMassUpLimit{"ConfStoreD0D0barInvMassUpLimit", 1.922, "Upper inv. mass limit of D0/D0bar candidate"};
   } ConfD0Selection;
 
   HfHelper hfHelper;
@@ -784,7 +787,7 @@ struct femtoUniverseProducerTask {
         continue;
       }
 
-      if (ConfD0Selection.ConfD0D0barCandEtaCut >= hfCand.eta()) {
+      if (std::abs(hfCand.eta()) < ConfD0Selection.ConfD0D0barCandEtaCut) {
         continue;
       }
 
@@ -815,6 +818,11 @@ struct femtoUniverseProducerTask {
         isD0D0bar = false;
       }
 
+      if (ConfD0Selection.ConfStoreD0D0barWithinTheMassRange) {
+        if ((invMassD0 > ConfD0Selection.ConfStoreD0D0barInvMassLowLimit && invMassD0 < ConfD0Selection.ConfStoreD0D0barInvMassUpLimit) || (invMassD0bar > ConfD0Selection.ConfStoreD0D0barInvMassLowLimit && invMassD0bar < ConfD0Selection.ConfStoreD0D0barInvMassUpLimit))
+          continue;
+      }
+
       if (isD0D0bar) {
         outputParts(outputCollision.lastIndex(),
                     hfCand.ptProng0(),