CMS 3D CMS Logo
ZDCTask.cc
Go to the documentation of this file.
1 
2 #include "DQM/HcalTasks/interface/ZDCTask.h"
3 #include <map>
4 
5 using namespace hcaldqm;
6 using namespace hcaldqm::constants;
7 ZDCTask::ZDCTask(edm::ParameterSet const& ps) {
8  // tags
9  _tagQIE10 = ps.getUntrackedParameter<edm::InputTag>("tagQIE10", edm::InputTag("hcalDigis"));
10  _tokQIE10 = consumes<ZDCDigiCollection>(_tagQIE10);
11 
12  // cuts
13  _cut = ps.getUntrackedParameter<double>("cut", 50.0);
14  _ped = ps.getUntrackedParameter<int>("ped", 4);
15 }
16 
17 /* virtual */ void ZDCTask::bookHistograms(DQMStore::IBooker& ib, edm::Run const& r, edm::EventSetup const& es) {
18  //############################## hardcode manually the zdc mapping #############################
19  //############################# this follows from https://github.com/cms-sw/cmssw/blob/CMSSW_8_0_X/EventFilter/CastorRawToDigi/src/ZdcUnpacker.cc#L118
20  //##############################################################################################
22  //PZDC
23  std::map<HcalElectronicsId, DetId> myEMap;
24  HcalElectronicsId eid = HcalElectronicsId(0, 1, 0, 3);
25  eid.setHTR(18, 8, 1);
26  myEMap[eid] = DetId(0x54000051); //PZDC EM1
27 
28  eid = HcalElectronicsId(1, 1, 0, 3);
29  eid.setHTR(18, 8, 1);
30  myEMap[eid] = DetId(0x54000052); //PZDC EM2
31 
32  eid = HcalElectronicsId(2, 1, 0, 3);
33  eid.setHTR(18, 8, 1);
34  myEMap[eid] = DetId(0x54000053); //PZDC EM3
35 
36  eid = HcalElectronicsId(0, 2, 0, 3);
37  eid.setHTR(18, 8, 1);
38  myEMap[eid] = DetId(0x54000061); //PZDC HAD1
39 
40  eid = HcalElectronicsId(1, 2, 0, 3);
41  eid.setHTR(18, 8, 1);
42  myEMap[eid] = DetId(0x54000054); //PZDC EM4
43 
44  eid = HcalElectronicsId(2, 2, 0, 3);
45  eid.setHTR(18, 8, 1);
46  myEMap[eid] = DetId(0x54000055); //PZDC EM5
47 
48  eid = HcalElectronicsId(0, 3, 0, 3);
49  eid.setHTR(18, 8, 1);
50  myEMap[eid] = DetId(0x54000062); //PZDC HAD2
51 
52  eid = HcalElectronicsId(1, 3, 0, 3);
53  eid.setHTR(18, 8, 1);
54  myEMap[eid] = DetId(0x54000063); //PZDC HAD3
55 
56  eid = HcalElectronicsId(2, 3, 0, 3);
57  eid.setHTR(18, 8, 1);
58  myEMap[eid] = DetId(0x54000064); //PZDC HAD4
59 
60  //NZDC
61  eid = HcalElectronicsId(0, 1, 1, 3);
62  eid.setHTR(18, 8, 0);
63  myEMap[eid] = DetId(0x54000011); //NZDC EM1
64 
65  eid = HcalElectronicsId(1, 1, 1, 3);
66  eid.setHTR(18, 8, 0);
67  myEMap[eid] = DetId(0x54000012); //NZDC EM2
68 
69  eid = HcalElectronicsId(2, 1, 1, 3);
70  eid.setHTR(18, 8, 0);
71  myEMap[eid] = DetId(0x54000013); //NZDC EM3
72 
73  eid = HcalElectronicsId(0, 2, 1, 3);
74  eid.setHTR(18, 8, 0);
75  myEMap[eid] = DetId(0x54000015); //NZDC EM5
76 
77  eid = HcalElectronicsId(1, 2, 1, 3);
78  eid.setHTR(18, 8, 0);
79  myEMap[eid] = DetId(0x54000021); //NZDC HAD1
80 
81  eid = HcalElectronicsId(2, 2, 1, 3);
82  eid.setHTR(18, 8, 0);
83  myEMap[eid] = DetId(0x54000014); //NZDC EM4
84 
85  eid = HcalElectronicsId(0, 3, 1, 3);
86  eid.setHTR(18, 8, 0);
87  myEMap[eid] = DetId(0x54000022); //NZDC HAD2
88 
89  eid = HcalElectronicsId(1, 3, 1, 3);
90  eid.setHTR(18, 8, 0);
91  myEMap[eid] = DetId(0x54000023); //NZDC HAD3
92 
93  eid = HcalElectronicsId(2, 3, 1, 3);
94  eid.setHTR(18, 8, 0);
95  myEMap[eid] = DetId(0x54000024); //NZDC HAD4
96  //##################################### end hardcoding ###################################
97 
98  ib.cd();
99 
100  //quantities for axis
101  hcaldqm::quantity::ValueQuantity xAxisShape(hcaldqm::quantity::fTiming_TS);
102  hcaldqm::quantity::ValueQuantity yAxisShape(hcaldqm::quantity::ffC_10000);
103 
104  hcaldqm::quantity::ValueQuantity xAxisADC(hcaldqm::quantity::fADC_128);
105 
106  //book histos per channel
107  for (std::map<HcalElectronicsId, DetId>::const_iterator itr = myEMap.begin(); itr != myEMap.end(); ++itr) {
108  char histoname[300];
109 
110  sprintf(histoname,
111  "%d_%d_%d_%d",
112  itr->first.fiberChanId(),
113  itr->first.fiberIndex(),
114  itr->first.spigot(),
115  itr->first.dccid());
116 
117  ib.setCurrentFolder("Hcal/ZDCTask/Shape_perChannel");
118  _cShape_EChannel[histoname] = ib.bookProfile(histoname,
119  histoname,
120  xAxisShape.nbins(),
121  xAxisShape.min(),
122  xAxisShape.max(),
123  yAxisShape.nbins(),
124  yAxisShape.min(),
125  yAxisShape.max());
126  _cShape_EChannel[histoname]->setAxisTitle("Timing", 1);
127  _cShape_EChannel[histoname]->setAxisTitle("fC QIE8", 2);
128 
129  ib.setCurrentFolder("Hcal/ZDCTask/ADC_perChannel");
130  _cADC_EChannel[histoname] = ib.book1D(histoname, histoname, xAxisADC.nbins(), xAxisADC.min(), xAxisADC.max());
131  _cADC_EChannel[histoname]->getRootObject()->SetBit(
132  BIT(hcaldqm::constants::BIT_OFFSET + hcaldqm::quantity::AxisType::fYAxis));
133  _cADC_EChannel[histoname]->setAxisTitle("ADC QIE8", 1);
134 
135  ib.setCurrentFolder("Hcal/ZDCTask/ADC_vs_TS_perChannel");
136  _cADC_vs_TS_EChannel[histoname] = ib.book2D(histoname,
137  histoname,
138  xAxisShape.nbins(),
139  xAxisShape.min(),
140  xAxisShape.max(),
141  xAxisADC.nbins(),
142  xAxisADC.min(),
143  xAxisADC.max());
144  _cADC_vs_TS_EChannel[histoname]->getRootObject()->SetBit(
145  BIT(hcaldqm::constants::BIT_OFFSET + hcaldqm::quantity::AxisType::fYAxis));
146  _cADC_vs_TS_EChannel[histoname]->setAxisTitle("Timing", 1);
147  _cADC_vs_TS_EChannel[histoname]->setAxisTitle("ADC QIE8", 2);
148  }
149 
150  //book global histos
151  ib.setCurrentFolder("Hcal/ZDCTask");
152 
153  _cShape = ib.bookProfile("Shape",
154  "Shape",
155  xAxisShape.nbins(),
156  xAxisShape.min(),
157  xAxisShape.max(),
158  yAxisShape.nbins(),
159  yAxisShape.min(),
160  yAxisShape.max());
161  _cShape->setAxisTitle("Timing", 1);
162  _cShape->setAxisTitle("fC QIE8", 2);
163 
164  _cADC = ib.book1D("ADC", "ADC", xAxisADC.nbins(), xAxisADC.min(), xAxisADC.max());
165  _cADC->getRootObject()->SetBit(BIT(hcaldqm::constants::BIT_OFFSET + hcaldqm::quantity::AxisType::fYAxis));
166  _cADC->setAxisTitle("ADC QIE8", 1);
167 
168  _cADC_vs_TS = ib.book2D("ADC_vs_TS",
169  "ADC_vs_TS",
170  xAxisShape.nbins(),
171  xAxisShape.min(),
172  xAxisShape.max(),
173  xAxisADC.nbins(),
174  xAxisADC.min(),
175  xAxisADC.max());
176  _cADC_vs_TS->getRootObject()->SetBit(BIT(hcaldqm::constants::BIT_OFFSET + hcaldqm::quantity::AxisType::fYAxis));
177  _cADC_vs_TS->setAxisTitle("Timing", 1);
178  _cADC_vs_TS->setAxisTitle("ADC QIE8", 2);
179 }
180 
181 /* virtual */ void ZDCTask::analyze(edm::Event const& e, edm::EventSetup const&) {
182  edm::Handle<ZDCDigiCollection> cqie10;
183  if (!e.getByToken(_tokQIE10, cqie10))
184  edm::LogError("Collection ZDCDigiCollection isn't available" + _tagQIE10.label() + " " + _tagQIE10.instance());
185 
186  for (uint32_t i = 0; i < cqie10->size(); i++) {
187  ZDCDataFrame frame = static_cast<ZDCDataFrame>((*cqie10)[i]);
188  HcalElectronicsId eid = frame.elecId();
189 
190  char histoname[300];
191  sprintf(histoname, "%d_%d_%d_%d", eid.fiberChanId(), eid.fiberIndex(), eid.spigot(), eid.dccid());
192 
193  // compute the signal, ped subracted
194  //double q = hcaldqm::utilities::sumQ_v10<ZDCDataFrame>(frame, constants::adc2fC[_ped], 0, frame.size()-1);
195 
196  // iterate thru all TS and fill
197  for (int j = 0; j < frame.size(); j++) {
198  _cShape_EChannel[histoname]->Fill(j, frame[j].nominal_fC());
199  _cShape->Fill(j, frame[j].nominal_fC());
200 
201  _cADC_EChannel[histoname]->Fill(frame[j].adc());
202  _cADC->Fill(frame[j].adc());
203 
204  _cADC_vs_TS_EChannel[histoname]->Fill(j, frame[j].adc());
205  _cADC_vs_TS->Fill(j, frame[j].adc());
206  }
207  }
208 }
209 
210 DEFINE_FWK_MODULE(ZDCTask);
HcalElectronicsId::setHTR
constexpr void setHTR(int crate, int slot, int tb)
Definition: HcalElectronicsId.h:72
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
mps_fire.i
i
Definition: mps_fire.py:338
ZDCTask::analyze
void analyze(edm::Event const &, edm::EventSetup const &) override
Definition: ZDCTask.cc:181
DQMStore::IBooker::bookProfile
MonitorElement * bookProfile(Args &&...args)
Definition: DQMStore.h:113
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:517
MonitorElement::getRootObject
TObject * getRootObject() const
Definition: MonitorElement.cc:1469
MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37
HcalElectronicsId::dccid
constexpr int dccid() const
get the (Hcal local) DCC id for VME, crate number for uTCA
Definition: HcalElectronicsId.h:98
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:16
DQMStore::IBooker::setCurrentFolder
void setCurrentFolder(std::string const &fullpath)
Definition: DQMStore.cc:268
ZDCTask::ZDCTask
ZDCTask(edm::ParameterSet const &)
Definition: ZDCTask.cc:7
DQMStore::IBooker::book1D
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:106
HcalElectronicsId::spigot
constexpr int spigot() const
get the spigot (input number on DCC, AMC card number for uTCA)
Definition: HcalElectronicsId.h:96
ecalMGPA::adc
constexpr int adc(sample_type sample)
get the ADC sample (12 bits)
Definition: EcalMGPASample.h:11
DetId
Definition: DetId.h:18
HcalElectronicsId::fiberIndex
constexpr int fiberIndex() const
get the fiber index. For VME 1-8 (which of eight fibers carried by a spigot), for uTCA fibers are zer...
Definition: HcalElectronicsId.h:83
ZDCTask::bookHistograms
void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
Definition: ZDCTask.cc:17
DQMStore::IBooker::book2D
MonitorElement * book2D(Args &&...args)
Definition: DQMStore.h:109
ZDCDataFrame::size
int size() const
total number of samples in the digi
Definition: ZDCDataFrame.h:26
edm::SortedCollection::size
size_type size() const
Definition: SortedCollection.h:236
ZDCDataFrame::elecId
const HcalElectronicsId & elecId() const
Definition: ZDCDataFrame.h:23
hcaldqm::constants::BIT_OFFSET
int const BIT_OFFSET
Definition: Constants.h:277
HcalElectronicsId::fiberChanId
constexpr int fiberChanId() const
get the fiber channel id (which of channels on a fiber)
Definition: HcalElectronicsId.h:81
MonitorElement::setAxisTitle
void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:973
HcalElectronicsId
Readout chain identification for Hcal.
Definition: HcalElectronicsId.h:32
edm::Run
Definition: Run.h:45
cuy.ib
ib
Definition: cuy.py:662
Generated for CMSSW Reference Manual by doxygen 1.8.11