MixerAlg.h File Reference

#include "GaudiKernel/Algorithm.h"
#include "GaudiKernel/MsgStream.h"
#include "GaudiKernel/NTuple.h"
#include "GaudiKernel/SmartDataPtr.h"
#include "BesRndmGenSvc/IBesRndmGenSvc.h"
#include "BesTimerSvc/IBesTimerSvc.h"
#include "DataInfoSvc/IDataInfoSvc.h"
#include "EmcRawEvent/EmcDigi.h"
#include "MdcRawEvent/MdcDigi.h"
#include "MucRawEvent/MucDigi.h"
#include "RawDataCnv/EventManagement/EmcBuilder.h"
#include "RawDataCnv/EventManagement/MdcBuilder.h"
#include "RawDataCnv/EventManagement/MucBuilder.h"
#include "RawDataCnv/EventManagement/RAWEVENT.h"
#include "RawDataCnv/EventManagement/TofBuilder.h"
#include "RawDataCnv/Util/EmcConverter.h"
#include "RawDataCnv/Util/MdcConverter.h"
#include "RawDataCnv/Util/MucConverter.h"
#include "RawDataCnv/Util/TofConverter.h"
#include "RawEvent/RawDataUtil.h"
#include "RawFile/RawFileReader.h"
#include "RealizationSvc/IRealizationSvc.h"
#include "TofRawEvent/TofDigi.h"

Go to the source code of this file.

Classes
class	MixerAlg

Macros
#define	MAX_LOOP_TIMES   1000000

Functions
template<class T1, class T2>
void	combineDigits (SmartDataPtr< T1 > &mcDigits, T1 &bgDigits, int verbosity)
void	combineMdcDigits (SmartDataPtr< MdcDigiCol > &mcDigits, MdcDigiCol &bgDigits, int verbosity)
void	combineTofDigits (SmartDataPtr< TofDigiCol > &mcDigits, TofDigiCol &bgDigits, int verbosity)

Macro Definition Documentation

MAX_LOOP_TIMES

#define MAX_LOOP_TIMES   1000000

Definition at line 34 of file MixerAlg.h.

Referenced by MixerAlg::execute().

Function Documentation

combineDigits()

template<class T1, class T2>

void combineDigits	(	SmartDataPtr< T1 > &	mcDigits,
		T1 &	bgDigits,
		int	verbosity )

Definition at line 972 of file MixerAlg.cxx.

                                                                              {
  vector<T2*>                 newDigiCol;
  typename T1::iterator       mc;
  typename T1::const_iterator bg;
  bool                        new_digi;
  for ( bg = bgDigits.begin(); bg != bgDigits.end(); bg++ )
  {
    new_digi = true;
    for ( mc = mcDigits->begin(); mc != mcDigits->end(); mc++ )
    {
      if ( ( *mc )->identify() == ( *bg )->identify() )
      {
        if ( verbosity < 2 )
        {
          cout << "****************************************" << endl;
          cout << "MC id " << ( *mc )->identify().get_value() << " BG Id "
               << ( *bg )->identify().get_value() << endl;
          cout << "==> MC Digi : ";
          ( *mc )->fillStream( cout );
          cout << "==> BG Digi : ";
          ( *bg )->fillStream( cout );
        }
 
        ( *mc )->setTrackIndex( ( *mc )->getTrackIndex() - 999 );
        *( *mc ) += *( *bg );
 
        new_digi = false;
        if ( verbosity < 2 )
        {
          cout << "==> New MC Digi: ";
          ( *mc )->fillStream( cout );
          cout << "****************************************" << endl;
        }
      }
    }
 
    // no signal digi in this channel. Create new digi with BG only
    if ( new_digi )
    {
      ( *bg )->setTrackIndex( -1000 );
      newDigiCol.push_back( *bg );
    }
  }
 
  for ( bg = newDigiCol.begin(); bg != newDigiCol.end(); bg++ ) mcDigits->push_back( *bg );
}

Referenced by MixerAlg::mixDigi().

combineMdcDigits()

void combineMdcDigits	(	SmartDataPtr< MdcDigiCol > &	mcDigits,
		MdcDigiCol &	bgDigits,
		int	verbosity )

Definition at line 1019 of file MixerAlg.cxx.

                                       {
  vector<MdcDigi*>           newDigiCol;
  MdcDigiCol::const_iterator mc;
  MdcDigiCol::const_iterator bg;
  bool                       new_digi;
  for ( bg = bgDigits.begin(); bg != bgDigits.end(); bg++ )
  {
    if ( ( *bg )->getChargeChannel() < 0x7FFFFFFF ) ( *bg )->setChargeChannel( 0 );
    new_digi = true;
    for ( mc = mcDigits->begin(); mc != mcDigits->end(); mc++ )
    {
      if ( ( *mc )->identify() == ( *bg )->identify() )
      {
        if ( verbosity < 2 )
        {
          cout << "****************************************" << endl;
          cout << "MC id " << ( *mc )->identify().get_value() << " BG Id "
               << ( *bg )->identify().get_value() << endl;
          cout << "==> MC Digi : ";
          ( *mc )->fillStream( cout );
          cout << "==> BG Digi : ";
          ( *bg )->fillStream( cout );
        }
 
        ( *mc )->setTrackIndex( ( *mc )->getTrackIndex() - 999 );
        *( *mc ) += *( *bg );
 
        new_digi = false;
        if ( verbosity < 2 )
        {
          cout << "==> New MC Digi: ";
          ( *mc )->fillStream( cout );
          cout << "****************************************" << endl;
        }
      }
    }
 
    // no signal digi in this channel. Create new digi with BG only
    if ( new_digi )
    {
      ( *bg )->setTrackIndex( -1000 );
      newDigiCol.push_back( *bg );
    }
  }
 
  for ( bg = newDigiCol.begin(); bg != newDigiCol.end(); bg++ ) mcDigits->push_back( *bg );
}

Referenced by MixerAlg::mixDigi().

combineTofDigits()

void combineTofDigits	(	SmartDataPtr< TofDigiCol > &	mcDigits,
		TofDigiCol &	bgDigits,
		int	verbosity )

Definition at line 1068 of file MixerAlg.cxx.

                                       {
  vector<TofDigi*> newDigiCol;
  // typename TofDigiCol::const_iterator bg;
  TofDigiCol::const_iterator bgTof = bgDigits.begin();
  for ( ; bgTof != bgDigits.end(); bgTof++ )
  {
    ( *bgTof )->setTrackIndex( -1000 );
    newDigiCol.push_back( *bgTof );
  }
  for ( bgTof = newDigiCol.begin(); bgTof != newDigiCol.end(); bgTof++ )
  { mcDigits->push_back( *bgTof ); }
}

Referenced by MixerAlg::mixDigi().