54 {
55 MsgStream log(
msgSvc(), name());
56 SmartDataPtr<RecExtTrackCol> extTracks(eventSvc(),"/Event/Recon/RecExtTrackCol");
57 if( ! extTracks )
58 {
59 log << MSG::ERROR << "Unable to retrieve RecExtTrackCol" << endreq;
60 return StatusCode::FAILURE;
61 } else {
62 log << MSG::DEBUG << "RecExtTrackCol retrieved of size "<< extTracks->size() << endreq;
63 for(RecExtTrackCol::iterator it=extTracks->begin(); it!=extTracks->end(); it++)
64 {
65
66 m_tof1=(*it)->tof1();
67 m_tof1Path=(*it)->tof1Path();
68 m_tof1PosSigmaAlongX=(*it)->tof1PosSigmaAlongX();
69 m_xTof1Pos = (*it)->tof1Position().x();
70 m_yTof1Pos = (*it)->tof1Position().y();
71 m_zTof1Pos = (*it)->tof1Position().z();
72 m_emTof1_11 = (*it)->tof1ErrorMatrix().fast(1,1);
73 m_emTof1_33 = (*it)->tof1ErrorMatrix().fast(3,3);
74 m_emTof1_55 = (*it)->tof1ErrorMatrix().fast(5,5);
75 m_emTof1_66 = (*it)->tof1ErrorMatrix().fast(6,6);
76 m_emTof1_23 = (*it)->tof1ErrorMatrix().fast(2,3);
77 m_emTof1_35 = (*it)->tof1ErrorMatrix().fast(3,5);
78
79 StatusCode status1 = m_tuple1->write();
80 if ( status1.isFailure() ) {
81 log << MSG::ERROR << "Cannot fill Ntuple1" << endreq;
82 }
83 }
84 }
85 return StatusCode::SUCCESS;
86}