ViewVC Help

View File | Revision Log | Show Annotations | Root Listing

root/cvsroot/UserCode/HbbAnalysis/src/HistosElecs.cc

Comparing UserCode/HbbAnalysis/src/HistosElecs.cc (file contents):
Revision 1.3 by amagnan, Mon May 18 14:59:04 2009 UTC vs.
Revision 1.5 by amagnan, Thu May 28 15:12:45 2009 UTC

#	Line 18 | Line 18 | namespace HbbAnalysis {//namespace
18
19		p_scSigmaEtaEta = aDir.make<TH1F>("p_scSigmaEtaEta",";;N_{entries}",100,0,0.1);
20		p_scSigmaIEtaIEta = aDir.make<TH1F>("p_scSigmaIEtaIEta",";;N_{entries}",100,0,0.1);
21	<	p_scE1x5 = aDir.make<TH1F>("p_scE1x5",";;N_{entries}",200,0,200);
22	<	p_scE2x5Max = aDir.make<TH1F>("p_scE2x5Max",";;N_{entries}",200,0,200);
21	>	p_scE1x5 = aDir.make<TH1F>("p_scE1x5",";;N_{entries}",500,0,500);
22	>	p_scE2x5Max = aDir.make<TH1F>("p_scE2x5Max",";;N_{entries}",500,0,500);
23		p_scE5x5 = aDir.make<TH1F>("p_scE5x5",";;N_{entries}",500,0,500);
24	+	p_scE1x5OverscE5x5 = aDir.make<TH1F>("p_scE1x5OverscE5x5",";;N_{entries}",500,0,1);
25	+	p_scE2x5MaxOverscE5x5 = aDir.make<TH1F>("p_scE2x5MaxOverscE5x5",";;N_{entries}",500,0,1);
26
27		p_trackIso = aDir.make<TH1F>("p_trackIso",";;N_{entries}",400,0,20);
28		p_caloIso = aDir.make<TH1F>("p_caloIso",";;N_{entries}",400,0,20);
#	Line 34 | Line 36 | namespace HbbAnalysis {//namespace
36
37		}
38
39	<	void HistosElecs::FillEventHistograms(edm::Handle<std::vector<pat::Electron> > aCol){
39	>	void HistosElecs::FillEventHistograms(const edm::Handle<std::vector<pat::Electron> > & aCol){
40
41		if (doGenMatched_) {
42		unsigned int nMatched = 0;
#	Line 50 | Line 52 | namespace HbbAnalysis {//namespace
52
53		}
54
55	<	void HistosElecs::FillHistograms(const pat::Electron & aElec, edm::Handle<std::vector<reco::Vertex> > aRecoVertices, bool isLead){//FillHistograms
55	>	void HistosElecs::FillHistograms(const pat::Electron & aElec, const edm::Handle<std::vector<reco::Vertex> >& aRecoVertices, bool isLead){//FillHistograms
56
57		bool lIsGenMatched = !doGenMatched_ || (doGenMatched_ && MatchesGenElectron(aElec));
58		if (lIsGenMatched) {//genMatched
#	Line 70 | Line 72 | namespace HbbAnalysis {//namespace
72		if (lIds.size() > 5) std::cout << "--- WARNING: histo will be out-of-range, number of electron IDs = " << lIds.size() << std::endl;
73		for (unsigned int iId(0); iId<lIds.size(); iId++){
74		const std::string & lName =  lIds.at(iId).first;
75	+	//std::cout << "--- ElectronID : Bin " << iId << ", id name = " << lName << std::endl;
76		if (aElec.isElectronIDAvailable(lName)){
77		p_electronID->Fill(2*iId+aElec.electronID(lName));
78		assert (aElec.electronID(lName) ==  lIds.at(iId).second);
#	Line 82 | Line 85 | namespace HbbAnalysis {//namespace
85		p_scE2x5Max->Fill(aElec.scE2x5Max());
86		p_scE5x5->Fill(aElec.scE5x5());
87
88	+	if (aElec.scE5x5() > 0) {
89	+	p_scE1x5OverscE5x5->Fill(aElec.scE1x5()/aElec.scE5x5());
90	+	p_scE2x5MaxOverscE5x5->Fill(aElec.scE2x5Max()/aElec.scE5x5());
91	+	}
92	+
93		p_trackIso->Fill(aElec.trackIso());
94		p_caloIso->Fill(aElec.caloIso());
95		p_ecalIso->Fill(aElec.ecalIso());

Diff Legend

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