ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/HbbAnalysis/plugins/HbbTreeMaker.cc
(Generate patch)

Comparing UserCode/HbbAnalysis/plugins/HbbTreeMaker.cc (file contents):
Revision 1.5 by amagnan, Wed Feb 10 13:50:17 2010 UTC vs.
Revision 1.9 by amagnan, Tue Mar 2 12:11:54 2010 UTC

# Line 111 | Line 111 | void HbbTreeMaker::analyze(const edm::Ev
111  
112    event_->run(aEvt.run());
113    event_->lumiBlock(aEvt.luminosityBlock());
114 +  event_->isRealData(aEvt.isRealData());
115  
116    edm::Handle<reco::GenParticleCollection> lGenParticles;
117    try {
# Line 213 | Line 214 | void HbbTreeMaker::analyze(const edm::Ev
214      std::cout << "AMM: Collection " << caloJetSrc_  << " not available! Exception : " << e.what() << ". " << std::endl;
215    }
216  
217 +  //std::cout << "Processing calo jets:" << std::endl;
218    HbbJets(lCaloJetCollection,jetFlav_,lGenParticles,event_->caloJets());
219    
220    edm::Handle<std::vector<pat::Jet> > lJptJetCollection;
# Line 227 | Line 229 | void HbbTreeMaker::analyze(const edm::Ev
229      std::cout << "AMM: Collection " << jptJetSrc_  << " not available! Exception : " << e.what() << ". " << std::endl;
230    }
231  
232 +  //std::cout << "Processing JPT jets:" << std::endl;
233    HbbJets(lJptJetCollection,jetFlav_,lGenParticles,event_->jptJets());
234    
235    edm::Handle<std::vector<pat::Jet> > lPfJetCollection;
# Line 241 | Line 244 | void HbbTreeMaker::analyze(const edm::Ev
244      std::cout << "AMM: Collection " << pfJetSrc_  << " not available! Exception : " << e.what() << ". " << std::endl;
245    }
246  
247 +  //std::cout << "Processing PF jets:" << std::endl;
248    HbbJets(lPfJetCollection,jetFlav_,lGenParticles,event_->pfJets());
249    
250    edm::Handle<std::vector<pat::MET> > lCaloMetCol;
# Line 955 | Line 959 | void HbbTreeMaker::HbbJets(const edm::Ha
959            lBtag.bProba = (*iter).bDiscriminator("jetBProbabilityBJetTags");
960            lBtag.probability = (*iter).bDiscriminator("jetProbabilityBJetTags");
961            lBtag.sSV = (*iter).bDiscriminator("simpleSecondaryVertexBJetTags");
962 <          lBtag.softElectronByPt = (*iter).bDiscriminator("softElectronByPtBJetTags");
963 <          lBtag.softElectronByIP3d = (*iter).bDiscriminator("softElectronByIP3dBJetTags");
964 <          lBtag.softMuon = (*iter).bDiscriminator("softMuonBJetTags");
965 <          lBtag.softMuonByPt = (*iter).bDiscriminator("softMuonByPtBJetTags");
966 <          lBtag.softMuonByIP3d = (*iter).bDiscriminator("softMuonByIP3dBJetTags");
962 >          lBtag.softElectronByPt = ((*iter).bDiscriminator("softElectronByPtBJetTags") < -1000) ? -1000 : (*iter).bDiscriminator("softElectronByPtBJetTags");
963 >          lBtag.softElectronByIP3d = ((*iter).bDiscriminator("softElectronByIP3dBJetTags") < -1000) ? -1000 : (*iter).bDiscriminator("softElectronByIP3dBJetTags");
964 >          lBtag.softMuon = ((*iter).bDiscriminator("softMuonBJetTags") < -1000) ? -1000 : (*iter).bDiscriminator("softMuonBJetTags");
965 >          lBtag.softMuonByPt = ((*iter).bDiscriminator("softMuonByPtBJetTags") < -1000) ? -1000 : (*iter).bDiscriminator("softMuonByPtBJetTags");
966 >          lBtag.softMuonByIP3d = ((*iter).bDiscriminator("softMuonByIP3dBJetTags") < -1000) ? -1000 : (*iter).bDiscriminator("softMuonByIP3dBJetTags");
967 >
968 >          //std::cout << " -- New values of b-discri for jet : " << iEle << std::endl
969 >          //        << " ---- softElecs: " << lBtag.softElectronByPt << " " << lBtag.softElectronByIP3d << std::endl
970 >          //        << " ---- softMus: " << lBtag.softMuon << " " << lBtag.softMuonByPt << " " << lBtag.softMuonByIP3d << std::endl;
971 >
972 >
973            lBtag.tCHE = (*iter).bDiscriminator("trackCountingHighEffBJetTags");
974            lBtag.tCHP = (*iter).bDiscriminator("trackCountingHighPurBJetTags");
975  
# Line 1158 | Line 1168 | void HbbTreeMaker::HbbParticles(const ed
1168  
1169  
1170  
1161
1162
1171   void HbbTreeMaker::HbbVertices(const edm::Handle<std::vector<reco::Vertex> > & aCol,
1172                                 std::vector<HbbAnalysis::Vertex> & aVec)
1173   {
# Line 1207 | Line 1215 | void HbbTreeMaker::HbbVertices(const edm
1215      }
1216  
1217   }
1218 +
1219 + #include "FWCore/Framework/interface/MakerMacros.h"
1220 + DEFINE_FWK_MODULE(HbbTreeMaker);
1221 +
1222 +

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines