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.2 by dasu, Wed Feb 9 07:20:36 2011 UTC vs.
Revision 1.4 by dasu, Thu Feb 10 14:29:41 2011 UTC

# Line 23 | Line 23 | int main(int argc, char **argv) {
23    int meanPileupEventCount = 0;
24    if(argc < 3 || argc > 5)
25      {
26 +      cerr << "Command syntax: " << argv[0] << " BBH[...]" << " <nEvents> [RandomSeed] [meanPileupEventCount]" << endl;
27        cerr << "Command syntax: " << argv[0] << " BBA[...]" << " <nEvents> [RandomSeed] [meanPileupEventCount]" << endl;
28        cerr << "or              " << argv[0] << " ZLL[...]" << " <nEvents> [RandomSeed] [meanPileupEventCount]" << endl;
29 +      cerr << "or              " << argv[0] << " ZMM[...]" << " <nEvents> [RandomSeed] [meanPileupEventCount]" << endl;
30        cerr << "or              " << argv[0] << " TLN[...]" << " <nEvents> [RandomSeed] [meanPileupEventCount]" << endl;
31        cerr << "or              " << argv[0] << " WLN[...]" << " <nEvents> [RandomSeed] [meanPileupEventCount]" << endl;
32        exit(1);
# Line 43 | Line 45 | int main(int argc, char **argv) {
45          pythia.readString("15:onMode = 2");       // tau- decays to electrons and muons - second tau decays in all modes
46          pythia.readString("15:onIfAny = 11 13");
47        }
48 +    else if(strncmp(argv[1], "ZMM", 3) == 0)
49 +      {
50 +        pythia.readString("WeakSingleBoson:ffbar2gmZ = on");
51 +        pythia.readString("PhaseSpace:mHatMin = 60.");    
52 +        pythia.readString("PhaseSpace:mHatMax = -1");    
53 +        pythia.readString("23:onMode = 0");       // Z decays to muons only
54 +        pythia.readString("23:onIfAny = 13");
55 +      }
56      else if(strncmp(argv[1], "WLN", 3) == 0)
57        {
58          pythia.readString("WeakSingleBoson:ffbar2W = on");
# Line 55 | Line 65 | int main(int argc, char **argv) {
65          pythia.readString("24:onMode = 2");       // W+ decays to light leptons only
66          pythia.readString("24:onIfAny = 11 13");
67        }
68 <    else if(strncmp(argv[1], "BBA", 3) == 0)
68 >    else if(strncmp(argv[1], "ZH", 2) == 0)
69 >      {
70 >        pythia.readString("HiggsSM:ffbar2HZ = on");
71 >        pythia.readString("25:m0 = 120");
72 >        pythia.readString("25:onMode = 0");      // Higgs decays to bbBar only
73 >        pythia.readString("25:onIfAny = 5");
74 >        pythia.readString("23:onMode = 0");      // Z decays to mu+,mu- only (for trigger)
75 >        pythia.readString("23:onIfAny = 13");
76 >      }
77 >    else if(strncmp(argv[1], "BB", 2) == 0)
78        {
79          pythia.readString("Higgs:useBSM = on");
80 <        pythia.readString("HiggsBSM:gg2A3bbbar = on");
80 >        if(strncmp(argv[1], "BBH", 3) == 0)
81 >          pythia.readString("HiggsBSM:gg2H2bbbar = on");
82 >        else
83 >          pythia.readString("HiggsBSM:gg2A3bbbar = on");
84 >        pythia.readString("25:m0 = 129");
85 >        pythia.readString("35:m0 = 499.7");
86 >        pythia.readString("36:m0 = 500");
87 >        pythia.readString("37:m0 = 506");
88 >        pythia.readString("HiggsHchg:tanBeta = 30");
89 >        pythia.readString("35:onMode = 0");      // H2(H_2) decays only to tau-pairs
90 >        pythia.readString("35:onIfAny = 15");
91          pythia.readString("36:onMode = 0");      // A0(H_3) decays only to tau-pairs
92          pythia.readString("36:onIfAny = 15");
93        }
# Line 93 | Line 122 | int main(int argc, char **argv) {
122      // Generate event. Skip if error. List first one.
123      if (!pythia.next()) continue;
124  
125 +    if(iEvent < 10) pythia.event.list();
126 +
127      // Add pileup
128  
129      int pileupEventCount = 0;

Diff Legend

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