|
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
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
|
/*
* \file L1TCSCTPG.cc
*
* \author J. Berryhill
*
*/
#include "DQM/L1TMonitor/interface/L1TCSCTPG.h"
using namespace std;
using namespace edm;
L1TCSCTPG::L1TCSCTPG(const ParameterSet& ps)
: csctpgSource_(ps.getParameter<InputTag>("csctpgSource")),
csctpgSource_token_(consumes<CSCCorrelatedLCTDigiCollection>(ps.getParameter<InputTag>("csctpgSource"))) {
// verbosity switch
verbose_ = ps.getUntrackedParameter<bool>("verbose", false);
if (verbose_)
cout << "L1TCSCTPG: constructor...." << endl;
outputFile_ = ps.getUntrackedParameter<string>("outputFile", "");
if (!outputFile_.empty()) {
cout << "L1T Monitoring histograms will be saved to " << outputFile_.c_str() << endl;
}
bool disable = ps.getUntrackedParameter<bool>("disableROOToutput", false);
if (disable) {
outputFile_ = "";
}
}
L1TCSCTPG::~L1TCSCTPG() {}
void L1TCSCTPG::bookHistograms(DQMStore::IBooker& ibooker, edm::Run const& iRun, edm::EventSetup const& iSetup) {
nev_ = 0;
ibooker.setCurrentFolder("L1T/L1TCSCTPG");
csctpgpattern = ibooker.book1D("CSC TPG hit pattern", "CSC TPG hit pattern", 8, -0.5, 7.5);
csctpgquality = ibooker.book1D("CSC TPG quality", "CSC TPG quality", 16, 0.5, 16.5);
csctpgwg = ibooker.book1D("CSC TPG wire group", "CSC TPG wire group", 116, -0.5, 115.5);
csctpgstrip = ibooker.book1D("CSC TPG strip", "CSC TPG strip", 160, -0.5, 159.5);
csctpgstriptype = ibooker.book1D("CSC TPG strip type", "CSC TPG strip type", 2, 0.5, 1.5);
csctpgbend = ibooker.book1D("CSC TPG bend", "CSC TPG bend", 3, 0.5, 2.5);
csctpgbx = ibooker.book1D("CSC TPG bx", "CSC TPG bx", 20, -0.5, 19.5);
}
void L1TCSCTPG::analyze(const Event& e, const EventSetup& c) {
nev_++;
if (verbose_)
cout << "L1TCSCTPG: analyze...." << endl;
Handle<CSCCorrelatedLCTDigiCollection> pCSCTPGcorrlcts;
e.getByToken(csctpgSource_token_, pCSCTPGcorrlcts);
if (!pCSCTPGcorrlcts.isValid()) {
edm::LogInfo("DataNotFound") << "can't find CSCCorrelatedLCTDigiCollection with label " << csctpgSource_.label();
return;
}
for (CSCCorrelatedLCTDigiCollection::DigiRangeIterator cscItr1 = pCSCTPGcorrlcts->begin();
cscItr1 != pCSCTPGcorrlcts->end();
cscItr1++) {
CSCCorrelatedLCTDigiCollection::Range range1 = pCSCTPGcorrlcts->get((*cscItr1).first);
for (CSCCorrelatedLCTDigiCollection::const_iterator lctItr1 = range1.first; lctItr1 != range1.second; lctItr1++) {
csctpgpattern->Fill(lctItr1->getCLCTPattern());
if (verbose_) {
std::cout << "CSC TPG CLCT pattern " << lctItr1->getCLCTPattern() << std::endl;
}
csctpgquality->Fill(lctItr1->getQuality());
if (verbose_) {
std::cout << "CSC LCT quality " << lctItr1->getQuality() << std::endl;
}
csctpgwg->Fill(lctItr1->getKeyWG());
if (verbose_) {
std::cout << "CSC LCT wire group " << lctItr1->getKeyWG() << std::endl;
}
csctpgstrip->Fill(lctItr1->getStrip());
if (verbose_) {
std::cout << "CSC LCT strip " << lctItr1->getStrip() << std::endl;
}
csctpgstriptype->Fill(lctItr1->getStripType());
if (verbose_) {
std::cout << "CSC LCT strip type" << lctItr1->getStripType() << std::endl;
}
csctpgbend->Fill(lctItr1->getBend());
if (verbose_) {
std::cout << "CSC LCT bend " << lctItr1->getBend() << std::endl;
}
csctpgbx->Fill(lctItr1->getBX());
if (verbose_) {
std::cout << "CSC LCT bx " << lctItr1->getBX() << std::endl;
}
}
}
}
|