CMS 3D CMS Logo

List of all members | Classes | Public Member Functions | Private Attributes
HcalConvertedPedestalWidthsGPU Class Reference

#include <HcalConvertedPedestalWidthsGPU.h>

Inheritance diagram for HcalConvertedPedestalWidthsGPU:
HcalConvertedEffectivePedestalWidthsGPU

Classes

struct  Product
 

Public Member Functions

Product const & getProduct (cudaStream_t) const
 
 HcalConvertedPedestalWidthsGPU (HcalPedestals const &, HcalPedestalWidths const &, HcalQIEData const &, HcalQIETypes const &)
 
 ~HcalConvertedPedestalWidthsGPU ()=default
 

Private Attributes

cms::cuda::ESProduct< Productproduct_
 
uint64_t totalChannels_
 
std::vector< float, cms::cuda::HostAllocator< float > > values_
 

Detailed Description

Definition at line 16 of file HcalConvertedPedestalWidthsGPU.h.

Constructor & Destructor Documentation

◆ HcalConvertedPedestalWidthsGPU()

HcalConvertedPedestalWidthsGPU::HcalConvertedPedestalWidthsGPU ( HcalPedestals const &  pedestals,
HcalPedestalWidths const &  pedestalWidths,
HcalQIEData const &  qieData,
HcalQIETypes const &  qieTypes 
)

Definition at line 22 of file HcalConvertedPedestalWidthsGPU.cc.

References HcalCondObjectContainer< Item >::getAllContainers().

26  : totalChannels_{pedestals.getAllContainers()[0].second.size() + pedestals.getAllContainers()[1].second.size()},
27  values_(totalChannels_ * 4) {
28 #ifdef HCAL_MAHI_CPUDEBUG
29  std::cout << "hello from converted pedestal widths" << std::endl;
30  std::cout << "pedestals HB values = " << pedestals.getAllContainers()[0].second.size()
31  << " HE values = " << pedestals.getAllContainers()[1].second.size() << std::endl;
32  std::cout << "qiedata HB values = " << qieData.getAllContainers()[0].second.size()
33  << " HE values = " << qieData.getAllContainers()[1].second.size() << std::endl;
34 #endif
35 
36  // retrieve all collections
37  auto const pedestalsAll = pedestals.getAllContainers();
38  auto const pedestalWidthsAll = pedestalWidths.getAllContainers();
39  auto const qieDataAll = qieData.getAllContainers();
40  auto const qieTypesAll = qieTypes.getAllContainers();
41 
42  // have to convert to fc if stored in adc
43  auto const unitIsADC = pedestals.isADC();
44 
45  // fill in barrel
46  auto const& pedestalBarrelValues = pedestalsAll[0].second;
47  auto const& pedestalWidthBarrelValues = pedestalWidthsAll[0].second;
48  auto const& qieDataBarrelValues = qieDataAll[0].second;
49  auto const& qieTypesBarrelValues = qieTypesAll[0].second;
50 
51 #ifdef HCAL_MAHI_CPUDEBUG
52  assert(pedestalWidthBarrelValues.size() == pedestalBarrelValues.size());
53  assert(pedestalBarrelValues.size() == qieDataBarrelValues.size());
54  assert(pedestalBarrelValues.size() == qieTypesBarrelValues.size());
55 #endif
56 
57  for (uint64_t i = 0; i < pedestalBarrelValues.size(); ++i) {
58  auto const& qieCoder = qieDataBarrelValues[i];
59  auto const qieType = qieTypesBarrelValues[i].getValue() > 1 ? 1 : 0;
60  auto const& qieShape = qieData.getShape(qieType);
61 
62  values_[i * 4] =
63  unitIsADC
64  ? convert(
65  pedestalBarrelValues[i].getValue(0), pedestalWidthBarrelValues[i].getWidth(0), 0, qieCoder, qieShape)
66  : pedestalWidthBarrelValues[i].getWidth(0);
67  values_[i * 4 + 1] =
68  unitIsADC
69  ? convert(
70  pedestalBarrelValues[i].getValue(1), pedestalWidthBarrelValues[i].getWidth(1), 1, qieCoder, qieShape)
71  : pedestalWidthBarrelValues[i].getWidth(1);
72  values_[i * 4 + 2] =
73  unitIsADC
74  ? convert(
75  pedestalBarrelValues[i].getValue(2), pedestalWidthBarrelValues[i].getWidth(2), 2, qieCoder, qieShape)
76  : pedestalWidthBarrelValues[i].getWidth(2);
77  values_[i * 4 + 3] =
78  unitIsADC
79  ? convert(
80  pedestalBarrelValues[i].getValue(3), pedestalWidthBarrelValues[i].getWidth(3), 3, qieCoder, qieShape)
81  : pedestalWidthBarrelValues[i].getWidth(3);
82  }
83 
84  // fill in endcap
85  auto const& pedestalEndcapValues = pedestalsAll[1].second;
86  auto const& pedestalWidthEndcapValues = pedestalWidthsAll[1].second;
87  auto const& qieDataEndcapValues = qieDataAll[1].second;
88  auto const& qieTypesEndcapValues = qieTypesAll[1].second;
89 
90 #ifdef HCAL_MAHI_CPUDEBUG
91  assert(pedestalWidthEndcapValues.size() == pedestalEndcapValues.size());
92  assert(pedestalEndcapValues.size() == qieDataEndcapValues.size());
93  assert(pedestalEndcapValues.size() == qieTypesEndcapValues.size());
94 #endif
95 
96  auto const offset = pedestalWidthBarrelValues.size();
97  for (uint64_t i = 0; i < pedestalEndcapValues.size(); ++i) {
98  auto const& qieCoder = qieDataEndcapValues[i];
99  auto const qieType = qieTypesEndcapValues[i].getValue() > 1 ? 1 : 0;
100  auto const& qieShape = qieData.getShape(qieType);
101  auto const off = offset + i;
102 
103  values_[off * 4] =
104  unitIsADC
105  ? convert(
106  pedestalEndcapValues[i].getValue(0), pedestalWidthEndcapValues[i].getWidth(0), 0, qieCoder, qieShape)
107  : pedestalWidthEndcapValues[i].getWidth(0);
108  values_[off * 4 + 1] =
109  unitIsADC
110  ? convert(
111  pedestalEndcapValues[i].getValue(1), pedestalWidthEndcapValues[i].getWidth(1), 1, qieCoder, qieShape)
112  : pedestalWidthEndcapValues[i].getWidth(1);
113  values_[off * 4 + 2] =
114  unitIsADC
115  ? convert(
116  pedestalEndcapValues[i].getValue(2), pedestalWidthEndcapValues[i].getWidth(2), 2, qieCoder, qieShape)
117  : pedestalWidthEndcapValues[i].getWidth(2);
118  values_[off * 4 + 3] =
119  unitIsADC
120  ? convert(
121  pedestalEndcapValues[i].getValue(3), pedestalWidthEndcapValues[i].getWidth(3), 3, qieCoder, qieShape)
122  : pedestalWidthEndcapValues[i].getWidth(3);
123 
124 #ifdef HCAL_MAHI_CPUDEBUG
125  if (pedestalEndcapValues[i].rawId() == DETID_TO_DEBUG) {
126  for (int i = 0; i < 4; i++)
127  printf("pedestalWidth(%d) = %f original pedestalWidth(%d) = %f\n",
128  i,
129  values_[off * 4 + i],
130  i,
131  pedestalWidthEndcapValues[i].getWidth(3));
132  }
133 #endif
134  }
135 }
std::vector< float, cms::cuda::HostAllocator< float > > values_
assert(be >=bs)
def convert(infile, ofile)
unsigned long long uint64_t
Definition: Time.h:13
const JetExtendedData & getValue(const Container &, const reco::JetBaseRef &)
get value for the association. Throw exception if no association found

◆ ~HcalConvertedPedestalWidthsGPU()

HcalConvertedPedestalWidthsGPU::~HcalConvertedPedestalWidthsGPU ( )
default

Member Function Documentation

◆ getProduct()

HcalConvertedPedestalWidthsGPU::Product const & HcalConvertedPedestalWidthsGPU::getProduct ( cudaStream_t  stream) const

Definition at line 137 of file HcalConvertedPedestalWidthsGPU.cc.

References cms::cuda::copyAsync(), product_, cms::cuda::stream, HcalConvertedPedestalWidthsGPU::Product::values, and values_.

137  {
138  auto const& product = product_.dataForCurrentDeviceAsync(
139  stream, [this](HcalConvertedPedestalWidthsGPU::Product& product, cudaStream_t stream) {
140  // allocate
141  product.values = cms::cuda::make_device_unique<float[]>(values_.size(), stream);
142 
143  // transfer
145  });
146 
147  return product;
148 }
std::vector< float, cms::cuda::HostAllocator< float > > values_
uint32_t T const *__restrict__ uint32_t const *__restrict__ int32_t int Histo::index_type cudaStream_t stream
cms::cuda::ESProduct< Product > product_
void copyAsync(device::unique_ptr< T > &dst, const host::unique_ptr< T > &src, cudaStream_t stream)
Definition: copyAsync.h:20
edm::propagate_const_array< cms::cuda::device::unique_ptr< float[]> > values

Member Data Documentation

◆ product_

cms::cuda::ESProduct<Product> HcalConvertedPedestalWidthsGPU::product_
private

Definition at line 39 of file HcalConvertedPedestalWidthsGPU.h.

Referenced by getProduct().

◆ totalChannels_

uint64_t HcalConvertedPedestalWidthsGPU::totalChannels_
private

Definition at line 36 of file HcalConvertedPedestalWidthsGPU.h.

◆ values_

std::vector<float, cms::cuda::HostAllocator<float> > HcalConvertedPedestalWidthsGPU::values_
private

Definition at line 37 of file HcalConvertedPedestalWidthsGPU.h.

Referenced by getProduct().