ParticleTowerProducer< T > Class Template Reference

#include <RecoHI/HiJetAlgos/plugins/ParticleTowerProducer.cc>

Inheritance diagram for ParticleTowerProducer< T >:

Public Member Functions
	ParticleTowerProducer (const edm::ParameterSet &)
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
int	eta2ieta (double eta) const
double	ieta2eta (int ieta) const
double	iphi2phi (int iphi, int ieta) const
int	phi2iphi (double phi, int ieta) const
void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
edm::EDGetTokenT< typename edm::View< T > >	src_
const bool	useHF_

Static Private Attributes
static constexpr int	ietaMax = 42

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

template<class T>
class ParticleTowerProducer< T >

Description: [one line class summary] Implementation: [Notes on implementation]

Definition at line 49 of file ParticleTowerProducer.cc.

Constructor & Destructor Documentation

template<class T >

ParticleTowerProducer< T >::ParticleTowerProducer ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 72 of file ParticleTowerProducer.cc.

    : src_(consumes<typename edm::View<T> >(iConfig.getParameter<edm::InputTag>("src"))),
      useHF_(iConfig.getParameter<bool>("useHF")) {
  produces<CaloTowerCollection>();
}

Member Function Documentation

template<class T >

int ParticleTowerProducer< T >::eta2ieta ( double  eta ) const

private

Definition at line 130 of file ParticleTowerProducer.cc.

References hi::etaedge.

                                                       {
  // binary search in the array of towers eta edges

  int ieta = 1;
  double xeta = fabs(eta);
  while (xeta > hi::etaedge[ieta] && ieta < ietaMax - 1) {
    ++ieta;
  }

  if (eta < 0)
    ieta = -ieta;
  return ieta;
}

template<class T >

void ParticleTowerProducer< T >::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

static

Definition at line 194 of file ParticleTowerProducer.cc.

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

                                                                                          {
  edm::ParameterSetDescription desc;
  desc.add<bool>("useHF", true);
  if (typeid(T) == typeid(reco::PFCandidate)) {
    desc.add<edm::InputTag>("src", edm::InputTag("particleFlow"));
    descriptions.add("PFTowers", desc);
  }
  if (typeid(T) == typeid(pat::PackedCandidate)) {
    desc.add<edm::InputTag>("src", edm::InputTag("packedPFCandidates"));
    descriptions.add("PackedPFTowers", desc);
  }
}

template<class T >

double ParticleTowerProducer< T >::ieta2eta ( int  ieta ) const

private

Definition at line 182 of file ParticleTowerProducer.cc.

References PVValHelper::eta, hi::etaedge, and jetcorrextractor::sign().

                                                        {
  int sign = 1;
  if (ieta < 0) {
    sign = -1;
    ieta = -ieta;
  }

  double eta = sign * (hi::etaedge[ieta] + hi::etaedge[ieta - 1]) / 2.;
  return eta;
}

template<class T >

double ParticleTowerProducer< T >::iphi2phi ( int  iphi, int  ieta  ) const

private

Definition at line 162 of file ParticleTowerProducer.cc.

References funct::abs(), M_PI, dqmiodumpmetadata::n, nphi, and PVValHelper::phi.

                                                                  {
  double phi = 0;
  int nphi = 72;

  int n = 1;
  if (abs(ieta) > 20)
    n = 2;
  if (abs(ieta) >= 40)
    n = 4;

  int myphi = (iphi - 1) / n + 1;

  phi = 2. * M_PI * (myphi - 0.5) / nphi * n;
  while (phi > M_PI)
    phi -= 2. * M_PI;

  return phi;
}

template<class T >

int ParticleTowerProducer< T >::phi2iphi ( double  phi, int  ieta  ) const

private

Definition at line 145 of file ParticleTowerProducer.cc.

References funct::abs(), reco::ceil(), M_PI, angle0to2pi::make0To2pi(), dqmiodumpmetadata::n, and nphi.

                                                                 {
  phi = angle0to2pi::make0To2pi(phi);
  int nphi = 72;
  int n = 1;
  if (abs(ieta) > 20)
    n = 2;
  if (abs(ieta) >= 40)
    n = 4;

  int iphi = (int)std::ceil(phi / 2.0 / M_PI * nphi / n);

  iphi = n * (iphi - 1) + 1;

  return iphi;
}

template<class T >

void ParticleTowerProducer< T >::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

overrideprivate

Definition at line 84 of file ParticleTowerProducer.cc.

References funct::abs(), PVValHelper::eta, edm::Event::get(), PixelMapPlotter::inputs, eostools::move(), point, edm::Event::put(), and hgcalTowerProducer_cfi::tower.

                                                                                    {
  using namespace edm;

  typedef std::pair<int, int> EtaPhi;
  typedef std::map<EtaPhi, double> EtaPhiMap;
  EtaPhiMap towers_;
  towers_.clear();

  auto const& inputs = iEvent.get(src_);
  for (auto const& particle : inputs) {
    double eta = particle.eta();

    int ieta = eta2ieta(eta);
    int iphi = phi2iphi(particle.phi(), ieta);

    if (!useHF_ && abs(ieta) > 29)
      continue;

    EtaPhi ep(ieta, iphi);
    towers_[ep] += particle.et();
  }

  auto prod = std::make_unique<CaloTowerCollection>();
  prod->reserve(towers_.size());
  for (auto const& tower : towers_) {
    EtaPhi ep = tower.first;
    double et = tower.second;

    int ieta = ep.first;
    int iphi = ep.second;

    CaloTowerDetId newTowerId(ieta, iphi);  // totally dummy id

    // currently sets et = pt, mass to zero
    // pt, eta, phi, mass
    reco::Particle::PolarLorentzVector p4(et, ieta2eta(ieta), iphi2phi(iphi, ieta), 0.);
    GlobalPoint point(p4.x(), p4.y(), p4.z());
    prod->emplace_back(newTowerId, p4.e(), 0, 0, 0, 0, p4, point, point);
  }

  iEvent.put(std::move(prod));
}

Member Data Documentation

template<class T >

constexpr int ParticleTowerProducer< T >::ietaMax = 42

staticprivate

Definition at line 66 of file ParticleTowerProducer.cc.

template<class T >

edm::EDGetTokenT<typename edm::View<T> > ParticleTowerProducer< T >::src_

private

Definition at line 62 of file ParticleTowerProducer.cc.

template<class T >

const bool ParticleTowerProducer< T >::useHF_

private

Definition at line 63 of file ParticleTowerProducer.cc.