************** MCFM - version 5.1 **************** * * * MCFM, v5.1 June 1st, 2006 * * * * Authors: John Campbell, Keith Ellis * * (J.Campbell@physics.gla.ac.uk, ellis@fnal.gov) * * * * For details see: * * N.P.B 726:109(2005), hep-ph/0506289 (W+t) * * Phys.Rev.D70:094012, hep-ph/0408158 (Sngl Top)* * (with Francesco Tramontano) * * * * Phys.Rev.D65:113007, hep-ph/0202176 (W,Z+2j) * * Phys.Rev.D62:114012, hep-ph/0006304 (W,Z+bb) * * Phys.Rev.D60:113006, hep-ph/9905386 (diboson) * * * * On the web: http://mcfm.fnal.gov/ * ************************************************** Using input file named input.DAT [Flags to specify the mode in which MCFM is run] evtgen = F creatent = T skipnt = F dswhisto = F [General options to specify the process and execution] nproc= 157 part=virt runstring=cteq6_m sqrts= 1960. ih1 1 ih2 -1 hmass 115. scale 200. facscale 200. dynamicscale F zerowidth T removebr T itmx1 10 ncall1 100000 itmx2 10 ncall2 100000 ij 1089 dryrun F Qflag F Gflag F [Pdf selection] pdlabel cteq6_m NGROUP= 4 NSET= 46 PDFname=cteq6mE.LHgrid PDFmember= -1 [Jet definition and event cuts] m34min 0.1 m34max 1960. m56min 0.1 m56max 1960. inclusive T algorithmcone ptjetmin 10. etajetmin 0. etajetmax 5. Rcut 0.4 makecuts F leptpt 20. leptrap 1.1 misspt 20. leptpt2 0. leptrap2 10. Rjlmin 0.7 Rllmin 0.7 delyjjmin 0. jetsopphem F lbjscheme 0 ptbjetmin 15. etabjetmax 2.8 gammpt 6.5 gammrap 2.5 gammcone 0.7 gammcut 1. [Anomalous couplings of the W and Z] delg1_z 0. delk_z 0. delk_g 0. lambda_z 0. lambda_g 0. tevscale 2. [How to resume/save a run] readin F writeout F ingridfile outgridfile Chooser:process chosen by nproc= 157 *************************** f(p1)+f(p2) --> ************************* * t t~ [for total Xsect] * ********************************************************************* ********************************** * MCFM is calling LHAPDF * * * * PDFname cteq6mE.LHgrid * * PDFmember -1 * ********************************** ************************************* * LHAPDF Version 5.3.0 * ************************************* >>>>>> PDF description: <<<<<< CTEQ6M Reference: J. Pumplin, D.R. Stump, J. Huston, H.L. Lai, P. Nadolsky, W.K. Tung hep-ph/0201195 This set has 41 member PDFs. mem=0 --> central value (cteq6m) mem=1-40 --> 20 sets of +ve/-ve error values >>>>>> <<<<<< **************************************** * Calculating errors using * * 40 sets of error PDFs * **************************************** ************** Electroweak parameters ************** * * * zmass(+) 91.1880000 wmass(+) 80.4190000 * * zwidth(+) 2.4900000 wwidth(+) 2.0600000 * * Gf(+) 0.116639E-04 1/aemmz 132.5069799 * * xw 0.2222465 mtop(+) 200.0000000 * * gwsq 0.4267132 esq 0.0948355 * * top width 2.4416107 * * mb 4.6200000 mc 1.5000000 * * * * Parameters marked (+) are input, others derived * **************************************************** ***************** CKM mixing matrix **************** * * * Vud=.97500 Vus=.22200 Vub=0.0000 * * Vcd=.22200 Vcs=.97500 Vcb=0.0000 * **************************************************** ************* Strong coupling, alpha_s ************ * * * alpha_s (scale) 0.10551725 * * alpha_s (zmass) 0.11798151 * * (using 2-loop running of alpha_s) * **************************************************** **************************************************** * Brn.Rat. removed = 1.0000 * **************************************************** **************************************** * Cross section in femtobarns * **************************************** ****************** Basic mass cuts ***************** * * * 0.10 < m34 < 1960.00 * * 0.10 < m56 < 1960.00 * **************************************************** ******* Cone algorithm additional parameter ******* * * * parton separation parameter, Rsep : 1.300 * * * **************************************************** Heavy Quark mass: 200. <----- Creating batch 0 of n-tuples -----> RZOPEN. record length: 8192 > maximum safe value (8191 words). RZOPEN. You may have problems transferring your file to other systems or writing it to tape. input parameters for vegas: ndim= 5 ncall= 98304. it= 0 itmx= 10 nprn= 0 alph= 1.50 mds= 1 nd= 50 xl( 1)= 0.000 xu( 1)= 1.000 xl( 2)= 0.000 xu( 2)= 1.000 xl( 3)= 0.000 xu( 3)= 1.000 xl( 4)= 0.000 xu( 4)= 1.000 xl( 5)= 0.000 xu( 5)= 1.000 Poles cancelled! ************* Integration by Vegas (iteration 1) ************** * * * integral = 3037.1729 accum. integral = 3037.1729 * * std. dev. = 5.7560195 accum. std. dev = 5.7560195 * * max. wt. = 0.778414 * * * ************** chi**2/iteration = 0.000 **************** ************* Integration by Vegas (iteration 2) ************** * * * integral = 3043.9448 accum. integral = 3042.9989 * * std. dev. = 2.3193112 accum. std. dev = 2.1512411 * * max. wt. = 0.778414 * * * ************** chi**2/iteration = 1.179 **************** ************* Integration by Vegas (iteration 3) ************** * * * integral = 3041.7629 accum. integral = 3042.3585 * * std. dev. = 2.0746319 accum. std. dev = 1.4933353 * * max. wt. = 0.778414 * * * ************** chi**2/iteration = 0.6775 **************** ************* Integration by Vegas (iteration 4) ************** * * * integral = 3044.1541 accum. integral = 3042.9386 * * std. dev. = 2.1615842 accum. std. dev = 1.2286449 * * max. wt. = 0.778414 * * * ************** chi**2/iteration = 0.6076 **************** ************* Integration by Vegas (iteration 5) ************** * * * integral = 3044.3776 accum. integral = 3043.2748 * * std. dev. = 2.2254288 accum. std. dev = 1.0756060 * * max. wt. = 0.778414 * * * ************** chi**2/iteration = 0.5360 **************** ************* Integration by Vegas (iteration 6) ************** * * * integral = 3040.4449 accum. integral = 3042.7557 * * std. dev. = 2.2695243 accum. std. dev = 0.97197166 * * max. wt. = 0.778414 * * * ************** chi**2/iteration = 0.6824 **************** ************* Integration by Vegas (iteration 7) ************** * * * integral = 3042.0719 accum. integral = 3042.6510 * * std. dev. = 2.2854759 accum. std. dev = 0.89444476 * * max. wt. = 0.778414 * * * ************** chi**2/iteration = 0.5815 **************** ************* Integration by Vegas (iteration 8) ************** * * * integral = 3039.7857 accum. integral = 3042.2636 * * std. dev. = 2.2622102 accum. std. dev = 0.83178811 * * max. wt. = 0.778414 * * * ************** chi**2/iteration = 0.6965 **************** ************* Integration by Vegas (iteration 9) ************** * * * integral = 3039.6036 accum. integral = 3041.9543 * * std. dev. = 2.2932538 accum. std. dev = 0.78194115 * * max. wt. = 0.778414 * * * ************** chi**2/iteration = 0.7579 **************** ************* Integration by Vegas (iteration 10) ************** * * * integral = 3044.1017 accum. integral = 3042.1706 * * std. dev. = 2.3367463 accum. std. dev = 0.74152596 * * max. wt. = 0.778414 * * * ************** chi**2/iteration = 0.7581 **************** input parameters for vegas: ndim= 5 ncall= 98304. it= 11 itmx= 10 nprn= 0 alph= 1.50 mds= 1 nd= 50 xl( 1)= 0.000 xu( 1)= 1.000 xl( 2)= 0.000 xu( 2)= 1.000 xl( 3)= 0.000 xu( 3)= 1.000 xl( 4)= 0.000 xu( 4)= 1.000 xl( 5)= 0.000 xu( 5)= 1.000 ************* Integration by Vegas (iteration 1) ************** * * * integral = 3040.9405 accum. integral = 3040.9405 * * std. dev. = 2.3266258 accum. std. dev = 2.3266258 * * max. wt. = 0.119011 * * * ************** chi**2/iteration = 0.000 **************** Writing out intermediate results for iteration 1 ************* Integration by Vegas (iteration 2) ************** * * * integral = 3041.8051 accum. integral = 3041.3655 * * std. dev. = 2.3663617 accum. std. dev = 1.6590432 * * max. wt. = 0.172830 * * * ************** chi**2/iteration = 0.6720E-01 **************** Writing out intermediate results for iteration 2 ************* Integration by Vegas (iteration 3) ************** * * * integral = 3041.6688 accum. integral = 3041.4663 * * std. dev. = 2.3516674 accum. std. dev = 1.3556449 * * max. wt. = 0.172830 * * * ************** chi**2/iteration = 0.3929E-01 **************** Writing out intermediate results for iteration 3 ************* Integration by Vegas (iteration 4) ************** * * * integral = 3044.8433 accum. integral = 3042.2735 * * std. dev. = 2.4188713 accum. std. dev = 1.1825838 * * max. wt. = 0.185903 * * * ************** chi**2/iteration = 0.5190 **************** Writing out intermediate results for iteration 4 ************* Integration by Vegas (iteration 5) ************** * * * integral = 3042.2490 accum. integral = 3042.2689 * * std. dev. = 2.4627466 accum. std. dev = 1.0660480 * * max. wt. = 0.185903 * * * ************** chi**2/iteration = 0.3896 **************** Writing out intermediate results for iteration 5 ************* Integration by Vegas (iteration 6) ************** * * * integral = 3040.7312 accum. integral = 3042.0396 * * std. dev. = 2.5466182 accum. std. dev = 0.98336335 * * max. wt. = 0.185903 * * * ************** chi**2/iteration = 0.3738 **************** Writing out intermediate results for iteration 6 ************* Integration by Vegas (iteration 7) ************** * * * integral = 3043.0499 accum. integral = 3042.1598 * * std. dev. = 2.6757457 accum. std. dev = 0.92300471 * * max. wt. = 0.185903 * * * ************** chi**2/iteration = 0.3325 **************** Writing out intermediate results for iteration 7 ************* Integration by Vegas (iteration 8) ************** * * * integral = 3044.4634 accum. integral = 3042.3825 * * std. dev. = 2.8215793 accum. std. dev = 0.87725973 * * max. wt. = 0.185903 * * * ************** chi**2/iteration = 0.3709 **************** Writing out intermediate results for iteration 8 ************* Integration by Vegas (iteration 9) ************** * * * integral = 3046.9992 accum. integral = 3042.7487 * * std. dev. = 2.9887316 accum. std. dev = 0.84174823 * * max. wt. = 0.185903 * * * ************** chi**2/iteration = 0.5989 **************** Writing out intermediate results for iteration 9 ************* Integration by Vegas (iteration 10) ************** * * * integral = 3051.3988 accum. integral = 3043.2893 * * std. dev. = 3.2602816 accum. std. dev = 0.81502242 * * max. wt. = 0.185903 * * * ************** chi**2/iteration = 1.265 **************** Value of final virt integral is 3043.289 +/- 0.815 fb Total number of shots : 1966081 Total no. failing cuts : 0 Number failing jet cuts : 0 Number failing process cuts : 0 Jet efficiency : 100.00% Cut efficiency : 100.00% Total efficiency : 100.00% Contribution from parton sub-processes: --------------------------------------- GG | 302.17410 9.93% GQ + GQB | -10.72495 -0.35% QG + QBG | -10.71313 -0.35% QQ + QBQB | 0.00000 0.00% QQB | 2762.55325 90.78% --------------------------------------- ************ PDF error analysis ************ * * * PDF error set 0 ---> 3043.815 fb * * PDF error set 1 ---> 3112.296 fb * * PDF error set 2 ---> 2976.227 fb * * PDF error set 3 ---> 3031.583 fb * * PDF error set 4 ---> 3055.851 fb * * PDF error set 5 ---> 3043.103 fb * * PDF error set 6 ---> 3045.791 fb * * PDF error set 7 ---> 2976.949 fb * * PDF error set 8 ---> 3113.209 fb * * PDF error set 9 ---> 3058.167 fb * * PDF error set 10 ---> 3034.258 fb * * PDF error set 11 ---> 3039.845 fb * * PDF error set 12 ---> 3046.188 fb * * PDF error set 13 ---> 3051.612 fb * * PDF error set 14 ---> 3034.752 fb * * PDF error set 15 ---> 3039.881 fb * * PDF error set 16 ---> 3078.017 fb * * PDF error set 17 ---> 3013.240 fb * * PDF error set 18 ---> 3103.656 fb * * PDF error set 19 ---> 3067.671 fb * * PDF error set 20 ---> 3029.021 fb * * PDF error set 21 ---> 3022.116 fb * * PDF error set 22 ---> 3053.440 fb * * PDF error set 23 ---> 2988.192 fb * * PDF error set 24 ---> 2975.705 fb * * PDF error set 25 ---> 3002.957 fb * * PDF error set 26 ---> 3036.960 fb * * PDF error set 27 ---> 3064.196 fb * * PDF error set 28 ---> 3072.570 fb * * PDF error set 29 ---> 2977.629 fb * * PDF error set 30 ---> 3240.614 fb * * PDF error set 31 ---> 2981.321 fb * * PDF error set 32 ---> 2971.901 fb * * PDF error set 33 ---> 2960.780 fb * * PDF error set 34 ---> 3017.454 fb * * PDF error set 35 ---> 3067.177 fb * * PDF error set 36 ---> 3068.015 fb * * PDF error set 37 ---> 3087.202 fb * * PDF error set 38 ---> 3023.204 fb * * PDF error set 39 ---> 3052.442 fb * * PDF error set 40 ---> 3036.462 fb * * * * --------------- SUMMARY ---------------- * * * * Minimum value 2960.780 fb * * Central value 3043.815 fb * * Maximum value 3240.614 fb * * * * Err estimate +/- 179.495 fb * * +ve direction 241.660 fb * * -ve direction 204.924 fb * ******************************************** **************************************************** output files tt_tot_virt_cteq6mE_200_cteq6_m **************************************************** RZOUT: current RZ file cannot support > 64K records or individual directories > 64K RZOUT: previous cycle(s) for this key ( 300) deleted RZOUT: please consult ZEBRA manual for further details <----- Completed a batch of n-tuples ----->