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	hasAbilityToProduceInLumis () const final
bool	hasAbilityToProduceInRuns () 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

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

explicit

Definition at line 49 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.

References iEvent, and produce().

{ }

Member Function Documentation

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

override

Definition at line 63 of file PATGenJetSlimmer.cc.

References reco::CompositePtrCandidate::addDaughter(), reco::CompositePtrCandidate::clearDaughters(), clearDaughters_, cut_, cutLoose_, reco::CompositePtrCandidate::daughterPtrVector(), DEFINE_FWK_MODULE, dropSpecific_, edm::Event::getByToken(), gp2pgp_, mps_fire::i, metsig::jet, crabWrapper::key, HPSPFTaus_cff::mapping, eostools::move(), nLoose_, MillePedeFileConverter_cfg::out, edm::Event::put(), edm::refToPtr(), reco::GenJet::setSpecific(), TrackRefitter_38T_cff::src, and src_.

Referenced by ~PATGenJetSlimmer().

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

Member Data Documentation

const bool pat::PATGenJetSlimmer::clearDaughters_

private

reset daughters to an empty vector

Definition at line 41 of file PATGenJetSlimmer.cc.

Referenced by produce().

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

private

Definition at line 36 of file PATGenJetSlimmer.cc.

Referenced by produce().

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

private

Definition at line 37 of file PATGenJetSlimmer.cc.

Referenced by produce().

const bool pat::PATGenJetSlimmer::dropSpecific_

private

drop the specific

Definition at line 43 of file PATGenJetSlimmer.cc.

Referenced by produce().

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

private

Definition at line 34 of file PATGenJetSlimmer.cc.

Referenced by produce().

const unsigned pat::PATGenJetSlimmer::nLoose_

private

Definition at line 38 of file PATGenJetSlimmer.cc.

Referenced by produce().

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

private

Definition at line 33 of file PATGenJetSlimmer.cc.

Referenced by produce().