************** 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 sqrts= 1960. ih1 1 ih2 -1 hmass 125. scale 175. facscale 175. dynamicscale F zerowidth T removebr T itmx1 10 ncall1 200000 itmx2 10 ncall2 200000 ij 1089 dryrun F Qflag F Gflag F [Pdf selection] pdlabel cteq6_m NGROUP= 4 NSET= 46 PDFname=cteq61.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 cteq61.LHgrid * * PDFmember -1 * ********************************** ************************************* * LHAPDF Version 5.3.0 * ************************************* >>>>>> PDF description: <<<<<< CTEQ6.1 Reference: J. Pumplin, D.R. Stump, J. Huston, H.L. Lai W.K. Tung, S. Kuhlmann, J. Owens hep-ph/0303013 This set has 41 member PDFs. mem=0 --> central value (cteq61m) 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(+) 172.5000000 * * gwsq 0.4267132 esq 0.0948355 * * top width 1.4803271 * * 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.10744339 * * 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: 172.5 <----- 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= 200000. 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 = 7041.8488 accum. integral = 7041.8488 * * std. dev. = 7.7303739 accum. std. dev = 7.7303739 * * max. wt. = 1.00388 * * * ************** chi**2/iteration = 0.4712E-08 **************** ************* Integration by Vegas (iteration 2) ************** * * * integral = 7035.8668 accum. integral = 7036.7910 * * std. dev. = 3.3045297 accum. std. dev = 3.0385487 * * max. wt. = 1.00388 * * * ************** chi**2/iteration = 0.5013 **************** ************* Integration by Vegas (iteration 3) ************** * * * integral = 7037.4092 accum. integral = 7037.0817 * * std. dev. = 3.2258928 accum. std. dev = 2.2118448 * * max. wt. = 1.00388 * * * ************** chi**2/iteration = 0.2616 **************** ************* Integration by Vegas (iteration 4) ************** * * * integral = 7041.7577 accum. integral = 7038.4690 * * std. dev. = 3.4054358 accum. std. dev = 1.8549273 * * max. wt. = 1.00388 * * * ************** chi**2/iteration = 0.6152 **************** ************* Integration by Vegas (iteration 5) ************** * * * integral = 7044.2204 accum. integral = 7039.7087 * * std. dev. = 3.5387243 accum. std. dev = 1.6429031 * * max. wt. = 1.00388 * * * ************** chi**2/iteration = 0.9785 **************** ************* Integration by Vegas (iteration 6) ************** * * * integral = 7042.7847 accum. integral = 7040.2472 * * std. dev. = 3.5662891 accum. std. dev = 1.4921787 * * max. wt. = 1.00388 * * * ************** chi**2/iteration = 0.9057 **************** ************* Integration by Vegas (iteration 7) ************** * * * integral = 7043.9263 accum. integral = 7040.7760 * * std. dev. = 3.6421301 accum. std. dev = 1.3807868 * * max. wt. = 1.00388 * * * ************** chi**2/iteration = 0.9004 **************** ************* Integration by Vegas (iteration 8) ************** * * * integral = 7038.4470 accum. integral = 7040.4822 * * std. dev. = 3.6344121 accum. std. dev = 1.2907710 * * max. wt. = 1.00388 * * * ************** chi**2/iteration = 0.8232 **************** ************* Integration by Vegas (iteration 9) ************** * * * integral = 7038.4985 accum. integral = 7040.2624 * * std. dev. = 3.6559822 accum. std. dev = 1.2171402 * * max. wt. = 1.00388 * * * ************** chi**2/iteration = 0.7531 **************** ************* Integration by Vegas (iteration 10) ************** * * * integral = 7046.8980 accum. integral = 7040.9117 * * std. dev. = 3.6954716 accum. std. dev = 1.1560512 * * max. wt. = 1.00388 * * * ************** chi**2/iteration = 0.9923 **************** input parameters for vegas: ndim= 5 ncall= 200000. 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 = 7039.7977 accum. integral = 7039.7977 * * std. dev. = 3.6621941 accum. std. dev = 3.6621941 * * max. wt. = 0.451372 * * * ************** chi**2/iteration = 0.000 **************** Writing out intermediate results for iteration 1 ************* Integration by Vegas (iteration 2) ************** * * * integral = 7037.8868 accum. integral = 7038.8536 * * std. dev. = 3.7059062 accum. std. dev = 2.6048793 * * max. wt. = 0.451372 * * * ************** chi**2/iteration = 0.1332 **************** Writing out intermediate results for iteration 2 ************* Integration by Vegas (iteration 3) ************** * * * integral = 7045.5643 accum. integral = 7041.0377 * * std. dev. = 3.7500161 accum. std. dev = 2.1393821 * * max. wt. = 0.451372 * * * ************** chi**2/iteration = 1.142 **************** Writing out intermediate results for iteration 3 ************* Integration by Vegas (iteration 4) ************** * * * integral = 7036.8994 accum. integral = 7040.0470 * * std. dev. = 3.8131646 accum. std. dev = 1.8657867 * * max. wt. = 0.451372 * * * ************** chi**2/iteration = 1.060 **************** Writing out intermediate results for iteration 4 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 -----> <----- Creating batch 1 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. ************* Integration by Vegas (iteration 5) ************** * * * integral = 7043.7927 accum. integral = 7040.7550 * * std. dev. = 3.8646180 accum. std. dev = 1.6802187 * * max. wt. = 0.451372 * * * ************** chi**2/iteration = 0.9856 **************** Writing out intermediate results for iteration 5 ************* Integration by Vegas (iteration 6) ************** * * * integral = 7052.9200 accum. integral = 7042.5794 * * std. dev. = 4.0002156 accum. std. dev = 1.5491139 * * max. wt. = 0.451372 * * * ************** chi**2/iteration = 2.358 **************** Writing out intermediate results for iteration 6 ************* Integration by Vegas (iteration 7) ************** * * * integral = 7044.8176 accum. integral = 7042.8585 * * std. dev. = 4.1038495 accum. std. dev = 1.4492963 * * max. wt. = 0.451372 * * * ************** chi**2/iteration = 2.009 **************** Writing out intermediate results for iteration 7 ************* Integration by Vegas (iteration 8) ************** * * * integral = 7036.7306 accum. integral = 7042.2420 * * std. dev. = 4.3333102 accum. std. dev = 1.3744601 * * max. wt. = 0.451372 * * * ************** chi**2/iteration = 1.979 **************** Writing out intermediate results for iteration 8 ************* Integration by Vegas (iteration 9) ************** * * * integral = 7039.4185 accum. integral = 7042.0119 * * std. dev. = 4.6144560 accum. std. dev = 1.3172674 * * max. wt. = 0.451372 * * * ************** chi**2/iteration = 1.775 **************** Writing out intermediate results for iteration 9 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 -----> <----- Creating batch 2 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. ************* Integration by Vegas (iteration 10) ************** * * * integral = 7043.8478 accum. integral = 7042.1323 * * std. dev. = 4.9732190 accum. std. dev = 1.2733569 * * max. wt. = 0.451372 * * * ************** chi**2/iteration = 1.592 **************** Value of final virt integral is 7042.132 +/- 1.273 fb Total number of shots : 4000001 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 | 1145.97866 16.27% GQ + GQB | -28.56225 -0.41% QG + QBG | -28.57611 -0.41% QQ + QBQB | 0.00000 0.00% QQB | 5953.29197 84.54% --------------------------------------- ************ PDF error analysis ************ * * * PDF error set 0 ---> 7042.168 fb * * PDF error set 1 ---> 7183.881 fb * * PDF error set 2 ---> 6905.933 fb * * PDF error set 3 ---> 7038.514 fb * * PDF error set 4 ---> 7046.938 fb * * PDF error set 5 ---> 7054.498 fb * * PDF error set 6 ---> 7033.232 fb * * PDF error set 7 ---> 6902.296 fb * * PDF error set 8 ---> 7177.233 fb * * PDF error set 9 ---> 7146.170 fb * * PDF error set 10 ---> 6948.147 fb * * PDF error set 11 ---> 7027.586 fb * * PDF error set 12 ---> 7055.124 fb * * PDF error set 13 ---> 7035.696 fb * * PDF error set 14 ---> 7044.491 fb * * PDF error set 15 ---> 7029.735 fb * * PDF error set 16 ---> 7113.224 fb * * PDF error set 17 ---> 6986.246 fb * * PDF error set 18 ---> 7147.620 fb * * PDF error set 19 ---> 7050.984 fb * * PDF error set 20 ---> 7057.532 fb * * PDF error set 21 ---> 7045.178 fb * * PDF error set 22 ---> 7119.875 fb * * PDF error set 23 ---> 6957.844 fb * * PDF error set 24 ---> 6913.607 fb * * PDF error set 25 ---> 6921.734 fb * * PDF error set 26 ---> 6979.722 fb * * PDF error set 27 ---> 7074.915 fb * * PDF error set 28 ---> 7145.634 fb * * PDF error set 29 ---> 6690.655 fb * * PDF error set 30 ---> 7579.616 fb * * PDF error set 31 ---> 7302.040 fb * * PDF error set 32 ---> 6900.780 fb * * PDF error set 33 ---> 6899.273 fb * * PDF error set 34 ---> 6910.899 fb * * PDF error set 35 ---> 7064.650 fb * * PDF error set 36 ---> 7114.125 fb * * PDF error set 37 ---> 7144.629 fb * * PDF error set 38 ---> 7062.821 fb * * PDF error set 39 ---> 7107.967 fb * * PDF error set 40 ---> 7032.896 fb * * * * --------------- SUMMARY ---------------- * * * * Minimum value 6690.655 fb * * Central value 7042.168 fb * * Maximum value 7579.616 fb * * * * Err estimate +/- 549.669 fb * * +ve direction 679.058 fb * * -ve direction 523.377 fb * ******************************************** **************************************************** output files tt_tot_virt_cteq61._175_cteq6 **************************************************** <----- Completed a batch of n-tuples ----->