In This Package:

TsSortTool.cc

Go to the documentation of this file.

#include "TsSortTool.h"

#include "Conventions/Site.h"
#include "Conventions/Trigger.h"

#include "Event/ElecCrateHeader.h"
#include "Event/ElecFeeCrate.h"

#include "Event/SimTrigCommand.h"
#include "Event/SimTrigCommandHeader.h"
#include "Event/SimTrigCommandCollection.h"
#include "CLHEP/Units/SystemOfUnits.h"

TsSortTool::TsSortTool(const std::string& type,
                                             const std::string& name, 
                                             const IInterface* parent)
    : GaudiTool(type,name,parent)
{
  declareInterface< ITsTriggerTool >(this) ;
  declareProperty("RecoveryTime",m_recoveryCycles=60,
                  "Number of 80MHz clock cycles for trigger dead time");
  declareProperty("SortTriggers",m_active=true,
                  "If set to true the tool will sort the trigger commands.");
}

TsSortTool::~TsSortTool(){}

StatusCode TsSortTool::initialize()
{
  debug() << "Initializing TsSortTool" << endreq;
  return StatusCode::SUCCESS;
}

StatusCode TsSortTool::finalize()
{
  debug() << "Finalizing TsSortTool" << endreq;
  return StatusCode::SUCCESS;
}

StatusCode TsSortTool::mutate(DayaBay::SimTrigHeader* trigHeader,
                                         const DayaBay::ElecHeader& elecHeader)
{     
  if ( !m_active ){
      debug() << "TsSortTool is not active, so doing nothing" << endreq;
      return StatusCode::SUCCESS;
  }
  
  debug() << "running mutate()" << endreq;
  const DayaBay::SimTrigCommandHeader *tch = trigHeader->commandHeader();
  const DayaBay::SimTrigCommandHeader::detCollMap trigCollMap 
                                                          = tch->collections();
  DayaBay::SimTrigCommandHeader::detCollMap::const_iterator trigIt;  
  //Loop over All trigger Collections
  for(trigIt  = trigCollMap.begin(); trigIt != trigCollMap.end(); ++trigIt)
  {
    DayaBay::Detector det(trigIt->first);
    debug() << "found " << (trigIt->second)->commands().size()
            << " commands for " << det.detName() << endreq;

    // Skip sorting if less than 2 triggers exist.
    if( (trigIt->second)->commands().size()<2 ) continue;
    DayaBay::SimTrigCommandCollection *triggerColl = 
               const_cast<DayaBay::SimTrigCommandCollection*>(trigIt->second);

    // sort the trigger commands by clock cycle
    triggerColl->sort();

    DayaBay::SimTrigCommandCollection::CommandContainer& detTriggers = 
        const_cast<DayaBay::SimTrigCommandCollection::CommandContainer&>(triggerColl->commands());

    DayaBay::SimTrigCommandCollection::CommandContainer reducedTriggerColl;
    DayaBay::SimTrigCommand* prevTrigger = detTriggers[0];
    reducedTriggerColl.push_back(prevTrigger);
    for(unsigned int trigIdx=1; trigIdx<detTriggers.size(); trigIdx++){
      unsigned int prevcycle = prevTrigger->clockCycle();
      unsigned int thiscycle = (detTriggers[trigIdx])->clockCycle();
      if( (thiscycle-prevcycle)<m_recoveryCycles ){
        // Remove this trigger since it occurs too soon after previous trigger
        verbose() << "Combining trigger at " << thiscycle << " with trigger at " 
                  << prevcycle << endreq;
        // Update the trigger type
        DayaBay::Trigger::TriggerType_t type = (detTriggers[trigIdx])->type();
        int newtype = (int) (type | prevTrigger->type()); // update its type
        prevTrigger->setType( (DayaBay::Trigger::TriggerType_t) newtype );
        delete detTriggers[trigIdx];
        detTriggers[trigIdx]=0;
      }else{
        prevTrigger = detTriggers[trigIdx];
        reducedTriggerColl.push_back(prevTrigger);
      }
    }
    triggerColl->setCommands( reducedTriggerColl );

    verbose() << "Done Sorting. "  << endreq;
  }

  return StatusCode::SUCCESS;                               
}

---