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.12 by amagnan, Fri Mar 26 17:52:50 2010 UTC

# Line 19 | Line 19
19   #include "DataFormats/PatCandidates/interface/MET.h"
20  
21   #include "FWCore/ServiceRegistry/interface/Service.h"
22 < #include "PhysicsTools/UtilAlgos/interface/TFileService.h"
22 > #include "CommonTools/UtilAlgos/interface/TFileService.h"
23  
24  
25   #include "UserCode/HbbAnalysis/interface/Electron.hh"
# 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 381 | Line 385 | void HbbTreeMaker::HbbElectrons(const ed
385            lID.hOverE = (*iter).hadronicOverEm();
386            lID.deltaPhiIn = (*iter).deltaPhiSuperClusterTrackAtVtx();
387            lID.deltaEtaIn = (*iter).deltaEtaSuperClusterTrackAtVtx();
388 +          lID.ecalDrivenSeed = (*iter).ecalDrivenSeed();
389 +          lID.trackerDrivenSeed = (*iter).trackerDrivenSeed();
390  
391            HbbAnalysis::Electron lObj(lGen,lReco,lSC,lIso,lID);
392            aVec.push_back(lObj);
# Line 944 | Line 950 | void HbbTreeMaker::HbbJets(const edm::Ha
950            lCommon.rawpT = (*iter).pt();
951            if ((*iter).hasCorrFactors())
952              lCommon.rawpT = (*iter).pt()/((*iter).corrFactor((*iter).corrStep()));
953 +
954 +          lCommon.etaMean = (*iter).etaPhiStatistics().etaMean;
955 +          lCommon.phiMean = (*iter).etaPhiStatistics().phiMean;
956 +          lCommon.etaEtaMoment = (*iter).etaPhiStatistics().etaEtaMoment;
957 +          lCommon.phiPhiMoment = (*iter).etaPhiStatistics().phiPhiMoment;
958 +          lCommon.etaPhiMoment = (*iter).etaPhiStatistics().etaPhiMoment;
959 +          
960            //lCommon.rawEta = (*iter).;
961            //lCommon.rawPhi = (*iter).;
962            //    p_hasJetCorrFactors->Fill(aJet.hasJetCorrFactors());
# Line 955 | Line 968 | void HbbTreeMaker::HbbJets(const edm::Ha
968            lBtag.bProba = (*iter).bDiscriminator("jetBProbabilityBJetTags");
969            lBtag.probability = (*iter).bDiscriminator("jetProbabilityBJetTags");
970            lBtag.sSV = (*iter).bDiscriminator("simpleSecondaryVertexBJetTags");
971 <          lBtag.softElectronByPt = (*iter).bDiscriminator("softElectronByPtBJetTags");
972 <          lBtag.softElectronByIP3d = (*iter).bDiscriminator("softElectronByIP3dBJetTags");
973 <          lBtag.softMuon = (*iter).bDiscriminator("softMuonBJetTags");
974 <          lBtag.softMuonByPt = (*iter).bDiscriminator("softMuonByPtBJetTags");
975 <          lBtag.softMuonByIP3d = (*iter).bDiscriminator("softMuonByIP3dBJetTags");
971 >          lBtag.softElectronByPt = ((*iter).bDiscriminator("softElectronByPtBJetTags") < -1000) ? -1000 : (*iter).bDiscriminator("softElectronByPtBJetTags");
972 >          lBtag.softElectronByIP3d = ((*iter).bDiscriminator("softElectronByIP3dBJetTags") < -1000) ? -1000 : (*iter).bDiscriminator("softElectronByIP3dBJetTags");
973 >          lBtag.softMuon = ((*iter).bDiscriminator("softMuonBJetTags") < -1000) ? -1000 : (*iter).bDiscriminator("softMuonBJetTags");
974 >          lBtag.softMuonByPt = ((*iter).bDiscriminator("softMuonByPtBJetTags") < -1000) ? -1000 : (*iter).bDiscriminator("softMuonByPtBJetTags");
975 >          lBtag.softMuonByIP3d = ((*iter).bDiscriminator("softMuonByIP3dBJetTags") < -1000) ? -1000 : (*iter).bDiscriminator("softMuonByIP3dBJetTags");
976 >
977 >          //std::cout << " -- New values of b-discri for jet : " << iEle << std::endl
978 >          //        << " ---- softElecs: " << lBtag.softElectronByPt << " " << lBtag.softElectronByIP3d << std::endl
979 >          //        << " ---- softMus: " << lBtag.softMuon << " " << lBtag.softMuonByPt << " " << lBtag.softMuonByIP3d << std::endl;
980 >
981 >
982            lBtag.tCHE = (*iter).bDiscriminator("trackCountingHighEffBJetTags");
983            lBtag.tCHP = (*iter).bDiscriminator("trackCountingHighPurBJetTags");
984  
# Line 1158 | Line 1177 | void HbbTreeMaker::HbbParticles(const ed
1177  
1178  
1179  
1161 –
1162 –
1180   void HbbTreeMaker::HbbVertices(const edm::Handle<std::vector<reco::Vertex> > & aCol,
1181                                 std::vector<HbbAnalysis::Vertex> & aVec)
1182   {
# Line 1207 | Line 1224 | void HbbTreeMaker::HbbVertices(const edm
1224      }
1225  
1226   }
1227 +
1228 + #include "FWCore/Framework/interface/MakerMacros.h"
1229 + DEFINE_FWK_MODULE(HbbTreeMaker);
1230 +
1231 +

Diff Legend

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