AOD: Add optional streamer to producer

Author: f3sch

Detectors/AOD/include/AODProducerWorkflow/AODProducerWorkflowSpec.h

@@ -30,6 +30,8 @@
 #include "TStopwatch.h"
 #include "ZDCBase/Constants.h"
 #include "GlobalTracking/MatchGlobalFwd.h"
+#include "CommonUtils/TreeStreamRedirector.h"
+#include "CommonUtils/EnumBitOperators.h"
 
 #include <cstdint>
 #include <limits>
@@ -205,7 +207,15 @@ class BunchCrossings
 
   std::vector<TimeWindow> mTimeWindows; // the time window structure covering the complete duration of mBCTimeVector
   double mWindowSize;                   // the size of a single time window
-};                                      // end internal class
+}; // end internal class
+
+// Steering bits for additional output during AOD production
+enum struct AODProducerStreamerMask : uint8_t {
+  None = 0,
+  TrackQA = O2_ENUM_SET_BIT(0),
+  All = std::numeric_limits<std::underlying_type_t<AODProducerStreamerMask>>::max(),
+};
+O2_DEFINE_ENUM_BIT_OPERATORS(AODProducerStreamerMask)
 
 class AODProducerWorkflowDPL : public Task
 {
@@ -243,6 +253,9 @@ class AODProducerWorkflowDPL : public Task
   std::unordered_set<GIndex> mGIDUsedBySVtx;
   std::unordered_set<GIndex> mGIDUsedByStr;
 
+  AODProducerStreamerMask mStreamerMask;
+  std::shared_ptr<o2::utils::TreeStreamRedirector> mStreamer;
+
   int mNThreads = 1;
   bool mUseMC = true;
   bool mEnableSV = true; // enable secondary vertices

Detectors/AOD/src/AODProducerWorkflowSpec.cxx

@@ -51,6 +51,7 @@
 #include "Framework/DataTypes.h"
 #include "Framework/TableBuilder.h"
 #include "Framework/CCDBParamSpec.h"
+#include "CommonUtils/TreeStreamRedirector.h"
 #include "FT0Base/Geometry.h"
 #include "GlobalTracking/MatchTOF.h"
 #include "ReconstructionDataFormats/Cascade.h"
@@ -1687,6 +1688,14 @@ void AODProducerWorkflowDPL::init(InitContext& ic)
   mThinTracks = ic.options().get<bool>("thin-tracks");
   mPropTracks = ic.options().get<bool>("propagate-tracks");
   mPropMuons = ic.options().get<bool>("propagate-muons");
+  if (auto s = ic.options().get<std::string>("with-streamers"); !s.empty()) {
+    mStreamerMask = static_cast<AODProducerStreamerMask>(std::stoul(s, nullptr, 2));
+    if (O2_ENUM_ANY_BIT(mStreamerMask)) {
+      LOGP(info, "Writing streamer data with mask {:0{}b}", static_cast<std::underlying_type_t<AODProducerStreamerMask>>(mStreamerMask), std::numeric_limits<std::underlying_type_t<AODProducerStreamerMask>>::digits);
+    } else {
+      LOGP(warn, "Specified non-default empty streamer mask!");
+    }
+  }
   mTrackQCFraction = ic.options().get<float>("trackqc-fraction");
   mTrackQCNTrCut = ic.options().get<int64_t>("trackqc-NTrCut");
   if (auto seed = ic.options().get<int>("seed"); seed == 0) {
@@ -1772,6 +1781,10 @@ void AODProducerWorkflowDPL::init(InitContext& ic)
   mHeavyIonUpdate = when;
 
   mTimer.Reset();
+
+  if (O2_ENUM_ANY_BIT(mStreamerMask)) {
+    mStreamer = std::make_unique<o2::utils::TreeStreamRedirector>("AO2DStreamer.root", "RECREATE");
+  }
 }
 
 void AODProducerWorkflowDPL::run(ProcessingContext& pc)
@@ -2643,6 +2656,51 @@ AODProducerWorkflowDPL::TrackQA AODProducerWorkflowDPL::processBarrelTrackQA(int
           trackQAHolder.dRefGloSnp = safeInt8Clamp(((itsCopy.getSnp() + tpcCopy.getSnp()) * 0.5f - gloCopy.getSnp()) * scaleGlo(2));
           trackQAHolder.dRefGloTgl = safeInt8Clamp(((itsCopy.getTgl() + tpcCopy.getTgl()) * 0.5f - gloCopy.getTgl()) * scaleGlo(3));
           trackQAHolder.dRefGloQ2Pt = safeInt8Clamp(((itsCopy.getQ2Pt() + tpcCopy.getQ2Pt()) * 0.5f - gloCopy.getQ2Pt()) * scaleGlo(4));
+
+          if (O2_ENUM_TEST_BIT(mStreamerMask, AODProducerStreamerMask::TrackQA)) {
+            (*mStreamer) << "trackQA"
+                         << "trackITSOrig=" << itsOrig
+                         << "trackTPCOrig=" << tpcOrig
+                         << "trackITSTPCOrig=" << trackPar
+                         << "trackITSProp=" << itsCopy
+                         << "trackTPCProp=" << tpcCopy
+                         << "trackITSTPCProp=" << gloCopy
+                         << "refRadius=" << o2::aod::track::trackQARefRadius
+                         << "scaleBins=" << o2::aod::track::trackQAScaleBins
+                         << "scaleCont0=" << scaleCont(0)
+                         << "scaleCont1=" << scaleCont(1)
+                         << "scaleCont2=" << scaleCont(2)
+                         << "scaleCont3=" << scaleCont(3)
+                         << "scaleCont4=" << scaleCont(4)
+                         << "scaleGlo0=" << scaleGlo(0)
+                         << "scaleGlo1=" << scaleGlo(1)
+                         << "scaleGlo2=" << scaleGlo(2)
+                         << "scaleGlo3=" << scaleGlo(3)
+                         << "scaleGlo4=" << scaleGlo(4)
+                         << "trackQAHolder.tpcTime0=" << trackQAHolder.tpcTime0
+                         << "trackQAHolder.tpcdcaR=" << trackQAHolder.tpcdcaR
+                         << "trackQAHolder.tpcdcaZ=" << trackQAHolder.tpcdcaZ
+                         << "trackQAHolder.tpcdcaClusterByteMask=" << trackQAHolder.tpcClusterByteMask
+                         << "trackQAHolder.tpcdEdxMax0R=" << trackQAHolder.tpcdEdxMax0R
+                         << "trackQAHolder.tpcdEdxMax1R=" << trackQAHolder.tpcdEdxMax1R
+                         << "trackQAHolder.tpcdEdxMax2R=" << trackQAHolder.tpcdEdxMax2R
+                         << "trackQAHolder.tpcdEdxMax3R=" << trackQAHolder.tpcdEdxMax3R
+                         << "trackQAHolder.tpcdEdxTot0R=" << trackQAHolder.tpcdEdxTot0R
+                         << "trackQAHolder.tpcdEdxTot1R=" << trackQAHolder.tpcdEdxTot1R
+                         << "trackQAHolder.tpcdEdxTot2R=" << trackQAHolder.tpcdEdxTot2R
+                         << "trackQAHolder.tpcdEdxTot3R=" << trackQAHolder.tpcdEdxTot3R
+                         << "trackQAHolder.dRefContY=" << trackQAHolder.dRefContY
+                         << "trackQAHolder.dRefContZ=" << trackQAHolder.dRefContZ
+                         << "trackQAHolder.dRefContSnp=" << trackQAHolder.dRefContSnp
+                         << "trackQAHolder.dRefContTgl=" << trackQAHolder.dRefContTgl
+                         << "trackQAHolder.dRefContQ2Pt=" << trackQAHolder.dRefContQ2Pt
+                         << "trackQAHolder.dRefGloY=" << trackQAHolder.dRefGloY
+                         << "trackQAHolder.dRefGloZ=" << trackQAHolder.dRefGloZ
+                         << "trackQAHolder.dRefGloSnp=" << trackQAHolder.dRefGloSnp
+                         << "trackQAHolder.dRefGloTgl=" << trackQAHolder.dRefGloTgl
+                         << "trackQAHolder.dRefGloQ2Pt=" << trackQAHolder.dRefGloQ2Pt
+                         << "\n";
+          }
         }
       }
     }
@@ -3012,6 +3070,8 @@ void AODProducerWorkflowDPL::endOfStream(EndOfStreamContext& /*ec*/)
 {
   LOGF(info, "aod producer dpl total timing: Cpu: %.3e Real: %.3e s in %d slots",
        mTimer.CpuTime(), mTimer.RealTime(), mTimer.Counter() - 1);
+
+  mStreamer.reset();
 }
 
 DataProcessorSpec getAODProducerWorkflowSpec(GID::mask_t src, bool enableSV, bool enableStrangenessTracking, bool useMC, bool CTPConfigPerRun)
@@ -3144,6 +3204,7 @@ DataProcessorSpec getAODProducerWorkflowSpec(GID::mask_t src, bool enableSV, boo
       ConfigParamSpec{"thin-tracks", VariantType::Bool, false, {"Produce thinned track tables"}},
       ConfigParamSpec{"trackqc-fraction", VariantType::Float, float(0.1), {"Fraction of tracks to QC"}},
       ConfigParamSpec{"trackqc-NTrCut", VariantType::Int64, 4L, {"Minimal length of the track - in amount of tracklets"}},
+      ConfigParamSpec{"with-streamers", VariantType::String, "", {"Bit-mask to steer writing of intermediate streamer files"}},
       ConfigParamSpec{"seed", VariantType::Int, 0, {"Set seed for random generator used for sampling (0 (default) means using a random_device)"}},
     }};
 }