213 |
|
std::cout << "AMM: Collection " << caloJetSrc_ << " not available! Exception : " << e.what() << ". " << std::endl; |
214 |
|
} |
215 |
|
|
216 |
+ |
//std::cout << "Processing calo jets:" << std::endl; |
217 |
|
HbbJets(lCaloJetCollection,jetFlav_,lGenParticles,event_->caloJets()); |
218 |
|
|
219 |
|
edm::Handle<std::vector<pat::Jet> > lJptJetCollection; |
228 |
|
std::cout << "AMM: Collection " << jptJetSrc_ << " not available! Exception : " << e.what() << ". " << std::endl; |
229 |
|
} |
230 |
|
|
231 |
+ |
//std::cout << "Processing JPT jets:" << std::endl; |
232 |
|
HbbJets(lJptJetCollection,jetFlav_,lGenParticles,event_->jptJets()); |
233 |
|
|
234 |
|
edm::Handle<std::vector<pat::Jet> > lPfJetCollection; |
243 |
|
std::cout << "AMM: Collection " << pfJetSrc_ << " not available! Exception : " << e.what() << ". " << std::endl; |
244 |
|
} |
245 |
|
|
246 |
+ |
//std::cout << "Processing PF jets:" << std::endl; |
247 |
|
HbbJets(lPfJetCollection,jetFlav_,lGenParticles,event_->pfJets()); |
248 |
|
|
249 |
|
edm::Handle<std::vector<pat::MET> > lCaloMetCol; |
958 |
|
lBtag.bProba = (*iter).bDiscriminator("jetBProbabilityBJetTags"); |
959 |
|
lBtag.probability = (*iter).bDiscriminator("jetProbabilityBJetTags"); |
960 |
|
lBtag.sSV = (*iter).bDiscriminator("simpleSecondaryVertexBJetTags"); |
961 |
< |
lBtag.softElectronByPt = (*iter).bDiscriminator("softElectronByPtBJetTags"); |
962 |
< |
lBtag.softElectronByIP3d = (*iter).bDiscriminator("softElectronByIP3dBJetTags"); |
963 |
< |
lBtag.softMuon = (*iter).bDiscriminator("softMuonBJetTags"); |
964 |
< |
lBtag.softMuonByPt = (*iter).bDiscriminator("softMuonByPtBJetTags"); |
965 |
< |
lBtag.softMuonByIP3d = (*iter).bDiscriminator("softMuonByIP3dBJetTags"); |
961 |
> |
lBtag.softElectronByPt = ((*iter).bDiscriminator("softElectronByPtBJetTags") < -1000) ? -1000 : (*iter).bDiscriminator("softElectronByPtBJetTags"); |
962 |
> |
lBtag.softElectronByIP3d = ((*iter).bDiscriminator("softElectronByIP3dBJetTags") < -1000) ? -1000 : (*iter).bDiscriminator("softElectronByIP3dBJetTags"); |
963 |
> |
lBtag.softMuon = ((*iter).bDiscriminator("softMuonBJetTags") < -1000) ? -1000 : (*iter).bDiscriminator("softMuonBJetTags"); |
964 |
> |
lBtag.softMuonByPt = ((*iter).bDiscriminator("softMuonByPtBJetTags") < -1000) ? -1000 : (*iter).bDiscriminator("softMuonByPtBJetTags"); |
965 |
> |
lBtag.softMuonByIP3d = ((*iter).bDiscriminator("softMuonByIP3dBJetTags") < -1000) ? -1000 : (*iter).bDiscriminator("softMuonByIP3dBJetTags"); |
966 |
> |
|
967 |
> |
//std::cout << " -- New values of b-discri for jet : " << iEle << std::endl |
968 |
> |
// << " ---- softElecs: " << lBtag.softElectronByPt << " " << lBtag.softElectronByIP3d << std::endl |
969 |
> |
// << " ---- softMus: " << lBtag.softMuon << " " << lBtag.softMuonByPt << " " << lBtag.softMuonByIP3d << std::endl; |
970 |
> |
|
971 |
> |
|
972 |
|
lBtag.tCHE = (*iter).bDiscriminator("trackCountingHighEffBJetTags"); |
973 |
|
lBtag.tCHP = (*iter).bDiscriminator("trackCountingHighPurBJetTags"); |
974 |
|
|
1167 |
|
|
1168 |
|
|
1169 |
|
|
1161 |
– |
|
1162 |
– |
|
1170 |
|
void HbbTreeMaker::HbbVertices(const edm::Handle<std::vector<reco::Vertex> > & aCol, |
1171 |
|
std::vector<HbbAnalysis::Vertex> & aVec) |
1172 |
|
{ |