|
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
103
104
105
106
107
108
109
110
111
|
/*----------------------------------------------------------------------
Toy EDAnalyzer for testing purposes only.
----------------------------------------------------------------------*/
#include <stdexcept>
#include <iostream>
#include <fstream>
#include <map>
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "CoralBase/TimeStamp.h"
#include "CondTools/DT/test/stubs/DTHVDump.h"
#include "CondFormats/DTObjects/interface/DTHVStatus.h"
#include "CondFormats/DataRecord/interface/DTHVStatusRcd.h"
#include "DataFormats/MuonDetId/interface/DTWireId.h"
#include "DataFormats/MuonDetId/interface/DTChamberId.h"
namespace edmtest {
DTHVDump::DTHVDump(edm::ParameterSet const& p) : dthvstatusToken_(esConsumes()) {
// parameters to setup
}
DTHVDump::DTHVDump(int i) : dthvstatusToken_(esConsumes()) {}
void DTHVDump::analyze(const edm::Event& e, const edm::EventSetup& context) {
using namespace edm::eventsetup;
// Context is not used.
std::cout << " I AM IN RUN NUMBER " << e.id().run() << std::endl;
std::cout << " ---EVENT NUMBER " << e.id().event() << std::endl;
// get configuration for current run
auto hv = context.getHandle(dthvstatusToken_);
std::cout << hv->version() << std::endl;
std::cout << std::distance(hv->begin(), hv->end()) << " data in the container" << std::endl;
edm::ValidityInterval iov(context.get<DTHVStatusRcd>().validityInterval());
/*
unsigned int currValidityStart = iov.first().eventID().run();
unsigned int currValidityEnd = iov.last( ).eventID().run();
std::cout << "valid since run " << currValidityStart
<< " to run " << currValidityEnd << std::endl;
*/
std::cout << "end of time: " << edm::Timestamp::endOfTime().value() << std::endl;
long long int currValidityStart = iov.first().time().value();
long long int currValidityEnd = iov.last().time().value();
std::cout << "valid since " << currValidityStart << " to " << currValidityEnd << std::endl;
long long int iTimeStart =
((((currValidityStart >> 32) & 0xFFFFFFFF) * 1000000000) + ((currValidityStart & 0xFFFFFFFF) * 1000));
long long int iTimeEnd =
((((currValidityEnd >> 32) & 0xFFFFFFFF) * 1000000000) + ((currValidityEnd & 0xFFFFFFFF) * 1000));
coral::TimeStamp cTimeStart(iTimeStart);
coral::TimeStamp cTimeEnd(iTimeEnd);
std::cout << " ( " << cTimeStart.year() << " , " << cTimeStart.month() << " , " << cTimeStart.day() << " h "
<< cTimeStart.hour() << ":" << cTimeStart.minute() << ":" << cTimeStart.second() << " to "
<< cTimeEnd.year() << " , " << cTimeEnd.month() << " , " << cTimeEnd.day() << " h " << cTimeEnd.hour()
<< ":" << cTimeEnd.minute() << ":" << cTimeEnd.second() << " ) " << std::endl;
/*
DTHVStatus::const_iterator iter = hv->begin();
DTHVStatus::const_iterator iend = hv->end();
while ( iter != iend ) {
const DTHVStatusId& hvId = iter->first;
const DTHVStatusData& hvData = iter->second;
std::cout << hvId.wheelId << " "
<< hvId.stationId << " "
<< hvId.sectorId << " "
<< hvId.slId << " "
<< hvId.layerId << " "
<< hvId.partId << " -> "
<< hvData.fCell << " "
<< hvData.lCell << " "
<< hvData.flagA << " "
<< hvData.flagC << " "
<< hvData.flagS << std::endl;
iter++;
}
*/
std::cout << "============" << std::endl;
int whe;
int sta;
int sec;
int qua;
int lay;
int cel;
int flagA;
int flagC;
int flagS;
std::ifstream fcel("cellList.txt");
while (fcel >> whe >> sta >> sec >> qua >> lay >> cel) {
DTWireId id(whe, sta, sec, qua, lay, cel);
hv->get(id, flagA, flagC, flagS);
std::cout << whe << " " << sta << " " << sec << " " << qua << " " << lay << " " << cel << " -> " << flagA << " "
<< flagC << " " << flagS << std::endl;
}
std::cout << "============" << std::endl;
std::ifstream fcha("chamList.txt");
while (fcha >> whe >> sta >> sec) {
DTChamberId id(whe, sta, sec);
std::cout << "chamber " << whe << " " << sta << " " << sec << " has " << hv->badChannelsNumber(id)
<< " bad cells and " << hv->offChannelsNumber(id) << " off cells" << std::endl;
}
std::cout << "============" << std::endl;
std::cout << "total " << hv->badChannelsNumber() << " bad cells and " << hv->offChannelsNumber() << " off cells"
<< std::endl;
}
DEFINE_FWK_MODULE(DTHVDump);
} // namespace edmtest
|