#include <SCRegressionCalculator.h>

Public Member Functions
float	getCorrection (const reco::SuperCluster &) const

std::pair< float, float >	getCorrectionWithErrors (const reco::SuperCluster &) const

	SCRegressionCalculator (const edm::ParameterSet &)

void	update (const edm::EventSetup &)

std::unique_ptr< VarCalc > &	varCalc ()

Private Attributes
edm::ESHandle< GBRForest >	eb_corr

std::string	eb_corr_name

edm::ESHandle< GBRForest >	eb_err

std::string	eb_err_name

edm::ESHandle< GBRForest >	ee_corr

std::string	ee_corr_name

edm::ESHandle< GBRForest >	ee_err

std::string	ee_err_name

const GBRWrapperRcd *	gbr_record

std::unique_ptr< VarCalc >	var_calc

Detailed Description

template<class VarCalc>
class SCRegressionCalculator< VarCalc >

Definition at line 21 of file SCRegressionCalculator.h.

Constructor & Destructor Documentation

template<class VarCalc >

SCRegressionCalculator< VarCalc >::SCRegressionCalculator ( const edm::ParameterSet & conf )

Definition at line 41 of file SCRegressionCalculator.h.

References SCRegressionCalculator< VarCalc >::eb_corr_name, SCRegressionCalculator< VarCalc >::eb_err_name, SCRegressionCalculator< VarCalc >::ee_corr_name, SCRegressionCalculator< VarCalc >::ee_err_name, edm::ParameterSet::existsAs(), edm::ParameterSet::getParameter(), AlCaHLTBitMon_QueryRunRegistry::string, and SCRegressionCalculator< VarCalc >::var_calc.

                                                     :
   gbr_record(NULL) {
   var_calc.reset(new VarCalc());
   eb_corr_name = conf.getParameter<std::string>("regressionKeyEB");
   ee_corr_name = conf.getParameter<std::string>("regressionKeyEE");
   if( conf.existsAs<std::string>("uncertaintyKeyEB") )
     eb_err_name = conf.getParameter<std::string>("uncertaintyKeyEB");
   if( conf.existsAs<std::string>("uncertaintyKeyEE") )
     ee_err_name = conf.getParameter<std::string>("uncertaintyKeyEE");
   
 }

Member Function Documentation

template<class VarCalc >

float SCRegressionCalculator< VarCalc >::getCorrection ( const reco::SuperCluster & sc ) const

Definition at line 74 of file SCRegressionCalculator.h.

References EcalBarrel, EcalEndcap, SiPixelRawToDigiRegional_cfi::inputs, and reco::SuperCluster::seed().

                                                 {
   std::vector<float> inputs;
   var_calc->set(sc,inputs);
   switch( sc.seed()->seed().subdetId() ) {
   case EcalSubdetector::EcalBarrel:
     return eb_corr->GetResponse(inputs.data());
     break;
   case EcalSubdetector::EcalEndcap:
     return ee_corr->GetResponse(inputs.data());
     break;
   }
   return -1.0f;
 }

template<class VarCalc >

std::pair< float, float > SCRegressionCalculator< VarCalc >::getCorrectionWithErrors ( const reco::SuperCluster & sc ) const

Definition at line 90 of file SCRegressionCalculator.h.

References EcalBarrel, EcalEndcap, f, SiPixelRawToDigiRegional_cfi::inputs, and reco::SuperCluster::seed().

                                                           {
   std::vector<float> inputs;
   var_calc->set(sc,inputs);
   switch( sc.seed()->seed().subdetId() ) {
   case EcalSubdetector::EcalBarrel:
     return std::make_pair( eb_corr->GetResponse(inputs.data()),
                eb_err->GetResponse(inputs.data())  );
     break;
   case EcalSubdetector::EcalEndcap:
     return std::make_pair( ee_corr->GetResponse(inputs.data()),
                ee_err->GetResponse(inputs.data())  );
     break;
   }
   return std::make_pair(-1.0f,-1.0f);
 }

template<class VarCalc >

void SCRegressionCalculator< VarCalc >::update ( const edm::EventSetup & es )

Definition at line 55 of file SCRegressionCalculator.h.

References edm::eventsetup::EventSetupRecord::cacheIdentifier(), edm::EventSetup::get(), and edm::eventsetup::EventSetupRecord::get().

                                 {
   var_calc->update(es);
   const GBRWrapperRcd& gbrfrom_es = es.get<GBRWrapperRcd>();
   if( !gbr_record || 
       gbrfrom_es.cacheIdentifier() != gbr_record->cacheIdentifier() ) {
     gbr_record = &gbrfrom_es;
     gbr_record->get(eb_corr_name.c_str(),eb_corr);    
     gbr_record->get(ee_corr_name.c_str(),ee_corr);
     if( eb_err_name.size() ) {
       gbr_record->get(eb_err_name.c_str(),eb_err);      
     }
     if( ee_err_name.size() ) {
       gbr_record->get(ee_err_name.c_str(),ee_err);
     }
   }
 }