ViewVC Help

View File | Revision Log | Show Annotations | Root Listing

root/cvsroot/UserCode/dasu/UltraFastSim/generateData.cc

Comparing UserCode/dasu/UltraFastSim/generateData.cc (file contents):
Revision 1.12 by mmhl, Thu Feb 24 20:17:12 2011 UTC vs.
Revision 1.15 by dasu, Mon Feb 28 11:13:18 2011 UTC

#	Line 13 | Line 13 | using namespace std;
13		#include "Pythia.h"
14		using namespace Pythia8;
15
16	–	#include "TROOT.h"
17	–	#include "TFile.h"
18	–
16		#include "UltraFastSim.h"
17	<	#include "bbHAnalysis.h"
21	<	//#include "ZHAnalysis.h"
22	<	#include "bbHAnalysis.cc"
23	<	//#include "ZHAnalysis.cc"
17	>	#include "UFSDataStore.h"
18
19		int main(int argc, char **argv) {
20		// Generator. Process selection. LHC initialization. Histogram.
#	Line 174 | Line 168 | int main(int argc, char **argv) {
168		}
169		}
170
177	–	string name(argv[1]);
178	–	name += ".root";
179	–	TFile *outFile = TFile::Open(name.c_str(), "recreate");
180	–
171		// Initialize pythia
172
173		pythia.init( 2212, 2212, 14000.);
#	Line 194 | Line 184 | int main(int argc, char **argv) {
184
185		// Ultra Fast Simulator
186
187	<	UltraFastSim ufs(rndmPtr);
188	<	bbHAnalysis bbH(outFile, &pythia, &ufs, true);
189	<	//ZHAnalysis ZH; // outFile, &pythia, &ufs, true);
190	<	//ZH.BookTree();
187	>	UltraFastSim ufs;
188	>	char jobName[256];
189	>	sprintf(jobName, "%s-%4.4d", argv[1], runNumber);
190	>	UFSDataStore dataStore(jobName, &ufs);
191
192		// Begin event loop
193		for (int iEvent = 0; iEvent < nEvents; ) {
#	Line 234 | Line 224 | int main(int argc, char **argv) {
224
225		// Ultra fast simulation
226
227	<	if(!ufs.run(pythia.event))
227	>	if(!ufs.run(pythia.event, rndmPtr))
228		{
229		cerr << "Ultra fast simulation failed - aborting" << endl;
230		exit(1);
231		}
232
233	<	// Run bbHAnalysis
233	>	// Store data
234
235	<	if(!bbH.run())
235	>	if(!dataStore.run())
236		{
237	<	cerr << "bbH Analysis failed - aborting" << endl;
248	<	exit(2);
249	<	}
250	<
251	<	// Run ZHAnalysis
252	<
253	<	/* if(!ZH.Run(&ufs))
254	<	{
255	<	cerr << "ZH Analysis failed - aborting" << endl;
237	>	cerr << "Failed to store data" << endl;
238		exit(3);
239		}
240	<	*/
240	>
241		// End of event loop. Statistics. Histogram. Done.
242
243		if(!(iEvent % 100)) cout << "Processed event " << iEvent << endl;
#	Line 268 | Line 250 | int main(int argc, char **argv) {
250		pythia.statistics();
251		pileupPythia.statistics();
252
271	–	bbH.end();
272	–
273	–	outFile->cd();
274	–	outFile->Write();
275	–	outFile->Close();
276	–
253		return 0;
254
255		}

Diff Legend

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