Inheritance diagram for PFCandidateRecalibrator:

Public Member Functions
	PFCandidateRecalibrator (const edm::ParameterSet &)

	~PFCandidateRecalibrator () override

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

Public Member Functions inherited from edm::stream::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

	~EDProducerBase () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

virtual	~ProducerBase () noexcept(false)

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

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

Static Public Member Functions inherited from edm::stream::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (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

Public Types inherited from edm::stream::EDProducerBase
typedef EDProducerAdaptorBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

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.

57 {};

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.

153 {

154 }

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().

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation