Additional QA

Author: ddobrigk

Common/DataModel/Centrality.h

@@ -43,8 +43,8 @@ DECLARE_SOA_TABLE(CentRun2SPDTrks, "AOD", "CENTRUN2SPDTRK", cent::CentRun2SPDTra
 DECLARE_SOA_TABLE(CentRun2SPDClss, "AOD", "CENTRUN2SPDCLS", cent::CentRun2SPDClusters);  //! Run 2 SPD clusters centrality table
 DECLARE_SOA_TABLE(CentRun2CL0s, "AOD", "CENTRUN2CL0", cent::CentRun2CL0);                //! Run 2 CL0 centrality table
 DECLARE_SOA_TABLE(CentRun2CL1s, "AOD", "CENTRUN2CL1", cent::CentRun2CL1);                //! Run 2 CL1 centrality table
-DECLARE_SOA_TABLE(CentRun2RefMult5, "AOD", "CENTRUN2REFMULT5", cent::CentRun2RefMult5);  //! Run 2, ref mult |eta| < 0.5
-DECLARE_SOA_TABLE(CentRun2RefMult8, "AOD", "CENTRUN2REFMULT8", cent::CentRun2RefMult8);  //! Run 2, ref mult |eta| < 0.8
+DECLARE_SOA_TABLE(CentRun2RefMult5s, "AOD", "CENTRUN2REFMULT5", cent::CentRun2RefMult5); //! Run 2, ref mult |eta| < 0.5
+DECLARE_SOA_TABLE(CentRun2RefMult8s, "AOD", "CENTRUN2REFMULT8", cent::CentRun2RefMult8); //! Run 2, ref mult |eta| < 0.8
 
 // Run 3 tables
 DECLARE_SOA_TABLE(CentFV0As, "AOD", "CENTFV0A", cent::CentFV0A);                         //! Run 3 FV0A centrality table
@@ -63,6 +63,8 @@ using CentRun2SPDTrk = CentRun2SPDTrks::iterator;
 using CentRun2SPDCls = CentRun2SPDClss::iterator;
 using CentRun2CL0 = CentRun2CL0s::iterator;
 using CentRun2CL1 = CentRun2CL1s::iterator;
+using CentRun2RefMult5 = CentRun2RefMult5s::iterator;
+using CentRun2RefMult8 = CentRun2RefMult8s::iterator;
 using CentFV0A = CentFV0As::iterator;
 using CentFT0M = CentFT0Ms::iterator;
 using CentFT0A = CentFT0As::iterator;

Common/LegacyDataQA/CMakeLists.txt

@@ -12,4 +12,9 @@
 o2physics_add_dpl_workflow(otfv0qa
                     SOURCES otfv0qa.cxx
                     PUBLIC_LINK_LIBRARIES O2::Framework O2Physics::AnalysisCore
+                    COMPONENT_NAME Analysis)
+
+o2physics_add_dpl_workflow(centqa
+                    SOURCES centqa.cxx
+                    PUBLIC_LINK_LIBRARIES O2::Framework O2Physics::AnalysisCore
                     COMPONENT_NAME Analysis)

Common/LegacyDataQA/centqa.cxx

@@ -0,0 +1,92 @@
+// Copyright 2019-2020 CERN and copyright holders of ALICE O2.
+// See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.
+// All rights not expressly granted are reserved.
+//
+// This software is distributed under the terms of the GNU General Public
+// License v3 (GPL Version 3), copied verbatim in the file "COPYING".
+//
+// In applying this license CERN does not waive the privileges and immunities
+// granted to it by virtue of its status as an Intergovernmental Organization
+// or submit itself to any jurisdiction.
+//
+// This code calculates output histograms for centrality calibration
+// as well as vertex-Z dependencies of raw variables (either for calibration
+// of vtx-Z dependencies or for the calibration of those).
+//
+// This task is not strictly necessary in a typical analysis workflow,
+// except for centrality calibration! The necessary task is the multiplicity
+// tables.
+//
+// Comments, suggestions, questions? Please write to:
+// - [email protected]
+// - [email protected]
+//
+
+#include "Framework/runDataProcessing.h"
+#include "Framework/AnalysisTask.h"
+#include "Framework/AnalysisDataModel.h"
+#include "Common/DataModel/McCollisionExtra.h"
+#include "Common/DataModel/Centrality.h"
+#include "Common/DataModel/EventSelection.h"
+#include "Framework/O2DatabasePDGPlugin.h"
+#include "TH1F.h"
+#include "TH2F.h"
+
+using namespace o2;
+using namespace o2::framework;
+
+struct CentQA {
+  // Raw multiplicities
+  HistogramRegistry histos{"Histos", {}, OutputObjHandlingPolicy::AnalysisObject};
+
+  void init(InitContext&)
+  {
+    const AxisSpec axisCentrality{205, 0, 205, "centrality"};
+
+    // Base histograms
+    histos.add("hV0M", "V0M centrality", kTH1F, {axisCentrality});
+    histos.add("hV0A", "V0A centrality", kTH1F, {axisCentrality});
+    histos.add("hCL0", "CL0 centrality", kTH1F, {axisCentrality});
+    histos.add("hCL1", "CL1 centrality", kTH1F, {axisCentrality});
+    histos.add("hRefMult5", "RefMult5 centrality", kTH1F, {axisCentrality});
+    histos.add("hRefMult8", "RefMult8 centrality", kTH1F, {axisCentrality});
+  }
+
+  void processV0M(soa::Join<aod::Collision, aod::CentRun2V0Ms>::iterator const& col)
+  {
+    histos.fill(HIST("hV0M"), col.centRun2V0M());
+  }
+  void processV0A(soa::Join<aod::Collision, aod::CentRun2V0As>::iterator const& col)
+  {
+    histos.fill(HIST("hV0A"), col.centRun2V0A());
+  }
+  void processCL0(soa::Join<aod::Collision, aod::CentRun2CL0s>::iterator const& col)
+  {
+    histos.fill(HIST("hCL0"), col.centRun2CL0());
+  }
+  void processCL1(soa::Join<aod::Collision, aod::CentRun2CL1s>::iterator const& col)
+  {
+    histos.fill(HIST("hCL1"), col.centRun2CL1());
+  }
+  void processRefMult5(soa::Join<aod::Collision, aod::CentRun2RefMult5s>::iterator const& col)
+  {
+    histos.fill(HIST("hRefMult5"), col.centRun2RefMult5());
+  }
+  void processRefMult8(soa::Join<aod::Collision, aod::CentRun2RefMult8s>::iterator const& col)
+  {
+    histos.fill(HIST("hRefMult8"), col.centRun2RefMult8());
+  }
+
+  PROCESS_SWITCH(CentQA, processV0M, "QA V0M centrality", true);
+  PROCESS_SWITCH(CentQA, processV0A, "QA V0A centrality", false);
+  PROCESS_SWITCH(CentQA, processCL0, "QA CL0 centrality", false);
+  PROCESS_SWITCH(CentQA, processCL1, "QA CL1 centrality", false);
+  PROCESS_SWITCH(CentQA, processRefMult5, "QA RefMult5 centrality", false);
+  PROCESS_SWITCH(CentQA, processRefMult8, "QA RefMult8 centrality", false);
+};
+
+WorkflowSpec defineDataProcessing(ConfigContext const& cfgc)
+{
+  return WorkflowSpec{
+    adaptAnalysisTask<CentQA>(cfgc)};
+}

Common/LegacyDataQA/otfv0qa.cxx

@@ -42,14 +42,17 @@ struct OTFV0Qa {
   void init(InitContext&)
   {
     const AxisSpec axisEvent{10, 0, 10, "Event counter"};
+    const AxisSpec axisNCandidates{500, 0, 500, "Number of OTF v0s"};
 
     // Base histograms
-    histos.add("hEventCounter", "Event counter", kTH1D, {axisEvent});
+    histos.add("hEventCounter", "Event counter", kTH1F, {axisEvent});
+    histos.add("hCandidates", "Number of OTF V0s", kTH1F, {axisNCandidates});
   }
 
-  void process(aod::Collision const& col)
+  void process(aod::Collision const& col, aod::OTFV0s const& v0s)
   {
     histos.fill(HIST("hEventCounter"), 0.5);
+    histos.fill(HIST("hCandidates"), v0s.size());
   }
 };