Inheritance diagram for BJetEnergyRegressionMVA:

Public Member Functions
	BJetEnergyRegressionMVA (const edm::ParameterSet &iConfig, const BaseMVACache *cache)

void	fillAdditionalVariables (const pat::Jet &j) override

void	readAdditionalCollections (edm::Event &iEvent, const edm::EventSetup &) override
	to be implemented in derived classes, filling values for additional variables More...

Public Member Functions inherited from BaseMVAValueMapProducer< pat::Jet >
	BaseMVAValueMapProducer (const edm::ParameterSet &iConfig, const BaseMVACache *cache)

void	setValue (const std::string var, float val)

	~BaseMVAValueMapProducer () override

Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< BaseMVACache > >
	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

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

Static Public Member Functions inherited from BaseMVAValueMapProducer< pat::Jet >
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

static edm::ParameterSetDescription	getDescription ()

static void	globalEndJob (const BaseMVACache *cache)

static std::unique_ptr < BaseMVACache >	initializeGlobalCache (const edm::ParameterSet &cfg)

Private Attributes
edm::Handle< std::vector < reco::Vertex > >	pvs_

const edm::EDGetTokenT < std::vector< reco::Vertex > >	pvsrc_

edm::Handle< double >	rho_

edm::EDGetTokenT< double >	rhosrc_

edm::Handle< edm::View < reco::VertexCompositePtrCandidate > >	svs_

const edm::EDGetTokenT < edm::View < reco::VertexCompositePtrCandidate > >	svsrc_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< BaseMVACache > >
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 30 of file BJetEnergyRegressionMVA.cc.

Constructor & Destructor Documentation

BJetEnergyRegressionMVA::BJetEnergyRegressionMVA	(	const edm::ParameterSet &	iConfig,
		const BaseMVACache *	cache
	)

inlineexplicit

Definition at line 32 of file BJetEnergyRegressionMVA.cc.

       : BaseMVAValueMapProducer<pat::Jet>(iConfig, cache),
         pvsrc_(consumes<std::vector<reco::Vertex>>(iConfig.getParameter<edm::InputTag>("pvsrc"))),
         svsrc_(consumes<edm::View<reco::VertexCompositePtrCandidate>>(iConfig.getParameter<edm::InputTag>("svsrc"))),
         rhosrc_(consumes<double>(iConfig.getParameter<edm::InputTag>("rhosrc"))) {}

Member Function Documentation

void BJetEnergyRegressionMVA::fillAdditionalVariables ( const pat::Jet & j )

inlineoverridevirtual

Reimplemented from BaseMVAValueMapProducer< pat::Jet >.

Definition at line 44 of file BJetEnergyRegressionMVA.cc.

References funct::abs(), ztail::d, pat::Jet::daughterPtrVector(), HLT_FULL_cff::deltaR, pat::Jet::jecFactor(), cuda_std::lower_bound(), reco::LeafCandidate::p4(), edm::Handle< T >::product(), pvs_, rho_, and BaseMVAValueMapProducer< pat::Jet >::setValue().

                                                          {
     this->setValue("nPVs", pvs_->size());
     this->setValue("rho", *(rho_.product()));
 
     float cone_boundaries[] = {0.05, 0.1, 0.2, 0.3, 0.4};
     size_t ncone_boundaries = sizeof(cone_boundaries) / sizeof(float);
     std::vector<float> emFractionEnergyRings(ncone_boundaries + 1);
     std::vector<float> chFractionEnergyRings(ncone_boundaries + 1);
     std::vector<float> neFractionEnergyRings(ncone_boundaries + 1);
     std::vector<float> muFractionEnergyRings(ncone_boundaries + 1);
     float jetRawEnergy = j.p4().E() * j.jecFactor("Uncorrected");
     int numDaughtersPt03 = 0;
     for (unsigned int ijcone = 0; ijcone < ncone_boundaries; ijcone++) {
       emFractionEnergyRings[ijcone] = 0;
       muFractionEnergyRings[ijcone] = 0;
       chFractionEnergyRings[ijcone] = 0;
       neFractionEnergyRings[ijcone] = 0;
     }
     for (const auto& d : j.daughterPtrVector()) {
       float candDr = Geom::deltaR(d->p4(), j.p4());
       size_t icone =
           std::lower_bound(&cone_boundaries[0], &cone_boundaries[ncone_boundaries], candDr) - &cone_boundaries[0];
       float candEnergy = d->energy() / jetRawEnergy;
       int pdgid = abs(d->pdgId());
       if (pdgid == 22 || pdgid == 11) {
         emFractionEnergyRings[icone] += candEnergy;
       } else if (pdgid == 13) {
         muFractionEnergyRings[icone] += candEnergy;
       } else if (d->charge() != 0) {
         chFractionEnergyRings[icone] += candEnergy;
       } else {
         neFractionEnergyRings[icone] += candEnergy;
       }
       if (d->pt() > 0.3)
         numDaughtersPt03 += 1;
     }  // end of jet daughters loop
 
     this->setValue("Jet_energyRing_dR0_em_Jet_rawEnergy", emFractionEnergyRings[0]);
     this->setValue("Jet_energyRing_dR1_em_Jet_rawEnergy", emFractionEnergyRings[1]);
     this->setValue("Jet_energyRing_dR2_em_Jet_rawEnergy", emFractionEnergyRings[2]);
     this->setValue("Jet_energyRing_dR3_em_Jet_rawEnergy", emFractionEnergyRings[3]);
     this->setValue("Jet_energyRing_dR4_em_Jet_rawEnergy", emFractionEnergyRings[4]);
     //  this->setValue("Jet_energyRing_dR5_em_Jet_rawEnergy", emFractionEnergyRings[5]);
 
     this->setValue("Jet_energyRing_dR0_ch_Jet_rawEnergy", chFractionEnergyRings[0]);
     this->setValue("Jet_energyRing_dR1_ch_Jet_rawEnergy", chFractionEnergyRings[1]);
     this->setValue("Jet_energyRing_dR2_ch_Jet_rawEnergy", chFractionEnergyRings[2]);
     this->setValue("Jet_energyRing_dR3_ch_Jet_rawEnergy", chFractionEnergyRings[3]);
     this->setValue("Jet_energyRing_dR4_ch_Jet_rawEnergy", chFractionEnergyRings[4]);
     //  this->setValue("Jet_energyRing_dR5_ch_Jet_rawEnergy", chFractionEnergyRings[5]);
 
     this->setValue("Jet_energyRing_dR0_mu_Jet_rawEnergy", muFractionEnergyRings[0]);
     this->setValue("Jet_energyRing_dR1_mu_Jet_rawEnergy", muFractionEnergyRings[1]);
     this->setValue("Jet_energyRing_dR2_mu_Jet_rawEnergy", muFractionEnergyRings[2]);
     this->setValue("Jet_energyRing_dR3_mu_Jet_rawEnergy", muFractionEnergyRings[3]);
     this->setValue("Jet_energyRing_dR4_mu_Jet_rawEnergy", muFractionEnergyRings[4]);
     //  this->setValue("Jet_energyRing_dR5_mu_Jet_rawEnergy", muFractionEnergyRings[5]);
 
     this->setValue("Jet_energyRing_dR0_neut_Jet_rawEnergy", neFractionEnergyRings[0]);
     this->setValue("Jet_energyRing_dR1_neut_Jet_rawEnergy", neFractionEnergyRings[1]);
     this->setValue("Jet_energyRing_dR2_neut_Jet_rawEnergy", neFractionEnergyRings[2]);
     this->setValue("Jet_energyRing_dR3_neut_Jet_rawEnergy", neFractionEnergyRings[3]);
     this->setValue("Jet_energyRing_dR4_neut_Jet_rawEnergy", neFractionEnergyRings[4]);
     //    this->setValue("Jet_energyRing_dR5_neut_Jet_rawEnergy", neFractionEnergyRings[5]);
 
     this->setValue("Jet_numDaughters_pt03", numDaughtersPt03);
   }

static void BJetEnergyRegressionMVA::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

inlinestatic

Definition at line 112 of file BJetEnergyRegressionMVA.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), submitPVResolutionJobs::desc, and BaseMVAValueMapProducer< T >::getDescription().

                                                                            {
     edm::ParameterSetDescription desc = BaseMVAValueMapProducer<pat::Jet>::getDescription();
     desc.add<edm::InputTag>("pvsrc")->setComment("primary vertices input collection");
     desc.add<edm::InputTag>("svsrc")->setComment("secondary vertices input collection");
     desc.add<edm::InputTag>("rhosrc")->setComment("rho  input collection");
     descriptions.add("BJetEnergyRegressionMVA", desc);
   }

void BJetEnergyRegressionMVA::readAdditionalCollections	(	edm::Event &	,
		const edm::EventSetup &
	)

inlineoverridevirtual

to be implemented in derived classes, filling values for additional variables

Reimplemented from BaseMVAValueMapProducer< pat::Jet >.

Definition at line 38 of file BJetEnergyRegressionMVA.cc.

References edm::Event::getByToken(), pvs_, pvsrc_, rho_, rhosrc_, svs_, and svsrc_.

                                                                                   {
     iEvent.getByToken(pvsrc_, pvs_);
     iEvent.getByToken(svsrc_, svs_);
     iEvent.getByToken(rhosrc_, rho_);
   }