PhotonEnergyCalibratorRun2 Class Reference

#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.

{}

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_(0),
  _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 23 of file PhotonEnergyCalibratorRun2.cc.

{}

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 36 of file PhotonEnergyCalibratorRun2.cc.

References _correctionRetriever, funct::abs(), assert(), 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 30 of file PhotonEnergyCalibratorRun2.cc.

References calibrate(), isMC_, PileupJetIDParams_cfi::simple, 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 59 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();
  }
}

void PhotonEnergyCalibratorRun2::initPrivateRng

(

TRandom *

rnd

)

Initialize with a random number generator (if not done, it will use the CMSSW service) Caller code owns the TRandom.

Definition at line 26 of file PhotonEnergyCalibratorRun2.cc.

References rng_.

                                                            {
     rng_ = rnd;   
}

Member Data Documentation

EnergyScaleCorrection_class PhotonEnergyCalibratorRun2::_correctionRetriever

protected

Definition at line 44 of file PhotonEnergyCalibratorRun2.h.

Referenced by calibrate(), and PhotonEnergyCalibratorRun2().

bool PhotonEnergyCalibratorRun2::isMC_

protected

Definition at line 35 of file PhotonEnergyCalibratorRun2.h.

Referenced by calibrate(), and PhotonEnergyCalibratorRun2().

TRandom* PhotonEnergyCalibratorRun2::rng_

protected

Definition at line 36 of file PhotonEnergyCalibratorRun2.h.

Referenced by gauss(), and initPrivateRng().

std::vector<double> PhotonEnergyCalibratorRun2::scales_

protected

Definition at line 38 of file PhotonEnergyCalibratorRun2.h.

std::vector<double> PhotonEnergyCalibratorRun2::smearings_

protected

Definition at line 37 of file PhotonEnergyCalibratorRun2.h.

bool PhotonEnergyCalibratorRun2::synchronization_

protected

Definition at line 35 of file PhotonEnergyCalibratorRun2.h.

Referenced by gauss().