IsolatedEcalPixelTrackCandidateProducer Class Reference

#include <Calibration/HcalIsolatedTrackReco/src/IsolatedEcalPixelTrackCandidateProducer.cc>

Inheritance diagram for IsolatedEcalPixelTrackCandidateProducer:

Public Member Functions
	IsolatedEcalPixelTrackCandidateProducer (const edm::ParameterSet &)
	~IsolatedEcalPixelTrackCandidateProducer ()
Public Member Functions inherited from edm::global::EDProducer<>
	EDProducer ()=default
Public Member Functions inherited from edm::global::EDProducerBase
	EDProducerBase ()
ModuleDescription const &	moduleDescription () const
virtual	~EDProducerBase ()
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 ()
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, std::unordered_multimap< std::string, edm::ProductResolverIndex > 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)

Private Member Functions
virtual void	produce (edm::StreamID, edm::Event &, const edm::EventSetup &) const override

Private Attributes
const double	coneSizeEta0_
const double	coneSizeEta1_
const double	hitCountEthr_
const double	hitEthr_
const edm::EDGetTokenT< EcalRecHitCollection >	tok_eb
const edm::EDGetTokenT< EcalRecHitCollection >	tok_ee
const edm::EDGetTokenT< trigger::TriggerFilterObjectWithRefs >	tok_trigcand

Additional Inherited Members
Public Types inherited from edm::global::EDProducerBase
typedef EDProducerBase	ModuleType
Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::global::EDProducerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
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

Description: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 20 of file IsolatedEcalPixelTrackCandidateProducer.h.

Constructor & Destructor Documentation

IsolatedEcalPixelTrackCandidateProducer::IsolatedEcalPixelTrackCandidateProducer ( const edm::ParameterSet &  conf )

explicit

Definition at line 37 of file IsolatedEcalPixelTrackCandidateProducer.cc.

                                                                                                            :
  tok_ee(        consumes<EcalRecHitCollection>(conf.getParameter<edm::InputTag>("EERecHitSource")) ),
  tok_eb(        consumes<EcalRecHitCollection>(conf.getParameter<edm::InputTag>("EBRecHitSource")) ),
  tok_trigcand(  consumes<trigger::TriggerFilterObjectWithRefs>(conf.getParameter<edm::InputTag>("filterLabel")) ),
  coneSizeEta0_( conf.getParameter<double>("EcalConeSizeEta0") ),
  coneSizeEta1_( conf.getParameter<double>("EcalConeSizeEta1") ),
  hitCountEthr_( conf.getParameter<double>("ECHitCountEnergyThreshold") ),
  hitEthr_(      conf.getParameter<double>("ECHitEnergyThreshold") )
{
  // register the products
  produces< reco::IsolatedPixelTrackCandidateCollection >();
}

IsolatedEcalPixelTrackCandidateProducer::~IsolatedEcalPixelTrackCandidateProducer

(

)

Definition at line 50 of file IsolatedEcalPixelTrackCandidateProducer.cc.

{ }

Member Function Documentation

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

overrideprivatevirtual

Definition at line 53 of file IsolatedEcalPixelTrackCandidateProducer.cc.

References funct::abs(), edm::SortedCollection< T, SORT >::begin(), coneSizeEta0_, coneSizeEta1_, reco::deltaR(), edm::SortedCollection< T, SORT >::end(), PVValHelper::eta, PV3DBase< T, PVType, FrameType >::eta(), edm::EventSetup::get(), edm::Event::getByToken(), trigger::TriggerRefsCollections::getObjects(), CaloGeometry::getPosition(), hitCountEthr_, hitEthr_, edm::HandleBase::isValid(), eostools::move(), AlCaHLTBitMon_ParallelJobs::p, phi, PV3DBase< T, PVType, FrameType >::phi(), edm::ESHandle< T >::product(), edm::Event::put(), dttmaxenums::R, reco::IsolatedPixelTrackCandidate::setEnergyIn(), reco::IsolatedPixelTrackCandidate::setEnergyOut(), reco::IsolatedPixelTrackCandidate::setNHitIn(), reco::IsolatedPixelTrackCandidate::setNHitOut(), edm::SortedCollection< T, SORT >::size(), tok_eb, tok_ee, tok_trigcand, HiIsolationCommonParameters_cff::track, and trigger::TriggerTrack.

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

                                                                                                                      {
#ifdef DebugLog
  edm::LogInfo("HcalIsoTrack") << "==============Inside IsolatedEcalPixelTrackCandidateProducer";
#endif
  edm::ESHandle<CaloGeometry> pG;
  iSetup.get<CaloGeometryRecord>().get(pG);
  const CaloGeometry* geo = pG.product();

  edm::Handle<EcalRecHitCollection> ecalEB;
  iEvent.getByToken(tok_eb,ecalEB);

  edm::Handle<EcalRecHitCollection> ecalEE;
  iEvent.getByToken(tok_ee,ecalEE);
#ifdef DebugLog
  edm::LogInfo("HcalIsoTrack") << "ecal Collections isValid: " << ecalEB.isValid() << "/" << ecalEE.isValid();
#endif

  edm::Handle<trigger::TriggerFilterObjectWithRefs> trigCand;
  iEvent.getByToken(tok_trigcand,trigCand);

  std::vector< edm::Ref<reco::IsolatedPixelTrackCandidateCollection> > isoPixTrackRefs;
  trigCand->getObjects(trigger::TriggerTrack, isoPixTrackRefs);
  int nCand=isoPixTrackRefs.size();  

  auto iptcCollection = std::make_unique<reco::IsolatedPixelTrackCandidateCollection>();
#ifdef DebugLog
  edm::LogInfo("HcalIsoTrack") << "coneSize_ " << coneSizeEta0_ << "/"<< coneSizeEta1_ << " hitCountEthr_ " << hitCountEthr_ << " hitEthr_ " << hitEthr_;
#endif
  for (int p=0; p<nCand; p++) {
    int    nhitIn(0), nhitOut(0);
    double inEnergy(0), outEnergy(0);
    std::pair<double,double> etaPhi(isoPixTrackRefs[p]->track()->eta(), isoPixTrackRefs[p]->track()->phi());
    if (isoPixTrackRefs[p]->etaPhiEcalValid()) etaPhi = isoPixTrackRefs[p]->etaPhiEcal();
    double etaAbs = std::abs(etaPhi.first);
    double coneSize_ = (etaAbs > 1.5) ? coneSizeEta1_ : (coneSizeEta0_*(1.5-etaAbs)+coneSizeEta1_*etaAbs)/1.5;
#ifdef DebugLog
    edm::LogInfo("HcalIsoTrack") << "Track: eta/phi " << etaPhi.first << "/" << etaPhi.second << " pt:" << isoPixTrackRefs[p]->track()->pt() << " cone " << coneSize_ << "\n" << "rechit size EB/EE : " << ecalEB->size() << "/" << ecalEE->size() << " coneSize_: " << coneSize_;
#endif
    if (etaAbs<1.7) {
      for (EcalRecHitCollection::const_iterator eItr=ecalEB->begin(); eItr!=ecalEB->end(); eItr++) {
    GlobalPoint pos = geo->getPosition(eItr->detid());
    double      R   = reco::deltaR(pos.eta(),pos.phi(),etaPhi.first,etaPhi.second);
    if (R < coneSize_) {
      nhitIn++;
      inEnergy += (eItr->energy());
      if (eItr->energy() > hitCountEthr_) nhitOut++;
      if (eItr->energy() > hitEthr_)      outEnergy += (eItr->energy());
#ifdef DebugLog
      edm::LogInfo("HcalIsoTrack") << "Rechit Close to the track has energy " << eItr->energy() << " eta/phi: " << pos.eta() << "/" << pos.phi() << " deltaR: " << R;
#endif
    }
      }
    }
    if (etaAbs>1.25) {
      for (EcalRecHitCollection::const_iterator eItr=ecalEE->begin(); eItr!=ecalEE->end(); eItr++) {
    GlobalPoint pos = geo->getPosition(eItr->detid());
    double      R   = reco::deltaR(pos.eta(),pos.phi(),etaPhi.first,etaPhi.second);
    if (R < coneSize_) {
      nhitIn++;
      inEnergy += (eItr->energy());
      if (eItr->energy() > hitCountEthr_) nhitOut++;
      if (eItr->energy() > hitEthr_)      outEnergy += (eItr->energy());
#ifdef DebugLog
      edm::LogInfo("HcalIsoTrack") << "Rechit Close to the track has energy " << eItr->energy() << " eta/phi: " << pos.eta() << "/" << pos.phi() << " deltaR: " << R;
#endif
    }
      }
    }
#ifdef DebugLog
    edm::LogInfo("HcalIsoTrack") << "nhitIn:" << nhitIn << " inEnergy:" << inEnergy << " nhitOut:" << nhitOut << " outEnergy:" << outEnergy;
#endif
    reco::IsolatedPixelTrackCandidate newca(*isoPixTrackRefs[p]);
    newca.setEnergyIn(inEnergy);
    newca.setEnergyOut(outEnergy);
    newca.setNHitIn(nhitIn);
    newca.setNHitOut(nhitOut);
    iptcCollection->push_back(newca);   
  }
#ifdef DebugLog
  edm::LogInfo("HcalIsoTrack") << "ncand:" << nCand << " outcollction size:" << iptcCollection->size();
#endif
  iEvent.put(std::move(iptcCollection));
}

Member Data Documentation

const double IsolatedEcalPixelTrackCandidateProducer::coneSizeEta0_

private

Definition at line 32 of file IsolatedEcalPixelTrackCandidateProducer.h.

Referenced by produce().

const double IsolatedEcalPixelTrackCandidateProducer::coneSizeEta1_

private

Definition at line 33 of file IsolatedEcalPixelTrackCandidateProducer.h.

Referenced by produce().

const double IsolatedEcalPixelTrackCandidateProducer::hitCountEthr_

private

Definition at line 34 of file IsolatedEcalPixelTrackCandidateProducer.h.

Referenced by produce().

const double IsolatedEcalPixelTrackCandidateProducer::hitEthr_

private

Definition at line 35 of file IsolatedEcalPixelTrackCandidateProducer.h.

Referenced by produce().

const edm::EDGetTokenT<EcalRecHitCollection> IsolatedEcalPixelTrackCandidateProducer::tok_eb

private

Definition at line 30 of file IsolatedEcalPixelTrackCandidateProducer.h.

Referenced by produce().

const edm::EDGetTokenT<EcalRecHitCollection> IsolatedEcalPixelTrackCandidateProducer::tok_ee

private

Definition at line 29 of file IsolatedEcalPixelTrackCandidateProducer.h.

Referenced by produce().

const edm::EDGetTokenT<trigger::TriggerFilterObjectWithRefs> IsolatedEcalPixelTrackCandidateProducer::tok_trigcand

private

Definition at line 31 of file IsolatedEcalPixelTrackCandidateProducer.h.

Referenced by produce().