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
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndRuns () const final

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

Matcher of reconstructed objects to L1 Muons.

Author

Giovanni Petrucciani

Definition at line 25 of file PATGenJetSlimmer.cc.

Constructor & Destructor Documentation

PATGenJetSlimmer()

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

~PATGenJetSlimmer()

pat::PATGenJetSlimmer::~PATGenJetSlimmer ( )

inlineoverride

Definition at line 28 of file PATGenJetSlimmer.cc.

{}

Member Function Documentation

produce()

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

override

Definition at line 61 of file PATGenJetSlimmer.cc.

                                                                               {
  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 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(it - src->begin());
 
    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");
}

References mps_fire::i, iEvent, metsig::jet, crabWrapper::key, taus_updatedMVAIds_cff::mapping, eostools::move(), MillePedeFileConverter_cfg::out, edm::refToPtr(), and TrackRefitter_38T_cff::src.

Member Data Documentation

clearDaughters_

const bool pat::PATGenJetSlimmer::clearDaughters_

private

reset daughters to an empty vector

Definition at line 41 of file PATGenJetSlimmer.cc.

cut_

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

private

Definition at line 36 of file PATGenJetSlimmer.cc.

cutLoose_

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

private

Definition at line 37 of file PATGenJetSlimmer.cc.

dropSpecific_

const bool pat::PATGenJetSlimmer::dropSpecific_

private

drop the specific

Definition at line 43 of file PATGenJetSlimmer.cc.

gp2pgp_

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

private

Definition at line 34 of file PATGenJetSlimmer.cc.

nLoose_

const unsigned pat::PATGenJetSlimmer::nLoose_

private

Definition at line 38 of file PATGenJetSlimmer.cc.

src_

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

private

Definition at line 33 of file PATGenJetSlimmer.cc.