PhotonEnergyCalibrator Class Reference

#include <PhotonEnergyCalibrator.h>

Public Types
enum	EventType { EventType::DATA, EventType::MC }

Public Member Functions
std::array< float, EGEnergySysIndex::kNrSysErrs >	calibrate (reco::Photon &photon, const unsigned int runNumber, const EcalRecHitCollection *recHits, const float smearNrSigma, const EventType eventType) const
std::array< float, EGEnergySysIndex::kNrSysErrs >	calibrate (reco::Photon &photon, const unsigned int runNumber, const EcalRecHitCollection *recHits, edm::StreamID const &id, const EventType eventType) const
void	initPrivateRng (TRandom *rnd)
	PhotonEnergyCalibrator ()
	PhotonEnergyCalibrator (const std::string &correctionFile)
void	setMinEt (float val)
	~PhotonEnergyCalibrator ()

Private Member Functions
double	gauss (edm::StreamID const &id) const
void	setEnergyAndSystVarations (const float scale, const float smearNrSigma, const float et, const EnergyScaleCorrection::ScaleCorrection &scaleCorr, const EnergyScaleCorrection::SmearCorrection &smearCorr, reco::Photon &photon, std::array< float, EGEnergySysIndex::kNrSysErrs > &energyData) const

Private Attributes
EnergyScaleCorrection	correctionRetriever_
float	minEt_
TRandom *	rng_

Static Private Attributes
static const EnergyScaleCorrection::ScaleCorrection	defaultScaleCorr_
static const EnergyScaleCorrection::SmearCorrection	defaultSmearCorr_

Detailed Description

Definition at line 20 of file PhotonEnergyCalibrator.h.

Member Enumeration Documentation

EventType

enum PhotonEnergyCalibrator::EventType

strong

Enumerator
DATA
MC

Definition at line 22 of file PhotonEnergyCalibrator.h.

                       {
    DATA,
    MC,
  };

Constructor & Destructor Documentation

PhotonEnergyCalibrator() [1/2]

PhotonEnergyCalibrator::PhotonEnergyCalibrator ( )

inline

Definition at line 27 of file PhotonEnergyCalibrator.h.

{}

PhotonEnergyCalibrator() [2/2]

PhotonEnergyCalibrator::PhotonEnergyCalibrator

(

const std::string &

correctionFile

)

Definition at line 12 of file PhotonEnergyCalibrator.cc.

    : correctionRetriever_(correctionFile), rng_(nullptr), minEt_(1.0) {}

~PhotonEnergyCalibrator()

PhotonEnergyCalibrator::~PhotonEnergyCalibrator ( )

inline

Definition at line 29 of file PhotonEnergyCalibrator.h.

{}

Member Function Documentation

calibrate() [1/2]

std::array< float, EGEnergySysIndex::kNrSysErrs > PhotonEnergyCalibrator::calibrate	(	reco::Photon &	photon,
		const unsigned int	runNumber,
		const EcalRecHitCollection *	recHits,
		const float	smearNrSigma,
		const EventType	eventType
	)		const

Definition at line 26 of file PhotonEnergyCalibrator.cc.

                                                           {
  const float scEtaAbs = std::abs(photon.superCluster()->eta());
  const float et = photon.getCorrectedEnergy(reco::Photon::P4type::regression2) / cosh(scEtaAbs);
 
  if (et < minEt_ || edm::isNotFinite(et)) {
    std::array<float, EGEnergySysIndex::kNrSysErrs> retVal;
    retVal.fill(photon.getCorrectedEnergy(reco::Photon::P4type::regression2));
    retVal[EGEnergySysIndex::kScaleValue] = 1.0;
    retVal[EGEnergySysIndex::kSmearValue] = 0.0;
    retVal[EGEnergySysIndex::kSmearNrSigma] = smearNrSigma;
    retVal[EGEnergySysIndex::kEcalErrPreCorr] = photon.getCorrectedEnergyError(reco::Photon::P4type::regression2);
    retVal[EGEnergySysIndex::kEcalErrPostCorr] = photon.getCorrectedEnergyError(reco::Photon::P4type::regression2);
    retVal[EGEnergySysIndex::kEcalTrkPreCorr] = 0.;
    retVal[EGEnergySysIndex::kEcalTrkErrPreCorr] = 0.;
    retVal[EGEnergySysIndex::kEcalTrkPostCorr] = 0.;
    retVal[EGEnergySysIndex::kEcalTrkErrPostCorr] = 0.;
 
    return retVal;
  }
 
  const DetId seedDetId = photon.superCluster()->seed()->seed();
  EcalRecHitCollection::const_iterator seedRecHit = recHits->find(seedDetId);
  unsigned int gainSeedSC = 12;
  if (seedRecHit != recHits->end()) {
    if (seedRecHit->checkFlag(EcalRecHit::kHasSwitchToGain6))
      gainSeedSC = 6;
    if (seedRecHit->checkFlag(EcalRecHit::kHasSwitchToGain1))
      gainSeedSC = 1;
  }
 
  const EnergyScaleCorrection::ScaleCorrection* scaleCorr =
      correctionRetriever_.getScaleCorr(runNumber, et, scEtaAbs, photon.full5x5_r9(), gainSeedSC);
  const EnergyScaleCorrection::SmearCorrection* smearCorr =
      correctionRetriever_.getSmearCorr(runNumber, et, scEtaAbs, photon.full5x5_r9(), gainSeedSC);
  if (scaleCorr == nullptr)
    scaleCorr = &defaultScaleCorr_;
  if (smearCorr == nullptr)
    smearCorr = &defaultSmearCorr_;
 
  std::array<float, EGEnergySysIndex::kNrSysErrs> uncertainties{};
 
  uncertainties[EGEnergySysIndex::kScaleValue] = scaleCorr->scale();
  uncertainties[EGEnergySysIndex::kSmearValue] =
      smearCorr->sigma(et);  //even though we use scale = 1.0, we still store the value returned for MC
  uncertainties[EGEnergySysIndex::kSmearNrSigma] = smearNrSigma;
  //MC central values are not scaled (scale = 1.0), data is not smeared (smearNrSigma = 0)
  //smearing still has a second order effect on data as it enters the E/p combination as an
  //extra uncertainty on the calo energy
  //MC gets all the scale systematics
  if (eventType == EventType::DATA) {
    setEnergyAndSystVarations(scaleCorr->scale(), 0., et, *scaleCorr, *smearCorr, photon, uncertainties);
  } else if (eventType == EventType::MC) {
    setEnergyAndSystVarations(1.0, smearNrSigma, et, *scaleCorr, *smearCorr, photon, uncertainties);
  }
 
  return uncertainties;
}

References funct::abs(), correctionRetriever_, DATA, defaultScaleCorr_, defaultSmearCorr_, EgHLTOffHistBins_cfi::et, hcalcalib_dqm_sourceclient-live_cfg::eventType, EnergyScaleCorrection::getScaleCorr(), EnergyScaleCorrection::getSmearCorr(), edm::isNotFinite(), EGEnergySysIndex::kEcalErrPostCorr, EGEnergySysIndex::kEcalErrPreCorr, EGEnergySysIndex::kEcalTrkErrPostCorr, EGEnergySysIndex::kEcalTrkErrPreCorr, EGEnergySysIndex::kEcalTrkPostCorr, EGEnergySysIndex::kEcalTrkPreCorr, EcalRecHit::kHasSwitchToGain1, EcalRecHit::kHasSwitchToGain6, EGEnergySysIndex::kScaleValue, EGEnergySysIndex::kSmearNrSigma, EGEnergySysIndex::kSmearValue, MC, minEt_, muons2muons_cfi::photon, FastTrackerRecHitMaskProducer_cfi::recHits, convertSQLiteXML::runNumber, EnergyScaleCorrection::ScaleCorrection::scale(), setEnergyAndSystVarations(), and EnergyScaleCorrection::SmearCorrection::sigma().

calibrate() [2/2]

std::array< float, EGEnergySysIndex::kNrSysErrs > PhotonEnergyCalibrator::calibrate	(	reco::Photon &	photon,
		const unsigned int	runNumber,
		const EcalRecHitCollection *	recHits,
		edm::StreamID const &	id,
		const EventType	eventType
	)		const

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

Definition at line 17 of file PhotonEnergyCalibrator.cc.

                                                           {
  return calibrate(photon, runNumber, recHits, gauss(id), eventType);
}

References hcalcalib_dqm_sourceclient-live_cfg::eventType, gauss(), muons2muons_cfi::photon, FastTrackerRecHitMaskProducer_cfi::recHits, and convertSQLiteXML::runNumber.

gauss()

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

private

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 144 of file PhotonEnergyCalibrator.cc.

                                                                {
  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();
  }
}

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

Referenced by calibrate().

initPrivateRng()

void PhotonEnergyCalibrator::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 15 of file PhotonEnergyCalibrator.cc.

{ rng_ = rnd; }

References rng_.

Referenced by CalibratedPhotonProducerT< T >::CalibratedPhotonProducerT().

setEnergyAndSystVarations()

void PhotonEnergyCalibrator::setEnergyAndSystVarations ( const float  scale, const float  smearNrSigma, const float  et, const EnergyScaleCorrection::ScaleCorrection &  scaleCorr, const EnergyScaleCorrection::SmearCorrection &  smearCorr, reco::Photon &  photon, std::array< float, EGEnergySysIndex::kNrSysErrs > &  energyData  ) const

private

Definition at line 89 of file PhotonEnergyCalibrator.cc.

                                                                   {
  const float smear = smearCorr.sigma(et);
  const float smearRhoUp = smearCorr.sigma(et, 1, 0);
  const float smearRhoDn = smearCorr.sigma(et, -1, 0);
  const float smearPhiUp = smearCorr.sigma(et, 0, 1);
  const float smearPhiDn = smearCorr.sigma(et, 0, -1);
 
  const float corr = scale + smear * smearNrSigma;
  const float corrRhoUp = scale + smearRhoUp * smearNrSigma;
  const float corrRhoDn = scale + smearRhoDn * smearNrSigma;
  const float corrPhiUp = scale + smearPhiUp * smearNrSigma;
  const float corrPhiDn = scale + smearPhiDn * smearNrSigma;
  const float corrUp = corrRhoUp;
  const float corrDn = corrRhoDn;
 
  const double oldEcalEnergy = photon.getCorrectedEnergy(reco::Photon::P4type::regression2);
  const double oldEcalEnergyError = photon.getCorrectedEnergyError(reco::Photon::P4type::regression2);
 
  energyData[EGEnergySysIndex::kEcalPreCorr] = oldEcalEnergy;
  energyData[EGEnergySysIndex::kEcalErrPreCorr] = oldEcalEnergyError;
 
  const double newEcalEnergy = oldEcalEnergy * corr;
  const double newEcalEnergyError = std::hypot(oldEcalEnergyError * corr, smear * newEcalEnergy);
  photon.setCorrectedEnergy(reco::Photon::P4type::regression2, newEcalEnergy, newEcalEnergyError, true);
 
  energyData[EGEnergySysIndex::kScaleStatUp] = oldEcalEnergy * (corr + scaleCorr.scaleErrStat());
  energyData[EGEnergySysIndex::kScaleStatDown] = oldEcalEnergy * (corr - scaleCorr.scaleErrStat());
  energyData[EGEnergySysIndex::kScaleSystUp] = oldEcalEnergy * (corr + scaleCorr.scaleErrSyst());
  energyData[EGEnergySysIndex::kScaleSystDown] = oldEcalEnergy * (corr - scaleCorr.scaleErrSyst());
  energyData[EGEnergySysIndex::kScaleGainUp] = oldEcalEnergy * (corr + scaleCorr.scaleErrGain());
  energyData[EGEnergySysIndex::kScaleGainDown] = oldEcalEnergy * (corr - scaleCorr.scaleErrGain());
  energyData[EGEnergySysIndex::kSmearRhoUp] = oldEcalEnergy * corrRhoUp;
  energyData[EGEnergySysIndex::kSmearRhoDown] = oldEcalEnergy * corrRhoDn;
  energyData[EGEnergySysIndex::kSmearPhiUp] = oldEcalEnergy * corrPhiUp;
  energyData[EGEnergySysIndex::kSmearPhiDown] = oldEcalEnergy * corrPhiDn;
 
  // The total variation
  energyData[EGEnergySysIndex::kScaleUp] =
      oldEcalEnergy * (corr + scaleCorr.scaleErr(EnergyScaleCorrection::kErrStatSystGain));
  energyData[EGEnergySysIndex::kScaleDown] =
      oldEcalEnergy * (corr - scaleCorr.scaleErr(EnergyScaleCorrection::kErrStatSystGain));
  energyData[EGEnergySysIndex::kSmearUp] = oldEcalEnergy * corrUp;
  energyData[EGEnergySysIndex::kSmearDown] = oldEcalEnergy * corrDn;
 
  energyData[EGEnergySysIndex::kEcalPostCorr] = photon.getCorrectedEnergy(reco::Photon::P4type::regression2);
  energyData[EGEnergySysIndex::kEcalErrPostCorr] = photon.getCorrectedEnergyError(reco::Photon::P4type::regression2);
}

References alignCSCRings::corr, EgHLTOffHistBins_cfi::et, EGEnergySysIndex::kEcalErrPostCorr, EGEnergySysIndex::kEcalErrPreCorr, EGEnergySysIndex::kEcalPostCorr, EGEnergySysIndex::kEcalPreCorr, EnergyScaleCorrection::kErrStatSystGain, EGEnergySysIndex::kScaleDown, EGEnergySysIndex::kScaleGainDown, EGEnergySysIndex::kScaleGainUp, EGEnergySysIndex::kScaleStatDown, EGEnergySysIndex::kScaleStatUp, EGEnergySysIndex::kScaleSystDown, EGEnergySysIndex::kScaleSystUp, EGEnergySysIndex::kScaleUp, EGEnergySysIndex::kSmearDown, EGEnergySysIndex::kSmearPhiDown, EGEnergySysIndex::kSmearPhiUp, EGEnergySysIndex::kSmearRhoDown, EGEnergySysIndex::kSmearRhoUp, EGEnergySysIndex::kSmearUp, muons2muons_cfi::photon, HLT_FULL_cff::scale, EnergyScaleCorrection::ScaleCorrection::scaleErr(), EnergyScaleCorrection::ScaleCorrection::scaleErrGain(), EnergyScaleCorrection::ScaleCorrection::scaleErrStat(), EnergyScaleCorrection::ScaleCorrection::scaleErrSyst(), and EnergyScaleCorrection::SmearCorrection::sigma().

Referenced by calibrate().

setMinEt()

void PhotonEnergyCalibrator::setMinEt ( float  val )

inline

Definition at line 36 of file PhotonEnergyCalibrator.h.

{ minEt_ = val; }

References minEt_, and heppy_batch::val.

Referenced by CalibratedPhotonProducerT< T >::CalibratedPhotonProducerT().

Member Data Documentation

correctionRetriever_

EnergyScaleCorrection PhotonEnergyCalibrator::correctionRetriever_

private

Definition at line 67 of file PhotonEnergyCalibrator.h.

Referenced by calibrate().

defaultScaleCorr_

const EnergyScaleCorrection::ScaleCorrection PhotonEnergyCalibrator::defaultScaleCorr_

staticprivate

Definition at line 72 of file PhotonEnergyCalibrator.h.

Referenced by calibrate().

defaultSmearCorr_

const EnergyScaleCorrection::SmearCorrection PhotonEnergyCalibrator::defaultSmearCorr_

staticprivate

Definition at line 73 of file PhotonEnergyCalibrator.h.

Referenced by calibrate().

minEt_

float PhotonEnergyCalibrator::minEt_

private

Definition at line 69 of file PhotonEnergyCalibrator.h.

Referenced by calibrate(), and setMinEt().

rng_

TRandom* PhotonEnergyCalibrator::rng_

private

Definition at line 68 of file PhotonEnergyCalibrator.h.

Referenced by gauss(), and initPrivateRng().