CgemAliRecTrk.cxx

Go to the documentation of this file.

#include "CgemAlignAlg/CgemAliRecTrk.h"
 
#include "GaudiKernel/MsgStream.h"
#include "GaudiKernel/IMessageSvc.h"
#include "GaudiKernel/StatusCode.h"
#include "GaudiKernel/ISvcLocator.h"
#include "GaudiKernel/Bootstrap.h"
#include "GaudiKernel/SmartDataPtr.h"
#include "GaudiKernel/IDataProviderSvc.h"
 
const double CgemAliRecTrk::m_df = 10.0;
 
CgemAliRecTrk::~CgemAliRecTrk(){
     unsigned int i;
     for(i=0; i<m_rechit.size(); i++){
      delete m_rechit[i];
     }
     m_rechit.clear();
}
 
void CgemAliRecTrk::setRecTrk(RecMdcTrackCol::iterator it_trk){
     IMessageSvc *msgSvc;
     Gaudi::svcLocator()->service("MessageSvc", msgSvc);
     MsgStream log(msgSvc, "CgemAliRecTrk");
     log << MSG::DEBUG << "CgemAliRecTrk::setRecTrk()" << endreq;
 
     m_stat = (*it_trk) -> stat();
     m_dr = (*it_trk) -> helix(0);
     m_phi0 = (*it_trk) -> helix(1);
     m_kappa = (*it_trk) -> helix(2);
     m_dz = (*it_trk) -> helix(3);
     m_tanl = (*it_trk) -> helix(4);
 
     m_helix = (*it_trk)->helix();
     m_helixerr = (*it_trk)->err();
 
     m_chisq = (*it_trk) -> chi2();
     m_nhits = (*it_trk) -> getNhits();
     m_ncluster = (*it_trk) -> getNcluster();
 
     if(m_kappa < 0.01){    // no magnetic field
      m_pt = 0.0;
      m_p = 0.0;
     }else{
      m_pt = 1.0 / m_kappa;
      m_p = m_pt * sqrt(m_tanl * m_tanl + 1.0);
     }
 
     ClusterRefVec gothits = (*it_trk) -> getVecClusters();
     ClusterRefVec::iterator it_hit = gothits.begin();
     CgemAliRecHit* rechit;
     for(; it_hit != gothits.end(); it_hit++){
      rechit = new CgemAliRecHit();
      rechit->setRecHit(it_hit);
      m_rechit.push_back(rechit);
     }
}
 
void CgemAliRecTrk::setKalTrk(RecMdcKalTrackCol::iterator it_trk){
     IMessageSvc *msgSvc;
     Gaudi::svcLocator()->service("MessageSvc", msgSvc);
     MsgStream log(msgSvc, "CgemAliRecTrk");
     log << MSG::DEBUG << "CgemAliRecTrk::setKalTrk()" << endreq;
 
     int m_pid = 1;
     if(0 == m_pid) RecMdcKalTrack::setPidType(RecMdcKalTrack::electron);
     else if(1 == m_pid) RecMdcKalTrack::setPidType(RecMdcKalTrack::muon);
     else RecMdcKalTrack::setPidType(RecMdcKalTrack::electron);// pion???
 
     m_stat = 0;
     m_dr = (*it_trk)->dr();
     m_phi0 = (*it_trk)->fi0();
     m_kappa = (*it_trk)->kappa();
     m_dz = (*it_trk)->dz();
     m_tanl = (*it_trk)->tanl();
     m_chisq = (*it_trk)->chi2();
 
     m_dr *= 10.0;      // cm -> mm
     m_dz *= 10.0;      // cm -> mm
 
     m_pt = 1.0 / fabs(m_kappa);
     m_p = m_pt * sqrt( m_tanl * m_tanl + 1.0 );
 
     HelixSegRefVec gothelixsegs = (*it_trk)->getVecHelixSegs();
     HelixSegRefVec::iterator it_hit = gothelixsegs.begin();
     CgemAliRecHit* rechit;
 
     int k = 0;
     m_nhits = k;
}