FT0 mc->raw->digits test added to full_system_test_ci_extra_tests

Author: shahor02

Detectors/FIT/FT0/reconstruction/CMakeLists.txt

@@ -30,3 +30,9 @@ o2_target_root_dictionary(
           include/FT0Reconstruction/CTFCoder.h
           include/FT0Reconstruction/InteractionTag.h
   )
+
+o2_add_executable(
+  test-raw-conversion
+  COMPONENT_NAME ft0
+  SOURCES src/test-raw-conversion.cxx
+  PUBLIC_LINK_LIBRARIES O2::FT0Reconstruction)

Detectors/FIT/FT0/reconstruction/src/test-raw-conversion.cxx

@@ -0,0 +1,85 @@
+// Copyright CERN and copyright holders of ALICE O2. This software is
+// distributed under the terms of the GNU General Public License v3 (GPL
+// Version 3), copied verbatim in the file "COPYING".
+//
+// See http://alice-o2.web.cern.ch/license for full licensing information.
+//
+// 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.
+
+#include <TSystem.h>
+#include <TTree.h>
+#include "Framework/Logger.h"
+#include "DataFormatsFT0/Digit.h"
+#include <TFile.h>
+#include <cstring>
+
+using namespace o2::ft0;
+
+int main(int argc, char* argv[])
+{
+  const std::string genDigDile{"ft0digits.root"};
+  const std::string decDigDile{"o2digit_ft0.root"};
+  const std::string branchBC{"FT0DIGITSBC"};
+  const std::string branchCH{"FT0DIGITSCH"};
+
+  if (gSystem->AccessPathName(genDigDile.c_str())) {
+    LOG(FATAL) << "Generated digits file " << genDigDile << " is absent";
+  }
+  TFile flIn(genDigDile.c_str());
+  std::unique_ptr<TTree> tree((TTree*)flIn.Get("o2sim"));
+  if (!flIn.IsOpen() || flIn.IsZombie() || !tree) {
+    LOG(FATAL) << "Failed to get tree from generated digits file " << genDigDile;
+  }
+  std::vector<o2::ft0::Digit> digitsBC, *ft0BCDataPtr = &digitsBC;
+  std::vector<o2::ft0::ChannelData> digitsCh, *ft0ChDataPtr = &digitsCh;
+  tree->SetBranchAddress("FT0DIGITSBC", &ft0BCDataPtr);
+  tree->SetBranchAddress("FT0DIGITSCH", &ft0ChDataPtr);
+
+  if (gSystem->AccessPathName(decDigDile.c_str())) {
+    LOG(FATAL) << "Decoded digits file " << genDigDile << " is absent";
+  }
+
+  TFile flIn2(decDigDile.c_str());
+  std::unique_ptr<TTree> tree2((TTree*)flIn2.Get("o2sim"));
+  if (!flIn2.IsOpen() || flIn2.IsZombie() || !tree2) {
+    LOG(FATAL) << "Failed to get tree from decoded digits file " << genDigDile;
+  }
+  std::vector<o2::ft0::Digit> digitsBC2, *ft0BCDataPtr2 = &digitsBC2;
+  std::vector<o2::ft0::ChannelData> digitsCh2, *ft0ChDataPtr2 = &digitsCh2;
+  tree2->SetBranchAddress("FT0DIGITSBC", &ft0BCDataPtr2);
+  tree2->SetBranchAddress("FT0DIGITSCH", &ft0ChDataPtr2);
+
+  int nbc = 0, nbc2 = 0, nch = 0, nch2 = 0;
+  for (int ient = 0; ient < tree->GetEntries(); ient++) {
+    tree->GetEntry(ient);
+    int nbcEntry = digitsBC.size();
+    nbc += nbcEntry;
+    for (int ibc = 0; ibc < nbcEntry; ibc++) {
+      auto& bcd = digitsBC[ibc];
+      int bc = bcd.getBC();
+      auto channels = bcd.getBunchChannelData(digitsCh);
+      nch += channels.size();
+    }
+  }
+
+  for (int ient = 0; ient < tree2->GetEntries(); ient++) {
+    tree2->GetEntry(ient);
+    int nbc2Entry = digitsBC2.size();
+    nbc2 += nbc2Entry;
+    for (int ibc = 0; ibc < nbc2Entry; ibc++) {
+      auto& bcd2 = digitsBC2[ibc];
+      int bc2 = bcd2.getBC();
+      auto channels2 = bcd2.getBunchChannelData(digitsCh2);
+      nch2 += channels2.size();
+    }
+  }
+  LOG(INFO) << "FT0 simulated: " << nbc << " triggers with " << nch << " channels";
+  LOG(INFO) << "FT0 decoded  : " << nbc2 << " triggers with " << nch2 << " channels";
+  if (nbc != nbc2 || nch != nch2) {
+    LOG(FATAL) << "Mismatch between the number of encoded and decoded objects";
+  }
+
+  return 0;
+}

prodtests/full_system_test_ci_extra_tests.sh

@@ -1,3 +1,32 @@
 #!/bin/bash
+#
+# Set of extra tests which may run after the full_system_test.sh
+# Particularly, they use the files generated by the full_system_test.sh
+#
+. ${O2_ROOT}/share/scripts/jobutils.sh
 
+if [ "0$O2_ROOT" == "0" ]; then
+  eval "`alienv shell-helper`"
+  alienv --no-refresh load O2/latest
+fi
 
+NHBPERTF=${NHBPERTF:-128} # number of HBFs per orbit to expect in pregenerated files
+
+# Set general arguments
+ARGS_ALL="--session default --shm-segment-size $SHMSIZE"
+
+# below we add individual tests
+execname=""
+logfile=""
+
+# FT0
+# run  raw->digits conversion from the raw data create elsewhere
+execname=o2-raw-file-reader-workflow
+taskwrapper ${execname}.log $execname "$ARGS_ALL --configKeyValues \"HBFUtils.nHBFPerTF=${NHBPERTF}\" --detect-tf0 --input-conf raw/FT0/FT0raw.cfg | o2-ft0-flp-dpl-workflow  $ARGS_ALL -b"
+# and check its correctness
+execname=o2-ft0-test-raw-conversion
+taskwrapper ${execname}.log $execname
+# don't proceed to other tests if error was found
+[ ! -f "${execname}.log_done" ] || [ "$?" != 0 ] && exit
+
+# other tests to add