CMS 3D CMS Logo

PartonJetCorrectionExample.cc
Go to the documentation of this file.
1 /* Example analizer to use L7 flavor JetCorrector services
2  Applies different parton corrections randomly
3  original from F.Ratnikov (UMd) Nov 16, 2007 - Adapted deom A.Santocchia Mar 01, 2008
4 */
5 
6 #include <string>
7 
10 
12  public:
13  explicit PartonJetCorrectionExample (const edm::ParameterSet& fParameters);
15  void analyze(const edm::Event&, const edm::EventSetup&) override;
16  private:
22 };
23 
24 
25 
26 
33 
34 using namespace std;
35 using namespace reco;
36 
38  : mInput (fConfig.getParameter <edm::InputTag> ("src")),
39  m_gJ_CorrectorName (fConfig.getParameter <std::string> ("gJetCorrector")),
40  m_qJ_CorrectorName (fConfig.getParameter <std::string> ("qJetCorrector")),
41  m_bJ_CorrectorName (fConfig.getParameter <std::string> ("bJetCorrector")),
42  m_bT_CorrectorName (fConfig.getParameter <std::string> ("bTopCorrector"))
43 {}
44 
46  // get all correctors
47  const JetCorrector* gJetCorrector = JetCorrector::getJetCorrector (m_gJ_CorrectorName, fSetup);
48  const JetCorrector* qJetCorrector = JetCorrector::getJetCorrector (m_qJ_CorrectorName, fSetup);
49  const JetCorrector* bJetCorrector = JetCorrector::getJetCorrector (m_bJ_CorrectorName, fSetup);
50  const JetCorrector* bTopCorrector = JetCorrector::getJetCorrector (m_bT_CorrectorName, fSetup);
51  const JetCorrector* corrector = nullptr;
52 
53  // get input jets (supposed to be MC corrected already)
55  fEvent.getByLabel (mInput, jets);
56  // loop over jets
57  for (unsigned ijet = 0; ijet < jets->size(); ++ijet) {
58  const CaloJet& jet = (*jets)[ijet];
59  std::cout << "PartonJetCorrectionExample::analize-> jet #" << ijet;
60  if (ijet%4 == 0) { // assume it is gluon from diJet
61  std::cout << ": use gJ corrections" << std::endl;
62  corrector = gJetCorrector;
63  }
64  else if (ijet%4 == 1) { // assume it is light quark from diJet
65  std::cout << ": use qJ corrections" << std::endl;
66  corrector = qJetCorrector;
67  }
68  else if (ijet%4 == 2) { // assume it is b quark from diJet
69  std::cout << ": use bJ corrections" << std::endl;
70  corrector = bJetCorrector;
71  }
72  else { // assume it is b quark from ttbar
73  std::cout << ": use bT corrections" << std::endl;
74  corrector = bTopCorrector;
75  }
76  // get selected correction for the jet
77  double correction = corrector->correction (jet);
78  // dump it
79  std::cout << " jet pt/eta/phi: " << jet.pt() << '/' << jet.eta() << '/' << jet.phi()
80  << " -> correction factor: " << correction
81  << ", corrected pt: " << jet.pt()*correction
82  << std::endl;
83  }
84 }
85 
88 
89 
void analyze(const edm::Event &, const edm::EventSetup &) override
double eta() const final
momentum pseudorapidity
PartonJetCorrectionExample(const edm::ParameterSet &fParameters)
Jets made from CaloTowers.
Definition: CaloJet.h:29
double pt() const final
transverse momentum
double correction(const LorentzVector &fJet) const
get correction using Jet information only
Definition: JetCorrector.h:49
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:16
vector< PseudoJet > jets
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:480
static const JetCorrector * getJetCorrector(const std::string &fName, const edm::EventSetup &fSetup)
retrieve corrector from the event setup. troughs exception if something is missing ...
Definition: JetCorrector.cc:50
fixed size matrix
HLT enums.
double phi() const final
momentum azimuthal angle