42 int extension_capacity,
45 int max_lvl1_emu_capacity,
46 int max_rpc_hits_capacity,
47 int max_tgc_hits_capacity,
48 int max_mdt_hits_capacity) :
50 m_id(
id), m_te_id(te_id), m_lvl1_id(lvl1_id), m_lumi_block(lumi_block),
51 m_muondetmask(muondetmask), m_RoIID(RoIID), m_RoISystem(RoISystem),
52 m_RoISubsystem(RoISubsystem), m_RoISector(RoISector), m_RoINumber(RoINumber),
53 m_RoIThreshold(RoIThreshold), m_RoIEta(RoIEta), m_RoIPhi(RoIPhi)
155 source.extension_capacity() : -1;
157 source.max_rob_capacity() : -1;
159 source.max_csm_capacity() : -1;
161 source.max_lvl1_emu_capacity() : -1;
163 source.max_rpc_hits_capacity() : -1;
165 source.max_tgc_hits_capacity() : -1;
167 source.max_mdt_hits_capacity() : -1;
345 std::vector< std::vector<int > > vec_int (2, std::vector<int >(4,99999));
346 std::vector< std::vector<float> > vec_float(2, std::vector<float>(4,99999.));
541 const std::vector<T>&
source,
545 if (resize>0)
v.reserve(resize);
591 float Mid1_r,
float Mid1_z) {
599 float Mid2_r,
float Mid2_z) {
607 float phi_chi2,
long phi_n) {
615 float Inn_r,
float Inn_z) {
623 float phi_std,
long phi_n) {
643 if (tr<0 || tr>1)
return;
644 if (st<0 || st>3)
return;
654 if (tr<0 || tr>1)
return;
655 if (st<0 || st>3)
return;
665 if (tr<0 || tr>1)
return;
666 if (st<0 || st>3)
return;
676 if (tr<0 || tr>1)
return;
677 if (st<0 || st>3)
return;
687 if (tr<0 || tr>1)
return;
688 if (st<0 || st>3)
return;
698 if (tr<0 || tr>1)
return;
699 if (st<0 || st>3)
return;
709 if (tr<0 || tr>1)
return;
710 if (st<0 || st>3)
return;
720 if (tr<0 || tr>1)
return;
721 if (st<0 || st>3)
return;
731 if (tr<0 || tr>1)
return;
732 if (st<0 || st>3)
return;
742 if (tr<0 || tr>1)
return;
743 if (st<0 || st>3)
return;
753 if (tr<0 || tr>1)
return;
754 if (st<0 || st>3)
return;
764 if (tr<0 || tr>1)
return;
765 if (st<0 || st>3)
return;
775 if (tr<0 || tr>1)
return;
776 if (st<0 || st>3)
return;
796 uint32_t myarray [] = { ly00,ly01,ly10,ly11 };
804 uint32_t myarray [] = { ly00,ly01,ly10,ly11 };
820 uint32_t myarray [] = { ly00,ly01,ly10,ly11 };
828 uint32_t myarray [] = { ly00,ly01,ly10,ly11 };
862 float x,
float y,
float z,
float r,
float p) {
1338 int nbit = (
n<=32)?
n : 32;
1339 std::stringstream
ss;
1341 for(
int i=nbit-1;
i>=0;--
i)
ss << ( ((
w>>
i)&0
x1)? 1 : 0 );
1354 case 0:
id =
"900 GEV Instance";
break;
1355 case 1:
id =
"MUON Instance";
break;
1356 case 2:
id =
"HALO Instance";
break;
1357 case 3:
id =
"COSMIC Instance";
break;
1358 default:
id =
"UNKNOWN Instance";
1361 std::stringstream
ss;
1364 ss <<
" *********************************************************" << std::endl;
1365 ss <<
" * Dump of the muFast ESD content *" << std::endl;
1366 ss <<
" *********************************************************" << std::endl;
1367 ss <<
" Instance : " <<
id << std::endl;
1368 ss <<
" TE id : " <<
d.te_id() << std::endl;
1369 ss <<
" Error code : 0x" << std::hex <<
d.error() << std::dec << std::endl;
1370 ss <<
" LVL1 ID : " <<
d.lvl1_id() << std::endl;
1371 ss <<
" Muon DetectorMask : " <<
d.muondetmask() << std::endl;
1372 ss <<
" RoI ID : " <<
d.roi_id() << std::endl;
1373 ss <<
" RoI System : " <<
d.roi_system() << std::endl;
1374 ss <<
" RoI Subsystem : " <<
d.roi_subsystem() << std::endl;
1375 ss <<
" RoI Sector : " <<
d.roi_sector() << std::endl;
1376 ss <<
" RoI Number : " <<
d.roi_number() << std::endl;
1377 ss <<
" RoI Threshold : " <<
d.roi_threshold() << std::endl;
1378 ss <<
" RoI Eta position : " <<
d.roi_eta() << std::endl;
1379 ss <<
" RoI Phi position : " <<
d.roi_phi() << std::endl;
1381 ss <<
" ********* Level-2 Data Access ********" << std::endl;
1382 ss <<
" Requested ROBs :";
1383 for (
unsigned int i=0;
i<
d.rob_id().
size();++
i)
1384 ss <<
" 0x" << std::hex <<
d.rob_id()[
i] << std::dec;
1387 ss <<
" Removed ROBs :";
1388 for (
unsigned int i=0;
i<
d.removed_rob_id().
size();++
i)
1389 ss <<
" 0x" << std::hex <<
d.removed_rob_id()[
i] << std::dec;
1392 ss <<
" RPC PAD access err: 0x" << std::hex <<
d.rpc_pad_error() << std::dec << std::endl;
1393 ss <<
" TGC RDO access err: 0x" << std::hex <<
d.tgc_rdo_error() << std::dec << std::endl;
1395 ss <<
" Requested CSMs :";
1396 for (
unsigned int i=0;
i<
d.csm_id().
size(); ++
i)
1397 ss <<
" " <<
d.csm_id()[
i] <<
"(" <<
d.csm_size()[
i] <<
")";
1400 ss <<
" Removed CSMs :";
1401 for (
unsigned int i=0;
i<
d.removed_csm_id().
size(); ++
i)
1402 ss <<
" " <<
d.removed_csm_id()[
i];
1405 ss <<
" CSMs access err :";
1406 for (
unsigned int i=0;
i<
d.csm_error().
size(); ++
i)
1407 ss <<
" 0x" << std::hex <<
d.csm_error()[
i] << std::dec;
1410 ss <<
" ******* RPC hit reconstruction *******" << std::endl;
1411 ss <<
" RPC point 1 : "
1412 <<
"x = " << std::setw(9) << std::setprecision(3) <<
d.rpc1_x()
1413 <<
", y = " << std::setw(9) << std::setprecision(3) <<
d.rpc1_y()
1414 <<
", z = " << std::setw(9) << std::setprecision(3) <<
d.rpc1_z() << std::endl;
1415 ss <<
" RPC point 2 : "
1416 <<
"x = " << std::setw(9) << std::setprecision(3) <<
d.rpc2_x()
1417 <<
", y = " << std::setw(9) << std::setprecision(3) <<
d.rpc2_y()
1418 <<
", z = " << std::setw(9) << std::setprecision(3) <<
d.rpc2_z() << std::endl;
1419 ss <<
" RPC point 3 : "
1420 <<
"x = " << std::setw(9) << std::setprecision(3) <<
d.rpc3_x()
1421 <<
", y = " << std::setw(9) << std::setprecision(3) <<
d.rpc3_y()
1422 <<
", z = " << std::setw(9) << std::setprecision(3) <<
d.rpc3_z() << std::endl;
1424 if(
d.lvl1_emulation().size()) {
1425 ss <<
" eta_pivot_lay0 : " <<
str_bit(
d.eta_pivot_lay0(),32) << std::endl;
1426 ss <<
" eta_pivot_lay1 : " <<
str_bit(
d.eta_pivot_lay1(),32) << std::endl;
1427 ss <<
" eta_low_0_lay0 : " <<
str_bit(
d.eta_low_0_lay0(),32) << std::endl;
1428 ss <<
" eta_low_1_lay0 : " <<
str_bit(
d.eta_low_1_lay0(),32) << std::endl;
1429 ss <<
" eta_low_0_lay1 : " <<
str_bit(
d.eta_low_0_lay1(),32) << std::endl;
1430 ss <<
" eta_low_1_lay1 : " <<
str_bit(
d.eta_low_1_lay1(),32) << std::endl;
1431 ss <<
" eta_high_0_lay0 : " <<
str_bit(
d.eta_high_0_lay0(),32) << std::endl;
1432 ss <<
" eta_high_1_lay0 : " <<
str_bit(
d.eta_high_1_lay0(),32) << std::endl;
1433 ss <<
" eta_high_0_lay1 : " <<
str_bit(
d.eta_high_0_lay1(),32) << std::endl;
1434 ss <<
" eta_high_1_lay1 : " <<
str_bit(
d.eta_high_1_lay1(),32) << std::endl;
1435 ss <<
" phi_pivot_lay0 : " <<
str_bit(
d.phi_pivot_lay0(),32) << std::endl;
1436 ss <<
" phi_pivot_lay1 : " <<
str_bit(
d.phi_pivot_lay1(),32) << std::endl;
1437 ss <<
" phi_low_0_lay0 : " <<
str_bit(
d.phi_low_0_lay0(),32) << std::endl;
1438 ss <<
" phi_low_1_lay0 : " <<
str_bit(
d.phi_low_1_lay0(),32) << std::endl;
1439 ss <<
" phi_low_0_lay1 : " <<
str_bit(
d.phi_low_0_lay1(),32) << std::endl;
1440 ss <<
" phi_low_1_lay1 : " <<
str_bit(
d.phi_low_1_lay1(),32) << std::endl;
1441 ss <<
" phi_high_0_lay0 : " <<
str_bit(
d.phi_high_0_lay0(),32) << std::endl;
1442 ss <<
" phi_high_1_lay0 : " <<
str_bit(
d.phi_high_1_lay0(),32) << std::endl;
1443 ss <<
" phi_high_0_lay1 : " <<
str_bit(
d.phi_high_0_lay1(),32) << std::endl;
1444 ss <<
" phi_high_1_lay1 : " <<
str_bit(
d.phi_high_1_lay1(),32) << std::endl;
1447 for (
unsigned int i=0;
i<
d.pad_hit_onlineId().
size();++
i)
1449 <<
"onlineId=0x" << std::hex <<
d.pad_hit_onlineId()[
i] << std::dec
1450 <<
", code=" <<
d.pad_hit_code()[
i]
1451 <<
", x=" << std::setw(9) << std::setprecision(3) <<
d.pad_hit_x()[
i]
1452 <<
", y=" << std::setw(9) << std::setprecision(3) <<
d.pad_hit_y()[
i]
1453 <<
", z=" << std::setw(9) << std::setprecision(3) <<
d.pad_hit_z()[
i]
1454 <<
", r=" << std::setw(10) << std::setprecision(3) <<
d.pad_hit_r()[
i]
1455 <<
", p=" << std::setw(6) << std::setprecision(3) <<
d.pad_hit_p()[
i]
1458 ss <<
" ******* TGC hit reconstruction *******" << std::endl;
1459 ss <<
" TGC Mid1 point :"
1460 <<
" eta=" << std::setw(9) << std::setprecision(4) <<
d.tgc_Mid1_eta()
1461 <<
", phi=" << std::setw(9) << std::setprecision(4) <<
d.tgc_Mid1_phi()
1462 <<
", r=" << std::setw(10) << std::setprecision(2) <<
d.tgc_Mid1_r()
1463 <<
", z=" << std::setw(10) << std::setprecision(2) <<
d.tgc_Mid1_z() << std::endl;
1464 ss <<
" TGC Mid2 point :"
1465 <<
" eta=" << std::setw(9) << std::setprecision(4) <<
d.tgc_Mid2_eta()
1466 <<
", phi=" << std::setw(9) << std::setprecision(4) <<
d.tgc_Mid2_phi()
1467 <<
", r=" << std::setw(10) << std::setprecision(2) <<
d.tgc_Mid2_r()
1468 <<
", z=" << std::setw(10) << std::setprecision(2) <<
d.tgc_Mid2_z() << std::endl;
1469 ss <<
" TGC Mid fit :"
1470 <<
" r_chi2=" << std::setw(9) << std::setprecision(1) <<
d.tgc_Mid_rho_chi2()
1471 <<
", r_num=" << std::setw(4) <<
d.tgc_Mid_rho_N()
1472 <<
", p_chi2=" << std::setw(9) << std::setprecision(1) <<
d.tgc_Mid_phi_chi2()
1473 <<
", p_num=" << std::setw(4) <<
d.tgc_Mid_phi_N() << std::endl;
1474 ss <<
" TGC Inner point :"
1475 <<
" eta=" << std::setw(9) << std::setprecision(4) <<
d.tgc_Inn_eta()
1476 <<
", phi=" << std::setw(9) << std::setprecision(4) <<
d.tgc_Inn_phi()
1477 <<
", r=" << std::setw(10) << std::setprecision(2) <<
d.tgc_Inn_r()
1478 <<
", z=" << std::setw(10) << std::setprecision(2) <<
d.tgc_Inn_z() << std::endl;
1479 ss <<
" TGC Inner fit :"
1480 <<
" r_std=" << std::setw(9) << std::setprecision(1) <<
d.tgc_Inn_rho_std()
1481 <<
", r_num=" << std::setw(4) <<
d.tgc_Inn_rho_N()
1482 <<
", p_std=" << std::setw(9) << std::setprecision(1) <<
d.tgc_Inn_phi_std()
1483 <<
", p_num=" << std::setw(4) <<
d.tgc_Inn_phi_N() << std::endl;
1485 ss <<
" TGC PT estimate : " <<
d.tgc_PT() << std::endl;
1487 for (
unsigned int i=0;
i<
d.tgc_Mid_rho_hit_phi().
size();++
i)
1488 ss <<
" TGC Mid rho hit :"
1489 <<
" phi=" << std::setw(9) << std::setprecision(4) <<
d.tgc_Mid_rho_hit_phi()[
i]
1490 <<
", r=" << std::setw(10) << std::setprecision(2) <<
d.tgc_Mid_rho_hit_r()[
i]
1491 <<
", z=" << std::setw(10) << std::setprecision(2) <<
d.tgc_Mid_rho_hit_z()[
i]
1492 <<
", width=" << std::setw(4) << std::setprecision(2) <<
d.tgc_Mid_rho_hit_width()[
i]
1493 <<
", in road=" <<
d.tgc_Mid_rho_hit_in_seg()[
i] << std::endl;
1495 for (
unsigned int i=0;
i<
d.tgc_Mid_phi_hit_phi().
size();++
i)
1496 ss <<
" TGC Mid phi hit :"
1497 <<
" phi=" << std::setw(9) << std::setprecision(4) <<
d.tgc_Mid_phi_hit_phi()[
i]
1498 <<
", r=" << std::setw(10) << std::setprecision(2) <<
d.tgc_Mid_phi_hit_r()[
i]
1499 <<
", z=" << std::setw(10) << std::setprecision(2) <<
d.tgc_Mid_phi_hit_z()[
i]
1500 <<
", width=" << std::setw(4) << std::setprecision(2) <<
d.tgc_Mid_phi_hit_width()[
i]
1501 <<
", in road=" <<
d.tgc_Mid_phi_hit_in_seg()[
i] << std::endl;
1503 for (
unsigned int i=0;
i<
d.tgc_Inn_rho_hit_phi().
size();++
i)
1504 ss <<
" TGC Inn rho hit :"
1505 <<
" phi=" << std::setw(9) << std::setprecision(4) <<
d.tgc_Inn_rho_hit_phi()[
i]
1506 <<
", r=" << std::setw(10) << std::setprecision(2) <<
d.tgc_Inn_rho_hit_r()[
i]
1507 <<
", z=" << std::setw(10) << std::setprecision(2) <<
d.tgc_Inn_rho_hit_z()[
i]
1508 <<
", width=" << std::setw(4) << std::setprecision(2) <<
d.tgc_Inn_rho_hit_width()[
i]
1509 <<
", in road=" <<
d.tgc_Inn_rho_hit_in_seg()[
i] << std::endl;
1511 for (
unsigned int i=0;
i<
d.tgc_Inn_phi_hit_phi().
size();++
i)
1512 ss <<
" TGC Inn rho hit :"
1513 <<
" phi=" << std::setw(9) << std::setprecision(4) <<
d.tgc_Inn_phi_hit_phi()[
i]
1514 <<
", r=" << std::setw(10) << std::setprecision(2) <<
d.tgc_Inn_phi_hit_r()[
i]
1515 <<
", z=" << std::setw(10) << std::setprecision(2) <<
d.tgc_Inn_phi_hit_z()[
i]
1516 <<
", width=" << std::setw(4) << std::setprecision(2) <<
d.tgc_Inn_phi_hit_width()[
i]
1517 <<
", in road=" <<
d.tgc_Inn_phi_hit_in_seg()[
i] << std::endl;
1519 ss <<
" ******** Track reconstruction ********" << std::endl;
1520 ss <<
" Track model : type=" <<
d.type()
1521 <<
", x=" << std::setw(9) << std::setprecision(3) <<
d.t_x_center()
1522 <<
", y=" << std::setw(9) << std::setprecision(3) <<
d.t_y_center()
1523 <<
", rad=" << std::setw(10) << std::setprecision(3) <<
d.t_radius()
1525 ss <<
" Phi position :"
1526 <<
" ap=" << std::setw(9) << std::setprecision(3) <<
d.t_ap()
1527 <<
", aphi=" << std::setw(9) << std::setprecision(3) <<
d.t_aphi()
1528 <<
", phi=" << std::setw(9) << std::setprecision(3) <<
d.t_phi()
1530 ss <<
" Extrapolation :"
1531 <<
" InnerEta=" << std::setw(9) << std::setprecision(3) <<
d.extrInnerEta()
1532 <<
" InnerPhi=" << std::setw(9) << std::setprecision(3) <<
d.extrInnerPhi()
1535 ss <<
" ************* Muon Roads *************" << std::endl;
1536 ss <<
" Inner trigger :"
1537 <<
" ty1=" << (
d.chamber_type_1())[0][0]
1538 <<
", ty2=" <<(
d.chamber_type_2())[0][0]
1539 <<
", pos=" <<(
d.pos())[0][0]
1540 <<
", aw=" << std::setw(9) << std::setprecision(3) << (
d.aw())[0][0]
1541 <<
", bw=" << std::setw(9) << std::setprecision(3) << (
d.bw())[0][0]
1542 <<
", phi=" << std::setw(8) << std::setprecision(3) << (
d.st_phi())[0][0]
1544 ss <<
" Inner trigger :"
1545 <<
" Zmin=" << std::setw(9) << std::setprecision(3) << (
d.zetaMin())[0][0]
1546 <<
", Zmax=" << std::setw(9) << std::setprecision(3) << (
d.zetaMax())[0][0]
1547 <<
", Rmin=" << std::setw(9) << std::setprecision(3) << (
d.radMin())[0][0]
1548 <<
", Rmax=" << std::setw(9) << std::setprecision(3) << (
d.radMax())[0][0]
1549 <<
", Emin=" << std::setw(7) << std::setprecision(3) << (
d.etaMin())[0][0]
1550 <<
", Emax=" << std::setw(7) << std::setprecision(3) << (
d.etaMax())[0][0]
1551 <<
", road=" << std::setw(7) << std::setprecision(3) << (
d.st_roads())[0][0]
1553 ss <<
" Middle trigger :"
1554 <<
" ty1=" << (
d.chamber_type_1())[0][1]
1555 <<
", ty2=" <<(
d.chamber_type_2())[0][1]
1556 <<
", pos=" <<(
d.pos())[0][1]
1557 <<
", aw=" << std::setw(9) << std::setprecision(3) << (
d.aw())[0][1]
1558 <<
", bw=" << std::setw(9) << std::setprecision(3) << (
d.bw())[0][1]
1559 <<
", phi=" << std::setw(8) << std::setprecision(3) << (
d.st_phi())[0][1]
1561 ss <<
" Middle trigger :"
1562 <<
" Zmin=" << std::setw(9) << std::setprecision(3) << (
d.zetaMin())[0][1]
1563 <<
", Zmax=" << std::setw(9) << std::setprecision(3) << (
d.zetaMax())[0][1]
1564 <<
", Rmin=" << std::setw(9) << std::setprecision(3) << (
d.radMin())[0][1]
1565 <<
", Rmax=" << std::setw(9) << std::setprecision(3) << (
d.radMax())[0][1]
1566 <<
", Emin=" << std::setw(7) << std::setprecision(3) << (
d.etaMin())[0][1]
1567 <<
", Emax=" << std::setw(7) << std::setprecision(3) << (
d.etaMax())[0][1]
1568 <<
", road=" << std::setw(7) << std::setprecision(3) << (
d.st_roads())[0][1]
1570 ss <<
" Outer trigger :"
1571 <<
" ty1=" << (
d.chamber_type_1())[0][2]
1572 <<
", ty2=" <<(
d.chamber_type_2())[0][2]
1573 <<
", pos=" <<(
d.pos())[0][2]
1574 <<
", aw=" << std::setw(9) << std::setprecision(3) << (
d.aw())[0][2]
1575 <<
", bw=" << std::setw(9) << std::setprecision(3) << (
d.bw())[0][2]
1576 <<
", phi=" << std::setw(8) << std::setprecision(3) << (
d.st_phi())[0][2]
1578 ss <<
" Outer trigger :"
1579 <<
" Zmin=" << std::setw(9) << std::setprecision(3) << (
d.zetaMin())[0][2]
1580 <<
", Zmax=" << std::setw(9) << std::setprecision(3) << (
d.zetaMax())[0][2]
1581 <<
", Rmin=" << std::setw(9) << std::setprecision(3) << (
d.radMin())[0][2]
1582 <<
", Rmax=" << std::setw(9) << std::setprecision(3) << (
d.radMax())[0][2]
1583 <<
", Emin=" << std::setw(7) << std::setprecision(3) << (
d.etaMin())[0][2]
1584 <<
", Emax=" << std::setw(7) << std::setprecision(3) << (
d.etaMax())[0][2]
1585 <<
", road=" << std::setw(7) << std::setprecision(3) << (
d.st_roads())[0][2]
1587 ss <<
" Ext. trigger :"
1588 <<
" ty1=" << (
d.chamber_type_1())[0][3]
1589 <<
", ty2=" <<(
d.chamber_type_2())[0][3]
1590 <<
", pos=" <<(
d.pos())[0][3]
1591 <<
", aw=" << std::setw(9) << std::setprecision(3) << (
d.aw())[0][3]
1592 <<
", bw=" << std::setw(9) << std::setprecision(3) << (
d.bw())[0][3]
1593 <<
", phi=" << std::setw(8) << std::setprecision(3) << (
d.st_phi())[0][3]
1595 ss <<
" Ext. trigger :"
1596 <<
" Zmin=" << std::setw(9) << std::setprecision(3) << (
d.zetaMin())[0][3]
1597 <<
", Zmax=" << std::setw(9) << std::setprecision(3) << (
d.zetaMax())[0][3]
1598 <<
", Rmin=" << std::setw(9) << std::setprecision(3) << (
d.radMin())[0][3]
1599 <<
", Rmax=" << std::setw(9) << std::setprecision(3) << (
d.radMax())[0][3]
1600 <<
", Emin=" << std::setw(7) << std::setprecision(3) << (
d.etaMin())[0][3]
1601 <<
", Emax=" << std::setw(7) << std::setprecision(3) << (
d.etaMax())[0][3]
1602 <<
", road=" << std::setw(7) << std::setprecision(3) << (
d.st_roads())[0][3]
1605 ss <<
" Inner overlap :"
1606 <<
" ty1=" << (
d.chamber_type_1())[1][0]
1607 <<
", ty2=" <<(
d.chamber_type_2())[1][0]
1608 <<
", pos=" <<(
d.pos())[1][0]
1609 <<
", aw=" << std::setw(9) << std::setprecision(3) << (
d.aw())[1][0]
1610 <<
", bw=" << std::setw(9) << std::setprecision(3) << (
d.bw())[1][0]
1611 <<
", phi=" << std::setw(8) << std::setprecision(3) << (
d.st_phi())[1][0]
1613 ss <<
" Inner overlap :"
1614 <<
" Zmin=" << std::setw(9) << std::setprecision(3) << (
d.zetaMin())[1][0]
1615 <<
", Zmax=" << std::setw(9) << std::setprecision(3) << (
d.zetaMax())[1][0]
1616 <<
", Rmin=" << std::setw(9) << std::setprecision(3) << (
d.radMin())[1][0]
1617 <<
", Rmax=" << std::setw(9) << std::setprecision(3) << (
d.radMax())[1][0]
1618 <<
", Emin=" << std::setw(7) << std::setprecision(3) << (
d.etaMin())[1][0]
1619 <<
", Emax=" << std::setw(7) << std::setprecision(3) << (
d.etaMax())[1][0]
1620 <<
", road=" << std::setw(7) << std::setprecision(3) << (
d.st_roads())[1][0]
1622 ss <<
" Middle overlap :"
1623 <<
" ty1=" << (
d.chamber_type_1())[1][1]
1624 <<
", ty2=" <<(
d.chamber_type_2())[1][1]
1625 <<
", pos=" <<(
d.pos())[1][1]
1626 <<
", aw=" << std::setw(9) << std::setprecision(3) << (
d.aw())[1][1]
1627 <<
", bw=" << std::setw(9) << std::setprecision(3) << (
d.bw())[1][1]
1628 <<
", phi=" << std::setw(8) << std::setprecision(3) << (
d.st_phi())[1][1]
1630 ss <<
" Middle overlap :"
1631 <<
" Zmin=" << std::setw(9) << std::setprecision(3) << (
d.zetaMin())[1][1]
1632 <<
", Zmax=" << std::setw(9) << std::setprecision(3) << (
d.zetaMax())[1][1]
1633 <<
", Rmin=" << std::setw(9) << std::setprecision(3) << (
d.radMin())[1][1]
1634 <<
", Rmax=" << std::setw(9) << std::setprecision(3) << (
d.radMax())[1][1]
1635 <<
", Emin=" << std::setw(7) << std::setprecision(3) << (
d.etaMin())[1][1]
1636 <<
", Emax=" << std::setw(7) << std::setprecision(3) << (
d.etaMax())[1][1]
1637 <<
", road=" << std::setw(7) << std::setprecision(3) << (
d.st_roads())[1][1]
1639 ss <<
" Outer overlap :"
1640 <<
" ty1=" << (
d.chamber_type_1())[1][2]
1641 <<
", ty2=" <<(
d.chamber_type_2())[1][2]
1642 <<
", pos=" <<(
d.pos())[1][2]
1643 <<
", aw=" << std::setw(9) << std::setprecision(3) << (
d.aw())[1][2]
1644 <<
", bw=" << std::setw(9) << std::setprecision(3) << (
d.bw())[1][2]
1645 <<
", phi=" << std::setw(8) << std::setprecision(3) << (
d.st_phi())[1][2]
1647 ss <<
" Outer overlap :"
1648 <<
" Zmin=" << std::setw(9) << std::setprecision(3) << (
d.zetaMin())[1][2]
1649 <<
", Zmax=" << std::setw(9) << std::setprecision(3) << (
d.zetaMax())[1][2]
1650 <<
", Rmin=" << std::setw(9) << std::setprecision(3) << (
d.radMin())[1][2]
1651 <<
", Rmax=" << std::setw(9) << std::setprecision(3) << (
d.radMax())[1][2]
1652 <<
", Emin=" << std::setw(7) << std::setprecision(3) << (
d.etaMin())[1][2]
1653 <<
", Emax=" << std::setw(7) << std::setprecision(3) << (
d.etaMax())[1][2]
1654 <<
", road=" << std::setw(7) << std::setprecision(3) << (
d.st_roads())[1][2]
1656 ss <<
" Ext. overlap :"
1657 <<
" ty1=" << (
d.chamber_type_1())[1][3]
1658 <<
", ty2=" <<(
d.chamber_type_2())[1][3]
1659 <<
", pos=" <<(
d.pos())[1][3]
1660 <<
", aw=" << std::setw(9) << std::setprecision(3) << (
d.aw())[1][3]
1661 <<
", bw=" << std::setw(9) << std::setprecision(3) << (
d.bw())[1][3]
1662 <<
", phi=" << std::setw(8) << std::setprecision(3) << (
d.st_phi())[1][3]
1664 ss <<
" Ext. overlap :"
1665 <<
" Zmin=" << std::setw(9) << std::setprecision(3) << (
d.zetaMin())[1][3]
1666 <<
", Zmax=" << std::setw(9) << std::setprecision(3) << (
d.zetaMax())[1][3]
1667 <<
", Rmin=" << std::setw(9) << std::setprecision(3) << (
d.radMin())[1][3]
1668 <<
", Rmax=" << std::setw(9) << std::setprecision(3) << (
d.radMax())[1][3]
1669 <<
", Emin=" << std::setw(7) << std::setprecision(3) << (
d.etaMin())[1][3]
1670 <<
", Emax=" << std::setw(7) << std::setprecision(3) << (
d.etaMax())[1][3]
1671 <<
", road=" << std::setw(7) << std::setprecision(3) << (
d.st_roads())[1][3]
1676 ss <<
" ******* MDT hit reconstruction *******" << std::endl;
1677 ss <<
" MDT Innermost fit :"
1678 <<
" slope=" << std::setw(9) << std::setprecision(3) <<
d.mdt_Inner_slope()
1679 <<
", interc.=" << std::setw(9) << std::setprecision(3) <<
d.mdt_Inner_intercept()
1680 <<
", z=" << std::setw(9) << std::setprecision(3) <<
d.mdt_Inner_Z()
1681 <<
", r=" << std::setw(9) << std::setprecision(3) <<
d.mdt_Inner_R()
1682 <<
", chi2=" << std::setw(9) << std::setprecision(2) <<
d.mdt_Inner_fit_chi() << std::endl;
1684 ss <<
" MDT Middle fit :"
1685 <<
" slope=" << std::setw(9) << std::setprecision(3) <<
d.mdt_Middle_slope()
1686 <<
", interc.=" << std::setw(9) << std::setprecision(3) <<
d.mdt_Middle_intercept()
1687 <<
", z=" << std::setw(9) << std::setprecision(3) <<
d.mdt_Middle_Z()
1688 <<
", r=" << std::setw(9) << std::setprecision(3) <<
d.mdt_Middle_R()
1689 <<
", chi2=" << std::setw(9) << std::setprecision(2) <<
d.mdt_Middle_fit_chi() << std::endl;
1691 ss <<
" MDT Outermost fit :"
1692 <<
" slope=" << std::setw(9) << std::setprecision(3) <<
d.mdt_Outer_slope()
1693 <<
", interc.=" << std::setw(9) << std::setprecision(3) <<
d.mdt_Outer_intercept()
1694 <<
", z=" << std::setw(9) << std::setprecision(3) <<
d.mdt_Outer_Z()
1695 <<
", r=" << std::setw(9) << std::setprecision(3) <<
d.mdt_Outer_R()
1696 <<
", chi2=" << std::setw(9) << std::setprecision(2) <<
d.mdt_Outer_fit_chi() << std::endl;
1698 for (
unsigned int i=0;
i<
d.mdt_onlineId().
size();++
i)
1700 <<
" onlineId=0x" << std::hex <<
d.mdt_onlineId()[
i] << std::dec
1701 <<
", offlineId=0x" << std::hex <<
d.mdt_offlineId()[
i] << std::dec
1702 <<
", r=" << std::setw(9) << std::setprecision(2) <<
d.mdt_tube_r()[
i]
1703 <<
", z=" << std::setw(9) << std::setprecision(2) <<
d.mdt_tube_z()[
i]
1704 <<
", res=" << std::setw(9) << std::setprecision(2) <<
d.mdt_tube_residual()[
i]
1705 <<
", time=" << std::setw(5) <<
d.mdt_tube_time()[
i]
1706 <<
", space=" <<
d.mdt_tube_space()[
i]
1707 <<
", sigma=" <<
d.mdt_tube_sigma()[
i] << std::endl;
1710 ss <<
" ********** Extended Variable *********" << std::endl;
1711 ss <<
" Extension 0 :";
1712 for (
unsigned int i=0;
i<
d.extension0().
size();++
i)
1713 ss <<
" " <<
d.extension0()[
i];
1715 ss <<
" Extension 1 :";
1716 for (
unsigned int i=0;
i<
d.extension1().
size();++
i)
1717 ss <<
" " <<
d.extension1()[
i];
1719 ss <<
" Extension 2 :";
1720 for (
unsigned int i=0;
i<
d.extension2().
size();++
i)
1721 ss <<
" " <<
d.extension2()[
i];
1723 ss <<
" Extension 3 :";
1724 for (
unsigned int i=0;
i<
d.extension3().
size();++
i)
1725 ss <<
" " <<
d.extension3()[
i];
1727 ss <<
" Extension 4 :";
1728 for (
unsigned int i=0;
i<
d.extension4().
size();++
i)
1729 ss <<
" " <<
d.extension4()[
i];
1731 ss <<
" Extension 5 :";
1732 for (
unsigned int i=0;
i<
d.extension5().
size();++
i)
1733 ss <<
" " <<
d.extension5()[
i];
1735 ss <<
" Extension 6 :";
1736 for (
unsigned int i=0;
i<
d.extension6().
size();++
i)
1737 ss <<
" " <<
d.extension6()[
i];
1739 ss <<
" Extension 7 :";
1740 for (
unsigned int i=0;
i<
d.extension7().
size();++
i)
1741 ss <<
" " <<
d.extension7()[
i];
1743 ss <<
" Extension 8 :";
1744 for (
unsigned int i=0;
i<
d.extension8().
size();++
i)
1745 ss <<
" " <<
d.extension8()[
i];
1747 ss <<
" Extension 9 :";
1748 for (
unsigned int i=0;
i<
d.extension9().
size();++
i)
1749 ss <<
" " <<
d.extension9()[
i];
1753 ss <<
" ********** Feature Extracted *********" << std::endl;
1754 ss <<
" Address : " <<
d.Address() << std::endl;
1755 ss <<
" Sagitta : " <<
d.Sagitta() << std::endl;
1756 ss <<
" Radius : " <<
d.Radius() << std::endl;
1757 ss <<
" Slope : " <<
d.Slope() << std::endl;
1758 ss <<
" Intercept : " <<
d.Intercept() << std::endl;
1759 ss <<
" Alpha : " <<
d.Alpha() << std::endl;
1760 ss <<
" Beta : " <<
d.Beta() << std::endl;
1761 ss <<
" DeltaR : " <<
d.DeltaR() << std::endl;
1762 ss <<
" Speed over c : " <<
d.Speed_over_c() << std::endl;
1763 ss <<
" Eta map : " <<
d.Eta() << std::endl;
1764 ss <<
" Phi map : " <<
d.PhiMap() << std::endl;
1765 ss <<
" Phi : " <<
d.Phi() << std::endl;
1766 ss <<
" Phi direction : " <<
d.PhiDir() << std::endl;
1767 ss <<
" PT : " <<
d.Pt() << std::endl;
1768 ss <<
" Charge : " <<
d.Charge() << std::endl;