38 :
mInput (fConfig.getParameter <
edm::InputTag> (
"src")),
57 for (
unsigned ijet = 0; ijet < jets->size(); ++ijet) {
59 std::cout <<
"PartonJetCorrectionExample::analize-> jet #" << ijet;
61 std::cout <<
": use gJ corrections" << std::endl;
62 corrector = gJetCorrector;
64 else if (ijet%4 == 1) {
65 std::cout <<
": use qJ corrections" << std::endl;
66 corrector = qJetCorrector;
68 else if (ijet%4 == 2) {
69 std::cout <<
": use bJ corrections" << std::endl;
70 corrector = bJetCorrector;
73 std::cout <<
": use bT corrections" << std::endl;
74 corrector = bTopCorrector;
77 double correction = corrector->
correction (jet);
79 std::cout <<
" jet pt/eta/phi: " << jet.
pt() <<
'/' << jet.
eta() <<
'/' << jet.
phi()
80 <<
" -> correction factor: " << correction
81 <<
", corrected pt: " << jet.
pt()*correction
std::string m_bT_CorrectorName
void analyze(const edm::Event &, const edm::EventSetup &) override
double eta() const final
momentum pseudorapidity
PartonJetCorrectionExample(const edm::ParameterSet &fParameters)
Jets made from CaloTowers.
~PartonJetCorrectionExample() override
double pt() const final
transverse momentum
double correction(const LorentzVector &fJet) const
get correction using Jet information only
std::string m_qJ_CorrectorName
std::string m_gJ_CorrectorName
#define DEFINE_FWK_MODULE(type)
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
std::string m_bJ_CorrectorName
static const JetCorrector * getJetCorrector(const std::string &fName, const edm::EventSetup &fSetup)
retrieve corrector from the event setup. troughs exception if something is missing ...
double phi() const final
momentum azimuthal angle