pat::PATGenJetSlimmer Class Reference

Matcher of reconstructed objects to L1 Muons. More...

Inheritance diagram for pat::PATGenJetSlimmer:

Public Member Functions
	PATGenJetSlimmer (const edm::ParameterSet &iConfig)
void	produce (edm::Event &iEvent, const edm::EventSetup &iSetup) override
	~PATGenJetSlimmer () 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

Private Attributes
const bool	clearDaughters_
	reset daughters to an empty vector More...
const StringCutObjectSelector < reco::GenJet >	cut_
const StringCutObjectSelector < reco::GenJet >	cutLoose_
const bool	dropSpecific_
	drop the specific More...
const edm::EDGetTokenT < edm::Association < std::vector < pat::PackedGenParticle > > >	gp2pgp_
const unsigned	nLoose_
const edm::EDGetTokenT < edm::View< reco::GenJet > >	src_

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

Matcher of reconstructed objects to L1 Muons.

Author

Giovanni Petrucciani

Definition at line 25 of file PATGenJetSlimmer.cc.

Constructor & Destructor Documentation

pat::PATGenJetSlimmer::PATGenJetSlimmer ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 48 of file PATGenJetSlimmer.cc.

    : src_(consumes<edm::View<reco::GenJet> >(iConfig.getParameter<edm::InputTag>("src"))),
      gp2pgp_(consumes<edm::Association<std::vector<pat::PackedGenParticle> > >(
          iConfig.getParameter<edm::InputTag>("packedGenParticles"))),
      cut_(iConfig.getParameter<std::string>("cut")),
      cutLoose_(iConfig.getParameter<std::string>("cutLoose")),
      nLoose_(iConfig.getParameter<unsigned>("nLoose")),
      clearDaughters_(iConfig.getParameter<bool>("clearDaughters")),
      dropSpecific_(iConfig.getParameter<bool>("dropSpecific")) {
  produces<std::vector<reco::GenJet> >();
  produces<edm::Association<std::vector<reco::GenJet> > >("slimmedGenJetAssociation");
}

pat::PATGenJetSlimmer::~PATGenJetSlimmer ( )

inlineoverride

Definition at line 28 of file PATGenJetSlimmer.cc.

{}

Member Function Documentation

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

override

Definition at line 61 of file PATGenJetSlimmer.cc.

References reco::CompositePtrCandidate::addDaughter(), reco::CompositePtrCandidate::clearDaughters(), reco::CompositePtrCandidate::daughterPtrVector(), edm::Event::getByToken(), mps_fire::i, metsig::jet, submitPVResolutionJobs::key, eostools::move(), submitPVResolutionJobs::out, edm::Event::put(), edm::refToPtr(), reco::GenJet::setSpecific(), and alcazmumu_cfi::src.

                                                                               {
  using namespace edm;
  using namespace std;

  Handle<View<reco::GenJet> > src;
  iEvent.getByToken(src_, src);

  auto out = std::make_unique<vector<reco::GenJet> >();
  out->reserve(src->size());

  Handle<edm::Association<std::vector<pat::PackedGenParticle> > > gp2pgp;
  iEvent.getByToken(gp2pgp_, gp2pgp);

  auto genJetSlimmedGenJetAssociation = make_unique<edm::Association<std::vector<reco::GenJet> > >();

  auto mapping = std::make_unique<std::vector<int> >();
  mapping->reserve(src->size());

  unsigned nm = 0;  // number of mapped jets
  unsigned nl = 0;  // number of loose jets
  for (View<reco::GenJet>::const_iterator it = src->begin(), ed = src->end(); it != ed; ++it) {
    bool selectedLoose = false;
    if (nLoose_ > 0 && nl < nLoose_ && cutLoose_(*it)) {
      selectedLoose = true;
      ++nl;
    }

    bool pass = cut_(*it) || selectedLoose;
    if (!pass) {
      mapping->push_back(-1);
      continue;
    }

    out->push_back(*it);
    reco::GenJet& jet = out->back();

    mapping->push_back(nm++);

    if (clearDaughters_) {
      jet.clearDaughters();
    } else  // rekey
    {
      //copy old
      reco::CompositePtrCandidate::daughters old = jet.daughterPtrVector();
      jet.clearDaughters();
      std::map<unsigned int, reco::CandidatePtr> ptrs;
      for (unsigned int i = 0; i < old.size(); i++) {
        //  if(! ((*gp2pgp)[old[i]]).isNonnull())   {
        //      std::cout << "Missing ref for key"  <<  old[i].key() << " pdgid " << old[i]->pdgId() << " st "<<   old[i]->status() <<  " pt " << old[i]->pt() << " eta " << old[i]->eta() << std::endl;
        //  }
        ptrs[((*gp2pgp)[old[i]]).key()] = refToPtr((*gp2pgp)[old[i]]);
      }
      for (std::map<unsigned int, reco::CandidatePtr>::iterator itp = ptrs.begin(); itp != ptrs.end();
           itp++)  //iterate on sorted items
      {
        jet.addDaughter(itp->second);
      }
    }
    if (dropSpecific_) {
      jet.setSpecific(reco::GenJet::Specific());
    }
  }

  edm::OrphanHandle<std::vector<reco::GenJet> > orphanHandle = iEvent.put(std::move(out));

  auto asso = std::make_unique<edm::Association<std::vector<reco::GenJet> > >(orphanHandle);
  edm::Association<std::vector<reco::GenJet> >::Filler slimmedAssoFiller(*asso);
  slimmedAssoFiller.insert(src, mapping->begin(), mapping->end());
  slimmedAssoFiller.fill();

  iEvent.put(std::move(asso), "slimmedGenJetAssociation");
}

Member Data Documentation

const bool pat::PATGenJetSlimmer::clearDaughters_

private

reset daughters to an empty vector

Definition at line 41 of file PATGenJetSlimmer.cc.

const StringCutObjectSelector<reco::GenJet> pat::PATGenJetSlimmer::cut_

private

Definition at line 36 of file PATGenJetSlimmer.cc.

const StringCutObjectSelector<reco::GenJet> pat::PATGenJetSlimmer::cutLoose_

private

Definition at line 37 of file PATGenJetSlimmer.cc.

const bool pat::PATGenJetSlimmer::dropSpecific_

private

drop the specific

Definition at line 43 of file PATGenJetSlimmer.cc.

const edm::EDGetTokenT<edm::Association<std::vector<pat::PackedGenParticle> > > pat::PATGenJetSlimmer::gp2pgp_

private

Definition at line 34 of file PATGenJetSlimmer.cc.

const unsigned pat::PATGenJetSlimmer::nLoose_

private

Definition at line 38 of file PATGenJetSlimmer.cc.

const edm::EDGetTokenT<edm::View<reco::GenJet> > pat::PATGenJetSlimmer::src_

private

Definition at line 33 of file PATGenJetSlimmer.cc.