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.6 by amagnan, Wed Feb 10 14:40:15 2010 UTC vs.
Revision 1.10 by amagnan, Tue Mar 2 16:03:34 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 944 | Line 948 | void HbbTreeMaker::HbbJets(const edm::Ha
948            lCommon.rawpT = (*iter).pt();
949            if ((*iter).hasCorrFactors())
950              lCommon.rawpT = (*iter).pt()/((*iter).corrFactor((*iter).corrStep()));
951 +
952 +          lCommon.etaMean = (*iter).etaPhiStatistics().etaMean;
953 +          lCommon.phiMean = (*iter).etaPhiStatistics().phiMean;
954 +          lCommon.etaEtaMoment = (*iter).etaPhiStatistics().etaEtaMoment;
955 +          lCommon.phiPhiMoment = (*iter).etaPhiStatistics().phiPhiMoment;
956 +          lCommon.etaPhiMoment = (*iter).etaPhiStatistics().etaPhiMoment;
957 +          
958            //lCommon.rawEta = (*iter).;
959            //lCommon.rawPhi = (*iter).;
960            //    p_hasJetCorrFactors->Fill(aJet.hasJetCorrFactors());
# Line 955 | Line 966 | void HbbTreeMaker::HbbJets(const edm::Ha
966            lBtag.bProba = (*iter).bDiscriminator("jetBProbabilityBJetTags");
967            lBtag.probability = (*iter).bDiscriminator("jetProbabilityBJetTags");
968            lBtag.sSV = (*iter).bDiscriminator("simpleSecondaryVertexBJetTags");
969 <          lBtag.softElectronByPt = (*iter).bDiscriminator("softElectronByPtBJetTags");
970 <          lBtag.softElectronByIP3d = (*iter).bDiscriminator("softElectronByIP3dBJetTags");
971 <          lBtag.softMuon = (*iter).bDiscriminator("softMuonBJetTags");
972 <          lBtag.softMuonByPt = (*iter).bDiscriminator("softMuonByPtBJetTags");
973 <          lBtag.softMuonByIP3d = (*iter).bDiscriminator("softMuonByIP3dBJetTags");
974 <
975 <          std::cout << " -- New values of b-discri for jet : " << iEle << std::endl
976 <                    << " ---- softElecs: " << lBtag.softElectronByPt << " " << lBtag.softElectronByIP3d << std::endl
977 <                    << " ---- softMus: " << lBtag.softMuon << " " << lBtag.softMuonByPt << " " << lBtag.softMuonByIP3d << std::endl;
969 >          lBtag.softElectronByPt = ((*iter).bDiscriminator("softElectronByPtBJetTags") < -1000) ? -1000 : (*iter).bDiscriminator("softElectronByPtBJetTags");
970 >          lBtag.softElectronByIP3d = ((*iter).bDiscriminator("softElectronByIP3dBJetTags") < -1000) ? -1000 : (*iter).bDiscriminator("softElectronByIP3dBJetTags");
971 >          lBtag.softMuon = ((*iter).bDiscriminator("softMuonBJetTags") < -1000) ? -1000 : (*iter).bDiscriminator("softMuonBJetTags");
972 >          lBtag.softMuonByPt = ((*iter).bDiscriminator("softMuonByPtBJetTags") < -1000) ? -1000 : (*iter).bDiscriminator("softMuonByPtBJetTags");
973 >          lBtag.softMuonByIP3d = ((*iter).bDiscriminator("softMuonByIP3dBJetTags") < -1000) ? -1000 : (*iter).bDiscriminator("softMuonByIP3dBJetTags");
974 >
975 >          //std::cout << " -- New values of b-discri for jet : " << iEle << std::endl
976 >          //        << " ---- softElecs: " << lBtag.softElectronByPt << " " << lBtag.softElectronByIP3d << std::endl
977 >          //        << " ---- softMus: " << lBtag.softMuon << " " << lBtag.softMuonByPt << " " << lBtag.softMuonByIP3d << std::endl;
978  
979  
980            lBtag.tCHE = (*iter).bDiscriminator("trackCountingHighEffBJetTags");
# Line 1211 | Line 1222 | void HbbTreeMaker::HbbVertices(const edm
1222      }
1223  
1224   }
1225 +
1226 + #include "FWCore/Framework/interface/MakerMacros.h"
1227 + DEFINE_FWK_MODULE(HbbTreeMaker);
1228 +
1229 +

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines