L1RCTParametersOnlineProd Class Reference

Inheritance diagram for L1RCTParametersOnlineProd:

Public Member Functions
void	fillScaleFactors (const l1t::OMDSReader::QueryResults &results, std::vector< double > &output, int nfactor=1)
	L1RCTParametersOnlineProd (const edm::ParameterSet &)
std::unique_ptr< L1RCTParameters >	newObject (const std::string &objectKey) override
	~L1RCTParametersOnlineProd () override
Public Member Functions inherited from L1ConfigOnlineProdBase< L1RCTParametersRcd, L1RCTParameters >
	L1ConfigOnlineProdBase (const edm::ParameterSet &)
virtual std::unique_ptr< L1RCTParameters >	produce (const L1RCTParametersRcd &iRecord)
	~L1ConfigOnlineProdBase () override
Public Member Functions inherited from edm::ESProducer
	ESProducer ()
	ESProducer (const ESProducer &)=delete
	ESProducer (ESProducer &&)=delete
ESProxyIndex const *	getTokenIndices (unsigned int iIndex) const
ESRecordIndex const *	getTokenRecordIndices (unsigned int iIndex) const
bool	hasMayConsumes () const noexcept
size_t	numberOfTokenIndices (unsigned int iIndex) const
ESProducer &	operator= (const ESProducer &)=delete
ESProducer &	operator= (ESProducer &&)=delete
SerialTaskQueueChain &	queue ()
template<typename Record >
std::optional< std::vector< ESProxyIndex > >	updateFromMayConsumes (unsigned int iIndex, const Record &iRecord) const
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &) final
	~ESProducer () noexcept(false) override
Public Member Functions inherited from edm::ESProxyFactoryProducer
	ESProxyFactoryProducer ()
	ESProxyFactoryProducer (const ESProxyFactoryProducer &)=delete
const ESProxyFactoryProducer &	operator= (const ESProxyFactoryProducer &)=delete
	~ESProxyFactoryProducer () noexcept(false) override
Public Member Functions inherited from edm::eventsetup::DataProxyProvider
void	createKeyedProxies (EventSetupRecordKey const &key, unsigned int nConcurrentIOVs)
	DataProxyProvider ()
	DataProxyProvider (const DataProxyProvider &)=delete
const ComponentDescription &	description () const
void	fillRecordsNotAllowingConcurrentIOVs (std::set< EventSetupRecordKey > &recordsNotAllowingConcurrentIOVs) const
virtual void	initConcurrentIOVs (EventSetupRecordKey const &key, unsigned int nConcurrentIOVs)
bool	isUsingRecord (const EventSetupRecordKey &key) const
KeyedProxies &	keyedProxies (const EventSetupRecordKey &iRecordKey, unsigned int iovIndex=0)
const DataProxyProvider &	operator= (const DataProxyProvider &)=delete
void	setAppendToDataLabel (const edm::ParameterSet &)
void	setDescription (const ComponentDescription &iDescription)
std::set< EventSetupRecordKey >	usingRecords () const
virtual	~DataProxyProvider () noexcept(false)

Additional Inherited Members
Static Public Member Functions inherited from edm::eventsetup::DataProxyProvider
static void	prevalidate (ConfigurationDescriptions &)
Protected Types inherited from edm::ESProxyFactoryProducer
using	EventSetupRecordKey = eventsetup::EventSetupRecordKey
Protected Types inherited from edm::eventsetup::DataProxyProvider
using	KeyedProxiesVector = std::vector< std::pair< DataKey, std::shared_ptr< DataProxy > >>
Protected Member Functions inherited from L1ConfigOnlineProdBase< L1RCTParametersRcd, L1RCTParameters >
bool	getObjectKey (const L1RCTParametersRcd &record, std::string &objectKey)
Protected Member Functions inherited from edm::ESProducer
ESConsumesInfo *	consumesInfoPushBackNew ()
unsigned int	consumesInfoSize () const
template<typename CallbackT , typename TList , typename TRecord >
void	registerProducts (std::shared_ptr< std::pair< unsigned int, std::shared_ptr< CallbackT >>> iCallback, const TList *, const TRecord *iRecord, const es::Label &iLabel)
template<typename CallbackT , typename TRecord >
void	registerProducts (std::shared_ptr< std::pair< unsigned int, std::shared_ptr< CallbackT >>>, const eventsetup::produce::Null *, const TRecord *, const es::Label &)
template<typename T >
auto	setWhatProduced (T *iThis, const es::Label &iLabel={})
template<typename T >
auto	setWhatProduced (T *iThis, const char *iLabel)
template<typename T >
auto	setWhatProduced (T *iThis, const std::string &iLabel)
template<typename T , typename TDecorator >
auto	setWhatProduced (T *iThis, const TDecorator &iDec, const es::Label &iLabel={})
template<typename T , typename TReturn , typename TRecord >
auto	setWhatProduced (T *iThis, TReturn(T::*iMethod)(const TRecord &), const es::Label &iLabel={})
template<typename T , typename TReturn , typename TRecord , typename TDecorator >
auto	setWhatProduced (T *iThis, TReturn(T ::*iMethod)(const TRecord &), const TDecorator &iDec, const es::Label &iLabel={})
template<typename TFunc >
auto	setWhatProduced (TFunc &&func, const es::Label &iLabel={})
template<typename TReturn , typename TRecord , typename TFunc , typename TDecorator >
ESConsumesCollectorT< TRecord >	setWhatProduced (TFunc &&func, TDecorator &&iDec, const es::Label &iLabel={})
void	usesResources (std::vector< std::string > const &)
Protected Member Functions inherited from edm::ESProxyFactoryProducer
template<class TFactory >
void	registerFactory (std::unique_ptr< TFactory > iFactory, const std::string &iLabel=std::string())
virtual void	registerFactoryWithKey (const EventSetupRecordKey &iRecord, std::unique_ptr< eventsetup::ProxyFactoryBase > iFactory, const std::string &iLabel=std::string())
KeyedProxiesVector	registerProxies (const EventSetupRecordKey &, unsigned int iovIndex) override
Protected Member Functions inherited from edm::eventsetup::DataProxyProvider
template<class T >
void	usingRecord ()
void	usingRecordWithKey (const EventSetupRecordKey &key)
Protected Attributes inherited from L1ConfigOnlineProdBase< L1RCTParametersRcd, L1RCTParameters >
std::optional< edm::ESConsumesCollectorT< L1RCTParametersRcd > >	m_consumesCollector
bool	m_copyFromCondDB
cond::persistency::Session	m_dbSession
bool	m_forceGeneration
l1t::OMDSReader	m_omdsReader

Detailed Description

Description: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 34 of file L1RCTParametersOnlineProd.cc.

Constructor & Destructor Documentation

L1RCTParametersOnlineProd()

L1RCTParametersOnlineProd::L1RCTParametersOnlineProd

(

const edm::ParameterSet &

iConfig

)

Definition at line 63 of file L1RCTParametersOnlineProd.cc.

    : L1ConfigOnlineProdBase<L1RCTParametersRcd, L1RCTParameters>(iConfig) {
  //the following line is needed to tell the framework what
  // data is being produced

  //now do what ever other initialization is needed
}

~L1RCTParametersOnlineProd()

L1RCTParametersOnlineProd::~L1RCTParametersOnlineProd ( )

override

Definition at line 71 of file L1RCTParametersOnlineProd.cc.

                                                      {
  // do anything here that needs to be done at desctruction time
  // (e.g. close files, deallocate resources etc.)
}

Member Function Documentation

fillScaleFactors()

void L1RCTParametersOnlineProd::fillScaleFactors	(	const l1t::OMDSReader::QueryResults &	results,
		std::vector< double > &	output,
		int	nfactor = 1
	)

Definition at line 428 of file L1RCTParametersOnlineProd.cc.

References mps_fire::i, LEDCalibrationChannels::ieta, cms::alpakatools::config::maxBin, convertSQLitetoXML_cfg::output, push_back(), and mysort::results.

Referenced by newObject().

                                                               {
  if ((nfactors < 1) || (nfactors > 6)) {
    edm::LogError("L1-O2O") << "invalid number of factors in scale factors fill";
    return;
  }

  std::vector<std::string> scaleFactorQuery6Strings;
  scaleFactorQuery6Strings.push_back("SCALEFACTOR");
  scaleFactorQuery6Strings.push_back("SF2");
  scaleFactorQuery6Strings.push_back("SF3");
  scaleFactorQuery6Strings.push_back("SF4");
  scaleFactorQuery6Strings.push_back("SF5");
  scaleFactorQuery6Strings.push_back("SF6");
  // Store scale factors in temporary array to get ordering right.
  // Reserve space for 100 bins.

  //  static const int reserve = 100 ;
  std::vector<double> sfTmp[100];
  /*  
  for( int i = 0 ; i < reserve ; ++i )
    {
      sfTmp[ i ] = 0. ;
    }
  */
  short maxBin = 0;
  for (int i = 0; i < results.numberRows(); ++i) {
    double sf[6];
    for (int nf = 0; nf < nfactors; nf++) {
      results.fillVariableFromRow(scaleFactorQuery6Strings.at(nf), i, sf[nf]);
    }
    short ieta = 0;
    results.fillVariableFromRow("FK_RCT_ETA", i, ieta);

    for (int nf = 0; nf < nfactors; nf++)
      //      sfTmp[ ieta-1 ] = sf ; // eta bins start at 1.
      sfTmp[ieta - 1].push_back(sf[nf]);

    if (ieta > maxBin) {
      maxBin = ieta;
    }
  }

  for (short i = 0; i < maxBin; ++i) {
    for (short nf = 0; nf < nfactors; nf++)
      output.push_back(sfTmp[i].at(nf));
  }
}

newObject()

std::unique_ptr< L1RCTParameters > L1RCTParametersOnlineProd::newObject ( const std::string &  objectKey )

overridevirtual

Implements L1ConfigOnlineProdBase< L1RCTParametersRcd, L1RCTParameters >.

Definition at line 76 of file L1RCTParametersOnlineProd.cc.

References l1t::OMDSReader::basicQuery(), caloStage1RCTLuts_cff::eActivityCut, L1StartupConfig_cff::eGammaLSB, caloStage1RCTLuts_cff::eicIsolationThreshold, caloStage1RCTLuts_cff::eMaxForFGCut, caloStage1RCTLuts_cff::eMaxForHoECut, caloStage1RCTLuts_cff::eMinForFGCut, caloStage1RCTLuts_cff::eMinForHoECut, fillScaleFactors(), l1t::OMDSReader::QueryResults::fillVariable(), caloStage1RCTLuts_cff::hActivityCut, caloStage1RCTLuts_cff::hMinForHoECut, caloStage1RCTLuts_cff::hOeCut, createfilelist::int, L1StartupConfig_cff::jetMETLSB, L1ConfigOnlineProdBase< L1RCTParametersRcd, L1RCTParameters >::m_omdsReader, L1RCTConfig_cfi::noiseVetoHB, L1RCTConfig_cfi::noiseVetoHEminus, L1RCTConfig_cfi::noiseVetoHEplus, l1t::OMDSReader::QueryResults::numberRows(), l1t::OMDSReader::QueryResults::queryFailed(), l1t::OMDSReader::singleAttribute(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                              {
  std::string rctSchema = "CMS_RCT";
  const l1t::OMDSReader::QueryResults paremKeyResults = m_omdsReader.singleAttribute(objectKey);

  // ~~~~~~~~~ Cut values ~~~~~~~~~

  // select egamma_lsb, jetmet_lsb, e_min_for_fg_cut, e_max_for_fg_cut,
  // h_ove r_e_cut, e_min_for_h_over_e_cut, e_max_for_h_over_e_cut,
  // h_min_for_h_over_e_cut, e_activity_cut, h_activity_cut,
  // eic_isolation_threshold, jsc_quiet_threshold_barrel,
  // jsc_quiet_threshold_endcap, noisevetohb, noisevetoheplus,
  // noisevetoheminus from parem_conf where parem_conf.parem_key =
  // (select rct_parameter from rct_conf where rct_conf.rct_key =
  // 'rct_cmssw_def');

  std::vector<std::string> queryStrings;
  queryStrings.push_back("EGAMMA_LSB");
  queryStrings.push_back("JETMET_LSB");
  queryStrings.push_back("E_MIN_FOR_FG_CUT");
  queryStrings.push_back("E_MAX_FOR_FG_CUT");
  queryStrings.push_back("H_OVER_E_CUT");
  queryStrings.push_back("E_MIN_FOR_H_OVER_E_CUT");
  queryStrings.push_back("E_MAX_FOR_H_OVER_E_CUT");
  queryStrings.push_back("H_MIN_FOR_H_OVER_E_CUT");
  queryStrings.push_back("E_ACTIVITY_CUT");
  queryStrings.push_back("H_ACTIVITY_CUT");
  queryStrings.push_back("EIC_ISOLATION_THRESHOLD");
  queryStrings.push_back("JSC_QUIET_THRESHOLD_BARREL");
  queryStrings.push_back("JSC_QUIET_THRESHOLD_ENDCAP");
  queryStrings.push_back("NOISEVETOHB");
  queryStrings.push_back("NOISEVETOHEPLUS");
  queryStrings.push_back("NOISEVETOHEMINUS");
  queryStrings.push_back("USECORR");
  l1t::OMDSReader::QueryResults paremResults =
      m_omdsReader.basicQuery(queryStrings, rctSchema, "PAREM_CONF", "PAREM_CONF.PAREM_KEY", paremKeyResults);

  if (paremResults.queryFailed() || paremResults.numberRows() != 1)  // check query successful
  {
    edm::LogError("L1-O2O") << "Problem with L1RCTParameters key.";
    return std::unique_ptr<L1RCTParameters>();
  }

  double eGammaLSB, jetMETLSB, eMinForFGCut, eMaxForFGCut, hOeCut;
  double eMinForHoECut, eMaxForHoECut, hMinForHoECut;
  double eActivityCut, hActivityCut;
  double jscQuietThreshBarrel, jscQuietThreshEndcap, eicIsolationThreshold;
  bool noiseVetoHB, noiseVetoHEplus, noiseVetoHEminus, useCorr;

  paremResults.fillVariable("EGAMMA_LSB", eGammaLSB);
  paremResults.fillVariable("JETMET_LSB", jetMETLSB);
  paremResults.fillVariable("E_MIN_FOR_FG_CUT", eMinForFGCut);
  paremResults.fillVariable("E_MAX_FOR_FG_CUT", eMaxForFGCut);
  paremResults.fillVariable("H_OVER_E_CUT", hOeCut);
  paremResults.fillVariable("E_MIN_FOR_H_OVER_E_CUT", eMinForHoECut);
  paremResults.fillVariable("E_MAX_FOR_H_OVER_E_CUT", eMaxForHoECut);
  paremResults.fillVariable("H_MIN_FOR_H_OVER_E_CUT", hMinForHoECut);
  paremResults.fillVariable("E_ACTIVITY_CUT", eActivityCut);
  paremResults.fillVariable("H_ACTIVITY_CUT", hActivityCut);
  paremResults.fillVariable("JSC_QUIET_THRESHOLD_BARREL", jscQuietThreshBarrel);
  paremResults.fillVariable("JSC_QUIET_THRESHOLD_ENDCAP", jscQuietThreshEndcap);
  paremResults.fillVariable("EIC_ISOLATION_THRESHOLD", eicIsolationThreshold);
  paremResults.fillVariable("NOISEVETOHB", noiseVetoHB);
  paremResults.fillVariable("NOISEVETOHEPLUS", noiseVetoHEplus);
  paremResults.fillVariable("NOISEVETOHEMINUS", noiseVetoHEminus);
  paremResults.fillVariable("USECORR", useCorr);
  //      std::cout << "eGammaLSB = " << eGammaLSB << std::endl ;
  //      std::cout << "jetMETLSB = " << jetMETLSB << std::endl ;
  //      std::cout << "eMinForFGCut = " << eMinForFGCut << std::endl ;
  //      std::cout << "eMaxForFGCut = " << eMaxForFGCut << std::endl ;
  //      std::cout << "hOeCut = " << hOeCut << std::endl ;
  //      std::cout << "eMinForHoECut = " << eMinForHoECut << std::endl ;
  //      std::cout << "eMaxForHoECut = " << eMaxForHoECut << std::endl ;
  //      std::cout << "hMinForHoECut = " << hMinForHoECut << std::endl ;
  //      std::cout << "eActivityCut = " << eActivityCut << std::endl ;
  //      std::cout << "hActivityCut = " << hActivityCut << std::endl ;
  //      std::cout << "eicIsolationThreshold = " << eicIsolationThreshold << std::endl ;
  //      std::cout << "jscQuietThreshBarrel = " << jscQuietThreshBarrel << std::endl ;
  //      std::cout << "jscQuietThreshEndcap = " << jscQuietThreshEndcap << std::endl ;
  //      std::cout << "noiseVetoHB = " << noiseVetoHB << std::endl ;
  //      std::cout << "noiseVetoHEplus = " << noiseVetoHEplus << std::endl ;
  //      std::cout << "noiseVetoHEminus = " << noiseVetoHEminus << std::endl ;

  // ~~~~~~~~~ EGamma ECAL scale factors ~~~~~~~~~

  // select scalefactor, fk_rct_eta from egamma_ecal_scalefactor where
  // egamma_ecal_scalefactor.fk_version =
  // (select egamma_ecal from parem_conf where parem_conf.parem_key =
  // (select rct_parameter from rct_conf where rct_conf.rct_key=
  // 'rct_cmssw_def'));

  std::vector<std::string> scaleFactorQueryStrings;
  scaleFactorQueryStrings.push_back("SCALEFACTOR");
  scaleFactorQueryStrings.push_back("FK_RCT_ETA");

  l1t::OMDSReader::QueryResults egammaEcalResults = m_omdsReader.basicQuery(
      scaleFactorQueryStrings,
      rctSchema,
      "EGAMMA_ECAL_SCALEFACTOR",
      "EGAMMA_ECAL_SCALEFACTOR.FK_VERSION",
      m_omdsReader.basicQuery("EGAMMA_ECAL", rctSchema, "PAREM_CONF", "PAREM_CONF.PAREM_KEY", paremKeyResults));

  if (egammaEcalResults.queryFailed())  // check query successful
  {
    edm::LogError("L1-O2O") << "Problem with EgammaEcal key.";
    return std::unique_ptr<L1RCTParameters>();
  }

  //      std::cout << "egammaEcal " ;
  std::vector<double> egammaEcalScaleFactors;
  fillScaleFactors(egammaEcalResults, egammaEcalScaleFactors);

  // ~~~~~~~~~ EGamma HCAL scale factors ~~~~~~~~~

  // select scalefactor, fk_rct_eta from egamma_hcal_scalefactor where
  // egamma_hcal_scalefactor.fk_version =
  // (select egamma_hcal from parem_conf where parem_conf.parem_key =
  // (select rct_parameter from rct_conf where rct_conf.rct_key=
  // 'rct_cmssw_def'));

  l1t::OMDSReader::QueryResults egammaHcalResults = m_omdsReader.basicQuery(
      scaleFactorQueryStrings,
      rctSchema,
      "EGAMMA_HCAL_SCALEFACTOR",
      "EGAMMA_HCAL_SCALEFACTOR.FK_VERSION",
      m_omdsReader.basicQuery("EGAMMA_HCAL", rctSchema, "PAREM_CONF", "PAREM_CONF.PAREM_KEY", paremKeyResults));

  if (egammaHcalResults.queryFailed())  // check query successful
  {
    edm::LogError("L1-O2O") << "Problem with EgammaHcal key.";
    return std::unique_ptr<L1RCTParameters>();
  }

  //      std::cout << "egammaHcal " ;
  std::vector<double> egammaHcalScaleFactors;
  fillScaleFactors(egammaHcalResults, egammaHcalScaleFactors);

  // ~~~~~~~~~ JetMET ECAL scale factors ~~~~~~~~~

  // select scalefactor, fk_rct_eta from jetmet_ecal_scalefactor where
  // jetmet_ecal_scalefactor.fk_version =
  // (select jetmet_ecal from parem_conf where parem_conf.parem_key =
  // (select rct_parameter from rct_conf where rct_conf.rct_key=
  // 'rct_cmssw_def'));

  l1t::OMDSReader::QueryResults jetmetEcalResults = m_omdsReader.basicQuery(
      scaleFactorQueryStrings,
      rctSchema,
      "JETMET_ECAL_SCALEFACTOR",
      "JETMET_ECAL_SCALEFACTOR.FK_VERSION",
      m_omdsReader.basicQuery("JETMET_ECAL", rctSchema, "PAREM_CONF", "PAREM_CONF.PAREM_KEY", paremKeyResults));

  if (jetmetEcalResults.queryFailed())  // check query successful
  {
    edm::LogError("L1-O2O") << "Problem with JetmetEcal key.";
    return std::unique_ptr<L1RCTParameters>();
  }

  //      std::cout << "jetmetEcal " ;
  std::vector<double> jetmetEcalScaleFactors;
  fillScaleFactors(jetmetEcalResults, jetmetEcalScaleFactors);

  // ~~~~~~~~~ JetMET HCAL scale factors ~~~~~~~~~

  // select scalefactor, fk_rct_eta from jetmet_hcal_scalefactor where
  // jetmet_hcal_scalefactor.fk_version =
  // (select jetmet_hcal from parem_conf where parem_conf.parem_key =
  // (select rct_parameter from rct_conf where rct_conf.rct_key=
  // 'rct_cmssw_def'));

  l1t::OMDSReader::QueryResults jetmetHcalResults = m_omdsReader.basicQuery(
      scaleFactorQueryStrings,
      rctSchema,
      "JETMET_HCAL_SCALEFACTOR",
      "JETMET_HCAL_SCALEFACTOR.FK_VERSION",
      m_omdsReader.basicQuery("JETMET_HCAL", rctSchema, "PAREM_CONF", "PAREM_CONF.PAREM_KEY", paremKeyResults));

  if (jetmetHcalResults.queryFailed())  // check query successful
  {
    edm::LogError("L1-O2O") << "Problem with JetmetHcal key.";
    return std::unique_ptr<L1RCTParameters>();
  }

  //      std::cout << "jetmetHcal " ;
  std::vector<double> jetmetHcalScaleFactors;
  fillScaleFactors(jetmetHcalResults, jetmetHcalScaleFactors);

  //Lindsay variables

  /*
    std::vector< std::vector< double > > hcalCalibScaleFactors ;
       std::vector< std::vector< double > > hcalCalibHighScaleFactors ;
       std::vector< std::vector< double > > ecalCalibScaleFactors ;
       std::vector< std::vector< double > > crossTermsScaleFactors ;
     */
  std::vector<double> lowHoverE_smear, highHoverE_smear;
  std::vector<double> hcalCalibScaleFactors, ecalCalibScaleFactors;
  std::vector<double> hcalCalibHighScaleFactors, crossTermsScaleFactors;

  if (useCorr) {  //lindsay corrections

    std::vector<std::string> scaleFactorQuery3Strings;
    scaleFactorQuery3Strings.push_back("SCALEFACTOR");
    scaleFactorQuery3Strings.push_back("SF2");
    scaleFactorQuery3Strings.push_back("SF3");
    scaleFactorQuery3Strings.push_back("FK_RCT_ETA");

    l1t::OMDSReader::QueryResults hcalCalibResults = m_omdsReader.basicQuery(
        scaleFactorQuery3Strings,
        rctSchema,
        "HCAL_CALIB_FACTOR",
        "HCAL_CALIB_FACTOR.VERSION",
        m_omdsReader.basicQuery(
            "HCAL_CALIB_VERSION", rctSchema, "PAREM_CONF", "PAREM_CONF.PAREM_KEY", paremKeyResults));

    if (hcalCalibResults.queryFailed()) {  // check query successful

      edm::LogError("L1-O2O") << "Problem with JetmetHcal key.";
      return std::unique_ptr<L1RCTParameters>();
    }

    //      std::cout << "jetmetHcal " ;

    fillScaleFactors(hcalCalibResults, hcalCalibScaleFactors, 3);

    l1t::OMDSReader::QueryResults hcalCalibHighResults = m_omdsReader.basicQuery(
        scaleFactorQuery3Strings,
        rctSchema,
        "HCAL_CALIB_HIGH_FACTOR",
        "HCAL_CALIB_HIGH_FACTOR.VERSION",
        m_omdsReader.basicQuery(
            "HCAL_CALIB_HIGH_VERSION", rctSchema, "PAREM_CONF", "PAREM_CONF.PAREM_KEY", paremKeyResults));

    if (hcalCalibHighResults.queryFailed())  // check query successful
    {
      edm::LogError("L1-O2O") << "Problem with hcalHigh key.";
      return std::unique_ptr<L1RCTParameters>();
    }

    fillScaleFactors(hcalCalibHighResults, hcalCalibHighScaleFactors, 3);

    l1t::OMDSReader::QueryResults ecalCalibResults = m_omdsReader.basicQuery(
        scaleFactorQuery3Strings,
        rctSchema,
        "ECAL_CALIB_FACTOR",
        "ECAL_CALIB_FACTOR.VERSION",
        m_omdsReader.basicQuery(
            "ECAL_CALIB_VERSION", rctSchema, "PAREM_CONF", "PAREM_CONF.PAREM_KEY", paremKeyResults));

    if (ecalCalibResults.queryFailed())  // check query successful
    {
      edm::LogError("L1-O2O") << "Problem with ecal calib key.";
      return std::unique_ptr<L1RCTParameters>();
    }

    fillScaleFactors(ecalCalibResults, ecalCalibScaleFactors, 3);

    std::vector<std::string> scaleFactorQuery6Strings;
    scaleFactorQuery6Strings.push_back("SCALEFACTOR");
    scaleFactorQuery6Strings.push_back("SF2");
    scaleFactorQuery6Strings.push_back("SF3");
    scaleFactorQuery6Strings.push_back("SF4");
    scaleFactorQuery6Strings.push_back("SF5");
    scaleFactorQuery6Strings.push_back("SF6");
    scaleFactorQuery6Strings.push_back("FK_RCT_ETA");
    l1t::OMDSReader::QueryResults crossTermResults = m_omdsReader.basicQuery(
        scaleFactorQuery6Strings,
        rctSchema,
        "CROSS_TERMS_FACTOR",
        "CROSS_TERMS_FACTOR.VERSION",
        m_omdsReader.basicQuery(
            "CROSS_TERMS_VERSION", rctSchema, "PAREM_CONF", "PAREM_CONF.PAREM_KEY", paremKeyResults));

    if (crossTermResults.queryFailed())  // check query successful
    {
      edm::LogError("L1-O2O") << "Problem with crossTerms key.";
      return std::unique_ptr<L1RCTParameters>();
    }

    fillScaleFactors(crossTermResults, crossTermsScaleFactors, 6);

    l1t::OMDSReader::QueryResults hoveresmearhighResults = m_omdsReader.basicQuery(
        scaleFactorQueryStrings,
        rctSchema,
        "H_OVER_E_SMEAR_HIGH_FACTOR",
        "H_OVER_E_SMEAR_HIGH_FACTOR.FK_VERSION",
        m_omdsReader.basicQuery(
            "H_OVER_E_SMEAR_HIGH_VERSION", rctSchema, "PAREM_CONF", "PAREM_CONF.PAREM_KEY", paremKeyResults));

    if (hoveresmearhighResults.queryFailed())  // check query successful
    {
      edm::LogError("L1-O2O") << "Problem with low h over e smear key.";
      return std::unique_ptr<L1RCTParameters>();
    }

    //      std::cout << "egammaEcal " ;
    fillScaleFactors(hoveresmearhighResults, highHoverE_smear);

    l1t::OMDSReader::QueryResults hoveresmearlowResults = m_omdsReader.basicQuery(
        scaleFactorQueryStrings,
        rctSchema,
        "H_OVER_E_SMEAR_LOW_FACTOR",
        "H_OVER_E_SMEAR_LOW_FACTOR.FK_VERSION",
        m_omdsReader.basicQuery(
            "H_OVER_E_SMEAR_LOW_VERSION", rctSchema, "PAREM_CONF", "PAREM_CONF.PAREM_KEY", paremKeyResults));

    if (hoveresmearlowResults.queryFailed())  // check query successful
    {
      edm::LogError("L1-O2O") << "Problem with low h over e smear key.";
      return std::unique_ptr<L1RCTParameters>();
    }

    //      std::cout << "egammaEcal " ;
    fillScaleFactors(hoveresmearlowResults, lowHoverE_smear);
  }

  //~~~~~~~~~ Instantiate new L1RCTParameters object. ~~~~~~~~~

  // Default objects for Lindsey

  return std::make_unique<L1RCTParameters>(eGammaLSB,
                                           jetMETLSB,
                                           eMinForFGCut,
                                           eMaxForFGCut,
                                           hOeCut,
                                           eMinForHoECut,
                                           eMaxForHoECut,
                                           hMinForHoECut,
                                           eActivityCut,
                                           hActivityCut,
                                           (unsigned int)eicIsolationThreshold,
                                           (int)jscQuietThreshBarrel,
                                           (int)jscQuietThreshEndcap,
                                           noiseVetoHB,
                                           noiseVetoHEplus,
                                           noiseVetoHEminus,
                                           useCorr,  // useLindsey
                                           egammaEcalScaleFactors,
                                           egammaHcalScaleFactors,
                                           jetmetEcalScaleFactors,
                                           jetmetHcalScaleFactors,
                                           ecalCalibScaleFactors,
                                           hcalCalibScaleFactors,
                                           hcalCalibHighScaleFactors,
                                           crossTermsScaleFactors,
                                           lowHoverE_smear,
                                           highHoverE_smear);
}