MixerAlg.cxx File Reference

#include "GaudiKernel/MsgStream.h"
#include "GaudiKernel/INTuple.h"
#include "GaudiKernel/INTupleSvc.h"
#include "BesEventMixer/MixerAlg.h"
#include "RawDataCnv/EventManagement/RAWEVENT.h"
#include "IRawFile/RawFileExceptions.h"
#include "EventModel/EventModel.h"
#include "EventModel/Event.h"
#include "EventModel/EventHeader.h"
#include "TrigEvent/TrigGTD.h"
#include "HltEvent/DstHltInf.h"
#include "BesKernel/IBesRndmGenSvc.h"
#include "CLHEP/Random/Random.h"
#include "CLHEP/Random/RandFlat.h"
#include "CLHEP/Random/RandExponential.h"
#include "RawDataCnv/Util/MdcConverter.h"
#include "RawDataCnv/Util/MucConverter.h"
#include "RawDataCnv/Util/TofConverter.h"
#include "RawDataCnv/Util/EmcConverter.h"
#include "RootCnvSvc/RootInterface.h"
#include <vector>
#include <map>
#include <algorithm>
#include <libgen.h>

Go to the source code of this file.

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)

Function Documentation

combineDigits()

template<class T1 , class T2 >

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

Definition at line 875 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 923 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 972 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().