/Users/jmonk/Physics/ForIA/src/TriggerDecision.cxx

Go to the documentation of this file.

#include "ForIA/TriggerDecision.hh"

#include "ForIA/Event.hh"

#include <stdexcept>
#include <boost/lexical_cast.hpp>

namespace ForIA {
  
  TriggerDecision::TriggerDecision(): m_ctpBits(0), m_filledBitSet(false), 
  m_efPassedPhysics(0), m_l2PassedPhysics(0), m_nTAVWords(8){}
  
  void TriggerDecision::setEvent(const Event &evt){
    m_trigConf = evt.triggerConfiguration();
    return;
  }
  
  bool TriggerDecision::passesL1(const string &name){
    fillBits();
    
    int pos = m_trigConf.l1Ctpid(name);
    
//    std::cout<<"pos = "<<pos<<std::endl;
//    std::cout<<m_l1Bits.size()<<std::endl;
    
    return m_l1Bits.test(pos);
    
  }
  
  void TriggerDecision::fillBits(){
    
    if(m_filledBitSet) return;
    
    if(m_ctpBits == 0){
      throw std::runtime_error("CTP TAV bits not set before requesting bitset");
    }
    
    if(m_ctpBits->size() != m_nTAVWords){
      throw std::runtime_error("Incorrect number of CTP TAV words: " + 
                               boost::lexical_cast<string>(m_ctpBits->size())+ ", should be "+
                               boost::lexical_cast<string>(m_nTAVWords));
    }
    
    m_l1Bits.reset();
    
    for(vector<unsigned int>::const_reverse_iterator word = m_ctpBits->rbegin();
        word != m_ctpBits->rend(); ++word){
      m_l1Bits <<= 32;
      m_l1Bits |= (*word);
    }
    
    m_filledBitSet = true;
    return;
  }
  
  bool TriggerDecision::passesEF(const string &name)const{
    if(m_efPassed.count(name) != 0) return true;
    if(m_efFailed.count(name) != 0) return false;
    
    if(m_efPassedPhysics == 0){
      throw std::runtime_error("TriggerDecision: Have not set the m_efPassedPhysics data!");
    }
    
    for(vector<short int>::const_iterator passed = m_efPassedPhysics->begin();
        passed != m_efPassedPhysics->end(); ++passed){
      string passedName = m_trigConf.hltChainName(*passed);
      m_efPassed.insert(passedName);
      if(passedName==name)return true;
    }
    
    m_efFailed.insert(name);
    
    return false;
  }
  bool TriggerDecision::passesL2(const string &name){
    if(m_l2Passed.count(name) != 0) return true;
    if(m_l2Failed.count(name) != 0) return false;
    
    if(m_l2PassedPhysics == 0){
      throw std::runtime_error("TriggerDecision: Have not set the m_l2PassedPhysics data!");
    }
    
    for(vector<short int>::const_iterator passed = m_l2PassedPhysics->begin();
        passed != m_l2PassedPhysics->end(); ++passed){
      string passedName = m_trigConf.hltChainName(*passed);
      m_l2Passed.insert(passedName);
      if(passedName==name) return true;
    }
    return false;
  }
  
  
}