EGEtScaleSysModifier Class Reference

Inheritance diagram for EGEtScaleSysModifier:

Classes
class	UncertFuncBase
class	UncertFuncV1

Public Member Functions
	EGEtScaleSysModifier (const edm::ParameterSet &conf, edm::ConsumesCollector &)
void	modifyObject (pat::Electron &ele) const final
void	modifyObject (pat::Photon &pho) const final
void	setEvent (const edm::Event &) final
void	setEventContent (const edm::EventSetup &) final
	~EGEtScaleSysModifier () override
Public Member Functions inherited from ModifyObjectValueBase
virtual void	modifyObject (reco::GsfElectron &) const
virtual void	modifyObject (reco::Photon &) const
virtual void	modifyObject (reco::Muon &) const
virtual void	modifyObject (reco::BaseTau &) const
virtual void	modifyObject (reco::Jet &) const
virtual void	modifyObject (pat::Muon &) const
virtual void	modifyObject (pat::Tau &) const
virtual void	modifyObject (pat::Jet &) const
	ModifyObjectValueBase (const edm::ParameterSet &conf)
const std::string &	name () const
virtual	~ModifyObjectValueBase ()

Private Member Functions
std::pair< float, float >	calCombinedMom (reco::GsfElectron &ele, const float scale, const float smear) const
void	setEcalEnergy (reco::GsfElectron &ele, const float scale, const float smear) const

Private Attributes
EpCombinationTool	epCombTool_
std::unique_ptr< UncertFuncBase >	uncertFunc_

Detailed Description

Definition at line 28 of file EGEtScaleSysModifier.cc.

Constructor & Destructor Documentation

EGEtScaleSysModifier()

EGEtScaleSysModifier::EGEtScaleSysModifier	(	const edm::ParameterSet &	conf,
		edm::ConsumesCollector &	cc
	)

Definition at line 89 of file EGEtScaleSysModifier.cc.

References gpuPixelDoublets::cc, edm::ParameterSet::getParameterSet(), and eostools::move().

    : ModifyObjectValueBase(conf), epCombTool_{conf.getParameterSet("epCombConfig"), std::move(cc)} {
  const edm::ParameterSet& funcPSet = conf.getParameterSet("uncertFunc");
  const std::string& funcName = funcPSet.getParameter<std::string>("name");
  if (funcName == "UncertFuncV1") {
    uncertFunc_ = std::make_unique<UncertFuncV1>(funcPSet);
  } else {
    throw cms::Exception("ConfigError") << "Error constructing EGEtScaleSysModifier, function name " << funcName
                                        << " not valid";
  }
}

~EGEtScaleSysModifier()

EGEtScaleSysModifier::~EGEtScaleSysModifier ( )

inlineoverride

Definition at line 31 of file EGEtScaleSysModifier.cc.

{}

Member Function Documentation

calCombinedMom()

std::pair< float, float > EGEtScaleSysModifier::calCombinedMom ( reco::GsfElectron &  ele, const float  scale, const float  smear  ) const

private

Definition at line 165 of file EGEtScaleSysModifier.cc.

References EpCombinationTool::combine(), reco::GsfElectron::correctMomentum(), reco::GsfElectron::ecalEnergy(), reco::GsfElectron::ecalEnergyError(), epCombTool_, reco::GsfElectron::p4(), reco::GsfElectron::P4_COMBINATION, reco::GsfElectron::p4Error(), l1tEGammaCrystalsEmulatorProducer_cfi::scale, reco::GsfElectron::setCorrectedEcalEnergy(), reco::GsfElectron::setCorrectedEcalEnergyError(), setEcalEnergy(), and reco::GsfElectron::trackMomentumError().

Referenced by modifyObject().

                                                                                      {
  const float oldEcalEnergy = ele.ecalEnergy();
  const float oldEcalEnergyErr = ele.ecalEnergyError();
  const auto oldP4 = ele.p4();
  const float oldP4Err = ele.p4Error(reco::GsfElectron::P4_COMBINATION);
  const float oldTrkMomErr = ele.trackMomentumError();

  setEcalEnergy(ele, scale, smear);
  const auto& combinedMomentum = epCombTool_.combine(ele);
  ele.setCorrectedEcalEnergy(oldEcalEnergy);
  ele.setCorrectedEcalEnergyError(oldEcalEnergyErr);
  ele.correctMomentum(oldP4, oldTrkMomErr, oldP4Err);

  return combinedMomentum;
}

modifyObject() [1/2]

void EGEtScaleSysModifier::modifyObject ( pat::Electron &  ele ) const

finalvirtual

Reimplemented from ModifyObjectValueBase.

Definition at line 105 of file EGEtScaleSysModifier.cc.

References pat::PATObject< ObjectType >::addUserFloat(), calCombinedMom(), alignCSCRings::corr, reco::GsfElectron::correctMomentum(), reco::GsfElectron::ecalEnergy(), reco::GsfElectron::ecalEnergyError(), reco::LeafCandidate::energy(), reco::LeafCandidate::et(), EGEnergySysIndex::kEcalErrPreCorr, EGEnergySysIndex::kEcalPostCorr, EGEnergySysIndex::kEcalPreCorr, EGEnergySysIndex::kScaleValue, EGEnergySysIndex::kSmearValue, EGEnergySysIndex::name(), reco::GsfElectron::p4(), reco::GsfElectron::P4_COMBINATION, reco::GsfElectron::p4Error(), reco::GsfElectron::setCorrectedEcalEnergy(), reco::GsfElectron::setCorrectedEcalEnergyError(), reco::GsfElectron::trackMomentumError(), uncertFunc_, and pat::PATObject< ObjectType >::userFloat().

                                                              {
  auto getVal = [](const pat::Electron& ele, EGEnergySysIndex::Index valIndex) {
    return ele.userFloat(EGEnergySysIndex::name(valIndex));
  };
  //so ele.energy() may be either pre or post corrections, we have no idea
  //so we explicity access the pre and post correction ecal energies
  //we need the pre corrected to properly do the e/p combination
  //we need the post corrected to get et uncertainty
  const float ecalEnergyPostCorr = getVal(ele, EGEnergySysIndex::kEcalPostCorr);
  const float ecalEnergyPreCorr = getVal(ele, EGEnergySysIndex::kEcalPreCorr);
  const float ecalEnergyErrPreCorr = getVal(ele, EGEnergySysIndex::kEcalErrPreCorr);

  //the et cut is in terms of ecal et using the track angle and post corr ecal energy
  const float etUncert = uncertFunc_->val(ele.et() / ele.energy() * ecalEnergyPostCorr);
  const float smear = getVal(ele, EGEnergySysIndex::kSmearValue);
  const float corr = getVal(ele, EGEnergySysIndex::kScaleValue);

  //get the values we have to reset back to
  const float oldEcalEnergy = ele.ecalEnergy();
  const float oldEcalEnergyErr = ele.ecalEnergyError();
  const auto oldP4 = ele.p4();
  const float oldP4Err = ele.p4Error(reco::GsfElectron::P4_COMBINATION);
  const float oldTrkMomErr = ele.trackMomentumError();

  ele.setCorrectedEcalEnergy(ecalEnergyPreCorr);
  ele.setCorrectedEcalEnergyError(ecalEnergyErrPreCorr);

  const float energyEtUncertUp = calCombinedMom(ele, corr + etUncert, smear).first;
  const float energyEtUncertDn = calCombinedMom(ele, corr - etUncert, smear).first;

  //reset it back to how it was
  ele.setCorrectedEcalEnergy(oldEcalEnergy);
  ele.setCorrectedEcalEnergyError(oldEcalEnergyErr);
  ele.correctMomentum(oldP4, oldTrkMomErr, oldP4Err);

  ele.addUserFloat("energyScaleEtUp", energyEtUncertUp);
  ele.addUserFloat("energyScaleEtDown", energyEtUncertDn);
}

modifyObject() [2/2]

void EGEtScaleSysModifier::modifyObject ( pat::Photon &  pho ) const

finalvirtual

Reimplemented from ModifyObjectValueBase.

Definition at line 144 of file EGEtScaleSysModifier.cc.

References pat::PATObject< ObjectType >::addUserFloat(), alignCSCRings::corr, reco::LeafCandidate::energy(), reco::LeafCandidate::et(), EGEnergySysIndex::kEcalPostCorr, EGEnergySysIndex::kEcalPreCorr, EGEnergySysIndex::kScaleValue, EGEnergySysIndex::name(), uncertFunc_, and pat::PATObject< ObjectType >::userFloat().

                                                            {
  auto getVal = [](const pat::Photon& pho, EGEnergySysIndex::Index valIndex) {
    return pho.userFloat(EGEnergySysIndex::name(valIndex));
  };
  //so pho.energy() may be either pre or post corrections, we have no idea
  //so we explicity access the pre and post correction ecal energies
  //post corr for the et value for the systematic, pre corr to apply them
  const float ecalEnergyPostCorr = getVal(pho, EGEnergySysIndex::kEcalPostCorr);
  const float ecalEnergyPreCorr = getVal(pho, EGEnergySysIndex::kEcalPreCorr);

  //the et cut is in terms of post corr ecal energy
  const float etUncert = uncertFunc_->val(pho.et() / pho.energy() * ecalEnergyPostCorr);
  const float corr = getVal(pho, EGEnergySysIndex::kScaleValue);

  const float energyEtUncertUp = ecalEnergyPreCorr * (corr + etUncert);
  const float energyEtUncertDn = ecalEnergyPreCorr * (corr - etUncert);

  pho.addUserFloat("energyScaleEtUp", energyEtUncertUp);
  pho.addUserFloat("energyScaleEtDown", energyEtUncertDn);
}

setEcalEnergy()

void EGEtScaleSysModifier::setEcalEnergy ( reco::GsfElectron &  ele, const float  scale, const float  smear  ) const

private

Definition at line 183 of file EGEtScaleSysModifier.cc.

References reco::GsfElectron::ecalEnergy(), reco::GsfElectron::ecalEnergyError(), Matriplex::hypot(), l1tEGammaCrystalsEmulatorProducer_cfi::scale, reco::GsfElectron::setCorrectedEcalEnergy(), and reco::GsfElectron::setCorrectedEcalEnergyError().

Referenced by calCombinedMom().

                                                                                                         {
  const float oldEcalEnergy = ele.ecalEnergy();
  const float oldEcalEnergyErr = ele.ecalEnergyError();
  ele.setCorrectedEcalEnergy(oldEcalEnergy * scale);
  ele.setCorrectedEcalEnergyError(std::hypot(oldEcalEnergyErr * scale, oldEcalEnergy * smear * scale));
}

setEvent()

void EGEtScaleSysModifier::setEvent ( const edm::Event &  iEvent )

finalvirtual

Reimplemented from ModifyObjectValueBase.

Definition at line 101 of file EGEtScaleSysModifier.cc.

{}

setEventContent()

void EGEtScaleSysModifier::setEventContent ( const edm::EventSetup &  iSetup )

finalvirtual

Reimplemented from ModifyObjectValueBase.

Definition at line 103 of file EGEtScaleSysModifier.cc.

References epCombTool_, and EpCombinationTool::setEventContent().

{ epCombTool_.setEventContent(iSetup); }

Member Data Documentation

epCombTool_

EpCombinationTool EGEtScaleSysModifier::epCombTool_

private

Definition at line 85 of file EGEtScaleSysModifier.cc.

Referenced by calCombinedMom(), and setEventContent().

uncertFunc_

std::unique_ptr<UncertFuncBase> EGEtScaleSysModifier::uncertFunc_

private

Definition at line 86 of file EGEtScaleSysModifier.cc.

Referenced by modifyObject().