testTPCCalDet.cxx

Go to the documentation of this file.

// 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.
 
#include <algorithm>
#define BOOST_TEST_MODULE Test TPC CalDet class
#define BOOST_TEST_MAIN
#define BOOST_TEST_DYN_LINK
#include <boost/range/combine.hpp>
#include <boost/test/unit_test.hpp>
#include <vector>
#include <limits>
#include <cstdlib>
 
#include "TMath.h"
#include "TPCBase/Mapper.h"
#include "TPCBase/CalArray.h"
#include "TPCBase/CalDet.h"
#include "TFile.h"
#include "Framework/TypeTraits.h"
#include "TPCBaseRecSim/DeadChannelMapCreator.h"
 
namespace o2::tpc
{
 
// templated euqality check
//  for integer one would need a specialisation to check for == instead of <
template <typename T>
bool isEqualAbs(T x, T y, int n = 1)
{
  // Since `epsilon()` is the gap size (ULP, unit in the last place)
  // of floating-point numbers in interval [1, 2), we can scale it to
  // the gap size in interval [2^e, 2^{e+1}), where `e` is the exponent
  // of `x` and `y`.
 
  // If `x` and `y` have different gap sizes (which means they have
  // different exponents), we take the smaller one. Taking the bigger
  // one is also reasonable, I guess.
  const T m = std::min(std::fabs(x), std::fabs(y));
 
  // Subnormal numbers have fixed exponent, which is `min_exponent - 1`.
  const int exp = m < std::numeric_limits<T>::min()
                    ? std::numeric_limits<T>::min_exponent - 1
                    : std::ilogb(m);
 
  // We consider `x` and `y` equal if the difference between them is
  // within `n` ULPs.
  return std::fabs(x - y) <= n * std::ldexp(std::numeric_limits<T>::epsilon(), exp);
}
 
template <typename T>
  requires(std::integral<T>)
bool isEqualAbs(T val1, T val2)
{
  return val1 == val2;
}
 
BOOST_AUTO_TEST_CASE(CalArray_ROOTIO)
{
  // CalROC roc(PadSubset::ROC, 10);
  CalArray<unsigned> roc(PadSubset::ROC, 10);
 
  int iter = 0;
  // unsigned iter=0;
  for (auto& val : roc.getData()) {
    val = iter++;
  }
 
  auto f = TFile::Open("CalArray_ROOTIO.root", "recreate");
  f->WriteObject(&roc, "roc");
  delete f;
 
  // CalROC *rocRead = nullptr;
  CalArray<unsigned>* rocRead = nullptr;
  f = TFile::Open("CalArray_ROOTIO.root");
  f->GetObject("roc", rocRead);
  delete f;
 
  BOOST_REQUIRE(rocRead != nullptr);
 
  float sumROC = 0;
  for (auto const& val : boost::combine(roc.getData(), rocRead->getData())) {
    sumROC += (val.get<0>() - val.get<1>());
  }
 
  BOOST_CHECK_CLOSE(sumROC, 0., 1.E-12);
}
 
BOOST_AUTO_TEST_CASE(CalDet_ROOTIO)
{
 
  auto& mapper = Mapper::instance();
  const auto numberOfPads = mapper.getPadsInSector() * 36;
 
  CalPad padROC(PadSubset::ROC);
  CalPad padPartition(PadSubset::Partition);
  CalPad padRegion(PadSubset::Region);
 
  // ===| Fill Data |===========================================================
  int iter = 0;
  // --- ROC type
  padROC.setName("ROCData");
  for (auto& calArray : padROC.getData()) {
    for (auto& value : calArray.getData()) {
      value = iter++;
    }
  }
 
  // --- Partition type
  padPartition.setName("PartitionData");
  for (auto& calArray : padPartition.getData()) {
    for (auto& value : calArray.getData()) {
      value = iter++;
    }
  }
 
  // --- Region type
  padRegion.setName("RegionData");
  for (auto& calArray : padRegion.getData()) {
    for (auto& value : calArray.getData()) {
      value = iter++;
    }
  }
 
  // ===| dump all objects to file |============================================
  auto f = TFile::Open("CalDet.root", "recreate");
  f->WriteObject(&padROC, "CalDetROC");
  f->WriteObject(&padPartition, "CalDetPartition");
  f->WriteObject(&padRegion, "CalDetRegion");
  f->Close();
  delete f;
 
  // ===| read back all values |================================================
  CalPad* padROCRead = nullptr;
  CalPad* padPartitionRead = nullptr;
  CalPad* padRegionRead = nullptr;
 
  f = TFile::Open("CalDet.root");
  f->GetObject("CalDetROC", padROCRead);
  f->GetObject("CalDetPartition", padPartitionRead);
  f->GetObject("CalDetRegion", padRegionRead);
 
  delete f;
 
  BOOST_REQUIRE(padROCRead != nullptr);
  BOOST_REQUIRE(padPartitionRead != nullptr);
  BOOST_REQUIRE(padRegionRead != nullptr);
 
  // ===| compare values before and after |=====================================
  float sumROC = 0.f;
  float sumPartition = 0.f;
  float sumRegion = 0.f;
 
  int numberOfPadsROC = 0;
  int numberOfPadsPartition = 0;
  int numberOfPadsRegion = 0;
 
  for (auto const& arrays : boost::combine(padROC.getData(), padROCRead->getData())) {
    for (auto const& val : boost::combine(arrays.get<0>().getData(), arrays.get<1>().getData())) {
      sumROC += (val.get<0>() - val.get<1>());
      ++numberOfPadsROC;
    }
  }
 
  for (auto const& arrays : boost::combine(padPartition.getData(), padPartitionRead->getData())) {
    for (auto const& val : boost::combine(arrays.get<0>().getData(), arrays.get<1>().getData())) {
      sumPartition += (val.get<0>() - val.get<1>());
      ++numberOfPadsPartition;
    }
  }
 
  for (auto const& arrays : boost::combine(padRegion.getData(), padRegionRead->getData())) {
    for (auto const& val : boost::combine(arrays.get<0>().getData(), arrays.get<1>().getData())) {
      sumRegion += (val.get<0>() - val.get<1>());
      ++numberOfPadsRegion;
    }
  }
 
  // ===| checks |==============================================================
  BOOST_CHECK_EQUAL(padROC.getName(), padROCRead->getName());
  BOOST_CHECK_CLOSE(sumROC, 0.f, 1.E-12);
  BOOST_CHECK_EQUAL(numberOfPadsROC, numberOfPads);
 
  BOOST_CHECK_EQUAL(padPartition.getName(), padPartitionRead->getName());
  BOOST_CHECK_CLOSE(sumPartition, 0.f, 1.E-12);
  BOOST_CHECK_EQUAL(numberOfPadsPartition, numberOfPads);
 
  BOOST_CHECK_EQUAL(padRegion.getName(), padRegionRead->getName());
  BOOST_CHECK_CLOSE(sumRegion, 0.f, 1.E-12);
  BOOST_CHECK_EQUAL(numberOfPadsRegion, numberOfPads);
} // BOOST_AUTO_TEST_CASE
 
BOOST_AUTO_TEST_CASE(CalDet_Arithmetics)
{
  // data
  CalPad pad(PadSubset::ROC);
 
  // data 2 for testing operators on objects
  CalPad pad2(PadSubset::ROC);
 
  // for applying the operators on
  CalPad padCmp(PadSubset::ROC);
 
  // ===| fill with data |======================================================
  int iter = 0;
  // --- ROC type
  for (auto& calArray : pad.getData()) {
    for (auto& value : calArray.getData()) {
      value = iter++;
    }
  }
 
  iter = 1;
  for (auto& calArray : pad2.getData()) {
    for (auto& value : calArray.getData()) {
      value = iter++;
    }
  }
 
  //
  // ===| test operators with simple numbers |==================================
  //
  const float number = 0.2f;
  bool isEqual = true;
 
  // + operator
  isEqual = true;
  padCmp = pad;
  padCmp += number;
 
  for (auto const& arrays : boost::combine(padCmp.getData(), pad.getData())) {
    for (auto const& val : boost::combine(arrays.get<0>().getData(), arrays.get<1>().getData())) {
      isEqual &= isEqualAbs(val.get<0>(), val.get<1>() + number);
    }
  }
  BOOST_CHECK_EQUAL(isEqual, true);
 
  // - operator
  isEqual = true;
  padCmp = pad;
  padCmp -= number;
 
  for (auto const& arrays : boost::combine(padCmp.getData(), pad.getData())) {
    for (auto const& val : boost::combine(arrays.get<0>().getData(), arrays.get<1>().getData())) {
      isEqual &= isEqualAbs(val.get<0>(), val.get<1>() - number);
    }
  }
  BOOST_CHECK_EQUAL(isEqual, true);
 
  // * operator
  isEqual = true;
  padCmp = pad;
  padCmp *= number;
 
  for (auto const& arrays : boost::combine(padCmp.getData(), pad.getData())) {
    for (auto const& val : boost::combine(arrays.get<0>().getData(), arrays.get<1>().getData())) {
      isEqual &= isEqualAbs(val.get<0>(), val.get<1>() * number);
    }
  }
  BOOST_CHECK_EQUAL(isEqual, true);
 
  // / operator
  isEqual = true;
  padCmp = pad;
  padCmp /= number;
 
  for (auto const& arrays : boost::combine(padCmp.getData(), pad.getData())) {
    for (auto const& val : boost::combine(arrays.get<0>().getData(), arrays.get<1>().getData())) {
      isEqual &= isEqualAbs(val.get<0>(), val.get<1>() / number);
    }
  }
  BOOST_CHECK_EQUAL(isEqual, true);
 
  //
  // ===| test operators with full object |=====================================
  //
  // + operator
  isEqual = true;
  padCmp = pad;
  padCmp += pad2;
 
  for (auto itpad = pad.getData().begin(), itpad2 = pad2.getData().begin(), itpadCmp = padCmp.getData().begin(); itpad != pad.getData().end(); ++itpad, ++itpad2, ++itpadCmp) {
    for (auto itval1 = (*itpad).getData().begin(), itval2 = (*itpad2).getData().begin(), itval3 = (*itpadCmp).getData().begin(); itval1 != (*itpad).getData().end(); ++itval1, ++itval2, ++itval3) {
      isEqual &= isEqualAbs(*itval3, *itval1 + *itval2);
    }
  }
  BOOST_CHECK_EQUAL(isEqual, true);
 
  // - operator
  isEqual = true;
  padCmp = pad;
  padCmp -= pad2;
 
  for (auto itpad = pad.getData().begin(), itpad2 = pad2.getData().begin(), itpadCmp = padCmp.getData().begin(); itpad != pad.getData().end(); ++itpad, ++itpad2, ++itpadCmp) {
    for (auto itval1 = (*itpad).getData().begin(), itval2 = (*itpad2).getData().begin(), itval3 = (*itpadCmp).getData().begin(); itval1 != (*itpad).getData().end(); ++itval1, ++itval2, ++itval3) {
      isEqual &= isEqualAbs(*itval3, *itval1 - *itval2);
    }
  }
  BOOST_CHECK_EQUAL(isEqual, true);
 
  // * operator
  isEqual = true;
  padCmp = pad;
  padCmp *= pad2;
 
  for (auto itpad = pad.getData().begin(), itpad2 = pad2.getData().begin(), itpadCmp = padCmp.getData().begin(); itpad != pad.getData().end(); ++itpad, ++itpad2, ++itpadCmp) {
    for (auto itval1 = (*itpad).getData().begin(), itval2 = (*itpad2).getData().begin(), itval3 = (*itpadCmp).getData().begin(); itval1 != (*itpad).getData().end(); ++itval1, ++itval2, ++itval3) {
      isEqual &= isEqualAbs(*itval3, *itval1 * *itval2);
    }
  }
  BOOST_CHECK_EQUAL(isEqual, true);
 
  // / operator
  isEqual = true;
  padCmp = pad;
  padCmp /= pad2;
 
  for (auto itpad = pad.getData().begin(), itpad2 = pad2.getData().begin(), itpadCmp = padCmp.getData().begin(); itpad != pad.getData().end(); ++itpad, ++itpad2, ++itpadCmp) {
    for (auto itval1 = (*itpad).getData().begin(), itval2 = (*itpad2).getData().begin(), itval3 = (*itpadCmp).getData().begin(); itval1 != (*itpad).getData().end(); ++itval1, ++itval2, ++itval3) {
      isEqual &= isEqualAbs(*itval3, *itval1 / *itval2);
    }
  }
  BOOST_CHECK_EQUAL(isEqual, true);
 
  // = operator
  isEqual = true;
  padCmp = 10.f;
  for (const auto& calArr : padCmp.getData()) {
    isEqual &= std::all_of(calArr.getData().begin(), calArr.getData().end(), [](const auto val) { return isEqualAbs(val, 10.f); });
  }
  BOOST_CHECK_EQUAL(isEqual, true);
}
 
BOOST_AUTO_TEST_CASE(CalDetTypeTest)
{
  using namespace o2::framework;
  BOOST_CHECK_EQUAL(has_root_dictionary<o2::tpc::CalDet<float>>::value, true);
  auto testDict = has_root_dictionary_mapped_type<std::unordered_map<std::string, o2::tpc::CalDet<float>>>::value;
  BOOST_CHECK(testDict == true);
}
 
BOOST_AUTO_TEST_CASE(CalDetStreamerTest)
{
  // simple code executing the TPC IDCPadFlags loading in a standalone env --> easy to valgrind
  //
  // Deliberately NOT ALICEO2_CCDB_HOST: that names the writable test instance,
  // which holds a *different* object at this path -- the timestamp below is
  // pinned to the production object's validity. The variable lets a build
  // container reach production through a broker; unset, behaviour is unchanged.
  const char* productionHost = std::getenv("ALICEO2_CCDB_PRODUCTION_HOST");
  o2::tpc::DeadChannelMapCreator creator{};
  creator.init((productionHost && *productionHost) ? productionHost : "https://alice-ccdb.cern.ch");
  creator.loadIDCPadFlags(1731274461770);
}
 
} // namespace o2::tpc