First full version of derived data consumer. Implemented preprocessor macro to avoid rewriting 20*4 histogram fills by hand (now fixed at compile time). Still not works with current JetPolarizationIons definition (collIdx has trouble with Preslice, as it is a non-monotonic. Possibly O2 is mixing the TimeFrames and combining multiple GlobalIndices in a same subfolder of the derived data)

Author: cmuncinelli

PWGLF/DataModel/lambdaJetPolarizationIons.h

@@ -63,7 +63,7 @@ DECLARE_SOA_TABLE(RingJets, "AOD", "RINGJETS", // Renamed to follow convention o
                   lambdajetpol::JetPhi,
                   lambdajetpol::JetNConstituents);
 
-DECLARE_SOA_TABLE(RingLaV0s, "AOD", "RINGLAV0S",
+DECLARE_SOA_TABLE(RingLaV0s, "AOD", "RINGLAV0S", // Had to write this in a shorter form because the derived data did not accept long names
                   lambdajetpol::CollIdx,
                   lambdajetpol::V0Pt,
                   lambdajetpol::V0Eta,

PWGLF/Tasks/Strangeness/CMakeLists.txt

@@ -171,6 +171,11 @@ SOURCES lambdaJetPolarizationIons.cxx
 PUBLIC_LINK_LIBRARIES O2::Framework O2Physics::AnalysisCore O2Physics::PWGJECore FastJet::FastJet FastJet::Contrib O2Physics::EventFilteringUtils O2Physics::AnalysisCCDB
 COMPONENT_NAME Analysis)
 
+o2physics_add_dpl_workflow(lambdajetpolarizationionsderived
+SOURCES lambdaJetPolarizationIonsDerived.cxx
+PUBLIC_LINK_LIBRARIES O2::Framework O2Physics::AnalysisCore
+COMPONENT_NAME Analysis)
+
 o2physics_add_dpl_workflow(lambdaspincorrderived
     SOURCES lambdaspincorrderived.cxx
     PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore

PWGLF/Tasks/Strangeness/lambdaJetPolarizationIons.cxx

@@ -1536,7 +1536,7 @@ struct lambdajetpolarizationions {
         const uint64_t collIdx = collision.globalIndex();
         if (v0Selections.rejectTPCsectorBoundary) initCCDB(bc); // Substituted call from collision to bc for raw data
 
-        bool hasValidV0 = false; // Bool to know if event information can be saved.
+        // bool hasValidV0 = false; // Bool to know if event information can be saved.
         for (auto const& v0 : fullV0s){
             V0SelCounter.resetForNewV0();
             V0SelCounter.fill(); // Fill for all v0 candidates
@@ -1552,7 +1552,7 @@ struct lambdajetpolarizationions {
             if (analyseAntiLambda) isAntiLambda = passesLambdaLambdaBarHypothesis(v0, collision, false);
 
             if (!isLambda && !isAntiLambda) continue; // Candidate is not considered to be a Lambda (TODO: expand this to a full if block with QA about rejections)
-            hasValidV0 = true;
+            // hasValidV0 = true;
 
             if (doArmenterosQA) histos.fill(HIST("GeneralQA/h2dArmenterosSelected"), v0.alpha(), v0.qtarm()); // cross-check
             if (isLambda && isAntiLambda) histos.fill(HIST("hAmbiguousLambdaCandidates"), 0);
@@ -1739,10 +1739,11 @@ struct lambdajetpolarizationions {
             } // end CompleteTopoQA
         } // end V0s loop
         // Only fills collision when there is a valid V0 in it: (TODO: could probably do the same for the jets table)
-        if (hasValidV0){
+        // if (hasValidV0){
             // LOG(INFO) << "Filling tableCollisions";
-            tableCollisions(collIdx, centrality); // (TODO: add InteractionRate info and other useful cuts for later on in the analysis!)
-        }
+        // Current logic now fills tables independently of collision having V0s, for the Jets table to match correctly
+        tableCollisions(collIdx, centrality); // (TODO: add InteractionRate info and other useful cuts for later on in the analysis!)
+        // }
     }
 
     PROCESS_SWITCH(lambdajetpolarizationions, processJetsData, "Process jets and produce derived data in Run 3 Data", true);