L1GtJetCountsCondition.cc

Go to the documentation of this file.

 
// this class header
#include "L1Trigger/GlobalTrigger/interface/L1GtJetCountsCondition.h"
 
// system include files
#include <iomanip>
#include <iostream>
 
#include <vector>
 
// user include files
//   base classes
#include "L1Trigger/GlobalTrigger/interface/L1GtConditionEvaluation.h"
 
#include "DataFormats/L1GlobalCaloTrigger/interface/L1GctJetCounts.h"
#include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerReadoutSetupFwd.h"
 
#include "CondFormats/L1TObjects/interface/L1GtJetCountsTemplate.h"
 
#include "L1Trigger/GlobalTrigger/interface/L1GlobalTriggerPSB.h"
 
// constructors
//     default
L1GtJetCountsCondition::L1GtJetCountsCondition() : L1GtConditionEvaluation() {
  // empty
}
 
//     from base template condition (from event setup usually)
L1GtJetCountsCondition::L1GtJetCountsCondition(const L1GtCondition *jcTemplate,
                                               const L1GlobalTriggerPSB *ptrPSB,
                                               const int nrL1JetCounts)
    : L1GtConditionEvaluation(),
      m_gtJetCountsTemplate(static_cast<const L1GtJetCountsTemplate *>(jcTemplate)),
      m_gtPSB(ptrPSB),
      m_numberL1JetCounts(nrL1JetCounts) {
  // maximum number of objects received for the evaluation of the condition
  // no objects, in fact, just a number
  m_condMaxNumberObjects = 1;
}
 
// copy constructor
void L1GtJetCountsCondition::copy(const L1GtJetCountsCondition &cp) {
  m_gtJetCountsTemplate = cp.gtJetCountsTemplate();
  m_gtPSB = cp.gtPSB();
 
  m_condMaxNumberObjects = cp.condMaxNumberObjects();
  m_condLastResult = cp.condLastResult();
  m_combinationsInCond = cp.getCombinationsInCond();
 
  m_verbosity = cp.m_verbosity;
}
 
L1GtJetCountsCondition::L1GtJetCountsCondition(const L1GtJetCountsCondition &cp) : L1GtConditionEvaluation() {
  copy(cp);
}
 
// destructor
L1GtJetCountsCondition::~L1GtJetCountsCondition() {
  // empty
}
 
// equal operator
L1GtJetCountsCondition &L1GtJetCountsCondition::operator=(const L1GtJetCountsCondition &cp) {
  copy(cp);
  return *this;
}
 
// methods
void L1GtJetCountsCondition::setGtJetCountsTemplate(const L1GtJetCountsTemplate *jcTemplate) {
  m_gtJetCountsTemplate = jcTemplate;
}
 
void L1GtJetCountsCondition::setGtPSB(const L1GlobalTriggerPSB *ptrPSB) { m_gtPSB = ptrPSB; }
 
// try all object permutations and check spatial correlations, if required
const bool L1GtJetCountsCondition::evaluateCondition() const {
  // number of trigger objects in the condition
  // in fact, there is only one object
  int iCondition = 0;
 
  // condition result condResult will be set to true if the jet counts
  // passes the requirement
  bool condResult = false;
 
  // store the index of the JetCount object
  // from the combination evaluated in the condition
  SingleCombInCond objectsInComb;
 
  // clear the m_combinationsInCond vector
  (combinationsInCond()).clear();
 
  // get the jet counts (event / condition)
  const L1GctJetCounts *jetCounts = m_gtPSB->getCandL1JetCounts();
 
  // protection against missing jet counts collection
  if (jetCounts == nullptr) {
    return false;
  }
 
  const L1GtJetCountsTemplate::ObjectParameter objPar = (*(m_gtJetCountsTemplate->objectParameter()))[iCondition];
 
  unsigned int cIndex = objPar.countIndex;
 
  if (cIndex >= m_numberL1JetCounts) {
    edm::LogError("L1GlobalTrigger") << "\nL1GtJetCountsCondition error: countIndex " << cIndex
                                     << "greater than maximum allowed count = " << m_numberL1JetCounts
                                     << "\n  ==> condResult = false " << std::endl;
    return false;
  }
 
  unsigned int countValue = jetCounts->count(cIndex);
 
  // check countThreshold
  if (!checkThreshold(objPar.countThreshold, countValue, m_gtJetCountsTemplate->condGEq())) {
    return false;
  }
 
  // index is always zero, as they are global quantities (there is only one
  // object)
  int indexObj = 0;
 
  objectsInComb.push_back(indexObj);
  (combinationsInCond()).push_back(objectsInComb);
 
  // if we get here all checks were successful for this combination
  // set the general result for evaluateCondition to "true"
 
  condResult = true;
  return condResult;
}
 
void L1GtJetCountsCondition::print(std::ostream &myCout) const {
  m_gtJetCountsTemplate->print(myCout);
  L1GtConditionEvaluation::print(myCout);
}