solution A

Author: Andreas Mikalsen

PWGDQ/Tasks/tableReader_withAssoc.cxx

@@ -280,6 +280,7 @@ struct AnalysisEventSelection {
 
   Produces<aod::EventCuts> eventSel;
   Produces<aod::MixingHashes> hash;
+  Produces<aod::StoredReducedEvents> JetEvents; 
   OutputObj<THashList> fOutputList{"output"};
 
   // TODO: Provide the mixing variables and binning directly via configurables (e.g. vectors of float)
@@ -290,6 +291,7 @@ struct AnalysisEventSelection {
   Configurable<std::string> fConfigAddEventHistogram{"cfgAddEventHistogram", "", "Comma separated list of histograms"};
   Configurable<std::string> fConfigAddJSONHistograms{"cfgAddJSONHistograms", "", "Add event histograms defined via JSON formatting (see HistogramsLibrary)"};
   Configurable<bool> fConfigQA{"cfgQA", true, "If true, QA histograms will be created and filled"};
+  Configurable<bool> fillJetEvents{"fillJetEvents", false, "Fill JetEvents table"};
 
   Configurable<int> fConfigITSROFrameStartBorderMargin{"cfgITSROFrameStartBorderMargin", -1, "Number of bcs at the start of ITS RO Frame border. Take from CCDB if -1"};
   Configurable<int> fConfigITSROFrameEndBorderMargin{"cfgITSROFrameEndBorderMargin", -1, "Number of bcs at the end of ITS RO Frame border. Take from CCDB if -1"};
@@ -547,7 +549,21 @@ struct AnalysisEventSelection {
   {
     runEventSelection<gkEventFillMap>(events);
     publishSelections<gkEventFillMap>(events);
+
+  if (!fillJetEvents) {
+    return;
+  }
+  for (auto& event : events) {
+    JetEvents(event.tag_raw(),
+              event.runNumber(),
+              event.posX(),
+              event.posY(),
+              event.posZ(),
+              event.numContrib(),
+              event.collisionTime(),
+              event.collisionTimeRes());
   }
+}
   void processSkimmedWithZdc(MyEventsZdc const& events)
   {
     runEventSelection<gkEventFillMapWithZdc>(events);