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
bool	hasAbilityToProduceInLumis () const final
bool	hasAbilityToProduceInRuns () const final

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_
edm::ESWatcher< HcalRecNumberingRecord >	hcalDbWatcher_
edm::ESWatcher< HcalRespCorrsRcd >	hcalRCWatcher_
float	longFibreThr_
edm::EDGetTokenT< reco::PFCandidateCollection >	pfcandidates_
float	shortFibreThr_

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

Detailed Description

Definition at line 54 of file PFCandidateRecalibrator.cc.

Constructor & Destructor Documentation

PFCandidateRecalibrator::PFCandidateRecalibrator

(

const edm::ParameterSet &

iConfig

)

Definition at line 78 of file PFCandidateRecalibrator.cc.

                                                                               :
  pfcandidates_(consumes<reco::PFCandidateCollection>(iConfig.getParameter<edm::InputTag>("pfcandidates"))),
  shortFibreThr_(iConfig.getParameter<double>("shortFibreThr")),
  longFibreThr_(iConfig.getParameter<double>("longFibreThr"))
{
    produces<reco::PFCandidateCollection>();
    produces<reco::PFCandidateCollection>("discarded");
    produces<edm::ValueMap<reco::PFCandidateRef>>();
}

PFCandidateRecalibrator::~PFCandidateRecalibrator ( )

inlineoverride

Definition at line 57 of file PFCandidateRecalibrator.cc.

References Ecal2004TBTDCRanges_v1_cff::endRun, and fillDescriptions().

{};

Member Function Documentation

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

overrideprivate

Definition at line 88 of file PFCandidateRecalibrator.cc.

References funct::abs(), badChHE_, badChHF_, edm::ESWatcher< T >::check(), HcalDetId::depth(), PV3DBase< T, PVType, FrameType >::eta(), f, edm::EventSetup::get(), HcalDbService::getHcalRespCorr(), HcalGeometry::getPosition(), CaloGeometry::getSubdetectorGeometry(), HcalGeometry::getValidDetIds(), HcalRespCorr::getValue(), DetId::Hcal, hcalDbWatcher_, HcalEndcap, HcalForward, hcalRCWatcher_, HcalDetId::ieta(), HcalDetId::iphi(), PV3DBase< T, PVType, FrameType >::phi(), edm::ESHandle< T >::product(), particleFlowDisplacedVertex_cfi::ratio, and HcalCondObjectContainerBase::setTopo().

{
    if (hcalDbWatcher_.check(iSetup) || hcalRCWatcher_.check(iSetup))
      {
    //Get Calib Constants from current GT
    edm::ESHandle<HcalDbService> gtCond;
    iSetup.get<HcalDbRecord>().get(gtCond);
    
    //Get Calib Constants from bugged tag
    edm::ESHandle<HcalTopology> htopo;
    iSetup.get<HcalRecNumberingRecord>().get(htopo);
    const HcalTopology* theHBHETopology = htopo.product();
    
    edm::ESHandle<HcalRespCorrs> buggedCond;
    iSetup.get<HcalRespCorrsRcd>().get("bugged", buggedCond);
    HcalRespCorrs buggedRespCorrs(*buggedCond.product());
    buggedRespCorrs.setTopo(theHBHETopology);
    
    //access calogeometry
    edm::ESHandle<CaloGeometry> calogeom;
    iSetup.get<CaloGeometryRecord>().get(calogeom);
    const CaloGeometry* cgeo = calogeom.product();
    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));
          }
      }
      }
}

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

overrideprivate

Definition at line 152 of file PFCandidateRecalibrator.cc.

{
}

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

static

Definition at line 341 of file PFCandidateRecalibrator.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), and DEFINE_FWK_MODULE.

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

}

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

overrideprivate

Definition at line 156 of file PFCandidateRecalibrator.cc.

References funct::abs(), badChHE_, badChHF_, popcon2dropbox::copy(), reco::deltaR2(), edm::helper::Filler< Map >::fill(), objects.autophobj::filler, edm::EventSetup::get(), edm::Event::getByToken(), HcalGeometry::getClosestCell(), CaloGeometry::getSubdetectorGeometry(), DetId::Hcal, HcalForward, mps_fire::i, edm::helper::Filler< Map >::insert(), LogDebug, longFibreThr_, eostools::move(), lumidbDDL::newToOld(), lumidbDDL::oldToNew(), packedPFCandidateRefMixer_cfi::pf, packedPFCandidateRefMixer_cfi::pf2pf, pfcandidates_, edm::ESHandle< T >::product(), edm::Event::put(), and shortFibreThr_.

Referenced by JSONExport.JsonExport::export(), HTMLExport.HTMLExport::export(), and HTMLExport.HTMLExportStatic::export().

{
    //access calogeometry
    edm::ESHandle<CaloGeometry> calogeom;
    iSetup.get<CaloGeometryRecord>().get(calogeom);
    const CaloGeometry* cgeo = calogeom.product();
    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

std::vector<HEChannel> PFCandidateRecalibrator::badChHE_

private

Definition at line 71 of file PFCandidateRecalibrator.cc.

Referenced by beginRun(), and produce().

std::vector<HFChannel> PFCandidateRecalibrator::badChHF_

private

Definition at line 72 of file PFCandidateRecalibrator.cc.

Referenced by beginRun(), and produce().

edm::ESWatcher<HcalRecNumberingRecord> PFCandidateRecalibrator::hcalDbWatcher_

private

Definition at line 66 of file PFCandidateRecalibrator.cc.

Referenced by beginRun().

edm::ESWatcher<HcalRespCorrsRcd> PFCandidateRecalibrator::hcalRCWatcher_

private

Definition at line 67 of file PFCandidateRecalibrator.cc.

Referenced by beginRun().

float PFCandidateRecalibrator::longFibreThr_

private

Definition at line 75 of file PFCandidateRecalibrator.cc.

Referenced by produce().

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

private

Definition at line 69 of file PFCandidateRecalibrator.cc.

Referenced by produce().

float PFCandidateRecalibrator::shortFibreThr_

private

Definition at line 74 of file PFCandidateRecalibrator.cc.

Referenced by produce().