MixerAlg.cxx File Reference

#include "CLHEP/Random/RandExponential.h"
#include "CLHEP/Random/RandFlat.h"
#include "CLHEP/Random/Random.h"
#include "GaudiKernel/INTuple.h"
#include "GaudiKernel/INTupleSvc.h"
#include "GaudiKernel/MsgStream.h"
#include "IRawFile/RawFileExceptions.h"
#include "MixerAlg.h"
#include "EventModel/Event.h"
#include "EventModel/EventHeader.h"
#include "EventModel/EventModel.h"
#include "HltEvent/DstHltInf.h"
#include "TrigEvent/TrigGTD.h"
#include "BesRndmGenSvc/IBesRndmGenSvc.h"
#include "RawDataCnv/EventManagement/RAWEVENT.h"
#include "RawDataCnv/SniperJSON.h"
#include "RawDataCnv/Util/EmcConverter.h"
#include "RawDataCnv/Util/MdcConverter.h"
#include "RawDataCnv/Util/MucConverter.h"
#include "RawDataCnv/Util/TofConverter.h"
#include "RootCnvSvc/RootInterface.h"
#include <algorithm>
#include <map>
#include <vector>
#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 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().