#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 12 of file ElectronEnergyCalibratorRun2.h.

Constructor & Destructor Documentation

ElectronEnergyCalibratorRun2::ElectronEnergyCalibratorRun2 ( )

inline

Definition at line 15 of file ElectronEnergyCalibratorRun2.h.

15 {}

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_(nullptr),
       _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.

26 {}

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 37 of file ElectronEnergyCalibratorRun2.cc.

References _correctionRetriever, funct::abs(), cms::cuda::assert(), EpCombinationTool::combine(), alignCSCRings::corr, epCombinationTool_, gauss(), SimpleElectron::getEta(), SimpleElectron::getNewEnergy(), SimpleElectron::getNewEnergyError(), SimpleElectron::getR9(), SimpleElectron::getRunNumber(), EnergyScaleCorrection_class::getSmearingSigma(), SimpleElectron::isEB(), SimpleElectron::isMC(), isMC_, pileupReCalc_HLTpaths::scale, EnergyScaleCorrection_class::ScaleCorrection(), SimpleElectron::setNewEnergy(), and SimpleElectron::setNewEnergyError().

Referenced by calibrate().

                                                                                                   {
   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 30 of file ElectronEnergyCalibratorRun2.cc.

References calibrate(), isMC_, 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 61 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();
   }
 }