#include <ElectronEnergyCalibratorRun2.h>

Public Member Functions
void	calibrate (SimpleElectron &electron, edm::StreamID const &id=edm::StreamID::invalidStreamID()) const

void	calibrate (reco::GsfElectron &electron, unsigned int runNumber, edm::StreamID const &id=edm::StreamID::invalidStreamID()) const

	ElectronEnergyCalibratorRun2 ()

	ElectronEnergyCalibratorRun2 (EpCombinationTool &combinator, bool isMC, bool synchronization, std::string)

void	initPrivateRng (TRandom *rnd)

	~ElectronEnergyCalibratorRun2 ()

Protected Member Functions
double	gauss (edm::StreamID const &id) const

Protected Attributes
EnergyScaleCorrection_class	_correctionRetriever

EpCombinationTool *	epCombinationTool_

bool	isMC_

TRandom *	rng_

bool	synchronization_

Detailed Description

Definition at line 13 of file ElectronEnergyCalibratorRun2.h.

Constructor & Destructor Documentation

ElectronEnergyCalibratorRun2::ElectronEnergyCalibratorRun2 ( )

inline

Definition at line 16 of file ElectronEnergyCalibratorRun2.h.

References calibrate(), metsig::electron, initPrivateRng(), edm::StreamID::invalidStreamID(), convertSQLiteXML::runNumber, AlCaHLTBitMon_QueryRunRegistry::string, and ~ElectronEnergyCalibratorRun2().

16 {}

ElectronEnergyCalibratorRun2::ElectronEnergyCalibratorRun2	(	EpCombinationTool &	combinator,
		bool	isMC,
		bool	synchronization,
		std::string	correctionFile
	)

Definition at line 7 of file ElectronEnergyCalibratorRun2.cc.

References _correctionRetriever, EnergyScaleCorrection_class::doScale, EnergyScaleCorrection_class::doSmearings, and isMC_.

                                  :
   epCombinationTool_(&combinator),
   isMC_(isMC), synchronization_(synchronization),
   rng_(0),
   _correctionRetriever(correctionFile) // here is opening the files and reading the corrections
 {
   if(isMC_) {
     _correctionRetriever.doScale = false; 
     _correctionRetriever.doSmearings = true;
   } else {
     _correctionRetriever.doScale = true; 
     _correctionRetriever.doSmearings = false;
   }
 }

ElectronEnergyCalibratorRun2::~ElectronEnergyCalibratorRun2 ( )

Definition at line 26 of file ElectronEnergyCalibratorRun2.cc.

Referenced by ElectronEnergyCalibratorRun2().

27 {}

Member Function Documentation

void ElectronEnergyCalibratorRun2::calibrate	(	SimpleElectron &	electron,
		edm::StreamID const &	id = `edm::StreamID::invalidStreamID()`
	)		const

Correct this electron. StreamID is needed when used with CMSSW Random Number Generator

Definition at line 40 of file ElectronEnergyCalibratorRun2.cc.

References _correctionRetriever, funct::abs(), EpCombinationTool::combine(), corr, epCombinationTool_, stringResolutionProvider_cfi::et, gauss(), SimpleElectron::getEta(), SimpleElectron::getNewEnergy(), SimpleElectron::getNewEnergyError(), SimpleElectron::getR9(), SimpleElectron::getRunNumber(), EnergyScaleCorrection_class::getSmearingSigma(), SimpleElectron::isEB(), SimpleElectron::isMC(), isMC_, Scenarios_cff::scale, EnergyScaleCorrection_class::ScaleCorrection(), SimpleElectron::setNewEnergy(), and SimpleElectron::setNewEnergyError().

Referenced by calibrate(), ElectronEnergyCalibratorRun2(), and CalibratedElectronProducerRun2T< T >::produce().

 {
   assert(isMC_ == electron.isMC());
   float smear = 0.0, scale = 1.0;
   float aeta = std::abs(electron.getEta()); //, r9 = electron.getR9();
   float et = electron.getNewEnergy()/cosh(aeta);
   
   scale = _correctionRetriever.ScaleCorrection(electron.getRunNumber(), electron.isEB(), electron.getR9(), aeta, et);
   smear = _correctionRetriever.getSmearingSigma(electron.getRunNumber(), electron.isEB(), electron.getR9(), aeta, et, 0., 0.); 
   
   double newEcalEnergy, newEcalEnergyError;
   if (isMC_) {
     double corr = 1.0 + smear * gauss(id);
     newEcalEnergy      = electron.getNewEnergy() * corr;
     newEcalEnergyError = std::hypot(electron.getNewEnergyError() * corr, smear * newEcalEnergy);
   } else {
     newEcalEnergy      = electron.getNewEnergy() * scale;
     newEcalEnergyError = std::hypot(electron.getNewEnergyError() * scale, smear * newEcalEnergy);
   }
   electron.setNewEnergy(newEcalEnergy); 
   electron.setNewEnergyError(newEcalEnergyError);
   epCombinationTool_->combine(electron);
 }

void ElectronEnergyCalibratorRun2::calibrate	(	reco::GsfElectron &	electron,
		unsigned int	runNumber,
		edm::StreamID const &	id = `edm::StreamID::invalidStreamID()`
	)		const

Correct this electron. StreamID is needed when used with CMSSW Random Number Generator

Definition at line 34 of file ElectronEnergyCalibratorRun2.cc.

References calibrate(), isMC_, PileupJetIDParams_cfi::simple, and SimpleElectron::writeTo().

 {
   SimpleElectron simple(electron, runNumber, isMC_);
   calibrate(simple, id);
   simple.writeTo(electron);
 }

double ElectronEnergyCalibratorRun2::gauss ( edm::StreamID const & id ) const

protected

Return a number distributed as a unit gaussian, drawn from the private RNG if initPrivateRng was called, or from the CMSSW RandomNumberGenerator service If synchronization is set to true, it returns a fixed number (1.0)

Definition at line 64 of file ElectronEnergyCalibratorRun2.cc.

References Exception, edm::RandomNumberGenerator::getEngine(), edm::Service< T >::isAvailable(), rng_, and synchronization_.

Referenced by calibrate().

 {
   if (synchronization_) return 1.0;
   if (rng_) {
     return rng_->Gaus();
   } else {
     edm::Service<edm::RandomNumberGenerator> rng;
     if ( !rng.isAvailable() ) {
       throw cms::Exception("Configuration")
     << "XXXXXXX requires the RandomNumberGeneratorService\n"
     "which is not present in the configuration file.  You must add the service\n"
     "in the configuration file or remove the modules that require it.";
     }
     CLHEP::RandGaussQ gaussDistribution(rng->getEngine(id), 0.0, 1.0);
     return gaussDistribution.fire();
   }
 }