Produces a ValueMap between JetCorrFactors and the to the originating reco jets. More...

#include "PhysicsTools/PatAlgos/interface/JetCorrFactorsProducer.h"

Inheritance diagram for pat::JetCorrFactorsProducer:

Public Types
typedef std::map < JetCorrFactors::Flavor, std::vector< std::string > >	FlavorCorrLevelMap
	map of correction levels to different flavors More...

typedef edm::ValueMap < pat::JetCorrFactors >	JetCorrFactorsMap
	value map for JetCorrFactors (to be written into the event) More...

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
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Public Member Functions
	JetCorrFactorsProducer (const edm::ParameterSet &cfg)
	default constructor More...

virtual void	produce (edm::Event &event, const edm::EventSetup &setup) override
	everything that needs to be done per event More...

	~JetCorrFactorsProducer ()
	default destructor More...

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

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

ModuleDescription const &	moduleDescription () const

virtual	~EDProducerBase ()

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

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

virtual	~ProducerBase ()

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

	EDConsumerBase ()

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

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

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

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

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

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)
	description of configuration file parameters More...

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 Member Functions
float	evaluate (edm::View< reco::Jet >::const_iterator &jet, const JetCorrFactors::Flavor &flavor, int level)
	evaluate jet correction factor up to a given level More...

std::vector< std::string >	expand (const std::vector< std::string > &levels, const JetCorrFactors::Flavor &flavor)

bool	flavorDependent () const
	return true if the jec levels contain at least one flavor dependent correction level More...

int	numberOf (const edm::Handle< std::vector< reco::Vertex > > &primaryVertices)
	determines the number of valid primary vertices for the standard L1Offset correction of JetMET More...

std::vector < JetCorrectorParameters >	params (const JetCorrectorParametersCollection &parameters, const std::vector< std::string > &levels) const
	return the jec parameters as input to the FactorizedJetCorrector for different flavors More...

std::string	payload ()
	map jet algorithm to payload in DB More...

Private Attributes
unsigned long long	cacheId_
	cache identifier for JetCorrectionsRecord More...

std::map < JetCorrFactors::Flavor, std::unique_ptr < FactorizedJetCorrector > >	correctors_
	cache container for jet corrections More...

bool	emf_
	use electromagnetic fraction for jet energy corrections or not (will only have an effect for jets CaloJets) More...

std::vector< std::string >	extraJPTOffset_

std::unique_ptr < FactorizedJetCorrector >	extraJPTOffsetCorrector_
	cache container for JPTOffset jet corrections More...

std::string	label_
	label of jec factors More...

FlavorCorrLevelMap	levels_

std::string	payload_
	label of payload More...

edm::InputTag	primaryVertices_
	label for L1Offset primaryVertex collection More...

edm::EDGetTokenT< std::vector < reco::Vertex > >	primaryVerticesToken_

edm::InputTag	rho_
	label for L1FastJet energy density parameter rho More...

edm::EDGetTokenT< double >	rhoToken_

edm::EDGetTokenT< edm::View < reco::Jet > >	srcToken_
	input jet collection More...

std::string	type_
	type of flavor dependent JEC factors (only 'J' and 'T' are allowed) More...

bool	useNPV_
	use the NPV and rho with the JEC? (used for L1Offset/L1FastJet and L1FastJet, resp.) More...

bool	useRho_

Additional Inherited Members
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

Produces a ValueMap between JetCorrFactors and the to the originating reco jets.

The JetCorrFactorsProducer produces a set of correction factors, defined in the class pat::JetCorrFactors. This vector is linked to the originating reco jets through an edm::ValueMap. The initializing parameters of the module can be found in the recoLayer1/jetCorrFactors_cfi.py of the PatAlgos package. In the standard PAT workflow the module has to be run before the creation of the pat::Jet. The edm::ValueMap will then be embedded into the pat::Jet.

Jets corrected up to a given correction level can then be accessed via the pat::Jet member function correctedJet. For more details have a look into the class description of the pat::Jet.

ATTENTION: available options for flavor corrections are L5Flavor_gJ L7Parton_gJ gluon from dijets L5Flavor_qJ/_qT L7Parton_qJ/_qT quark from dijets/top L5Flavor_cJ/_cT L7Parton_cJ/_cT charm from dijets/top L5Flavor_bJ/_bT L7Parton_bJ/_bT beauty from dijets/top L7Parton_jJ/_tT mixture from dijets/top

where mixture refers to the flavor mixture as determined from the MC sample the flavor dependent corrections have been derived from. 'J' and 'T' stand for a typical dijet (ttbar) sample.

L1Offset corrections require the collection of offlinePrimaryVertices, which are supposed to be added as an additional optional parameter primaryVertices in the jetCorrFactors_cfi.py file.

L1FastJet corrections, which are an alternative to the standard L1Offset correction as recommended by the JetMET PAG the energy density parameter rho is supposed to be added as an additional optional parameter rho in the jetCorrFactors_cfi.py file.

NOTE: the mixed mode (mc input mixture from dijets/ttbar) only exists for parton level corrections. jJ and tT are not covered in this implementation of the JetCorrFactorsProducer there are no gluon corrections available from the top sample on the L7Parton level.

Definition at line 59 of file JetCorrFactorsProducer.h.

Member Typedef Documentation

typedef std::map<JetCorrFactors::Flavor, std::vector<std::string> > pat::JetCorrFactorsProducer::FlavorCorrLevelMap

map of correction levels to different flavors

Definition at line 64 of file JetCorrFactorsProducer.h.

typedef edm::ValueMap<pat::JetCorrFactors> pat::JetCorrFactorsProducer::JetCorrFactorsMap

value map for JetCorrFactors (to be written into the event)

Definition at line 62 of file JetCorrFactorsProducer.h.

Constructor & Destructor Documentation

JetCorrFactorsProducer::JetCorrFactorsProducer ( const edm::ParameterSet & cfg )

explicit

default constructor

Definition at line 22 of file JetCorrFactorsProducer.cc.

References pat::JetCorrFactors::BOTTOM, pat::JetCorrFactors::CHARM, Exception, edm::ParameterSet::existsAs(), expand(), extraJPTOffset_, spr::find(), edm::ParameterSet::getParameter(), pat::JetCorrFactors::GLUON, levels_, python.rootplot.argparse::message, pat::JetCorrFactors::NONE, primaryVertices_, primaryVerticesToken_, rho_, rhoToken_, AlCaHLTBitMon_QueryRunRegistry::string, pat::JetCorrFactors::UDS, useNPV_, and useRho_.

                                                                         :
   emf_(cfg.getParameter<bool>( "emf" )),
   srcToken_(consumes<edm::View<reco::Jet> >(cfg.getParameter<edm::InputTag>( "src" ))),
   type_ (cfg.getParameter<std::string>("flavorType")),
   label_(cfg.getParameter<std::string>( "@module_label" )),
   payload_( cfg.getParameter<std::string>("payload") ),
   useNPV_(cfg.getParameter<bool>("useNPV")),
   useRho_(cfg.getParameter<bool>("useRho")),
   cacheId_(0)
 {
   std::vector<std::string> levels = cfg.getParameter<std::vector<std::string> >("levels");
   // fill the std::map for levels_, which might be flavor dependent or not;
   // flavor dependency is determined from the fact whether the std::string
   // L5Flavor or L7Parton can be found in levels; if flavor dependent four
   // vectors of strings will be filled into the map corresponding to GLUON,
   // UDS, CHARM and BOTTOM (according to JetCorrFactors::Flavor), 'L5Flavor'
   // and 'L7Parton' will be expanded accordingly; if not levels_ is filled
   // with only one vector of strings according to NONE. This vector will be
   // equivalent to the original vector of strings.
   if(std::find(levels.begin(), levels.end(), "L5Flavor")!=levels.end() || std::find(levels.begin(), levels.end(), "L7Parton")!=levels.end()){
     levels_[JetCorrFactors::GLUON ] = expand(levels, JetCorrFactors::GLUON );
     levels_[JetCorrFactors::UDS   ] = expand(levels, JetCorrFactors::UDS   );
     levels_[JetCorrFactors::CHARM ] = expand(levels, JetCorrFactors::CHARM );
     levels_[JetCorrFactors::BOTTOM] = expand(levels, JetCorrFactors::BOTTOM);
   }
   else{
     levels_[JetCorrFactors::NONE  ] = levels;
   }
   // if the std::string L1JPTOffset can be found in levels an additional
   // parameter extraJPTOffset is needed, which should pass on the the usual
   // L1Offset correction, which is an additional input to the L1JPTOffset
   // corrector
   if(std::find(levels.begin(), levels.end(), "L1JPTOffset")!=levels.end()){
     if(cfg.existsAs<std::string>("extraJPTOffset")){
       extraJPTOffset_.push_back(cfg.getParameter<std::string>("extraJPTOffset"));
     }
     else{
       throw cms::Exception("No parameter extraJPTOffset specified")
     << "The configured correction levels contain a L1JPTOffset correction, which re- \n"
     << "quires the additional parameter extraJPTOffset or type std::string. This     \n"
     << "string should correspond to the L1Offset corrections that should be applied  \n"
     << "together with the JPTL1Offset corrections. These corrections can be of type  \n"
     << "L1Offset or L1FastJet.                                                       \n";
     }
   }
   // if the std::string L1Offset can be found in levels an additional para-
   // meter primaryVertices is needed, which should pass on the offline pri-
   // mary vertex collection. The size of this collection is needed for the
   // L1Offset correction.
   if(useNPV_){
     if(cfg.existsAs<edm::InputTag>("primaryVertices")){
       primaryVertices_=cfg.getParameter<edm::InputTag>("primaryVertices");
       primaryVerticesToken_=mayConsume<std::vector<reco::Vertex> >(primaryVertices_);
     }
     else{
       throw cms::Exception("No primaryVertices specified")
     << "The configured correction levels contain an L1Offset or L1FastJet correction, \n"
     << "which requires the number of offlinePrimaryVertices. Please specify this col- \n"
     << "lection as additional optional parameter primaryVertices of type edm::InputTag\n"
     << "in the jetCorrFactors module.                                                 \n";
     }
   }
   // if the std::string L1FastJet can be found in levels an additional
   // parameter rho is needed, which should pass on the energy density
   // parameter for the corresponding jet collection.
   if(useRho_){
     if((!extraJPTOffset_.empty() && extraJPTOffset_.front()==std::string("L1FastJet")) || std::find(levels.begin(), levels.end(), "L1FastJet")!=levels.end()){
       if(cfg.existsAs<edm::InputTag>("rho")){
     rho_=cfg.getParameter<edm::InputTag>("rho");
     rhoToken_=mayConsume<double>(rho_);
       }
       else{
     throw cms::Exception("No parameter rho specified")
       << "The configured correction levels contain a L1FastJet correction, which re- \n"
       << "quires the energy density parameter rho. Please specify this collection as \n"
       << "additional optional parameter rho of type edm::InputTag in the jetCorrFac- \n"
       << "tors module.                                                               \n";
       }
     }
     else{
       edm::LogInfo message( "Parameter rho not used" );
       message << "Module is configured to use the parameter rho, but rho is only used     \n"
           << "for L1FastJet corrections. The configuration of levels does not contain \n"
           << "L1FastJet corrections though, so rho will not be used by this module.   \n";
     }
   }
   produces<JetCorrFactorsMap>();
 }

pat::JetCorrFactorsProducer::~JetCorrFactorsProducer ( )

inline

default destructor

Definition at line 70 of file JetCorrFactorsProducer.h.

70 {};

Member Function Documentation

float JetCorrFactorsProducer::evaluate	(	edm::View< reco::Jet >::const_iterator &	jet,
		const JetCorrFactors::Flavor &	flavor,
		int	level
	)

private

evaluate jet correction factor up to a given level

Definition at line 149 of file JetCorrFactorsProducer.cc.

References pat::Jet::correctedP4(), mvaPFMET_cff::corrector, correctors_, emf_, reco::LeafCandidate::energy(), reco::LeafCandidate::eta(), extraJPTOffsetCorrector_, reco::JPTJet::getCaloJetRef(), metsig::jet, testEve_cfg::level, p4, reco::LeafCandidate::phi(), and reco::LeafCandidate::pt().

Referenced by produce().

 {
   std::unique_ptr<FactorizedJetCorrector>& corrector = correctors_.find(flavor)->second;
   // add parameters for JPT corrections
   const reco::JPTJet* jpt = dynamic_cast<reco::JPTJet const *>( &*jet );
   if( jpt ){
     TLorentzVector p4; p4.SetPtEtaPhiE(jpt->getCaloJetRef()->pt(), jpt->getCaloJetRef()->eta(), jpt->getCaloJetRef()->phi(), jpt->getCaloJetRef()->energy());
     if( extraJPTOffsetCorrector_ ){
       extraJPTOffsetCorrector_->setJPTrawP4(p4);
       corrector->setJPTrawOff(extraJPTOffsetCorrector_->getSubCorrections()[0]);
     }
     corrector->setJPTrawP4(p4);
   }
   //For PAT jets undo previous jet energy corrections
   const Jet* patjet = dynamic_cast<Jet const *>( &*jet );
   if( patjet ){
     corrector->setJetEta(patjet->correctedP4(0).eta()); corrector->setJetPt(patjet->correctedP4(0).pt()); corrector->setJetE(patjet->correctedP4(0).energy());
   } else {
     corrector->setJetEta(jet->eta()); corrector->setJetPt(jet->pt()); corrector->setJetE(jet->energy());
   }
   if( emf_ && dynamic_cast<const reco::CaloJet*>(&*jet)){
     corrector->setJetEMF(dynamic_cast<const reco::CaloJet*>(&*jet)->emEnergyFraction());
   }
   return corrector->getSubCorrections()[level];
 }

std::vector< std::string > JetCorrFactorsProducer::expand	(	const std::vector< std::string > &	levels,
		const JetCorrFactors::Flavor &	flavor
	)

private

return an expanded version of correction levels for different flavors; the result should be of type ['L2Relative', 'L3Absolute', 'L5FLavor_gJ', 'L7Parton_gJ']; L7Parton_gT will result in an empty string as this correction level is not available

Definition at line 112 of file JetCorrFactorsProducer.cc.

References pat::JetCorrFactors::BOTTOM, pat::JetCorrFactors::CHARM, pat::JetCorrFactors::GLUON, testEve_cfg::level, python.rootplot.argparse::message, AlCaHLTBitMon_QueryRunRegistry::string, type_, and pat::JetCorrFactors::UDS.

Referenced by JetCorrFactorsProducer().

 {
   std::vector<std::string> expand;
   for(std::vector<std::string>::const_iterator level=levels.begin(); level!=levels.end(); ++level){
     if((*level)=="L5Flavor" || (*level)=="L7Parton"){
       if(flavor==JetCorrFactors::GLUON ){
     if(*level=="L7Parton" && type_=="T"){
       edm::LogWarning message( "L7Parton::GLUON not available" );
       message << "Jet energy corrections requested for level: L7Parton and type: 'T'. \n"
           << "For this combination there is no GLUON correction available. The    \n"
           << "correction for this flavor type will be taken from 'J'.";
     }
     expand.push_back(std::string(*level).append("_").append("g").append("J"));
       }
       if(flavor==JetCorrFactors::UDS   ) expand.push_back(std::string(*level).append("_").append("q").append(type_));
       if(flavor==JetCorrFactors::CHARM ) expand.push_back(std::string(*level).append("_").append("c").append(type_));
       if(flavor==JetCorrFactors::BOTTOM) expand.push_back(std::string(*level).append("_").append("b").append(type_));
     }
     else{
       expand.push_back(*level);
     }
   }
   return expand;
 }

void JetCorrFactorsProducer::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

description of configuration file parameters

Definition at line 297 of file JetCorrFactorsProducer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), HLT_25ns14e33_v1_cff::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

 {
   edm::ParameterSetDescription iDesc;
   iDesc.add<bool>("emf", false);
   iDesc.add<std::string>("flavorType", "J");
   iDesc.add<edm::InputTag>("src", edm::InputTag("ak5CaloJets"));
   iDesc.add<std::string>("payload", "AK5Calo");
   iDesc.add<bool>("useNPV", true);
   iDesc.add<edm::InputTag>("primaryVertices", edm::InputTag("offlinePrimaryVertices"));
   iDesc.add<bool>("useRho", true);
   iDesc.add<edm::InputTag>("rho", edm::InputTag("fixedGridRhoFastjetAllCalo"));
   iDesc.add<std::string>("extraJPTOffset", "L1Offset");
 
   std::vector<std::string> levels;
   levels.push_back(std::string("L1Offset"  ));
   levels.push_back(std::string("L2Relative"));
   levels.push_back(std::string("L3Absolute"));
   levels.push_back(std::string("L2L3Residual"));
   levels.push_back(std::string("L5Flavor"  ));
   levels.push_back(std::string("L7Parton"  ));
   iDesc.add<std::vector<std::string> >("levels", levels);
   descriptions.add("JetCorrFactorsProducer", iDesc);
 }

bool pat::JetCorrFactorsProducer::flavorDependent ( ) const

inlineprivate

return true if the jec levels contain at least one flavor dependent correction level

Definition at line 78 of file JetCorrFactorsProducer.h.

References levels_.

Referenced by produce().

78 { return (levels_.size()>1); };

pat::JetCorrFactorsProducer::levels_

FlavorCorrLevelMap levels_

Definition: JetCorrFactorsProducer.h:125

int pat::JetCorrFactorsProducer::numberOf ( const edm::Handle< std::vector< reco::Vertex > > & primaryVertices )

inlineprivate

determines the number of valid primary vertices for the standard L1Offset correction of JetMET

Definition at line 135 of file JetCorrFactorsProducer.h.

References HLT_25ns14e33_v1_cff::primaryVertices, and MetAnalyzer::pv().

Referenced by produce().

   {
     int npv=0;
     for(std::vector<reco::Vertex>::const_iterator pv=primaryVertices->begin(); pv!=primaryVertices->end(); ++pv){
       if(pv->ndof()>=4) ++npv;
     }
     return npv;
   }

std::vector< JetCorrectorParameters > JetCorrFactorsProducer::params	(	const JetCorrectorParametersCollection &	parameters,
		const std::vector< std::string > &	levels
	)		const

private

return the jec parameters as input to the FactorizedJetCorrector for different flavors

Definition at line 138 of file JetCorrFactorsProducer.cc.

References testEve_cfg::level, and JetCorrectorParametersCollection::push_back().

Referenced by produce().

 {
   std::vector<JetCorrectorParameters> params;
   for(std::vector<std::string>::const_iterator level=levels.begin(); level!=levels.end(); ++level){
     const JetCorrectorParameters& ip = parameters[*level]; //ip.printScreen();
     params.push_back(ip);
   }
   return params;
 }

std::string JetCorrFactorsProducer::payload ( )

private

map jet algorithm to payload in DB

Definition at line 176 of file JetCorrFactorsProducer.cc.

References payload_.

Referenced by produce().

 {
   return payload_;
 }

void JetCorrFactorsProducer::produce	(	edm::Event &	event,
		const edm::EventSetup &	setup
	)

overridevirtual

everything that needs to be done per event

Implements edm::stream::EDProducerBase.

Definition at line 182 of file JetCorrFactorsProducer.cc.

References cacheId_, correctors_, evaluate(), Exception, extraJPTOffset_, extraJPTOffsetCorrector_, edm::helper::Filler< Map >::fill(), flavorDependent(), edm::EventSetup::get(), edm::Event::getByToken(), customizeTrackingMonitorSeedNumber::idx, edm::helper::Filler< Map >::insert(), patTestJEC_cfi::jec, metsig::jet, fwrapper::jets, edm::InputTag::label(), label_, levels_, numberOf(), Parameters::parameters, params(), payload(), HLT_25ns14e33_v1_cff::primaryVertices, primaryVertices_, primaryVerticesToken_, rho, rho_, rhoToken_, srcToken_, and AlCaHLTBitMon_QueryRunRegistry::string.

 {
   // get jet collection from the event
   edm::Handle<edm::View<reco::Jet> > jets;
   event.getByToken(srcToken_, jets);
 
   // get primary vertices for L1Offset correction level if needed
   edm::Handle<std::vector<reco::Vertex> > primaryVertices;
   if(!primaryVertices_.label().empty()) event.getByToken(primaryVerticesToken_, primaryVertices);
 
   // get parameter rho for L1FastJet correction level if needed
   edm::Handle<double> rho;
   if(!rho_.label().empty()) event.getByToken(rhoToken_, rho);
 
   auto const& rec = setup.get<JetCorrectionsRecord>();
   if (cacheId_ != rec.cacheIdentifier()) {
     // retreive parameters from the DB
     edm::ESHandle<JetCorrectorParametersCollection> parameters;
     setup.get<JetCorrectionsRecord>().get(payload(), parameters);
     // initialize jet correctors
     for(FlavorCorrLevelMap::const_iterator flavor=levels_.begin(); flavor!=levels_.end(); ++flavor){
       correctors_[flavor->first].reset( new FactorizedJetCorrector(params(*parameters, flavor->second)) );
     }
     // initialize extra jet corrector for jpt if needed
     if(!extraJPTOffset_.empty()){
       extraJPTOffsetCorrector_.reset( new FactorizedJetCorrector(params(*parameters, extraJPTOffset_)) );
     }
     cacheId_ = rec.cacheIdentifier();
   }
 
   // fill the jetCorrFactors
   std::vector<JetCorrFactors> jcfs;
   for(edm::View<reco::Jet>::const_iterator jet = jets->begin(); jet!=jets->end(); ++jet){
     // the JetCorrFactors::CorrectionFactor is a std::pair<std::string, std::vector<float> >
     // the string corresponds to the label of the correction level, the vector contains four
     // floats if flavor dependent and one float else. Per construction jet energy corrections
     // will be flavor independent up to the first flavor dependent correction and flavor de-
     // pendent afterwards. The first correction level is predefined with label 'Uncorrected'.
     // Per definition it is flavor independent. The correction factor is 1.
     std::vector<JetCorrFactors::CorrectionFactor> jec;
     jec.push_back(std::make_pair<std::string, std::vector<float> >(std::string("Uncorrected"), std::vector<float>(1, 1)));
 
     // pick the first element in the map (which could be the only one) and loop all jec
     // levels listed for that element. If this is not the only element all jec levels, which
     // are flavor independent will give the same correction factors until the first flavor
     // dependent correction level is reached. So the first element is still a good choice.
     FlavorCorrLevelMap::const_iterator corrLevel=levels_.begin();
     if(corrLevel==levels_.end()){
       throw cms::Exception("No JECFactors") << "You request to create a jetCorrFactors object with no JEC Levels indicated. \n"
                         << "This makes no sense, either you should correct this or drop the module from \n"
                         << "the sequence.";
     }
     for(unsigned int idx=0; idx<corrLevel->second.size(); ++idx){
       bool flavorDependent=false;
       std::vector<float> factors;
       if(flavorDependent ||
      corrLevel->second[idx].find("L5Flavor")!=std::string::npos ||
      corrLevel->second[idx].find("L7Parton")!=std::string::npos){
     flavorDependent=true;
     // after the first encounter all subsequent correction levels are flavor dependent
     for(FlavorCorrLevelMap::const_iterator flavor=corrLevel; flavor!=levels_.end(); ++flavor){
       if(!primaryVertices_.label().empty()){
         // if primaryVerticesToken_ has a value the number of primary vertices needs to be
         // specified
         correctors_.find(flavor->first)->second->setNPV(numberOf(primaryVertices));
       }
       if(!rho_.label().empty()){
         // if rhoToken_ has a value the energy density parameter rho and the jet area need
         //  to be specified
         correctors_.find(flavor->first)->second->setRho(*rho);
         correctors_.find(flavor->first)->second->setJetA(jet->jetArea());
       }
       factors.push_back(evaluate(jet, flavor->first, idx));
     }
       }
       else{
     if(!primaryVertices_.label().empty()){
       // if primaryVerticesToken_ has a value the number of primary vertices needs to be
       // specified
       correctors_.find(corrLevel->first)->second->setNPV(numberOf(primaryVertices));
     }
     if(!rho_.label().empty()){
       // if rhoToken_ has a value the energy density parameter rho and the jet area need
       // to be specified
       correctors_.find(corrLevel->first)->second->setRho(*rho);
       correctors_.find(corrLevel->first)->second->setJetA(jet->jetArea());
     }
     factors.push_back(evaluate(jet, corrLevel->first, idx));
       }
       // push back the set of JetCorrFactors: the first entry corresponds to the label
       // of the correction level, which is taken from the first element in levels_. For
       // L5Flavor and L7Parton the part including the first '_' indicating the flavor
       // of the first element in levels_ is chopped of from the label to avoid confusion
       // of the correction levels. The second parameter corresponds to the set of jec
       // factors, which might be flavor dependent or not. In the default configuration
       // the CorrectionFactor will look like this: 'Uncorrected': 1 ; 'L2Relative': x ;
       // 'L3Absolute': x ; 'L5Flavor': v, x, y, z ; 'L7Parton': v, x, y, z
       jec.push_back(std::make_pair((corrLevel->second[idx]).substr(0, (corrLevel->second[idx]).find("_")), factors));
     }
     // create the actual object with the scale factors we want the valuemap to refer to
     // label_ corresponds to the label of the module instance
     JetCorrFactors corrFactors(label_, jec);
     jcfs.push_back(corrFactors);
   }
   // build the value map
   std::auto_ptr<JetCorrFactorsMap> jetCorrsMap(new JetCorrFactorsMap());
   JetCorrFactorsMap::Filler filler(*jetCorrsMap);
   // jets and jetCorrs have their indices aligned by construction
   filler.insert(jets, jcfs.begin(), jcfs.end());
   filler.fill(); // do the actual filling
   // put our produced stuff in the event
   event.put(jetCorrsMap);
 }

Member Data Documentation

unsigned long long pat::JetCorrFactorsProducer::cacheId_

private

cache identifier for JetCorrectionsRecord

Definition at line 127 of file JetCorrFactorsProducer.h.

Referenced by produce().

std::map<JetCorrFactors::Flavor, std::unique_ptr<FactorizedJetCorrector> > pat::JetCorrFactorsProducer::correctors_

private

cache container for jet corrections

Definition at line 129 of file JetCorrFactorsProducer.h.

Referenced by evaluate(), and produce().

bool pat::JetCorrFactorsProducer::emf_

private

use electromagnetic fraction for jet energy corrections or not (will only have an effect for jets CaloJets)

Definition at line 94 of file JetCorrFactorsProducer.h.

Referenced by evaluate().

std::vector<std::string> pat::JetCorrFactorsProducer::extraJPTOffset_

private

label of additional L1Offset corrector for JPT jets; for format reasons this string is kept in a vector of strings

Definition at line 105 of file JetCorrFactorsProducer.h.

Referenced by JetCorrFactorsProducer(), and produce().

std::unique_ptr<FactorizedJetCorrector> pat::JetCorrFactorsProducer::extraJPTOffsetCorrector_

private

cache container for JPTOffset jet corrections

Definition at line 131 of file JetCorrFactorsProducer.h.

Referenced by evaluate(), and produce().

std::string pat::JetCorrFactorsProducer::label_

private

label of jec factors

Definition at line 100 of file JetCorrFactorsProducer.h.

Referenced by Mixins._Labelable::_findDependencies(), produce(), and Mixins._Labelable::setLabel().

FlavorCorrLevelMap pat::JetCorrFactorsProducer::levels_

private

jec levels for different flavors. In the default configuration this map would look like this: GLUON : 'L2Relative', 'L3Absolute', 'L5FLavor_jg', L7Parton_jg' UDS : 'L2Relative', 'L3Absolute', 'L5FLavor_jq', L7Parton_jq' CHARM : 'L2Relative', 'L3Absolute', 'L5FLavor_jc', L7Parton_jc' BOTTOM : 'L2Relative', 'L3Absolute', 'L5FLavor_jb', L7Parton_jb' or just like this: NONE : 'L2Relative', 'L3Absolute', 'L2L3Residual' per definition the vectors for all elements in this map should have the same size

Definition at line 125 of file JetCorrFactorsProducer.h.

Referenced by flavorDependent(), JetCorrFactorsProducer(), and produce().

std::string pat::JetCorrFactorsProducer::payload_

private

label of payload

Definition at line 102 of file JetCorrFactorsProducer.h.

Referenced by payload().

edm::InputTag pat::JetCorrFactorsProducer::primaryVertices_

private

label for L1Offset primaryVertex collection

Definition at line 107 of file JetCorrFactorsProducer.h.

Referenced by JetCorrFactorsProducer(), and produce().