GenPartIsoProducer Class Reference

Inheritance diagram for GenPartIsoProducer:

Public Member Functions
	GenPartIsoProducer (const edm::ParameterSet &iConfig)
	~GenPartIsoProducer () 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 Member Functions
float	computeIso (TLorentzVector thisPart, edm::Handle< pat::PackedGenParticleCollection > packedGenParticles, std::set< int > gen_fsrset, bool skip_leptons)
void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
int	additionalPdgId_
edm::EDGetTokenT< reco::GenParticleCollection >	finalGenParticleToken
std::vector< float >	Lepts_RelIso
edm::EDGetTokenT< pat::PackedGenParticleCollection >	packedGenParticlesToken

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

Definition at line 31 of file GenPartIsoProducer.cc.

Constructor & Destructor Documentation

GenPartIsoProducer()

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

inlineexplicit

Definition at line 33 of file GenPartIsoProducer.cc.

      : finalGenParticleToken(consumes<reco::GenParticleCollection>(iConfig.getParameter<edm::InputTag>("genPart"))),
        packedGenParticlesToken(
            consumes<pat::PackedGenParticleCollection>(iConfig.getParameter<edm::InputTag>("packedGenPart"))),
        additionalPdgId_(iConfig.getParameter<int>("additionalPdgId")) {
    produces<edm::ValueMap<float>>();
  }

~GenPartIsoProducer()

GenPartIsoProducer::~GenPartIsoProducer ( )

override

Definition at line 56 of file GenPartIsoProducer.cc.

{}

Member Function Documentation

computeIso()

float GenPartIsoProducer::computeIso ( TLorentzVector  thisPart, edm::Handle< pat::PackedGenParticleCollection >  packedGenParticles, std::set< int >  gen_fsrset, bool  skip_leptons  )

private

Definition at line 125 of file GenPartIsoProducer.cc.

References funct::abs(), reco::deltaR(), dqmdumpme::k, packedGenParticles_cfi::packedGenParticles, and EgammaValidation_cff::pdgId.

                                                        {
  double this_GENiso = 0.0;
  for (size_t k = 0; k < packedGenParticles->size(); k++) {
    if ((*packedGenParticles)[k].status() != 1)
      continue;
    if (abs((*packedGenParticles)[k].pdgId()) == 12 || abs((*packedGenParticles)[k].pdgId()) == 14 ||
        abs((*packedGenParticles)[k].pdgId()) == 16)
      continue;
    if (reco::deltaR(thisPart.Eta(), thisPart.Phi(), (*packedGenParticles)[k].eta(), (*packedGenParticles)[k].phi()) <
        0.001)
      continue;
    if (skip_leptons == true) {
      if ((abs((*packedGenParticles)[k].pdgId()) == 11 || abs((*packedGenParticles)[k].pdgId()) == 13))
        continue;
      if (gen_fsrset.find(k) != gen_fsrset.end())
        continue;
    }
    double this_dRvL_nolep =
        reco::deltaR(thisPart.Eta(), thisPart.Phi(), (*packedGenParticles)[k].eta(), (*packedGenParticles)[k].phi());
    if (this_dRvL_nolep < 0.3) {
      this_GENiso = this_GENiso + (*packedGenParticles)[k].pt();
    }
  }
  this_GENiso = this_GENiso / thisPart.Pt();
  return this_GENiso;
}

fillDescriptions()

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

static

Definition at line 155 of file GenPartIsoProducer.cc.

References edm::ConfigurationDescriptions::addDefault(), and submitPVResolutionJobs::desc.

                                                                                    {
  // Description of external inputs requered by this module
  // genPart: collection of gen particles for which to compute Iso
  // packedGenPart: collection of particles to be used for leptons dressing
  // additionalPdgId: additional particle (besides leptons) for which Iso is computed
  edm::ParameterSetDescription desc;
  desc.add<edm::InputTag>("genPart")->setComment("input physics object collection");
  desc.add<edm::InputTag>("packedGenPart")->setComment("input stable hadrons collection");
  desc.add<int>("additionalPdgId")->setComment("additional pdgId for Iso computation");
  descriptions.addDefault(desc);
}

produce()

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

overrideprivate

Definition at line 58 of file GenPartIsoProducer.cc.

References funct::abs(), reco::deltaR(), hcalRecHitTable_cff::energy, PVValHelper::eta, CustomPhysics_cfi::gamma, genparticles_cff::genPart, iEvent, dqmdumpme::k, visualization-live-secondInstance_cfg::m, reco::Candidate::mother(), eostools::move(), reco::Candidate::numberOfMothers(), packedGenParticles_cfi::packedGenParticles, and DiDispStaMuonMonitor_cfi::pt.

                                                                              {
  using namespace edm;

  auto finalParticles = iEvent.getHandle(finalGenParticleToken);
  auto packedGenParticles = iEvent.getHandle(packedGenParticlesToken);

  reco::GenParticleCollection::const_iterator genPart;

  Lepts_RelIso.clear();

  for (genPart = finalParticles->begin(); genPart != finalParticles->end(); genPart++) {
    if (abs(genPart->pdgId()) == 11 || abs(genPart->pdgId()) == 13 || abs(genPart->pdgId()) == 15) {
      TLorentzVector lep_dressed;
      lep_dressed.SetPtEtaPhiE(genPart->pt(), genPart->eta(), genPart->phi(), genPart->energy());
      std::set<int> gen_fsrset;
      for (size_t k = 0; k < packedGenParticles->size(); k++) {
        if ((*packedGenParticles)[k].status() != 1)
          continue;
        if ((*packedGenParticles)[k].pdgId() != 22)
          continue;
        double this_dR_lgamma = reco::deltaR(
            genPart->eta(), genPart->phi(), (*packedGenParticles)[k].eta(), (*packedGenParticles)[k].phi());
        bool idmatch = false;
        if ((*packedGenParticles)[k].mother(0)->pdgId() == genPart->pdgId())
          idmatch = true;
        const reco::Candidate* mother = (*packedGenParticles)[k].mother(0);
        for (size_t m = 0; m < mother->numberOfMothers(); m++) {
          if ((*packedGenParticles)[k].mother(m)->pdgId() == genPart->pdgId())
            idmatch = true;
        }
        if (!idmatch)
          continue;
        if (this_dR_lgamma < 0.3) {
          gen_fsrset.insert(k);
          TLorentzVector gamma;
          gamma.SetPtEtaPhiE((*packedGenParticles)[k].pt(),
                             (*packedGenParticles)[k].eta(),
                             (*packedGenParticles)[k].phi(),
                             (*packedGenParticles)[k].energy());
          lep_dressed = lep_dressed + gamma;
        }
      }
      float this_GENiso = 0.0;
      TLorentzVector thisLep;
      thisLep.SetPtEtaPhiM(lep_dressed.Pt(), lep_dressed.Eta(), lep_dressed.Phi(), lep_dressed.M());
      this_GENiso = computeIso(thisLep, packedGenParticles, gen_fsrset, true);
      Lepts_RelIso.push_back(this_GENiso);
    } else if (abs(genPart->pdgId()) == additionalPdgId_) {
      float this_GENiso = 0.0;
      std::set<int> gen_fsrset_nolep;
      TLorentzVector thisPart;
      thisPart.SetPtEtaPhiE(genPart->pt(), genPart->eta(), genPart->phi(), genPart->energy());
      this_GENiso = computeIso(thisPart, packedGenParticles, gen_fsrset_nolep, false);
      Lepts_RelIso.push_back(this_GENiso);
    } else {
      float this_GENiso = 0.0;
      Lepts_RelIso.push_back(this_GENiso);
    }
  }

  auto isoV = std::make_unique<edm::ValueMap<float>>();
  edm::ValueMap<float>::Filler fillerIsoMap(*isoV);
  fillerIsoMap.insert(finalParticles, Lepts_RelIso.begin(), Lepts_RelIso.end());
  fillerIsoMap.fill();
  iEvent.put(std::move(isoV));
}

Member Data Documentation

additionalPdgId_

int GenPartIsoProducer::additionalPdgId_

private

Definition at line 53 of file GenPartIsoProducer.cc.

finalGenParticleToken

edm::EDGetTokenT<reco::GenParticleCollection> GenPartIsoProducer::finalGenParticleToken

private

Definition at line 51 of file GenPartIsoProducer.cc.

Lepts_RelIso

std::vector<float> GenPartIsoProducer::Lepts_RelIso

private

Definition at line 50 of file GenPartIsoProducer.cc.

packedGenParticlesToken

edm::EDGetTokenT<pat::PackedGenParticleCollection> GenPartIsoProducer::packedGenParticlesToken

private

Definition at line 52 of file GenPartIsoProducer.cc.