JetPlusTrackAddonSeedProducer Class Reference

#include <JetPlusTrackAddonSeedProducer.cc>

Inheritance diagram for JetPlusTrackAddonSeedProducer:

Public Member Functions
	JetPlusTrackAddonSeedProducer (const edm::ParameterSet &)
void	produce (edm::Event &, const edm::EventSetup &) override
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 Attributes
const double	dRcone_
const edm::EDGetTokenT< CaloTowerCollection >	input_ctw_token_
const edm::EDGetTokenT< edm::View< reco::CaloJet > >	input_jets_token_
const edm::EDGetTokenT< edm::View< reco::TrackJet > >	input_trackjets_token_
const edm::EDGetTokenT< reco::VertexCollection >	input_vertex_token_
const edm::EDGetTokenT< std::vector< pat::PackedCandidate > >	tokenPFCandidates_
const bool	usePAT_

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 48 of file JetPlusTrackAddonSeedProducer.cc.

Constructor & Destructor Documentation

JetPlusTrackAddonSeedProducer()

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

explicit

Definition at line 64 of file JetPlusTrackAddonSeedProducer.cc.

    : dRcone_(iConfig.getParameter<double>("dRcone")),
      usePAT_(iConfig.getParameter<bool>("UsePAT")),
      input_jets_token_(consumes<edm::View<reco::CaloJet> >(iConfig.getParameter<edm::InputTag>("srcCaloJets"))),
      input_trackjets_token_(consumes<edm::View<reco::TrackJet> >(iConfig.getParameter<edm::InputTag>("srcTrackJets"))),
      input_vertex_token_(consumes<reco::VertexCollection>(iConfig.getParameter<edm::InputTag>("srcPVs"))),
      tokenPFCandidates_(consumes<pat::PackedCandidateCollection>(iConfig.getParameter<edm::InputTag>("PFCandidates"))),
      input_ctw_token_(consumes<CaloTowerCollection>(iConfig.getParameter<edm::InputTag>("towerMaker"))) {
  //register your products
  produces<reco::CaloJetCollection>();
}

Member Function Documentation

fillDescriptions()

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

static

Definition at line 76 of file JetPlusTrackAddonSeedProducer.cc.

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

                                                                                                {
  edm::ParameterSetDescription desc;
  desc.add<double>("dRcone", 0.4);
  desc.add<bool>("UsePAT", false);
  desc.add<edm::InputTag>("srcCaloJets", edm::InputTag("ak4CaloJets"));
  desc.add<edm::InputTag>("srcTrackJets", edm::InputTag("ak4TrackJets"));
  desc.add<edm::InputTag>("srcPVs", edm::InputTag("primaryVertex"));
  desc.add<edm::InputTag>("PFCandidates", edm::InputTag("PFCandidates"));
  desc.add<edm::InputTag>("towerMaker", edm::InputTag("towerMaker"));
  iDescriptions.addWithDefaultLabel(desc);
}

produce()

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

override

Definition at line 89 of file JetPlusTrackAddonSeedProducer.cc.

References funct::abs(), HLTMuonOfflineAnalyzer_cfi::deltaR2, dRcone_, iEvent, input_ctw_token_, input_jets_token_, input_trackjets_token_, metsig::jet, PDWG_EXODelayedJetMET_cff::jets, M_PI, eostools::move(), reco::CaloJet::Specific::mTowersArea, packedPFCandidateRefMixer_cfi::pf, zmumugammaAnalyzer_cfi::pfCandidates, reco::Jet::setJetArea(), tokenPFCandidates_, and usePAT_.

                                                                                         {
  using namespace edm;
  // get stuff from Event
  auto const& jets = iEvent.get(input_jets_token_);
  auto const& jetsTrackJets = iEvent.get(input_trackjets_token_);
  auto pCaloOut = std::make_unique<reco::CaloJetCollection>();

  for (auto const& jet : jetsTrackJets) {
    int iflag = 0;
    for (auto const& oldjet : jets) {
      double dr2 = deltaR2(jet, oldjet);
      if (dr2 < dRcone_ * dRcone_) {
        iflag = 1;
      }
    }  // Calojets

    if (iflag == 1)
      continue;
    double caloen = 0.;
    double hadinho = 0.;
    double hadinhb = 0.;
    double hadinhe = 0.;
    double hadinhf = 0.;
    double emineb = 0.;
    double eminee = 0.;
    double eminhf = 0.;
    double eefraction = 0.;
    double hhfraction = 0.;
    int ncand = 0;

    if (usePAT_) {
      auto const& pfCandidates = iEvent.get(tokenPFCandidates_);
      for (auto const& pf : pfCandidates) {
        double dr2 = deltaR2(jet, pf);
        if (dr2 > dRcone_ * dRcone_)
          continue;
        // jetconstit
        caloen = caloen + pf.energy() * pf.caloFraction();
        hadinho += 0.;
        if (std::abs(pf.eta()) <= 1.4) {
          hadinhb += pf.energy() * pf.caloFraction() * pf.hcalFraction();
          emineb += pf.energy() * pf.caloFraction() * (1. - pf.hcalFraction());
        } else if (std::abs(pf.eta()) < 3.) {
          hadinhe += pf.energy() * pf.caloFraction() * pf.hcalFraction();
          eminee += pf.energy() * pf.caloFraction() * (1. - pf.hcalFraction());
        } else {
          hadinhf += pf.energy() * pf.caloFraction() * pf.hcalFraction();
          eminhf += pf.energy() * pf.caloFraction() * (1. - pf.hcalFraction());
        }
        ncand++;
      }  // pfcandidates
    } else {
      auto const& cts = iEvent.get(input_ctw_token_);
      for (auto const& ct : cts) {
        double dr2 = deltaR2(jet, ct);
        if (dr2 > dRcone_ * dRcone_)
          continue;
        caloen = caloen + ct.energy();
        hadinho += ct.energyInHO();
        hadinhb += ct.energyInHB();
        hadinhe += ct.energyInHE();
        hadinhf += 0.5 * ct.energyInHF();
        emineb += ct.energy() - ct.energyInHB() - ct.energyInHO();
        eminee += ct.energy() - ct.energyInHE();
        eminhf += 0.5 * ct.energyInHF();
        ncand++;
      }
    }
    eefraction = (emineb + eminee) / caloen;
    hhfraction = (hadinhb + hadinhe + hadinhf + hadinho) / caloen;

    double trackp = jet.p();
    if (caloen <= 0.)
      caloen = 0.001;
    math::XYZTLorentzVector pcalo4(caloen * jet.p4() / trackp);
    reco::CaloJet::Specific calospe;
    calospe.mTowersArea = -1 * ncand;
    calospe.mHadEnergyInHO = hadinho;
    calospe.mHadEnergyInHB = hadinhb;
    calospe.mHadEnergyInHE = hadinhe;
    calospe.mHadEnergyInHF = hadinhf;
    calospe.mEmEnergyInEB = emineb;
    calospe.mEmEnergyInEE = eminee;
    calospe.mEmEnergyInHF = eminhf;
    calospe.mEnergyFractionEm = eefraction / caloen;
    calospe.mEnergyFractionHadronic = hhfraction / caloen;

    reco::CaloJet mycalo(pcalo4, jet.primaryVertex()->position(), calospe);
    mycalo.setJetArea(M_PI * dRcone_ * dRcone_);
    pCaloOut->push_back(mycalo);

  }  // trackjets
  iEvent.put(std::move(pCaloOut));
}

Member Data Documentation

dRcone_

const double JetPlusTrackAddonSeedProducer::dRcone_

private

Definition at line 55 of file JetPlusTrackAddonSeedProducer.cc.

Referenced by produce().

input_ctw_token_

const edm::EDGetTokenT<CaloTowerCollection> JetPlusTrackAddonSeedProducer::input_ctw_token_

private

Definition at line 61 of file JetPlusTrackAddonSeedProducer.cc.

Referenced by produce().

input_jets_token_

const edm::EDGetTokenT<edm::View<reco::CaloJet> > JetPlusTrackAddonSeedProducer::input_jets_token_

private

Definition at line 57 of file JetPlusTrackAddonSeedProducer.cc.

Referenced by produce().

input_trackjets_token_

const edm::EDGetTokenT<edm::View<reco::TrackJet> > JetPlusTrackAddonSeedProducer::input_trackjets_token_

private

Definition at line 58 of file JetPlusTrackAddonSeedProducer.cc.

Referenced by produce().

input_vertex_token_

const edm::EDGetTokenT<reco::VertexCollection> JetPlusTrackAddonSeedProducer::input_vertex_token_

private

Definition at line 59 of file JetPlusTrackAddonSeedProducer.cc.

tokenPFCandidates_

const edm::EDGetTokenT<std::vector<pat::PackedCandidate> > JetPlusTrackAddonSeedProducer::tokenPFCandidates_

private

Definition at line 60 of file JetPlusTrackAddonSeedProducer.cc.

Referenced by produce().

usePAT_

const bool JetPlusTrackAddonSeedProducer::usePAT_

private

Definition at line 56 of file JetPlusTrackAddonSeedProducer.cc.

Referenced by produce().