#include <L1Trigger/L1GctConfigProducers/interface/L1GctConfigProducers.h>

Inheritance diagram for L1GctConfigProducers:

Public Types
typedef boost::shared_ptr < L1GctChannelMask >	ChanMaskReturnType

typedef boost::shared_ptr < L1GctJetFinderParams >	JfParamsReturnType

Public Types inherited from edm::eventsetup::DataProxyProvider
typedef std::vector< std::pair < DataKey, boost::shared_ptr < DataProxy > > >	KeyedProxies

typedef std::vector < EventSetupRecordKey >	Keys

typedef std::map < EventSetupRecordKey, KeyedProxies >	RecordProxies

Public Member Functions
std::vector< double >	etToEnergyConversion (const L1CaloGeometry *geom) const
	Legacy nonsense. More...

	L1GctConfigProducers (const edm::ParameterSet &)

ChanMaskReturnType	produceChanMask (const L1GctChannelMaskRcd &)

JfParamsReturnType	produceJfParams (const L1GctJetFinderParamsRcd &)

	~L1GctConfigProducers ()

Public Member Functions inherited from edm::ESProducer
	ESProducer ()

virtual	~ESProducer ()

Public Member Functions inherited from edm::ESProxyFactoryProducer
	ESProxyFactoryProducer ()

virtual void	newInterval (const eventsetup::EventSetupRecordKey &iRecordType, const ValidityInterval &iInterval)
	overrides DataProxyProvider method More...

virtual	~ESProxyFactoryProducer ()

Public Member Functions inherited from edm::eventsetup::DataProxyProvider
	DataProxyProvider ()

const ComponentDescription &	description () const

bool	isUsingRecord (const EventSetupRecordKey &) const

const KeyedProxies &	keyedProxies (const EventSetupRecordKey &iRecordKey) const

void	resetProxies (const EventSetupRecordKey &iRecordType)

void	resetProxiesIfTransient (const EventSetupRecordKey &iRecordType)

void	setAppendToDataLabel (const edm::ParameterSet &)

void	setDescription (const ComponentDescription &iDescription)

std::set< EventSetupRecordKey >	usingRecords () const

virtual	~DataProxyProvider ()

Private Attributes
double	m_CenJetSeed

bool	m_convertToEnergy

unsigned	m_corrFunType

unsigned	m_EtaBoundry

double	m_FwdJetSeed

double	m_htJetThresh

double	m_htLsb

std::vector< std::vector < double > >	m_jetCalibFunc

unsigned	m_metEtaMask

unsigned	m_mhtEtaMask

double	m_mhtJetThresh

double	m_rgnEtLsb

std::vector< std::vector < double > >	m_tauCalibFunc

double	m_tauIsoThresh

double	m_TauJetSeed

unsigned	m_tetEtaMask

unsigned	m_thtEtaMask

Additional Inherited Members
Static Public Member Functions inherited from edm::eventsetup::DataProxyProvider
static void	prevalidate (ConfigurationDescriptions &)

Protected Member Functions inherited from edm::ESProducer
template<typename T >
void	setWhatProduced (T *iThis, const es::Label &iLabel=es::Label())

template<typename T >
void	setWhatProduced (T iThis, const char iLabel)

template<typename T >
void	setWhatProduced (T *iThis, const std::string &iLabel)

template<typename T , typename TDecorator >
void	setWhatProduced (T *iThis, const TDecorator &iDec, const es::Label &iLabel=es::Label())

template<typename T , typename TReturn , typename TRecord >
void	setWhatProduced (T iThis, TReturn(T::iMethod)(const TRecord &), const es::Label &iLabel=es::Label())

template<typename T , typename TReturn , typename TRecord , typename TArg >
void	setWhatProduced (T iThis, TReturn(T::iMethod)(const TRecord &), const TArg &iDec, const es::Label &iLabel=es::Label())

Protected Member Functions inherited from edm::ESProxyFactoryProducer
template<class TFactory >
void	registerFactory (std::auto_ptr< TFactory > iFactory, const std::string &iLabel=std::string())

virtual void	registerFactoryWithKey (const eventsetup::EventSetupRecordKey &iRecord, std::auto_ptr< eventsetup::ProxyFactoryBase > &iFactory, const std::string &iLabel=std::string())

virtual void	registerProxies (const eventsetup::EventSetupRecordKey &iRecord, KeyedProxies &aProxyList)
	override DataProxyProvider method More...

Protected Member Functions inherited from edm::eventsetup::DataProxyProvider
void	eraseAll (const EventSetupRecordKey &iRecordKey)
	deletes all the Proxies in aStream More...

void	invalidateProxies (const EventSetupRecordKey &iRecordKey)

template<class T >
void	usingRecord ()

void	usingRecordWithKey (const EventSetupRecordKey &)

Detailed Description

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

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

Definition at line 51 of file L1GctConfigProducers.h.

Member Typedef Documentation

typedef boost::shared_ptr<L1GctChannelMask> L1GctConfigProducers::ChanMaskReturnType

Definition at line 57 of file L1GctConfigProducers.h.

typedef boost::shared_ptr<L1GctJetFinderParams> L1GctConfigProducers::JfParamsReturnType

Definition at line 56 of file L1GctConfigProducers.h.

Constructor & Destructor Documentation

L1GctConfigProducers::L1GctConfigProducers ( const edm::ParameterSet & iConfig )

Definition at line 25 of file L1GctConfigProducers.cc.

References edm::ParameterSet::getParameter(), i, m_corrFunType, m_jetCalibFunc, m_tauCalibFunc, L1GctJetFinderParams::N_CENTRAL_ETA_VALUES, L1GctJetFinderParams::NUMBER_ETA_VALUES, produceChanMask(), produceJfParams(), and edm::ESProducer::setWhatProduced().

                                                                          :
   m_rgnEtLsb(iConfig.getParameter<double>("RctRegionEtLSB")),
   m_htLsb(iConfig.getParameter<double>("GctHtLSB")),
   m_CenJetSeed(iConfig.getParameter<double>("JetFinderCentralJetSeed")),
   m_FwdJetSeed(iConfig.getParameter<double>("JetFinderForwardJetSeed")),
   m_TauJetSeed(iConfig.getParameter<double>("JetFinderCentralJetSeed")), // no separate tau jet seed yet
   m_tauIsoThresh(iConfig.getParameter<double>("TauIsoEtThreshold")),
   m_htJetThresh(iConfig.getParameter<double>("HtJetEtThreshold")),
   m_mhtJetThresh(iConfig.getParameter<double>("MHtJetEtThreshold")),
   m_EtaBoundry(7), // not programmable!
   m_corrFunType(0),
   m_convertToEnergy (iConfig.getParameter<bool>("ConvertEtValuesToEnergy")),
   m_jetCalibFunc(),
   m_tauCalibFunc(),
   m_metEtaMask(iConfig.getParameter<unsigned>("MEtEtaMask")),
   m_tetEtaMask(iConfig.getParameter<unsigned>("TEtEtaMask")),
   m_mhtEtaMask(iConfig.getParameter<unsigned>("MHtEtaMask")),
   m_thtEtaMask(iConfig.getParameter<unsigned>("HtEtaMask"))
 {
 
    //the following lines are needed to tell the framework what
    // data is being produced
    setWhatProduced(this,&L1GctConfigProducers::produceJfParams);
    setWhatProduced(this,&L1GctConfigProducers::produceChanMask);
 
    //now do what ever other initialization is needed
    std::string CalibStyle = iConfig.getParameter<std::string>("CalibrationStyle");
    edm::ParameterSet calibCoeffs;
    
    if (CalibStyle == "PowerSeries") {
      m_corrFunType = 1;
      calibCoeffs = iConfig.getParameter<edm::ParameterSet>("PowerSeriesCoefficients");
    }
    
    if (CalibStyle == "ORCAStyle") {
      m_corrFunType = 2;
      calibCoeffs = iConfig.getParameter<edm::ParameterSet>("OrcaStyleCoefficients");
    }
    
    if (CalibStyle == "Simple") {
      m_corrFunType = 3;
      calibCoeffs = iConfig.getParameter<edm::ParameterSet>("SimpleCoefficients");
    }
    
    if (CalibStyle == "PiecewiseCubic") {
      m_corrFunType = 4;
      calibCoeffs = iConfig.getParameter<edm::ParameterSet>("PiecewiseCubicCoefficients");
    }
    
    if (CalibStyle == "PF") {
      m_corrFunType = 5;
      calibCoeffs = iConfig.getParameter<edm::ParameterSet>("PFCoefficients");
    }
 
    // check 
    if (CalibStyle != "None") {
      
      // Read the coefficients from file
      // coefficients for non-tau jet corrections
      for (unsigned i=0; i<L1GctJetFinderParams::NUMBER_ETA_VALUES; ++i) {
        std::stringstream ss;
        std::string str;
        ss << "nonTauJetCalib" << i;
        ss >> str;
        m_jetCalibFunc.push_back(calibCoeffs.getParameter< std::vector<double> >(str));
      }
      // coefficients for tau jet corrections
      for (unsigned i=0; i<L1GctJetFinderParams::N_CENTRAL_ETA_VALUES; ++i) {
        std::stringstream ss;
        std::string str;
        ss << "tauJetCalib" << i;
        ss >> str;
        m_tauCalibFunc.push_back(calibCoeffs.getParameter< std::vector<double> >(str));
      }
      
    } else {
      // No corrections to be applied
      m_corrFunType = 0;  // no correction
      // Set the vector sizes to those expected by the CalibrationFunction
      m_jetCalibFunc.resize(L1GctJetFinderParams::NUMBER_ETA_VALUES);
      m_tauCalibFunc.resize(L1GctJetFinderParams::N_CENTRAL_ETA_VALUES);
    }
 
    edm::LogWarning("L1GctConfig") << "Calibration Style option " << CalibStyle << std::endl;
    
 }

L1GctConfigProducers::~L1GctConfigProducers ( )

Definition at line 113 of file L1GctConfigProducers.cc.

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

Member Function Documentation

std::vector< double > L1GctConfigProducers::etToEnergyConversion ( const L1CaloGeometry * geom ) const

Legacy nonsense.

Calculate Et-to-energy conversion factors for eta bins

Definition at line 173 of file L1GctConfigProducers.cc.

References L1CaloGeometry::etaBinCenter(), create_public_lumi_plots::exp, and query::result.

Referenced by produceJfParams().

                                      {
   //  L1CaloGeometry* geom = new L1CaloGeometry();
   std::vector<double> result;
   // Factors for central eta bins
   for (unsigned ieta=0; ieta<7; ieta++) {
     double bineta = geom->etaBinCenter(ieta, true);
     double factor = 0.5*(exp(bineta)+exp(-bineta)); // Conversion from eta to cosec(theta)
     result.push_back(factor);
   }
   // Factors for forward eta bins
   for (unsigned ieta=0; ieta<4; ieta++) {
     double bineta = geom->etaBinCenter(ieta, false);
     double factor = 0.5*(exp(bineta)+exp(-bineta)); // Conversion from eta to cosec(theta)
     result.push_back(factor);
   }
   return result;
 }

L1GctConfigProducers::ChanMaskReturnType L1GctConfigProducers::produceChanMask ( const L1GctChannelMaskRcd & )

Definition at line 153 of file L1GctConfigProducers.cc.

References m_metEtaMask, m_mhtEtaMask, m_tetEtaMask, m_thtEtaMask, L1GctChannelMask::maskMissingEt(), L1GctChannelMask::maskMissingHt(), L1GctChannelMask::maskTotalEt(), and L1GctChannelMask::maskTotalHt().

Referenced by L1GctConfigProducers().

                                                                 {
 
   L1GctChannelMask* mask = new L1GctChannelMask;
 
   for (unsigned ieta=0; ieta<22; ++ieta) {
     if (((m_metEtaMask>>ieta)&0x1)==1) mask->maskMissingEt(ieta);
     if (((m_tetEtaMask>>ieta)&0x1)==1) mask->maskTotalEt(ieta);
     if (((m_mhtEtaMask>>ieta)&0x1)==1) mask->maskMissingHt(ieta);
     if (((m_thtEtaMask>>ieta)&0x1)==1) mask->maskTotalHt(ieta);
   }
 
   return boost::shared_ptr<L1GctChannelMask>(mask);
 
 }

L1GctConfigProducers::JfParamsReturnType L1GctConfigProducers::produceJfParams ( const L1GctJetFinderParamsRcd & aRcd )

Definition at line 124 of file L1GctConfigProducers.cc.

References etToEnergyConversion(), relativeConstraints::geom, edm::eventsetup::EventSetupRecord::get(), edm::eventsetup::DependentRecordImplementation< RecordT, ListT >::getRecord(), m_CenJetSeed, m_convertToEnergy, m_corrFunType, m_EtaBoundry, m_FwdJetSeed, m_htJetThresh, m_htLsb, m_jetCalibFunc, m_mhtJetThresh, m_rgnEtLsb, m_tauCalibFunc, m_tauIsoThresh, and m_TauJetSeed.

Referenced by L1GctConfigProducers().

 {
   // get geometry
   const L1CaloGeometryRecord& geomRcd = aRcd.getRecord< L1CaloGeometryRecord >() ;
   edm::ESHandle< L1CaloGeometry > geom ;
   geomRcd.get( geom ) ;
   
   // construct jet finder params object
   boost::shared_ptr<L1GctJetFinderParams> pL1GctJetFinderParams =
     boost::shared_ptr<L1GctJetFinderParams> (new L1GctJetFinderParams(m_rgnEtLsb,
                                       m_htLsb,
                                       m_CenJetSeed,
                                       m_FwdJetSeed,
                                       m_TauJetSeed,
                                       m_tauIsoThresh,
                                       m_htJetThresh,
                                       m_mhtJetThresh,
                                       m_EtaBoundry,
                                       m_corrFunType,
                                       m_jetCalibFunc,
                                       m_tauCalibFunc,
                                       m_convertToEnergy,
                                       etToEnergyConversion(geom.product())) );
   
   return pL1GctJetFinderParams ;
 
 }