HcalConvertedEffectivePedestalsGPU Class Reference

#include <HcalConvertedEffectivePedestalsGPU.h>

Inheritance diagram for HcalConvertedEffectivePedestalsGPU:

Public Member Functions
	HcalConvertedPedestalsGPU (HcalPedestals const &, HcalQIEData const &, HcalQIETypes const &)
Public Member Functions inherited from HcalConvertedPedestalsGPU
Product const &	getProduct (cudaStream_t) const
	HcalConvertedPedestalsGPU (HcalPedestals const &, HcalQIEData const &, HcalQIETypes const &)
uint32_t	offsetForHashes () const
	~HcalConvertedPedestalsGPU ()=default

Additional Inherited Members
Protected Attributes inherited from HcalConvertedPedestalsGPU
uint32_t	offsetForHashes_
cms::cuda::ESProduct< Product >	product_
uint64_t	totalChannels_
std::vector< float, cms::cuda::HostAllocator< float > >	values_

Detailed Description

Definition at line 9 of file HcalConvertedEffectivePedestalsGPU.h.

Member Function Documentation

HcalConvertedPedestalsGPU()

HcalConvertedPedestalsGPU::HcalConvertedPedestalsGPU

Definition at line 18 of file HcalConvertedPedestalsGPU.cc.

    : totalChannels_{pedestals.getAllContainers()[0].second.size() + pedestals.getAllContainers()[1].second.size()},
      offsetForHashes_{static_cast<uint32_t>(pedestals.getAllContainers()[0].second.size())},
      values_(totalChannels_ * 4) {
#ifdef HCAL_MAHI_CPUDEBUG
  std::cout << "hello from converted pedestals" << std::endl;
  std::cout << "pedestals HB values = " << pedestals.getAllContainers()[0].second.size()
            << "  HE values = " << pedestals.getAllContainers()[1].second.size() << std::endl;
  std::cout << "qiedata HB values = " << qieData.getAllContainers()[0].second.size()
            << "  HE values = " << qieData.getAllContainers()[1].second.size() << std::endl;
#endif

  // retrieve all collections
  auto const pedestalsAll = pedestals.getAllContainers();
  auto const qieDataAll = qieData.getAllContainers();
  auto const qieTypesAll = qieTypes.getAllContainers();

  // have to convert to fc if stored in adc
  auto const unitIsADC = pedestals.isADC();

  // fill in barrel
  auto const& pedestalBarrelValues = pedestalsAll[0].second;
  auto const& qieDataBarrelValues = qieDataAll[0].second;
  auto const& qieTypesBarrelValues = qieTypesAll[0].second;

#ifdef HCAL_MAHI_CPUDEBUG
  assert(pedestalBarrelValues.size() == qieDataBarrelValues.size());
  assert(pedestalBarrelValues.size() == qieTypesBarrelValues.size());
#endif

  for (uint64_t i = 0; i < pedestalBarrelValues.size(); ++i) {
    auto const& qieCoder = qieDataBarrelValues[i];
    auto const qieType = qieTypesBarrelValues[i].getValue() > 1 ? 1 : 0;
    auto const& qieShape = qieData.getShape(qieType);

    values_[i * 4] = unitIsADC ? convert(pedestalBarrelValues[i].getValue(0), 0, qieCoder, qieShape)
                               : pedestalBarrelValues[i].getValue(0);
    values_[i * 4 + 1] = unitIsADC ? convert(pedestalBarrelValues[i].getValue(1), 1, qieCoder, qieShape)
                                   : pedestalBarrelValues[i].getValue(1);
    values_[i * 4 + 2] = unitIsADC ? convert(pedestalBarrelValues[i].getValue(2), 2, qieCoder, qieShape)
                                   : pedestalBarrelValues[i].getValue(2);
    values_[i * 4 + 3] = unitIsADC ? convert(pedestalBarrelValues[i].getValue(3), 3, qieCoder, qieShape)
                                   : pedestalBarrelValues[i].getValue(3);
  }

  // fill in endcap
  auto const& pedestalEndcapValues = pedestalsAll[1].second;
  auto const& qieDataEndcapValues = qieDataAll[1].second;
  auto const& qieTypesEndcapValues = qieTypesAll[1].second;

#ifdef HCAL_MAHI_CPUDEBUG
  assert(pedestalEndcapValues.size() == qieDataEndcapValues.size());
  assert(pedestalEndcapValues.size() == qieTypesEndcapValues.size());
#endif

  auto const offset = pedestalBarrelValues.size();
  for (uint64_t i = 0; i < pedestalEndcapValues.size(); ++i) {
    auto const& qieCoder = qieDataEndcapValues[i];
    auto const qieType = qieTypesEndcapValues[i].getValue() > 1 ? 1 : 0;
    auto const& qieShape = qieData.getShape(qieType);
    auto const off = offset + i;

    values_[off * 4] = unitIsADC ? convert(pedestalEndcapValues[i].getValue(0), 0, qieCoder, qieShape)
                                 : pedestalEndcapValues[i].getValue(0);
    values_[off * 4 + 1] = unitIsADC ? convert(pedestalEndcapValues[i].getValue(1), 1, qieCoder, qieShape)
                                     : pedestalEndcapValues[i].getValue(1);
    values_[off * 4 + 2] = unitIsADC ? convert(pedestalEndcapValues[i].getValue(2), 2, qieCoder, qieShape)
                                     : pedestalEndcapValues[i].getValue(2);
    values_[off * 4 + 3] = unitIsADC ? convert(pedestalEndcapValues[i].getValue(3), 3, qieCoder, qieShape)
                                     : pedestalEndcapValues[i].getValue(3);

#ifdef HCAL_MAHI_CPUDEBUG
    if (pedestalEndcapValues[i].rawId() == DETID_TO_DEBUG) {
      printf("qietype = %d\n", qieType);
      printf("ped0 = %f ped1 = %f ped2 = %f ped3 = %f\n",
             pedestalEndcapValues[i].getValue(0),
             pedestalEndcapValues[i].getValue(1),
             pedestalEndcapValues[i].getValue(2),
             pedestalEndcapValues[i].getValue(3));
      printf("converted: ped0 = %f ped1 = %f ped2 = %f ped3 = %f\n",
             values_[off * 4],
             values_[off * 4 + 1],
             values_[off * 4 + 2],
             values_[off * 4 + 3]);
    }
#endif
  }
}