EventNavigator.cxx

Go to the documentation of this file.

#include <iostream>
 
#include "EventNavigator/EventNavigator.h"
 
using namespace std;
using namespace Event;
 
// methods to retrieve physics objects corresponding given detector id
 
RecEmcShowerVector EventNavigator::getEmcRecShowers( int id ) {
  RecEmcShowerVector                                    tmp;
  multimap<unsigned int, const RecEmcShower*>::iterator it;
  for ( it = m_emcRecShowerId.equal_range( id ).first;
        it != m_emcRecShowerId.equal_range( id ).second; it++ )
    tmp.push_back( ( *it ).second );
  return tmp;
}
 
void EventNavigator::Print() {
  cout << "EventNavigator object:" << std::endl;
  cout << "  Mdc hits:    " << m_mcMdcMcHitsIdx.size() << endl;
  cout << "  Mdc tracks:  " << m_mcMdcTracksIdx.size() << endl;
  cout << "  Emc hits:    " << m_mcEmcMcHitsIdx.size() << endl;
  cout << "  Emc showers: " << m_mcEmcRecShowersIdx.size() << endl;
  cout << "  TABLES:" << endl;
  cout << "   m_mdcTracksWrtMcParticle " << m_mdcTracksWrtMcParticle.size() << endl;
  cout << "   m_mcParticlesWrtMdcTrack " << m_mcParticlesWrtMdcTrack.size() << endl;
  cout << "   m_mdcKalTracksWrtMcParticle " << m_mdcKalTracksWrtMcParticle.size() << endl;
  cout << "   m_mcParticlesWrtMdcKalTrack " << m_mcParticlesWrtMdcKalTrack.size() << endl;
  cout << "   m_mdcTracksWrtMcParticleP " << m_mdcTracksWrtMcParticleP.size() << endl;
  cout << "   m_mcParticlesWrtMdcTrackP " << m_mcParticlesWrtMdcTrackP.size() << endl;
  cout << "   m_mdcKalTracksWrtMcParticle " << m_mdcKalTracksWrtMcParticleP.size() << endl;
  cout << "   m_mcParticlesWrtMdcKalTrack " << m_mcParticlesWrtMdcKalTrackP.size() << endl;
  cout << "   m_mcParticlesWrtRecShower " << m_mcParticlesWrtRecShower.size() << endl;
  cout << "   m_emcRecShowersWrtMcParticle " << m_emcRecShowersWrtMcParticle.size() << endl;
  cout << "   m_mcParticlesWrtEmcShowerP " << m_mcParticlesWrtRecEmcShowerP.size() << endl;
  cout << "   m_emcRecShowersWrtMcParticle " << m_emcRecShowersWrtMcParticleP.size() << endl;
}
 
//=====================================
// MDC
 
McParticleVector& EventNavigator::getMcParticles( const RecMdcTrack* key ) {
  return m_mcParticlesWrtMdcTrack[key->trackId()];
}
 
McParticleVector& EventNavigator::getMcParticles( const RecMdcKalTrack* key ) {
  return m_mcParticlesWrtMdcKalTrack[key->trackId()];
}
 
McParticleVector& EventNavigator::getMcParticles( const RecEmcShower* key ) {
  return m_mcParticlesWrtRecShower[key->getShowerId().get_value()];
}
 
RecMdcTrackVector& EventNavigator::getMdcTracks( const McParticle* key ) {
  return m_mdcTracksWrtMcParticle[key->trackIndex()];
}
 
RecMdcKalTrackVector& EventNavigator::getMdcKalTracks( const McParticle* key ) {
  return m_mdcKalTracksWrtMcParticle[key->trackIndex()];
}
 
// Return number of hits contributed by the McParticle to given MdcTrack
int EventNavigator::getMcParticleRelevance( const RecMdcTrack* key, const McParticle* value ) {
  McParticleVectorP&          vpair = m_mcParticlesWrtMdcTrackP[key->trackId()];
  McParticleVectorP::iterator it;
  for ( it = vpair.begin(); it != vpair.end(); it++ )
  {
    if ( ( *it ).first == value ) { return ( *it ).second; }
  }
 
  std::cerr << __FILE__ << ":" << __LINE__ << " Should not reach here!" << std::endl;
  exit( 1 );
}
 
int EventNavigator::getMcParticleRelevance( const RecMdcKalTrack* key,
                                            const McParticle*     value ) {
  McParticleVectorP&          vpair = m_mcParticlesWrtMdcKalTrackP[key->trackId()];
  McParticleVectorP::iterator it;
  for ( it = vpair.begin(); it != vpair.end(); it++ )
  {
    if ( ( *it ).first == value ) { return ( *it ).second; }
  }
 
  std::cerr << __FILE__ << ":" << __LINE__ << " Should not reach here!" << std::endl;
  exit( 1 );
}
 
int EventNavigator::getMcParticleRelevance( const RecEmcShower* key,
                                            const McParticle*   value ) {
  McParticleVectorP& vpair = m_mcParticlesWrtRecEmcShowerP[key->getShowerId().get_value()];
  McParticleVectorP::iterator it;
  for ( it = vpair.begin(); it != vpair.end(); it++ )
  {
    if ( ( *it ).first == value ) { return ( *it ).second; }
  }
 
  std::cerr << __FILE__ << ":" << __LINE__ << " Should not reach here!" << std::endl;
  exit( 1 );
}
 
RecEmcShowerVector& EventNavigator::getEmcRecShowers( const McParticle* key ) {
  return m_emcRecShowersWrtMcParticle[key->trackIndex()];
}
 
// ==== Methods to store relations ====
// store list of MdcTracks corresponding given McParticle
void EventNavigator::addLink( const McParticle* key, const RecMdcTrack* value ) {
  RecMdcTrackVectorP&          vpair = m_mdcTracksWrtMcParticleP[key->trackIndex()];
  RecMdcTrackVectorP::iterator it;
  bool                         found = false;
  for ( it = vpair.begin(); it != vpair.end(); it++ )
  {
    if ( ( *it ).first == value )
    {
      ( *it ).second++;
      found = true;
      break;
    }
  }
 
  if ( !found ) vpair.push_back( pair<const RecMdcTrack*, int>( value, 1 ) );
 
  RecMdcTrackVector& tmp = m_mdcTracksWrtMcParticle[key->trackIndex()];
  for ( it = vpair.begin(); it != vpair.end(); it++ )
  {
    if ( ( *it ).second > m_mdcCut )
    {
      if ( find( tmp.begin(), tmp.end(), ( *it ).first ) == tmp.end() )
      { tmp.push_back( ( *it ).first ); }
    }
  }
}
 
// store list of MdcKalTracks corresponding given McParticle
void EventNavigator::addLink( const McParticle* key, const RecMdcKalTrack* value ) {
  RecMdcKalTrackVectorP&          vpair = m_mdcKalTracksWrtMcParticleP[key->trackIndex()];
  RecMdcKalTrackVectorP::iterator it;
  bool                            found = false;
  for ( it = vpair.begin(); it != vpair.end(); it++ )
  {
    if ( ( *it ).first == value )
    {
      ( *it ).second++;
      found = true;
      break;
    }
  }
 
  if ( !found ) vpair.push_back( pair<const RecMdcKalTrack*, int>( value, 1 ) );
 
  RecMdcKalTrackVector& tmp = m_mdcKalTracksWrtMcParticle[key->trackIndex()];
  for ( it = vpair.begin(); it != vpair.end(); it++ )
  {
    if ( ( *it ).second > m_mdcCut )
    {
      if ( find( tmp.begin(), tmp.end(), ( *it ).first ) == tmp.end() )
      { tmp.push_back( ( *it ).first ); }
    }
  }
}
 
// store list of McParticles corresponding given MdcTrack
void EventNavigator::addLink( const RecMdcTrack* key, const McParticle* value ) {
  McParticleVectorP&          vpair = m_mcParticlesWrtMdcTrackP[key->trackId()];
  McParticleVectorP::iterator it;
  bool                        found = false;
  for ( it = vpair.begin(); it != vpair.end(); it++ )
  {
    if ( ( *it ).first == value )
    {
      ( *it ).second++;
      found = true;
      break;
    }
  }
 
  if ( !found ) vpair.push_back( pair<const McParticle*, int>( value, 1 ) );
 
  McParticleVector& tmp = m_mcParticlesWrtMdcTrack[key->trackId()];
  for ( it = vpair.begin(); it != vpair.end(); it++ )
  {
    if ( ( *it ).second > m_mdcCut )
    {
      if ( find( tmp.begin(), tmp.end(), ( *it ).first ) == tmp.end() )
      { tmp.push_back( ( *it ).first ); }
    }
  }
}
 
// store list of McParticles corresponding given MdcKalTrack
void EventNavigator::addLink( const RecMdcKalTrack* key, const McParticle* value ) {
  McParticleVectorP&          vpair = m_mcParticlesWrtMdcKalTrackP[key->trackId()];
  McParticleVectorP::iterator it;
  bool                        found = false;
  for ( it = vpair.begin(); it != vpair.end(); it++ )
  {
    if ( ( *it ).first == value )
    {
      ( *it ).second++;
      found = true;
      break;
    }
  }
 
  if ( !found ) vpair.push_back( pair<const McParticle*, int>( value, 1 ) );
 
  McParticleVector& tmp = m_mcParticlesWrtMdcKalTrack[key->trackId()];
  for ( it = vpair.begin(); it != vpair.end(); it++ )
  {
    if ( ( *it ).second > m_mdcCut )
    {
      if ( find( tmp.begin(), tmp.end(), ( *it ).first ) == tmp.end() )
      { tmp.push_back( ( *it ).first ); }
    }
  }
}
 
// store list of McParticles corresponding given EmcRecShower
void EventNavigator::addLink( const RecEmcShower* key, const McParticle* value ) {
  //  McParticleVector& tmp =
  //  m_mcParticlesWrtRecShower[((RecEmcShower*)key)->getShowerId().get_value()]; if
  //  (find(tmp.begin(),tmp.end(),value) == tmp.end() )
  //    tmp.push_back(value);
 
  McParticleVectorP& vpair =
      m_mcParticlesWrtRecEmcShowerP[( (RecEmcShower*)key )->getShowerId().get_value()];
  McParticleVectorP::iterator it;
  bool                        found = false;
  for ( it = vpair.begin(); it != vpair.end(); it++ )
  {
    if ( ( *it ).first == value )
    {
      ( *it ).second++;
      found = true;
      break;
    }
  }
 
  if ( !found ) vpair.push_back( pair<const McParticle*, int>( value, 1 ) );
 
  McParticleVector& tmp =
      m_mcParticlesWrtRecShower[( (RecEmcShower*)key )->getShowerId().get_value()];
  for ( it = vpair.begin(); it != vpair.end(); it++ )
  {
    if ( find( tmp.begin(), tmp.end(), ( *it ).first ) == tmp.end() )
    { tmp.push_back( ( *it ).first ); }
  }
}
 
// store list of EmcRecShowers corresponding given McParticle
void EventNavigator::addLink( const McParticle* key, const RecEmcShower* value ) {
  //  RecEmcShowerVector& tmp = m_emcRecShowersWrtMcParticle[key->trackIndex()];
  //  if (find(tmp.begin(),tmp.end(),value) == tmp.end() )
  //    tmp.push_back(value);
  RecEmcShowerVectorP&          vpair = m_emcRecShowersWrtMcParticleP[key->trackIndex()];
  RecEmcShowerVectorP::iterator it;
  bool                          found = false;
  for ( it = vpair.begin(); it != vpair.end(); it++ )
  {
    if ( ( *it ).first == value )
    {
      ( *it ).second++;
      found = true;
      break;
    }
  }
 
  if ( !found ) vpair.push_back( pair<const RecEmcShower*, int>( value, 1 ) );
 
  RecEmcShowerVector& tmp = m_emcRecShowersWrtMcParticle[key->trackIndex()];
  for ( it = vpair.begin(); it != vpair.end(); it++ )
  {
    if ( find( tmp.begin(), tmp.end(), ( *it ).first ) == tmp.end() )
    { tmp.push_back( ( *it ).first ); }
  }
}