L1TCaloParamsOnlineProd.cc File Reference

#include <iostream>
#include <fstream>
#include <stdexcept>
#include <sstream>
#include "CondTools/L1TriggerExt/interface/L1ConfigOnlineProdBaseExt.h"
#include "CondFormats/L1TObjects/interface/CaloParams.h"
#include "CondFormats/DataRecord/interface/L1TCaloStage2ParamsRcd.h"
#include "CondFormats/DataRecord/interface/L1TCaloParamsO2ORcd.h"
#include "L1Trigger/L1TCommon/interface/TrigSystem.h"
#include "L1Trigger/L1TCommon/interface/XmlConfigReader.h"
#include "L1Trigger/L1TCalorimeter/interface/CaloParamsHelper.h"
#include "xercesc/util/PlatformUtils.hpp"

Go to the source code of this file.

Classes
class	L1TCaloParamsOnlineProd

Functions
bool	readCaloLayer1OnlineSettings (l1t::CaloParamsHelper &paramsHelper, std::map< std::string, l1t::Setting > &conf, std::map< std::string, l1t::Mask > &)
bool	readCaloLayer2OnlineSettings (l1t::CaloParamsHelper &paramsHelper, std::map< std::string, l1t::Setting > &conf, std::map< std::string, l1t::Mask > &)
static const edm::ParameterSetDescriptionFillerPluginFactory::PMaker < edm::DescriptionFillerForESProducers < L1TCaloParamsOnlineProd > >	s_filler__LINE__ ("L1TCaloParamsOnlineProd")
static const edm::eventsetup::ModulePluginFactory::PMaker < edm::eventsetup::ModuleMaker < L1TCaloParamsOnlineProd > >	s_maker__LINE__ ("L1TCaloParamsOnlineProd")

Function Documentation

bool readCaloLayer1OnlineSettings	(	l1t::CaloParamsHelper &	paramsHelper,
		std::map< std::string, l1t::Setting > &	conf,
		std::map< std::string, l1t::Mask > &
	)

Definition at line 27 of file L1TCaloParamsOnlineProd.cc.

References ecal_dqm_sourceclient-live_cfg::cerr, l1t::CaloParamsHelper::setLayer1ECalScaleETBins(), l1t::CaloParamsHelper::setLayer1ECalScaleFactors(), l1t::CaloParamsHelper::setLayer1HCalScaleETBins(), l1t::CaloParamsHelper::setLayer1HCalScaleFactors(), l1t::CaloParamsHelper::setLayer1HFScaleETBins(), and l1t::CaloParamsHelper::setLayer1HFScaleFactors().

Referenced by L1TCaloParamsOnlineProd::newObject().

                                                                                                                                   {
  const char * expectedParams[] = {
    "layer1ECalScaleFactors",
    "layer1HCalScaleFactors",
    "layer1HFScaleFactors",
    "layer1ECalScaleETBins",
    "layer1HCalScaleETBins",
    "layer1HFScaleETBins"
  };
  for (const auto param : expectedParams) {
    if ( conf.find(param) == conf.end() ) {
      std::cerr << "Unable to locate expected CaloLayer1 parameter: " << param << " in L1 settings payload!";
      return false;
    }
  }
  // Layer 1 LUT specification
  paramsHelper.setLayer1ECalScaleFactors((conf["layer1ECalScaleFactors"].getVector<double>()));
  paramsHelper.setLayer1HCalScaleFactors((conf["layer1HCalScaleFactors"].getVector<double>()));
  paramsHelper.setLayer1HFScaleFactors  ((conf["layer1HFScaleFactors"]  .getVector<double>()));
  paramsHelper.setLayer1ECalScaleETBins(conf["layer1ECalScaleETBins"].getVector<int>());
  paramsHelper.setLayer1HCalScaleETBins(conf["layer1HCalScaleETBins"].getVector<int>());
  paramsHelper.setLayer1HFScaleETBins  (conf["layer1HFScaleETBins"]  .getVector<int>());

  return true;
}

bool readCaloLayer2OnlineSettings	(	l1t::CaloParamsHelper &	paramsHelper,
		std::map< std::string, l1t::Setting > &	conf,
		std::map< std::string, l1t::Mask > &
	)

Definition at line 54 of file L1TCaloParamsOnlineProd.cc.

References ecal_dqm_sourceclient-live_cfg::cerr, i, l1t::CaloParamsHelper::setEgCalibrationLUT(), l1t::CaloParamsHelper::setEgEtaCut(), l1t::CaloParamsHelper::setEgIsolationLUT(), l1t::CaloParamsHelper::setEgMaxPtHOverE(), l1t::CaloParamsHelper::setEgNeighbourThreshold(), l1t::CaloParamsHelper::setEgSeedThreshold(), l1t::CaloParamsHelper::setEtSumEcalSumPUSLUT(), l1t::CaloParamsHelper::setEtSumEtaMax(), l1t::CaloParamsHelper::setEtSumEtThreshold(), l1t::CaloParamsHelper::setEtSumEttPUSLUT(), l1t::CaloParamsHelper::setEtSumXPUSLUT(), l1t::CaloParamsHelper::setIsoTauEtaMax(), l1t::CaloParamsHelper::setJetCalibrationLUT(), l1t::CaloParamsHelper::setJetSeedThreshold(), l1t::CaloParamsHelper::setTauCalibrationLUT(), l1t::CaloParamsHelper::setTauIsolationLUT(), l1t::CaloParamsHelper::setTauIsolationLUT2(), l1t::CaloParamsHelper::setTauNeighbourThreshold(), and l1t::CaloParamsHelper::setTauSeedThreshold().

Referenced by L1TCaloParamsOnlineProd::newObject().

                                                                                                                                   {
  const char * expectedParams[] = {
    "leptonSeedThreshold",
    "leptonTowerThreshold",
    "pileUpTowerThreshold",
    "jetSeedThreshold",
    "jetMaxEta",
    "HTMHT_maxJetEta",
    "HT_jetThreshold",
    "MHT_jetThreshold",
    "jetEnergyCalibLUT",
    "ETMET_maxTowerEta",
    "ET_energyCalibLUT",
    "ecalET_energyCalibLUT",
    "METX_energyCalibLUT",
    "METY_energyCalibLUT",
    "egammaRelaxationThreshold",
    "egammaMaxEta",
    "egammaEnergyCalibLUT",
    "egammaIsoLUT",
    "tauMaxEta",
    "tauEnergyCalibLUT",
    "tauIsoLUT1",
    "tauIsoLUT2"
  };
  for (const auto param : expectedParams) {
    if ( conf.find(param) == conf.end() ) {
      std::cerr << "Unable to locate expected CaloLayer2 parameter: " << param << " in L1 settings payload!";
      return false;
    }
  }
  // Layer 2 params specification
  paramsHelper.setEgSeedThreshold((conf["leptonSeedThreshold"].getValue<int>())/2);
  paramsHelper.setTauSeedThreshold((conf["leptonSeedThreshold"].getValue<int>())/2);
  paramsHelper.setEgNeighbourThreshold((conf["leptonTowerThreshold"].getValue<int>())/2);
  paramsHelper.setTauNeighbourThreshold((conf["leptonTowerThreshold"].getValue<int>())/2);
  paramsHelper.setJetSeedThreshold((conf["jetSeedThreshold"].getValue<int>())/2);

  // Currently not used // paramsHelper.setEgPileupTowerThresh((conf["pileUpTowerThreshold"].getValue<int>())); 
  // Currently not used // paramsHelper.setTauPileupTowerThresh((conf["pileUpTowerThreshold"].getValue<int>())); 
  // Currently not used // paramsHelper.setJetMaxEta((conf["jetMaxEta"].getValue<int>()));
  
  std::vector<int> etSumEtaMax;
  std::vector<int> etSumEtThresh;
  
  etSumEtaMax.push_back(conf["ETMET_maxTowerEta"].getValue<int>());
  etSumEtaMax.push_back(conf["HTMHT_maxJetEta"].getValue<int>());
  etSumEtaMax.push_back(conf["ETMET_maxTowerEta"].getValue<int>());
  etSumEtaMax.push_back(conf["HTMHT_maxJetEta"].getValue<int>());
  
  etSumEtThresh.push_back(0); // ETT tower threshold
  etSumEtThresh.push_back(conf["HT_jetThreshold"].getValue<int>());
  etSumEtThresh.push_back(0); // ETM tower threshold
  etSumEtThresh.push_back(conf["MHT_jetThreshold"].getValue<int>());

  for (uint i=0; i<4; ++i) {
    paramsHelper.setEtSumEtaMax(i, etSumEtaMax.at(i));
    paramsHelper.setEtSumEtThreshold(i, etSumEtThresh.at(i));
  }

  std::stringstream oss;

  std::vector<uint32_t> jetEnergyCalibLUT = conf["jetEnergyCalibLUT"].getVector<uint32_t>();
  oss <<"#<header> V1 "<< ( 32 - __builtin_clz( uint32_t(jetEnergyCalibLUT.size()-1) ) ) <<" 63 </header> "<<std::endl; // hardcode max bits for data
  for(unsigned int i=0; i<jetEnergyCalibLUT.size(); i++) oss << i << " " << jetEnergyCalibLUT[i] << std::endl;

  std::istringstream iss1( oss.str() );
  paramsHelper.setJetCalibrationLUT( l1t::LUT( (std::istream&)iss1 ) );
  oss.str("");

  std::vector<int> etSumEttPUSLUT = conf["ET_energyCalibLUT"].getVector<int>();
  oss <<"#<header> V1 "<< ( 32 - __builtin_clz( uint32_t(etSumEttPUSLUT.size()-1) ) ) <<" 63 </header> "<<std::endl; // hardcode max bits for data
  for(unsigned int i=0; i<etSumEttPUSLUT.size(); i++) oss << i << " " << etSumEttPUSLUT[i] << std::endl;

  std::istringstream iss2( oss.str() );
  paramsHelper.setEtSumEttPUSLUT( l1t::LUT( (std::istream&)iss2 ) );
  oss.str("");

  std::vector<int> etSumEcalSumPUTLUT = conf["ecalET_energyCalibLUT"].getVector<int>();
  oss <<"#<header> V1 "<< ( 32 - __builtin_clz( uint32_t(etSumEcalSumPUTLUT.size()-1) ) ) <<" 63 </header> "<<std::endl; // hardcode max bits for data
  for(unsigned int i=0; i<etSumEcalSumPUTLUT.size(); i++) oss << i << " " << etSumEcalSumPUTLUT[i] << std::endl;

  std::istringstream iss3( oss.str() );
  paramsHelper.setEtSumEcalSumPUSLUT( l1t::LUT( (std::istream&)iss3 ) );
  oss.str("");

  std::vector<int> etSumXPUSLUT = conf["METX_energyCalibLUT"].getVector<int>();
  oss <<"#<header> V1 "<< ( 32 - __builtin_clz( uint32_t(etSumXPUSLUT.size()-1) ) ) <<" 63 </header> "<<std::endl; // hardcode max bits for data
  for(unsigned int i=0; i<etSumXPUSLUT.size(); i++) oss << i << " " << etSumXPUSLUT[i] << std::endl;

  std::istringstream iss4( oss.str() );
  paramsHelper.setEtSumXPUSLUT( l1t::LUT( (std::istream&)iss4 ) );
  oss.str("");

  paramsHelper.setEgMaxPtHOverE((conf["egammaRelaxationThreshold"].getValue<int>()));
  paramsHelper.setEgEtaCut((conf["egammaMaxEta"].getValue<int>()));


  std::vector<int> egCalibrationLUT = conf["egammaEnergyCalibLUT"].getVector<int>();
  oss <<"#<header> V1 "<< ( 32 - __builtin_clz( uint32_t(egCalibrationLUT.size()-1) ) ) <<" 63 </header> "<<std::endl; // hardcode max bits for data
  for(unsigned int i=0; i<egCalibrationLUT.size(); i++) oss << i << " " << egCalibrationLUT[i] << std::endl;

  std::istringstream iss5( oss.str() );
  paramsHelper.setEgCalibrationLUT( l1t::LUT( (std::istream&)iss5 ) );
  oss.str("");

  std::vector<int> egIsolationLUT = conf["egammaIsoLUT"].getVector<int>();
  oss <<"#<header> V1 "<< ( 32 - __builtin_clz( uint32_t(egIsolationLUT.size()-1) ) ) <<" 63 </header> "<<std::endl; // hardcode max bits for data
  for(unsigned int i=0; i<egIsolationLUT.size(); i++) oss << i << " " << egIsolationLUT[i] << std::endl;

  std::istringstream iss6( oss.str() );
  paramsHelper.setEgIsolationLUT( l1t::LUT( (std::istream&)iss6 ) );
  oss.str("");

  //std::cout<<"egammaIsoLUT: "<<std::endl;
  //for(unsigned int i=0; i<paramsHelper.egIsolationLUT()->maxSize(); i++) std::cout << std::setprecision(14) << paramsHelper.egIsolationLUT()->data(i) <<", ";
  //std::cout << std::endl;

  paramsHelper.setIsoTauEtaMax((conf["tauMaxEta"].getValue<int>()));

  std::vector<int> tauCalibrationLUT = conf["tauEnergyCalibLUT"].getVector<int>();
  oss <<"#<header> V1 "<< ( 32 - __builtin_clz( uint32_t(tauCalibrationLUT.size()-1) ) ) <<" 63 </header> "<<std::endl; // hardcode max bits for data
  for(unsigned int i=0; i<tauCalibrationLUT.size(); i++) oss << i << " " << tauCalibrationLUT[i] << std::endl;

  std::istringstream iss7( oss.str() );
  paramsHelper.setTauCalibrationLUT(l1t::LUT( (std::istream&)iss7 ));
  oss.str("");

  std::vector<int> tauIsolationLUT = conf["tauIsoLUT1"].getVector<int>();
  oss <<"#<header> V1 "<< ( 32 - __builtin_clz( uint32_t(tauIsolationLUT.size()-1) ) ) <<" 63 </header> "<<std::endl; // hardcode max bits for data
  for(unsigned int i=0; i<tauIsolationLUT.size(); i++) oss << i << " " << tauIsolationLUT[i] << std::endl;

  std::istringstream iss8( oss.str() );
  paramsHelper.setTauIsolationLUT( l1t::LUT((std::istream&)iss8 ) );
  oss.str("");

  std::vector<int> tauIsolationLUT2 = conf["tauIsoLUT2"].getVector<int>();
  oss <<"#<header> V1 "<< ( 32 - __builtin_clz( uint32_t(tauIsolationLUT2.size()-1) ) ) <<" 63 </header> "<<std::endl; // hardcode max bits for data
  for(unsigned int i=0; i<tauIsolationLUT2.size(); i++) oss << i << " " << tauIsolationLUT2[i] << std::endl;

  std::istringstream iss9( oss.str() );
  paramsHelper.setTauIsolationLUT2( l1t::LUT( (std::istream&)iss9 ) );
  oss.str("");

  return true;
}

static const edm::ParameterSetDescriptionFillerPluginFactory::PMaker<edm::DescriptionFillerForESProducers< L1TCaloParamsOnlineProd > > s_filler__LINE__ ( "L1TCaloParamsOnlineProd"  )

static

static const edm::eventsetup::ModulePluginFactory ::PMaker< edm::eventsetup::ModuleMaker< L1TCaloParamsOnlineProd > > s_maker__LINE__ ( "L1TCaloParamsOnlineProd"  )

static