14 |
|
using namespace Pythia8; |
15 |
|
|
16 |
|
#include "UltraFastSim.h" |
17 |
< |
#include "Analysis.h" |
17 |
> |
#include "bbHAnalysis.h" |
18 |
|
|
19 |
|
int main(int argc, char **argv) { |
20 |
|
// Generator. Process selection. LHC initialization. Histogram. |
178 |
|
// Ultra Fast Simulator |
179 |
|
|
180 |
|
UltraFastSim ufs(rndmPtr); |
181 |
< |
Analysis analysis(argv[1], &ufs); |
181 |
> |
bbHAnalysis analysis(argv[1], &pythia, &ufs, true); |
182 |
|
|
183 |
|
// Begin event loop |
184 |
|
for (int iEvent = 0; iEvent < nEvents; ) { |
215 |
|
exit(1); |
216 |
|
} |
217 |
|
|
218 |
– |
cout << "Number of Gen Chrgd = " << ufs.chargedHadronList().size() << endl; |
219 |
– |
cout << "Number of Gen Neutr = " << ufs.neutralHadronList().size() << endl; |
220 |
– |
cout << "Number of Gen Phtns = " << ufs.photonList().size() << endl; |
221 |
– |
cout << "Number of Gen Elcns = " << ufs.electronList().size() << endl; |
222 |
– |
cout << "Number of Gen Muons = " << ufs.muonList().size() << endl; |
223 |
– |
cout << "Number of Gen Taus = " << ufs.genTauList().size() << endl; |
224 |
– |
cout << "Number of taus = " << ufs.tauList().size() << endl; |
225 |
– |
cout << "Number of c Quarks = " << ufs.cQuarkList().size() << endl; |
226 |
– |
cout << "Number of b Quarks = " << ufs.bQuarkList().size() << endl; |
227 |
– |
cout << "Number of Jets = " << ufs.jetList().size() << endl; |
228 |
– |
cout << "Number of bJets = " << ufs.bJetList().size() << endl; |
229 |
– |
cout << endl; |
230 |
– |
|
218 |
|
// Run analysis |
219 |
|
|
220 |
|
if(!analysis.run()) |