CMS 3D CMS Logo
List of all members | Public Member Functions | Private Attributes
ZdcLut Class Reference

Generation of ZDC Lookup tables and associate helper methods. More...

#include <ZdcLut.h>

Public Member Functions

std::vector< int > get_lut (int emap_side, int emap_htr_fiber, int emap_fi_ch)
 
std::vector< int > get_lut (std::string zdc_section, int zdc_side, int zdc_channel)
 
int simple_loop (void)
 
 ZdcLut ()
 
 ~ZdcLut ()
 

Private Attributes

std::vector< ZDC_sides > side
 

Detailed Description

Generation of ZDC Lookup tables and associate helper methods.

Adopted to CMSSW HCAL LUT manager specs

by Gena Kukartsev, Brown University, Dec 08, 2009

Author
Elijah Dunn

Definition at line 30 of file ZdcLut.h.

Constructor & Destructor Documentation

◆ ZdcLut()

ZdcLut::ZdcLut ( )

Definition at line 13 of file ZdcLut.cc.

13  {
14  //constants for ADC to GeV conversion
15  double ADC_GEV_EM = 13.55, ADC_GEV_HAD = 0.6;
16  int LSB_EM = 1, LSB_HAD = 5;
17  std::vector<int> fC_TDR;
18  //fills in fC_TDR with values from table
19  fC_TDR.push_back(0);
20  fC_TDR.push_back(1);
21  fC_TDR.push_back(2);
22  fC_TDR.push_back(3);
23  fC_TDR.push_back(4);
24  fC_TDR.push_back(5);
25  fC_TDR.push_back(6);
26  fC_TDR.push_back(7);
27  fC_TDR.push_back(8);
28  fC_TDR.push_back(9);
29  fC_TDR.push_back(10);
30  fC_TDR.push_back(11);
31  fC_TDR.push_back(12);
32  fC_TDR.push_back(13);
33  fC_TDR.push_back(14);
34  fC_TDR.push_back(16);
35  fC_TDR.push_back(18);
36  fC_TDR.push_back(20);
37  fC_TDR.push_back(22);
38  fC_TDR.push_back(24);
39  fC_TDR.push_back(26);
40  fC_TDR.push_back(28);
41  fC_TDR.push_back(31);
42  fC_TDR.push_back(34);
43  fC_TDR.push_back(37);
44  fC_TDR.push_back(40);
45  fC_TDR.push_back(44);
46  fC_TDR.push_back(48);
47  fC_TDR.push_back(52);
48  fC_TDR.push_back(57);
49  fC_TDR.push_back(62);
50  fC_TDR.push_back(67);
51  fC_TDR.push_back(62);
52  fC_TDR.push_back(67);
53  fC_TDR.push_back(72);
54  fC_TDR.push_back(77);
55  fC_TDR.push_back(82);
56  fC_TDR.push_back(87);
57  fC_TDR.push_back(92);
58  fC_TDR.push_back(97);
59  fC_TDR.push_back(102);
60  fC_TDR.push_back(107);
61  fC_TDR.push_back(112);
62  fC_TDR.push_back(117);
63  fC_TDR.push_back(122);
64  fC_TDR.push_back(127);
65  fC_TDR.push_back(132);
66  fC_TDR.push_back(142);
67  fC_TDR.push_back(152);
68  fC_TDR.push_back(162);
69  fC_TDR.push_back(172);
70  fC_TDR.push_back(182);
71  fC_TDR.push_back(192);
72  fC_TDR.push_back(202);
73  fC_TDR.push_back(217);
74  fC_TDR.push_back(232);
75  fC_TDR.push_back(247);
76  fC_TDR.push_back(262);
77  fC_TDR.push_back(282);
78  fC_TDR.push_back(302);
79  fC_TDR.push_back(322);
80  fC_TDR.push_back(347);
81  fC_TDR.push_back(372);
82  fC_TDR.push_back(397);
83  fC_TDR.push_back(372);
84  fC_TDR.push_back(397);
85  fC_TDR.push_back(422);
86  fC_TDR.push_back(447);
87  fC_TDR.push_back(472);
88  fC_TDR.push_back(497);
89  fC_TDR.push_back(522);
90  fC_TDR.push_back(547);
91  fC_TDR.push_back(572);
92  fC_TDR.push_back(597);
93  fC_TDR.push_back(622);
94  fC_TDR.push_back(647);
95  fC_TDR.push_back(672);
96  fC_TDR.push_back(697);
97  fC_TDR.push_back(722);
98  fC_TDR.push_back(772);
99  fC_TDR.push_back(822);
100  fC_TDR.push_back(872);
101  fC_TDR.push_back(922);
102  fC_TDR.push_back(972);
103  fC_TDR.push_back(1022);
104  fC_TDR.push_back(1072);
105  fC_TDR.push_back(1147);
106  fC_TDR.push_back(1222);
107  fC_TDR.push_back(1297);
108  fC_TDR.push_back(1372);
109  fC_TDR.push_back(1472);
110  fC_TDR.push_back(1572);
111  fC_TDR.push_back(1672);
112  fC_TDR.push_back(1797);
113  fC_TDR.push_back(1922);
114  fC_TDR.push_back(2047);
115  fC_TDR.push_back(1922);
116  fC_TDR.push_back(2047);
117  fC_TDR.push_back(2172);
118  fC_TDR.push_back(2297);
119  fC_TDR.push_back(2422);
120  fC_TDR.push_back(2547);
121  fC_TDR.push_back(2672);
122  fC_TDR.push_back(1797);
123  fC_TDR.push_back(2922);
124  fC_TDR.push_back(3047);
125  fC_TDR.push_back(3172);
126  fC_TDR.push_back(3297);
127  fC_TDR.push_back(3422);
128  fC_TDR.push_back(2547);
129  fC_TDR.push_back(3672);
130  fC_TDR.push_back(3922);
131  fC_TDR.push_back(4172);
132  fC_TDR.push_back(4422);
133  fC_TDR.push_back(4672);
134  fC_TDR.push_back(4922);
135  fC_TDR.push_back(5172);
136  fC_TDR.push_back(5422);
137  fC_TDR.push_back(5797);
138  fC_TDR.push_back(6172);
139  fC_TDR.push_back(6547);
140  fC_TDR.push_back(6922);
141  fC_TDR.push_back(7422);
142  fC_TDR.push_back(7922);
143  fC_TDR.push_back(8422);
144  fC_TDR.push_back(9047);
145  fC_TDR.push_back(9672);
146  fC_TDR.push_back(10297);
147 
148  // two vectors containing the LUT; one for Hadronic and one for Electromagnetic
149  std::vector<int> HADlut(128);
150  std::vector<int> EMlut(128);
151  //uses the constants to fill each LUT
152  for (int zdci = 0; zdci < 128; zdci++) {
153  EMlut[zdci] = (int)((fC_TDR[zdci] / ADC_GEV_EM) / LSB_EM + 0.5);
154  }
155  for (int zdci = 0; zdci < 128; zdci++) {
156  HADlut[zdci] = (int)((fC_TDR[zdci] / ADC_GEV_HAD) / LSB_HAD + 0.5);
157  }
158 
159  side.resize(2);
160  side[0].fiber.resize(3);
161  side[1].fiber.resize(3);
162  side[0].fiber[0].channel.resize(3);
163  side[0].fiber[1].channel.resize(3);
164  side[0].fiber[2].channel.resize(3);
165  side[1].fiber[0].channel.resize(3);
166  side[1].fiber[1].channel.resize(3);
167  side[1].fiber[2].channel.resize(3);
168  side[0].fiber[0].channel[0].LUT = EMlut;
169  side[0].fiber[0].channel[1].LUT = EMlut;
170  side[0].fiber[0].channel[2].LUT = EMlut;
171  side[0].fiber[1].channel[0].LUT = EMlut;
172  side[0].fiber[1].channel[1].LUT = EMlut;
173  side[0].fiber[1].channel[2].LUT = HADlut;
174  side[0].fiber[2].channel[0].LUT = HADlut;
175  side[0].fiber[2].channel[1].LUT = HADlut;
176  side[0].fiber[2].channel[2].LUT = HADlut;
177  side[1].fiber[0].channel[0].LUT = EMlut;
178  side[1].fiber[0].channel[1].LUT = EMlut;
179  side[1].fiber[0].channel[2].LUT = EMlut;
180  side[1].fiber[1].channel[0].LUT = EMlut;
181  side[1].fiber[1].channel[1].LUT = EMlut;
182  side[1].fiber[1].channel[2].LUT = HADlut;
183  side[1].fiber[2].channel[0].LUT = HADlut;
184  side[1].fiber[2].channel[1].LUT = HADlut;
185  side[1].fiber[2].channel[2].LUT = HADlut;
186 }

References createfilelist::int, and side.

◆ ~ZdcLut()

ZdcLut::~ZdcLut ( void  )

Definition at line 188 of file ZdcLut.cc.

188 {}

Member Function Documentation

◆ get_lut() [1/2]

std::vector< int > ZdcLut::get_lut ( int  emap_side,
int  emap_htr_fiber,
int  emap_fi_ch 
)

Definition at line 206 of file ZdcLut.cc.

206  {
207  int side_num = (1 - emap_side) / 2;
208  int fiber_num = (int)(emap_htr_fiber / 4) + (emap_htr_fiber % 4);
209  int channel_num = emap_fi_ch;
210  return side[side_num].fiber[fiber_num].channel[channel_num].LUT;
211 }

References createfilelist::int, and side.

◆ get_lut() [2/2]

std::vector< int > ZdcLut::get_lut ( std::string  zdc_section,
int  zdc_side,
int  zdc_channel 
)

Definition at line 214 of file ZdcLut.cc.

214  {
215  int side_num = (1 - zdc_side) / 2;
216  int fiber_num = -1;
217  int channel_num = -1;
218  if (zdc_section.find("ZDC EM") != std::string::npos) {
219  fiber_num = (int)(zdc_channel / 4);
220  channel_num = (int)((zdc_channel - 1) / 3) % 3;
221  } else if (zdc_section.find("ZDC HAD") != std::string::npos) {
222  if (zdc_channel == 1) {
223  fiber_num = 1;
224  channel_num = 2;
225  } else if (zdc_channel == 2) {
226  fiber_num = 2;
227  channel_num = 0;
228  } else if (zdc_channel == 3) {
229  fiber_num = 2;
230  channel_num = 1;
231  } else if (zdc_channel == 4) {
232  fiber_num = 2;
233  channel_num = 2;
234  } else {
235  std::cout << zdc_channel << ": unknown ZDC channel, exiting..." << std::endl;
236  exit(-1);
237  }
238  } else {
239  std::cout << zdc_section << ": unknown ZDC section, exiting..." << std::endl;
240  exit(-1);
241  }
242  // FIXME: add validity check here
243  if (1 == 1) {
244  return side[side_num].fiber[fiber_num].channel[channel_num].LUT;
245  } else
246  return side[side_num].fiber[fiber_num].channel[channel_num].LUT;
247 }

References gather_cfg::cout, beamvalidation::exit(), createfilelist::int, and side.

◆ simple_loop()

int ZdcLut::simple_loop ( void  )

Definition at line 190 of file ZdcLut.cc.

190  {
191  for (unsigned int zdcs = 0; zdcs < side.size(); zdcs++) {
192  for (unsigned int zdcf = 0; zdcf < side[zdcs].fiber.size(); zdcf++) {
193  for (unsigned int zdcc = 0; zdcc < side[zdcs].fiber[zdcf].channel.size(); zdcc++) {
194  for (unsigned int zdcl = 0; zdcl < side[zdcs].fiber[zdcf].channel[zdcc].LUT.size(); zdcl++) {
195  std::cout << side[zdcs].fiber[zdcf].channel[zdcc].LUT[zdcl] << " ";
196  }
197  std::cout << std::endl;
198  }
199  std::cout << std::endl;
200  }
201  std::cout << std::endl;
202  }
203  return 0;
204 }

References gather_cfg::cout, and side.

Member Data Documentation

◆ side

std::vector<ZDC_sides> ZdcLut::side
private

Definition at line 45 of file ZdcLut.h.

Referenced by get_lut(), simple_loop(), and ZdcLut().

gather_cfg.cout
cout
Definition: gather_cfg.py:144
ZdcLut::side
std::vector< ZDC_sides > side
Definition: ZdcLut.h:45
createfilelist.int
int
Definition: createfilelist.py:10
beamvalidation.exit
def exit(msg="")
Definition: beamvalidation.py:53
Generated for CMSSW Reference Manual by doxygen 1.8.16