HFJetShowerShape Class Reference

#include <HFJetShowerShape/HFJetShowerShape/plugins/HFJetShowerShape.cc>

Inheritance diagram for HFJetShowerShape:

Public Member Functions
	HFJetShowerShape (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
void	produce (edm::Event &, const edm::EventSetup &) override
template<typename T >
void	putInEvent (const std::string &, const edm::Handle< edm::View< reco::Jet >> &, std::vector< T >, edm::Event &) const
	Function to put product into event. More...

Private Attributes
const double	hfTowerEtaWidth_
const double	hfTowerPhiWidth_
const double	jetEtaThreshold_
const double	jetPtThreshold_
const double	jetReferenceRadius_
const edm::EDGetTokenT< edm::View< reco::Jet > >	jets_token_
const double	offsetPerPU_
const double	stripPtThreshold_
const double	vertexRecoEffcy_
const edm::EDGetTokenT< std::vector< reco::Vertex > >	vertices_token_
const double	widthPtThreshold_

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

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 43 of file HFJetShowerShape.cc.

Constructor & Destructor Documentation

HFJetShowerShape()

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

explicit

Definition at line 78 of file HFJetShowerShape.cc.

    : jets_token_(consumes<edm::View<reco::Jet>>(iConfig.getParameter<edm::InputTag>("jets"))),
      vertices_token_(consumes<std::vector<reco::Vertex>>(iConfig.getParameter<edm::InputTag>("vertices"))),
      jetPtThreshold_(iConfig.getParameter<double>("jetPtThreshold")),
      jetEtaThreshold_(iConfig.getParameter<double>("jetEtaThreshold")),
      hfTowerEtaWidth_(iConfig.getParameter<double>("hfTowerEtaWidth")),
      hfTowerPhiWidth_(iConfig.getParameter<double>("hfTowerPhiWidth")),
      vertexRecoEffcy_(iConfig.getParameter<double>("vertexRecoEffcy")),
      offsetPerPU_(iConfig.getParameter<double>("offsetPerPU")),
      jetReferenceRadius_(iConfig.getParameter<double>("jetReferenceRadius")),
      stripPtThreshold_(iConfig.getParameter<double>("stripPtThreshold")),
      widthPtThreshold_(iConfig.getParameter<double>("widthPtThreshold")) {
  produces<edm::ValueMap<float>>("sigmaEtaEta");
  produces<edm::ValueMap<float>>("sigmaPhiPhi");
  produces<edm::ValueMap<int>>("centralEtaStripSize");
  produces<edm::ValueMap<int>>("adjacentEtaStripsSize");
}

Member Function Documentation

fillDescriptions()

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

static

Definition at line 212 of file HFJetShowerShape.cc.

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

                                                                                  {
  edm::ParameterSetDescription desc;
  desc.add<edm::InputTag>("jets", edm::InputTag("ak4PFJetsCHS"));
  desc.add<edm::InputTag>("vertices", edm::InputTag("offlineSlimmedPrimaryVertices"));
  desc.add<double>("jetPtThreshold", 25.);
  desc.add<double>("jetEtaThreshold", 2.9);
  desc.add<double>("hfTowerEtaWidth", 0.175);
  desc.add<double>("hfTowerPhiWidth", 0.175);
  desc.add<double>("vertexRecoEffcy", 0.7);
  desc.add<double>("offsetPerPU", 0.4);
  desc.add<double>("jetReferenceRadius", 0.4);
  desc.add<double>("stripPtThreshold", 10.);
  desc.add<double>("widthPtThreshold", 3.);
  descriptions.add("hfJetShowerShape", desc);
}

produce()

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

overrideprivate

Definition at line 101 of file HFJetShowerShape.cc.

References funct::abs(), SiPixelRawToDigiRegional_cfi::deltaPhi, reco::Candidate::eta(), hfTowerEtaWidth_, hfTowerPhiWidth_, mps_fire::i, iEvent, metsig::jet, jetEtaThreshold_, jetPtThreshold_, jetReferenceRadius_, jets_token_, M_PI, nPV, offsetPerPU_, reco::Candidate::pdgId(), reco::Candidate::phi(), reco::Candidate::pt(), putInEvent(), mathSSE::sqrt(), stripPtThreshold_, vertexRecoEffcy_, vertices_token_, mps_merge::weight, and widthPtThreshold_.

                                                                            {
  using namespace edm;

  //Jets
  auto theJets = iEvent.getHandle(jets_token_);

  //Vertices
  int nPV = iEvent.get(vertices_token_).size();

  //Products
  std::vector<float> v_sigmaEtaEta, v_sigmaPhiPhi;
  v_sigmaEtaEta.reserve(theJets->size());
  v_sigmaPhiPhi.reserve(theJets->size());
  std::vector<int> v_size_CentralEtaStrip, v_size_AdjacentEtaStrips;
  v_size_CentralEtaStrip.reserve(theJets->size());
  v_size_AdjacentEtaStrips.reserve(theJets->size());

  //Et offset for HF PF candidates
  double puoffset = offsetPerPU_ / (M_PI * jetReferenceRadius_ * jetReferenceRadius_) * nPV / vertexRecoEffcy_ *
                    (hfTowerEtaWidth_ * hfTowerPhiWidth_);

  for (auto const& jet : *theJets) {
    double pt_jet = jet.pt();
    double eta_jet = jet.eta();

    //If central or low pt jets, fill with dummy variables
    if (pt_jet <= jetPtThreshold_ || std::abs(eta_jet) <= jetEtaThreshold_) {
      v_sigmaEtaEta.push_back(-1.);
      v_sigmaPhiPhi.push_back(-1.);
      v_size_CentralEtaStrip.push_back(0);
      v_size_AdjacentEtaStrips.push_back(0);
    } else {
      //First loop over PF candidates to compute some global variables needed for shower shape calculations
      double sumptPFcands = 0.;

      for (unsigned i = 0; i < jet.numberOfSourceCandidatePtrs(); ++i) {
        const reco::Candidate* icand = jet.sourceCandidatePtr(i).get();
        //Only look at pdgId =1,2 (HF PF cands)
        if (std::abs(icand->pdgId()) > 2)
          continue;
        double pt_PUsub = icand->pt() - puoffset;
        if (pt_PUsub < widthPtThreshold_)
          continue;
        sumptPFcands += pt_PUsub;
      }

      //Second loop to compute the various shower shape variables
      int size_CentralEtaStrip(0.), size_AdjacentEtaStrips(0.);
      double sigmaEtaEtaSq(0.), sigmaPhiPhiSq(0.);
      double sumweightsPFcands = 0;
      for (unsigned i = 0; i < jet.numberOfSourceCandidatePtrs(); ++i) {
        const reco::Candidate* icand = jet.sourceCandidatePtr(i).get();
        //Only look at pdgId =1,2 (HF PF cands)
        if (std::abs(icand->pdgId()) > 2)
          continue;

        double deta = std::abs(icand->eta() - jet.eta());
        double dphi = std::abs(deltaPhi(icand->phi(), jet.phi()));
        double pt_PUsub = icand->pt() - puoffset;

        //This is simply the size of the central eta strip and the adjacent strips
        if (pt_PUsub >= stripPtThreshold_) {
          if (dphi <= hfTowerPhiWidth_ * 0.5)
            size_CentralEtaStrip++;
          else if (dphi <= hfTowerPhiWidth_ * 1.5)
            size_AdjacentEtaStrips++;
        }

        //Now computing sigmaEtaEta/PhiPhi
        if (pt_PUsub >= widthPtThreshold_ && sumptPFcands > 0) {
          double weight = pt_PUsub / sumptPFcands;
          sigmaEtaEtaSq += deta * deta * weight;
          sigmaPhiPhiSq += dphi * dphi * weight;
          sumweightsPFcands += weight;
        }
      }

      v_size_CentralEtaStrip.push_back(size_CentralEtaStrip);
      v_size_AdjacentEtaStrips.push_back(size_AdjacentEtaStrips);

      if (sumweightsPFcands > 0 && sigmaEtaEtaSq > 0 && sigmaPhiPhiSq > 0) {
        v_sigmaEtaEta.push_back(sqrt(sigmaEtaEtaSq / sumweightsPFcands));
        v_sigmaPhiPhi.push_back(sqrt(sigmaPhiPhiSq / sumweightsPFcands));
      } else {
        v_sigmaEtaEta.push_back(-1.);
        v_sigmaPhiPhi.push_back(-1.);
      }

    }  //End loop over jets
  }

  putInEvent("sigmaEtaEta", theJets, v_sigmaEtaEta, iEvent);
  putInEvent("sigmaPhiPhi", theJets, v_sigmaPhiPhi, iEvent);
  putInEvent("centralEtaStripSize", theJets, v_size_CentralEtaStrip, iEvent);
  putInEvent("adjacentEtaStripsSize", theJets, v_size_AdjacentEtaStrips, iEvent);
}

putInEvent()

template<typename T >

void HFJetShowerShape::putInEvent ( const std::string &  name, const edm::Handle< edm::View< reco::Jet >> &  jets, std::vector< T >  product, edm::Event &  iEvent  ) const

private

Function to put product into event.

Definition at line 200 of file HFJetShowerShape.cc.

References trigObjTnPSource_cfi::filler, iEvent, singleTopDQM_cfi::jets, eostools::move(), Skims_PA_cff::name, and MillePedeFileConverter_cfg::out.

Referenced by produce().

                                                          {
  auto out = std::make_unique<edm::ValueMap<T>>();
  typename edm::ValueMap<T>::Filler filler(*out);
  filler.insert(jets, product.begin(), product.end());
  filler.fill();
  iEvent.put(std::move(out), name);
}

Member Data Documentation

hfTowerEtaWidth_

const double HFJetShowerShape::hfTowerEtaWidth_

private

Definition at line 60 of file HFJetShowerShape.cc.

Referenced by produce().

hfTowerPhiWidth_

const double HFJetShowerShape::hfTowerPhiWidth_

private

Definition at line 60 of file HFJetShowerShape.cc.

Referenced by produce().

jetEtaThreshold_

const double HFJetShowerShape::jetEtaThreshold_

private

Definition at line 58 of file HFJetShowerShape.cc.

Referenced by produce().

jetPtThreshold_

const double HFJetShowerShape::jetPtThreshold_

private

Definition at line 58 of file HFJetShowerShape.cc.

Referenced by produce().

jetReferenceRadius_

const double HFJetShowerShape::jetReferenceRadius_

private

Definition at line 62 of file HFJetShowerShape.cc.

Referenced by produce().

jets_token_

const edm::EDGetTokenT<edm::View<reco::Jet> > HFJetShowerShape::jets_token_

private

Definition at line 54 of file HFJetShowerShape.cc.

Referenced by produce().

offsetPerPU_

const double HFJetShowerShape::offsetPerPU_

private

Definition at line 62 of file HFJetShowerShape.cc.

Referenced by produce().

stripPtThreshold_

const double HFJetShowerShape::stripPtThreshold_

private

Definition at line 64 of file HFJetShowerShape.cc.

Referenced by produce().

vertexRecoEffcy_

const double HFJetShowerShape::vertexRecoEffcy_

private

Definition at line 62 of file HFJetShowerShape.cc.

Referenced by produce().

vertices_token_

const edm::EDGetTokenT<std::vector<reco::Vertex> > HFJetShowerShape::vertices_token_

private

Definition at line 55 of file HFJetShowerShape.cc.

Referenced by produce().

widthPtThreshold_

const double HFJetShowerShape::widthPtThreshold_

private

Definition at line 64 of file HFJetShowerShape.cc.

Referenced by produce().