#include "PhysicsTools/PatAlgos/interface/PATJetProducer.h"

Inheritance diagram for pat::PATJetProducer:

Public Member Functions
	PATJetProducer (const edm::ParameterSet &iConfig)

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

	~PATJetProducer () override

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Private Attributes
bool	addAssociatedTracks_

bool	addBTagInfo_

bool	addDiscriminators_

bool	addEfficiencies_

bool	addGenJetMatch_

bool	addGenPartonMatch_

bool	addJetCharge_

bool	addJetCorrFactors_

bool	addJetFlavourInfo_

bool	addJetID_

bool	addPartonJetMatch_

bool	addResolutions_

bool	addTagInfos_

GreaterByPt< CaloTower >	caloPTComparator_

std::vector< std::string >	discriminatorLabels_

std::vector< edm::InputTag >	discriminatorTags_

std::vector< edm::EDGetTokenT< reco::JetFloatAssociation::Container > >	discriminatorTokens_

pat::helper::EfficiencyLoader	efficiencyLoader_

bool	embedCaloTowers_

bool	embedGenJetMatch_

bool	embedGenPartonMatch_

bool	embedPFCandidates_

edm::EDGetTokenT< edm::Association< reco::GenJetCollection > >	genJetToken_

edm::EDGetTokenT< edm::Association< reco::GenParticleCollection > >	genPartonToken_

bool	getJetMCFlavour_

edm::EDGetTokenT< reco::JetFloatAssociation::Container >	jetChargeToken_

std::vector< edm::EDGetTokenT< edm::ValueMap< JetCorrFactors > > >	jetCorrFactorsTokens_

edm::EDGetTokenT< reco::JetFlavourInfoMatchingCollection >	jetFlavourInfoToken_

edm::EDGetTokenT< reco::JetIDValueMap >	jetIDMapToken_

edm::EDGetTokenT< reco::JetFlavourMatchingCollection >	jetPartonMapToken_

edm::EDGetTokenT< edm::View< reco::Jet > >	jetsToken_

bool	printWarning_

GreaterByPt< Jet >	pTComparator_

pat::helper::KinResolutionsLoader	resolutionLoader_

std::vector< std::string >	tagInfoLabels_

std::vector< edm::InputTag >	tagInfoTags_

std::vector< edm::EDGetTokenT< edm::View< reco::BaseTagInfo > > >	tagInfoTokens_

edm::EDGetTokenT< reco::JetTracksAssociation::Container >	trackAssociationToken_

bool	useLegacyJetMCFlavour_

pat::PATUserDataHelper< pat::Jet >	userDataHelper_

bool	useUserData_

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

Produces pat::Jet's.

The PATJetProducer produces analysis-level pat::Jet's starting from a collection of objects of JetType.

Author: Steven Lowette, Jeremy Andrea

Version

Id: PATJetProducer.h,v 1.26 2010/08/09 18:13:54 srappocc Exp

Definition at line 43 of file PATJetProducer.h.

Constructor & Destructor Documentation

◆ PATJetProducer()

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

explicit

Definition at line 40 of file PATJetProducer.cc.

     : useUserData_(iConfig.exists("userData")), printWarning_(true) {
   // initialize configurables
   jetsToken_ = consumes<edm::View<reco::Jet>>(iConfig.getParameter<edm::InputTag>("jetSource"));
   embedCaloTowers_ = false;  // parameter is optional
   if (iConfig.exists("embedCaloTowers")) {
     embedCaloTowers_ = iConfig.getParameter<bool>("embedCaloTowers");
   }
   embedPFCandidates_ = iConfig.getParameter<bool>("embedPFCandidates");
   getJetMCFlavour_ = iConfig.getParameter<bool>("getJetMCFlavour");
   useLegacyJetMCFlavour_ = iConfig.getParameter<bool>("useLegacyJetMCFlavour");
   addJetFlavourInfo_ = (useLegacyJetMCFlavour_ ? false : iConfig.getParameter<bool>("addJetFlavourInfo"));
   if (getJetMCFlavour_ && useLegacyJetMCFlavour_)
     jetPartonMapToken_ =
         consumes<reco::JetFlavourMatchingCollection>(iConfig.getParameter<edm::InputTag>("JetPartonMapSource"));
   else if (getJetMCFlavour_ && !useLegacyJetMCFlavour_)
     jetFlavourInfoToken_ =
         consumes<reco::JetFlavourInfoMatchingCollection>(iConfig.getParameter<edm::InputTag>("JetFlavourInfoSource"));
   addGenPartonMatch_ = iConfig.getParameter<bool>("addGenPartonMatch");
   embedGenPartonMatch_ = iConfig.getParameter<bool>("embedGenPartonMatch");
   if (addGenPartonMatch_)
     genPartonToken_ =
         consumes<edm::Association<reco::GenParticleCollection>>(iConfig.getParameter<edm::InputTag>("genPartonMatch"));
   addGenJetMatch_ = iConfig.getParameter<bool>("addGenJetMatch");
   embedGenJetMatch_ = iConfig.getParameter<bool>("embedGenJetMatch");
   if (addGenJetMatch_)
     genJetToken_ =
         consumes<edm::Association<reco::GenJetCollection>>(iConfig.getParameter<edm::InputTag>("genJetMatch"));
   addPartonJetMatch_ = iConfig.getParameter<bool>("addPartonJetMatch");
   //   partonJetToken_ = mayConsume<reco::SomePartonJetType>(iConfig.getParameter<edm::InputTag>( "partonJetSource" ));
   addJetCorrFactors_ = iConfig.getParameter<bool>("addJetCorrFactors");
   if (addJetCorrFactors_) {
     jetCorrFactorsTokens_ = edm::vector_transform(
         iConfig.getParameter<std::vector<edm::InputTag>>("jetCorrFactorsSource"),
         [this](edm::InputTag const &tag) { return consumes<edm::ValueMap<JetCorrFactors>>(tag); });
   }
   addBTagInfo_ = iConfig.getParameter<bool>("addBTagInfo");
   addDiscriminators_ = iConfig.getParameter<bool>("addDiscriminators");
   discriminatorTags_ = iConfig.getParameter<std::vector<edm::InputTag>>("discriminatorSources");
   discriminatorTokens_ = edm::vector_transform(discriminatorTags_, [this](edm::InputTag const &tag) {
     return mayConsume<reco::JetFloatAssociation::Container>(tag);
   });
   addTagInfos_ = iConfig.getParameter<bool>("addTagInfos");
   tagInfoTags_ = iConfig.getParameter<std::vector<edm::InputTag>>("tagInfoSources");
   tagInfoTokens_ = edm::vector_transform(
       tagInfoTags_, [this](edm::InputTag const &tag) { return mayConsume<edm::View<reco::BaseTagInfo>>(tag); });
   addAssociatedTracks_ = iConfig.getParameter<bool>("addAssociatedTracks");
   if (addAssociatedTracks_)
     trackAssociationToken_ =
         consumes<reco::JetTracksAssociation::Container>(iConfig.getParameter<edm::InputTag>("trackAssociationSource"));
   addJetCharge_ = iConfig.getParameter<bool>("addJetCharge");
   if (addJetCharge_)
     jetChargeToken_ =
         consumes<reco::JetFloatAssociation::Container>(iConfig.getParameter<edm::InputTag>("jetChargeSource"));
   addJetID_ = iConfig.getParameter<bool>("addJetID");
   if (addJetID_)
     jetIDMapToken_ = consumes<reco::JetIDValueMap>(iConfig.getParameter<edm::InputTag>("jetIDMap"));
   // Efficiency configurables
   addEfficiencies_ = iConfig.getParameter<bool>("addEfficiencies");
   if (addEfficiencies_) {
     efficiencyLoader_ =
         pat::helper::EfficiencyLoader(iConfig.getParameter<edm::ParameterSet>("efficiencies"), consumesCollector());
   }
   // Resolution configurables
   addResolutions_ = iConfig.getParameter<bool>("addResolutions");
   if (addResolutions_) {
     resolutionLoader_ = pat::helper::KinResolutionsLoader(iConfig.getParameter<edm::ParameterSet>("resolutions"));
   }
   if (discriminatorTags_.empty()) {
     addDiscriminators_ = false;
   } else {
     for (std::vector<edm::InputTag>::const_iterator it = discriminatorTags_.begin(), ed = discriminatorTags_.end();
          it != ed;
          ++it) {
       std::string label = it->label();
       std::string::size_type pos = label.find("JetTags");
       if ((pos != std::string::npos) && (pos != label.length() - 7)) {
         label.erase(pos + 7);  // trim a tail after "JetTags"
       }
       if (!it->instance().empty()) {
         label = (label + std::string(":") + it->instance());
       }
       discriminatorLabels_.push_back(label);
     }
   }
   if (tagInfoTags_.empty()) {
     addTagInfos_ = false;
   } else {
     for (std::vector<edm::InputTag>::const_iterator it = tagInfoTags_.begin(), ed = tagInfoTags_.end(); it != ed;
          ++it) {
       std::string label = it->label();
       std::string::size_type pos = label.find("TagInfos");
       if ((pos != std::string::npos) && (pos != label.length() - 8)) {
         label.erase(pos + 8);  // trim a tail after "TagInfos"
       }
       tagInfoLabels_.push_back(label);
     }
   }
   if (!addBTagInfo_) {
     addDiscriminators_ = false;
     addTagInfos_ = false;
   }
   // Check to see if the user wants to add user data
   if (useUserData_) {
     userDataHelper_ = PATUserDataHelper<Jet>(iConfig.getParameter<edm::ParameterSet>("userData"), consumesCollector());
   }
   // produces vector of jets
   produces<std::vector<Jet>>();
   produces<reco::GenJetCollection>("genJets");
   produces<std::vector<CaloTower>>("caloTowers");
   produces<reco::PFCandidateCollection>("pfCandidates");
   produces<edm::OwnVector<reco::BaseTagInfo>>("tagInfos");
 }

References addAssociatedTracks_, addBTagInfo_, addDiscriminators_, addEfficiencies_, addGenJetMatch_, addGenPartonMatch_, addJetCharge_, addJetCorrFactors_, addJetFlavourInfo_, addJetID_, addPartonJetMatch_, addResolutions_, addTagInfos_, discriminatorLabels_, discriminatorTags_, discriminatorTokens_, efficiencyLoader_, embedCaloTowers_, embedGenJetMatch_, embedGenPartonMatch_, embedPFCandidates_, edm::ParameterSet::exists(), funct::false, genJetToken_, genPartonToken_, getJetMCFlavour_, edm::ParameterSet::getParameter(), jetChargeToken_, jetCorrFactorsTokens_, jetFlavourInfoToken_, jetIDMapToken_, jetPartonMapToken_, jetsToken_, label, resolutionLoader_, AlCaHLTBitMon_QueryRunRegistry::string, GlobalPosition_Frontier_DevDB_cff::tag, tagInfoLabels_, tagInfoTags_, tagInfoTokens_, trackAssociationToken_, useLegacyJetMCFlavour_, userDataHelper_, useUserData_, and edm::vector_transform().

◆ ~PATJetProducer()

PATJetProducer::~PATJetProducer ( )

override

Definition at line 154 of file PATJetProducer.cc.

154 {}

Member Function Documentation

◆ fillDescriptions()

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

static

Definition at line 462 of file PATJetProducer.cc.

                                                                                 {
   edm::ParameterSetDescription iDesc;
   iDesc.setComment("PAT jet producer module");
  
   // input source
   iDesc.add<edm::InputTag>("jetSource", edm::InputTag("no default"))->setComment("input collection");
  
   // embedding
   iDesc.addOptional<bool>("embedCaloTowers", false)
       ->setComment("embed external CaloTowers (not to be used on AOD input)");
   iDesc.add<bool>("embedPFCandidates", true)->setComment("embed external PFCandidates");
  
   // MC matching configurables
   iDesc.add<bool>("addGenPartonMatch", true)->setComment("add MC matching");
   iDesc.add<bool>("embedGenPartonMatch", false)->setComment("embed MC matched MC information");
   iDesc.add<edm::InputTag>("genPartonMatch", edm::InputTag())->setComment("input with MC match information");
  
   iDesc.add<bool>("addGenJetMatch", true)->setComment("add MC matching");
   iDesc.add<bool>("embedGenJetMatch", false)->setComment("embed MC matched MC information");
   iDesc.add<edm::InputTag>("genJetMatch", edm::InputTag())->setComment("input with MC match information");
  
   iDesc.add<bool>("addJetCharge", true);
   iDesc.add<edm::InputTag>("jetChargeSource", edm::InputTag("patJetCharge"));
  
   // jet id
   iDesc.add<bool>("addJetID", true)->setComment("Add jet ID information");
   iDesc.add<edm::InputTag>("jetIDMap", edm::InputTag())->setComment("jet id map");
  
   iDesc.add<bool>("addPartonJetMatch", false);
   iDesc.add<edm::InputTag>("partonJetSource", edm::InputTag("NOT IMPLEMENTED"));
  
   // track association
   iDesc.add<bool>("addAssociatedTracks", true);
   iDesc.add<edm::InputTag>("trackAssociationSource", edm::InputTag("ak4JTA"));
  
   // tag info
   iDesc.add<bool>("addTagInfos", true);
   std::vector<edm::InputTag> emptyVInputTags;
   iDesc.add<std::vector<edm::InputTag>>("tagInfoSources", emptyVInputTags);
  
   // jet energy corrections
   iDesc.add<bool>("addJetCorrFactors", true);
   iDesc.add<std::vector<edm::InputTag>>("jetCorrFactorsSource", emptyVInputTags);
  
   // btag discriminator tags
   iDesc.add<bool>("addBTagInfo", true);
   iDesc.add<bool>("addDiscriminators", true);
   iDesc.add<std::vector<edm::InputTag>>("discriminatorSources", emptyVInputTags);
  
   // jet flavour idetification configurables
   iDesc.add<bool>("getJetMCFlavour", true);
   iDesc.add<bool>("useLegacyJetMCFlavour", false);
   iDesc.add<bool>("addJetFlavourInfo", false);
   iDesc.add<edm::InputTag>("JetPartonMapSource", edm::InputTag("jetFlavourAssociationLegacy"));
   iDesc.add<edm::InputTag>("JetFlavourInfoSource", edm::InputTag("jetFlavourAssociation"));
  
   pat::helper::KinResolutionsLoader::fillDescription(iDesc);
  
   // Efficiency configurables
   edm::ParameterSetDescription efficienciesPSet;
   efficienciesPSet.setAllowAnything();  // TODO: the pat helper needs to implement a description.
   iDesc.add("efficiencies", efficienciesPSet);
   iDesc.add<bool>("addEfficiencies", false);
  
   // Check to see if the user wants to add user data
   edm::ParameterSetDescription userDataPSet;
   PATUserDataHelper<Jet>::fillDescription(userDataPSet);
   iDesc.addOptional("userData", userDataPSet);
  
   descriptions.add("PATJetProducer", iDesc);
 }

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), edm::ParameterSetDescription::addOptional(), pat::helper::KinResolutionsLoader::fillDescription(), pat::PATUserDataHelper< ObjectType >::fillDescription(), HLT_FULL_cff::InputTag, edm::ParameterSetDescription::setAllowAnything(), and edm::ParameterSetDescription::setComment().

◆ produce()

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

override

Definition at line 156 of file PATJetProducer.cc.

                                                                           {
   // check whether dealing with MC or real data
   if (iEvent.isRealData()) {
     getJetMCFlavour_ = false;
     useLegacyJetMCFlavour_ = false;
     addJetFlavourInfo_ = false;
     addGenPartonMatch_ = false;
     addGenJetMatch_ = false;
     addPartonJetMatch_ = false;
   }
  
   // Get the vector of jets
   edm::Handle<edm::View<reco::Jet>> jets;
   iEvent.getByToken(jetsToken_, jets);
  
   if (efficiencyLoader_.enabled())
     efficiencyLoader_.newEvent(iEvent);
   if (resolutionLoader_.enabled())
     resolutionLoader_.newEvent(iEvent, iSetup);
  
   // for jet flavour
   edm::Handle<reco::JetFlavourMatchingCollection> jetFlavMatch;
   edm::Handle<reco::JetFlavourInfoMatchingCollection> jetFlavInfoMatch;
   if (getJetMCFlavour_ && useLegacyJetMCFlavour_)
     iEvent.getByToken(jetPartonMapToken_, jetFlavMatch);
   else if (getJetMCFlavour_ && !useLegacyJetMCFlavour_)
     iEvent.getByToken(jetFlavourInfoToken_, jetFlavInfoMatch);
  
   // Get the vector of generated particles from the event if needed
   edm::Handle<edm::Association<reco::GenParticleCollection>> partonMatch;
   if (addGenPartonMatch_)
     iEvent.getByToken(genPartonToken_, partonMatch);
   // Get the vector of GenJets from the event if needed
   edm::Handle<edm::Association<reco::GenJetCollection>> genJetMatch;
   if (addGenJetMatch_)
     iEvent.getByToken(genJetToken_, genJetMatch);
   /* TO BE IMPLEMENTED FOR >= 1_5_X
   // Get the vector of PartonJets from the event if needed
   edm::Handle<edm::View<reco::SomePartonJetType> > partonJets;
   if (addPartonJetMatch_) iEvent.getByToken(partonJetToken_, partonJets);
 */
  
   // read in the jet correction factors ValueMap
   std::vector<edm::ValueMap<JetCorrFactors>> jetCorrs;
   if (addJetCorrFactors_) {
     for (size_t i = 0; i < jetCorrFactorsTokens_.size(); ++i) {
       edm::Handle<edm::ValueMap<JetCorrFactors>> jetCorr;
       iEvent.getByToken(jetCorrFactorsTokens_[i], jetCorr);
       jetCorrs.push_back(*jetCorr);
     }
   }
  
   // Get the vector of jet tags with b-tagging info
   std::vector<edm::Handle<reco::JetFloatAssociation::Container>> jetDiscriminators;
   if (addBTagInfo_ && addDiscriminators_) {
     jetDiscriminators.resize(discriminatorTokens_.size());
     for (size_t i = 0; i < discriminatorTokens_.size(); ++i) {
       iEvent.getByToken(discriminatorTokens_[i], jetDiscriminators[i]);
     }
   }
   std::vector<edm::Handle<edm::View<reco::BaseTagInfo>>> jetTagInfos;
   if (addBTagInfo_ && addTagInfos_) {
     jetTagInfos.resize(tagInfoTokens_.size());
     for (size_t i = 0; i < tagInfoTokens_.size(); ++i) {
       iEvent.getByToken(tagInfoTokens_[i], jetTagInfos[i]);
     }
   }
  
   // tracks Jet Track Association
   edm::Handle<reco::JetTracksAssociation::Container> hTrackAss;
   if (addAssociatedTracks_)
     iEvent.getByToken(trackAssociationToken_, hTrackAss);
   edm::Handle<reco::JetFloatAssociation::Container> hJetChargeAss;
   if (addJetCharge_)
     iEvent.getByToken(jetChargeToken_, hJetChargeAss);
  
   // jet ID handle
   edm::Handle<reco::JetIDValueMap> hJetIDMap;
   if (addJetID_)
     iEvent.getByToken(jetIDMapToken_, hJetIDMap);
  
   // loop over jets
   auto patJets = std::make_unique<std::vector<Jet>>();
  
   auto genJetsOut = std::make_unique<reco::GenJetCollection>();
   auto caloTowersOut = std::make_unique<std::vector<CaloTower>>();
   auto pfCandidatesOut = std::make_unique<reco::PFCandidateCollection>();
   auto tagInfosOut = std::make_unique<edm::OwnVector<reco::BaseTagInfo>>();
  
   edm::RefProd<reco::GenJetCollection> h_genJetsOut = iEvent.getRefBeforePut<reco::GenJetCollection>("genJets");
   edm::RefProd<std::vector<CaloTower>> h_caloTowersOut = iEvent.getRefBeforePut<std::vector<CaloTower>>("caloTowers");
   edm::RefProd<reco::PFCandidateCollection> h_pfCandidatesOut =
       iEvent.getRefBeforePut<reco::PFCandidateCollection>("pfCandidates");
   edm::RefProd<edm::OwnVector<reco::BaseTagInfo>> h_tagInfosOut =
       iEvent.getRefBeforePut<edm::OwnVector<reco::BaseTagInfo>>("tagInfos");
  
   for (edm::View<reco::Jet>::const_iterator itJet = jets->begin(); itJet != jets->end(); itJet++) {
     // construct the Jet from the ref -> save ref to original object
     unsigned int idx = itJet - jets->begin();
     edm::RefToBase<reco::Jet> jetRef = jets->refAt(idx);
     edm::Ptr<reco::Jet> jetPtr = jets->ptrAt(idx);
     Jet ajet(jetRef);
  
     // add the FwdPtrs to the CaloTowers
     if ((ajet.isCaloJet() || ajet.isJPTJet()) && embedCaloTowers_) {
       const reco::CaloJet *cj = nullptr;
       const reco::JPTJet *jptj = nullptr;
       if (ajet.isCaloJet())
         cj = dynamic_cast<const reco::CaloJet *>(jetRef.get());
       else {
         jptj = dynamic_cast<const reco::JPTJet *>(jetRef.get());
         cj = dynamic_cast<const reco::CaloJet *>(jptj->getCaloJetRef().get());
       }
       pat::CaloTowerFwdPtrCollection itowersRef;
       std::vector<CaloTowerPtr> itowers = cj->getCaloConstituents();
       for (std::vector<CaloTowerPtr>::const_iterator towBegin = itowers.begin(),
                                                      towEnd = itowers.end(),
                                                      itow = towBegin;
            itow != towEnd;
            ++itow) {
         if (itow->isAvailable() && itow->isNonnull()) {
           caloTowersOut->push_back(**itow);
           // set the "forward" ref to the thinned collection
           edm::Ref<std::vector<CaloTower>> caloTowerRef(h_caloTowersOut, caloTowersOut->size() - 1);
           edm::Ptr<CaloTower> caloForwardRef(h_caloTowersOut.id(), caloTowerRef.key(), h_caloTowersOut.productGetter());
           // set the "backward" ref to the original collection for association
           edm::Ptr<CaloTower> caloBackRef(*itow);
           // add to the list of FwdPtr's
           itowersRef.push_back(pat::CaloTowerFwdPtrCollection::value_type(caloForwardRef, caloBackRef));
         }
       }
       ajet.setCaloTowers(itowersRef);
     }
  
     // add the FwdPtrs to the PFCandidates
     if (ajet.isPFJet() && embedPFCandidates_) {
       const reco::PFJet *cj = dynamic_cast<const reco::PFJet *>(jetRef.get());
       pat::PFCandidateFwdPtrCollection iparticlesRef;
       std::vector<reco::PFCandidatePtr> iparticles = cj->getPFConstituents();
       for (std::vector<reco::PFCandidatePtr>::const_iterator partBegin = iparticles.begin(),
                                                              partEnd = iparticles.end(),
                                                              ipart = partBegin;
            ipart != partEnd;
            ++ipart) {
         pfCandidatesOut->push_back(**ipart);
         // set the "forward" ref to the thinned collection
         edm::Ref<reco::PFCandidateCollection> pfCollectionRef(h_pfCandidatesOut, pfCandidatesOut->size() - 1);
         edm::Ptr<reco::PFCandidate> pfForwardRef(
             h_pfCandidatesOut.id(), pfCollectionRef.key(), h_pfCandidatesOut.productGetter());
         // set the "backward" ref to the original collection for association
         edm::Ptr<reco::PFCandidate> pfBackRef(*ipart);
         // add to the list of FwdPtr's
         iparticlesRef.push_back(pat::PFCandidateFwdPtrCollection::value_type(pfForwardRef, pfBackRef));
       }
       ajet.setPFCandidates(iparticlesRef);
     }
  
     if (addJetCorrFactors_) {
       // add additional JetCorrs to the jet
       for (unsigned int i = 0; i < jetCorrFactorsTokens_.size(); ++i) {
         const JetCorrFactors &jcf = jetCorrs[i][jetRef];
         // uncomment for debugging
         // jcf.print();
         ajet.addJECFactors(jcf);
       }
       std::vector<std::string> levels = jetCorrs[0][jetRef].correctionLabels();
       if (std::find(levels.begin(), levels.end(), "L2L3Residual") != levels.end()) {
         ajet.initializeJEC(jetCorrs[0][jetRef].jecLevel("L2L3Residual"));
       } else if (std::find(levels.begin(), levels.end(), "L3Absolute") != levels.end()) {
         ajet.initializeJEC(jetCorrs[0][jetRef].jecLevel("L3Absolute"));
       } else {
         ajet.initializeJEC(jetCorrs[0][jetRef].jecLevel("Uncorrected"));
         if (printWarning_) {
           edm::LogWarning("L3Absolute not found")
               << "L2L3Residual and L3Absolute are not part of the jetCorrFactors\n"
               << "of module " << jetCorrs[0][jetRef].jecSet() << ". Jets will remain"
               << " uncorrected.";
           printWarning_ = false;
         }
       }
     }
  
     // get the MC flavour information for this jet
     if (getJetMCFlavour_ && useLegacyJetMCFlavour_) {
       ajet.setPartonFlavour((*jetFlavMatch)[edm::RefToBase<reco::Jet>(jetRef)].getFlavour());
     } else if (getJetMCFlavour_ && !useLegacyJetMCFlavour_) {
       if (addJetFlavourInfo_)
         ajet.setJetFlavourInfo((*jetFlavInfoMatch)[edm::RefToBase<reco::Jet>(jetRef)]);
       else {
         ajet.setPartonFlavour((*jetFlavInfoMatch)[edm::RefToBase<reco::Jet>(jetRef)].getPartonFlavour());
         ajet.setHadronFlavour((*jetFlavInfoMatch)[edm::RefToBase<reco::Jet>(jetRef)].getHadronFlavour());
       }
     }
     // store the match to the generated partons
     if (addGenPartonMatch_) {
       reco::GenParticleRef parton = (*partonMatch)[jetRef];
       if (parton.isNonnull() && parton.isAvailable()) {
         ajet.setGenParton(parton, embedGenPartonMatch_);
       }  // leave empty if no match found
     }
     // store the match to the GenJets
     if (addGenJetMatch_) {
       reco::GenJetRef genjet = (*genJetMatch)[jetRef];
       if (genjet.isNonnull() && genjet.isAvailable()) {
         genJetsOut->push_back(*genjet);
         // set the "forward" ref to the thinned collection
         edm::Ref<reco::GenJetCollection> genForwardRef(h_genJetsOut, genJetsOut->size() - 1);
         // set the "backward" ref to the original collection
         const edm::Ref<reco::GenJetCollection> &genBackRef(genjet);
         // make the FwdPtr
         edm::FwdRef<reco::GenJetCollection> genjetFwdRef(genForwardRef, genBackRef);
         ajet.setGenJetRef(genjetFwdRef);
       }  // leave empty if no match found
     }
  
     if (efficiencyLoader_.enabled()) {
       efficiencyLoader_.setEfficiencies(ajet, jetRef);
     }
  
     // IMPORTANT: DO THIS AFTER JES CORRECTIONS
     if (resolutionLoader_.enabled()) {
       resolutionLoader_.setResolutions(ajet);
     }
  
     // TO BE IMPLEMENTED FOR >=1_5_X: do the PartonJet matching
     if (addPartonJetMatch_) {
     }
  
     // add b-tag info if available & required
     if (addBTagInfo_) {
       if (addDiscriminators_) {
         for (size_t k = 0; k < jetDiscriminators.size(); ++k) {
           float value = (*jetDiscriminators[k])[jetRef];
           ajet.addBDiscriminatorPair(std::make_pair(discriminatorLabels_[k], value));
         }
       }
       if (addTagInfos_) {
         for (size_t k = 0; k < jetTagInfos.size(); ++k) {
           const edm::View<reco::BaseTagInfo> &taginfos = *jetTagInfos[k];
           // This is not associative, so we have to search the jet
           edm::Ptr<reco::BaseTagInfo> match;
           // Try first by 'same index'
           if ((idx < taginfos.size()) && (taginfos[idx].jet() == jetRef)) {
             match = taginfos.ptrAt(idx);
           } else {
             // otherwise fail back to a simple search
             for (edm::View<reco::BaseTagInfo>::const_iterator itTI = taginfos.begin(), edTI = taginfos.end();
                  itTI != edTI;
                  ++itTI) {
               if (itTI->jet() == jetRef) {
                 match = taginfos.ptrAt(itTI - taginfos.begin());
                 break;
               }
             }
           }
           if (match.isNonnull()) {
             tagInfosOut->push_back(match->clone());
             // set the "forward" ptr to the thinned collection
             edm::Ptr<reco::BaseTagInfo> tagInfoForwardPtr(
                 h_tagInfosOut.id(), &tagInfosOut->back(), tagInfosOut->size() - 1);
             // set the "backward" ptr to the original collection for association
             const edm::Ptr<reco::BaseTagInfo> &tagInfoBackPtr(match);
             // make FwdPtr
             TagInfoFwdPtrCollection::value_type tagInfoFwdPtr(tagInfoForwardPtr, tagInfoBackPtr);
             ajet.addTagInfo(tagInfoLabels_[k], tagInfoFwdPtr);
           }
         }
       }
     }
  
     if (addAssociatedTracks_)
       ajet.setAssociatedTracks((*hTrackAss)[jetRef]);
  
     if (addJetCharge_)
       ajet.setJetCharge((*hJetChargeAss)[jetRef]);
  
     // add jet ID for calo jets
     if (addJetID_ && ajet.isCaloJet()) {
       reco::JetID jetId = (*hJetIDMap)[jetRef];
       ajet.setJetID(jetId);
     }
     // add jet ID jpt jets
     else if (addJetID_ && ajet.isJPTJet()) {
       const reco::JPTJet *jptj = dynamic_cast<const reco::JPTJet *>(jetRef.get());
       reco::JetID jetId = (*hJetIDMap)[jptj->getCaloJetRef()];
       ajet.setJetID(jetId);
     }
     if (useUserData_) {
       userDataHelper_.add(ajet, iEvent, iSetup);
     }
     patJets->push_back(ajet);
   }
  
   // sort jets in pt
   std::sort(patJets->begin(), patJets->end(), pTComparator_);
  
   // put genEvt  in Event
   iEvent.put(std::move(patJets));
  
   iEvent.put(std::move(genJetsOut), "genJets");
   iEvent.put(std::move(caloTowersOut), "caloTowers");
   iEvent.put(std::move(pfCandidatesOut), "pfCandidates");
   iEvent.put(std::move(tagInfosOut), "tagInfos");
 }