Changes to implement phi particle in FemtoDream

modified/ added new files for particle reconstruction via the inv mass
due to short life time. Selection is implemented in the ResoSelection
header. For this new types in femtoDreived.h were implemented.
The Producer task was changed to reconstruct the phi and filltables.
In order to get pid-cuts for further partitioning(inside debug-reso and
pair-task) the cutculator was
modified.

Author: Christopher Klumm

PWGCF/FemtoDream/Core/femtoDreamObjectSelection.h

@@ -60,12 +60,12 @@ class FemtoDreamObjectSelection
     }
   }
 
-  /// assigns value from configurbale to protected class member
-  ///  \param selVals o2 configurable containing the values employed for the selection
-  template <typename V>
+  ///assigns value from configurbale to protected class member 
+  /// \param selVals o2 configurable containing the values employed for the selection
+  template<typename V>
   void assign(V& selVals)
   {
-    assignedValue = static_cast<selValDataType>(selVals);
+    assignedValue = static_cast<selValDataType>(selVals);   
   }
 
   /// Pass the Configurable of selection values in the analysis task to the selection class
@@ -201,7 +201,7 @@ class FemtoDreamObjectSelection
   HistogramRegistry* mHistogramRegistry;                                     ///< For Analysis QA output
   HistogramRegistry* mQAHistogramRegistry;                                   ///< For QA output
   std::vector<FemtoDreamSelection<selValDataType, selVariable>> mSelections; ///< Vector containing all selections
-  selValDataType assignedValue;
+  selValDataType assignedValue; 
 };
 
 } // namespace femtoDream

PWGCF/FemtoDream/Core/femtoDreamTrackSelection.h

@@ -33,6 +33,9 @@
 #include <string>
 #include <vector>
 
+
+
+
 using namespace o2::framework;
 using namespace o2::analysis::femtoDream::femtoDreamSelection;
 
@@ -54,7 +57,7 @@ enum TrackSel { kSign,         ///< Sign of the track
                 kDCAxyMax,     ///< Max. DCA_xy (cm)
                 kDCAzMax,      ///< Max. DCA_z (cm)
                 kDCAMin,       ///< Min. DCA_xyz (cm)
-                kPIDnSigmaMax  ///< Max. |n_sigma| for PID
+                kPIDnSigmaMax ///< Max. |n_sigma| for PID
 };
 
 enum TrackContainerPosition {
@@ -425,7 +428,7 @@ auto FemtoDreamTrackSelection::getNsigmaITS(T const& track, o2::track::PID pid)
 }
 
 template <typename T>
-bool FemtoDreamTrackSelection::isSelectedMinimal(T const& track, bool UseThreshold)
+bool FemtoDreamTrackSelection::isSelectedMinimal(T const& track, bool UseThreshold ) 
 {
   const auto pT = track.pt();
   const auto eta = track.eta();
@@ -502,19 +505,18 @@ bool FemtoDreamTrackSelection::isSelectedMinimal(T const& track, bool UseThresho
     bool pass = false;
     for (size_t i = 0; i < pidTPC.size(); ++i) {
 
-      auto pidTPCVal = pidTPC.at(i);
-      if (pT < nPTPCThr) {
-        pass = std::fabs(pidTPCVal) < nSigmaPIDMax;
-      } else if (pT >= nPTPCThr) {
-        auto pidTOFVal = pidTOF.at(i);
-        pass = std::sqrt(pidTPCVal * pidTPCVal + pidTOFVal * pidTOFVal) < nSigmaPIDMax;
-      }
-      if (pass)
-        break; // early exit if any condition is satisfied
+        auto pidTPCVal = pidTPC.at(i);
+        if(pT < nPTPCThr) {
+            pass = std::fabs(pidTPCVal) < nSigmaPIDMax;
+        }else if (pT >= nPTPCThr) {
+            auto pidTOFVal = pidTOF.at(i);
+            pass = std::sqrt(pidTPCVal * pidTPCVal + pidTOFVal * pidTOFVal) < nSigmaPIDMax;
+        }
+        if (pass) break; // early exit if any condition is satisfied
     }
 
     if (!pass) {
-      return false;
+        return false;
     }
   }
   return true;