MdcPrintSvc Class Reference

#include <MdcPrintSvc.h>

Inheritance diagram for MdcPrintSvc:

Public Member Functions
	MdcPrintSvc (const std::string &name, ISvcLocator *svcloc)
	~MdcPrintSvc ()
StatusCode	initialize ()
StatusCode	finalize ()
void	printRecMdcTrackCol () const
void	printRecMdcTrack (const RecMdcTrack *tk) const
void	printDigi (uint32_t getDigiFlag=0) const
void	printMdcMcHitCol () const

Detailed Description

Definition at line 20 of file MdcPrintSvc.h.

Constructor & Destructor Documentation

MdcPrintSvc()

MdcPrintSvc::MdcPrintSvc	(	const std::string &	name,
		ISvcLocator *	svcloc )

Definition at line 30 of file MdcPrintSvc.cxx.

    : base_class( name, svcloc ) {
 
  // declare properties
  declareProperty( "printHit", m_printHit = 1 );
  declareProperty( "printTrackDetail", m_printTrackDetail = 1 );
  declareProperty( "printErrMat", m_printErrMat = 0 );
  declareProperty( "printTdcAdc", m_printTdcAdc = 0 );
  declareProperty( "printMcTk", m_printMcTk = 0 );
}

Referenced by MdcPrintSvc().

~MdcPrintSvc()

MdcPrintSvc::~MdcPrintSvc

(

)

Definition at line 41 of file MdcPrintSvc.cxx.

{}

Member Function Documentation

finalize()

StatusCode MdcPrintSvc::finalize

(

)

Definition at line 70 of file MdcPrintSvc.cxx.

                                 {
  MsgStream log( msgSvc(), name() );
  log << MSG::INFO << "MdcPrintSvc::finalize()" << endmsg;
 
  return StatusCode::SUCCESS;
}

initialize()

StatusCode MdcPrintSvc::initialize

(

)

Definition at line 43 of file MdcPrintSvc.cxx.

                                   {
  MsgStream log( msgSvc(), name() );
  log << MSG::INFO << "MdcPrintSvc::initialize()" << endmsg;
 
  StatusCode sc = Service::initialize();
  if ( sc.isFailure() ) return sc;
 
  ISvcLocator* svcLocator = Gaudi::svcLocator();
  svcLocator->service( "EventDataSvc", eventSvc, true );
  if ( NULL == eventSvc )
  {
    std::cout << "ERROR Could not load service DataSvc" << std::endl;
    return StatusCode::FAILURE;
  }
 
  IRawDataProviderSvc* irawDataProviderSvc;
  sc                   = service( "RawDataProviderSvc", irawDataProviderSvc );
  m_rawDataProviderSvc = irawDataProviderSvc;
  if ( sc.isFailure() )
  {
    log << MSG::FATAL << "Could not load RawDataProviderSvc!" << endmsg;
    return StatusCode::FAILURE;
  }
 
  return StatusCode::SUCCESS;
}

printDigi()

void MdcPrintSvc::printDigi

(

uint32_t

getDigiFlag = 0

)

const

Definition at line 86 of file MdcPrintSvc.cxx.

                                                        {
  MdcDigiVec           mdcDigiVec = m_rawDataProviderSvc->getMdcDigiVec( getDigiFlag );
  MdcDigiCol::iterator iter       = mdcDigiVec.begin();
  std::cout << "//==== " << name() << " print MdcDigiVec, nDigi=" << mdcDigiVec.size() << " :"
            << std::endl;
  for ( int iDigi = 0; iter != mdcDigiVec.end(); iter++, iDigi++ )
  {
    long l = MdcID::layer( ( *iter )->identify() );
    long w = MdcID::wire( ( *iter )->identify() );
    cout << "(" << l << "," << w;
    if ( m_printMcTk )
    {
      int tkTruth = ( *iter )->getTrackIndex();
      cout << ";" << tkTruth;
    }
    if ( m_printTdcAdc )
    {
      double tdc = RawDataUtil::MdcTime( ( *iter )->getTimeChannel() );
      double adc = RawDataUtil::MdcCharge( ( *iter )->getChargeChannel() );
      cout << ",t " << tdc << ",c " << adc;
    }
    cout << ")";
    if ( iDigi % 4 == 0 ) std::cout << std::endl;
  }
  std::cout << std::endl;
}

printMdcMcHitCol()

void MdcPrintSvc::printMdcMcHitCol

(

)

const

Definition at line 192 of file MdcPrintSvc.cxx.

                                         {
  MsgStream log( msgSvc(), name() );
 
  SmartDataPtr<MdcMcHitCol> mdcMcHitCol( eventSvc, "/Event/MC/MdcMcHitCol" );
  if ( !mdcMcHitCol )
  {
    log << MSG::WARNING << "Could not find MdcMcHitCol" << endmsg;
    return;
  }
 
  int ihit = 0;
  cout << "//==== " << name() << " print MdcMcHitCol nMcHit=" << mdcMcHitCol->size() << endl;
  if ( mdcMcHitCol->size() > 0 )
  { cout << "No. (layer,wire) mcTk  (x,y,z)   lr  driftdist(cm)  energy" << endl; }
  MdcMcHitCol::iterator iter_mchit = mdcMcHitCol->begin();
  for ( ; iter_mchit != mdcMcHitCol->end(); iter_mchit++, ihit++ )
  {
    const Identifier id       = ( *iter_mchit )->identify();
    int              layer    = MdcID::layer( id );
    int              wire     = MdcID::wire( id );
    int              iMcTk    = ( *iter_mchit )->getTrackIndex();
    double           mcX      = ( *iter_mchit )->getPositionX() / 10.;
    double           mcY      = ( *iter_mchit )->getPositionY() / 10.;
    double           mcZ      = ( *iter_mchit )->getPositionZ() / 10.;
    double           mcLR     = ( *iter_mchit )->getPositionFlag();
    double           mcDrift  = ( *iter_mchit )->getDriftDistance() / 10.;
    double           mcEnergy = ( *iter_mchit )->getDepositEnergy() / 10.;
    if ( mcLR == 0 ) mcLR = -1;
    cout << ihit << " (" << layer << "," << wire << ") " << iMcTk << " (" << mcX << "," << mcY
         << "," << mcZ << ")  " << mcLR << "  " << mcDrift << "  " << mcEnergy << endl;
  }
}

printRecMdcTrack()

void MdcPrintSvc::printRecMdcTrack

(

const RecMdcTrack *

tk

)

const

Definition at line 129 of file MdcPrintSvc.cxx.

                                                                {
  std::cout << "//==== " << name() << " print RecMdcTrack No." << tk->trackId() << " :"
            << std::endl;
  cout << " dr " << tk->helix( 0 ) << " phi0 " << tk->helix( 1 ) << " cpa " << tk->helix( 2 )
       << " z0 " << tk->helix( 3 ) << " tanl " << tk->helix( 4 ) << endl;
  if ( m_printTrackDetail )
  {
    std::cout << " q " << tk->charge() << " theta " << tk->theta() << " phi " << tk->phi()
              << " x0 " << tk->x() << " y0 " << tk->y() << " z0 " << tk->z() << " r0 "
              << tk->r() << endl;
    std::cout << " p " << tk->p() << " pt " << tk->pxy() << " pxyz(" << tk->px() << ","
              << tk->py() << "," << tk->pz() << ")" << endl;
    std::cout << " tkStat " << tk->stat() << " chi2 " << tk->chi2() << " ndof " << tk->ndof()
              << " nhit " << tk->getNhits() << " nst " << tk->nster() << endl;
    if ( m_printErrMat )
    {
      std::cout << "errmat   " << std::endl;
      for ( int i = 0; i < 15; i++ ) { std::cout << " " << tk->err( i ); }
      std::cout << "   " << std::endl;
    }
  }
 
  if ( m_printHit )
  {
    int haveDigi[43][288];
    if ( m_printMcTk )
    {
      for ( int ii = 0; ii < 43; ii++ )
      {
        for ( int jj = 0; jj < 43; jj++ ) { haveDigi[ii][jj] = -9999; }
      }
      MdcDigiVec           mdcDigiVec = m_rawDataProviderSvc->getMdcDigiVec();
      MdcDigiCol::iterator iter       = mdcDigiVec.begin();
      std::cout << name() << "//==== " << name()
                << " print MdcDigiVec, nDigi=" << mdcDigiVec.size() << " :" << std::endl;
      for ( int iDigi = 0; iter != mdcDigiVec.end(); iter++, iDigi++ )
      {
        long l         = MdcID::layer( ( *iter )->identify() );
        long w         = MdcID::wire( ( *iter )->identify() );
        haveDigi[l][w] = ( *iter )->getTrackIndex();
      }
    }
 
    int nhits = tk->getVecHits().size();
    std::cout << "nHits=" << nhits << std::endl;
    cout << "No. ";
    if ( m_printMcTk ) cout << "mcTkId ";
    cout << "(layer,wire,lr) stat z" << endl;
    for ( int ii = 0; ii < nhits; ii++ )
    {
      Identifier id( tk->getVecHits()[ii]->getMdcId() );
      int        layer = MdcID::layer( id );
      int        wire  = MdcID::wire( id );
      cout << ii << ":";
      if ( m_printMcTk ) { cout << haveDigi[layer][wire]; }
      cout << " (" << layer << "," << wire << "," << tk->getVecHits()[ii]->getFlagLR() << ")  "
           << tk->getVecHits()[ii]->getStat() << "  " << tk->getVecHits()[ii]->getZhit() << " "
           << std::endl;
    } // end of hit list
    std::cout << "  " << std::endl;
  }
}

Referenced by printRecMdcTrackCol().

printRecMdcTrackCol()

void MdcPrintSvc::printRecMdcTrackCol

(

)

const

Definition at line 113 of file MdcPrintSvc.cxx.

                                            {
  MsgStream log( msgSvc(), name() );
 
  SmartDataPtr<RecMdcTrackCol> recMdcTrackCol( eventSvc, "/Event/Recon/RecMdcTrackCol" );
  if ( !recMdcTrackCol )
  {
    log << MSG::WARNING << "Could not find RecMdcTrackCol" << endmsg;
    return;
  }
 
  std::cout << "//======= " << name() << " RecMdcTrackCol size = " << recMdcTrackCol->size()
            << std::endl;
  RecMdcTrackCol::iterator it = recMdcTrackCol->begin();
  for ( ; it != recMdcTrackCol->end(); it++ ) { printRecMdcTrack( *it ); }
}

---

The documentation for this class was generated from the following files:

• MdcPrintSvc.h
• MdcPrintSvc.cxx