GenericTriggerEventFlag.cc

Go to the documentation of this file.

//
// $Id: GenericTriggerEventFlag.cc,v 1.13 2012/04/22 15:09:29 vadler Exp $
//

#include "CommonTools/TriggerUtils/interface/GenericTriggerEventFlag.h"

#include "DataFormats/L1GlobalTrigger/interface/L1GtLogicParser.h"

#include <memory>

#include <vector>

#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/Utilities/interface/Transition.h"

// Constants' definitions
static const bool useL1EventSetup(true);
static const bool useL1GtTriggerMenuLite(false);

GenericTriggerEventFlag::GenericTriggerEventFlag(const edm::ParameterSet& config,
                                                 edm::ConsumesCollector& iC,
                                                 l1t::UseEventSetupIn use)
    : GenericTriggerEventFlag(config, iC, false) {
  if (config.exists("andOrL1")) {
    if (stage2_) {
      l1uGt_ = std::make_unique<l1t::L1TGlobalUtil>(config, iC, use);
    }
  }
}

GenericTriggerEventFlag::GenericTriggerEventFlag(const edm::ParameterSet& config,
                                                 edm::ConsumesCollector& iC,
                                                 bool stage1Valid)
    : watchDB_(),
      hltConfigInit_(false),
      andOr_(false),
      dbLabel_(""),
      verbose_(0),
      andOrDcs_(false),
      errorReplyDcs_(false),
      andOrGt_(false),
      gtInputTag_(""),
      gtEvmInputTag_(""),
      gtDBKey_(""),
      errorReplyGt_(false),
      andOrL1_(false),
      stage2_(false),
      l1BeforeMask_(true),
      l1DBKey_(""),
      errorReplyL1_(false),
      andOrHlt_(false),
      hltDBKey_(""),
      errorReplyHlt_(false),
      on_(true),
      onDcs_(true),
      onGt_(true),
      onL1_(true),
      onHlt_(true),
      configError_("CONFIG_ERROR"),
      emptyKeyError_("EMPTY_KEY_ERROR") {
  // General switch(es)
  if (config.exists("andOr")) {
    andOr_ = config.getParameter<bool>("andOr");
    if (config.exists("verbosityLevel"))
      verbose_ = config.getParameter<unsigned>("verbosityLevel");
  } else {
    on_ = false;
    onDcs_ = false;
    onGt_ = false;
    onL1_ = false;
    onHlt_ = false;
  }

  if (on_) {
    if (config.exists("andOrDcs")) {
      andOrDcs_ = config.getParameter<bool>("andOrDcs");
      dcsInputTag_ = config.getParameter<edm::InputTag>("dcsInputTag");
      dcsInputToken_ = iC.mayConsume<DcsStatusCollection>(dcsInputTag_);
      dcsRecordInputTag_ = config.getParameter<edm::InputTag>("dcsRecordInputTag");
      dcsRecordToken_ = iC.mayConsume<DCSRecord>(dcsRecordInputTag_);
      dcsPartitions_ = config.getParameter<std::vector<int>>("dcsPartitions");
      errorReplyDcs_ = config.getParameter<bool>("errorReplyDcs");
    } else {
      onDcs_ = false;
    }
    if (config.exists("andOrGt")) {
      andOrGt_ = config.getParameter<bool>("andOrGt");
      gtInputTag_ = config.getParameter<edm::InputTag>("gtInputTag");
      gtInputToken_ = iC.mayConsume<L1GlobalTriggerReadoutRecord>(gtInputTag_);
      gtLogicalExpressions_ = config.getParameter<std::vector<std::string>>("gtStatusBits");
      errorReplyGt_ = config.getParameter<bool>("errorReplyGt");
      if (config.exists("gtEvmInputTag")) {
        gtEvmInputTag_ = config.getParameter<edm::InputTag>("gtEvmInputTag");
        gtEvmInputToken_ = iC.mayConsume<L1GlobalTriggerEvmReadoutRecord>(gtEvmInputTag_);
      }
      if (config.exists("gtDBKey"))
        gtDBKey_ = config.getParameter<std::string>("gtDBKey");
    } else {
      onGt_ = false;
    }
    if (config.exists("andOrL1")) {
      andOrL1_ = config.getParameter<bool>("andOrL1");
      if (config.exists("stage2"))
        stage2_ = config.getParameter<bool>("stage2");
      else
        stage2_ = false;
      l1LogicalExpressionsCache_ = config.getParameter<std::vector<std::string>>("l1Algorithms");
      errorReplyL1_ = config.getParameter<bool>("errorReplyL1");
      if (config.exists("l1DBKey"))
        l1DBKey_ = config.getParameter<std::string>("l1DBKey");
      if (config.exists("l1BeforeMask"))
        l1BeforeMask_ = config.getParameter<bool>("l1BeforeMask");
    } else {
      onL1_ = false;
    }
    if (config.exists("andOrHlt")) {
      andOrHlt_ = config.getParameter<bool>("andOrHlt");
      hltInputTag_ = config.getParameter<edm::InputTag>("hltInputTag");
      hltInputToken_ = iC.mayConsume<edm::TriggerResults>(hltInputTag_);
      hltLogicalExpressionsCache_ = config.getParameter<std::vector<std::string>>("hltPaths");
      errorReplyHlt_ = config.getParameter<bool>("errorReplyHlt");
      if (config.exists("hltDBKey"))
        hltDBKey_ = config.getParameter<std::string>("hltDBKey");
    } else {
      onHlt_ = false;
    }
    if (!onDcs_ && !onGt_ && !onL1_ && !onHlt_)
      on_ = false;
    else {
      if (config.exists("dbLabel"))
        dbLabel_ = config.getParameter<std::string>("dbLabel");
      watchDB_ = std::make_unique<edm::ESWatcher<AlCaRecoTriggerBitsRcd>>();
    }
  }

  if (onL1_ && !stage2_) {
    l1GtTriggerMenuToken_ = iC.esConsumes<L1GtTriggerMenu, L1GtTriggerMenuRcd, edm::Transition::BeginRun>();
  }
  if ((onGt_ && !gtDBKey_.empty()) || (onL1_ && !l1DBKey_.empty()) || (onHlt_ && !hltDBKey_.empty())) {
    alCaRecoTriggerBitsToken_ = iC.esConsumes<AlCaRecoTriggerBits, AlCaRecoTriggerBitsRcd, edm::Transition::BeginRun>(
        edm::ESInputTag{"", dbLabel_});
  }

  //check to see we arent trying to setup legacy / stage-1 from a constructor call
  //that does not support it
  if (config.exists("andOrL1") && stage2_ == false) {  //stage-1 setup
    if (stage1Valid == false)
      throw cms::Exception("ConfigError") << " Error when constructing GenericTriggerEventFlag, legacy/stage-1 is "
                                             "requested but the constructor called is stage2 only";
  }
}

void GenericTriggerEventFlag::initRun(const edm::Run& run, const edm::EventSetup& setup) {
  if (watchDB_->check(setup)) {
    if (onGt_ && !gtDBKey_.empty()) {
      const std::vector<std::string> exprs(expressionsFromDB(gtDBKey_, setup));
      if (exprs.empty() || exprs.at(0) != configError_)
        gtLogicalExpressions_ = exprs;
    }
    if (onL1_ && !l1DBKey_.empty()) {
      const std::vector<std::string> exprs(expressionsFromDB(l1DBKey_, setup));
      if (exprs.empty() || exprs.at(0) != configError_)
        l1LogicalExpressionsCache_ = exprs;
    }
    if (onHlt_ && !hltDBKey_.empty()) {
      const std::vector<std::string> exprs(expressionsFromDB(hltDBKey_, setup));
      if (exprs.empty() || exprs.at(0) != configError_)
        hltLogicalExpressionsCache_ = exprs;
    }
  }

  // Re-initialise starting valuse before wild-card expansion
  l1LogicalExpressions_ = l1LogicalExpressionsCache_;
  hltLogicalExpressions_ = hltLogicalExpressionsCache_;

  hltConfigInit_ = false;
  if (onHlt_) {
    if (hltInputTag_.process().empty()) {
      if (verbose_ > 0)
        edm::LogError("GenericTriggerEventFlag")
            << "HLT TriggerResults InputTag \"" << hltInputTag_.encode() << "\" specifies no process";
    } else {
      bool hltChanged(false);
      if (!hltConfig_.init(run, setup, hltInputTag_.process(), hltChanged)) {
        if (verbose_ > 0)
          edm::LogError("GenericTriggerEventFlag")
              << "HLT config initialization error with process name \"" << hltInputTag_.process() << "\"";
      } else if (hltConfig_.size() <= 0) {
        if (verbose_ > 0)
          edm::LogError("GenericTriggerEventFlag") << "HLT config size error";
      } else
        hltConfigInit_ = true;
    }
  }

  // Expand version wild-cards in HLT logical expressions
  // L1
  if (onL1_) {
    // build vector of algo names

    std::vector<std::string> algoNames;

    if (stage2_) {
      l1uGt_->retrieveL1Setup(setup);

      const std::vector<std::pair<std::string, double>> prescales = l1uGt_->prescales();
      for (const auto& ip : prescales)
        algoNames.push_back(ip.first);
    } else {
      l1Gt_->getL1GtRunCache(run, setup, useL1EventSetup, useL1GtTriggerMenuLite);
      L1GtTriggerMenu const& l1GtTriggerMenu = setup.get<L1GtTriggerMenuRcd>().get(l1GtTriggerMenuToken_);

      const AlgorithmMap& l1GtPhys(l1GtTriggerMenu.gtAlgorithmMap());
      for (CItAlgo iAlgo = l1GtPhys.begin(); iAlgo != l1GtPhys.end(); ++iAlgo) {
        algoNames.push_back(iAlgo->second.algoName());
      }
      const AlgorithmMap& l1GtTech(l1GtTriggerMenu.gtTechnicalTriggerMap());
      for (CItAlgo iAlgo = l1GtTech.begin(); iAlgo != l1GtTech.end(); ++iAlgo) {
        algoNames.push_back(iAlgo->second.algoName());
      }
    }

    for (unsigned iExpr = 0; iExpr < l1LogicalExpressions_.size(); ++iExpr) {
      std::string l1LogicalExpression(l1LogicalExpressions_.at(iExpr));
      L1GtLogicParser l1AlgoLogicParser(l1LogicalExpression);
      // Loop over algorithms
      for (size_t iAlgo = 0; iAlgo < l1AlgoLogicParser.operandTokenVector().size(); ++iAlgo) {
        const std::string l1AlgoName(l1AlgoLogicParser.operandTokenVector().at(iAlgo).tokenName);
        if (l1AlgoName.find('*') != std::string::npos) {
          l1LogicalExpression.replace(
              l1LogicalExpression.find(l1AlgoName), l1AlgoName.size(), expandLogicalExpression(algoNames, l1AlgoName));
        }
      }
      l1LogicalExpressions_[iExpr] = l1LogicalExpression;
    }
    std::vector<std::string> tmp = l1LogicalExpressions_;
    for (unsigned iExpr = 0; iExpr < tmp.size(); ++iExpr)
      if (std::find(algoNames.begin(), algoNames.end(), tmp[iExpr]) == algoNames.end()) {
        l1LogicalExpressions_.erase(l1LogicalExpressions_.begin() + iExpr);
        if (verbose_ > 1)
          edm::LogWarning("GenericTriggerEventFlag")
              << "L1 algorithm \"" << tmp[iExpr]
              << "\" does not exist in the L1 menu ==> drop it from the list of l1LogicalExpressions";
      }
  }
  // HLT
  if (hltConfigInit_) {
    for (unsigned iExpr = 0; iExpr < hltLogicalExpressions_.size(); ++iExpr) {
      std::string hltLogicalExpression(hltLogicalExpressions_.at(iExpr));
      L1GtLogicParser hltAlgoLogicParser(hltLogicalExpression);
      // Loop over paths
      for (size_t iPath = 0; iPath < hltAlgoLogicParser.operandTokenVector().size(); ++iPath) {
        const std::string hltPathName(hltAlgoLogicParser.operandTokenVector().at(iPath).tokenName);
        if (hltPathName.find('*') != std::string::npos) {
          hltLogicalExpression.replace(hltLogicalExpression.find(hltPathName),
                                       hltPathName.size(),
                                       expandLogicalExpression(hltConfig_.triggerNames(), hltPathName));
        }
      }
      hltLogicalExpressions_[iExpr] = hltLogicalExpression;
    }
  }
}

bool GenericTriggerEventFlag::accept(const edm::Event& event, const edm::EventSetup& setup) {
  if (!on_)
    return true;

  // Determine decision
  if (andOr_)
    return (acceptDcs(event) || acceptGt(event) || acceptL1(event, setup) || acceptHlt(event));
  return (acceptDcs(event) && acceptGt(event) && acceptL1(event, setup) && acceptHlt(event));
}

bool GenericTriggerEventFlag::acceptDcs(const edm::Event& event) {
  // An empty DCS partitions list acts as switch.
  if (!onDcs_ || dcsPartitions_.empty())
    return (!andOr_);  // logically neutral, depending on base logical connective

  bool useDCSRecord(false);

  // Accessing the DcsStatusCollection
  edm::Handle<DcsStatusCollection> dcsStatus;
  event.getByToken(dcsInputToken_, dcsStatus);

  edm::Handle<DCSRecord> dcsRecord;
  event.getByToken(dcsRecordToken_, dcsRecord);

  // none of the DCS products is valid
  if (!dcsStatus.isValid() && !dcsRecord.isValid()) {
    if (verbose_ > 1)
      edm::LogWarning("GenericTriggerEventFlag")
          << "DcsStatusCollection product with InputTag \"" << dcsInputTag_.encode() << "\" not in event \n"
          << "DCSRecord product with InputTag \"" << dcsRecordInputTag_.encode() << "\" not in event \n"
          << " ==> decision: " << errorReplyDcs_;
    return errorReplyDcs_;
  }
  if (dcsStatus.isValid() && (*dcsStatus).empty()) {
    if (event.eventAuxiliary().isRealData()) {
      // this is the Data case for >= Run3, DCSStatus is available (unpacked), but empty
      // becasue SCAL is not in data-taking. In this case we fall back to s/w FED 1022
      if (dcsRecord.isValid()) {
        useDCSRecord = true;
      } else {
        if (verbose_ > 1)
          edm::LogWarning("GenericTriggerEventFlag")
              << "DCSRecord product with InputTag \"" << dcsRecordInputTag_.encode()
              << "\" empty ==> decision: " << errorReplyDcs_;
        return errorReplyDcs_;
      }
    } else {
      // this is the case in which the DCS status is empty, but it's not real data.
      if (verbose_ > 1)
        edm::LogWarning("GenericTriggerEventFlag")
            << "DcsStatusCollection product with InputTag \"" << dcsInputTag_.encode()
            << "\" empty ==> decision: " << errorReplyDcs_;
      return errorReplyDcs_;
    }
  }

  // Determine decision of DCS partition combination and return
  if (andOrDcs_) {  // OR combination
    for (std::vector<int>::const_iterator partitionNumber = dcsPartitions_.begin();
         partitionNumber != dcsPartitions_.end();
         ++partitionNumber) {
      if (acceptDcsPartition(dcsStatus, dcsRecord, useDCSRecord, *partitionNumber))
        return true;
    }
    return false;
  }
  for (std::vector<int>::const_iterator partitionNumber = dcsPartitions_.begin();
       partitionNumber != dcsPartitions_.end();
       ++partitionNumber) {
    if (!acceptDcsPartition(dcsStatus, dcsRecord, useDCSRecord, *partitionNumber))
      return false;
  }
  return true;
}

bool GenericTriggerEventFlag::acceptDcsPartition(const edm::Handle<DcsStatusCollection>& dcsStatus,
                                                 const edm::Handle<DCSRecord>& dcsRecord,
                                                 bool useDCSRecord,
                                                 int dcsPartition) const {
  int theDCSRecordPartition;
  // Error checks
  switch (dcsPartition) {
    case DcsStatus::EBp:
      theDCSRecordPartition = DCSRecord::EBp;
      break;
    case DcsStatus::EBm:
      theDCSRecordPartition = DCSRecord::EBm;
      break;
    case DcsStatus::EEp:
      theDCSRecordPartition = DCSRecord::EEp;
      break;
    case DcsStatus::EEm:
      theDCSRecordPartition = DCSRecord::EBm;
      break;
    case DcsStatus::HBHEa:
      theDCSRecordPartition = DCSRecord::HBHEa;
      break;
    case DcsStatus::HBHEb:
      theDCSRecordPartition = DCSRecord::HBHEb;
      break;
    case DcsStatus::HBHEc:
      theDCSRecordPartition = DCSRecord::HBHEc;
      break;
    case DcsStatus::HF:
      theDCSRecordPartition = DCSRecord::HF;
      break;
    case DcsStatus::HO:
      theDCSRecordPartition = DCSRecord::HO;
      break;
    case DcsStatus::RPC:
      theDCSRecordPartition = DCSRecord::RPC;
      break;
    case DcsStatus::DT0:
      theDCSRecordPartition = DCSRecord::DT0;
      break;
    case DcsStatus::DTp:
      theDCSRecordPartition = DCSRecord::DTp;
      break;
    case DcsStatus::DTm:
      theDCSRecordPartition = DCSRecord::DTm;
      break;
    case DcsStatus::CSCp:
      theDCSRecordPartition = DCSRecord::CSCp;
      break;
    case DcsStatus::CSCm:
      theDCSRecordPartition = DCSRecord::CSCm;
      break;
    case DcsStatus::CASTOR:
      theDCSRecordPartition = DCSRecord::CASTOR;
      break;
    case DcsStatus::TIBTID:
      theDCSRecordPartition = DCSRecord::TIBTID;
      break;
    case DcsStatus::TOB:
      theDCSRecordPartition = DCSRecord::TOB;
      break;
    case DcsStatus::TECp:
      theDCSRecordPartition = DCSRecord::TECp;
      break;
    case DcsStatus::TECm:
      theDCSRecordPartition = DCSRecord::TECm;
      break;
    case DcsStatus::BPIX:
      theDCSRecordPartition = DCSRecord::BPIX;
      break;
    case DcsStatus::FPIX:
      theDCSRecordPartition = DCSRecord::FPIX;
      break;
    case DcsStatus::ESp:
      theDCSRecordPartition = DCSRecord::ESp;
      break;
    case DcsStatus::ESm:
      theDCSRecordPartition = DCSRecord::ESm;
      break;
    default:
      if (verbose_ > 1)
        edm::LogWarning("GenericTriggerEventFlag")
            << "DCS partition number \"" << dcsPartition << "\" does not exist ==> decision: " << errorReplyDcs_;
      return errorReplyDcs_;
  }

  // Determine decision
  if (!useDCSRecord) {
    return dcsStatus->at(0).ready(dcsPartition);
  } else {
    if (verbose_ > 2) {
      LogDebug("GenericTriggerEventFlag")
          << "using dcs record, dcsPartition:" << dcsPartition << " " << theDCSRecordPartition << " "
          << (*dcsRecord).partitionName(theDCSRecordPartition) << " "
          << (*dcsRecord).highVoltageReady(theDCSRecordPartition) << std::endl;
    }
    return (*dcsRecord).highVoltageReady(theDCSRecordPartition);
  }
}

bool GenericTriggerEventFlag::acceptGt(const edm::Event& event) {
  // An empty GT status bits logical expressions list acts as switch.
  if (!onGt_ || gtLogicalExpressions_.empty())
    return (!andOr_);  // logically neutral, depending on base logical connective

  // Determine decision of GT status bits logical expression combination and return
  if (andOrGt_) {  // OR combination
    for (std::vector<std::string>::const_iterator gtLogicalExpression = gtLogicalExpressions_.begin();
         gtLogicalExpression != gtLogicalExpressions_.end();
         ++gtLogicalExpression) {
      if (acceptGtLogicalExpression(event, *gtLogicalExpression))
        return true;
    }
    return false;
  }
  for (std::vector<std::string>::const_iterator gtLogicalExpression = gtLogicalExpressions_.begin();
       gtLogicalExpression != gtLogicalExpressions_.end();
       ++gtLogicalExpression) {
    if (!acceptGtLogicalExpression(event, *gtLogicalExpression))
      return false;
  }
  return true;
}

bool GenericTriggerEventFlag::acceptGtLogicalExpression(const edm::Event& event, std::string gtLogicalExpression) {
  // Check empty std::strings
  if (gtLogicalExpression.empty()) {
    if (verbose_ > 1)
      edm::LogWarning("GenericTriggerEventFlag") << "Empty logical expression ==> decision: " << errorReplyGt_;
    return errorReplyGt_;
  }

  // Negated paths
  bool negExpr(negate(gtLogicalExpression));
  if (negExpr && gtLogicalExpression.empty()) {
    if (verbose_ > 1)
      edm::LogWarning("GenericTriggerEventFlag")
          << "Empty (negated) logical expression ==> decision: " << errorReplyGt_;
    return errorReplyGt_;
  }

  // Parse logical expression and determine GT status bit decision
  L1GtLogicParser gtAlgoLogicParser(gtLogicalExpression);
  // Loop over status bits
  for (size_t iStatusBit = 0; iStatusBit < gtAlgoLogicParser.operandTokenVector().size(); ++iStatusBit) {
    const std::string gtStatusBit(gtAlgoLogicParser.operandTokenVector().at(iStatusBit).tokenName);
    // Manipulate status bit decision as stored in the parser
    bool decision(errorReplyDcs_);
    // Hard-coded status bits!!!
    if (gtStatusBit == "PhysDecl" || gtStatusBit == "PhysicsDeclared") {
      edm::Handle<L1GlobalTriggerReadoutRecord> gtReadoutRecord;
      event.getByToken(gtInputToken_, gtReadoutRecord);
      if (!gtReadoutRecord.isValid()) {
        if (verbose_ > 1)
          edm::LogWarning("GenericTriggerEventFlag")
              << "L1GlobalTriggerReadoutRecord product with InputTag \"" << gtInputTag_.encode()
              << "\" not in event ==> decision: " << errorReplyGt_;
        gtAlgoLogicParser.operandTokenVector().at(iStatusBit).tokenResult = errorReplyDcs_;
        continue;
      }
      decision = (gtReadoutRecord->gtFdlWord().physicsDeclared() == 1);
    } else if (gtStatusBit == "Stable" || gtStatusBit == "StableBeam" || gtStatusBit == "Adjust" ||
               gtStatusBit == "Sqeeze" || gtStatusBit == "Flat" || gtStatusBit == "FlatTop" || gtStatusBit == "7TeV" ||
               gtStatusBit == "8TeV" || gtStatusBit == "13TeV" || gtStatusBit == "2360GeV" || gtStatusBit == "900GeV") {
      edm::Handle<L1GlobalTriggerEvmReadoutRecord> gtEvmReadoutRecord;
      event.getByToken(gtEvmInputToken_, gtEvmReadoutRecord);
      if (!gtEvmReadoutRecord.isValid()) {
        if (verbose_ > 1)
          edm::LogWarning("GenericTriggerEventFlag")
              << "L1GlobalTriggerEvmReadoutRecord product with InputTag \"" << gtEvmInputTag_.encode()
              << "\" not in event ==> decision: " << errorReplyGt_;
        gtAlgoLogicParser.operandTokenVector().at(iStatusBit).tokenResult = errorReplyDcs_;
        continue;
      }
      if (gtStatusBit == "Stable" || gtStatusBit == "StableBeam") {
        decision = (gtEvmReadoutRecord->gtfeWord().beamMode() == 11);
      } else if (gtStatusBit == "Adjust") {
        decision = (10 <= gtEvmReadoutRecord->gtfeWord().beamMode() && gtEvmReadoutRecord->gtfeWord().beamMode() <= 11);
      } else if (gtStatusBit == "Sqeeze") {
        decision = (9 <= gtEvmReadoutRecord->gtfeWord().beamMode() && gtEvmReadoutRecord->gtfeWord().beamMode() <= 11);
      } else if (gtStatusBit == "Flat" || gtStatusBit == "FlatTop") {
        decision = (8 <= gtEvmReadoutRecord->gtfeWord().beamMode() && gtEvmReadoutRecord->gtfeWord().beamMode() <= 11);
      } else if (gtStatusBit == "7TeV") {
        decision = (gtEvmReadoutRecord->gtfeWord().beamMomentum() == 3500);
      } else if (gtStatusBit == "8TeV") {
        decision = (gtEvmReadoutRecord->gtfeWord().beamMomentum() == 4000);
      } else if (gtStatusBit == "13TeV") {
        decision = (gtEvmReadoutRecord->gtfeWord().beamMomentum() == 6500);
      } else if (gtStatusBit == "2360GeV") {
        decision = (gtEvmReadoutRecord->gtfeWord().beamMomentum() == 1180);
      } else if (gtStatusBit == "900GeV") {
        decision = (gtEvmReadoutRecord->gtfeWord().beamMomentum() == 450);
      }
    } else {
      if (verbose_ > 1)
        edm::LogWarning("GenericTriggerEventFlag")
            << "GT status bit \"" << gtStatusBit << "\" is not defined ==> decision: " << errorReplyGt_;
    }
    gtAlgoLogicParser.operandTokenVector().at(iStatusBit).tokenResult = decision;
  }

  // Determine decision
  const bool gtDecision(gtAlgoLogicParser.expressionResult());
  return negExpr ? (!gtDecision) : gtDecision;
}

bool GenericTriggerEventFlag::acceptL1(const edm::Event& event, const edm::EventSetup& setup) {
  // An empty L1 logical expressions list acts as switch.
  if (!onL1_ || l1LogicalExpressions_.empty())
    return (!andOr_);  // logically neutral, depending on base logical connective

  // Determine decision of L1 logical expression combination and return
  if (andOrL1_) {  // OR combination
    for (std::vector<std::string>::const_iterator l1LogicalExpression = l1LogicalExpressions_.begin();
         l1LogicalExpression != l1LogicalExpressions_.end();
         ++l1LogicalExpression) {
      if (acceptL1LogicalExpression(event, setup, *l1LogicalExpression))
        return true;
    }
    return false;
  }
  for (std::vector<std::string>::const_iterator l1LogicalExpression = l1LogicalExpressions_.begin();
       l1LogicalExpression != l1LogicalExpressions_.end();
       ++l1LogicalExpression) {
    if (!acceptL1LogicalExpression(event, setup, *l1LogicalExpression))
      return false;
  }
  return true;
}

bool GenericTriggerEventFlag::acceptL1LogicalExpression(const edm::Event& event,
                                                        const edm::EventSetup& setup,
                                                        std::string l1LogicalExpression) {
  // Getting the L1 event setup
  if (stage2_)
    l1uGt_->retrieveL1(event, setup);
  else
    // FIXME This can possibly go to initRun()
    l1Gt_->getL1GtRunCache(event, setup, useL1EventSetup, useL1GtTriggerMenuLite);

  // Check empty std::strings
  if (l1LogicalExpression.empty()) {
    if (verbose_ > 1)
      edm::LogWarning("GenericTriggerEventFlag") << "Empty logical expression ==> decision: " << errorReplyL1_;
    return errorReplyL1_;
  }

  // Negated logical expression
  bool negExpr(negate(l1LogicalExpression));
  if (negExpr && l1LogicalExpression.empty()) {
    if (verbose_ > 1)
      edm::LogWarning("GenericTriggerEventFlag")
          << "Empty (negated) logical expression ==> decision: " << errorReplyL1_;
    return errorReplyL1_;
  }

  // Parse logical expression and determine L1 decision
  L1GtLogicParser l1AlgoLogicParser(l1LogicalExpression);
  // Loop over algorithms
  for (size_t iAlgorithm = 0; iAlgorithm < l1AlgoLogicParser.operandTokenVector().size(); ++iAlgorithm) {
    const std::string l1AlgoName(l1AlgoLogicParser.operandTokenVector().at(iAlgorithm).tokenName);

    bool decision = false;
    bool error = false;
    if (stage2_) {
      bool errorBOOL = (l1BeforeMask_ ? l1uGt_->getInitialDecisionByName(l1AlgoName, decision)
                                      : l1uGt_->getFinalDecisionByName(l1AlgoName, decision));
      error = !errorBOOL;
    } else {
      int errorINT(-1);
      //      const bool decision( l1BeforeMask_ ? l1Gt_->decisionBeforeMask( event, l1AlgoName, errorINT ) : l1Gt_->decisionAfterMask( event, l1AlgoName, errorINT ) );
      decision = (l1BeforeMask_ ? l1Gt_->decisionBeforeMask(event, l1AlgoName, errorINT)
                                : l1Gt_->decisionAfterMask(event, l1AlgoName, errorINT));
      error = (errorINT != 0);
      if (errorINT > 1)
        if (verbose_ > 1)
          edm::LogWarning("GenericTriggerEventFlag")
              << "L1 algorithm \"" << l1AlgoName << "\" received error code " << error
              << " from L1GtUtils::decisionBeforeMask ==> decision: " << errorReplyL1_;
    }

    // Error checks
    if (error) {
      if (verbose_ > 1)
        edm::LogWarning("GenericTriggerEventFlag")
            << "L1 algorithm \"" << l1AlgoName << "\" does not exist in the L1 menu ==> decision: " << errorReplyL1_;
      l1AlgoLogicParser.operandTokenVector().at(iAlgorithm).tokenResult = errorReplyL1_;
      continue;
    }
    // Manipulate algo decision as stored in the parser
    l1AlgoLogicParser.operandTokenVector().at(iAlgorithm).tokenResult = decision;
  }

  // Return decision
  const bool l1Decision(l1AlgoLogicParser.expressionResult());
  return negExpr ? (!l1Decision) : l1Decision;
}

bool GenericTriggerEventFlag::acceptHlt(const edm::Event& event) {
  // An empty HLT logical expressions list acts as switch.
  if (!onHlt_ || hltLogicalExpressions_.empty())
    return (!andOr_);  // logically neutral, depending on base logical connective

  // Checking the HLT configuration,
  if (!hltConfigInit_) {
    if (verbose_ > 1)
      edm::LogWarning("GenericTriggerEventFlag") << "HLT config error ==> decision: " << errorReplyHlt_;
    return errorReplyHlt_;
  }

  // Accessing the TriggerResults
  edm::Handle<edm::TriggerResults> hltTriggerResults;
  event.getByToken(hltInputToken_, hltTriggerResults);
  if (!hltTriggerResults.isValid()) {
    if (verbose_ > 1)
      edm::LogWarning("GenericTriggerEventFlag") << "TriggerResults product with InputTag \"" << hltInputTag_.encode()
                                                 << "\" not in event ==> decision: " << errorReplyHlt_;
    return errorReplyHlt_;
  }
  if ((*hltTriggerResults).size() == 0) {
    if (verbose_ > 1)
      edm::LogWarning("GenericTriggerEventFlag") << "TriggerResults product with InputTag \"" << hltInputTag_.encode()
                                                 << "\" empty ==> decision: " << errorReplyHlt_;
    return errorReplyDcs_;
  }

  // Determine decision of HLT logical expression combination and return
  if (andOrHlt_) {  // OR combination
    for (std::vector<std::string>::const_iterator hltLogicalExpression = hltLogicalExpressions_.begin();
         hltLogicalExpression != hltLogicalExpressions_.end();
         ++hltLogicalExpression) {
      if (acceptHltLogicalExpression(hltTriggerResults, *hltLogicalExpression))
        return true;
    }
    return false;
  }
  for (std::vector<std::string>::const_iterator hltLogicalExpression = hltLogicalExpressions_.begin();
       hltLogicalExpression != hltLogicalExpressions_.end();
       ++hltLogicalExpression) {
    if (!acceptHltLogicalExpression(hltTriggerResults, *hltLogicalExpression))
      return false;
  }
  return true;
}

bool GenericTriggerEventFlag::acceptHltLogicalExpression(const edm::Handle<edm::TriggerResults>& hltTriggerResults,
                                                         std::string hltLogicalExpression) const {
  // Check empty std::strings
  if (hltLogicalExpression.empty()) {
    if (verbose_ > 1)
      edm::LogWarning("GenericTriggerEventFlag") << "Empty logical expression ==> decision: " << errorReplyHlt_;
    return errorReplyHlt_;
  }

  // Negated paths
  bool negExpr(negate(hltLogicalExpression));
  if (negExpr && hltLogicalExpression.empty()) {
    if (verbose_ > 1)
      edm::LogWarning("GenericTriggerEventFlag")
          << "Empty (negated) logical expression ==> decision: " << errorReplyHlt_;
    return errorReplyHlt_;
  }

  // Parse logical expression and determine HLT decision
  L1GtLogicParser hltAlgoLogicParser(hltLogicalExpression);
  // Loop over paths
  for (size_t iPath = 0; iPath < hltAlgoLogicParser.operandTokenVector().size(); ++iPath) {
    const std::string hltPathName(hltAlgoLogicParser.operandTokenVector().at(iPath).tokenName);
    const unsigned indexPath(hltConfig_.triggerIndex(hltPathName));
    // Further error checks
    if (indexPath == hltConfig_.size()) {
      if (verbose_ > 1)
        edm::LogWarning("GenericTriggerEventFlag") << "HLT path \"" << hltPathName << "\" is not found in process "
                                                   << hltInputTag_.process() << " ==> decision: " << errorReplyHlt_;
      hltAlgoLogicParser.operandTokenVector().at(iPath).tokenResult = errorReplyHlt_;
      continue;
    }
    if (hltTriggerResults->error(indexPath)) {
      if (verbose_ > 1)
        edm::LogWarning("GenericTriggerEventFlag")
            << "HLT path \"" << hltPathName << "\" in error ==> decision: " << errorReplyHlt_;
      hltAlgoLogicParser.operandTokenVector().at(iPath).tokenResult = errorReplyHlt_;
      continue;
    }
    // Manipulate algo decision as stored in the parser
    const bool decision(hltTriggerResults->accept(indexPath));
    hltAlgoLogicParser.operandTokenVector().at(iPath).tokenResult = decision;
  }

  // Determine decision
  const bool hltDecision(hltAlgoLogicParser.expressionResult());
  return negExpr ? (!hltDecision) : hltDecision;
}

std::string GenericTriggerEventFlag::expandLogicalExpression(const std::vector<std::string>& targets,
                                                             const std::string& expr,
                                                             bool useAnd) const {
  // Find matching entries in the menu
  std::vector<std::string> matched;
  const std::string versionWildcard("_v*");
  if (expr.substr(expr.size() - versionWildcard.size()) == versionWildcard) {
    const std::string exprBase(expr.substr(0, expr.size() - versionWildcard.size()));
    matched = hltConfig_.restoreVersion(targets, exprBase);
  } else {
    matched = hltConfig_.matched(targets, expr);
  }

  // Return input, if no match is found
  if (matched.empty()) {
    if (verbose_ > 1)
      edm::LogWarning("GenericTriggerEventFlag") << "Logical expression: \"" << expr << "\" could not be resolved";
    return expr;
  }

  // Compose logical expression
  std::string expanded("(");
  for (unsigned iVers = 0; iVers < matched.size(); ++iVers) {
    if (iVers > 0)
      expanded.append(useAnd ? " AND " : " OR ");
    expanded.append(matched.at(iVers));
  }
  expanded.append(")");
  if (verbose_ > 1)
    edm::LogInfo("GenericTriggerEventFlag") << "Logical expression: \"" << expr << "\"\n"
                                            << "   --> expanded to  \"" << expanded << "\"";

  return expanded;
}

bool GenericTriggerEventFlag::negate(std::string& word) const {
  bool negate(false);
  if (word.at(0) == '~') {
    negate = true;
    word.erase(0, 1);
  }
  return negate;
}

std::vector<std::string> GenericTriggerEventFlag::expressionsFromDB(const std::string& key,
                                                                    const edm::EventSetup& setup) {
  if (key.empty())
    return std::vector<std::string>(1, emptyKeyError_);
  std::vector<edm::eventsetup::DataKey> labels;
  setup.get<AlCaRecoTriggerBitsRcd>().fillRegisteredDataKeys(labels);
  std::vector<edm::eventsetup::DataKey>::const_iterator iKey = labels.begin();
  while (iKey != labels.end() && iKey->name().value() != dbLabel_)
    ++iKey;
  if (iKey == labels.end()) {
    if (verbose_ > 0)
      edm::LogWarning("GenericTriggerEventFlag")
          << "Label " << dbLabel_ << " not found in DB for 'AlCaRecoTriggerBitsRcd'";
    return std::vector<std::string>(1, configError_);
  }
  auto const& alCaRecoTriggerBits = setup.getData(alCaRecoTriggerBitsToken_);
  const std::map<std::string, std::string>& expressionMap = alCaRecoTriggerBits.m_alcarecoToTrig;
  std::map<std::string, std::string>::const_iterator listIter = expressionMap.find(key);
  if (listIter == expressionMap.end()) {
    if (verbose_ > 0)
      edm::LogWarning("GenericTriggerEventFlag")
          << "No logical expressions found under key " << key << " in 'AlCaRecoTriggerBitsRcd'";
    return std::vector<std::string>(1, configError_);
  }
  return alCaRecoTriggerBits.decompose(listIter->second);
}

bool GenericTriggerEventFlag::allHLTPathsAreValid() const {
  if (not onHlt_) {
    return true;
  }

  if (not hltConfigInit_) {
    if (verbose_ > 0) {
      edm::LogWarning("GenericTriggerEventFlag::allHLTPathsAreValid()")
          << "HLTConfigProvider is not initialized, method will return \"false\"";
    }

    return false;
  }

  for (unsigned iExpr = 0; iExpr < hltLogicalExpressions_.size(); ++iExpr) {
    std::string hltLogicalExpression = hltLogicalExpressions_.at(iExpr);

    L1GtLogicParser hltAlgoLogicParser(hltLogicalExpression);

    if (hltAlgoLogicParser.operandTokenVector().empty()) {
      return false;
    }

    for (size_t iPath = 0; iPath < hltAlgoLogicParser.operandTokenVector().size(); ++iPath) {
      const std::string hltPathName(hltAlgoLogicParser.operandTokenVector().at(iPath).tokenName);

      const unsigned indexPath(hltConfig_.triggerIndex(hltPathName));

      if (indexPath == hltConfig_.size()) {
        if (verbose_ > 1) {
          edm::LogWarning("GenericTriggerEventFlag::allHLTPathsAreValid()")
              << "HLT path \"" << hltPathName << "\" is not found in process " << hltInputTag_.process();
        }

        return false;
      }
    }
  }

  return true;
}

void GenericTriggerEventFlag::fillPSetDescription(edm::ParameterSetDescription& desc) {
  desc.add<bool>("ReadPrescalesFromFile", false);
  desc.add<bool>("andOr", false);
  desc.add<bool>("andOrDcs", false);
  desc.add<bool>("andOrHlt", false);
  desc.add<bool>("andOrL1", false);
  desc.add<bool>("errorReplyDcs", false);
  desc.add<bool>("errorReplyHlt", false);
  desc.add<bool>("errorReplyL1", false);
  desc.add<bool>("l1BeforeMask", false);
  desc.add<bool>("stage2", false);
  desc.add<edm::InputTag>("dcsInputTag", edm::InputTag("scalersRawToDigi"));
  desc.add<edm::InputTag>("dcsRecordInputTag", edm::InputTag("onlineMetaDataDigis"));
  desc.add<edm::InputTag>("hltInputTag", edm::InputTag("TriggerResults::HLT"));
  desc.add<edm::InputTag>("l1tAlgBlkInputTag", edm::InputTag("gtStage2Digis"));
  desc.add<edm::InputTag>("l1tExtBlkInputTag", edm::InputTag("gtStage2Digis"));
  desc.add<std::string>("dbLabel", "");
  desc.add<std::string>("hltDBKey", "");
  desc.add<std::vector<int>>("dcsPartitions", {});
  desc.add<std::vector<std::string>>("hltPaths", {});
  desc.add<std::vector<std::string>>("l1Algorithms", {});
  desc.add<unsigned int>("verbosityLevel", 0);
}