File indexing completed on 2024-04-06 12:20:36
0001
0002
0003
0004
0005
0006 #include "L1Trigger/L1TGlobal/interface/L1TGlobalUtil.h"
0007
0008 #include <fstream>
0009 #include <iostream>
0010 #include <memory>
0011
0012 #include "FWCore/Framework/interface/Event.h"
0013 #include "FWCore/ParameterSet/interface/ParameterSet.h"
0014 #include "FWCore/Utilities/interface/EDGetToken.h"
0015 #include "FWCore/Utilities/interface/InputTag.h"
0016 #include "FWCore/Utilities/interface/Transition.h"
0017
0018 #include "FWCore/MessageLogger/interface/MessageLogger.h"
0019 #include "FWCore/MessageLogger/interface/MessageDrop.h"
0020
0021
0022 l1t::L1TGlobalUtil::L1TGlobalUtil() : m_l1GtMenu(nullptr) {
0023
0024 m_l1GtMenuCacheID = 0ULL;
0025 m_l1GtPfAlgoCacheID = 0ULL;
0026 m_filledPrescales = false;
0027
0028 edm::FileInPath f1("L1Trigger/L1TGlobal/data/Luminosity/startup/prescale_L1TGlobal.csv");
0029 m_preScaleFileName = f1.fullPath();
0030 m_numberPhysTriggers = 512;
0031 m_PreScaleColumn = 0;
0032 m_readPrescalesFromFile = false;
0033
0034 m_prescaleFactorsAlgoTrig = nullptr;
0035 m_triggerMaskAlgoTrig = nullptr;
0036 }
0037
0038 l1t::L1TGlobalUtil::L1TGlobalUtil(edm::ParameterSet const& pset,
0039 edm::ConsumesCollector&& iC,
0040 UseEventSetupIn useEventSetupIn)
0041 : L1TGlobalUtil(pset, iC, useEventSetupIn) {}
0042
0043 l1t::L1TGlobalUtil::L1TGlobalUtil(edm::ParameterSet const& pset,
0044 edm::ConsumesCollector& iC,
0045 UseEventSetupIn useEventSetupIn)
0046 : L1TGlobalUtil() {
0047 m_l1tGlobalUtilHelper = std::make_unique<L1TGlobalUtilHelper>(pset, iC);
0048 m_readPrescalesFromFile = m_l1tGlobalUtilHelper->readPrescalesFromFile();
0049 eventSetupConsumes(iC, useEventSetupIn);
0050 }
0051
0052
0053 l1t::L1TGlobalUtil::~L1TGlobalUtil() {
0054
0055 }
0056
0057
0058 bool l1t::L1TGlobalUtil::valid() const { return m_l1GtMenuCacheID != 0ULL and m_l1GtMenu != nullptr; }
0059
0060 void l1t::L1TGlobalUtil::OverridePrescalesAndMasks(std::string filename, unsigned int psColumn) {
0061 edm::FileInPath f1("L1Trigger/L1TGlobal/data/Luminosity/startup/" + filename);
0062 m_preScaleFileName = f1.fullPath();
0063 m_PreScaleColumn = psColumn;
0064 }
0065
0066 void l1t::L1TGlobalUtil::retrieveL1(const edm::Event& iEvent, const edm::EventSetup& evSetup) {
0067
0068 bool isRun = false;
0069 retrieveL1Setup(evSetup, isRun);
0070
0071 retrieveL1Event(iEvent, evSetup);
0072 }
0073
0074 void l1t::L1TGlobalUtil::retrieveL1(const edm::Event& iEvent,
0075 const edm::EventSetup& evSetup,
0076 edm::EDGetToken gtAlgToken) {
0077
0078 bool isRun = false;
0079 retrieveL1Setup(evSetup, isRun);
0080
0081 retrieveL1Event(iEvent, evSetup, gtAlgToken);
0082 }
0083
0084 void l1t::L1TGlobalUtil::retrieveL1Setup(const edm::EventSetup& evSetup) {
0085 bool isRun = true;
0086 retrieveL1Setup(evSetup, isRun);
0087 }
0088
0089 void l1t::L1TGlobalUtil::retrieveL1Setup(const edm::EventSetup& evSetup, bool isRun) {
0090
0091
0092 auto menuRcd = evSetup.get<L1TUtmTriggerMenuRcd>();
0093 unsigned long long l1GtMenuCacheID = menuRcd.cacheIdentifier();
0094
0095 if (m_l1GtMenuCacheID != l1GtMenuCacheID) {
0096 if (isRun) {
0097 m_l1GtMenu = &menuRcd.get(m_L1TUtmTriggerMenuRunToken);
0098 } else {
0099 m_l1GtMenu = &menuRcd.get(m_L1TUtmTriggerMenuEventToken);
0100 }
0101
0102
0103 m_algorithmMap = &(m_l1GtMenu->getAlgorithmMap());
0104
0105
0106 resetDecisionVectors();
0107
0108 m_l1GtMenuCacheID = l1GtMenuCacheID;
0109 }
0110
0111 if (!m_readPrescalesFromFile) {
0112 auto vetosRcd = evSetup.get<L1TGlobalPrescalesVetosFractRcd>();
0113 unsigned long long l1GtPfAlgoCacheID = vetosRcd.cacheIdentifier();
0114
0115 if (m_l1GtPfAlgoCacheID != l1GtPfAlgoCacheID) {
0116
0117
0118
0119 resetPrescaleVectors();
0120 resetMaskVectors();
0121 m_PreScaleColumn = 0;
0122 m_numberOfPreScaleColumns = 0;
0123 m_numberPhysTriggers = 0;
0124
0125 const L1TGlobalPrescalesVetosFract* es = nullptr;
0126 if (isRun) {
0127 es = &vetosRcd.get(m_L1TGlobalPrescalesVetosFractRunToken);
0128 } else {
0129 es = &vetosRcd.get(m_L1TGlobalPrescalesVetosFractEventToken);
0130 }
0131 m_l1GtPrescalesVetoes = PrescalesVetosFractHelper::readFromEventSetup(es);
0132
0133 m_prescaleFactorsAlgoTrig = &(m_l1GtPrescalesVetoes->prescaleTable());
0134 m_numberOfPreScaleColumns = m_prescaleFactorsAlgoTrig->size();
0135 m_numberPhysTriggers =
0136 (*m_prescaleFactorsAlgoTrig)[0].size();
0137
0138 m_triggerMaskAlgoTrig = &(m_l1GtPrescalesVetoes->triggerAlgoBxMask());
0139
0140 m_l1GtPfAlgoCacheID = l1GtPfAlgoCacheID;
0141 }
0142 } else {
0143
0144
0145 if (!m_filledPrescales) {
0146 resetPrescaleVectors();
0147 resetMaskVectors();
0148
0149 loadPrescalesAndMasks();
0150
0151
0152 m_prescaleFactorsAlgoTrig = &m_initialPrescaleFactorsAlgoTrig;
0153 m_triggerMaskAlgoTrig = &m_initialTriggerMaskAlgoTrig;
0154 m_filledPrescales = true;
0155 }
0156 }
0157
0158
0159 if (m_PreScaleColumn > m_prescaleFactorsAlgoTrig->size()) {
0160 LogTrace("l1t|Global") << "\nNo Prescale Set: " << m_PreScaleColumn
0161 << "\nMax Prescale Set value : " << m_prescaleFactorsAlgoTrig->size()
0162 << "\nSetting prescale column to 0" << std::endl;
0163 m_PreScaleColumn = 0;
0164 }
0165
0166 const std::vector<double>& prescaleSet = (*m_prescaleFactorsAlgoTrig)[m_PreScaleColumn];
0167
0168 for (std::map<std::string, L1TUtmAlgorithm>::const_iterator itAlgo = m_algorithmMap->begin();
0169 itAlgo != m_algorithmMap->end();
0170 itAlgo++) {
0171
0172 std::string algName = itAlgo->first;
0173 int algBit = (itAlgo->second).getIndex();
0174
0175 (m_prescales[algBit]).first = algName;
0176 if (size_t(algBit) < prescaleSet.size()) {
0177 (m_prescales[algBit]).second = prescaleSet[algBit];
0178 }
0179 LogDebug("l1t|Global") << "Number of bunch crossings stored: " << (*m_triggerMaskAlgoTrig).size() << endl;
0180
0181 const std::map<int, std::vector<int> >* triggerAlgoMaskAlgoTrig = m_triggerMaskAlgoTrig;
0182 std::map<int, std::vector<int> >::const_iterator it = triggerAlgoMaskAlgoTrig->begin();
0183
0184 std::vector<int> maskedBxs;
0185 (m_masks[algBit]).first = algName;
0186 (m_masks[algBit]).second = maskedBxs;
0187 while (it != triggerAlgoMaskAlgoTrig->end()) {
0188 std::vector<int> masks = it->second;
0189
0190
0191 for (unsigned int imask = 0; imask < masks.size(); imask++) {
0192 if (masks.at(imask) == algBit)
0193 maskedBxs.push_back(it->first);
0194
0195 }
0196
0197 it++;
0198 }
0199
0200 if (!maskedBxs.empty()) {
0201 LogDebug("l1t|Global") << "i Algo: " << algBit << "\t" << algName << " masked\n";
0202 for (unsigned int ibx = 0; ibx < maskedBxs.size(); ibx++) {
0203
0204 (m_masks[algBit]).second = maskedBxs;
0205 }
0206 }
0207 }
0208 }
0209
0210 void l1t::L1TGlobalUtil::retrieveL1Event(const edm::Event& iEvent, const edm::EventSetup& evSetup) {
0211 retrieveL1Event(iEvent, evSetup, m_l1tGlobalUtilHelper->l1tAlgBlkToken());
0212 }
0213
0214 void l1t::L1TGlobalUtil::retrieveL1Event(const edm::Event& iEvent,
0215 const edm::EventSetup& evSetup,
0216 edm::EDGetToken gtAlgToken) {
0217
0218 iEvent.getByToken(gtAlgToken, m_uGtAlgBlk);
0219 m_finalOR = false;
0220
0221
0222 if (m_uGtAlgBlk.isValid()) {
0223
0224 std::vector<GlobalAlgBlk>::const_iterator algBlk = m_uGtAlgBlk->begin(0);
0225 if (algBlk != m_uGtAlgBlk->end(0)) {
0226 if (!m_readPrescalesFromFile) {
0227 m_PreScaleColumn = static_cast<unsigned int>(algBlk->getPreScColumn());
0228
0229
0230 if (iEvent.run() == 1) {
0231 if (m_prescaleFactorsAlgoTrig->size() == 1 && m_PreScaleColumn == 1)
0232 m_PreScaleColumn = 0;
0233 }
0234
0235
0236 if (m_PreScaleColumn >= m_prescaleFactorsAlgoTrig->size()) {
0237 LogDebug("l1t|Global") << "Prescale column extracted from GlobalAlgBlk too large: " << m_PreScaleColumn
0238 << "\tMaximum value allowed: " << m_prescaleFactorsAlgoTrig->size() - 1
0239 << "\tResetting prescale column to 0" << std::endl;
0240 m_PreScaleColumn = 0;
0241 }
0242 }
0243 const std::vector<double>& prescaleSet = (*m_prescaleFactorsAlgoTrig)[m_PreScaleColumn];
0244
0245
0246 m_finalOR = algBlk->getFinalOR();
0247
0248
0249
0250 for (std::map<std::string, L1TUtmAlgorithm>::const_iterator itAlgo = m_algorithmMap->begin();
0251 itAlgo != m_algorithmMap->end();
0252 itAlgo++) {
0253
0254 std::string algName = itAlgo->first;
0255 int algBit = (itAlgo->second).getIndex();
0256
0257 bool decisionInitial = algBlk->getAlgoDecisionInitial(algBit);
0258 (m_decisionsInitial[algBit]).first = algName;
0259 (m_decisionsInitial[algBit]).second = decisionInitial;
0260
0261 bool decisionInterm = algBlk->getAlgoDecisionInterm(algBit);
0262 (m_decisionsInterm[algBit]).first = algName;
0263 (m_decisionsInterm[algBit]).second = decisionInterm;
0264
0265 bool decisionFinal = algBlk->getAlgoDecisionFinal(algBit);
0266 (m_decisionsFinal[algBit]).first = algName;
0267 (m_decisionsFinal[algBit]).second = decisionFinal;
0268
0269 (m_prescales[algBit]).first = algName;
0270 if (size_t(algBit) < prescaleSet.size()) {
0271 (m_prescales[algBit]).second = prescaleSet[algBit];
0272 }
0273 }
0274 } else {
0275
0276 }
0277 } else {
0278
0279 }
0280 }
0281
0282 void l1t::L1TGlobalUtil::loadPrescalesAndMasks() {
0283 std::ifstream inputPrescaleFile;
0284
0285 inputPrescaleFile.open(m_preScaleFileName);
0286
0287 std::vector<std::vector<int> > vec;
0288 std::vector<std::vector<double> > prescale_vec;
0289
0290 if (inputPrescaleFile) {
0291 std::string prefix1("#");
0292 std::string prefix2("-1");
0293
0294 std::string line;
0295
0296 bool first = true;
0297
0298 while (getline(inputPrescaleFile, line)) {
0299 if (!line.compare(0, prefix1.size(), prefix1))
0300 continue;
0301
0302
0303 istringstream split(line);
0304 int value;
0305 int col = 0;
0306 char sep;
0307
0308 while (split >> value) {
0309 if (first) {
0310
0311 vec.push_back(std::vector<int>());
0312 }
0313
0314 vec[col].push_back(value);
0315 ++col;
0316
0317
0318 split >> sep;
0319 }
0320
0321
0322
0323 first = false;
0324 }
0325
0326 int NumPrescaleSets = 0;
0327 for (int iCol = 0; iCol < int(vec.size()); iCol++) {
0328 if (!vec[iCol].empty()) {
0329 int firstRow = vec[iCol][0];
0330
0331 if (firstRow >= 0)
0332 NumPrescaleSets++;
0333
0334
0335 }
0336 }
0337
0338
0339 if (NumPrescaleSets > 0) {
0340
0341 for (int iSet = 0; iSet < NumPrescaleSets; iSet++) {
0342 prescale_vec.push_back(std::vector<double>());
0343 for (unsigned int iBit = 0; iBit < m_numberPhysTriggers; ++iBit) {
0344 int inputDefaultPrescale = 1;
0345 prescale_vec[iSet].push_back(inputDefaultPrescale);
0346 }
0347 }
0348
0349
0350 for (int iBit = 1; iBit < int(vec[0].size()); iBit++) {
0351 unsigned int algoBit = vec[0][iBit];
0352
0353 if (algoBit < m_numberPhysTriggers) {
0354 for (int iSet = 0; iSet < int(vec.size()); iSet++) {
0355 int useSet = -1;
0356 if (!vec[iSet].empty()) {
0357 useSet = vec[iSet][0];
0358 }
0359 useSet -= 1;
0360
0361 if (useSet < 0)
0362 continue;
0363
0364 int prescale = vec[iSet][iBit];
0365 prescale_vec[useSet][algoBit] = prescale;
0366 }
0367 } else {
0368 LogTrace("l1t|Global") << "\nPrescale file has algo bit: " << algoBit
0369 << "\nThis is larger than the number of triggers: " << m_numberPhysTriggers
0370 << "\nSomething is wrong. Ignoring." << std::endl;
0371 }
0372 }
0373 }
0374
0375 } else {
0376 LogTrace("l1t|Global") << "\nCould not find file: " << m_preScaleFileName
0377 << "\nFilling the prescale vectors with prescale 1"
0378 << "\nSetting prescale set to 0" << std::endl;
0379
0380 m_PreScaleColumn = 0;
0381
0382 for (int col = 0; col < 1; col++) {
0383 prescale_vec.push_back(std::vector<double>());
0384 for (unsigned int iBit = 0; iBit < m_numberPhysTriggers; ++iBit) {
0385 int inputDefaultPrescale = 0;
0386 prescale_vec[col].push_back(inputDefaultPrescale);
0387 }
0388 }
0389 }
0390
0391 inputPrescaleFile.close();
0392
0393 m_initialPrescaleFactorsAlgoTrig = prescale_vec;
0394
0395 std::map<int, std::vector<int> > m_initialTriggerMaskAlgoTrig;
0396 }
0397
0398 void l1t::L1TGlobalUtil::eventSetupConsumes(edm::ConsumesCollector& iC, UseEventSetupIn useEventSetupIn) {
0399 if (useEventSetupIn == UseEventSetupIn::Run || useEventSetupIn == UseEventSetupIn::RunAndEvent) {
0400 m_L1TUtmTriggerMenuRunToken = iC.esConsumes<L1TUtmTriggerMenu, L1TUtmTriggerMenuRcd, edm::Transition::BeginRun>();
0401 if (!m_readPrescalesFromFile) {
0402 m_L1TGlobalPrescalesVetosFractRunToken =
0403 iC.esConsumes<L1TGlobalPrescalesVetosFract, L1TGlobalPrescalesVetosFractRcd, edm::Transition::BeginRun>();
0404 }
0405 }
0406 if (useEventSetupIn == UseEventSetupIn::Event || useEventSetupIn == UseEventSetupIn::RunAndEvent) {
0407 m_L1TUtmTriggerMenuEventToken = iC.esConsumes<L1TUtmTriggerMenu, L1TUtmTriggerMenuRcd>();
0408 if (!m_readPrescalesFromFile) {
0409 m_L1TGlobalPrescalesVetosFractEventToken =
0410 iC.esConsumes<L1TGlobalPrescalesVetosFract, L1TGlobalPrescalesVetosFractRcd>();
0411 }
0412 }
0413 }
0414
0415 void l1t::L1TGlobalUtil::resetDecisionVectors() {
0416
0417 m_decisionsInitial.clear();
0418 m_decisionsInitial.resize(m_maxNumberPhysTriggers);
0419 m_decisionsInterm.clear();
0420 m_decisionsInterm.resize(m_maxNumberPhysTriggers);
0421 m_decisionsFinal.clear();
0422 m_decisionsFinal.resize(m_maxNumberPhysTriggers);
0423
0424 for (unsigned int algBit = 0; algBit < m_maxNumberPhysTriggers; algBit++) {
0425 (m_decisionsInitial[algBit]).first = "NULL";
0426 (m_decisionsInitial[algBit]).second = false;
0427
0428 (m_decisionsInterm[algBit]).first = "NULL";
0429 (m_decisionsInterm[algBit]).second = false;
0430
0431 (m_decisionsFinal[algBit]).first = "NULL";
0432 (m_decisionsFinal[algBit]).second = false;
0433 }
0434 }
0435
0436 void l1t::L1TGlobalUtil::resetPrescaleVectors() {
0437
0438 m_prescales.clear();
0439 m_prescales.resize(m_maxNumberPhysTriggers);
0440
0441 for (unsigned int algBit = 0; algBit < m_maxNumberPhysTriggers; algBit++) {
0442 (m_prescales[algBit]).first = "NULL";
0443 (m_prescales[algBit]).second = 1;
0444 }
0445 }
0446
0447 void l1t::L1TGlobalUtil::resetMaskVectors() {
0448
0449 m_masks.clear();
0450 m_masks.resize(m_maxNumberPhysTriggers);
0451
0452 for (unsigned int algBit = 0; algBit < m_maxNumberPhysTriggers; algBit++) {
0453 (m_masks[algBit]).first = "NULL";
0454
0455 }
0456 }
0457
0458 const bool l1t::L1TGlobalUtil::getAlgBitFromName(const std::string& algName, int& bit) const {
0459 std::map<std::string, L1TUtmAlgorithm>::const_iterator itAlgo = m_algorithmMap->find(algName);
0460 if (itAlgo != m_algorithmMap->end()) {
0461 bit = (itAlgo->second).getIndex();
0462 return true;
0463 }
0464
0465 return false;
0466 }
0467
0468 const bool l1t::L1TGlobalUtil::getAlgNameFromBit(int& bit, std::string& algName) const {
0469
0470 if ((m_decisionsInitial[bit]).first != "NULL") {
0471 algName = (m_decisionsInitial[bit]).first;
0472 return true;
0473 }
0474 return false;
0475 }
0476
0477 const bool l1t::L1TGlobalUtil::getInitialDecisionByBit(int& bit, bool& decision) const {
0478
0479
0480
0481
0482
0483
0484 if ((m_decisionsInitial[bit]).first != "NULL") {
0485 decision = (m_decisionsInitial[bit]).second;
0486 return true;
0487 }
0488
0489 return false;
0490 }
0491 const bool l1t::L1TGlobalUtil::getIntermDecisionByBit(int& bit, bool& decision) const {
0492
0493 if ((m_decisionsInterm[bit]).first != "NULL") {
0494 decision = (m_decisionsInterm[bit]).second;
0495 return true;
0496 }
0497
0498 return false;
0499 }
0500 const bool l1t::L1TGlobalUtil::getFinalDecisionByBit(int& bit, bool& decision) const {
0501
0502 if ((m_decisionsFinal[bit]).first != "NULL") {
0503 decision = (m_decisionsFinal[bit]).second;
0504 return true;
0505 }
0506
0507 return false;
0508 }
0509 const bool l1t::L1TGlobalUtil::getPrescaleByBit(int& bit, double& prescale) const {
0510
0511 if ((m_prescales[bit]).first != "NULL") {
0512 prescale = (m_prescales[bit]).second;
0513 return true;
0514 }
0515
0516 return false;
0517 }
0518 const bool l1t::L1TGlobalUtil::getMaskByBit(int& bit, std::vector<int>& mask) const {
0519
0520 if ((m_masks[bit]).first != "NULL") {
0521 mask = (m_masks[bit]).second;
0522 return true;
0523 }
0524
0525 return false;
0526 }
0527
0528 const bool l1t::L1TGlobalUtil::getInitialDecisionByName(const std::string& algName, bool& decision) const {
0529 int bit = -1;
0530 if (getAlgBitFromName(algName, bit)) {
0531 decision = (m_decisionsInitial[bit]).second;
0532 return true;
0533 }
0534
0535 return false;
0536 }
0537
0538 const bool l1t::L1TGlobalUtil::getIntermDecisionByName(const std::string& algName, bool& decision) const {
0539 int bit = -1;
0540 if (getAlgBitFromName(algName, bit)) {
0541 decision = (m_decisionsInterm[bit]).second;
0542 return true;
0543 }
0544
0545 return false;
0546 }
0547
0548 const bool l1t::L1TGlobalUtil::getFinalDecisionByName(const std::string& algName, bool& decision) const {
0549 int bit = -1;
0550 if (getAlgBitFromName(algName, bit)) {
0551 decision = (m_decisionsFinal[bit]).second;
0552 return true;
0553 }
0554
0555 return false;
0556 }
0557 const bool l1t::L1TGlobalUtil::getPrescaleByName(const std::string& algName, double& prescale) const {
0558 int bit = -1;
0559 if (getAlgBitFromName(algName, bit)) {
0560 prescale = (m_prescales[bit]).second;
0561 return true;
0562 }
0563
0564 return false;
0565 }
0566 const bool l1t::L1TGlobalUtil::getMaskByName(const std::string& algName, std::vector<int>& mask) const {
0567 int bit = -1;
0568 if (getAlgBitFromName(algName, bit)) {
0569 mask = (m_masks[bit]).second;
0570 return true;
0571 }
0572
0573 return false;
0574 }