38 void endJob()
override ;
57 mVertexToken = consumes<reco::VertexCollection>(
edm::InputTag(
"offlinePrimaryVertices"));
58 mMinRawJetPt = iConfig.
getParameter<
double> (
"MinRawJetPt");
79 for(
unsigned int ind=0;ind<recVtxs->size();ind++) {
80 if (!((*recVtxs)[ind].isFake())) {
84 typename JetCollection::const_iterator i_jet;
86 for(i_jet = jets->begin(); i_jet != jets->end(); i_jet++) {
87 if (i_jet->pt() < mMinRawJetPt)
continue;
91 std::cout<<
"energy = "<<i_jet->energy()<<
", " 92 <<
"eta = "<<i_jet->eta()<<
", " 93 <<
"raw pt = "<<i_jet->pt()<<
", " 95 <<
"correction = "<<scale<<
", " 96 <<
"cor pt = "<<scale*i_jet->pt()<<endl;
98 mRawPt->Fill(i_jet->pt());
99 mCorPt->Fill(scale*i_jet->pt());
108 mRawPt = fs->make<TH1F>(
"RawJetPt",
"RawJetPt",1000,0,1000);
109 mCorPt = fs->make<TH1F>(
"CorJetPt",
"CorJetPt",1000,0,1000);
T getParameter(std::string const &) const
edm::EDGetTokenT< JetCollection > mJetToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
~JetCorrectorOnTheFly() override
std::vector< Jet > JetCollection
double correction(const LorentzVector &fJet) const
get correction using Jet information only
JetCorrectorOnTheFly(const edm::ParameterSet &)
edm::EDGetTokenT< reco::VertexCollection > mVertexToken
example_stream void analyze(const edm::Event &, const edm::EventSetup &) override
edm::EDGetTokenT< reco::JetCorrector > mJetCorrector
#define DEFINE_FWK_MODULE(type)
JetCorrectorOnTheFly< PFJet > PFJetCorrectorOnTheFly
edm::Service< TFileService > fs
JetCorrectorOnTheFly< JPTJet > JPTJetCorrectorOnTheFly
void analyze(const edm::Event &, const edm::EventSetup &) override
JetCorrectorOnTheFly< CaloJet > CaloJetCorrectorOnTheFly