MucTrigHit.cxx

Go to the documentation of this file.

//---------------------------------------------------------------------------//
////      BOOST --- BESIII Object_Oriented Simulation Tool                     /
////
////---------------------------------------------------------------------------/
////
////Description:
////Author:   Caogf
////Created:  March, 2007
////Modified:
////Comment:
////
//
#include "Trigger/MucTrigHit.h"
 
#include "Identifier/Identifier.h"
#include "Identifier/MucID.h"
#include "RawEvent/DigiEvent.h"
#include "RawEvent/RawDataUtil.h"
#include <math.h>
using namespace std;
 
MucTrigHit* MucTrigHit::muc_Pointer = 0;
 
MucTrigHit* MucTrigHit::get_Muc( void ) {
  if ( !muc_Pointer ) muc_Pointer = new MucTrigHit();
  return muc_Pointer;
}
 
MucTrigHit::MucTrigHit() {}
 
MucTrigHit::~MucTrigHit() {}
 
void MucTrigHit::getMucDigi( MucDigiCol* mucDigiCol ) {
  // reset values
  for ( int i = 0; i < 8; i++ )
    for ( int j = 0; j < 4; j++ )
      for ( int k = 0; k < 112; k++ ) { hitBR[i][j][k] = 0; }
  for ( int i = 0; i < 2; i++ )
    for ( int j = 0; j < 4; j++ )
      for ( int k = 0; k < 4; k++ )
        for ( int l = 0; l < 64; l++ ) { hitEC[i][j][k][l] = 0; }
 
  int                  partId, segId, gapId, stripId;
  MucDigiCol::iterator iDigiCol;
  for ( iDigiCol = mucDigiCol->begin(); iDigiCol != mucDigiCol->end(); iDigiCol++ )
  {
    const Identifier ident = ( *iDigiCol )->identify();
    partId                 = MucID::part( ident );
    segId                  = MucID::seg( ident );
    gapId                  = MucID::gap( ident );
    stripId                = MucID::strip( ident );
    if ( partId == 1 && gapId < 4 ) { hitBR[segId][gapId][stripId] = 1; }
    if ( partId == 0 && gapId < 4 ) hitEC[0][segId][gapId][stripId] = 1; // east endcap
    if ( partId == 2 && gapId < 4 ) hitEC[1][segId][gapId][stripId] = 1; // west endcap
  }
}
 
void MucTrigHit::getFastOr() {
  // reset values
  for ( int i = 0; i < 8; i++ )
    for ( int j = 0; j < 4; j++ )
      for ( int k = 0; k < 7; k++ ) { foBR[i][j][k] = 0; }
  for ( int i = 0; i < 2; i++ )
    for ( int j = 0; j < 4; j++ )
      for ( int k = 0; k < 4; k++ )
        for ( int l = 0; l < 4; l++ ) { foEC[i][j][k][l] = 0; }
  // get fastOr signal
  for ( int i = 0; i < 8; i++ )
    for ( int j = 0; j < 4; j++ )
      for ( int k = 0; k < 112; k++ )
      {
        if ( hitBR[i][j][k] == 1 )
        {
          int nbox         = (int)k / 16;
          foBR[i][j][nbox] = 1;
        }
      }
  for ( int i = 0; i < 2; i++ )
    for ( int j = 0; j < 4; j++ )
      for ( int k = 0; k < 4; k++ )
        for ( int l = 0; l < 64; l++ )
        {
          if ( hitEC[i][j][k][l] == 1 )
          {
            int nbox            = (int)l / 16;
            foEC[i][j][k][nbox] = 1;
          }
        }
}