ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/dasu/UltraFastSim/generateData.cc
(Generate patch)

Comparing UserCode/dasu/UltraFastSim/generateData.cc (file contents):
Revision 1.5 by dasu, Tue Feb 15 06:49:11 2011 UTC vs.
Revision 1.9 by dasu, Thu Feb 17 14:57:43 2011 UTC

# Line 14 | Line 14 | using namespace std;
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.
# Line 24 | Line 24 | int main(int argc, char **argv) {
24    int meanPileupEventCount = 0;
25    if(argc < 3 || argc > 5)
26      {
27 <      cerr << "Command syntax: " << argv[0] << " BBH[...]" << " <nEvents> [RandomSeed] [meanPileupEventCount]" << endl;
28 <      cerr << "Command syntax: " << argv[0] << " BBA[...]" << " <nEvents> [RandomSeed] [meanPileupEventCount]" << endl;
29 <      cerr << "or              " << argv[0] << " ZLL[...]" << " <nEvents> [RandomSeed] [meanPileupEventCount]" << endl;
30 <      cerr << "or              " << argv[0] << " ZMM[...]" << " <nEvents> [RandomSeed] [meanPileupEventCount]" << endl;
31 <      cerr << "or              " << argv[0] << " TLN[...]" << " <nEvents> [RandomSeed] [meanPileupEventCount]" << endl;
32 <      cerr << "or              " << argv[0] << " WLN[...]" << " <nEvents> [RandomSeed] [meanPileupEventCount]" << endl;
27 >      cerr << "Command syntax: " << argv[0] << " ZHmmbb[...]" << " <nEvents> [RandomSeed] [meanPileupEventCount]" << endl;
28 >      cerr << "or              " << argv[0] << " ZHeebb[...]" << " <nEvents> [RandomSeed] [meanPileupEventCount]" << endl;
29 >      cerr << "or              " << argv[0] << " BBHmh[...] " << " <nEvents> [RandomSeed] [meanPileupEventCount]" << endl;
30 >      cerr << "or              " << argv[0] << " BBAmh[...] " << " <nEvents> [RandomSeed] [meanPileupEventCount]" << endl;
31 >      cerr << "or              " << argv[0] << " Zee[...]   " << " <nEvents> [RandomSeed] [meanPileupEventCount]" << endl;
32 >      cerr << "or              " << argv[0] << " Zmm[...]   " << " <nEvents> [RandomSeed] [meanPileupEventCount]" << endl;
33 >      cerr << "or              " << argv[0] << " Zeh[...]   " << " <nEvents> [RandomSeed] [meanPileupEventCount]" << endl;
34 >      cerr << "or              " << argv[0] << " Zmh[...]   " << " <nEvents> [RandomSeed] [meanPileupEventCount]" << endl;
35 >      cerr << "or              " << argv[0] << " Zhh[...]   " << " <nEvents> [RandomSeed] [meanPileupEventCount]" << endl;
36 >      cerr << "or              " << argv[0] << " Wen[...]   " << " <nEvents> [RandomSeed] [meanPileupEventCount]" << endl;
37 >      cerr << "or              " << argv[0] << " Wmn[...]   " << " <nEvents> [RandomSeed] [meanPileupEventCount]" << endl;
38 >      cerr << "or              " << argv[0] << " Ten[...]   " << " <nEvents> [RandomSeed] [meanPileupEventCount]" << endl;
39 >      cerr << "or              " << argv[0] << " Tmn[...]   " << " <nEvents> [RandomSeed] [meanPileupEventCount]" << endl;
40 >      cerr << "In above e-electron, m-muon, h-hadronic tau. n-neutrino" << endl;
41        exit(1);
42      }
43    else {
44      nEvents = atoi(argv[2]);
45      if(argc >= 4) runNumber = atoi(argv[3]);
46      if(argc >= 5) meanPileupEventCount = atoi(argv[4]);
47 <    if(strncmp(argv[1], "ZLL", 3) == 0)
47 >    if(strncmp(argv[1], "Zee", 3) == 0)
48        {
49          pythia.readString("WeakSingleBoson:ffbar2gmZ = on");
50          pythia.readString("PhaseSpace:mHatMin = 60.");    
51          pythia.readString("PhaseSpace:mHatMax = -1");    
52 <        pythia.readString("23:onMode = 0");       // Z decays to leptons only
53 <        pythia.readString("23:onIfAny = 11 13 15");
46 <        pythia.readString("15:onMode = 2");       // tau- decays to electrons and muons - second tau decays in all modes
47 <        pythia.readString("15:onIfAny = 11 13");
52 >        pythia.readString("23:onMode = 0");       // Z decays to electrons only
53 >        pythia.readString("23:onIfAny = 11");
54        }
55 <    else if(strncmp(argv[1], "ZMM", 3) == 0)
55 >    else if(strncmp(argv[1], "Zmm", 3) == 0)
56        {
57          pythia.readString("WeakSingleBoson:ffbar2gmZ = on");
58          pythia.readString("PhaseSpace:mHatMin = 60.");    
# Line 54 | Line 60 | int main(int argc, char **argv) {
60          pythia.readString("23:onMode = 0");       // Z decays to muons only
61          pythia.readString("23:onIfAny = 13");
62        }
63 <    else if(strncmp(argv[1], "WLN", 3) == 0)
63 >    else if(strncmp(argv[1], "Zeh", 3) == 0)
64 >      {
65 >        pythia.readString("WeakSingleBoson:ffbar2gmZ = on");
66 >        pythia.readString("PhaseSpace:mHatMin = 60.");    
67 >        pythia.readString("PhaseSpace:mHatMax = -1");    
68 >        pythia.readString("23:onMode = 0");       // Z decays to taus only
69 >        pythia.readString("23:onIfAny = 15");
70 >        pythia.readString("15:onMode = 2");       // tau- decays to electrons - second tau decays in all modes
71 >        pythia.readString("15:onIfAny = 11");
72 >      }
73 >    else if(strncmp(argv[1], "Zmh", 3) == 0)
74 >      {
75 >        pythia.readString("WeakSingleBoson:ffbar2gmZ = on");
76 >        pythia.readString("PhaseSpace:mHatMin = 60.");    
77 >        pythia.readString("PhaseSpace:mHatMax = -1");    
78 >        pythia.readString("23:onMode = 0");       // Z decays to taus only
79 >        pythia.readString("23:onIfAny = 15");
80 >        pythia.readString("15:onMode = 2");       // tau- decays to muons - second tau decays in all modes
81 >        pythia.readString("15:onIfAny = 13");
82 >      }
83 >    else if(strncmp(argv[1], "Zhh", 3) == 0)
84 >      {
85 >        pythia.readString("WeakSingleBoson:ffbar2gmZ = on");
86 >        pythia.readString("PhaseSpace:mHatMin = 60.");    
87 >        pythia.readString("PhaseSpace:mHatMax = -1");    
88 >        pythia.readString("23:onMode = 0");       // Z decays to taus only
89 >        pythia.readString("23:onIfAny = 15");
90 >      }
91 >    else if(strncmp(argv[1], "Wen", 3) == 0)
92 >      {
93 >        pythia.readString("WeakSingleBoson:ffbar2W = on");
94 >        pythia.readString("24:onMode = 0");       // W decays to electron
95 >        pythia.readString("24:onIfAny = 11");
96 >      }
97 >    else if(strncmp(argv[1], "Wmn", 3) == 0)
98        {
99          pythia.readString("WeakSingleBoson:ffbar2W = on");
100 <        pythia.readString("24:onMode = 0");       // W decays to light leptons only
101 <        pythia.readString("24:onIfAny = 11 13");
100 >        pythia.readString("24:onMode = 0");       // W decays to muon
101 >        pythia.readString("24:onIfAny = 13");
102        }
103 <    else if(strncmp(argv[1], "TLN", 3) == 0)
103 >    else if(strncmp(argv[1], "Ten", 3) == 0)
104        {
105          pythia.readString("Top:all = on");
106 <        pythia.readString("24:onMode = 2");       // W+ decays to light leptons only
107 <        pythia.readString("24:onIfAny = 11 13");
106 >        pythia.readString("24:onMode = 2");       // W+ decays to electron
107 >        pythia.readString("24:onIfAny = 11");
108        }
109 <    else if(strncmp(argv[1], "ZH", 2) == 0)
109 >    else if(strncmp(argv[1], "Tmn", 3) == 0)
110 >      {
111 >        pythia.readString("Top:all = on");
112 >        pythia.readString("24:onMode = 2");       // W+ decays to muon
113 >        pythia.readString("24:onIfAny = 13");
114 >      }
115 >    else if(strncmp(argv[1], "ZHmmbb", 6) == 0)
116        {
117          pythia.readString("HiggsSM:ffbar2HZ = on");
118          pythia.readString("25:m0 = 120");
# Line 75 | Line 121 | int main(int argc, char **argv) {
121          pythia.readString("23:onMode = 0");      // Z decays to mu+,mu- only (for trigger)
122          pythia.readString("23:onIfAny = 13");
123        }
124 +    else if(strncmp(argv[1], "ZHeebb", 6) == 0)
125 +      {
126 +        pythia.readString("HiggsSM:ffbar2HZ = on");
127 +        pythia.readString("25:m0 = 120");
128 +        pythia.readString("25:onMode = 0");      // Higgs decays to bbBar only
129 +        pythia.readString("25:onIfAny = 5");
130 +        pythia.readString("23:onMode = 0");      // Z decays to e+,e- only (for trigger)
131 +        pythia.readString("23:onIfAny = 11");
132 +      }
133      else if(strncmp(argv[1], "BB", 2) == 0)
134        {
135          pythia.readString("Higgs:useBSM = on");
81 –        if(strncmp(argv[1], "BBH", 3) == 0)
82 –          pythia.readString("HiggsBSM:gg2H2bbbar = on");
83 –        else
84 –          pythia.readString("HiggsBSM:gg2A3bbbar = on");
136          pythia.readString("25:m0 = 129");
137          pythia.readString("35:m0 = 499.7");
138          pythia.readString("36:m0 = 500");
139          pythia.readString("37:m0 = 506");
140          pythia.readString("HiggsHchg:tanBeta = 30");
141 <        pythia.readString("35:onMode = 0");      // H2(H_2) decays only to tau-pairs
142 <        pythia.readString("35:onIfAny = 15");
143 <        pythia.readString("36:onMode = 0");      // A0(H_3) decays only to tau-pairs
144 <        pythia.readString("36:onIfAny = 15");
141 >        if(strncmp(argv[1], "BBH", 3) == 0) {
142 >          pythia.readString("HiggsBSM:gg2H2bbbar = on");
143 >          pythia.readString("35:onMode = 0");      // H2(H_2) decays only to tau-pairs
144 >          pythia.readString("35:onIfAny = 15");
145 >        }
146 >        else {
147 >          pythia.readString("HiggsBSM:gg2A3bbbar = on");
148 >          pythia.readString("36:onMode = 0");      // A0(H_3) decays only to tau-pairs
149 >          pythia.readString("36:onIfAny = 15");
150 >        }
151 >        if(strncmp(argv[1], "BBHmh", 5) == 0 ||
152 >           strncmp(argv[1], "BBAmh", 5) == 0) {
153 >          pythia.readString("15:onMode = 2");       // tau- decays to muons - second tau decays in all modes
154 >          pythia.readString("15:onIfAny = 13");
155 >        }
156        }
157 +     else if(strncmp(argv[1], "ZZmmbb", 6) == 0)
158 +       {
159 +        pythia.readString("WeakDoubleBoson:ffbar2gmZgmZ = on");
160 +        pythia.readString("23:onMode = 2");      
161 +        pythia.readString("23:onIfAny = 13");
162 +        pythia.readString("23:onMode = 3");
163 +        pythia.readString("23:onIfAny = 5");
164 +       }
165      else
166        {
167          cerr << "Unknown process type " << argv[1] << " - aborting" << endl;
# Line 116 | Line 186 | int main(int argc, char **argv) {
186    // Ultra Fast Simulator
187  
188    UltraFastSim ufs(rndmPtr);
189 <  Analysis analysis(argv[1], &ufs);
189 >  bbHAnalysis analysis(argv[1], &pythia, &ufs, true);
190  
191    // Begin event loop
192    for (int iEvent = 0; iEvent < nEvents; ) {
# Line 153 | Line 223 | int main(int argc, char **argv) {
223          exit(1);
224        }
225  
156 –    cout << "Number of Gen Chrgd = " << ufs.chargedHadronList().size() << endl;
157 –    cout << "Number of Gen Neutr = " << ufs.neutralHadronList().size() << endl;
158 –    cout << "Number of Gen Phtns = " << ufs.photonList().size() << endl;
159 –    cout << "Number of Gen Elcns = " << ufs.electronList().size() << endl;
160 –    cout << "Number of Gen Muons = " << ufs.muonList().size() << endl;
161 –    cout << "Number of Gen Taus =  " << ufs.genTauList().size() << endl;
162 –    cout << "Number of taus =      " << ufs.tauList().size() << endl;
163 –    cout << "Number of c Quarks =  " << ufs.cQuarkList().size() << endl;
164 –    cout << "Number of b Quarks =  " << ufs.bQuarkList().size() << endl;
165 –    cout << "Number of Jets =      " << ufs.jetList().size() << endl;
166 –    cout << "Number of bJets =     " << ufs.bJetList().size() << endl;
167 –    cout << endl;
168 –
226      // Run analysis
227  
228      if(!analysis.run())

Diff Legend

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