|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
|
#include "DQM/CastorMonitor/interface/CastorLEDMonitor.h"
#include "DQMServices/Core/interface/DQMStore.h"
//***************************************************//
//********** CastorLEDMonitor ***********************//
//********** Author: Dmytro Volyanskyy ************//
//********** Date : 20.11.2008 (first version) ******//
//---------- last revision: 31.05.2011 (Panos Katsas)
//***************************************************//
//---- critical revision 26.06.2014 (Vladimir Popov)
//==================================================================//
CastorLEDMonitor::CastorLEDMonitor(const edm::ParameterSet &ps) {
fVerbosity = ps.getUntrackedParameter<int>("debug", 0);
if (fVerbosity > 0)
std::cout << "CastorLEDMonitor Constructor: " << this << std::endl;
subsystemname = ps.getUntrackedParameter<std::string>("subSystemFolder", "Castor");
ievt_ = 0;
}
CastorLEDMonitor::~CastorLEDMonitor() {}
void CastorLEDMonitor::bookHistograms(DQMStore::IBooker &ibooker, const edm::Run &iRun) {
char s[60];
ibooker.setCurrentFolder(subsystemname + "/CastorLEDMonitor");
sprintf(s, "CastorLEDqMap(cumulative)");
h2qMap = ibooker.book2D(s, s, 14, 0, 14, 16, 0, 16);
h2qMap->setOption("colz");
sprintf(s, "CastorLED_QmeanMap");
h2meanMap = ibooker.book2D(s, s, 14, 0, 14, 16, 0, 16);
h2meanMap->setAxisTitle("moduleZ", /* axis */ 1);
h2meanMap->setAxisTitle("sectorPhi", /* axis */ 2);
h2meanMap->setOption("colz");
ievt_ = 0;
return;
}
void CastorLEDMonitor::processEvent(const CastorDigiCollection &castorDigis, const CastorDbService &cond) {
if (fVerbosity > 0)
std::cout << "CastorLEDMonitor::processEvent (start)" << std::endl;
/* be implemented
edm::Handle<HcalTBTriggerData> trigger_data;
iEvent.getByToken(tok_tb_, trigger_data);
if(trigger_data.isValid())
if(trigger_data->triggerWord()==6) LEDevent=true;
*/
if (castorDigis.empty()) {
if (fVerbosity > 0)
std::cout << "CastorLEDMonitor::processEvent NO Castor Digis" << std::endl;
return;
}
for (CastorDigiCollection::const_iterator j = castorDigis.begin(); j != castorDigis.end(); j++) {
const CastorDataFrame digi = (const CastorDataFrame)(*j);
int module = digi.id().module() - 1;
int sector = digi.id().sector() - 1;
double qsum = 0.;
for (int i = 0; i < digi.size(); i++) {
int dig = digi.sample(i).adc() & 0x7f;
float ets = LedMonAdc2fc[dig] + 0.5;
// h2qts->Fill(i,ets);
qsum += ets;
}
// int ind = sector*14 + module;
// h2QvsPMT->Fill(ind,qsum);
h2qMap->Fill(module, sector, qsum);
} // end for(CastorDigiCollection::const_iterator j=castorDigis...
ievt_++;
if (ievt_ % 100 == 0) {
for (int mod = 1; mod <= 14; mod++)
for (int sec = 1; sec <= 16; sec++) {
double a = h2qMap->getBinContent(mod, sec);
h2meanMap->setBinContent(mod, sec, a / double(ievt_));
}
}
if (fVerbosity > 0)
std::cout << "CastorLEDMonitor::processEvent(end)" << std::endl;
return;
}
|