CMS 3D CMS Logo

All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Modules Pages
List of all members | Public Member Functions
HcalConvertedEffectivePedestalWidthsGPU Class Referencefinal

#include <HcalConvertedEffectivePedestalWidthsGPU.h>

Inheritance diagram for HcalConvertedEffectivePedestalWidthsGPU:
HcalConvertedPedestalWidthsGPU

Public Member Functions

 HcalConvertedPedestalWidthsGPU (HcalPedestals const &, HcalPedestalWidths const &, HcalQIEData const &, HcalQIETypes const &)
 
- Public Member Functions inherited from HcalConvertedPedestalWidthsGPU
Product const & getProduct (cudaStream_t) const
 
 HcalConvertedPedestalWidthsGPU (HcalPedestals const &, HcalPedestalWidths const &, HcalQIEData const &, HcalQIETypes const &)
 
 ~HcalConvertedPedestalWidthsGPU ()=default
 

Detailed Description

Definition at line 7 of file HcalConvertedEffectivePedestalWidthsGPU.h.

Member Function Documentation

◆ HcalConvertedPedestalWidthsGPU()

HcalConvertedPedestalWidthsGPU::HcalConvertedPedestalWidthsGPU

Definition at line 22 of file HcalConvertedPedestalWidthsGPU.cc.

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