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

Public Member Functions inherited from edm::stream::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

	~EDProducerBase () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

virtual	~ProducerBase () noexcept(false)

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

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 Public Member Functions inherited from edm::stream::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

Public Types inherited from edm::stream::EDProducerBase
typedef EDProducerAdaptorBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 31 of file BJetEnergyRegressionMVA.cc.

Constructor & Destructor Documentation

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

inlineexplicit

Definition at line 33 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 50 of file BJetEnergyRegressionMVA.cc.

References funct::abs(), edmIntegrityCheck::d, pat::Jet::daughterPtrVector(), deltaR(), pat::Jet::jecFactor(), reco::LeafCandidate::p4(), BPhysicsValidation_cfi::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 119 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 43 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_);
 
       }