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

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

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

produce()

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

overrideprivate

Definition at line 101 of file HFJetShowerShape.cc.

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

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

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.

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

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

Referenced by produce().

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