TrigEventMaker Class Reference

#include <TrigEventMaker.h>

List of all members.

Public Member Functions
	TrigEventMaker (const std::string &name, ISvcLocator *pSvcLocator)
StatusCode	initialize ()
StatusCode	execute ()
StatusCode	finalize ()

---

Detailed Description

Definition at line 6 of file TrigEventMaker.h.

---

Constructor & Destructor Documentation

TrigEventMaker::TrigEventMaker	(	const std::string &	name,
		ISvcLocator *	pSvcLocator
	)

Definition at line 17 of file TrigEventMaker.cxx.

                                                                              : Algorithm(name, pSvcLocator)
{
  //Declare the properties
}

---

Member Function Documentation

StatusCode TrigEventMaker::execute

(

)

Definition at line 28 of file TrigEventMaker.cxx.

References Bes_Common::DEBUG, Bes_Common::FATAL, genRecEmupikp::i, Bes_Common::INFO, iter(), ganga-rec::j, msgSvc(), delete_small_size::size, and EventModel::Trig::TrigData.

                                  {
  MsgStream log(msgSvc(),name());
  log<<MSG::DEBUG<< "TrigEventMaker: in execute()" <<endreq;

  // Get the event header, print out event and run number
  
  SmartDataPtr<Event::EventHeader> eventHeader(eventSvc(),"/Event/EventHeader");
  if (!eventHeader) {
    log << MSG::FATAL << "Could not find Event Header" << endreq;
    return StatusCode::FAILURE;
  }

  int eventId = eventHeader->eventNumber() ;
  int runId = eventHeader->runNumber();
  log << MSG::INFO << "TrigEventMaker: retrieved event: " <<"  Event Number  "<< eventId << "  run: " <<runId  << endreq;

  //define elements in TrigData structure
  int trgCond[48];
  int trgChan[16];
  int window = 0;
  int timing = 0;
  bool preScale = false;

  //initialize the elements defined above
  for(int i = 0; i < 48; i++) {
    if(i < 16) trgChan[i] = 0;
    trgCond[i] = 0;
  }

  //Get trigger information from TDS
  SmartDataPtr<TrigGTDCol> trigGTDCol(eventSvc(),"/Event/Trig/TrigGTDCol");
  if (!trigGTDCol) {
    log << MSG::FATAL << "Could not find TrigGTDCol!" << endreq;
    return StatusCode::FAILURE;
  }
  
  TrigGTDCol::iterator iter = trigGTDCol->begin();
  for (;iter != trigGTDCol->end(); iter++ ) {
    const uint32_t  boardId = (*iter)->getId();   //The board Id 0xd3: GTL, 0xD2: SAF1, 0xD4: SAF2, 0xD6: SAF3
    const uint32_t  timeWindow = (*iter)->getTimeWindow(); //Time window, bit8 to bit13, total: 0--31
    //const uint32_t  dataType = (*iter)->getDataType(); //bit3 to bit7, total: 0--31
    const uint32_t  size = (*iter)->getDataSize();    //The size of trigger data, not include head
    const uint32_t* trigData = (*iter)->getDataPtr(); //Trigger data

    window = timeWindow;

    //Get data group 5 in GTL, including trigger channel, timing and prescale.
    if(boardId == 0xd3) {
      if(size%timeWindow != 0) {
        log << MSG::FATAL << "GTL data is NOT completed" << endreq;
        return StatusCode::FAILURE;
      }
      for(uint32_t j = 0; j < size; j++) {
        uint32_t dataId = ((trigData[j] >> 24) & 0x7);
        if(dataId != 5) continue; //find data group 5
        for(uint32_t i = 1, loop = 0; loop < 24; i <<= 1, loop++) {
          if(loop < 16) {
            if(trigData[j] & i) trgChan[loop] = 1;
          }
          if((loop == 16) && (trigData[j] & i)) timing = 1;
          if((loop == 17) && (trigData[j] & i) && (timing != 1)) timing = 2;
          if((loop == 18) && (trigData[j] & i) && (timing == 0)) timing = 3;
          if((loop == 21) && (trigData[j] & i)) preScale = true;
        }
      }
    }
    //Get "or 4" in SAF
    if(boardId == 0xd2 || boardId == 0xd4 || boardId == 0xd6) {
      for(uint32_t j = 0; j < size; j++) {
        uint32_t dataId = ((trigData[j] >> 16) & 0xFF);
        if(dataId != 4) continue; //find data "or 4" in SAF
        for(uint32_t i = 1, loop = 0; loop < 16; i <<= 1, loop++) {
          if((boardId == 0xd2) && (trigData[j] & i)) trgCond[32+loop] = 1;
          if((boardId == 0xd4) && (trigData[j] & i)) trgCond[16+loop] = 1;
          if((boardId == 0xd6) && (trigData[j] & i)) trgCond[loop] = 1;
        }
      }
    }
  }

  //Register related trigger information to TDS for physics analysis   
  TrigData* aTrigData = new TrigData(window, timing, trgCond, trgChan, preScale);

  StatusCode sc = StatusCode::SUCCESS ;
  sc = eventSvc()->registerObject(EventModel::Trig::TrigData,aTrigData);
  if(sc!=StatusCode::SUCCESS) {
    log<<MSG::DEBUG<< "Could not register TrigData" <<endreq;
    return StatusCode::FAILURE;
  }
  
  return StatusCode::SUCCESS;
}

StatusCode TrigEventMaker::finalize

(

)

Definition at line 121 of file TrigEventMaker.cxx.

References Bes_Common::DEBUG, and msgSvc().

                                    {
  MsgStream log(msgSvc(), name());
  log << MSG::DEBUG << "==> Finalize TrigEventMaker" << endreq;
  return StatusCode::SUCCESS;
}

StatusCode TrigEventMaker::initialize

(

)

Definition at line 22 of file TrigEventMaker.cxx.

References Bes_Common::DEBUG, and msgSvc().

                                     {
  MsgStream log(msgSvc(), name());
  log << MSG::DEBUG << "TrigEventMaker: in initialize()" << endreq;
  return StatusCode::SUCCESS;
}

---