In This Package:

TsESumTriggerTool.cc

Go to the documentation of this file.

#include "TsESumTriggerTool.h"

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

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

#include "Event/SimTrigCommand.h"
#include "Event/SimTrigCommandHeader.h"
#include "Event/SimTrigCommandCollection.h"


TsESumTriggerTool::TsESumTriggerTool(const std::string& type,
                   const std::string& name, 
                   const IInterface* parent)
    : GaudiTool(type,name,parent)
{
  declareInterface< ITsTriggerTool >(this) ;
  // for now add all but RPC's as default.
  m_detectorsToProcess.push_back("DayaBayAD1");
  m_detectorsToProcess.push_back("DayaBayAD2");
  m_detectorsToProcess.push_back("DayaBayIWS");
  m_detectorsToProcess.push_back("DayaBayOWS");
  m_detectorsToProcess.push_back("LingAoAD1");
  m_detectorsToProcess.push_back("LingAoAD2");
  m_detectorsToProcess.push_back("LingAoIWS");
  m_detectorsToProcess.push_back("LingAoOWS"); 
  m_detectorsToProcess.push_back("FarAD1");
  m_detectorsToProcess.push_back("FarAD2");
  m_detectorsToProcess.push_back("FarAD3");
  m_detectorsToProcess.push_back("FarAD4");
  m_detectorsToProcess.push_back("FarIWS");
  m_detectorsToProcess.push_back("FarOWS");
  
  m_UpperLowerThresholdMap["DayaBayAD1"]=DayaBay::Trigger::kADESumThreshold;
  m_UpperLowerThresholdMap["DayaBayAD2"]=DayaBay::Trigger::kADESumThreshold;
  m_UpperLowerThresholdMap["DayaBayIWS"]=DayaBay::Trigger::kNearIWSESumThreshold;
  m_UpperLowerThresholdMap["DayaBayOWS"]=DayaBay::Trigger::kNearOWSESumThreshold;
  m_UpperLowerThresholdMap["LingAoAD1"]=DayaBay::Trigger::kADESumThreshold;
  m_UpperLowerThresholdMap["LingAoAD2"]=DayaBay::Trigger::kADESumThreshold;
  m_UpperLowerThresholdMap["LingAoIWS"]=DayaBay::Trigger::kNearIWSESumThreshold;
  m_UpperLowerThresholdMap["LingAoOWS"]=DayaBay::Trigger::kNearOWSESumThreshold;
  m_UpperLowerThresholdMap["FarAD1"]=DayaBay::Trigger::kADESumThreshold;
  m_UpperLowerThresholdMap["FarAD2"]=DayaBay::Trigger::kADESumThreshold;
  m_UpperLowerThresholdMap["FarAD3"]=DayaBay::Trigger::kADESumThreshold;
  m_UpperLowerThresholdMap["FarAD4"]=DayaBay::Trigger::kADESumThreshold;
  m_UpperLowerThresholdMap["FarIWS"]=DayaBay::Trigger::kFarIWSESumThreshold;
  m_UpperLowerThresholdMap["FarOWS"]=DayaBay::Trigger::kFarOWSESumThreshold;
  
  m_TotalDigitalThresholdMap["DayaBayAD1"]=DayaBay::Trigger::kADESumThreshold;
  m_TotalDigitalThresholdMap["DayaBayAD2"]=DayaBay::Trigger::kADESumThreshold;
  m_TotalDigitalThresholdMap["DayaBayIWS"]=DayaBay::Trigger::kNearIWSESumThreshold;
  m_TotalDigitalThresholdMap["DayaBayOWS"]=DayaBay::Trigger::kNearOWSESumThreshold;
  m_TotalDigitalThresholdMap["LingAoAD1"]=DayaBay::Trigger::kADESumThreshold;
  m_TotalDigitalThresholdMap["LingAoAD2"]=DayaBay::Trigger::kADESumThreshold;
  m_TotalDigitalThresholdMap["LingAoIWS"]=DayaBay::Trigger::kNearIWSESumThreshold;
  m_TotalDigitalThresholdMap["LingAoOWS"]=DayaBay::Trigger::kNearOWSESumThreshold;
  m_TotalDigitalThresholdMap["FarAD1"]=DayaBay::Trigger::kADESumThreshold;
  m_TotalDigitalThresholdMap["FarAD2"]=DayaBay::Trigger::kADESumThreshold;
  m_TotalDigitalThresholdMap["FarAD3"]=DayaBay::Trigger::kADESumThreshold;
  m_TotalDigitalThresholdMap["FarAD4"]=DayaBay::Trigger::kADESumThreshold;
  m_TotalDigitalThresholdMap["FarIWS"]=DayaBay::Trigger::kFarIWSESumThreshold;
  m_TotalDigitalThresholdMap["FarOWS"]=DayaBay::Trigger::kFarOWSESumThreshold;

  
  declareProperty("DetectorsToProcess",m_detectorsToProcess,
                  "List of detectors to process with this tool");
  declareProperty("RecoveryTime",m_recoveryCycles=20,
               "Number of esum (200 MHz) clock cycles for trigger dead time");
  declareProperty("UpperLowerThresholdMap",m_UpperLowerThresholdMap,
                  "Map indexed by detector name with value of Thrshold in volts");
  declareProperty("DigitalTotalThresholdMap",m_TotalDigitalThresholdMap,
                  "Map indexed by detector name with value of Thrshold in volts");
}

TsESumTriggerTool::~TsESumTriggerTool(){}

StatusCode TsESumTriggerTool::initialize()
{
  std::vector<std::string>::iterator it;

  // loop over detectors 
  for(it = m_detectorsToProcess.begin(); 
      it != m_detectorsToProcess.end(); 
      ++it){    
        short int detId =  DayaBay::Detector::siteDetPackedFromString(*it);
        DayaBay::Detector det(detId);
        info() << det 
               << " is setup with Upper Lower threshold " 
               << getULThreshold(det) 
               << endreq;
        info() << det 
               << " is setup with Digital Total threshold " 
               << getDTThreshold(det) 
               << endreq;   
    }
      
    return StatusCode::SUCCESS;
}

StatusCode TsESumTriggerTool::finalize()
{
  return StatusCode::SUCCESS;
}

StatusCode TsESumTriggerTool::fillTriggers(DayaBay::SimTrigHeader* trigHeader,
                                        const DayaBay::ElecHeader& elecHeader)
{
  debug() << "running fillTriggers()" << endreq;
  DayaBay::SimTrigCommandHeader* trigCH 
    = const_cast<DayaBay::SimTrigCommandHeader*>(trigHeader->commandHeader());
  
  const DayaBay::ElecCrateHeader* crateHeader = elecHeader.crateHeader();
  DayaBay::ElecCrateHeader::CrateMap cmap = crateHeader->crates();
  DayaBay::ElecCrateHeader::CrateMap::iterator detIterator;
  
  std::vector<std::string>::iterator it;
  debug() << "processing " << m_detectorsToProcess.size() 
          << " detectors." << endreq;
  
  // loop over detectors 
  for(it = m_detectorsToProcess.begin(); 
      it != m_detectorsToProcess.end(); 
      ++it){
    short int detId =  DayaBay::Detector::siteDetPackedFromString(*it);
    DayaBay::Detector det(detId);
    verbose() << "Processing " << *it << " which is " << det << endreq;
    
    const DayaBay::SimTrigCommandHeader::detCollMap trigCollMap = trigCH->collections();
    

    detIterator = cmap.find(det);
    if(detIterator != cmap.end()){
      // A crate was found for this detector
      debug() << "A crate was found for " << det 
              << " processing..." << endreq;
      
      DayaBay::SimTrigCommandCollection *tcc = 0;
      if( trigCollMap.find(det) == trigCollMap.end() ){
          tcc = new DayaBay::SimTrigCommandCollection();
          tcc->setHeader(trigCH);
          tcc->setDetector(det);
          trigCH->addCollection(tcc);
      }
      
      
      // Get the crate from the header object
      DayaBay::ElecFeeCrate *crate 
                        = dynamic_cast<DayaBay::ElecFeeCrate*>(cmap[det]);
      
      // get esum hits...
      DayaBay::AnalogSignal& upperEsum = crate->esumUpper();
      DayaBay::AnalogSignal& lowerEsum = crate->esumLower();
      DayaBay::AnalogSignal& totalEsum = crate->esumTotal();
      DayaBay::DigitalSignal& totalAdc = crate->esumADC();    
                    
      // threshold for upper and lower esums.
      double ULthreshold = getULThreshold(det);
      
      // threshold for digital and total esums.
      double DTthreshold = getDTThreshold(det);
      
      verbose() << "UL Trigger Threshold " << ULthreshold  << endreq;
      verbose() << "DT Trigger Threshold " << DTthreshold  << endreq;
      
      //
      // IN ALL CASES BELOW THE FOLLOWING APPLIES 
      //
      // convert the esum clock to the next Nhit clock cycle 
      // for the trigger command.  by default this will report the 
      // 200MHz esum clock to a 80 MHz adc clock cycle so it's 
      // properly handeled by the readout.
      // 2*(1 + (previous ADC cycle)) --> trigger reported at next 80 MHz
      //                                 clock cycle
      // The Esum Triggers are checked at 40MHz but ReadoutSim expectes them
      // to be reported at 80MHz, so its calculate at 40 then doubled
      // to 80.
      //
      
      // Do upperEsumTriggers
      
      for( unsigned int i = 0; i < upperEsum.size() ; ++i )
      {
          if ( upperEsum[i] >= ULthreshold )
          {
              DayaBay::Trigger::TriggerType_t type = getTriggerType(det);
              type = DayaBay::Trigger::TriggerType_t(type | DayaBay::Trigger::kESumLow);
              
              unsigned int trigcycle = 2* ( 1 + 
                           (i * DayaBay::AdcCycles) / (DayaBay::EsumCycles) );
              DayaBay::SimTrigCommand *tc = 
                             new DayaBay::SimTrigCommand(det,type,trigcycle,0,totalAdc[i],DayaBay::ESumComp::kESumHigh);
              debug() << "Adding Upper Esum Trigger at cycle\t" <<  i
                      << " With type " << type 
                      << " \n\t Nhit clock cycle: " << trigcycle
                      << " \n\t Esum clock cycle: " << i << endreq;
              verbose() << "Waiting " << m_recoveryCycles 
                        << " cycles to recover" << endreq;           
              i += m_recoveryCycles;
              trigCH->addCommand(tc);
              
          }
      }
      
      // Do lowerEsumTriggers
      
      for( unsigned int i = 0; i < lowerEsum.size() ; ++i )
      {
          if ( lowerEsum[i] >= ULthreshold )
          {
              DayaBay::Trigger::TriggerType_t type = getTriggerType(det);
              type = DayaBay::Trigger::TriggerType_t( DayaBay::Trigger::kESumHigh);
              
              unsigned int trigcycle = 2* ( 1 + 
                           (i * DayaBay::AdcCycles) / (DayaBay::EsumCycles) );
              DayaBay::SimTrigCommand *tc = 
                             new DayaBay::SimTrigCommand(det,type,trigcycle,0,totalAdc[i],DayaBay::ESumComp::kESumLow);
              debug() << "Adding Lower Esum Trigger at cycle\t" <<  i
                      << " With type " << type 
                      << " \n\t Nhit clock cycle: " << trigcycle
                      << " \n\t Esum clock cycle: " << i << endreq;
              verbose() << "Waiting " << m_recoveryCycles 
                        << " cycles to recover" << endreq;           
              i += m_recoveryCycles;
              trigCH->addCommand(tc);
              
          }
      }
      
      // Do TotalEsumTriggers
      for( unsigned int i = 0; i < totalEsum.size() ; ++i )
      {
          if ( totalEsum[i] >= DTthreshold )
          {
              DayaBay::Trigger::TriggerType_t type = getTriggerType(det);
              type = DayaBay::Trigger::TriggerType_t(DayaBay::Trigger::kESumAll);

              unsigned int trigcycle = 2* ( 1 + 
                           (i * DayaBay::AdcCycles) / (DayaBay::EsumCycles) );
              DayaBay::SimTrigCommand *tc = 
                             new DayaBay::SimTrigCommand(det,type,trigcycle,0,totalAdc[i],DayaBay::ESumComp::kESumTotal);
              debug() << "Adding Total Esum Trigger at cycle\t" <<  i
                      << " With type " << type 
                      << " \n\t Nhit clock cycle: " << trigcycle
                      << " \n\t Esum clock cycle: " << i << endreq;
              verbose() << "Waiting " << m_recoveryCycles 
                        << " cycles to recover" << endreq;           
              i += m_recoveryCycles;
              trigCH->addCommand(tc);
              
          }
      }
      
      // Do DigitalEsumTriggers
      for( unsigned int i = 0; i < totalAdc.size() ; ++i )
      {
          if ( totalAdc[i] >= getDigitalValue(DTthreshold) )
          {
              DayaBay::Trigger::TriggerType_t type = getTriggerType(det);
              type = DayaBay::Trigger::TriggerType_t(DayaBay::Trigger::kESumADC);

              unsigned int trigcycle = 2* ( 1 + 
                           (i * DayaBay::AdcCycles) / (DayaBay::EsumCycles) );
              DayaBay::SimTrigCommand *tc = 
                             new DayaBay::SimTrigCommand(det,type,trigcycle,0,totalAdc[i],0);
              debug() << "Adding ADC Esum Trigger at cycle\t" <<  i
                      << " With type " << type 
                      << " \n\t Nhit clock cycle: " << trigcycle
                      << " \n\t Esum clock cycle: " << i << endreq;
              verbose() << "Waiting " << m_recoveryCycles 
                        << " cycles to recover" << endreq;           
              i += m_recoveryCycles;
              trigCH->addCommand(tc);
              
          }
      }
        
    }// end if crate exists
  }
  return StatusCode::SUCCESS;                               
}

// provide mutate to satisfy interface ITsTriggerTool.
StatusCode TsESumTriggerTool::mutate(DayaBay::SimTrigHeader* trigHeader,
                                     const DayaBay::ElecHeader& elecHeader)
{     
  return fillTriggers(trigHeader,elecHeader);
}

// provide interfaces so the defualt way of assigning threshods / 
// types can be transparantly changed.

double TsESumTriggerTool::getULThreshold(DayaBay::Detector det){
    if( m_UpperLowerThresholdMap.find(det.detName()) 
                   != m_UpperLowerThresholdMap.end() ){
        return m_UpperLowerThresholdMap[det.detName()];
    }
    warning() << "Requesting a esum threshold for " << det.detName() 
              << " but none was specified in imput map " << endreq;
    warning() << "Using default value as defined in Conventions/Trigger.h"
              << endreq;
  if(det.detectorId() == DetectorId::kAD1 || 
     det.detectorId() == DetectorId::kAD2 ||
     det.detectorId() == DetectorId::kAD3 || 
     det.detectorId() == DetectorId::kAD4 ){
    return DayaBay::Trigger::kADESumThreshold;
  }
  if(det.detectorId() == DetectorId::kIWS){
    if(det.site() == Site::kFar){
       return DayaBay::Trigger::kFarIWSESumThreshold;
    }else{
       return DayaBay::Trigger::kNearIWSESumThreshold;
    }
  }
  if(det.detectorId() == DetectorId::kOWS){
    if(det.site() == Site::kFar){
       return DayaBay::Trigger::kFarOWSESumThreshold;
    }else{
       return DayaBay::Trigger::kNearOWSESumThreshold;
    }
  }
  return 0.0;
}

double TsESumTriggerTool::getDTThreshold(DayaBay::Detector det){
  if( m_TotalDigitalThresholdMap.find(det.detName()) 
                 != m_TotalDigitalThresholdMap.end() ){
      return m_TotalDigitalThresholdMap[det.detName()];
  }
  
  warning() << "Requesting a esum threshold for " << det.detName() 
            << " but none was specified in imput map " << endreq;
  warning() << "Using default value as defined in Conventions/Trigger.h"
            << endreq;
            
  if(det.detectorId() == DetectorId::kAD1 || 
     det.detectorId() == DetectorId::kAD2 ||
     det.detectorId() == DetectorId::kAD3 || 
     det.detectorId() == DetectorId::kAD4 ){
    return DayaBay::Trigger::kADESumThreshold;
  }
  if(det.detectorId() == DetectorId::kIWS){
    if(det.site() == Site::kFar){
       return DayaBay::Trigger::kFarIWSESumThreshold;
    }else{
       return DayaBay::Trigger::kNearIWSESumThreshold;
    }
  }
  if(det.detectorId() == DetectorId::kOWS){
    if(det.site() == Site::kFar){
       return DayaBay::Trigger::kFarOWSESumThreshold;
    }else{
       return DayaBay::Trigger::kNearOWSESumThreshold;
    }
  }
  return 0.0;
}

DayaBay::Trigger::TriggerType_t TsESumTriggerTool::getTriggerType(DayaBay::Detector det){
  DayaBay::Trigger::TriggerType_t type = DayaBay::Trigger::kUnknown;
  if(det.detectorId() == DetectorId::kAD1) type = DayaBay::Trigger::kAD1;
  if(det.detectorId() == DetectorId::kAD2) type = DayaBay::Trigger::kAD2;
  if(det.detectorId() == DetectorId::kAD3) type = DayaBay::Trigger::kAD3;
  if(det.detectorId() == DetectorId::kAD4) type = DayaBay::Trigger::kAD4;
  if(det.detectorId() == DetectorId::kIWS) type = DayaBay::Trigger::kIWS;
  if(det.detectorId() == DetectorId::kOWS) type = DayaBay::Trigger::kOWS;
  if( type == DayaBay::Trigger::kUnknown) return type;
  return DayaBay::Trigger::TriggerType_t(DayaBay::Trigger::kESum | type); 
}

int TsESumTriggerTool::getDigitalValue(double analogThresh){
    analogThresh=0; // yes compliler I am used, please don't warn me.
    return 4095; // essentially turn it off right now
    //return (unsigned int) analogThresh;
}

---