CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
HcalPedestalsAnalysis.cc
Go to the documentation of this file.
1 // Original Author: Steven Won
2 // Created: Fri May 2 15:34:43 CEST 2008
3 // Written to replace the combination of HcalPedestalAnalyzer and HcalPedestalAnalysis
4 // This code runs 1000x faster and produces all outputs from a single run
5 // (ADC, fC in .txt plus an .xml file)
6 //
7 // $Id: HcalPedestalsAnalysis.cc,v 1.25 2012/11/13 03:30:20 dlange Exp $
8 
9 #include <memory>
10 #include "CalibCalorimetry/HcalStandardModules/interface/HcalPedestalsAnalysis.h"
11 #include "Geometry/CaloTopology/interface/HcalTopology.h"
12 
13 HcalPedestalsAnalysis::HcalPedestalsAnalysis(const edm::ParameterSet& ps) :
14  hbheDigiCollectionTag_(ps.getParameter<edm::InputTag>("hbheDigiCollectionTag")),
15  hoDigiCollectionTag_(ps.getParameter<edm::InputTag>("hoDigiCollectionTag")),
16  hfDigiCollectionTag_(ps.getParameter<edm::InputTag>("hfDigiCollectionTag")) {
17 
18  std::cout << "Code version 10.6\n";
19  hiSaveFlag = ps.getUntrackedParameter<bool>("hiSaveFlag", false);
20  dumpXML = ps.getUntrackedParameter<bool>("dumpXML", true);
21  verboseflag = ps.getUntrackedParameter<bool>("verbose", false);
22  firstTS = ps.getUntrackedParameter<int>("firstTS", 0);
23  lastTS = ps.getUntrackedParameter<int>("lastTS", 9);
24  firsttime = true;
25  ievt = 0;
26 
27  rawPedsItem = 0;
28  rawWidthsItem = 0;
29  rawPedsItemfc = 0;
30  rawWidthsItemfc = 0;
31 }
32 
33 
34 HcalPedestalsAnalysis::~HcalPedestalsAnalysis()
35 {}
36 
37 void HcalPedestalsAnalysis::endJob()
38 {
39 // std::cout << "NEvents " << ievt << std::endl;
40  if(ievt < 1000) return;
41  std::string ievt_string;
42  std::stringstream tempstringout;
43  tempstringout << ievt;
44  ievt_string = tempstringout.str();
45  pedsADCfilename += ievt_string;
46  pedsADCfilename += ".txt";
47  //Calculate pedestal constants
48  std::cout << "Calculating Pedestal constants...\n";
49  std::vector<NewPedBunch>::iterator bunch_it;
50  for(bunch_it=Bunches.begin(); bunch_it != Bunches.end(); bunch_it++)
51  {
52  if(bunch_it->usedflag){
53 
54  if(verboseflag) std::cout << "Analyzing channel " << bunch_it->detid << std::endl;
55  //pedestal constant is the mean
56  if(bunch_it->num[0][0]!=0) bunch_it->cap[0] /= bunch_it->num[0][0];
57  if(bunch_it->num[1][1]!=0) bunch_it->cap[1] /= bunch_it->num[1][1];
58  if(bunch_it->num[2][2]!=0) bunch_it->cap[2] /= bunch_it->num[2][2];
59  if(bunch_it->num[3][3]!=0) bunch_it->cap[3] /= bunch_it->num[3][3];
60  if(bunch_it->num[0][0]!=0) bunch_it->capfc[0] /= bunch_it->num[0][0];
61  if(bunch_it->num[1][1]!=0) bunch_it->capfc[1] /= bunch_it->num[1][1];
62  if(bunch_it->num[2][2]!=0) bunch_it->capfc[2] /= bunch_it->num[2][2];
63  if(bunch_it->num[3][3]!=0) bunch_it->capfc[3] /= bunch_it->num[3][3];
64  bunch_it->sig[0][0] = (bunch_it->prod[0][0]/bunch_it->num[0][0])-(bunch_it->cap[0]*bunch_it->cap[0]);
65  bunch_it->sig[0][1] = (bunch_it->prod[0][1]/bunch_it->num[0][1])-(bunch_it->cap[0]*bunch_it->cap[1]);
66  bunch_it->sig[0][2] = (bunch_it->prod[0][2]/bunch_it->num[0][2])-(bunch_it->cap[0]*bunch_it->cap[2]);
67  bunch_it->sig[0][3] = (bunch_it->prod[0][3]/bunch_it->num[0][3])-(bunch_it->cap[0]*bunch_it->cap[3]);
68  bunch_it->sig[1][0] = (bunch_it->prod[1][0]/bunch_it->num[1][0])-(bunch_it->cap[1]*bunch_it->cap[0]);
69  bunch_it->sig[1][1] = (bunch_it->prod[1][1]/bunch_it->num[1][1])-(bunch_it->cap[1]*bunch_it->cap[1]);
70  bunch_it->sig[1][2] = (bunch_it->prod[1][2]/bunch_it->num[1][2])-(bunch_it->cap[1]*bunch_it->cap[2]);
71  bunch_it->sig[1][3] = (bunch_it->prod[1][3]/bunch_it->num[1][3])-(bunch_it->cap[1]*bunch_it->cap[3]);
72  bunch_it->sig[2][0] = (bunch_it->prod[2][0]/bunch_it->num[2][0])-(bunch_it->cap[2]*bunch_it->cap[0]);
73  bunch_it->sig[2][1] = (bunch_it->prod[2][1]/bunch_it->num[2][1])-(bunch_it->cap[2]*bunch_it->cap[1]);
74  bunch_it->sig[2][2] = (bunch_it->prod[2][2]/bunch_it->num[2][2])-(bunch_it->cap[2]*bunch_it->cap[2]);
75  bunch_it->sig[2][3] = (bunch_it->prod[2][3]/bunch_it->num[2][3])-(bunch_it->cap[2]*bunch_it->cap[3]);
76  bunch_it->sig[3][0] = (bunch_it->prod[3][0]/bunch_it->num[3][0])-(bunch_it->cap[3]*bunch_it->cap[0]);
77  bunch_it->sig[3][1] = (bunch_it->prod[3][1]/bunch_it->num[3][1])-(bunch_it->cap[3]*bunch_it->cap[1]);
78  bunch_it->sig[3][2] = (bunch_it->prod[3][2]/bunch_it->num[3][2])-(bunch_it->cap[3]*bunch_it->cap[2]);
79  bunch_it->sig[3][3] = (bunch_it->prod[3][3]/bunch_it->num[3][3])-(bunch_it->cap[3]*bunch_it->cap[3]);
80 
81  bunch_it->sigfc[0][0] = (bunch_it->prodfc[0][0]/bunch_it->num[0][0])-(bunch_it->capfc[0]*bunch_it->capfc[0]);
82  bunch_it->sigfc[0][1] = (bunch_it->prodfc[0][1]/bunch_it->num[0][1])-(bunch_it->capfc[0]*bunch_it->capfc[1]);
83  bunch_it->sigfc[0][2] = (bunch_it->prodfc[0][2]/bunch_it->num[0][2])-(bunch_it->capfc[0]*bunch_it->capfc[2]);
84  bunch_it->sigfc[0][3] = (bunch_it->prodfc[0][3]/bunch_it->num[0][3])-(bunch_it->capfc[0]*bunch_it->capfc[3]);
85  bunch_it->sigfc[1][0] = (bunch_it->prodfc[1][0]/bunch_it->num[1][0])-(bunch_it->capfc[1]*bunch_it->capfc[0]);
86  bunch_it->sigfc[1][1] = (bunch_it->prodfc[1][1]/bunch_it->num[1][1])-(bunch_it->capfc[1]*bunch_it->capfc[1]);
87  bunch_it->sigfc[1][2] = (bunch_it->prodfc[1][2]/bunch_it->num[1][2])-(bunch_it->capfc[1]*bunch_it->capfc[2]);
88  bunch_it->sigfc[1][3] = (bunch_it->prodfc[1][3]/bunch_it->num[1][3])-(bunch_it->capfc[1]*bunch_it->capfc[3]);
89  bunch_it->sigfc[2][0] = (bunch_it->prodfc[2][0]/bunch_it->num[2][0])-(bunch_it->capfc[2]*bunch_it->capfc[0]);
90  bunch_it->sigfc[2][1] = (bunch_it->prodfc[2][1]/bunch_it->num[2][1])-(bunch_it->capfc[2]*bunch_it->capfc[1]);
91  bunch_it->sigfc[2][2] = (bunch_it->prodfc[2][2]/bunch_it->num[2][2])-(bunch_it->capfc[2]*bunch_it->capfc[2]);
92  bunch_it->sigfc[2][3] = (bunch_it->prodfc[2][3]/bunch_it->num[2][3])-(bunch_it->capfc[2]*bunch_it->capfc[3]);
93  bunch_it->sigfc[3][0] = (bunch_it->prodfc[3][0]/bunch_it->num[3][0])-(bunch_it->capfc[3]*bunch_it->capfc[0]);
94  bunch_it->sigfc[3][1] = (bunch_it->prodfc[3][1]/bunch_it->num[3][1])-(bunch_it->capfc[3]*bunch_it->capfc[1]);
95  bunch_it->sigfc[3][2] = (bunch_it->prodfc[3][2]/bunch_it->num[3][2])-(bunch_it->capfc[3]*bunch_it->capfc[2]);
96  bunch_it->sigfc[3][3] = (bunch_it->prodfc[3][3]/bunch_it->num[3][3])-(bunch_it->capfc[3]*bunch_it->capfc[3]);
97 
98  if(bunch_it->detid.subdet() == 1){
99  for(int i = 0; i != 4; i++){
100  HBMeans->Fill(bunch_it->cap[i]);
101  HBWidths->Fill(bunch_it->sig[i][i]);
102  }
103  }
104  if(bunch_it->detid.subdet() == 2){
105  for(int i = 0; i != 4; i++){
106  HEMeans->Fill(bunch_it->cap[i]);
107  HEWidths->Fill(bunch_it->sig[i][i]);
108  }
109  }
110  if(bunch_it->detid.subdet() == 3){
111  for(int i = 0; i != 4; i++){
112  HOMeans->Fill(bunch_it->cap[i]);
113  HOWidths->Fill(bunch_it->sig[i][i]);
114  }
115  }
116  if(bunch_it->detid.subdet() == 4){
117  for(int i = 0; i != 4; i++){
118  HFMeans->Fill(bunch_it->cap[i]);
119  HFWidths->Fill(bunch_it->sig[i][i]);
120  }
121  }
122 
123  const HcalPedestal item(bunch_it->detid, bunch_it->cap[0], bunch_it->cap[1], bunch_it->cap[2], bunch_it->cap[3],
124  bunch_it->sig[0][0], bunch_it->sig[1][1], bunch_it->sig[2][2], bunch_it->sig[3][3]);
125  rawPedsItem->addValues(item);
126  HcalPedestalWidth widthsp(bunch_it->detid);
127  widthsp.setSigma(0,0,bunch_it->sig[0][0]);
128  widthsp.setSigma(0,1,bunch_it->sig[0][1]);
129  widthsp.setSigma(0,2,bunch_it->sig[0][2]);
130  widthsp.setSigma(0,3,bunch_it->sig[0][3]);
131  widthsp.setSigma(1,1,bunch_it->sig[1][1]);
132  widthsp.setSigma(1,2,bunch_it->sig[1][2]);
133  widthsp.setSigma(1,3,bunch_it->sig[1][3]);
134  widthsp.setSigma(2,2,bunch_it->sig[2][2]);
135  widthsp.setSigma(2,3,bunch_it->sig[2][3]);
136  widthsp.setSigma(3,3,bunch_it->sig[3][3]);
137  rawWidthsItem->addValues(widthsp);
138 
139  const HcalPedestal itemfc(bunch_it->detid, bunch_it->capfc[0], bunch_it->capfc[1], bunch_it->capfc[2], bunch_it->capfc[3],
140  bunch_it->sigfc[0][0], bunch_it->sigfc[1][1], bunch_it->sigfc[2][2], bunch_it->sigfc[3][3]);
141  rawPedsItemfc->addValues(itemfc);
142  HcalPedestalWidth widthspfc(bunch_it->detid);
143  widthspfc.setSigma(0,0,bunch_it->sigfc[0][0]);
144  widthspfc.setSigma(0,1,bunch_it->sigfc[0][1]);
145  widthspfc.setSigma(0,2,bunch_it->sigfc[0][2]);
146  widthspfc.setSigma(0,3,bunch_it->sigfc[0][3]);
147  widthspfc.setSigma(1,1,bunch_it->sigfc[1][1]);
148  widthspfc.setSigma(1,2,bunch_it->sigfc[1][2]);
149  widthspfc.setSigma(1,3,bunch_it->sigfc[1][3]);
150  widthspfc.setSigma(2,2,bunch_it->sigfc[2][2]);
151  widthspfc.setSigma(2,3,bunch_it->sigfc[2][3]);
152  widthspfc.setSigma(3,3,bunch_it->sigfc[3][3]);
153  rawWidthsItemfc->addValues(widthspfc);
154  }
155  }
156 
157  // dump the resulting list of pedestals into a file
158  std::ofstream outStream1(pedsADCfilename.c_str());
159  HcalDbASCIIIO::dumpObject (outStream1, (*rawPedsItem) );
160 // std::ofstream outStream2(widthsADCfilename.c_str());
161 // HcalDbASCIIIO::dumpObject (outStream2, (*rawWidthsItem) );
162 //
163 // std::ofstream outStream3(pedsfCfilename.c_str());
164 // HcalDbASCIIIO::dumpObject (outStream3, (*rawPedsItemfc) );
165 // std::ofstream outStream4(widthsfCfilename.c_str());
166 // HcalDbASCIIIO::dumpObject (outStream4, (*rawWidthsItemfc) );
167 //
168 // if(dumpXML){
169 // std::ofstream outStream5(XMLfilename.c_str());
170 // HcalDbXml::dumpObject (outStream5, runnum, 0, 2147483647, XMLtag, 1, (*rawPedsItem), (*rawWidthsItem));
171 // }
172 //
173  if(hiSaveFlag){
174  theFile->Write();
175  }else{
176  theFile->cd();
177  theFile->cd("HB");
178  HBMeans->Write();
179  HBWidths->Write();
180  theFile->cd();
181  theFile->cd("HF");
182  HFMeans->Write();
183  HFWidths->Write();
184  theFile->cd();
185  theFile->cd("HE");
186  HEMeans->Write();
187  HEWidths->Write();
188  theFile->cd();
189  theFile->cd("HO");
190  HOMeans->Write();
191  HOWidths->Write();
192  }
193 
194  std::cout << "Writing ROOT file... ";
195  theFile->Close();
196  std::cout << "ROOT file closed.\n";
197 
198  delete rawPedsItem;
199  delete rawWidthsItem;
200  delete rawPedsItemfc;
201  delete rawWidthsItemfc;
202 }
203 
204 // ------------ method called to for each event ------------
205 void
206 HcalPedestalsAnalysis::analyze(const edm::Event& e, const edm::EventSetup& iSetup)
207 {
208  using namespace edm;
209  using namespace std;
210 
211  edm::Handle<HBHEDigiCollection> hbhe; e.getByLabel(hbheDigiCollectionTag_, hbhe);
212  edm::Handle<HODigiCollection> ho; e.getByLabel(hoDigiCollectionTag_, ho);
213  edm::Handle<HFDigiCollection> hf; e.getByLabel(hfDigiCollectionTag_, hf);
214  edm::ESHandle<HcalDbService> conditions;
215  iSetup.get<HcalDbRecord>().get(conditions);
216 
217 
218  if(firsttime)
219  {
220 
221  edm::ESHandle<HcalTopology> topology;
222  iSetup.get<IdealGeometryRecord>().get( topology );
223  theTopology=new HcalTopology(*topology);
224 
225  rawPedsItem = new HcalPedestals(theTopology,true);
226  rawWidthsItem = new HcalPedestalWidths(theTopology,true);
227  rawPedsItemfc = new HcalPedestals(theTopology,false);
228  rawWidthsItemfc = new HcalPedestalWidths(theTopology,false);
229 
230  runnum = e.id().run();
231  std::string runnum_string;
232  std::stringstream tempstringout;
233  tempstringout << runnum;
234  runnum_string = tempstringout.str();
235  ROOTfilename = runnum_string + "-peds_ADC.root";
236  pedsADCfilename = runnum_string + "-peds_ADC_";
237  pedsfCfilename = runnum_string + "-peds_fC.txt";
238  widthsADCfilename = runnum_string + "-widths_ADC.txt";
239  widthsfCfilename = runnum_string + "-widths_fC.txt";
240  XMLfilename = runnum_string + "-peds_ADC_complete.xml";
241  XMLtag = "Hcal_pedestals_" + runnum_string;
242 
243  theFile = new TFile(ROOTfilename.c_str(), "RECREATE");
244  theFile->cd();
245  // Create sub-directories
246  theFile->mkdir("HB");
247  theFile->mkdir("HE");
248  theFile->mkdir("HF");
249  theFile->mkdir("HO");
250  theFile->cd();
251 
252  HBMeans = new TH1F("All Ped Means HB","All Ped Means HB", 100, 0, 9);
253  HBWidths = new TH1F("All Ped Widths HB","All Ped Widths HB", 100, 0, 3);
254  HEMeans = new TH1F("All Ped Means HE","All Ped Means HE", 100, 0, 9);
255  HEWidths = new TH1F("All Ped Widths HE","All Ped Widths HE", 100, 0, 3);
256  HFMeans = new TH1F("All Ped Means HF","All Ped Means HF", 100, 0, 9);
257  HFWidths = new TH1F("All Ped Widths HF","All Ped Widths HF", 100, 0, 3);
258  HOMeans = new TH1F("All Ped Means HO","All Ped Means HO", 100, 0, 9);
259  HOWidths = new TH1F("All Ped Widths HO","All Ped Widths HO", 100, 0, 3);
260 
261  edm::ESHandle<HcalElectronicsMap> refEMap;
262  iSetup.get<HcalElectronicsMapRcd>().get(refEMap);
263  const HcalElectronicsMap* myRefEMap = refEMap.product();
264  std::vector<HcalGenericDetId> listEMap = myRefEMap->allPrecisionId();
265  for (std::vector<HcalGenericDetId>::const_iterator it = listEMap.begin(); it != listEMap.end(); it++)
266  {
267  HcalGenericDetId mygenid(it->rawId());
268  if(mygenid.isHcalDetId())
269  {
270  NewPedBunch a;
271  HcalDetId chanid(mygenid.rawId());
272  a.detid = chanid;
273  a.usedflag = false;
274  for(int i = 0; i != 4; i++)
275  {
276  a.cap[i] = 0;
277  a.capfc[i] = 0;
278  for(int j = 0; j != 4; j++)
279  {
280  a.sig[i][j] = 0;
281  a.sigfc[i][j] = 0;
282  a.prod[i][j] = 0;
283  a.prodfc[i][j] = 0;
284  a.num[i][j] = 0;
285  }
286  }
287  Bunches.push_back(a);
288  }
289  }
290  firsttime = false;
291  }
292 
293  std::vector<NewPedBunch>::iterator bunch_it;
294 
295  for(HBHEDigiCollection::const_iterator j = hbhe->begin(); j != hbhe->end(); j++)
296  {
297  const HBHEDataFrame digi = (const HBHEDataFrame)(*j);
298  for(bunch_it = Bunches.begin(); bunch_it != Bunches.end(); bunch_it++)
299  if(bunch_it->detid.rawId() == digi.id().rawId()) break;
300  bunch_it->usedflag = true;
301  for(int ts = firstTS; ts != lastTS+1; ts++)
302  {
303  const HcalQIECoder* coder = conditions->getHcalCoder(digi.id().rawId());
304  const HcalQIEShape* shape = conditions->getHcalShape(coder);
305  bunch_it->num[digi.sample(ts).capid()][digi.sample(ts).capid()] += 1;
306  bunch_it->cap[digi.sample(ts).capid()] += digi.sample(ts).adc();
307  double charge1 = coder->charge(*shape, digi.sample(ts).adc(), digi.sample(ts).capid());
308  bunch_it->capfc[digi.sample(ts).capid()] += charge1;
309  bunch_it->prod[digi.sample(ts).capid()][digi.sample(ts).capid()] += (digi.sample(ts).adc() * digi.sample(ts).adc());
310  bunch_it->prodfc[digi.sample(ts).capid()][digi.sample(ts).capid()] += charge1 * charge1;
311  if((ts+1 < digi.size()) && (ts+1 < lastTS)){
312  bunch_it->prod[digi.sample(ts).capid()][digi.sample(ts+1).capid()] += digi.sample(ts).adc()*digi.sample(ts+1).adc();
313  double charge2 = coder->charge(*shape, digi.sample(ts+1).adc(), digi.sample(ts+1).capid());
314  bunch_it->prodfc[digi.sample(ts).capid()][digi.sample(ts+1).capid()] += charge1*charge2;
315  bunch_it->num[digi.sample(ts).capid()][digi.sample(ts+1).capid()] += 1;
316  }
317  if((ts+2 < digi.size()) && (ts+2 < lastTS)){
318  bunch_it->prod[digi.sample(ts).capid()][digi.sample(ts+2).capid()] += digi.sample(ts).adc()*digi.sample(ts+2).adc();
319  double charge2 = coder->charge(*shape, digi.sample(ts+2).adc(), digi.sample(ts+2).capid());
320  bunch_it->prodfc[digi.sample(ts).capid()][digi.sample(ts+2).capid()] += charge1*charge2;
321  bunch_it->num[digi.sample(ts).capid()][digi.sample(ts+2).capid()] += 1;
322  }
323  if((ts+3 < digi.size()) && (ts+3 < lastTS)){
324  bunch_it->prod[digi.sample(ts).capid()][digi.sample(ts+3).capid()] += digi.sample(ts).adc()*digi.sample(ts+3).adc();
325  double charge2 = coder->charge(*shape, digi.sample(ts+3).adc(), digi.sample(ts+3).capid());
326  bunch_it->prodfc[digi.sample(ts).capid()][digi.sample(ts+3).capid()] += charge1*charge2;
327  bunch_it->num[digi.sample(ts).capid()][digi.sample(ts+3).capid()] += 1;
328  }
329  }
330  }
331 
332  for(HODigiCollection::const_iterator j = ho->begin(); j != ho->end(); j++)
333  {
334  const HODataFrame digi = (const HODataFrame)(*j);
335  for(bunch_it = Bunches.begin(); bunch_it != Bunches.end(); bunch_it++)
336  if(bunch_it->detid.rawId() == digi.id().rawId()) break;
337  bunch_it->usedflag = true;
338  for(int ts = firstTS; ts <= lastTS; ts++)
339  {
340  const HcalQIECoder* coder = conditions->getHcalCoder(digi.id().rawId());
341  const HcalQIEShape* shape = conditions->getHcalShape(coder);
342  bunch_it->num[digi.sample(ts).capid()][digi.sample(ts).capid()] += 1;
343  bunch_it->cap[digi.sample(ts).capid()] += digi.sample(ts).adc();
344  double charge1 = coder->charge(*shape, digi.sample(ts).adc(), digi.sample(ts).capid());
345  bunch_it->capfc[digi.sample(ts).capid()] += charge1;
346  bunch_it->prod[digi.sample(ts).capid()][digi.sample(ts).capid()] += (digi.sample(ts).adc() * digi.sample(ts).adc());
347  bunch_it->prodfc[digi.sample(ts).capid()][digi.sample(ts).capid()] += charge1 * charge1;
348  if((ts+1 < digi.size()) && (ts+1 < lastTS)){
349  bunch_it->prod[digi.sample(ts).capid()][digi.sample(ts+1).capid()] += digi.sample(ts).adc()*digi.sample(ts+1).adc();
350  double charge2 = coder->charge(*shape, digi.sample(ts+1).adc(), digi.sample(ts+1).capid());
351  bunch_it->prodfc[digi.sample(ts).capid()][digi.sample(ts+1).capid()] += charge1*charge2;
352  bunch_it->num[digi.sample(ts).capid()][digi.sample(ts+1).capid()] += 1;
353  }
354  if((ts+2 < digi.size()) && (ts+2 < lastTS)){
355  bunch_it->prod[digi.sample(ts).capid()][digi.sample(ts+2).capid()] += digi.sample(ts).adc()*digi.sample(ts+2).adc();
356  double charge2 = coder->charge(*shape, digi.sample(ts+2).adc(), digi.sample(ts+2).capid());
357  bunch_it->prodfc[digi.sample(ts).capid()][digi.sample(ts+2).capid()] += charge1*charge2;
358  bunch_it->num[digi.sample(ts).capid()][digi.sample(ts+2).capid()] += 1;
359  }
360  if((ts+3 < digi.size()) && (ts+3 < lastTS)){
361  bunch_it->prod[digi.sample(ts).capid()][digi.sample(ts+3).capid()] += digi.sample(ts).adc()*digi.sample(ts+3).adc();
362  double charge2 = coder->charge(*shape, digi.sample(ts+3).adc(), digi.sample(ts+3).capid());
363  bunch_it->prodfc[digi.sample(ts).capid()][digi.sample(ts+3).capid()] += charge1*charge2;
364  bunch_it->num[digi.sample(ts).capid()][digi.sample(ts+3).capid()] += 1;
365  }
366  }
367  }
368 
369  for(HFDigiCollection::const_iterator j = hf->begin(); j != hf->end(); j++)
370  {
371  const HFDataFrame digi = (const HFDataFrame)(*j);
372  for(bunch_it = Bunches.begin(); bunch_it != Bunches.end(); bunch_it++)
373  if(bunch_it->detid.rawId() == digi.id().rawId()) break;
374  bunch_it->usedflag = true;
375  for(int ts = firstTS; ts <= lastTS; ts++)
376  {
377  const HcalQIECoder* coder = conditions->getHcalCoder(digi.id().rawId());
378  const HcalQIEShape* shape = conditions->getHcalShape(coder);
379  bunch_it->num[digi.sample(ts).capid()][digi.sample(ts).capid()] += 1;
380  bunch_it->cap[digi.sample(ts).capid()] += digi.sample(ts).adc();
381  double charge1 = coder->charge(*shape, digi.sample(ts).adc(), digi.sample(ts).capid());
382  bunch_it->capfc[digi.sample(ts).capid()] += charge1;
383  bunch_it->prod[digi.sample(ts).capid()][digi.sample(ts).capid()] += (digi.sample(ts).adc() * digi.sample(ts).adc());
384  bunch_it->prodfc[digi.sample(ts).capid()][digi.sample(ts).capid()] += charge1 * charge1;
385  if((ts+1 < digi.size()) && (ts+1 < lastTS)){
386  bunch_it->prod[digi.sample(ts).capid()][digi.sample(ts+1).capid()] += digi.sample(ts).adc()*digi.sample(ts+1).adc();
387  double charge2 = coder->charge(*shape, digi.sample(ts+1).adc(), digi.sample(ts+1).capid());
388  bunch_it->prodfc[digi.sample(ts).capid()][digi.sample(ts+1).capid()] += charge1*charge2;
389  bunch_it->num[digi.sample(ts).capid()][digi.sample(ts+1).capid()] += 1;
390  }
391  if((ts+2 < digi.size()) && (ts+2 < lastTS)){
392  bunch_it->prod[digi.sample(ts).capid()][digi.sample(ts+2).capid()] += digi.sample(ts).adc()*digi.sample(ts+2).adc();
393  double charge2 = coder->charge(*shape, digi.sample(ts+2).adc(), digi.sample(ts+2).capid());
394  bunch_it->prodfc[digi.sample(ts).capid()][digi.sample(ts+2).capid()] += charge1*charge2;
395  bunch_it->num[digi.sample(ts).capid()][digi.sample(ts+2).capid()] += 1;
396  }
397  if((ts+3 < digi.size()) && (ts+3 < lastTS)){
398  bunch_it->prod[digi.sample(ts).capid()][digi.sample(ts+3).capid()] += digi.sample(ts).adc()*digi.sample(ts+3).adc();
399  double charge2 = coder->charge(*shape, digi.sample(ts+3).adc(), digi.sample(ts+3).capid());
400  bunch_it->prodfc[digi.sample(ts).capid()][digi.sample(ts+3).capid()] += charge1*charge2;
401  bunch_it->num[digi.sample(ts).capid()][digi.sample(ts+3).capid()] += 1;
402  }
403  }
404  }
405 
406 //this is the last brace
407 ievt++;
408 }
409 
410 //define this as a plug-in
411 DEFINE_FWK_MODULE(HcalPedestalsAnalysis);
edm::EventID::run
RunNumber_t run() const
Definition: EventID.h:42
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
i
int i
Definition: DBlmapReader.cc:9
HcalPedestalsAnalysis::hbheDigiCollectionTag_
edm::InputTag hbheDigiCollectionTag_
Definition: HcalPedestalsAnalysis.h:120
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:17
HcalQIESample::adc
int adc() const
get the ADC sample
Definition: HcalQIESample.h:24
HBHEDataFrame::size
int size() const
total number of samples in the digi
Definition: HBHEDataFrame.h:26
edm::SortedCollection::const_iterator
std::vector< T >::const_iterator const_iterator
Definition: SortedCollection.h:80
HODataFrame::sample
const HcalQIESample & sample(int i) const
access a sample
Definition: HODataFrame.h:40
HcalPedestalsAnalysis::analyze
void analyze(const edm::Event &event, const edm::EventSetup &eventSetup)
Definition: HcalPedestalsAnalysis.cc:206
HcalPedestalsAnalysis::rawWidthsItem
HcalPedestalWidths * rawWidthsItem
Definition: HcalPedestalsAnalysis.h:115
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
HODataFrame::id
const HcalDetId & id() const
Definition: HODataFrame.h:23
DetId::rawId
uint32_t rawId() const
get the raw id
Definition: DetId.h:45
HODataFrame::size
int size() const
total number of samples in the digi
Definition: HODataFrame.h:27
HFDataFrame::sample
const HcalQIESample & sample(int i) const
access a sample
Definition: HFDataFrame.h:39
j
int j
Definition: DBlmapReader.cc:9
HcalElectronicsMap::allPrecisionId
std::vector< HcalGenericDetId > allPrecisionId() const
Definition: HcalElectronicsMap.cc:137
edm::Event::getByLabel
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:361
NewPedBunch::detid
HcalCastorDetId detid
Definition: CastorPedestalsAnalysis.h:62
HFDataFrame::size
int size() const
total number of samples in the digi
Definition: HFDataFrame.h:26
HcalQIESample::capid
int capid() const
get the Capacitor id
Definition: HcalQIESample.h:28
HBHEDataFrame::sample
const HcalQIESample & sample(int i) const
access a sample
Definition: HBHEDataFrame.h:39
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:62
edm::EventBase::id
edm::EventID id() const
Definition: EventBase.h:56
HcalPedestalsAnalysis::HcalPedestalsAnalysis
HcalPedestalsAnalysis(const edm::ParameterSet &ps)
Definition: HcalPedestalsAnalysis.cc:13
HcalDbASCIIIO::dumpObject
bool dumpObject(std::ostream &fOutput, const HcalPedestals &fObject)
a
double a
Definition: hdecay.h:121
HBHEDataFrame::id
const HcalDetId & id() const
Definition: HBHEDataFrame.h:22
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
HcalPedestalsAnalysis::rawWidthsItemfc
HcalPedestalWidths * rawWidthsItemfc
Definition: HcalPedestalsAnalysis.h:117
HFDataFrame::id
const HcalDetId & id() const
Definition: HFDataFrame.h:22
HcalCondObjectContainer::addValues
bool addValues(const Item &myItem)
Definition: HcalCondObjectContainer.h:182
HcalPedestalWidth::setSigma
void setSigma(int fCapId1, int fCapId2, float fSigma)
Definition: HcalPedestalWidth.cc:38
HcalPedestalsAnalysis::Bunches
std::vector< NewPedBunch > Bunches
Definition: HcalPedestalsAnalysis.h:86
HcalQIECoder::charge
float charge(const HcalQIEShape &fShape, unsigned fAdc, unsigned fCapId) const
ADC [0..127] + capid [0..3] -&gt; fC conversion.
Definition: HcalQIECoder.cc:22
Generated for CMSSW Reference Manual by doxygen 1.8.5