K*(892) optimized (#6978)

* adding new task for k*(892) using filters on resonancetables

* modified cmakelists

* adding partitions for tpc PID at low pt, process functions

* adding Occupancy and ImpactParameter to ResonanceTables

* adding Maximum occupancy cut to PWGLF/Utils/collisionCuts.h

Author: murioni

PWGLF/DataModel/LFResonanceTables.h

@@ -66,6 +66,7 @@ DECLARE_SOA_COLUMN(IsINELgt0, isINELgt0, bool);               //! INEL>0
 DECLARE_SOA_COLUMN(IsTriggerTVX, isTriggerTVX, bool);         //! TriggerTVX
 DECLARE_SOA_COLUMN(IsInSel8, isInSel8, bool);                 //! InSel8
 DECLARE_SOA_COLUMN(IsInAfterAllCuts, isInAfterAllCuts, bool); //! InAfterAllCuts
+DECLARE_SOA_COLUMN(ImpactParameter, impactParameter, float);  //! ImpactParameter
 
 } // namespace resocollision
 DECLARE_SOA_TABLE(ResoCollisions, "AOD", "RESOCOLLISION",
@@ -89,7 +90,8 @@ DECLARE_SOA_TABLE(ResoMCCollisions, "AOD", "RESOMCCOL",
                   resocollision::IsINELgt0,
                   resocollision::IsTriggerTVX,
                   resocollision::IsInSel8,
-                  resocollision::IsInAfterAllCuts);
+                  resocollision::IsInAfterAllCuts,
+                  resocollision::ImpactParameter);
 using ResoMCCollision = ResoMCCollisions::iterator;
 
 // Resonance Daughters

PWGLF/TableProducer/Resonances/LFResonanceInitializer.cxx

@@ -89,6 +89,7 @@ struct reso2initializer {
   /// Event cuts
   o2::analysis::CollisonCuts colCuts;
   Configurable<float> ConfEvtZvtx{"ConfEvtZvtx", 10.f, "Evt sel: Max. z-Vertex (cm)"};
+  Configurable<int> ConfEvtOccupancyInTimeRange{"ConfEvtOccupancyInTimeRange", -1, "Evt sel: maximum track occupancy"};
   Configurable<bool> ConfEvtTriggerCheck{"ConfEvtTriggerCheck", false, "Evt sel: check for trigger"};
   Configurable<int> ConfEvtTriggerSel{"ConfEvtTriggerSel", 8, "Evt sel: trigger"};
   Configurable<bool> ConfEvtOfflineCheck{"ConfEvtOfflineCheck", true, "Evt sel: check for offline selection"};
@@ -837,7 +838,7 @@ struct reso2initializer {
   }
 
   template <bool isRun2, typename MCCol, typename MCPart>
-  void fillMCCollision(MCCol const& mccol, MCPart const& mcparts)
+  void fillMCCollision(MCCol const& mccol, MCPart const& mcparts, float impactpar = -999.0)
   {
     auto centrality = 0.0;
     if constexpr (!isRun2)
@@ -849,7 +850,7 @@ struct reso2initializer {
     bool isTriggerTVX = mccol.selection_bit(aod::evsel::kIsTriggerTVX);
     bool isSel8 = mccol.sel8();
     bool isSelected = colCuts.isSelected(mccol);
-    resoMCCollisions(inVtx10, isTrueINELgt0, isTriggerTVX, isSel8, isSelected);
+    resoMCCollisions(inVtx10, isTrueINELgt0, isTriggerTVX, isSel8, isSelected, impactpar);
 
     // QA for Trigger efficiency
     qaRegistry.fill(HIST("Event/hMCEventIndices"), centrality, aod::resocollision::kINEL);
@@ -923,7 +924,7 @@ struct reso2initializer {
     if (doprocessTrackDataRun2 || doprocessTrackV0DataRun2 || doprocessTrackV0CascDataRun2 || doprocessTrackMCRun2 || doprocessTrackV0MCRun2 || doprocessTrackV0CascMCRun2) {
       colCuts.setCuts(ConfEvtZvtx, ConfEvtTriggerCheck, ConfEvtTriggerSel, ConfEvtOfflineCheck, false);
     } else if (doprocessTrackData || doprocessTrackV0Data || doprocessTrackV0CascData || doprocessTrackMC || doprocessTrackV0MC || doprocessTrackV0CascMC || doprocessTrackEPData) {
-      colCuts.setCuts(ConfEvtZvtx, ConfEvtTriggerCheck, ConfEvtTriggerSel, ConfEvtOfflineCheck, true);
+      colCuts.setCuts(ConfEvtZvtx, ConfEvtTriggerCheck, ConfEvtTriggerSel, ConfEvtOfflineCheck, true, false, ConfEvtOccupancyInTimeRange);
     }
     colCuts.init(&qaRegistry);
     colCuts.setTriggerTVX(ConfEvtTriggerTVXSel);
@@ -1144,7 +1145,10 @@ struct reso2initializer {
     colCuts.fillQA(collision);
 
     resoCollisions(0, collision.posX(), collision.posY(), collision.posZ(), CentEst(collision), ComputeSpherocity(tracks, trackSphMin, trackSphDef), 0., 0., 0., 0., d_bz, bc.timestamp());
-    fillMCCollision<false>(collision, mcParticles);
+
+    auto mccollision = collision.mcCollision_as<aod::McCollisions>();
+    float impactpar = mccollision.impactParameter();
+    fillMCCollision<false>(collision, mcParticles, impactpar);
 
     // Loop over tracks
     fillTracks<true>(collision, tracks);

PWGLF/Tasks/Resonances/CMakeLists.txt

@@ -24,6 +24,11 @@ o2physics_add_dpl_workflow(k892analysis
     PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore
     COMPONENT_NAME Analysis)
 
+o2physics_add_dpl_workflow(kstar892analysis
+    SOURCES kstar892analysis.cxx
+    PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore
+    COMPONENT_NAME Analysis)
+
 o2physics_add_dpl_workflow(k892pmanalysis
     SOURCES k892pmanalysis.cxx
     PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore