Inheritance diagram for BJetEnergyRegressionMVA:

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

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)

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

	~BaseMVAValueMapProducer () override

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

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

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 ()

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<>
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

Definition at line 30 of file BJetEnergyRegressionMVA.cc.

Constructor & Destructor Documentation

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

inlineexplicit

Definition at line 32 of file BJetEnergyRegressionMVA.cc.

       : BaseMVAValueMapProducer<pat::Jet>(iConfig),
         pvsrc_(edm::stream::EDProducer<>::consumes<std::vector<reco::Vertex>>(
             iConfig.getParameter<edm::InputTag>("pvsrc"))),
         svsrc_(edm::stream::EDProducer<>::consumes<edm::View<reco::VertexCompositePtrCandidate>>(
             iConfig.getParameter<edm::InputTag>("svsrc"))),
         rhosrc_(edm::stream::EDProducer<>::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 47 of file BJetEnergyRegressionMVA.cc.

References funct::abs(), ztail::d, pat::Jet::daughterPtrVector(), PbPb_ZMuSkimMuonDPG_cff::deltaR, pat::Jet::jecFactor(), reco::LeafCandidate::p4(), EgammaValidation_cff::pdgid, 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 115 of file BJetEnergyRegressionMVA.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), 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 41 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_);
   }