BesGlobalTrigSvc.cxx

Go to the documentation of this file.

//---------------------------------------------------------------------------//
////      BOOST --- BESIII Object_Oriented Simulation Tool                     /
////
////---------------------------------------------------------------------------/
////
////Description:
////Author:   Caogf
////Created:  Feb, 2006
////Modified:
////Comment:
////
//
#include "GaudiKernel/AlgFactory.h"
#include "GaudiKernel/Bootstrap.h"
#include "GaudiKernel/DataObject.h"
#include "GaudiKernel/IDataProviderSvc.h"
#include "GaudiKernel/IIncidentSvc.h"
#include "GaudiKernel/ISvcLocator.h"
#include "GaudiKernel/Incident.h"
#include "GaudiKernel/MsgStream.h"
#include "GaudiKernel/SmartDataPtr.h"
#include "GaudiKernel/SvcFactory.h"
 
#include "CLHEP/Units/PhysicalConstants.h"
#include "Trigger/BesEmcTrig.h"
#include "Trigger/BesGlobalTrigSvc.h"
#include "Trigger/BesMdcTrig.h"
#include "Trigger/BesMucTrig.h"
#include "Trigger/BesTMTrig.h"
#include "Trigger/BesTofTrig.h"
#include <fstream>
#include <iostream>
#include <string>
using namespace CLHEP;
using namespace std;
 
const std::string BesGlobalTrigSvc::COND_NAME[] = {
    "NClus1",     "NClus2",     "BClusBB",    "EClusBB",  "Clus_Z",     "BClus_PHI",
    "EClus_PHI",  "BEtot_H",    "EEtot_H",    "Etot_L",   "Etot_M",     "BL_Z",
    "NBClus1",    "NEClus1",    "BL_BBLK",    "BL_EBLK",  "ETofBB",     "BTofBB",
    "NETof2",     "NETof1",     "NBTof2",     "NBTof1",   "NTof1",      "CTrk_BB",
    "NCTrk2",     "NCTrk1",     "BTrk_BB",    "NBTrk2",   "NBTrk1",     "ATrk_BB",
    "NATrk2",     "NATrk1",     "MUONBB3478", "MUONBB37", "NotDefined", "NotDefined",
    "NotDefined", "NotDefined", "STrk_BB",    "NStrkN",   "NStrk2",     "NStrk1",
    "LTrk_BB",    "NLtrkN",     "NLtrk2",     "NLtrk1",   "NItrk2",     "NItrk1" };
 
// Instantiation of a static factory class used by clients to create
// instances of this service
// static SvcFactory<BesGlobalTrigSvc> s_factory;
// const ISvcFactory& BesGlobalTrigSvcFactory = s_factory;
 
BesGlobalTrigSvc::BesGlobalTrigSvc( const std::string& name, ISvcLocator* svc )
    : Service( name, svc )
    , m_MdcTrig( 0 )
    , m_TofTrig( 0 )
    , m_EmcTrig( 0 )
    , m_MucTrig( 0 )
    , m_TMTrig( 0 ) {
  trigTableFile = std::string( getenv( "TRIGGERROOT" ) );
  trigTableFile += std::string( "/share/TriggerTable.txt" );
  declareProperty( "TrigTableCon", trigTableFile );
  declareProperty( "TofLayerControl", m_tofLayer = 2 );
  declareProperty( "L1TC_GATE", L1TC_GATE = 70 );
  declareProperty( "L1TC_THRESH", L1TC_THRESH = 130 );
  declareProperty( "L1ETOT_L", L1ETOT_L = 0x56 );
  declareProperty( "L1ETOT_M", L1ETOT_M = 0xba );    // 0xba(run10584),0x190,0x10a(run9571)
  declareProperty( "L1ETOT_BR", L1ETOT_BR = 0x19a ); // 0x216,0x19a
  declareProperty( "L1ETOT_EC", L1ETOT_EC = 0x19a ); // 0x1bd,0x19a
  declareProperty( "L1EBL_BR", L1EBL_BR = 0x30 );
  declareProperty( "L1EBL_EC", L1EBL_EC = 0x1ec );
  declareProperty( "L1EDIFF_BR", L1EDIFF_BR = 0xf6 );
  declareProperty( "L1EDIFF_EC", L1EDIFF_EC = 0xe1 );
  declareProperty( "L1BLK_GATE", L1BLK_GATE = 0x30 );
  declareProperty( "L1EBL_Z", L1EBL_Z = 0xa4 );
  declareProperty( "MAXMDCTRACK", MAXMDCTRACK = 4 );
  declareProperty( "TM_ENDCAP_MODE", tm_endcap = 0 );
  declareProperty( "Energy_Ratio", Energy_Ratio = 1. );
}
BesGlobalTrigSvc::~BesGlobalTrigSvc() {}
//=============================================================================
// Initialize
//=============================================================================
StatusCode BesGlobalTrigSvc::initialize() {
  MsgStream log( msgSvc(), name() );
 
  StatusCode sc = Service::initialize();
  if ( sc.isFailure() ) return sc;
 
  StatusCode status;
  // Get the references to the services that are needed by the ApplicationMgr itself
  IIncidentSvc* incsvc;
  status       = service( "IncidentSvc", incsvc );
  int priority = 100;
  if ( status.isSuccess() ) { incsvc->addListener( this, "NewRun", priority ); }
 
  IRealizationSvc* tmpReal;
  status = service( "RealizationSvc", tmpReal );
  if ( !status.isSuccess() )
  { log << MSG::FATAL << " Could not initialize Realization Service" << endmsg; }
  else { m_RealizationSvc = dynamic_cast<RealizationSvc*>( tmpReal ); }
 
  m_MdcTrig = new BesMdcTrig();
  m_TofTrig = new BesTofTrig();
  m_EmcTrig = new BesEmcTrig();
  m_MucTrig = new BesMucTrig();
  m_TMTrig  = new BesTMTrig();
 
  if ( !m_MdcTrig ) log << MSG::FATAL << "can not get BesMdcTrig pointer " << endmsg;
  if ( !m_TofTrig ) log << MSG::FATAL << "can not get BesTofTrig pointer " << endmsg;
  if ( !m_EmcTrig ) log << MSG::FATAL << "can not get BesEmcTrig pointer " << endmsg;
  if ( !m_MucTrig ) log << MSG::FATAL << "can not get BesMucTrig pointer " << endmsg;
  if ( !m_TMTrig ) log << MSG::FATAL << "can not get BesTMTrig pointer " << endmsg;
 
  trigChannel.clear();
  trigCondName.clear();
  trigTable.clear();
  channelNo.clear();
 
  if ( m_RealizationSvc->UseDBFlag() == false )
  {
    ifstream infile;
    infile.open( trigTableFile.c_str(), ios_base::in );
    if ( !infile )
      log << MSG::FATAL << "can not open Trigger Table file: " << trigTableFile << endmsg;
    char line[255];
    // Skip the header
    do {
      infile.getline( line, 255 );
    } while ( line[0] != 'C' );
    // the name of each trigger channel
    std::string chan;
    int         num;
    char*       token = strtok( line, " " );
    do {
      if ( token )
      {
        chan  = token;
        token = strtok( NULL, " " );
        if ( chan != "CHANNEL" ) trigChannel.push_back( chan );
      }
      else continue;
    } while ( token != NULL );
    // use this trigger channel, yes or no ?
    infile.getline( line, 255 );
    token = strtok( line, " " );
    do {
      if ( token )
      {
        chan  = token;
        token = strtok( NULL, " " );
        num   = atoi( chan.c_str() );
        channelNo.push_back( num );
      }
      else continue;
    } while ( token != NULL );
 
    while ( infile )
    {
      infile.getline( line, 255 );
      if ( line[0] == ' ' ) continue;
      if ( line[0] == '#' ) break;
      token = strtok( line, " " );
      do {
        if ( token )
        {
          chan  = token;
          token = strtok( NULL, " " );
          if ( chan == "1" || chan == "0" || chan == "-" ) { trigTable.push_back( chan ); }
          else { trigCondName.push_back( chan ); }
        }
        else continue;
      } while ( token != NULL );
    }
 
    for ( unsigned int i = 0; i < trigChannel.size(); i++ )
    { cout << " " << trigChannel[i] << " "; }
    cout << endl;
    for ( unsigned int i = 0; i < channelNo.size(); i++ )
    { cout << " " << channelNo[i] << " "; }
    cout << endl;
    for ( unsigned int j = 0; j < trigCondName.size(); j++ )
    {
      cout << trigCondName[j] << " ";
      for ( unsigned int i = 0 + j * trigChannel.size();
            i < trigChannel.size() + j * trigChannel.size(); i++ )
      { cout << trigTable[i] << " "; }
      cout << endl;
    }
    infile.close();
  }
 
  return StatusCode::SUCCESS;
}
 
StatusCode BesGlobalTrigSvc::finalize() {
  MsgStream log( msgSvc(), name() );
  log << MSG::INFO << "Service finalized successfully" << endmsg;
  return StatusCode::SUCCESS;
}
 
//=============================================================================
// QueryInterface
//=============================================================================
StatusCode BesGlobalTrigSvc::queryInterface( const InterfaceID& riid, void** ppvInterface ) {
  if ( IBesGlobalTrigSvc::interfaceID().versionMatch( riid ) )
  { *ppvInterface = (IBesGlobalTrigSvc*)this; }
  else { return Service::queryInterface( riid, ppvInterface ); }
  addRef();
  return StatusCode::SUCCESS;
}
 
void BesGlobalTrigSvc::handle( const Incident& inc ) {
  MsgStream log( msgSvc(), name() );
  log << MSG::DEBUG << "handle: " << inc.type() << endmsg;
  if ( inc.type() != "NewRun" ) { return; }
  log << MSG::DEBUG << "Begin New Run" << endmsg;
  if ( m_RealizationSvc->UseDBFlag() == true && m_RealizationSvc->ifReadTrg() == true )
  {
 
    std::vector<uint32_t> trgTable_DB = m_RealizationSvc->getTrgTable();
 
    trigChannel.clear();
    trigCondName.clear();
    trigTable.clear();
    channelNo.clear();
 
    for ( unsigned int i = 0, mask = 1; i < 16; i++, mask <<= 1 )
    {
      if ( trgTable_DB[0] & mask )
      {
        if ( i == 9 ) channelNo.push_back( 0 );
        else channelNo.push_back( 1 );
      }
      else channelNo.push_back( 0 );
      std::ostringstream osname;
      osname << "Channel_" << i;
      std::string name = osname.str();
      trigChannel.push_back( name );
    }
 
    string table_tmp[16][64];
    for ( unsigned int i = 1; i < 65; i++ )
    {
      for ( unsigned int j = 0; j < 16; j++ )
      {
        int bit2   = ( trgTable_DB[i] >> 2 * j ) & 0x3;
        int row    = (int)( i - 1 ) / 4;
        int column = ( ( i - 1 ) % 4 ) * 16 + j;
        if ( bit2 == 0 ) table_tmp[row][column] = "1";
        else if ( bit2 == 1 ) table_tmp[row][column] = "0";
        else table_tmp[row][column] = "-";
      }
    }
 
    for ( int j = 0; j < 48; j++ )
    {
      trigCondName.push_back( COND_NAME[j] );
      for ( int i = 0; i < 16; i++ ) { trigTable.push_back( table_tmp[i][j] ); }
    }
    // print
    for ( unsigned int i = 0; i < trigChannel.size(); i++ )
    { cout << " " << trigChannel[i] << " "; }
    cout << endl;
    for ( unsigned int i = 0; i < channelNo.size(); i++ )
    { cout << " " << channelNo[i] << " "; }
    cout << endl;
    for ( unsigned int j = 0; j < trigCondName.size(); j++ )
    {
      cout << setw( 15 ) << trigCondName[j] << " ";
      for ( unsigned int i = 0 + j * trigChannel.size();
            i < trigChannel.size() + j * trigChannel.size(); i++ )
      { cout << trigTable[i] << " "; }
      cout << endl;
    }
  }
}
 
int BesGlobalTrigSvc::getL1ETOT_L() {
  L1ETOT_L = m_RealizationSvc->getVthEtotL();
  return L1ETOT_L;
}
 
int BesGlobalTrigSvc::getL1ETOT_M() {
  L1ETOT_M = m_RealizationSvc->getVthEtotM();
  return L1ETOT_M;
}
 
int BesGlobalTrigSvc::getL1ETOT_BR() {
  L1ETOT_BR = m_RealizationSvc->getVthBEtotH();
  return L1ETOT_BR;
}
 
int BesGlobalTrigSvc::getL1ETOT_EC() {
  L1ETOT_EC = m_RealizationSvc->getVthEEtotH();
  return L1ETOT_EC;
}
 
int BesGlobalTrigSvc::getL1EBL_BR() {
  L1EBL_BR = m_RealizationSvc->getVthBalBLK();
  return L1EBL_BR;
}
 
int BesGlobalTrigSvc::getL1EBL_EC() {
  L1EBL_EC = m_RealizationSvc->getVthBalEEMC();
  return L1EBL_EC;
}
 
int BesGlobalTrigSvc::getL1EDIFF_BR() {
  L1EDIFF_BR = m_RealizationSvc->getVthDiffB();
  return L1EDIFF_BR;
}
 
int BesGlobalTrigSvc::getL1EDIFF_EC() {
  L1EDIFF_EC = m_RealizationSvc->getVthDiffE();
  return L1EDIFF_EC;
}
 
int BesGlobalTrigSvc::getL1BLK_GATE() {
  L1BLK_GATE = m_RealizationSvc->getVthBalBLK();
  return L1BLK_GATE;
}
 
int BesGlobalTrigSvc::getL1EBL_Z() {
  L1EBL_Z = m_RealizationSvc->getVthBLZ();
  return L1EBL_Z;
}
 
void BesGlobalTrigSvc::startMdcTrig() {
  if ( m_MdcTrig ) m_MdcTrig->startMdcTrig();
}
 
void BesGlobalTrigSvc::startTofTrig() {
  if ( m_TofTrig ) m_TofTrig->startTofTrig();
}
 
void BesGlobalTrigSvc::startEmcTrig() {
  if ( m_EmcTrig ) m_EmcTrig->startEmcTrig();
}
 
void BesGlobalTrigSvc::startTMTrig() {
  if ( m_TMTrig ) m_TMTrig->startTMTrig();
}
 
StatusCode BesGlobalTrigSvc::GlobalTrig() {
  // initialize
  ifpass = false;
  trigOut.clear();
 
  for ( int i = 0; i < 16; i++ )
  {
    if ( i < 16 ) m_trigChannel[i] = 0;
  }
 
  MsgStream log( msgSvc(), name() );
 
  // get Muc trigger information
  muc_vlayerSeg = m_MucTrig->getNlayerSeg();
 
  muc_nlayerEE  = m_MucTrig->getNlayerPart( 0 );
  muc_nlayerBR  = m_MucTrig->getNlayerPart( 1 );
  muc_nlayerWE  = m_MucTrig->getNlayerPart( 2 );
  muc_vhitLayer = m_MucTrig->getNhitLayer();
  muc_vhitSeg   = m_MucTrig->getNhitSeg();
  muc_nhitEE    = m_MucTrig->getNhitPart( 0 );
  muc_nhitBR    = m_MucTrig->getNhitPart( 1 );
  muc_nhitWE    = m_MucTrig->getNhitPart( 2 );
  muc_nhitTotal = m_MucTrig->getNhitTotal();
 
  // set trigger condition id
  StatusCode status;
  // StatusCode status = setTrigCondition();
 
  // make out trigger conditon table;
  for ( unsigned int i = 0; i < 48; i++ )
  {
    for ( unsigned int j = 0; j < trigChannel.size(); j++ )
    {
      if ( trigCond[i] ) trigOut.push_back( 1 );
      else trigOut.push_back( 0 );
    }
  }
 
  // compare trigger table and trigger output
 
  for ( unsigned int i = 0; i < channelNo.size(); i++ )
  {
    if ( channelNo.size() > 16 )
    {
      log << MSG::FATAL << "Trigger Channel is greater than 16!!!!!" << endmsg;
      return StatusCode::FAILURE;
    }
    int condition   = 0;
    int conditionNo = 0;
    if ( channelNo[i] == 0 )
    {
      m_trigChannel[i] = 0; // data structure member in TDS
      continue;
    }
    else
    {
      log << MSG::INFO << trigChannel[i] << ": ";
      for ( unsigned int j = 0; j < trigCondName.size(); j++ )
      {
        if ( trigTable[i + j * channelNo.size()] == "1" &&
             trigOut[i + j * channelNo.size()] == 1 )
          condition++;
        if ( trigTable[i + j * channelNo.size()] == "0" &&
             trigOut[i + j * channelNo.size()] == 0 )
          condition++;
        if ( trigTable[i + j * channelNo.size()] == "1" )
          log << MSG::INFO << trigCondName[j] << " ";
        if ( trigTable[i + j * channelNo.size()] == "0" )
          log << MSG::INFO << "anti_" << trigCondName[j] << " ";
      }
      for ( unsigned int k = 0; k < trigCondName.size(); k++ )
      {
        if ( trigTable[i + k * channelNo.size()] == "1" ||
             trigTable[i + k * channelNo.size()] == "0" )
        {
          log << MSG::INFO << trigOut[i + k * channelNo.size()] << " ";
          conditionNo++;
        }
      }
      log << MSG::INFO << endmsg;
      if ( ( condition == conditionNo ) && conditionNo != 0 )
      {
        ifpass           = true;
        m_trigChannel[i] = 1; // data structure member in TDS
      }
      else
      {
        m_trigChannel[i] = 0; // data structure member in TDS
      }
    }
  }
 
  return status;
}
 
StatusCode BesGlobalTrigSvc::setTrigCondition() {
  for ( int i = 0; i < 48; i++ )
  {
    trigCond[i]        = false;
    m_trigCondition[i] = 0;
  }
 
  trigCond[0]  = emc_NClus1;
  trigCond[1]  = emc_NClus2;
  trigCond[2]  = emc_BClusBB;
  trigCond[3]  = emc_EClusBB;
  trigCond[4]  = emc_Clus_Z;
  trigCond[5]  = emc_BClus_PHI;
  trigCond[6]  = emc_EClus_PHI;
  trigCond[7]  = emc_BEtot_H;
  trigCond[8]  = emc_EEtot_H;
  trigCond[9]  = emc_Etot_L;
  trigCond[10] = emc_Etot_M;
  trigCond[11] = emc_BL_Z;
  trigCond[12] = emc_NBClus1;
  trigCond[13] = emc_NEClus1;
  trigCond[14] = emc_BL_BBLK;
  trigCond[15] = emc_BL_EBLK;
  trigCond[16] = tof_ETofBB;
  trigCond[17] = tof_BTofBB;
  trigCond[18] = tof_NETof2;
  trigCond[19] = tof_NETof1;
  trigCond[20] = tof_NBTof2;
  trigCond[21] = tof_NBTof1;
  trigCond[22] = tof_NTof1;
  trigCond[23] = tm_CTrk_BB;
  trigCond[24] = tm_NCTrk2;
  trigCond[25] = tm_NCTrk1;
  trigCond[26] = tm_BTrk_BB;
  trigCond[27] = tm_NBTrk2;
  trigCond[28] = tm_NBTrk1;
  trigCond[29] = tm_ATrk_BB;
  trigCond[30] = tm_NATrk2;
  trigCond[31] = tm_NATrk1;
  trigCond[32] = muc_BB3478;
  trigCond[33] = muc_BB37;
  trigCond[34] = false;
  trigCond[35] = false;
  trigCond[36] = false;
  trigCond[37] = false;
  trigCond[38] = mdc_STrk_BB;
  trigCond[39] = mdc_NStrkN;
  trigCond[40] = mdc_NStrk2;
  trigCond[41] = mdc_NStrk1;
  trigCond[42] = mdc_LTrk_BB;
  trigCond[43] = mdc_NLtrkN;
  trigCond[44] = mdc_NLtrk2;
  trigCond[45] = mdc_NLtrk1;
  trigCond[46] = mdc_NItrk2;
  trigCond[47] = mdc_NItrk1;
 
  MsgStream log( msgSvc(), name() );
  if ( trigCondName.size() != 48 )
  {
    log << MSG::FATAL << "The trigger condition CAN NOT be successfully read!!!" << endmsg;
    return StatusCode::FAILURE;
  }
 
  for ( unsigned int i = 0; i < 48; i++ )
  {
    if ( trigCond[i] ) m_trigCondition[i] = 1; // data structure member in TDS
  }
 
  return StatusCode::SUCCESS;
}