PileupJPTJetIdProducer Class Reference

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::PileupJPTJetIdProducer ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 36 of file PileupJPTJetIdProducer.cc.

References allowMissingInputs_, cms::PileupJPTJetIdAlgo::bookMVAReader(), 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::~PileupJPTJetIdProducer ( )

override

Definition at line 48 of file PileupJPTJetIdProducer.cc.

{}

Member Function Documentation

produce()

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, cms::PileupJPTJetIdAlgo::fillJPTBlock(), mps_fire::i, iEvent, input_token_, metsig::jet, PDWG_EXODelayedJetMET_cff::jets, eostools::move(), beam_dqm_sourceclient-live_cfg::mva, pualgo, and verbosity.

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

Member Data Documentation

allowMissingInputs_

bool PileupJPTJetIdProducer::allowMissingInputs_

private

Definition at line 30 of file PileupJPTJetIdProducer.cc.

Referenced by PileupJPTJetIdProducer().

input_token_

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

private

Definition at line 29 of file PileupJPTJetIdProducer.cc.

Referenced by PileupJPTJetIdProducer(), and produce().

jets_

edm::InputTag PileupJPTJetIdProducer::jets_

private

Definition at line 28 of file PileupJPTJetIdProducer.cc.

Referenced by PileupJPTJetIdProducer().

pualgo

cms::PileupJPTJetIdAlgo* PileupJPTJetIdProducer::pualgo

private

Definition at line 32 of file PileupJPTJetIdProducer.cc.

Referenced by PileupJPTJetIdProducer(), and produce().

verbosity

int PileupJPTJetIdProducer::verbosity

private

Definition at line 31 of file PileupJPTJetIdProducer.cc.

Referenced by PileupJPTJetIdProducer(), and produce().