|
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
112
113
114
115
116
117
118
119
120
121
|
/*----------------------------------------------------------------------
Toy EDAnalyzer for testing purposes only.
----------------------------------------------------------------------*/
#include "FWCore/Framework/interface/one/EDAnalyzer.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include <cmath>
#include <fstream>
#include <iostream>
#include <map>
#include <sstream>
#include <stdexcept>
#include <string>
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "CondTools/DT/test/validate/DTMtimeValidateDBRead.h"
#include "CondFormats/DTObjects/interface/DTMtime.h"
#include "CondFormats/DataRecord/interface/DTMtimeRcd.h"
#include "CondFormats/DTObjects/interface/DTRecoConditions.h"
#include "CondFormats/DataRecord/interface/DTRecoConditionsVdriftRcd.h"
DTMtimeValidateDBRead::DTMtimeValidateDBRead(edm::ParameterSet const& p)
: dataFileName(p.getParameter<std::string>("chkFile")),
elogFileName(p.getParameter<std::string>("logFile")),
readLegacyVDriftDB(p.getParameter<bool>("readLegacyVDriftDB")),
dtmtTimeToken_(esConsumes()),
dtrecoCondToken_(esConsumes()) {}
DTMtimeValidateDBRead::DTMtimeValidateDBRead(int i) : dtmtTimeToken_(esConsumes()), dtrecoCondToken_(esConsumes()) {}
void DTMtimeValidateDBRead::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;
std::stringstream run_fn;
run_fn << "run" << e.id().run() << dataFileName;
std::ifstream chkFile(run_fn.str().c_str());
std::ofstream logFile(elogFileName.c_str(), std::ios_base::app);
int whe;
int sta;
int sec;
int qua;
float mTime;
float mTrms;
float ckmt;
float ckrms;
if (readLegacyVDriftDB) { //legacy format
auto mT = context.getHandle(dtmtTimeToken_);
std::cout << mT->version() << std::endl;
std::cout << std::distance(mT->begin(), mT->end()) << " data in the container" << std::endl;
int status;
DTMtime::const_iterator iter = mT->begin();
DTMtime::const_iterator iend = mT->end();
while (iter != iend) {
const DTMtimeId& mTId = iter->first;
const DTMtimeData& mTData = iter->second;
status = mT->get(mTId.wheelId, mTId.stationId, mTId.sectorId, mTId.slId, mTime, mTrms, DTTimeUnits::counts);
if (status)
logFile << "ERROR while getting sl Mtime " << mTId.wheelId << " " << mTId.stationId << " " << mTId.sectorId
<< " " << mTId.slId << " , status = " << status << std::endl;
if ((fabs(mTData.mTime - mTime) > 0.01) || (fabs(mTData.mTrms - mTrms) > 0.0001))
logFile << "MISMATCH WHEN READING sl Mtime " << mTId.wheelId << " " << mTId.stationId << " " << mTId.sectorId
<< " " << mTId.slId << " : " << mTData.mTime << " " << mTData.mTrms << " -> " << mTime << " " << mTrms
<< std::endl;
iter++;
}
while (chkFile >> whe >> sta >> sec >> qua >> ckmt >> ckrms) {
status = mT->get(whe, sta, sec, qua, mTime, mTrms, DTTimeUnits::counts);
if ((fabs(ckmt - mTime) > 0.01) || (fabs(ckrms - mTrms) > 0.0001))
logFile << "MISMATCH IN WRITING AND READING sl Mtime " << whe << " " << sta << " " << sec << " " << qua << " : "
<< ckmt << " " << ckrms << " -> " << mTime << " " << mTrms << std::endl;
}
} else {
//hVdrift;
//context.get<DTRecoConditionsVdriftRcd>().get(hVdrift);
const DTRecoConditions* vDriftMap_ = &context.getData(dtrecoCondToken_);
// Consistency check: no parametrization is implemented for the time being
int version = vDriftMap_->version();
if (version != 1) {
throw cms::Exception("Configuration") << "only version 1 is presently supported for VDriftDB";
}
while (chkFile >> whe >> sta >> sec >> qua >> ckmt >> ckrms) {
mTime = vDriftMap_->get(DTWireId(whe, sta, sec, 1, 0, 0));
if ((fabs(ckmt - mTime) > 0.01))
logFile << "MISMATCH IN WRITING AND READING sl Mtime " << whe << " " << sta << " " << sec << " " << qua << " : "
<< ckmt << " " << ckrms << " -> " << mTime << " --- " << std::endl;
}
}
}
void DTMtimeValidateDBRead::endJob() {
std::ifstream logFile(elogFileName.c_str());
char* line = new char[1000];
int errors = 0;
std::cout << "Mtime validation result:" << std::endl;
while (logFile.getline(line, 1000)) {
std::cout << line << std::endl;
errors++;
}
if (!errors) {
std::cout << " ********************************* " << std::endl;
std::cout << " *** *** " << std::endl;
std::cout << " *** NO ERRORS FOUND *** " << std::endl;
std::cout << " *** *** " << std::endl;
std::cout << " ********************************* " << std::endl;
}
return;
}
DEFINE_FWK_MODULE(DTMtimeValidateDBRead);
|