ViewVC Help

View File | Revision Log | Show Annotations | Root Listing

root/cvsroot/UserCode/MitHzz4l/Selection/src/applySelection.cc

Comparing UserCode/MitHzz4l/Selection/src/applySelection.cc (file contents):
Revision 1.1 by khahn, Thu Sep 8 13:33:18 2011 UTC vs.
Revision 1.5 by dkralph, Wed Oct 12 17:25:15 2011 UTC

#	Line 32 | Line 32 | using namespace std;
32		//
33		// ntuple format headers
34		//
35	<	#include "EWKAnaDefs.hh"
35	>	#include "HiggsAnaDefs.hh"
36		#include "TEventInfo.hh"
37		#include "TElectron.hh"
38		#include "TMuon.hh"
#	Line 43 | Line 43 | using namespace std;
43		// utility headers
44		//
45		#include "ParseArgs.h"
46	+	#include "SampleWeight.h"
47		#include "Selection.h"
48		#include "HZZCiCElectronSelection.h"
49	+	#include "HZZLikelihoodElectronSelection.h"
50	+	#include "HZZBDTElectronSelection.h"
51	+	#include "SampleWeight.h"
52
53		//#define THEIR_EVENTS
54
#	Line 70 | Line 74 | int main(int argc, char** argv)
74		}
75		ctrl.dump();
76
77	+	map<string,double> entrymap; // number of unskimmed entries in each file
78	+
79		//
80		// File I/O
81		//--------------------------------------------------------------------------------------------------------------
#	Line 77 | Line 83 | int main(int argc, char** argv)
83		ifstream f(ctrl.inputfile.c_str());
84		string fname;
85		while (f >> fname) {
86	<	if( !(strncmp( fname.c_str(), "#", 1 ) ) ) continue; // skip commented lines
86	>	if( !(strncmp( fname.c_str(), "#", 1 ) ) ) continue; // skip commented lines
87		cout << "adding inputfile : " << fname.c_str() << endl;
88	+	entrymap[string(fname.c_str())] = unskimmedEntries(fname.c_str());
89		chain->AddFile(fname.c_str());
90		}
91
92	+	// table of cross section values
93	+	SimpleTable xstab("./data/xs.dat");
94	+
95		//
96		// Setup
97		//--------------------------------------------------------------------------------------------------------------
98		TH1F * h_evt = new TH1F("hevt", "hevt", 2, 0, 2 );
99	<	TH1F * h_evtfail = new TH1F("hevtfail", "hevtfail", 100, 0, 100 );
99	>	TH1F * h_evtfail = new TH1F("hevtfail", "hevtfail", 1024, 0, 1024 );
100		TNtuple * passtuple = new TNtuple( "passtuple", "passtuple", "run:evt:lumi:channel:mZ1:mZ2:m4l:pt4l:w" );
101		initCiCSelection();
102	+	if(ctrl.eleSele=="lik") initLikSelection();
103	+	if(ctrl.eleSele=="bdt") initBDTSelection();
104		initRunLumiRangeMap();
105
94	–
106		//
107		// Loop
108		//--------------------------------------------------------------------------------------------------------------
#	Line 123 | Line 134 | int main(int argc, char** argv)
134
135
136		cout << "nEntries: " << chain->GetEntries() << endl;
137	<	for(UInt_t ientry=0; ientry<chain->GetEntries(); ientry++) {
138	<
137	>	for(UInt_t ientry=0; ientry<chain->GetEntries(); ientry++) {
138	>
139		chain->GetEntry(ientry);
140	+	if(!(ientry%100000)) cout << "entry: " << ientry << endl;
141	+
142	+	// get event weight for this file
143	+	double eventweight=1;
144	+	if(ctrl.mc) eventweight = getWeight(xstab,entrymap,chain);
145
146		#ifdef THEIR_EVENTS
147		if( !( info->evtNum == 504867308  ||
#	Line 152 | Line 168 | int main(int argc, char** argv)
168		info->evtNum == 107360878  ) ) continue;
169		#endif
170
171	<	unsigned evtfail = fails_HZZ4L_selection(ctrl, info, electronArr, muonArr, passtuple );
156	<	double eventweight = info->eventweight;
171	>	unsigned evtfail = fails_HZZ4L_selection(ctrl, info, electronArr, muonArr, eventweight, passtuple);
172		h_evtfail->Fill( evtfail, eventweight );
173	<	cout << endl << endl;
173	>	if( ctrl.debug ) cout << endl << endl;
174
175		} //end loop over data
176
#	Line 184 | Line 199 | int main(int argc, char** argv)
199		file->Close();
200		delete file;
201
202	+	// map<TString,TMVA::Reader*>::iterator it;
203	+	// for(it=readers.begin(); it!=readers.end(); it++) delete (*it).second;
204	+
205		cerr << "done!" << endl;
206		}
207

Diff Legend

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