p p > Q Q~ > d u~ Z W+ 2 2 Enter process to generate (pp > e+ ve jj) Enter maximum number of QCD vertices for each diagram. Enter maximum number of QED vertices for each diagram. Started Sun Mar 13 18:16:26 CST 2005 Running...... Finished Sun Mar 13 18:16:28 CST 2005 Number of unweighted events Enter 1 for parallel 0 for serial run Compiling libraries Using random number seed offset = 0 Working on subprocess: Pddx_duxzw+ 8.000u 0.030s 0:08.06 99.6% 0+0k 0+0io 1605pf+0w ddx_duxzw+ 0.47385E-01 Iteration 1 0.048897 0.000132254515 Iteration 2 0.0485570939 0.000101377277 Iteration 3 0.0484847005 8.00265187E-05 Iteration 4 0.0472443872 6.45865414E-05 1 0.048897 0.002543 0.425308653 2 0.048148 0.0027471 0.515841544 3 0.048353 0.0025067 0.357038239 4 0.044143 0.0021552 1.02678953 Results 0.047385 0.0012485 100000 8.33196896 Updated results in file results.html Pdxd_duxzw+ 7.910u 0.030s 0:07.96 99.7% 0+0k 0+0io 1605pf+0w dxd_duxzw+ 0.36674E-04 ddx_duxzw+ 0.47385E-01 Iteration 1 0.048933646 0.000132155591 Iteration 2 0.0485927266 0.000101302935 Iteration 3 0.0485205975 7.99672994E-05 Iteration 4 0.0472809873 6.45353565E-05 1 0.048933646 0.00254300116 0.425303522 2 0.048182384 0.00274710127 0.515380439 3 0.048389383 0.00250670133 0.356762732 4 0.044182283 0.00215520196 1.02540821 Results 0.047421674 0.0012485007 200000 11.7740779 Updated results in file results.html Pgg_duxzw+ 10.630u 0.070s 0:10.77 99.3% 0+0k 0+0io 3457pf+0w dxd_duxzw+ 0.36674E-04 ddx_duxzw+ 0.47385E-01 gg_duxzw+ 0.46538E-03 Iteration 1 0.049450134 0.000130820033 Iteration 2 0.0490844285 0.000100317963 Iteration 3 0.0489785039 7.92420126E-05 Iteration 4 0.0477449014 6.39165305E-05 1 0.049450134 0.00254343629 0.453217178 2 0.04864341 0.00274751478 0.542634149 3 0.048785112 0.00250684425 0.358286252 4 0.044670576 0.00215566849 1.01542867 Results 0.047887054 0.00124868122 299996 14.2820719 Updated results in file results.html Placing events in Events/events.dat Using random number seed offset = 3 Pddx_duxzw+ dxd_duxzw+ 0.36674E-04 ddx_duxzw+ 0.50972E-01 gg_duxzw+ 0.46538E-03 Iteration 1 0.051552634 1.37707605E-06 Iteration 2 0.0514319138 1.01655256E-06 Iteration 3 0.0512485278 1.0129151E-06 Iteration 4 2.04138752E+10 9.89336968E-20 1 0.051552634 0.000266443048 4260.04153 2 0.05129941 0.000277773083 7831.07391 3 0.000432112 2.6892458E-05 5357.55384 4 2.04138752E+10 4.4940184E-05 3571.70256 Results 0.051474054 0.000239267573 1449996 5.59730389 Updated results in file results.html Pdxd_duxzw+ dxd_duxzw+ 0.36674E-04 ddx_duxzw+ 0.50972E-01 gg_duxzw+ 0.46538E-03 Iteration 1 0.051552634 1.37707605E-06 Iteration 2 0.0514319138 1.01655256E-06 Iteration 3 0.0512485278 1.0129151E-06 Iteration 4 2.04138752E+10 9.89336968E-20 1 0.051552634 0.000266443048 4260.04153 2 0.05129941 0.000277773083 7831.07391 3 0.000432112 2.6892458E-05 5357.55384 4 2.04138752E+10 4.4940184E-05 3571.70256 Results 0.051474054 0.000239267573 1449996 5.59730389 Updated results in file results.html Pgg_duxzw+ dxd_duxzw+ 0.36674E-04 ddx_duxzw+ 0.50972E-01 gg_duxzw+ 0.45850E-03 Iteration 1 0.0344275453 1.39518923E-06 Iteration 2 0.0313000086 1.26672912E-06 1 0.0344275453 0.000219164186 369.744367 2 0.000494232 9.87436683E-06 1794491.7 Results 0.051467174 0.000238443298 1459452 5.59692756 Updated results in file results.html Pddx_duxzw+ Pdxd_duxzw+ Pgg_duxzw+ input the event file (e.g. Events/events.dat) input the banner file name (e.g. banner-events.dat) input the process file (e.g. Process/input.dat) Sorting 38554 Found 38554 events Integrated weight 0.0514638983 Maximum wgt 0.0002382441 1.38675603E-06 Average wgt 1.33485237E-06 Unweight Efficiency 0.962571891 Number of events 38554 0.0002382441 1.38675603E-06 Enter number of events desired ( < 37110) Wrote 19824 events Integrated weight 0.0514638983 Truncated events 17 0.000406908987 0.790668801% c -rw-r--r-- 1 frisch cdf 40k Mar 13 18:32 get.kumac.~1~ -rw-r--r-- 1 frisch cdf 40k Mar 13 18:32 get.kumac drwxr-xr-x 2 frisch cdf 4.0k Mar 13 18:32 CVS lrwxrwxrwx 1 frisch cdf 61 Mar 13 18:34 unweighted_events.dat -> ../../../MadGraph_BJ/QQ300_udW+Z/Events/unweighted_events.dat cdf:MadGraph_Feb6_05/QQ300_udW+Z/Events>>../bin/decay ***************************************************** * DECAY v1.0 * * a MadEvent program * * for decaying unstable particles * * -------------------------------- * * http://madgraph.physics.uiuc.edu * * * * * * * * * * today : 13-Mar-05 * * last update: 07-Nov-03 * ***************************************************** Input run mode: --------------- 0 = calculates decay widths 1 = decay events in file 1 ***************************************************** * run_mode=1 => decaying events * ***************************************************** input event file: (e.g. events.dat) unweighted_events.dat name for output file: (e.g. dec-events.dat) dec-events_W.dat Warning: file input-card.dat is not in the main directory Warning: parameter pdlabel not found setting it to default value cteq6l1 Warning: parameter hgf not found setting it to default value 1.166E-05 Warning: parameter 1/hale not found setting it to default value 128.9 Warning: parameter hss not found setting it to default value 0.2312 Warning: parameter wm not found setting it to default value 80.419 Warning: parameter zm not found setting it to default value 91.188 Warning: parameter lzdecay not found setting it to default value T Warning: parameter zwidth not found setting it to default value 2.44139 Warning: parameter lwdecay not found setting it to default value T Warning: parameter wwidth not found setting it to default value 2.04759 Warning: parameter hmass not found setting it to default value 165. Warning: parameter lhdecay not found setting it to default value F Warning: parameter hwidth not found setting it to default value 0.0037 Warning: parameter tmass not found setting it to default value 174.3 Warning: parameter ltdecay not found setting it to default value T Warning: parameter twidth not found setting it to default value 1.5083 Warning: parameter bmass not found setting it to default value 4.7 Warning: parameter lmass not found setting it to default value 1.777 Warning: parameter mtms not found setting it to default value 164.5 Warning: parameter mbms not found setting it to default value 4.2 Warning: parameter mcms not found setting it to default value 1.42 Warning: parameter iyuk not found setting it to default value 1 Warning: parameter vud not found setting it to default value 0.975 Warning: parameter idef not found setting it to default value 1 Warning: parameter fixed_ren_scale not found setting it to default value T Warning: parameter fixed_fac_scale not found setting it to default value T Warning: parameter scale not found setting it to default value 91.188 Warning: parameter dsqrt(q2fact(1)) not found setting it to default value 91.188 Warning: parameter dsqrt(q2fact(1)) not found setting it to default value 91.188 Warning: parameter lpp(1) not found setting it to default value 1 Warning: parameter lpp(2) not found setting it to default value 1 Warning: parameter ebeam(1) not found setting it to default value 7000. Warning: parameter ebeam(2) not found setting it to default value 7000. ***************************************************** * Input parameters read from banner * ***************************************************** ================================================ | COUPSM called with option 1: | | New MadEvent default "G_mu scheme" | ================================================ EW Params Input Output ----------------------------------------------- GF (10^-5*GeV^-2) 1.1660 * 1.1660 1/alpha 128.9000 132.5513 sin^2(W) 0.2312 0.2222 M_Z (GeV) 91.1880 * 91.1880 M_W (GeV) 80.4190 * 80.4190 Electroweak coupling: --------------------- gw = 0.6531 Boson masses and widths: ------------------------ Z mass = 91.188 GeV Z width = 2.4409 GeV (calc @ LO) W mass = 80.419 GeV W width = 2.0472 GeV (calc @ LO) H mass = 165.000 GeV H width = 0.0037 GeV Fermion masses and widths: -------------------------- top mass = 174.300 GeV top width = 1.5076 GeV (calc @ LO) bottom mass = 4.700 GeV bottom width = 0.0000 GeV charm mass = 0.000 GeV charm width = 0.0000 GeV tau mass = 1.777 GeV tau width = 0.0000 GeV all other quark and lepton masses set to zero Quark running MSbar masses @ 1 loops: -------------------------------------- m_top(hmass)_MSbar = 164.4622 GeV m_bot(hmass)_MSbar = 2.8210 GeV m_cha(hmass)_MSbar = 0.7517 GeV Z,W LO partial widths and effective BRs: ------------------------------------- width Z->nn = 0.1658 GeV, BR(Z->nn) = 0.2038 width Z->ll = 0.0839 GeV, BR(Z->ll) = 0.0688 width Z->TT = 0.0837 GeV, BR(Z->TT) = 0.0343 width Z->uu = 0.2900 GeV, BR(Z->uu) = 0.1188 width Z->dd = 0.3719 GeV, BR(Z->dd) = 0.3047 width Z->cc = 0.2900 GeV, BR(Z->cc) = 0.1188 width Z->bb = 0.3680 GeV, BR(Z->bb) = 0.1507 width W->nl = 0.2275 GeV, BR(W->nl) = 0.2222 width W->nT = 0.2273 GeV, BR(W->nT) = 0.1110 width W->ud = 0.6825 GeV, BR(W->ud) = 0.3334 width W->cs = 0.6825 GeV, BR(W->cs) = 0.3334 Higgs total width = 0.00370 GeV (set by hand) Higgs LO partial widths and BRs: -------------------------------- width H->tt~ = 0.00000000 GeV, BR_LO(H->tt~) = 0.0000 width H->bb~ = 0.00257203 GeV, BR_LO(H->bb~) = 0.6951 width H->cc~ = 0.00018350 GeV, BR_LO(H->cc~) = 0.0496 width H->ta ta~ = 0.00034161 GeV, BR_LO(H->ta ta~) = 0.0923 width H->zz = 0.00000000 GeV, BR_LO(H->zz) = 0.0000 width H->w- w+ = 0.23909523 GeV, BR_LO(H->w- w+) = ****** Boson couplings: ---------------- gwwa = 0.30790 gwwz = 0.57599 gwwh = 52.52355 0.00000 gzzh = 67.53239 0.00000 ghhh = -331.66215 0.00000 gwwhh = 0.21329 0.00000 gzzhh = 0.27423 0.00000 ghhhh = -1.34680 0.00000 FFV couplings: -------------- gal(L) = 0.3079023 0.0000000 gal(R) = 0.3079023 0.0000000 gau(L) = -0.2052682 0.0000000 gau(R) = -0.2052682 0.0000000 gad(L) = 0.1026341 0.0000000 gad(R) = 0.1026341 0.0000000 gwf(L) = -0.4618282 0.0000000 gwf(R) = 0.0000000 0.0000000 gzn(L) = -0.3702921 0.0000000 gzn(R) = 0.0000000 0.0000000 gzl(L) = 0.2056998 0.0000000 gzl(R) = -0.1645923 0.0000000 gzu(L) = -0.2605639 0.0000000 gzu(R) = 0.1097282 0.0000000 gzd(L) = 0.3154280 0.0000000 gzd(R) = -0.0548641 0.0000000 FFH couplings: -------------- gHtop(L) = -0.6678408, 0.0000000 gHtop(R) = -0.6678408, 0.0000000 gHbot(L) = -0.0114555, 0.0000000 gHbot(R) = -0.0114555, 0.0000000 gHcha(L) = -0.0030523, 0.0000000 gHcha(R) = -0.0030523, 0.0000000 gHtau(L) = -0.0072160, 0.0000000 gHtau(R) = -0.0072160, 0.0000000 Strong couplings: ----------------- gg(1) = -1.2781346, 0.0000000 gg(2) = -1.2781346, 0.0000000 Collider parameters: -------------------- Running at P P machine @ 14000. GeV PDF set = cteq6l1 alpha_s(Mz)= 0.1300 running at 1 loops. Renormalization scale fixed @ 91.188 Factorization scales fixed @ 91.188 91.188 ***************************************************** * seed not present in banner * * no iseed.dat => iseed=-1 * ***************************************************** ***************************************************** * first rnd number= 0.390043 * ***************************************************** Implemented decays are for: --------------------------- Leptons: ta- ta+ Quarks : t t~ Bosons : z w+ w- h Input particle to be decayed (e.g. t~): w+ particle to decay :w+ Implemented decay modes: ------------------------ 1 w+ -> e+ ve 2 w+ -> mu+ vm 3 w+ -> ta+ vt 4 w+ -> l+ vl (e+mu) 5 w+ -> l+ vl (e+mu+ta) 6 w+ -> j j (ud+cs) 7 w+ -> anything (e+mu+ta+ud+cs) your choice is: 7 ---------------------------------------------- Relevant branching ratios ---------------------------------------------- BR (w ->anything) = 1 (=e+mu+ta+ud+cs) ---------------------------------------------- Decay information ---------------------------------------------- particle name : w+ decay mode : w+ -> anything (e+mu+ta+ud+cs) ---------------------------------------------- MC Evaluation of the Partial Width ---------------------------------------------- input param. for vegas: ndim= 2 calls= 96800. it= 0 itmx= 5 nprn= 0 alph= 1.50 mds= -1 nd= 44 iter. no. 1: integral = 2.039774 +/- 0.98E-02 iter. no. 2: integral = 2.054559 +/- 0.97E-02 iter. no. 3: integral = 2.049870 +/- 0.97E-02 iter. no. 4: integral = 2.046835 +/- 0.98E-02 iter. no. 5: integral = 2.056194 +/- 0.97E-02 input param. for vegas: ndim= 2 calls= 490050. it= 6 itmx= 1 nprn= 0 alph= 1.50 mds= -1 nd= 45 iter. no. 1: integral = 2.047368 +/- 0.43E-02 LO partial width = 2.04722511 GeV MC partial width = 2.04736781 +- 0.00431996398 GeV Max. wgt = 14.0728418 Unwgt. Eff.(%) = 14.5 Wait....Reading In Event File ---------------------------------------------- Input events information ---------------------------------------------- Input Event file : unweighted_events.dat Number of Events : 19824 Integrated weight (pb) : 0.051463897 Max wgt : 2.59604E-06 Average wgt : 2.59604E-06 Events with 1 w+ : 19824 ---------------------------------------------- Decay events in file : in progress ---------------------------------------------- 10 % done 20 % done 30 % done 40 % done 50 % done 60 % done 70 % done 80 % done 90 % done 100% done ---------------------------------------------- Statistics of the unweighting ---------------------------------------------- Events weighted = 137127 Events unweighted = 19824 Width (wgt) = 2.0485986 GeV Width (unwgt) = 2.03446451 GeV Events over = 14 Integral over = 2.27920067E-05 GeV Wait....Writing Out Decayed Event file ---------------------------------------------- output events information ---------------------------------------------- Output Event file : dec-events_W.dat Number of Events : 19824 Integrated weight (pb) : 0.051463897 Max wgt : 2.59604E-06 Average wgt : 2.59604E-06 cdf:MadGraph_Feb6_05/QQ300_udW+Z/Events>>edit decayW.log &