ParticleTowerProducer Class Reference

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

Inheritance diagram for ParticleTowerProducer:

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

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< reco::PFCandidateCollection >	src_
const bool	useHF_

Static Private Attributes
static constexpr int	ietaMax = 42

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

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

Definition at line 48 of file ParticleTowerProducer.cc.

Constructor & Destructor Documentation

ParticleTowerProducer()

ParticleTowerProducer::ParticleTowerProducer ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 70 of file ParticleTowerProducer.cc.

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

Member Function Documentation

eta2ieta()

int ParticleTowerProducer::eta2ieta ( double  eta ) const

private

Definition at line 126 of file ParticleTowerProducer.cc.

                                                    {
  // 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;
}

References PVValHelper::eta, hi::etaedge, LEDCalibrationChannels::ieta, and ietaMax.

Referenced by produce().

fillDescriptions()

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

static

Definition at line 186 of file ParticleTowerProducer.cc.

                                                                                       {
  // particleTowerProducer
  edm::ParameterSetDescription desc;
  desc.add<edm::InputTag>("src", edm::InputTag("particleFlow"));
  desc.add<bool>("useHF", true);
  descriptions.add("particleTowerProducer", desc);
}

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

ieta2eta()

double ParticleTowerProducer::ieta2eta ( int  ieta ) const

private

Definition at line 175 of file ParticleTowerProducer.cc.

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

References PVValHelper::eta, hi::etaedge, LEDCalibrationChannels::ieta, and Validation_hcalonly_cfi::sign.

Referenced by produce().

iphi2phi()

double ParticleTowerProducer::iphi2phi ( int  iphi, int  ieta  ) const

private

Definition at line 156 of file ParticleTowerProducer.cc.

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

References funct::abs(), LEDCalibrationChannels::ieta, LEDCalibrationChannels::iphi, M_PI, dqmiodumpmetadata::n, nphi, and phi.

Referenced by produce().

phi2iphi()

int ParticleTowerProducer::phi2iphi ( double  phi, int  ieta  ) const

private

Definition at line 140 of file ParticleTowerProducer.cc.

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

References funct::abs(), reco::ceil(), LEDCalibrationChannels::ieta, createfilelist::int, LEDCalibrationChannels::iphi, M_PI, angle0to2pi::make0To2pi(), dqmiodumpmetadata::n, nphi, and phi.

Referenced by produce().

produce()

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

overrideprivate

Definition at line 81 of file ParticleTowerProducer.cc.

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

References funct::abs(), SiStripBadComponentsDQMServiceTemplate_cfg::ep, EgHLTOffHistBins_cfi::et, PVValHelper::eta, eta2ieta(), LEDCalibrationChannels::ieta, ieta2eta(), iEvent, PixelMapPlotter::inputs, LEDCalibrationChannels::iphi, iphi2phi(), eostools::move(), p4, phi2iphi(), point, dumpMFGeometry_cfg::prod, src_, hgcalTowerProducer_cfi::tower, and useHF_.

Member Data Documentation

ietaMax

constexpr int ParticleTowerProducer::ietaMax = 42

staticconstexprprivate

Definition at line 65 of file ParticleTowerProducer.cc.

Referenced by eta2ieta().

src_

edm::EDGetTokenT<reco::PFCandidateCollection> ParticleTowerProducer::src_

private

Definition at line 61 of file ParticleTowerProducer.cc.

Referenced by produce().

useHF_

const bool ParticleTowerProducer::useHF_

private

Definition at line 62 of file ParticleTowerProducer.cc.

Referenced by produce().