75 |
|
pfTauSrc_(pset.getParameter<edm::InputTag>("PFTaus")), |
76 |
|
caloJetSrc_(pset.getParameter<edm::InputTag>("CaloJets")), |
77 |
|
jptJetSrc_(pset.getParameter<edm::InputTag>("JPTJets")), |
78 |
+ |
ak7JetSrc_(pset.getParameter<edm::InputTag>("AK7Jets")), |
79 |
|
pfJetSrc_(pset.getParameter<edm::InputTag>("PFJets")), |
80 |
|
caloMetSrc_(pset.getParameter<edm::InputTag>("CaloMET")), |
81 |
|
tcMetSrc_(pset.getParameter<edm::InputTag>("TCMET")), |
340 |
|
//std::cout << "Processing JPT jets:" << std::endl; |
341 |
|
HbbJets(lJptJetCollection,jetFlav_,lGenParticles,event_->jptJets()); |
342 |
|
|
343 |
+ |
edm::Handle<std::vector<pat::Jet> > lAk7JetCollection; |
344 |
+ |
|
345 |
+ |
try { |
346 |
+ |
aEvt.getByLabel(ak7JetSrc_,lAk7JetCollection); |
347 |
+ |
if (!lAkJetCollection.isValid()){ |
348 |
+ |
edm::LogInfo("ERROR")<< "Error! Can't get ak7Jets by label. "; |
349 |
+ |
} |
350 |
+ |
if (debug_) std::cout << "** ak7Jetcollection = " << lAk7JetCollection->size() << " elements." << std::endl; |
351 |
+ |
} catch(cms::Exception& e) { |
352 |
+ |
std::cout << "AMM: Collection " << ak7JetSrc_ << " not available! Exception : " << e.what() << ". " << std::endl; |
353 |
+ |
} |
354 |
+ |
|
355 |
+ |
//std::cout << "Processing JPT jets:" << std::endl; |
356 |
+ |
HbbJets(lAk7JetCollection,jetFlav_,lGenParticles,event_->ak7Jets()); |
357 |
+ |
|
358 |
|
edm::Handle<std::vector<pat::Jet> > lPfJetCollection; |
359 |
|
|
360 |
|
try { |
1514 |
|
lVtx.x = (*iter).x(); |
1515 |
|
lVtx.y = (*iter).y(); |
1516 |
|
lVtx.z = (*iter).z(); |
1517 |
+ |
lVtx.rho = (*iter).position().Rho(); |
1518 |
|
lVtx.xError = (*iter).xError(); |
1519 |
|
lVtx.yError = (*iter).yError(); |
1520 |
|
lVtx.zError = (*iter).zError(); |