InstallArea/include/BesEventMixer/BesEventMixer/MixerAlg.h File Reference

#include "RawDataCnv/EventManagement/RAWEVENT.h"
#include "GaudiKernel/Algorithm.h"
#include "GaudiKernel/SmartDataPtr.h"
#include "RawFile/RawFileReader.h"
#include "RawDataCnv/EventManagement/EmcBuilder.h"
#include "RawDataCnv/EventManagement/MdcBuilder.h"
#include "RawDataCnv/EventManagement/MucBuilder.h"
#include "RawDataCnv/EventManagement/TofBuilder.h"
#include "EmcRawEvent/EmcDigi.h"
#include "MdcRawEvent/MdcDigi.h"
#include "TofRawEvent/TofDigi.h"
#include "MucRawEvent/MucDigi.h"
#include "RawEvent/RawDataUtil.h"
#include "RealizationSvc/RealizationSvc.h"
#include "RealizationSvc/IRealizationSvc.h"
#include "DataInfoSvc/DataInfoSvc.h"
#include "DataInfoSvc/IDataInfoSvc.h"
#include "BesTimerSvc/IBesTimerSvc.h"
#include "BesTimerSvc/BesTimerSvc.h"
#include "GaudiKernel/NTuple.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 30 of file InstallArea/include/BesEventMixer/BesEventMixer/MixerAlg.h.

Function Documentation

combineDigits()

template<class T1 , class T2 >

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

Definition at line 879 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);
}

combineMdcDigits()

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

Definition at line 927 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);
}

combineTofDigits()

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

Definition at line 976 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);
  }
}