930 std::vector < std::vector<double> > amplwgtv(
nTDCbins_);
934 std::vector<double> vampl;
936 vampl.push_back( -0.33333 );
937 vampl.push_back( -0.33333 );
938 vampl.push_back( -0.33333 );
939 vampl.push_back( 0. );
940 vampl.push_back( 0. );
941 vampl.push_back( 1. );
942 vampl.push_back( 0. );
943 vampl.push_back( 0. );
944 vampl.push_back( 0. );
945 vampl.push_back( 0. );
953 while (!amplFile.eof() && tdcBin <
nTDCbins_)
955 for(
int j = 0;
j < 10; ++
j) {
958 amplwgtv[tdcBin].push_back(ww);
968 edm::LogError(
"EcalTrivialConditionRetriever") <<
"Configuration not supported. Exception is raised ";
975 assert(amplwgtv[
i].
size() == 10);
977 for(std::vector<double>::const_iterator it = amplwgtv[
i].
begin(); it != amplwgtv[
i].end(); ++it)
985 std::vector < std::vector<double> > amplwgtvAftGain(nTDCbins_);
989 std::vector<double> vamplAftGain;
990 vamplAftGain.push_back( 0. );
991 vamplAftGain.push_back( 0. );
992 vamplAftGain.push_back( 0. );
993 vamplAftGain.push_back( 0. );
994 vamplAftGain.push_back( 0. );
995 vamplAftGain.push_back( 1. );
996 vamplAftGain.push_back( 0. );
997 vamplAftGain.push_back( 0. );
998 vamplAftGain.push_back( 0. );
999 vamplAftGain.push_back( 0. );
1000 amplwgtvAftGain[0] = ps.
getUntrackedParameter< std::vector<double> >(
"amplWeightsAftGain", vamplAftGain);
1008 while (!amplFile.eof() && tdcBin <
nTDCbins_)
1010 for(
int j = 0; j < 10; ++
j) {
1013 amplwgtvAftGain[tdcBin].push_back(ww);
1017 assert (tdcBin == nTDCbins_);
1022 edm::LogError(
"EcalTrivialConditionRetriever") <<
"Configuration not supported. Exception is raised ";
1028 assert(amplwgtvAftGain[
i].
size() == 10);
1030 for(std::vector<double>::const_iterator it = amplwgtvAftGain[
i].
begin(); it != amplwgtvAftGain[
i].end(); ++it) {
1038 std::vector< std::vector<double> > pedwgtv(nTDCbins_);
1042 std::vector<double> vped;
1043 vped.push_back( 0.33333 );
1044 vped.push_back( 0.33333 );
1045 vped.push_back( 0.33333 );
1046 vped.push_back( 0. );
1047 vped.push_back( 0. );
1048 vped.push_back( 0. );
1049 vped.push_back( 0. );
1050 vped.push_back( 0. );
1051 vped.push_back( 0. );
1052 vped.push_back( 0. );
1061 while (!pedFile.eof() && tdcBin <
nTDCbins_)
1063 for(
int j = 0; j < 10; ++
j) {
1066 pedwgtv[tdcBin].push_back(ww);
1070 assert (tdcBin == nTDCbins_);
1075 edm::LogError(
"EcalTrivialConditionRetriever") <<
"Configuration not supported. Exception is raised ";
1081 assert(pedwgtv[
i].
size() == 10);
1083 for(std::vector<double>::const_iterator it = pedwgtv[
i].
begin(); it != pedwgtv[
i].end(); ++it) {
1089 std::vector< std::vector<double> > pedwgtvaft(nTDCbins_);
1093 std::vector<double> vped;
1094 vped.push_back( 0. );
1095 vped.push_back( 0. );
1096 vped.push_back( 0. );
1097 vped.push_back( 0. );
1098 vped.push_back( 0. );
1099 vped.push_back( 0. );
1100 vped.push_back( 0. );
1101 vped.push_back( 0. );
1102 vped.push_back( 0. );
1103 vped.push_back( 0. );
1112 while (!pedFile.eof() && tdcBin <
nTDCbins_)
1114 for(
int j = 0; j < 10; ++
j) {
1117 pedwgtvaft[tdcBin].push_back(ww);
1121 assert (tdcBin == nTDCbins_);
1126 edm::LogError(
"EcalTrivialConditionRetriever") <<
"Configuration not supported. Exception is raised ";
1132 assert(pedwgtvaft[
i].
size() == 10);
1134 for(std::vector<double>::const_iterator it = pedwgtvaft[
i].
begin(); it != pedwgtvaft[
i].end(); ++it) {
1144 std::vector< std::vector<double> > jittwgtv(nTDCbins_);
1148 std::vector<double> vjitt;
1149 vjitt.push_back( 0.04066309 );
1150 vjitt.push_back( 0.04066309 );
1151 vjitt.push_back( 0.04066309 );
1152 vjitt.push_back( 0.000 );
1153 vjitt.push_back( 1.325176 );
1154 vjitt.push_back( -0.04997078 );
1155 vjitt.push_back( -0.504338 );
1156 vjitt.push_back( -0.5024844 );
1157 vjitt.push_back( -0.3903718 );
1158 vjitt.push_back( 0.000 );
1167 while (!jittFile.eof() && tdcBin <
nTDCbins_)
1169 for(
int j = 0; j < 10; ++
j) {
1172 jittwgtv[tdcBin].push_back(ww);
1176 assert (tdcBin == nTDCbins_);
1181 edm::LogError(
"EcalTrivialConditionRetriever") <<
"Configuration not supported. Exception is raised ";
1187 assert(jittwgtv[
i].
size() == 10);
1189 for(std::vector<double>::const_iterator it = jittwgtv[
i].
begin(); it != jittwgtv[
i].end(); ++it) {
1195 std::vector< std::vector<double> > jittwgtvaft(nTDCbins_);
1199 std::vector<double> vjitt;
1200 vjitt.push_back( 0. );
1201 vjitt.push_back( 0. );
1202 vjitt.push_back( 0. );
1203 vjitt.push_back( 0. );
1204 vjitt.push_back( 1.097871 );
1205 vjitt.push_back( -0.04551035 );
1206 vjitt.push_back( -0.4159156 );
1207 vjitt.push_back( -0.4185352 );
1208 vjitt.push_back( -0.3367127 );
1209 vjitt.push_back( 0. );
1218 while (!jittFile.eof() && tdcBin <
nTDCbins_)
1220 for(
int j = 0; j < 10; ++
j) {
1223 jittwgtvaft[tdcBin].push_back(ww);
1227 assert (tdcBin == nTDCbins_);
1232 edm::LogError(
"EcalTrivialConditionRetriever") <<
"Configuration not supported. Exception is raised ";
1238 assert(jittwgtvaft[
i].
size() == 10);
1240 for(std::vector<double>::const_iterator it = jittwgtvaft[
i].
begin(); it != jittwgtvaft[
i].end(); ++it) {
1247 std::vector< EcalWeightSet::EcalChi2WeightMatrix > chi2Matrix(nTDCbins_);
1254 chi2Matrix[0](0,0) = 0.694371;
1255 chi2Matrix[0](0,1) = -0.305629;
1256 chi2Matrix[0](0,2) = -0.305629;
1257 chi2Matrix[0](0,3) = 0.;
1258 chi2Matrix[0](0,4) = 0.;
1259 chi2Matrix[0](0,5) = 0.;
1260 chi2Matrix[0](0,6) = 0.;
1261 chi2Matrix[0](0,7) = 0.;
1262 chi2Matrix[0](0,8) = 0.;
1263 chi2Matrix[0](0,9) = 0.;
1264 chi2Matrix[0](1,0) = -0.305629;
1265 chi2Matrix[0](1,1) = 0.694371;
1266 chi2Matrix[0](1,2) = -0.305629;
1267 chi2Matrix[0](1,3) = 0.;
1268 chi2Matrix[0](1,4) = 0.;
1269 chi2Matrix[0](1,5) = 0.;
1270 chi2Matrix[0](1,6) = 0.;
1271 chi2Matrix[0](1,7) = 0.;
1272 chi2Matrix[0](1,8) = 0.;
1273 chi2Matrix[0](1,9) = 0.;
1274 chi2Matrix[0](2,0) = -0.305629;
1275 chi2Matrix[0](2,1) = -0.305629;
1276 chi2Matrix[0](2,2) = 0.694371;
1277 chi2Matrix[0](2,3) = 0.;
1278 chi2Matrix[0](2,4) = 0.;
1279 chi2Matrix[0](2,5) = 0.;
1280 chi2Matrix[0](2,6) = 0.;
1281 chi2Matrix[0](2,7) = 0.;
1282 chi2Matrix[0](2,8) = 0.;
1283 chi2Matrix[0](2,9) = 0.;
1284 chi2Matrix[0](3,0) = 0.;
1285 chi2Matrix[0](3,1) = 0.;
1286 chi2Matrix[0](3,2) = 0.;
1287 chi2Matrix[0](3,3) = 0.;
1288 chi2Matrix[0](3,4) = 0.;
1289 chi2Matrix[0](3,5) = 0.;
1290 chi2Matrix[0](3,6) = 0.;
1291 chi2Matrix[0](3,7) = 0.;
1292 chi2Matrix[0](3,8) = 0.;
1293 chi2Matrix[0](3,9) = 0.;
1294 chi2Matrix[0](4,0) = 0.;
1295 chi2Matrix[0](4,1) = 0.;
1296 chi2Matrix[0](4,2) = 0.;
1297 chi2Matrix[0](4,3) = 0.;
1298 chi2Matrix[0](4,4) = 0.8027116;
1299 chi2Matrix[0](4,5) = -0.2517103;
1300 chi2Matrix[0](4,6) = -0.2232882;
1301 chi2Matrix[0](4,7) = -0.1716192;
1302 chi2Matrix[0](4,8) = -0.1239006;
1303 chi2Matrix[0](4,9) = 0.;
1304 chi2Matrix[0](5,0) = 0.;
1305 chi2Matrix[0](5,1) = 0.;
1306 chi2Matrix[0](5,2) = 0.;
1307 chi2Matrix[0](5,3) = 0.;
1308 chi2Matrix[0](5,4) = -0.2517103;
1309 chi2Matrix[0](5,5) = 0.6528964;
1310 chi2Matrix[0](5,6) = -0.2972839;
1311 chi2Matrix[0](5,7) = -0.2067162;
1312 chi2Matrix[0](5,8) = -0.1230729;
1313 chi2Matrix[0](5,9) = 0.;
1314 chi2Matrix[0](6,0) = 0.;
1315 chi2Matrix[0](6,1) = 0.;
1316 chi2Matrix[0](6,2) = 0.;
1317 chi2Matrix[0](6,3) = 0.;
1318 chi2Matrix[0](6,4) = -0.2232882;
1319 chi2Matrix[0](6,5) = -0.2972839;
1320 chi2Matrix[0](6,6) = 0.7413607;
1321 chi2Matrix[0](6,7) = -0.1883866;
1322 chi2Matrix[0](6,8) = -0.1235052;
1323 chi2Matrix[0](6,9) = 0.;
1324 chi2Matrix[0](7,0) = 0.;
1325 chi2Matrix[0](7,1) = 0.;
1326 chi2Matrix[0](7,2) = 0.;
1327 chi2Matrix[0](7,3) = 0.;
1328 chi2Matrix[0](7,4) = -0.1716192;
1329 chi2Matrix[0](7,5) = -0.2067162;
1330 chi2Matrix[0](7,6) = -0.1883866;
1331 chi2Matrix[0](7,7) = 0.844935;
1332 chi2Matrix[0](7,8) = -0.124291;
1333 chi2Matrix[0](7,9) = 0.;
1334 chi2Matrix[0](8,0) = 0.;
1335 chi2Matrix[0](8,1) = 0.;
1336 chi2Matrix[0](8,2) = 0.;
1337 chi2Matrix[0](8,3) = 0.;
1338 chi2Matrix[0](8,4) = -0.1239006;
1339 chi2Matrix[0](8,5) = -0.1230729;
1340 chi2Matrix[0](8,6) = -0.1235052;
1341 chi2Matrix[0](8,7) = -0.124291;
1342 chi2Matrix[0](8,8) = 0.8749833;
1343 chi2Matrix[0](8,9) = 0.;
1344 chi2Matrix[0](9,0) = 0.;
1345 chi2Matrix[0](9,1) = 0.;
1346 chi2Matrix[0](9,2) = 0.;
1347 chi2Matrix[0](9,3) = 0.;
1348 chi2Matrix[0](9,4) = 0.;
1349 chi2Matrix[0](9,5) = 0.;
1350 chi2Matrix[0](9,6) = 0.;
1351 chi2Matrix[0](9,7) = 0.;
1352 chi2Matrix[0](9,8) = 0.;
1353 chi2Matrix[0](9,9) = 0.;
1361 while (!chi2MatrixFile.eof() && tdcBin <
nTDCbins_)
1364 for(
int j = 0; j < 10; ++
j) {
1365 for(
int l = 0;
l < 10; ++
l) {
1367 chi2MatrixFile >> ww;
1368 chi2Matrix[tdcBin](
j,
l)=ww;
1373 assert (tdcBin == nTDCbins_);
1378 edm::LogError(
"EcalTrivialConditionRetriever") <<
"Configuration not supported. Exception is raised ";
1388 std::vector< EcalWeightSet::EcalChi2WeightMatrix > chi2MatrixAft(nTDCbins_);
1395 chi2MatrixAft[0](0,0) = 0.;
1396 chi2MatrixAft[0](0,1) = 0.;
1397 chi2MatrixAft[0](0,2) = 0.;
1398 chi2MatrixAft[0](0,3) = 0.;
1399 chi2MatrixAft[0](0,4) = 0.;
1400 chi2MatrixAft[0](0,5) = 0.;
1401 chi2MatrixAft[0](0,6) = 0.;
1402 chi2MatrixAft[0](0,7) = 0.;
1403 chi2MatrixAft[0](0,8) = 0.;
1404 chi2MatrixAft[0](0,9) = 0.;
1405 chi2MatrixAft[0](1,0) = 0.;
1406 chi2MatrixAft[0](1,1) = 0.;
1407 chi2MatrixAft[0](1,2) = 0.;
1408 chi2MatrixAft[0](1,3) = 0.;
1409 chi2MatrixAft[0](1,4) = 0.;
1410 chi2MatrixAft[0](1,5) = 0.;
1411 chi2MatrixAft[0](1,6) = 0.;
1412 chi2MatrixAft[0](1,7) = 0.;
1413 chi2MatrixAft[0](1,8) = 0.;
1414 chi2MatrixAft[0](1,9) = 0.;
1415 chi2MatrixAft[0](2,0) = 0.;
1416 chi2MatrixAft[0](2,1) = 0.;
1417 chi2MatrixAft[0](2,2) = 0.;
1418 chi2MatrixAft[0](2,3) = 0.;
1419 chi2MatrixAft[0](2,4) = 0.;
1420 chi2MatrixAft[0](2,5) = 0.;
1421 chi2MatrixAft[0](2,6) = 0.;
1422 chi2MatrixAft[0](2,7) = 0.;
1423 chi2MatrixAft[0](2,8) = 0.;
1424 chi2MatrixAft[0](2,9) = 0.;
1425 chi2MatrixAft[0](3,0) = 0.;
1426 chi2MatrixAft[0](3,1) = 0.;
1427 chi2MatrixAft[0](3,2) = 0.;
1428 chi2MatrixAft[0](3,3) = 0.;
1429 chi2MatrixAft[0](3,4) = 0.;
1430 chi2MatrixAft[0](3,5) = 0.;
1431 chi2MatrixAft[0](3,6) = 0.;
1432 chi2MatrixAft[0](3,7) = 0.;
1433 chi2MatrixAft[0](3,8) = 0.;
1434 chi2MatrixAft[0](3,9) = 0.;
1435 chi2MatrixAft[0](4,0) = 0.;
1436 chi2MatrixAft[0](4,1) = 0.;
1437 chi2MatrixAft[0](4,2) = 0.;
1438 chi2MatrixAft[0](4,3) = 0.;
1439 chi2MatrixAft[0](4,4) = 0.8030884;
1440 chi2MatrixAft[0](4,5) = -0.2543541;
1441 chi2MatrixAft[0](4,6) = -0.2243544;
1442 chi2MatrixAft[0](4,7) = -0.1698177;
1443 chi2MatrixAft[0](4,8) = -0.1194506;
1444 chi2MatrixAft[0](4,9) = 0.;
1445 chi2MatrixAft[0](5,0) = 0.;
1446 chi2MatrixAft[0](5,1) = 0.;
1447 chi2MatrixAft[0](5,2) = 0.;
1448 chi2MatrixAft[0](5,3) = 0.;
1449 chi2MatrixAft[0](5,4) = -0.2543541;
1450 chi2MatrixAft[0](5,5) = 0.6714465;
1451 chi2MatrixAft[0](5,6) = -0.2898025;
1452 chi2MatrixAft[0](5,7) = -0.2193564;
1453 chi2MatrixAft[0](5,8) = -0.1542964;
1454 chi2MatrixAft[0](5,9) = 0.;
1455 chi2MatrixAft[0](6,0) = 0.;
1456 chi2MatrixAft[0](6,1) = 0.;
1457 chi2MatrixAft[0](6,2) = 0.;
1458 chi2MatrixAft[0](6,3) = 0.;
1459 chi2MatrixAft[0](6,4) = -0.2243544;
1460 chi2MatrixAft[0](6,5) = -0.2898025;
1461 chi2MatrixAft[0](6,6) = 0.7443781;
1462 chi2MatrixAft[0](6,7) = -0.1934846;
1463 chi2MatrixAft[0](6,8) = -0.136098;
1464 chi2MatrixAft[0](6,9) = 0.;
1465 chi2MatrixAft[0](7,0) = 0.;
1466 chi2MatrixAft[0](7,1) = 0.;
1467 chi2MatrixAft[0](7,2) = 0.;
1468 chi2MatrixAft[0](7,3) = 0.;
1469 chi2MatrixAft[0](7,4) = -0.1698177;
1470 chi2MatrixAft[0](7,5) = -0.2193564;
1471 chi2MatrixAft[0](7,6) = -0.1934846;
1472 chi2MatrixAft[0](7,7) = 0.8535482;
1473 chi2MatrixAft[0](7,8) = -0.1030149;
1474 chi2MatrixAft[0](7,9) = 0.;
1475 chi2MatrixAft[0](8,0) = 0.;
1476 chi2MatrixAft[0](8,1) = 0.;
1477 chi2MatrixAft[0](8,2) = 0.;
1478 chi2MatrixAft[0](8,3) = 0.;
1479 chi2MatrixAft[0](8,4) = -0.1194506;
1480 chi2MatrixAft[0](8,5) = -0.1542964;
1481 chi2MatrixAft[0](8,6) = -0.136098;
1482 chi2MatrixAft[0](8,7) = -0.1030149;
1483 chi2MatrixAft[0](8,8) = 0.9275388;
1484 chi2MatrixAft[0](8,9) = 0.;
1485 chi2MatrixAft[0](9,0) = 0.;
1486 chi2MatrixAft[0](9,1) = 0.;
1487 chi2MatrixAft[0](9,2) = 0.;
1488 chi2MatrixAft[0](9,3) = 0.;
1489 chi2MatrixAft[0](9,4) = 0.;
1490 chi2MatrixAft[0](9,5) = 0.;
1491 chi2MatrixAft[0](9,6) = 0.;
1492 chi2MatrixAft[0](9,7) = 0.;
1493 chi2MatrixAft[0](9,8) = 0.;
1494 chi2MatrixAft[0](9,9) = 0.;
1502 while (!chi2MatrixAftFile.eof() && tdcBin <
nTDCbins_)
1505 for(
int j = 0; j < 10; ++
j) {
1506 for(
int l = 0;
l < 10; ++
l) {
1508 chi2MatrixAftFile >> ww;
1509 chi2MatrixAft[tdcBin](
j,
l)=ww;
1514 assert (tdcBin == nTDCbins_);
1519 edm::LogError(
"EcalTrivialConditionRetriever") <<
"Configuration not supported. Exception is raised ";
T getUntrackedParameter(std::string const &, T const &) const
std::vector< ROOT::Math::SVector< double, EcalDataFrame::MAXSAMPLES > > jittWeightsAft_
std::vector< ROOT::Math::SVector< double, EcalDataFrame::MAXSAMPLES > > pedWeightsAft_
std::string jittWeightsAftFile_
std::string amplWeightsFile_
std::vector< ROOT::Math::SVector< double, EcalDataFrame::MAXSAMPLES > > pedWeights_
std::string pedWeightsAftFile_
std::vector< ROOT::Math::SVector< double, EcalDataFrame::MAXSAMPLES > > amplWeights_
std::vector< EcalWeightSet::EcalChi2WeightMatrix > chi2MatrixAft_
std::string chi2MatrixAftFile_
std::vector< ROOT::Math::SVector< double, EcalDataFrame::MAXSAMPLES > > amplWeightsAft_
std::string amplWeightsAftFile_
std::string fullPath() const
std::string chi2MatrixFile_
tuple size
Write out results.
std::string pedWeightsFile_
std::string jittWeightsFile_
std::vector< EcalWeightSet::EcalChi2WeightMatrix > chi2Matrix_
std::vector< ROOT::Math::SVector< double, EcalDataFrame::MAXSAMPLES > > jittWeights_