p p > Q Q~ > d u~ H 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 6 15:05:09 CST 2005 Running...... Finished Sun Mar 6 15:05:10 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_duxhw+ 8.080u 0.000s 0:08.09 99.8% 0+0k 0+0io 1605pf+0w ddx_duxhw+ 0.23886E-01 Iteration 1 0.024792 3.88356047E-05 Iteration 2 0.0244955087 3.67173239E-05 Iteration 3 0.0241648264 2.70060834E-05 Iteration 4 0.0238720626 2.12153721E-05 1 0.024792 0.00098123 0.996513732 2 0.023862 0.0013805 0.996541923 3 0.023656 0.00095855 0.525748611 4 0.023231 0.00087423 0.57526966 Results 0.023886 0.00053336 100000 7.06117534 Updated results in file results.html Pdxd_duxhw+ 8.110u 0.020s 0:08.17 99.5% 0+0k 0+0io 1603pf+0w dxd_duxhw+ 0.19990E-04 ddx_duxhw+ 0.23886E-01 Iteration 1 0.024810337 3.88071219E-05 Iteration 2 0.0245147036 3.66894765E-05 Iteration 3 0.0241846156 2.69842295E-05 Iteration 4 0.0238917854 2.11979712E-05 1 0.024810337 0.000981232781 0.993326405 2 0.02388322 0.00138050048 0.993346833 3 0.023676797 0.000958551678 0.523881086 4 0.023250605 0.00087423068 0.574113531 Results 0.02390599 0.000533360656 200000 9.97767233 Updated results in file results.html Pgg_duxhw+ 10.710u 0.040s 0:10.85 99.0% 0+0k 0+0io 3445pf+0w dxd_duxhw+ 0.19990E-04 ddx_duxhw+ 0.23886E-01 gg_duxhw+ 0.23214E-03 Iteration 1 0.025040466 3.84594696E-05 Iteration 2 0.0247468563 3.63557796E-05 Iteration 3 0.0244203983 2.67298165E-05 Iteration 4 0.0241224519 2.10007808E-05 1 0.025040466 0.000981347564 0.991163894 2 0.024120587 0.0013805865 0.991164863 3 0.023918847 0.000958693903 0.520003253 4 0.023469621 0.000874325276 0.579613865 Results 0.02413813 0.000533413437 299996 12.1036981 Updated results in file results.html Placing events in Events/events.dat Using random number seed offset = 3 Pddx_duxhw+ dxd_duxhw+ 0.19990E-04 ddx_duxhw+ 0.25305E-01 gg_duxhw+ 0.23214E-03 Iteration 1 0.0255615265 3.02922906E-07 Iteration 2 0.0255567474 3.13946266E-07 Iteration 3 -1.97045759 4.19262962E-09 1 0.0255615265 8.79952947E-05 277028.932 2 0.025543621 0.000145732701 313801.032 3 -1.97094636 9.0907E-05 156943.867 Results 0.02555713 6.45262404E-05 949996 2.46085011 Updated results in file results.html Pdxd_duxhw+ dxd_duxhw+ 0.19990E-04 ddx_duxhw+ 0.25305E-01 gg_duxhw+ 0.23214E-03 Iteration 1 0.0255615265 3.02922906E-07 Iteration 2 0.0255567474 3.13946266E-07 Iteration 3 -1.96960288 4.19444782E-09 1 0.0255615265 8.79952947E-05 277031.832 2 0.025543621 0.000145732701 313804.317 3 -1.97009187 9.0907E-05 156945.547 Results 0.02555713 6.45262404E-05 949996 2.46085011 Updated results in file results.html Pgg_duxhw+ dxd_duxhw+ 0.19990E-04 ddx_duxhw+ 0.25305E-01 gg_duxhw+ 0.22913E-03 Iteration 1 0.0255584434 5.62256549E-07 Iteration 2 0.0255400968 2.91079849E-07 1 0.0255584434 0.000119876612 0.0171585061 2 0.025529493 9.10325068E-05 0.0343559497 Results 0.02555412 6.41425898E-05 948098 2.44406178 Updated results in file results.html Pddx_duxhw+ Pdxd_duxhw+ Pgg_duxhw+ 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 30910 Found 30910 events Integrated weight 0.0255527414 Maximum wgt 5.723355E-05 8.61352817E-07 Average wgt 8.26682025E-07 Unweight Efficiency 0.959748444 Number of events 30910 5.723355E-05 8.61352817E-07 Enter number of events desired ( < 29665) Wrote 19899 events Integrated weight 0.0255527414 Truncated events 21 0.000104103387 0.407405941% cdf:MadGraph_Feb6_05/QQ300_W+H165/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_H.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~): h particle to decay :h Implemented decay modes: ------------------------ 1 h -> b b~ 2 h -> ta- ta+ 3 h -> mu- mu+ 4 h -> c c~ 5 h -> t t~ (when m_h>2*m_t) 6 h -> g g 7 h -> a a 8 h -> z a (when m_h> m_z) 9 h -> w+ w- (when m_h>2*m_w) 10 h -> w* w -> l vl l' vl' (l,l'=e,mu) 11 h -> w* w -> l vl l' vl' (l,l'=e,mu,ta) 12 h -> w* w -> j j l vl (jj=ud,cs;l=e,mu) 13 h -> w* w -> j j l vl (jj=ud,cs;l=e,mu,ta) 14 h -> z z (when m_h>2*m_z) 15 h -> z* z -> l- l+ l-' l+'(l,l'=e,mu) 16 h -> z* z -> l- l+ l-' l+'(l,l'=e,mu,ta) 17 h -> z* z -> j j~ l- l+ (j=u,d,c,s;l=e,mu ) 18 h -> z* z -> j j~ l- l+ (j=u,d,c,s;l=e,mu,ta) 19 h -> z* z -> b b~ l- l+ (l=e,mu ) 20 h -> z* z -> b b~ l- l+ (l=e,mu,ta ) 21 h -> z* z -> vl vl~ l-' l+'(l=e,mu,ta;l'=e,mu) 22 h -> z* z -> vl vl~ l-' l+'(l,l'=e,mu,ta) your choice is: 9 ---------------------------------------------- Relevant branching ratios ---------------------------------------------- Higgs mass = 165.00000 GeV Higgs tot width = 0.00000 GeV BR (h -> w+ w- ) = 0.00000 Best tot Br Ratio = 0.00000 ---------------------------------------------- Decay information ---------------------------------------------- particle name : h decay mode : h -> w+ w- (when m_h>2*m_w) ---------------------------------------------- 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 = 0.2381594 +/- 0.91E-03 iter. no. 2: integral = 0.2383139 +/- 0.90E-03 iter. no. 3: integral = 0.2386205 +/- 0.91E-03 iter. no. 4: integral = 0.2382190 +/- 0.91E-03 iter. no. 5: integral = 0.2383377 +/- 0.91E-03 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 = 0.2383720 +/- 0.40E-03 LO partial width = 0.23909523 GeV MC partial width = 0.238371998 +- 0.000403452286 GeV Max. wgt = 0.860983896 Unwgt. Eff.(%) = 27.7 Wait....Reading In Event File ---------------------------------------------- Input events information ---------------------------------------------- Input Event file : unweighted_events.dat Number of Events : 19899 Integrated weight (pb) : 0.0255527437 Max wgt : 1.284122E-06 Average wgt : 1.284122E-06 Events with 1 h : 19899 ---------------------------------------------- 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 = 71765 Events unweighted = 19899 Width (wgt) = 0.239313931 GeV Width (unwgt) = 0.238733624 GeV Events over = 0 Integral over = 0. GeV Wait....Writing Out Decayed Event file ---------------------------------------------- output events information ---------------------------------------------- Output Event file : dec-events_H.dat Number of Events : 19899 Integrated weight (pb) : 0. Max wgt : 0. Average wgt : 0. cdf:MadGraph_Feb6_05/QQ300_W+H165/Events>>edit decay.log &