Inheritance diagram for PileupJPTJetIdProducer:

Public Member Functions
	PileupJPTJetIdProducer (const edm::ParameterSet &)

	~PileupJPTJetIdProducer () 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 Member Functions
void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
bool	allowMissingInputs_

edm::EDGetTokenT< edm::View < reco::JPTJet > >	input_token_

edm::InputTag	jets_

cms::PileupJPTJetIdAlgo *	pualgo

int	verbosity

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 20 of file PileupJPTJetIdProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 36 of file PileupJPTJetIdProducer.cc.

References allowMissingInputs_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), input_token_, jets_, pualgo, and verbosity.

                                                                              {
   jets_ = iConfig.getParameter<edm::InputTag>("jets");
   input_token_ = consumes<edm::View<reco::JPTJet>>(jets_);
   verbosity = iConfig.getParameter<int>("Verbosity");
   allowMissingInputs_ = iConfig.getUntrackedParameter<bool>("AllowMissingInputs", false);
   pualgo = new cms::PileupJPTJetIdAlgo(iConfig);
   pualgo->bookMVAReader();
   produces<edm::ValueMap<float>>("JPTPUDiscriminant");
   produces<edm::ValueMap<int>>("JPTPUId");
 }

PileupJPTJetIdProducer::~PileupJPTJetIdProducer ( )

override

Definition at line 48 of file PileupJPTJetIdProducer.cc.

48 {}

Member Function Documentation

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

overrideprivate

|eta|<2.6

Definition at line 51 of file PileupJPTJetIdProducer.cc.

References b, gather_cfg::cout, reco::LeafCandidate::eta(), cms::PileupJPTJetIdAlgo::fillJPTBlock(), edm::Event::getByToken(), reco::JPTJet::getSpecific(), mps_fire::i, input_token_, metsig::jet, fwrapper::jets, eostools::move(), beam_dqm_sourceclient-live_cfg::mva, reco::LeafCandidate::phi(), reco::LeafCandidate::pt(), pualgo, edm::Event::put(), dt_dqm_sourceclient_common_cff::reco, verbosity, and reco::JPTJet::Specific::Zch.

                                                                                   {
   using namespace edm;
   using namespace std;
   using namespace reco;
   edm::Handle<View<JPTJet>> jets;
   iEvent.getByToken(input_token_, jets);
   vector<float> mva;
   vector<int> idflag;
   for (unsigned int i = 0; i < jets->size(); ++i) {
     int b = -1;
     const JPTJet& jet = jets->at(i);
 
     float mvapu = pualgo->fillJPTBlock(&jet);
 
     mva.push_back(mvapu);
 
     // Get PUid type
     //WP 95% JPT PUID > 0.3
     //WP 90% JPT PUID > 0.7
     //WP 80% JPT PUID > 0.9
 
     //|eta|>=2.6
     //WP 90% JPT PUID > -0.55
     //WP 80% JPT PUID > -0.3
     //WP 70% JPT PUID > -0.1
 
     if (fabs(jet.eta()) < 2.6) {
       if (mvapu > 0.3)
         b = 0;
       if (mvapu > 0.7)
         b = 1;
       if (mvapu > 0.9)
         b = 2;
     } else {
       if (mvapu > -0.55)
         b = 0;
       if (mvapu > -0.3)
         b = 1;
       if (mvapu > -0.1)
         b = 2;
     }
 
     idflag.push_back(b);
 
     if (verbosity > 0)
       std::cout << " PUID producer::Corrected JPT Jet is " << jet.pt() << " " << jet.eta() << " " << jet.phi() << " "
                 << jet.getSpecific().Zch << std::endl;
   }
 
   auto mvaout = std::make_unique<ValueMap<float>>();
   ValueMap<float>::Filler mvafiller(*mvaout);
   mvafiller.insert(jets, mva.begin(), mva.end());
   mvafiller.fill();
   iEvent.put(std::move(mvaout), "JPTPUDiscriminant");
 
   auto idflagout = std::make_unique<ValueMap<int>>();
   ValueMap<int>::Filler idflagfiller(*idflagout);
   idflagfiller.insert(jets, idflag.begin(), idflag.end());
   idflagfiller.fill();
   iEvent.put(std::move(idflagout), "JPTPUId");
 }