1098 std::vector < std::vector<double> > amplwgtv(
nTDCbins_);
1102 std::vector<double> vampl;
1104 vampl.push_back( -0.33333 );
1105 vampl.push_back( -0.33333 );
1106 vampl.push_back( -0.33333 );
1107 vampl.push_back( 0. );
1108 vampl.push_back( 0. );
1109 vampl.push_back( 1. );
1110 vampl.push_back( 0. );
1111 vampl.push_back( 0. );
1112 vampl.push_back( 0. );
1113 vampl.push_back( 0. );
1121 while (!amplFile.eof() && tdcBin <
nTDCbins_)
1123 for(
int j = 0;
j < 10; ++
j) {
1126 amplwgtv[tdcBin].push_back(ww);
1136 edm::LogError(
"EcalTrivialConditionRetriever") <<
"Configuration not supported. Exception is raised ";
1143 assert(amplwgtv[
i].
size() == 10);
1145 for(std::vector<double>::const_iterator it = amplwgtv[
i].
begin(); it != amplwgtv[
i].end(); ++it)
1153 std::vector < std::vector<double> > amplwgtvAftGain(nTDCbins_);
1157 std::vector<double> vamplAftGain;
1158 vamplAftGain.push_back( 0. );
1159 vamplAftGain.push_back( 0. );
1160 vamplAftGain.push_back( 0. );
1161 vamplAftGain.push_back( 0. );
1162 vamplAftGain.push_back( 0. );
1163 vamplAftGain.push_back( 1. );
1164 vamplAftGain.push_back( 0. );
1165 vamplAftGain.push_back( 0. );
1166 vamplAftGain.push_back( 0. );
1167 vamplAftGain.push_back( 0. );
1168 amplwgtvAftGain[0] = ps.
getUntrackedParameter< std::vector<double> >(
"amplWeightsAftGain", vamplAftGain);
1176 while (!amplFile.eof() && tdcBin <
nTDCbins_)
1178 for(
int j = 0; j < 10; ++
j) {
1181 amplwgtvAftGain[tdcBin].push_back(ww);
1185 assert (tdcBin == nTDCbins_);
1190 edm::LogError(
"EcalTrivialConditionRetriever") <<
"Configuration not supported. Exception is raised ";
1196 assert(amplwgtvAftGain[
i].
size() == 10);
1198 for(std::vector<double>::const_iterator it = amplwgtvAftGain[
i].
begin(); it != amplwgtvAftGain[
i].end(); ++it) {
1206 std::vector< std::vector<double> > pedwgtv(nTDCbins_);
1210 std::vector<double> vped;
1211 vped.push_back( 0.33333 );
1212 vped.push_back( 0.33333 );
1213 vped.push_back( 0.33333 );
1214 vped.push_back( 0. );
1215 vped.push_back( 0. );
1216 vped.push_back( 0. );
1217 vped.push_back( 0. );
1218 vped.push_back( 0. );
1219 vped.push_back( 0. );
1220 vped.push_back( 0. );
1229 while (!pedFile.eof() && tdcBin <
nTDCbins_)
1231 for(
int j = 0; j < 10; ++
j) {
1234 pedwgtv[tdcBin].push_back(ww);
1238 assert (tdcBin == nTDCbins_);
1243 edm::LogError(
"EcalTrivialConditionRetriever") <<
"Configuration not supported. Exception is raised ";
1249 assert(pedwgtv[
i].
size() == 10);
1251 for(std::vector<double>::const_iterator it = pedwgtv[
i].
begin(); it != pedwgtv[
i].end(); ++it) {
1257 std::vector< std::vector<double> > pedwgtvaft(nTDCbins_);
1261 std::vector<double> vped;
1262 vped.push_back( 0. );
1263 vped.push_back( 0. );
1264 vped.push_back( 0. );
1265 vped.push_back( 0. );
1266 vped.push_back( 0. );
1267 vped.push_back( 0. );
1268 vped.push_back( 0. );
1269 vped.push_back( 0. );
1270 vped.push_back( 0. );
1271 vped.push_back( 0. );
1280 while (!pedFile.eof() && tdcBin <
nTDCbins_)
1282 for(
int j = 0; j < 10; ++
j) {
1285 pedwgtvaft[tdcBin].push_back(ww);
1289 assert (tdcBin == nTDCbins_);
1294 edm::LogError(
"EcalTrivialConditionRetriever") <<
"Configuration not supported. Exception is raised ";
1300 assert(pedwgtvaft[
i].
size() == 10);
1302 for(std::vector<double>::const_iterator it = pedwgtvaft[
i].
begin(); it != pedwgtvaft[
i].end(); ++it) {
1312 std::vector< std::vector<double> > jittwgtv(nTDCbins_);
1316 std::vector<double> vjitt;
1317 vjitt.push_back( 0.04066309 );
1318 vjitt.push_back( 0.04066309 );
1319 vjitt.push_back( 0.04066309 );
1320 vjitt.push_back( 0.000 );
1321 vjitt.push_back( 1.325176 );
1322 vjitt.push_back( -0.04997078 );
1323 vjitt.push_back( -0.504338 );
1324 vjitt.push_back( -0.5024844 );
1325 vjitt.push_back( -0.3903718 );
1326 vjitt.push_back( 0.000 );
1335 while (!jittFile.eof() && tdcBin <
nTDCbins_)
1337 for(
int j = 0; j < 10; ++
j) {
1340 jittwgtv[tdcBin].push_back(ww);
1344 assert (tdcBin == nTDCbins_);
1349 edm::LogError(
"EcalTrivialConditionRetriever") <<
"Configuration not supported. Exception is raised ";
1355 assert(jittwgtv[
i].
size() == 10);
1357 for(std::vector<double>::const_iterator it = jittwgtv[
i].
begin(); it != jittwgtv[
i].end(); ++it) {
1363 std::vector< std::vector<double> > jittwgtvaft(nTDCbins_);
1367 std::vector<double> vjitt;
1368 vjitt.push_back( 0. );
1369 vjitt.push_back( 0. );
1370 vjitt.push_back( 0. );
1371 vjitt.push_back( 0. );
1372 vjitt.push_back( 1.097871 );
1373 vjitt.push_back( -0.04551035 );
1374 vjitt.push_back( -0.4159156 );
1375 vjitt.push_back( -0.4185352 );
1376 vjitt.push_back( -0.3367127 );
1377 vjitt.push_back( 0. );
1386 while (!jittFile.eof() && tdcBin <
nTDCbins_)
1388 for(
int j = 0; j < 10; ++
j) {
1391 jittwgtvaft[tdcBin].push_back(ww);
1395 assert (tdcBin == nTDCbins_);
1400 edm::LogError(
"EcalTrivialConditionRetriever") <<
"Configuration not supported. Exception is raised ";
1406 assert(jittwgtvaft[
i].
size() == 10);
1408 for(std::vector<double>::const_iterator it = jittwgtvaft[
i].
begin(); it != jittwgtvaft[
i].end(); ++it) {
1415 std::vector< EcalWeightSet::EcalChi2WeightMatrix > chi2Matrix(nTDCbins_);
1422 chi2Matrix[0](0,0) = 0.694371;
1423 chi2Matrix[0](0,1) = -0.305629;
1424 chi2Matrix[0](0,2) = -0.305629;
1425 chi2Matrix[0](0,3) = 0.;
1426 chi2Matrix[0](0,4) = 0.;
1427 chi2Matrix[0](0,5) = 0.;
1428 chi2Matrix[0](0,6) = 0.;
1429 chi2Matrix[0](0,7) = 0.;
1430 chi2Matrix[0](0,8) = 0.;
1431 chi2Matrix[0](0,9) = 0.;
1432 chi2Matrix[0](1,0) = -0.305629;
1433 chi2Matrix[0](1,1) = 0.694371;
1434 chi2Matrix[0](1,2) = -0.305629;
1435 chi2Matrix[0](1,3) = 0.;
1436 chi2Matrix[0](1,4) = 0.;
1437 chi2Matrix[0](1,5) = 0.;
1438 chi2Matrix[0](1,6) = 0.;
1439 chi2Matrix[0](1,7) = 0.;
1440 chi2Matrix[0](1,8) = 0.;
1441 chi2Matrix[0](1,9) = 0.;
1442 chi2Matrix[0](2,0) = -0.305629;
1443 chi2Matrix[0](2,1) = -0.305629;
1444 chi2Matrix[0](2,2) = 0.694371;
1445 chi2Matrix[0](2,3) = 0.;
1446 chi2Matrix[0](2,4) = 0.;
1447 chi2Matrix[0](2,5) = 0.;
1448 chi2Matrix[0](2,6) = 0.;
1449 chi2Matrix[0](2,7) = 0.;
1450 chi2Matrix[0](2,8) = 0.;
1451 chi2Matrix[0](2,9) = 0.;
1452 chi2Matrix[0](3,0) = 0.;
1453 chi2Matrix[0](3,1) = 0.;
1454 chi2Matrix[0](3,2) = 0.;
1455 chi2Matrix[0](3,3) = 0.;
1456 chi2Matrix[0](3,4) = 0.;
1457 chi2Matrix[0](3,5) = 0.;
1458 chi2Matrix[0](3,6) = 0.;
1459 chi2Matrix[0](3,7) = 0.;
1460 chi2Matrix[0](3,8) = 0.;
1461 chi2Matrix[0](3,9) = 0.;
1462 chi2Matrix[0](4,0) = 0.;
1463 chi2Matrix[0](4,1) = 0.;
1464 chi2Matrix[0](4,2) = 0.;
1465 chi2Matrix[0](4,3) = 0.;
1466 chi2Matrix[0](4,4) = 0.8027116;
1467 chi2Matrix[0](4,5) = -0.2517103;
1468 chi2Matrix[0](4,6) = -0.2232882;
1469 chi2Matrix[0](4,7) = -0.1716192;
1470 chi2Matrix[0](4,8) = -0.1239006;
1471 chi2Matrix[0](4,9) = 0.;
1472 chi2Matrix[0](5,0) = 0.;
1473 chi2Matrix[0](5,1) = 0.;
1474 chi2Matrix[0](5,2) = 0.;
1475 chi2Matrix[0](5,3) = 0.;
1476 chi2Matrix[0](5,4) = -0.2517103;
1477 chi2Matrix[0](5,5) = 0.6528964;
1478 chi2Matrix[0](5,6) = -0.2972839;
1479 chi2Matrix[0](5,7) = -0.2067162;
1480 chi2Matrix[0](5,8) = -0.1230729;
1481 chi2Matrix[0](5,9) = 0.;
1482 chi2Matrix[0](6,0) = 0.;
1483 chi2Matrix[0](6,1) = 0.;
1484 chi2Matrix[0](6,2) = 0.;
1485 chi2Matrix[0](6,3) = 0.;
1486 chi2Matrix[0](6,4) = -0.2232882;
1487 chi2Matrix[0](6,5) = -0.2972839;
1488 chi2Matrix[0](6,6) = 0.7413607;
1489 chi2Matrix[0](6,7) = -0.1883866;
1490 chi2Matrix[0](6,8) = -0.1235052;
1491 chi2Matrix[0](6,9) = 0.;
1492 chi2Matrix[0](7,0) = 0.;
1493 chi2Matrix[0](7,1) = 0.;
1494 chi2Matrix[0](7,2) = 0.;
1495 chi2Matrix[0](7,3) = 0.;
1496 chi2Matrix[0](7,4) = -0.1716192;
1497 chi2Matrix[0](7,5) = -0.2067162;
1498 chi2Matrix[0](7,6) = -0.1883866;
1499 chi2Matrix[0](7,7) = 0.844935;
1500 chi2Matrix[0](7,8) = -0.124291;
1501 chi2Matrix[0](7,9) = 0.;
1502 chi2Matrix[0](8,0) = 0.;
1503 chi2Matrix[0](8,1) = 0.;
1504 chi2Matrix[0](8,2) = 0.;
1505 chi2Matrix[0](8,3) = 0.;
1506 chi2Matrix[0](8,4) = -0.1239006;
1507 chi2Matrix[0](8,5) = -0.1230729;
1508 chi2Matrix[0](8,6) = -0.1235052;
1509 chi2Matrix[0](8,7) = -0.124291;
1510 chi2Matrix[0](8,8) = 0.8749833;
1511 chi2Matrix[0](8,9) = 0.;
1512 chi2Matrix[0](9,0) = 0.;
1513 chi2Matrix[0](9,1) = 0.;
1514 chi2Matrix[0](9,2) = 0.;
1515 chi2Matrix[0](9,3) = 0.;
1516 chi2Matrix[0](9,4) = 0.;
1517 chi2Matrix[0](9,5) = 0.;
1518 chi2Matrix[0](9,6) = 0.;
1519 chi2Matrix[0](9,7) = 0.;
1520 chi2Matrix[0](9,8) = 0.;
1521 chi2Matrix[0](9,9) = 0.;
1529 while (!chi2MatrixFile.eof() && tdcBin <
nTDCbins_)
1532 for(
int j = 0; j < 10; ++
j) {
1533 for(
int l = 0;
l < 10; ++
l) {
1535 chi2MatrixFile >> ww;
1536 chi2Matrix[tdcBin](
j,
l)=ww;
1541 assert (tdcBin == nTDCbins_);
1546 edm::LogError(
"EcalTrivialConditionRetriever") <<
"Configuration not supported. Exception is raised ";
1556 std::vector< EcalWeightSet::EcalChi2WeightMatrix > chi2MatrixAft(nTDCbins_);
1563 chi2MatrixAft[0](0,0) = 0.;
1564 chi2MatrixAft[0](0,1) = 0.;
1565 chi2MatrixAft[0](0,2) = 0.;
1566 chi2MatrixAft[0](0,3) = 0.;
1567 chi2MatrixAft[0](0,4) = 0.;
1568 chi2MatrixAft[0](0,5) = 0.;
1569 chi2MatrixAft[0](0,6) = 0.;
1570 chi2MatrixAft[0](0,7) = 0.;
1571 chi2MatrixAft[0](0,8) = 0.;
1572 chi2MatrixAft[0](0,9) = 0.;
1573 chi2MatrixAft[0](1,0) = 0.;
1574 chi2MatrixAft[0](1,1) = 0.;
1575 chi2MatrixAft[0](1,2) = 0.;
1576 chi2MatrixAft[0](1,3) = 0.;
1577 chi2MatrixAft[0](1,4) = 0.;
1578 chi2MatrixAft[0](1,5) = 0.;
1579 chi2MatrixAft[0](1,6) = 0.;
1580 chi2MatrixAft[0](1,7) = 0.;
1581 chi2MatrixAft[0](1,8) = 0.;
1582 chi2MatrixAft[0](1,9) = 0.;
1583 chi2MatrixAft[0](2,0) = 0.;
1584 chi2MatrixAft[0](2,1) = 0.;
1585 chi2MatrixAft[0](2,2) = 0.;
1586 chi2MatrixAft[0](2,3) = 0.;
1587 chi2MatrixAft[0](2,4) = 0.;
1588 chi2MatrixAft[0](2,5) = 0.;
1589 chi2MatrixAft[0](2,6) = 0.;
1590 chi2MatrixAft[0](2,7) = 0.;
1591 chi2MatrixAft[0](2,8) = 0.;
1592 chi2MatrixAft[0](2,9) = 0.;
1593 chi2MatrixAft[0](3,0) = 0.;
1594 chi2MatrixAft[0](3,1) = 0.;
1595 chi2MatrixAft[0](3,2) = 0.;
1596 chi2MatrixAft[0](3,3) = 0.;
1597 chi2MatrixAft[0](3,4) = 0.;
1598 chi2MatrixAft[0](3,5) = 0.;
1599 chi2MatrixAft[0](3,6) = 0.;
1600 chi2MatrixAft[0](3,7) = 0.;
1601 chi2MatrixAft[0](3,8) = 0.;
1602 chi2MatrixAft[0](3,9) = 0.;
1603 chi2MatrixAft[0](4,0) = 0.;
1604 chi2MatrixAft[0](4,1) = 0.;
1605 chi2MatrixAft[0](4,2) = 0.;
1606 chi2MatrixAft[0](4,3) = 0.;
1607 chi2MatrixAft[0](4,4) = 0.8030884;
1608 chi2MatrixAft[0](4,5) = -0.2543541;
1609 chi2MatrixAft[0](4,6) = -0.2243544;
1610 chi2MatrixAft[0](4,7) = -0.1698177;
1611 chi2MatrixAft[0](4,8) = -0.1194506;
1612 chi2MatrixAft[0](4,9) = 0.;
1613 chi2MatrixAft[0](5,0) = 0.;
1614 chi2MatrixAft[0](5,1) = 0.;
1615 chi2MatrixAft[0](5,2) = 0.;
1616 chi2MatrixAft[0](5,3) = 0.;
1617 chi2MatrixAft[0](5,4) = -0.2543541;
1618 chi2MatrixAft[0](5,5) = 0.6714465;
1619 chi2MatrixAft[0](5,6) = -0.2898025;
1620 chi2MatrixAft[0](5,7) = -0.2193564;
1621 chi2MatrixAft[0](5,8) = -0.1542964;
1622 chi2MatrixAft[0](5,9) = 0.;
1623 chi2MatrixAft[0](6,0) = 0.;
1624 chi2MatrixAft[0](6,1) = 0.;
1625 chi2MatrixAft[0](6,2) = 0.;
1626 chi2MatrixAft[0](6,3) = 0.;
1627 chi2MatrixAft[0](6,4) = -0.2243544;
1628 chi2MatrixAft[0](6,5) = -0.2898025;
1629 chi2MatrixAft[0](6,6) = 0.7443781;
1630 chi2MatrixAft[0](6,7) = -0.1934846;
1631 chi2MatrixAft[0](6,8) = -0.136098;
1632 chi2MatrixAft[0](6,9) = 0.;
1633 chi2MatrixAft[0](7,0) = 0.;
1634 chi2MatrixAft[0](7,1) = 0.;
1635 chi2MatrixAft[0](7,2) = 0.;
1636 chi2MatrixAft[0](7,3) = 0.;
1637 chi2MatrixAft[0](7,4) = -0.1698177;
1638 chi2MatrixAft[0](7,5) = -0.2193564;
1639 chi2MatrixAft[0](7,6) = -0.1934846;
1640 chi2MatrixAft[0](7,7) = 0.8535482;
1641 chi2MatrixAft[0](7,8) = -0.1030149;
1642 chi2MatrixAft[0](7,9) = 0.;
1643 chi2MatrixAft[0](8,0) = 0.;
1644 chi2MatrixAft[0](8,1) = 0.;
1645 chi2MatrixAft[0](8,2) = 0.;
1646 chi2MatrixAft[0](8,3) = 0.;
1647 chi2MatrixAft[0](8,4) = -0.1194506;
1648 chi2MatrixAft[0](8,5) = -0.1542964;
1649 chi2MatrixAft[0](8,6) = -0.136098;
1650 chi2MatrixAft[0](8,7) = -0.1030149;
1651 chi2MatrixAft[0](8,8) = 0.9275388;
1652 chi2MatrixAft[0](8,9) = 0.;
1653 chi2MatrixAft[0](9,0) = 0.;
1654 chi2MatrixAft[0](9,1) = 0.;
1655 chi2MatrixAft[0](9,2) = 0.;
1656 chi2MatrixAft[0](9,3) = 0.;
1657 chi2MatrixAft[0](9,4) = 0.;
1658 chi2MatrixAft[0](9,5) = 0.;
1659 chi2MatrixAft[0](9,6) = 0.;
1660 chi2MatrixAft[0](9,7) = 0.;
1661 chi2MatrixAft[0](9,8) = 0.;
1662 chi2MatrixAft[0](9,9) = 0.;
1670 while (!chi2MatrixAftFile.eof() && tdcBin <
nTDCbins_)
1673 for(
int j = 0; j < 10; ++
j) {
1674 for(
int l = 0;
l < 10; ++
l) {
1676 chi2MatrixAftFile >> ww;
1677 chi2MatrixAft[tdcBin](
j,
l)=ww;
1682 assert (tdcBin == nTDCbins_);
1687 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_