MucMappingAlg.cxx

Go to the documentation of this file.

//------------------------------------------------------------------------------|
//      [File  ]:                       MucMappingAlg.cxx                       |
//      [Brief ]:       Muon detector electronics mapping algrithom for BESIII  |
//      [Author]:       Xie Yuguang, <ygxie@mail.ihep.ac.cn>                    |
//      [Date  ]:       Mar 28, 2006                                            |
//      [Log   ]:       See ChangLog                                            |
//------------------------------------------------------------------------------|
 
#include "GaudiKernel/MsgStream.h"
 
#include "MucMappingAlg/MucMappingAlg.h"
DECLARE_COMPONENT( MucMappingAlg )
// -----------------------------------Declaration----------------------------------------
MucMappingAlg::MucMappingAlg( const std::string& name, ISvcLocator* pSvcLocator )
    : Algorithm( name, pSvcLocator ), m_fFecOrder( 1 ) {
  // Declear the properties
  declareProperty( "FecOrder", m_fFecOrder );
  declareProperty( "ChainMap", m_sFileName[0] );
  declareProperty( "FecMap", m_sFileName[1] );
  declareProperty( "OutputOption", m_sOutputOption );
}
 
// ----------------------------------Initialize-----------------------------------------
StatusCode MucMappingAlg::initialize() {
  MsgStream log( msgSvc(), name() );
  log << MSG::INFO << endmsg << "In initialize()" << endmsg;
  log << MSG::INFO << "Fec order:\t\t" << m_fFecOrder << endmsg;
  log << MSG::INFO << "Chain map:\t" << m_sFileName[0] << endmsg;
  log << MSG::INFO << "Output map:\t\t" << m_sFileName[1] << endmsg;
  log << MSG::INFO << "Output option:\t\t" << m_sOutputOption << endmsg;
 
  ifstream fin( m_sFileName[0].c_str(), ios::in );
 
  log << MSG::INFO << "-------Chain map--------" << endmsg;
  if ( m_sFileName[0].size() == 0 || fin.bad() )
  {
    log << MSG::INFO << endmsg << "Load default map!" << endmsg;
 
    for ( int i = 0; i < MODULE_MAX; i++ )
    {
      for ( int j = 0; j < SOCKET_MAX; j++ )
      {
        m_sChainMap[i][j] = DEFAULT_MAP[i][j];
        // strcpy( m_ChainMap[i][j], DEFAULT_MAP[i][j] );
        log << MSG::INFO << m_sChainMap[i][j] << "\t";
      }
 
      log << MSG::INFO << endmsg;
    }
  }
  else
  {
    string tempMark;
    for ( int i = 0; i < MODULE_MAX; i++ )
    {
      for ( int j = 0; j < SOCKET_MAX; j++ )
      {
        fin >> tempMark;
        // if( strlen( tempMark ) > MARK_NUM )
        if ( tempMark.size() > MARK_NUM )
        {
          log << MSG::INFO << endmsg << "Mark:\t"
              << "Module\t" << i << "Socket\t" << j << "\t" << tempMark << "\terror!"
              << endmsg;
          continue;
        }
 
        m_sChainMap[i][j] = tempMark;
        // strcpy( m_sChainMap[i][j], tempMark );
        log << MSG::INFO << m_sChainMap[i][j] << "\t";
      }
      log << MSG::INFO << endmsg;
    }
 
    fin.close();
  }
  log << MSG::INFO << "------------------------" << endmsg;
 
  return StatusCode::SUCCESS;
}
 
// ----------------------------------Execute--------------------------------------------
StatusCode MucMappingAlg::execute() {
  MsgStream log( msgSvc(), name() );
  log << MSG::INFO << endmsg << "In execute()" << endmsg;
 
  return StatusCode::SUCCESS;
}
 
// ----------------------------------Finalize------------------------------------------
StatusCode MucMappingAlg::finalize() {
  MsgStream log( msgSvc(), name() );
  log << MSG::INFO << endmsg << "In finalize()" << endmsg << endmsg;
 
  ofstream fout( m_sFileName[1].c_str(), ios::out );
  if ( fout.bad() )
  {
    log << MSG::INFO << "Output file created error!" << endmsg;
    exit( -1 );
  }
 
  for ( int i = 0; i < DATA_NUM; i++ )
  {
    if ( m_sOutputOption[i] == '1' ) fout << DATA_NAME[i] << "\t";
  }
  fout << endl;
 
  for ( int i = 0; i < MODULE_MAX; i++ )
  {
    log << MSG::INFO << "----------------------Module [ " << i << " ]----------------------"
        << endmsg;
    for ( int j = 0; j < SOCKET_MAX; j++ )
    {
      int id         = i * SOCKET_MAX + j;
      m_mucChain[id] = new MucChain( id, m_sChainMap[i][j], i, j, m_fFecOrder );
 
      for ( int k = 0; k < m_mucChain[id]->GetFecTotal(); k++ )
      {
        MucFec* aFec = m_mucChain[id]->GetFec( k );
        if ( m_sOutputOption[0] == '1' ) fout << aFec->GetChainName() << "\t";
        if ( m_sOutputOption[1] == '1' ) fout << aFec->DecodeModule() << "\t";
        if ( m_sOutputOption[2] == '1' ) fout << aFec->DecodeSocket() << "\t";
        if ( m_sOutputOption[3] == '1' ) fout << aFec->DecodeFecID() << "\t";
        if ( m_sOutputOption[4] == '1' ) fout << aFec->GetVmeAddress() << "\t";
        if ( m_sOutputOption[5] == '1' ) fout << aFec->GetVmeBin() << "\t";
        if ( m_sOutputOption[6] == '1' ) fout << aFec->GetPart() << "\t";
        if ( m_sOutputOption[7] == '1' ) fout << aFec->GetSegment() << "\t";
        if ( m_sOutputOption[8] == '1' ) fout << aFec->GetLayer() << "\t";
        if ( m_sOutputOption[9] == '1' ) fout << aFec->GetFirstStripID() << "\t";
        if ( m_sOutputOption[10] == '1' ) fout << aFec->GetStripOrder() << "\t";
        if ( m_sOutputOption[11] == '1' ) fout << aFec->GetGeoAddress() << "\t";
        if ( m_sOutputOption[12] == '1' ) fout << aFec->GetGeoBin() << "\t";
        if ( m_sOutputOption[13] == '1' ) fout << aFec->GetDataBin() << "\t";
        fout << endl;
      }
 
      log << MSG::INFO << "Chain ID: " << id << "\tName: " << m_sChainMap[i][j];
      log << MSG::INFO << "\tSocket: " << j << "\tFecTotal: " << m_mucChain[id]->GetFecTotal()
          << "\tFecOrder: " << m_fFecOrder << "\tdone!" << endmsg;
    } // end socket
  }   // end module
 
  fout.close();
 
  return StatusCode::SUCCESS;
}