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.6 by dasu, Tue Feb 15 21:56:55 2011 UTC

# 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 light leptons only
95 <        pythia.readString("24:onIfAny = 11 13");
94 >        pythia.readString("24:onMode = 0");       // W decays to electron
95 >        pythia.readString("24:onIfAny = 11");
96        }
97 <    else if(strncmp(argv[1], "TLN", 3) == 0)
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 muon
101 >        pythia.readString("24:onIfAny = 13");
102 >      }
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 electron
107 >        pythia.readString("24:onIfAny = 11");
108 >      }
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 light leptons only
113 <        pythia.readString("24:onIfAny = 11 13");
112 >        pythia.readString("24:onMode = 2");       // W+ decays to muon
113 >        pythia.readString("24:onIfAny = 13");
114        }
115 <    else if(strncmp(argv[1], "ZH", 2) == 0)
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
158        {

Diff Legend

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