Digi.cxx

Go to the documentation of this file.

#include "PrintMcInfo/PrintMcInfo.h"
using namespace std;
 
void PrintMcInfo::printDigi( ofstream& os ) {
 
  os << '\t' << setw( 40 ) << "MdcDigi" << setw( 30 ) << "TofDigi" << setw( 29 ) << "EmcDigi"
     << "MucDig" << endl;
  os << '\t' << "(layer , wire, trkIndex)"
     << "   "
     << "(B/EC, layer, phi_module, trkIndex)"
     << "  "
     << "(B/EC, theta, phi, trkIndex)"
     << "  "
     << "(B/EC, segment, layer, strip, trkIndex)" << endl;
 
  SmartDataPtr<MdcDigiCol> mdcCol( eventSvc(), "/Event/Digi/MdcDigiCol" );
  SmartDataPtr<TofDigiCol> tofCol( eventSvc(), "/Event/Digi/TofDigiCol" );
  SmartDataPtr<EmcDigiCol> emcCol( eventSvc(), "/Event/Digi/EmcDigiCol" );
  SmartDataPtr<MucDigiCol> mucCol( eventSvc(), "/Event/Digi/MucDigiCol" );
 
  MdcDigiCol::const_iterator it_mdc = mdcCol->begin();
  TofDigiCol::const_iterator it_tof = tofCol->begin();
  EmcDigiCol::const_iterator it_emc = emcCol->begin();
  MucDigiCol::const_iterator it_muc = mucCol->begin();
 
  vector<MdcDigiCol::const_iterator> vmdc;
  vector<TofDigiCol::const_iterator> vtof;
  vector<EmcDigiCol::const_iterator> vemc;
  vector<MucDigiCol::const_iterator> vmuc;
 
  if ( mdcCol->size() != 0 )
  {
    for ( ; it_mdc != mdcCol->end(); it_mdc++ ) { vmdc.push_back( it_mdc ); }
  }
  if ( tofCol->size() != 0 )
  {
    for ( ; it_tof != tofCol->end(); it_tof++ ) { vtof.push_back( it_tof ); }
  }
 
  if ( emcCol->size() != 0 )
  {
    for ( ; it_emc != emcCol->end(); it_emc++ ) { vemc.push_back( it_emc ); }
  }
 
  if ( mucCol->size() != 0 )
  {
    for ( ; it_muc != mucCol->end(); it_muc++ ) { vmuc.push_back( it_muc ); }
  }
 
  for ( int i = 0;; i++ )
  {
    bool bmdc = ( i > vmdc.size() ) ? true : false;
    bool btof = ( i > vtof.size() ) ? true : false;
    bool bemc = ( i > vemc.size() ) ? true : false;
    bool bmuc = ( i > vmuc.size() ) ? true : false;
    if ( ( i >= vmdc.size() ) && ( i >= vtof.size() ) && ( i >= vemc.size() ) &&
         ( i >= vmuc.size() ) )
      break;
    for ( int m = 0; m < 8; m++ ) os << " ";
    if ( vmdc.size() > 0 )
    {
      if ( i < vmdc.size() )
      {
        const Identifier ident_mdc = ( *vmdc[i] )->identify();
        os << setw( 6 ) << MdcID::layer( ident_mdc );
        os << setw( 6 ) << MdcID::wire( ident_mdc );
        os << setw( 15 ) << ( *vmdc[i] )->getTrackIndex();
      }
      else
      {
        for ( int m = 0; m < 27; m++ ) os << " ";
      }
    }
    else
    {
      for ( int m = 0; m < 27; m++ ) os << " ";
    }
 
    if ( vtof.size() > 0 )
    {
      if ( i < vtof.size() )
      {
        const Identifier ident_tof = ( *vtof[i] )->identify();
        if ( TofID::barrel_ec( ident_tof ) == 1 )
        {
          os << "B  "
             << "  ";
        }
        if ( TofID::barrel_ec( ident_tof ) == 0 )
        {
          os << "E_E"
             << "  ";
        }
        if ( TofID::barrel_ec( ident_tof ) == 2 )
        {
          os << "W_E"
             << "  ";
        }
        os << setw( 6 ) << TofID::layer( ident_tof );
        os << setw( 6 ) << TofID::phi_module( ident_tof );
        os << setw( 22 ) << ( *vtof[i] )->getTrackIndex();
      }
      else
      {
        for ( int m = 0; m < 39; m++ ) os << " ";
      }
    }
    else
    {
      for ( int m = 0; m < 39; m++ ) os << " ";
    }
 
    if ( vemc.size() > 0 )
    {
      if ( i < vemc.size() )
      {
        const Identifier ident_emc = ( *vemc[i] )->identify();
        if ( EmcID::barrel_ec( ident_emc ) == 1 )
        {
          os << "B  "
             << "  ";
        }
        if ( EmcID::barrel_ec( ident_emc ) == 0 )
        {
          os << "E_E"
             << "  ";
        }
        if ( EmcID::barrel_ec( ident_emc ) == 2 )
        {
          os << "W_E"
             << "  ";
        }
        os << setw( 6 ) << EmcID::theta_module( ident_emc );
        os << setw( 6 ) << EmcID::phi_module( ident_emc );
        os << setw( 12 ) << ( *vemc[i] )->getTrackIndex();
      }
      else
      {
        for ( int m = 0; m < 29; m++ ) os << " ";
      }
    }
    else
    {
      for ( int m = 0; m < 29; m++ ) os << " ";
    }
 
    if ( vmuc.size() > 0 )
    {
      if ( i < vmuc.size() )
      {
        const Identifier ident_muc = ( *vmuc[i] )->identify();
        if ( MucID::barrel_ec( ident_muc ) == 1 )
        {
          os << "B  "
             << "  ";
        }
        if ( MucID::barrel_ec( ident_muc ) == 0 )
        {
          os << "E_E"
             << "  ";
        }
        if ( MucID::barrel_ec( ident_muc ) == 2 )
        {
          os << "W_E"
             << "  ";
        }
        os << setw( 6 ) << MucID::segment( ident_muc );
        os << setw( 6 ) << MucID::layer( ident_muc );
        os << setw( 6 ) << MucID::strip( ident_muc );
        os << setw( 6 ) << ( *vmuc[i] )->getTrackIndex();
        os << endl;
      }
      else { os << endl; }
    }
    else { os << endl; }
  }
  vmdc.clear();
  vtof.clear();
  vemc.clear();
  vmuc.clear();
}