PFCandidateRecalibrator Class Reference

Inheritance diagram for PFCandidateRecalibrator:

Public Member Functions
	PFCandidateRecalibrator (const edm::ParameterSet &)
	~PFCandidateRecalibrator () override
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Private Member Functions
void	beginRun (const edm::Run &iRun, edm::EventSetup const &iSetup) override
void	endRun (const edm::Run &iRun, edm::EventSetup const &iSetup) override
void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
std::vector< HEChannel >	badChHE_
std::vector< HFChannel >	badChHF_
const edm::ESGetToken< HcalRespCorrs, HcalRespCorrsRcd >	buggedCondToken_
const edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	calogeomTokenEvent_
const edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	calogeomTokenRun_
const edm::ESGetToken< HcalDbService, HcalDbRecord >	gtCondToken_
edm::ESWatcher< HcalRecNumberingRecord >	hcalDbWatcher_
edm::ESWatcher< HcalRespCorrsRcd >	hcalRCWatcher_
const edm::ESGetToken< HcalTopology, HcalRecNumberingRecord >	htopoToken_
float	longFibreThr_
edm::EDGetTokenT< reco::PFCandidateCollection >	pfcandidates_
float	shortFibreThr_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
using	CacheTypes = CacheContexts< T... >
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T... >
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

Definition at line 50 of file PFCandidateRecalibrator.cc.

Constructor & Destructor Documentation

PFCandidateRecalibrator()

PFCandidateRecalibrator::PFCandidateRecalibrator

(

const edm::ParameterSet &

iConfig

)

Definition at line 80 of file PFCandidateRecalibrator.cc.

    : pfcandidates_(consumes<reco::PFCandidateCollection>(iConfig.getParameter<edm::InputTag>("pfcandidates"))),
      gtCondToken_(esConsumes<edm::Transition::BeginRun>()),
      htopoToken_(esConsumes<edm::Transition::BeginRun>()),
      buggedCondToken_(esConsumes<edm::Transition::BeginRun>(edm::ESInputTag("", "bugged"))),
      calogeomTokenRun_(esConsumes<edm::Transition::BeginRun>()),
      calogeomTokenEvent_(esConsumes()),
      shortFibreThr_(iConfig.getParameter<double>("shortFibreThr")),
      longFibreThr_(iConfig.getParameter<double>("longFibreThr")) {
  produces<reco::PFCandidateCollection>();
  produces<reco::PFCandidateCollection>("discarded");
  produces<edm::ValueMap<reco::PFCandidateRef>>();
}

~PFCandidateRecalibrator()

PFCandidateRecalibrator::~PFCandidateRecalibrator ( )

inlineoverride

Definition at line 53 of file PFCandidateRecalibrator.cc.

{};

Member Function Documentation

beginRun()

void PFCandidateRecalibrator::beginRun ( const edm::Run &  iRun, edm::EventSetup const &  iSetup  )

overrideprivate

Definition at line 94 of file PFCandidateRecalibrator.cc.

References funct::abs(), badChHE_, badChHF_, buggedCondToken_, calogeomTokenRun_, edm::ESWatcher< T >::check(), HcalDetId::depth(), f, edm::EventSetup::getData(), HcalDbService::getHcalRespCorr(), HcalGeometry::getPosition(), CaloGeometry::getSubdetectorGeometry(), HcalGeometry::getValidDetIds(), HcalRespCorr::getValue(), HcalCondObjectContainer< Item >::getValues(), gtCondToken_, DetId::Hcal, hcalDbWatcher_, HcalEndcap, HcalForward, hcalRCWatcher_, htopoToken_, HcalDetId::ieta(), HcalDetId::iphi(), particleFlowDisplacedVertex_cfi::ratio, and HcalCondObjectContainerBase::setTopo().

                                                                                      {
  auto newDb = hcalDbWatcher_.check(iSetup);
  auto newRC = hcalRCWatcher_.check(iSetup);
  if (newDb || newRC) {
    //Get Calib Constants from current GT
    HcalDbService const& gtCond = iSetup.getData(gtCondToken_);

    //Get Calib Constants from bugged tag
    const HcalTopology& theHBHETopology = iSetup.getData(htopoToken_);

    HcalRespCorrs buggedRespCorrs = iSetup.getData(buggedCondToken_);
    buggedRespCorrs.setTopo(&theHBHETopology);

    //access calogeometry
    const CaloGeometry& cgeo = iSetup.getData(calogeomTokenRun_);
    const HcalGeometry* hgeom = static_cast<const HcalGeometry*>(cgeo.getSubdetectorGeometry(DetId::Hcal, HcalForward));

    //reset the bad channel containers
    badChHE_.clear();
    badChHF_.clear();

    //fill bad cells HE (use eta, phi)
    const std::vector<DetId>& cellsHE = hgeom->getValidDetIds(DetId::Detector::Hcal, HcalEndcap);
    for (auto id : cellsHE) {
      float currentRespCorr = gtCond.getHcalRespCorr(id)->getValue();
      float buggedRespCorr = buggedRespCorrs.getValues(id)->getValue();
      if (buggedRespCorr == 0.)
        continue;

      float ratio = currentRespCorr / buggedRespCorr;
      if (std::abs(ratio - 1.f) > 0.001) {
        GlobalPoint pos = hgeom->getPosition(id);
        badChHE_.push_back(HEChannel(pos.eta(), pos.phi(), ratio));
      }
    }

    //fill bad cells HF (use ieta, iphi)
    auto const& cellsHF = hgeom->getValidDetIds(DetId::Detector::Hcal, HcalForward);
    for (auto id : cellsHF) {
      float currentRespCorr = gtCond.getHcalRespCorr(id)->getValue();
      float buggedRespCorr = buggedRespCorrs.getValues(id)->getValue();
      if (buggedRespCorr == 0.)
        continue;

      float ratio = currentRespCorr / buggedRespCorr;
      if (std::abs(ratio - 1.f) > 0.001) {
        HcalDetId dummyId(id);
        badChHF_.push_back(HFChannel(dummyId.ieta(), dummyId.iphi(), dummyId.depth(), ratio));
      }
    }
  }
}

endRun()

void PFCandidateRecalibrator::endRun ( const edm::Run &  iRun, edm::EventSetup const &  iSetup  )

overrideprivate

Definition at line 147 of file PFCandidateRecalibrator.cc.

{}

fillDescriptions()

void PFCandidateRecalibrator::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

static

Definition at line 312 of file PFCandidateRecalibrator.cc.

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, and ProducerED_cfi::InputTag.

                                                                                         {
  edm::ParameterSetDescription desc;

  desc.add<edm::InputTag>("pfcandidates", edm::InputTag("particleFlow"));
  desc.add<double>("shortFibreThr", 1.4);
  desc.add<double>("longFibreThr", 1.4);

  descriptions.add("pfCandidateRecalibrator", desc);
}

produce()

void PFCandidateRecalibrator::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

overrideprivate

Definition at line 149 of file PFCandidateRecalibrator.cc.

References funct::abs(), badChHE_, badChHF_, calogeomTokenEvent_, filterCSVwithJSON::copy, reco::deltaR2(), trigObjTnPSource_cfi::filler, HcalGeometry::getClosestCell(), edm::EventSetup::getData(), CaloGeometry::getSubdetectorGeometry(), DetId::Hcal, HcalForward, mps_fire::i, iEvent, LogDebug, longFibreThr_, eostools::move(), packedPFCandidateRefMixer_cfi::pf, packedPFCandidateRefMixer_cfi::pf2pf, OniaPhotonConversionProducer_cfi::pfcandidates, pfcandidates_, and shortFibreThr_.

                                                                                   {
  //access calogeometry
  const CaloGeometry& cgeo = iSetup.getData(calogeomTokenEvent_);
  const HcalGeometry* hgeom = static_cast<const HcalGeometry*>(cgeo.getSubdetectorGeometry(DetId::Hcal, HcalForward));

  //access PFCandidates
  edm::Handle<reco::PFCandidateCollection> pfcandidates;
  iEvent.getByToken(pfcandidates_, pfcandidates);

  int nPfCand = pfcandidates->size();
  std::unique_ptr<reco::PFCandidateCollection> copy(new reco::PFCandidateCollection());
  std::unique_ptr<reco::PFCandidateCollection> discarded(new reco::PFCandidateCollection());
  copy->reserve(nPfCand);
  std::vector<int> oldToNew(nPfCand), newToOld, badToOld;
  newToOld.reserve(nPfCand);

  LogDebug("PFCandidateRecalibrator") << "NEW EV:";

  //loop over PFCandidates
  int i = -1;
  for (const reco::PFCandidate& pf : *pfcandidates) {
    ++i;
    float absEta = std::abs(pf.eta());

    //deal with HE
    if (pf.particleId() == reco::PFCandidate::ParticleType::h0 &&
        !badChHE_.empty() &&  //don't touch if no miscalibration is found
        absEta > 1.4 && absEta < 3.) {
      bool toKill = false;
      for (auto const& badIt : badChHE_)
        if (reco::deltaR2(pf.eta(), pf.phi(), badIt.eta, badIt.phi) < 0.07)
          toKill = true;

      if (toKill) {
        discarded->push_back(pf);
        oldToNew[i] = (-discarded->size());
        badToOld.push_back(i);
        continue;
      } else {
        copy->push_back(pf);
        oldToNew[i] = (copy->size());
        newToOld.push_back(i);
      }
    }
    //deal with HF
    else if ((pf.particleId() == reco::PFCandidate::ParticleType::h_HF ||
              pf.particleId() == reco::PFCandidate::ParticleType::egamma_HF) &&
             !badChHF_.empty() &&  //don't touch if no miscalibration is found
             absEta >= 3.) {
      const math::XYZPointF& ecalPoint = pf.positionAtECALEntrance();
      GlobalPoint ecalGPoint(ecalPoint.X(), ecalPoint.Y(), ecalPoint.Z());
      HcalDetId closestDetId(hgeom->getClosestCell(ecalGPoint));

      if (closestDetId.subdet() == HcalForward) {
        HcalDetId hDetId(closestDetId.subdet(), closestDetId.ieta(), closestDetId.iphi(), 1);  //depth1

        //raw*calEnergy() is the same as *calEnergy() - no corrections are done for HF
        float longE = pf.rawEcalEnergy() + pf.rawHcalEnergy() / 2.;  //depth1
        float shortE = pf.rawHcalEnergy() / 2.;                      //depth2

        float ecalEnergy = pf.rawEcalEnergy();
        float hcalEnergy = pf.rawHcalEnergy();
        float totEnergy = ecalEnergy + hcalEnergy;
        float totEnergyOrig = totEnergy;

        bool toKill = false;

        for (auto const& badIt : badChHF_) {
          if (hDetId.ieta() == badIt.ieta && hDetId.iphi() == badIt.iphi) {
            LogDebug("PFCandidateRecalibrator")
                << "==> orig en (tot,H,E): " << pf.energy() << " " << pf.rawHcalEnergy() << " " << pf.rawEcalEnergy();
            if (badIt.depth == 1)  //depth1
            {
              longE *= badIt.ratio;
              ecalEnergy = ((longE - shortE) > 0.) ? (longE - shortE) : 0.;
              totEnergy = ecalEnergy + hcalEnergy;
            } else  //depth2
            {
              shortE *= badIt.ratio;
              hcalEnergy = 2 * shortE;
              ecalEnergy = ((longE - shortE) > 0.) ? (longE - shortE) : 0.;
              totEnergy = ecalEnergy + hcalEnergy;
            }
            //kill candidate if goes below thr
            if ((pf.particleId() == reco::PFCandidate::ParticleType::h_HF && shortE < shortFibreThr_) ||
                (pf.particleId() == reco::PFCandidate::ParticleType::egamma_HF && longE < longFibreThr_))
              toKill = true;

            LogDebug("PFCandidateRecalibrator") << "====> ieta,iphi,depth: " << badIt.ieta << " " << badIt.iphi << " "
                                                << badIt.depth << " corr: " << badIt.ratio;
            LogDebug("PFCandidateRecalibrator")
                << "====> recal en (tot,H,E): " << totEnergy << " " << hcalEnergy << " " << ecalEnergy;
          }
        }

        if (toKill) {
          discarded->push_back(pf);
          oldToNew[i] = (-discarded->size());
          badToOld.push_back(i);

          LogDebug("PFCandidateRecalibrator") << "==> KILLED ";
        } else {
          copy->push_back(pf);
          oldToNew[i] = (copy->size());
          newToOld.push_back(i);

          copy->back().setHcalEnergy(hcalEnergy, hcalEnergy);
          copy->back().setEcalEnergy(ecalEnergy, ecalEnergy);

          float scalingFactor = totEnergy / totEnergyOrig;
          math::XYZTLorentzVector recalibP4 = pf.p4() * scalingFactor;
          copy->back().setP4(recalibP4);

          LogDebug("PFCandidateRecalibrator") << "====> stored en (tot,H,E): " << copy->back().energy() << " "
                                              << copy->back().hcalEnergy() << " " << copy->back().ecalEnergy();
        }
      } else {
        copy->push_back(pf);
        oldToNew[i] = (copy->size());
        newToOld.push_back(i);
      }
    } else {
      copy->push_back(pf);
      oldToNew[i] = (copy->size());
      newToOld.push_back(i);
    }
  }

  // Now we put things in the event
  edm::OrphanHandle<reco::PFCandidateCollection> newpf = iEvent.put(std::move(copy));
  edm::OrphanHandle<reco::PFCandidateCollection> badpf = iEvent.put(std::move(discarded), "discarded");

  std::unique_ptr<edm::ValueMap<reco::PFCandidateRef>> pf2pf(new edm::ValueMap<reco::PFCandidateRef>());
  edm::ValueMap<reco::PFCandidateRef>::Filler filler(*pf2pf);
  std::vector<reco::PFCandidateRef> refs;
  refs.reserve(nPfCand);

  // old to new
  for (auto iOldToNew : oldToNew) {
    if (iOldToNew > 0) {
      refs.push_back(reco::PFCandidateRef(newpf, iOldToNew - 1));
    } else {
      refs.push_back(reco::PFCandidateRef(badpf, -iOldToNew - 1));
    }
  }
  filler.insert(pfcandidates, refs.begin(), refs.end());
  // new good to old
  refs.clear();
  for (int i : newToOld) {
    refs.push_back(reco::PFCandidateRef(pfcandidates, i));
  }
  filler.insert(newpf, refs.begin(), refs.end());
  // new bad to old
  refs.clear();
  for (int i : badToOld) {
    refs.push_back(reco::PFCandidateRef(pfcandidates, i));
  }
  filler.insert(badpf, refs.begin(), refs.end());
  // done
  filler.fill();
  iEvent.put(std::move(pf2pf));
}

Member Data Documentation

badChHE_

std::vector<HEChannel> PFCandidateRecalibrator::badChHE_

private

Definition at line 73 of file PFCandidateRecalibrator.cc.

Referenced by beginRun(), and produce().

badChHF_

std::vector<HFChannel> PFCandidateRecalibrator::badChHF_

private

Definition at line 74 of file PFCandidateRecalibrator.cc.

Referenced by beginRun(), and produce().

buggedCondToken_

const edm::ESGetToken<HcalRespCorrs, HcalRespCorrsRcd> PFCandidateRecalibrator::buggedCondToken_

private

Definition at line 69 of file PFCandidateRecalibrator.cc.

Referenced by beginRun().

calogeomTokenEvent_

const edm::ESGetToken<CaloGeometry, CaloGeometryRecord> PFCandidateRecalibrator::calogeomTokenEvent_

private

Definition at line 71 of file PFCandidateRecalibrator.cc.

Referenced by produce().

calogeomTokenRun_

const edm::ESGetToken<CaloGeometry, CaloGeometryRecord> PFCandidateRecalibrator::calogeomTokenRun_

private

Definition at line 70 of file PFCandidateRecalibrator.cc.

Referenced by beginRun().

gtCondToken_

const edm::ESGetToken<HcalDbService, HcalDbRecord> PFCandidateRecalibrator::gtCondToken_

private

Definition at line 67 of file PFCandidateRecalibrator.cc.

Referenced by beginRun().

hcalDbWatcher_

edm::ESWatcher<HcalRecNumberingRecord> PFCandidateRecalibrator::hcalDbWatcher_

private

Definition at line 62 of file PFCandidateRecalibrator.cc.

Referenced by beginRun().

hcalRCWatcher_

edm::ESWatcher<HcalRespCorrsRcd> PFCandidateRecalibrator::hcalRCWatcher_

private

Definition at line 63 of file PFCandidateRecalibrator.cc.

Referenced by beginRun().

htopoToken_

const edm::ESGetToken<HcalTopology, HcalRecNumberingRecord> PFCandidateRecalibrator::htopoToken_

private

Definition at line 68 of file PFCandidateRecalibrator.cc.

Referenced by beginRun().

longFibreThr_

float PFCandidateRecalibrator::longFibreThr_

private

Definition at line 77 of file PFCandidateRecalibrator.cc.

Referenced by produce().

pfcandidates_

edm::EDGetTokenT<reco::PFCandidateCollection> PFCandidateRecalibrator::pfcandidates_

private

Definition at line 65 of file PFCandidateRecalibrator.cc.

Referenced by produce().

shortFibreThr_

float PFCandidateRecalibrator::shortFibreThr_

private

Definition at line 76 of file PFCandidateRecalibrator.cc.

Referenced by produce().