|
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
122
123
124
125
126
127
128
129
130
131
|
/*
* see file for a description of this class.
*
* $Date: 2012/02/07 18:35:00 $
* $Revision: 1.6.2.1 $
* \author Paolo Ronchese INFN Padova
*
*/
//-----------------------
// This Class' Header --
//-----------------------
#include "CondTools/DT/test/validate/DTT0ValidateHandler.h"
//-------------------------------
// Collaborating Class Headers --
//-------------------------------
#include "CondFormats/DTObjects/interface/DTT0.h"
//---------------
// C++ Headers --
//---------------
#include <cmath>
#include <iostream>
#include <fstream>
#include <sstream>
//-------------------
// Initializations --
//-------------------
//----------------
// Constructors --
//----------------
DTT0ValidateHandler::DTT0ValidateHandler(const edm::ParameterSet& ps)
: firstRun(ps.getParameter<unsigned int>("firstRun")),
lastRun(ps.getParameter<unsigned int>("lastRun")),
dataVersion(ps.getParameter<std::string>("version")),
dataFileName(ps.getParameter<std::string>("outFile")),
elogFileName(ps.getParameter<std::string>("logFile")) {
std::ofstream logFile(elogFileName.c_str());
}
//--------------
// Destructor --
//--------------
DTT0ValidateHandler::~DTT0ValidateHandler() {}
//--------------
// Operations --
//--------------
void DTT0ValidateHandler::getNewObjects() {
int runNumber = firstRun;
while (runNumber <= lastRun)
addNewObject(runNumber++);
return;
}
void DTT0ValidateHandler::addNewObject(int runNumber) {
DTT0* t0 = new DTT0(dataVersion);
std::stringstream run_fn;
run_fn << "run" << runNumber << dataFileName;
int status = 0;
std::ofstream outFile(run_fn.str().c_str());
std::ofstream logFile(elogFileName.c_str(), std::ios_base::app);
int whe;
int sta;
int sec;
int qua;
int lay;
int cel;
int cur;
int ndt = 20;
float t0mean;
float t0rms;
float ckmean;
float ckrms;
int ckrun = runNumber % 3;
whe = 3;
while (--whe >= -2) {
sta = 5;
while (--sta) {
if (sta == 4)
sec = 15;
else
sec = 13;
while (--sec) {
qua = 4;
while (--qua) {
if ((sta == 4) && (qua == 2))
continue;
lay = 5;
while (--lay) {
cur = ndt;
while (--cur) {
cel = (ckrun ? cur : ndt - cur);
t0mean = random() * 1.0 / 0x0fffffff;
t0rms = random() * 0.2 / 0x7fffffff;
t0mean -= 4.0;
// t0rms /= 4.0;
status = t0->set(whe, sta, sec, qua, lay, cel, t0mean, t0rms, DTTimeUnits::counts);
outFile << whe << " " << sta << " " << sec << " " << qua << " " << lay << " " << cel << " " << t0mean
<< " " << t0rms << std::endl;
if (status)
logFile << "ERROR while setting cell T0 " << whe << " " << sta << " " << sec << " " << qua << " " << lay
<< " " << cel << " , status = " << status << std::endl;
status = t0->get(whe, sta, sec, qua, lay, cel, ckmean, ckrms, DTTimeUnits::counts);
if (status)
logFile << "ERROR while checking cell T0 " << whe << " " << sta << " " << sec << " " << qua << " "
<< lay << " " << cel << " , status = " << status << std::endl;
if ((fabs(ckmean - t0mean) > 0.0001) || (fabs(ckrms - t0rms) > 0.0001))
logFile << "MISMATCH WHEN WRITING cell T0 " << whe << " " << sta << " " << sec << " " << qua << " "
<< lay << " " << cel << " : " << t0mean << " " << t0rms << " -> " << ckmean << " " << ckrms
<< std::endl;
}
}
}
}
}
}
cond::Time_t snc = runNumber;
m_to_transfer.push_back(std::make_pair(t0, snc));
return;
}
std::string DTT0ValidateHandler::id() const { return dataVersion; }
|