#include <PhotonEnergyCalibratorRun2.h>

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

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

void	initPrivateRng (TRandom *rnd)

	PhotonEnergyCalibratorRun2 ()

	PhotonEnergyCalibratorRun2 (bool isMC, bool synchronization, std::string correctionFile)

	~PhotonEnergyCalibratorRun2 ()

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

Protected Attributes
EnergyScaleCorrection_class	_correctionRetriever

bool	isMC_

TRandom *	rng_

std::vector< double >	scales_

std::vector< double >	smearings_

bool	synchronization_

Detailed Description

Definition at line 11 of file PhotonEnergyCalibratorRun2.h.

Constructor & Destructor Documentation

PhotonEnergyCalibratorRun2::PhotonEnergyCalibratorRun2 ( )

inline

Definition at line 14 of file PhotonEnergyCalibratorRun2.h.

14 {}

PhotonEnergyCalibratorRun2::PhotonEnergyCalibratorRun2	(	bool	isMC,
		bool	synchronization,
		std::string	correctionFile
	)

Definition at line 7 of file PhotonEnergyCalibratorRun2.cc.

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

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

PhotonEnergyCalibratorRun2::~PhotonEnergyCalibratorRun2 ( )

Definition at line 22 of file PhotonEnergyCalibratorRun2.cc.

22 {}

Member Function Documentation

void PhotonEnergyCalibratorRun2::calibrate	(	SimplePhoton &	photon,
		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 PhotonEnergyCalibratorRun2.cc.

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

Referenced by calibrate().

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

void PhotonEnergyCalibratorRun2::calibrate	(	reco::Photon &	photon,
		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 26 of file PhotonEnergyCalibratorRun2.cc.

References calibrate(), isMC_, and SimplePhoton::writeTo().

                                                                         {
   SimplePhoton simple(photon, runNumber, isMC_);
   calibrate(simple, id);
   simple.writeTo(photon);
 }

double PhotonEnergyCalibratorRun2::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 56 of file PhotonEnergyCalibratorRun2.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();
   }
 }