ios main ios main ios main ios main ios main ios main ios main ios main ios main ios main ios main ios main ios main ios main RANLUX LUXURY LEVEL SET BY RLUXGO : 4 P= 389 RANLUX INITIALIZED BY RLUXGO FROM SEEDS 44788029 0 0 NEW BLOCK DATA for RAPGAP initialised HERACLES 4.6.6 BLOCK DATA initialized 1 ****************************************************************************** ****************************************************************************** ** ** ** ** ** *......* Welcome to the Lund Monte Carlo! ** ** *:::!!:::::::::::* ** ** *::::::!!::::::::::::::* PPP Y Y TTTTT H H III A ** ** *::::::::!!::::::::::::::::* P P Y Y T H H I A A ** ** *:::::::::!!:::::::::::::::::* PPP Y T HHHHH I AAAAA ** ** *:::::::::!!:::::::::::::::::* P Y T H H I A A ** ** *::::::::!!::::::::::::::::*! P Y T H H III A A ** ** *::::::!!::::::::::::::* !! ** ** !! *:::!!:::::::::::* !! This is PYTHIA version 6.413 ** ** !! !* -><- * !! Last date of change: 12 Sep 2007 ** ** !! !! !! ** ** !! !! !! Now is 0 Jan 2000 at 0:00:00 ** ** !! !! ** ** !! lh !! Disclaimer: this program comes ** ** !! !! without any guarantees. Beware ** ** !! hh !! of errors and use common sense ** ** !! ll !! when interpreting results. ** ** !! !! ** ** !! Copyright T. Sjostrand (2007) ** ** ** ** An archive of program versions and documentation is found on the web: ** ** http://www.thep.lu.se/~torbjorn/Pythia.html ** ** ** ** When you cite this program, the official reference is to the 6.4 manual: ** ** T. Sjostrand, S. Mrenna and P. Skands, JHEP05 (2006) 026 ** ** (LU TP 06-13, FERMILAB-PUB-06-052-CD-T) [hep-ph/0603175]. ** ** ** ** Also remember that the program, to a large extent, represents original ** ** physics research. Other publications of special relevance to your ** ** studies may therefore deserve separate mention. ** ** ** ** Main author: Torbjorn Sjostrand; CERN/PH, CH-1211 Geneva, Switzerland, ** ** and Department of Theoretical Physics, Lund University, Lund, Sweden; ** ** phone: + 41 - 22 - 767 82 27; e-mail: torbjorn@thep.lu.se ** ** Author: Stephen Mrenna; Computing Division, Simulations Group, ** ** Fermi National Accelerator Laboratory, MS 234, Batavia, IL 60510, USA; ** ** phone: + 1 - 630 - 840 - 2556; e-mail: mrenna@fnal.gov ** ** Author: Peter Skands; Theoretical Physics Department, ** ** Fermi National Accelerator Laboratory, MS 106, Batavia, IL 60510, USA; ** ** phone: + 1 - 630 - 840 - 2270; e-mail: skands@fnal.gov ** ** ** ** ** ****************************************************************************** ****************************************************************************** PMAS(4,1) changed from 1.50000 to 1.50000 GRAINI: initialisation of RAPGAP Nevents = 2500 COMPUTED MACHINE ACCURACY = 0.11E-18 rapcha:PLEP( 0) changed from the default -0.300E+02 TO-0.276E+02 rapcha:PPIN( 0) changed from the default 0.820E+03 TO 0.920E+03 rapcha (WARNING): LEPI( 0) not changed from the default 11 TO 11 rapcha (WARNING): HADI( 0) not changed from the default 2212 TO 2212 rapcha:QMIN( 0) changed from the default 0.500E+01 TO 0.000E+00 rapcha:QMAX( 0) changed from the default 0.100E+10 TO 0.100E+01 rapcha (WARNING): NFLA( 0) not changed from the default 5 TO 5 rapcha (WARNING): INTE( 0) not changed from the default 0 TO 0 rapcha (WARNING): IPRO( 0) not changed from the default 12 TO 12 rapcha (WARNING): IDIR( 0) not changed from the default 0 TO 0 rapcha (WARNING):LQ2S( 0) not changed from the default 0.337E+01 TO 0.337E+01 rapcha: NFQC( 0) changed from the default 3 TO 5 rapcha (WARNING): IBS ( 0) not changed from the default 0 TO 0 rapcha (WARNING): IRPA( 0) not changed from the default 1 TO 1 rapcha (WARNING): IRPB( 0) not changed from the default 1 TO 1 rapcha (WARNING): IRPC( 0) not changed from the default 1 TO 1 rapcha (WARNING): IRPD( 0) not changed from the default 1 TO 1 rapcha (WARNING): IRPE( 0) not changed from the default 1 TO 1 rapcha (WARNING): IRPF( 0) not changed from the default 1 TO 1 rapcha (WARNING): ISEH( 0) not changed from the default 0 TO 0 rapcha (WARNING): IHFL( 0) not changed from the default 3 TO 3 rapcha (WARNING): IDIF( 0) not changed from the default 0 TO 0 rapcha (WARNING): IFUL( 0) not changed from the default 1 TO 1 rapcha (WARNING): IGRI( 0) not changed from the default 1 TO 1 rapcha (WARNING): IVME( 0) not changed from the default 0 TO 0 rapcha (WARNING): NFRA( 0) not changed from the default 1 TO 1 rapcha (WARNING): IJET( 0) not changed from the default 0 TO 0 rapcha: IRGL( 0) changed from the default 0 TO 1 rapcha (WARNING): IFPS( 0) not changed from the default 3 TO 3 rapcha (WARNING): IPST( 0) not changed from the default 0 TO 0 rapcha (WARNING): IORD( 0) not changed from the default 1 TO 1 rapcha (WARNING): IALP( 0) not changed from the default 1 TO 1 rapcha (WARNING): ITIM( 0) not changed from the default 1 TO 1 rapcha (WARNING): ISOG( 0) not changed from the default 1 TO 1 rapcha (WARNING): Igkt( 0) not changed from the default 1 TO 1 rapcha:KT1 ( 0) changed from the default 0.770E+00 TO 0.700E+00 rapcha (WARNING):KT2 ( 0) not changed from the default 0.440E+00 TO 0.440E+00 rapcha (WARNING): IRSP( 0) not changed from the default 4 TO 4 rapcha: IAKT( 0) changed from the default 0 TO 1 rapcha (WARNING): IREM( 0) not changed from the default 1 TO 1 rapcha (WARNING): IRAS( 0) not changed from the default 1 TO 1 rapcha: IRAM( 0) changed from the default 0 TO 1 rapcha: IQ2S( 0) changed from the default 5 TO 3 rapcha (WARNING):SCAL( 0) not changed from the default 0.100E+01 TO 0.100E+01 rapcha: INPR( 0) changed from the default 9 TO 10150 rapcha: INGA( 0) changed from the default 1 TO 3005002 rapcha (WARNING): ISET( 0) not changed from the default 3 TO 3 rapcha (WARNING): IP2 ( 0) not changed from the default 0 TO 0 rapcha:OMEG( 0) changed from the default 0.100E-01 TO-0.100E-01 rapcha (WARNING):SCQ2( 0) not changed from the default 0.100E+01 TO 0.100E+01 rapcha: INGG( 0) changed from the default -14 TO -31 rapcha: IPOM( 0) changed from the default -10 TO -30 rapcha:EESP( 0) changed from the default 0.850E-01 TO 0.139E+00 rapcha:RRN2( 0) changed from the default 0.470E+01 TO 0.400E+01 rapcha:RALP( 0) changed from the default 0.250E+00 TO 0.300E+00 rapcha:RRXF( 0) changed from the default 0.900E+00 TO 0.800E+00 rapcha (WARNING):RT2M( 0) not changed from the default 0.500E+01 TO 0.500E+01 rapcha: Iqqg( 0) changed from the default 0 TO 1 rapcha (WARNING):C1 ( 0) not changed from the default 0.100E+01 TO 0.100E+01 rapcha (WARNING):Cg ( 0) not changed from the default 0.100E+01 TO 0.100E+01 rapcha: INPI( 0) changed from the default 1 TO 2005001 rapcha: NCAL( 0) changed from the default 60000 TO 100000 rapcha:ACC1( 0) changed from the default 0.750E+00 TO 0.200E+00 rapcha:ACC2( 0) changed from the default 0.500E+00 TO 0.100E+00 COMPUTED MACHINE ACCURACY = 0.11E-18 hercha:PARAMETER XMIN changed from default 0.100000E-02 TO 0.100000E-04 hercha:Parameter INT2( 1) changed from default 0 TO 1 hercha:Parameter INT3( 1) changed from default 0 TO 5 hercha:Parameter INT3( 2) changed from default 0 TO 5 hercha:Parameter INT3( 3) changed from default 0 TO 5 hercha:Parameter NPVE changed from default 10000 TO 50000 hercha:Parameter ISM2( 1) changed from default 0 TO 1 hercha:Parameter ISM3( 1) changed from default 0 TO 1 hercha:Parameter ISM3( 2) changed from default 0 TO 1 hercha:Parameter ISM3( 3) changed from default 0 TO 1 COMPUTED MACHINE ACCURACY = 0.11E-18 pytcha (waring):MSTJ( 22) not changed from the default 1 TO 1 pytcha (waring):MSTU( 111) not changed from the default 1 TO 1 pytcha (waring):MSTU( 112) not changed from the default 5 TO 5 pytcha (waring):PARU( 112) not changed from the default 0.250E+00 TO 0.250E+00 pytcha: PMAS( 21) changed from the default 0.000E+00 TO 0.700E+00 pytcha: PARJ( 82) changed from the default 0.100E+01 TO 0.120E+01 pytcha: PMAS( 134) changed from the default 0.528E+01 TO 0.528E+01 pytcha: PMAS( 137) changed from the default 0.528E+01 TO 0.528E+01 pytcha: PMAS( 140) changed from the default 0.537E+01 TO 0.537E+01 pytcha: PMAS( 210) changed from the default 0.564E+01 TO 0.562E+01 pytcha: PMAS( 134) changed from the default 0.468E+00 TO 0.460E+00 pytcha: PMAS( 137) changed from the default 0.462E+00 TO 0.501E+00 pytcha: PMAS( 140) changed from the default 0.483E+00 TO 0.438E+00 pytcha: PMAS( 210) changed from the default 0.342E+00 TO 0.368E+00 pytcha: PARJ( 76) changed from the default 0.700E+00 TO 0.730E+00 pytcha: PARJ( 77) changed from the default 0.100E+02 TO 0.180E+02 pytcha: MSTJ( 11) changed from the default 4 TO 3 pytcha (waring):MSTJ( 12) not changed from the default 2 TO 2 pytcha (waring):MSTJ( 24) not changed from the default 2 TO 2 pytcha: MSTJ( 46) changed from the default 3 TO 0 pytcha: MSTJ( 51) changed from the default 0 TO 2 pytcha: MSTJ( 52) changed from the default 3 TO 9 pytcha (waring):MSTJ( 53) not changed from the default 0 TO 0 pytcha (waring):MSTJ( 54) not changed from the default 2 TO 2 pytcha: PARJ( 92) changed from the default 0.100E+01 TO 0.111E+01 pytcha: PARJ( 93) changed from the default 0.200E+00 TO 0.341E+00 pytcha: PARJ( 1) changed from the default 0.100E+00 TO 0.108E+00 pytcha: PARJ( 2) changed from the default 0.300E+00 TO 0.286E+00 pytcha: PARJ( 3) changed from the default 0.400E+00 TO 0.690E+00 pytcha (waring):PARJ( 4) not changed from the default 0.500E-01 TO 0.500E-01 pytcha: PARJ( 11) changed from the default 0.500E+00 TO 0.553E+00 pytcha: PARJ( 12) changed from the default 0.600E+00 TO 0.470E+00 pytcha: PARJ( 13) changed from the default 0.750E+00 TO 0.650E+00 pytcha: PARJ( 14) changed from the default 0.000E+00 TO 0.120E+00 pytcha: PARJ( 15) changed from the default 0.000E+00 TO 0.400E-01 pytcha: PARJ( 16) changed from the default 0.000E+00 TO 0.120E+00 pytcha: PARJ( 17) changed from the default 0.000E+00 TO 0.200E+00 pytcha: PARJ( 19) changed from the default 0.100E+01 TO 0.550E+00 pytcha: PARJ( 21) changed from the default 0.360E+00 TO 0.382E+00 pytcha (waring):PARJ( 25) not changed from the default 0.100E+01 TO 0.100E+01 pytcha: PARJ( 26) changed from the default 0.400E+00 TO 0.276E+00 pytcha: PARJ( 41) changed from the default 0.300E+00 TO 0.437E+00 pytcha: PARJ( 42) changed from the default 0.580E+00 TO 0.850E+00 pytcha: PARJ( 54) changed from the default -0.500E-01 TO-0.400E-01 pytcha: PARJ( 55) changed from the default -0.500E-02 TO-0.200E-02 pytcha: PARJ( 81) changed from the default 0.290E+00 TO 0.299E+00 pytcha: PARJ( 82) changed from the default 0.120E+01 TO 0.154E+01 pytcha (waring):PARP( 72) not changed from the default 0.250E+00 TO 0.250E+00 PARU(11) changed from 0.00100 to 0.01000 Event listing (summary) I particle/jet KS KF orig p_x p_y p_z E m 1 !e-! 21 11 0 0.000 0.000 -27.560 27.560 0.001 2 !p+! 21 2212 0 0.000 0.000 920.000 920.000 0.938 sum: 0.00 0.000 0.000 0.000 0.000 0.000 ***************************************************** * * * You are using the RAPGAP MC generator * * version 3.02/00 * * neutral current interaction selected * * gamma + q_pomeron --> q (DIS diffractive) * * mixing of QPM and O(alpha_s) processes selected * * QCD weights calculated in grid * * no cut on max angle of scattered electron * * no cut on min angle of scattered electron * * Q^2 _min according to kinematics * * Q^2 _max = 1.000 * * y_min according to kinematics * * y_max according to kinematics * * Nr of flavours in target sea = 5 * * Nr of flavours for QCDC = 5 * * no p_t ^2 cut imposed for IPRO = 12 * * p_t^2 cut:16.00 for IPRO=15 16.00 for IPRO=13 * *** diffractive hard scattering *** * cuts for diffractive events: * * t_max = 5.000 * * cut on x_f: 0.800 * ***************************************************** ##################################################### # parton shower selection: # # initial state parton parameters: # # ordering in Q2 of emitted partons # # alphas first order with scale Q2 # # time like partons may shower # # soft gluons are resummed # # final state parton shower # ##################################################### ################################################## # proton remnant parameters # # treated inside RAPGAP # # energy sharing IREM = 4 # ################################################## cm energy of ep system 318.468 GEV running alpha_em selected scale for alpha_s: 4*m_q**2 + p_t **2 scale is multiplied by: 1. Lambda_qcd = 0.250 at NF= 5 flavours RASTFU: NG = -31 NPOM = -30 ############################################# # H1QCD 2006 fits Selected # ############################################# # FIT is -31 # fit B # Pomeron Part Only ############################################# [QCD_2006] Initialising H1 2006 DPDF Fit B [QCD_2006] Warning: using DPDFs outside fitted range of 0.0043 < z < 0.8 ; 8.5 < Q2 < 1600 extrapolation used for z,q2: 0.000177622045 1.88031754E-05 charm density defined from F2c*9/8 used only within RAPGAP [QCD_2006] Warning: using DPDFs outside fitted range of 0.0043 < z < 0.8 ; 8.5 < Q2 < 1600 extrapolation used for z,q2: 6.89423487E-08 1.36598297E-08 charm density defined from F2c*9/8 used only within RAPGAP [QCD_2006] Warning: using DPDFs outside fitted range of 0.0043 < z < 0.8 ; 8.5 < Q2 < 1600 extrapolation used for z,q2: 6.08087103E-05 0.00413867272 charm density defined from F2c*9/8 used only within RAPGAP Date: 2011/10/ 6 09:42 ********************************************************** * * * BBBBBBB AAAA SSSSSS EEEEEE SSSSSS * * BB BB AA AA SS SS EE SS SS * * BB BB AA AA SS EE SS * * BBBBBBB AAAAAAAA SSSSSS EEEEEE SSSSSS * * BB BB AA AA SS EE SS * * BB BB AA AA SS SS EE SS SS * * BBBB BB AA AA SSSSSS EEEEEE SSSSSS * * * * BASES Version 5.1 * * coded by S.Kawabata KEK, March 1994 * ********************************************************** << Parameters for BASES >> (1) Dimensions of integration etc. # of dimensions : Ndim = 5 ( 50 at max.) # of Wilds : Nwild = 5 ( 15 at max.) # of sample points : Ncall = 98304(real) 100000(given) # of subregions : Ng = 48 / variable # of regions : Nregion = 8 / variable # of Hypercubes : Ncube = 32768 (2) About the integration variables ------+---------------+---------------+-------+------- i XL(i) XU(i) IG(i) Wild ------+---------------+---------------+-------+------- 1 1.000000E-12 1.000000E+00 1 yes 2 1.000000E-12 1.000000E+00 1 yes 3 1.000000E-12 1.000000E+00 1 yes 4 1.000000E-12 1.000000E+00 1 yes 5 1.000000E-12 1.000000E+00 1 yes ------+---------------+---------------+-------+------- (3) Parameters for the grid optimization step Max.# of iterations: ITMX1 = 30 Expected accuracy : Acc1 = 0.2000 % (4) Parameters for the integration step Max.# of iterations: ITMX2 = 100 Expected accuracy : Acc2 = 0.1000 % Date: 2011/10/ 6 09:42 Convergency Behavior for the Grid Optimization Step ------------------------------------------------------------------------------ <- Result of each iteration -> <- Cumulative Result -> < CPU time > IT Eff R_Neg Estimate Acc % Estimate(+- Error )order Acc % ( H: M: Sec ) ------------------------------------------------------------------------------ [QCD_2006] Warning: using DPDFs outside fitted range of 0.0043 < z < 0.8 ; 8.5 < Q2 < 1600 extrapolation used for z,q2: 0.209498674 0.0211123936 charm density defined from F2c*9/8 used only within RAPGAP [QCD_2006] Warning: using DPDFs outside fitted range of 0.0043 < z < 0.8 ; 8.5 < Q2 < 1600 extrapolation used for z,q2: 0.268383622 0.0495216474 charm density defined from F2c*9/8 used only within RAPGAP 1 19 0.00 7.864E+11 16.181 7.864317(+-1.272517)E 11 16.181 0: 0: 2.58 2 56 0.00 9.117E+11 0.741 9.113461(+-0.067436)E 11 0.740 0: 0: 6.01 3 59 0.00 9.580E+11 5.388 9.121296(+-0.066868)E 11 0.733 0: 0: 9.56 4 51 0.00 9.116E+11 1.122 9.119712(+-0.055971)E 11 0.614 0: 0:12.69 5 60 0.00 9.109E+11 0.655 9.114527(+-0.040817)E 11 0.448 0: 0:16.29 6 58 0.00 9.081E+11 0.491 9.099266(+-0.030120)E 11 0.331 0: 0:19.86 7 61 0.00 9.067E+11 0.370 9.084805(+-0.022421)E 11 0.247 0: 0:23.80 8 62 0.00 9.025E+11 0.355 9.065172(+-0.018368)E 11 0.203 0: 0:27.90 9 61 0.00 9.036E+11 0.395 9.059076(+-0.016330)E 11 0.180 0: 0:31.90 ------------------------------------------------------------------------------ Date: 2011/10/ 6 09:42 Convergency Behavior for the Integration Step ------------------------------------------------------------------------------ <- Result of each iteration -> <- Cumulative Result -> < CPU time > IT Eff R_Neg Estimate Acc % Estimate(+- Error )order Acc % ( H: M: Sec ) ------------------------------------------------------------------------------ 1 60 0.00 9.166E+11 1.178 9.165857(+-0.107982)E 11 1.178 0: 0:35.50 2 60 0.00 9.010E+11 0.396 9.024945(+-0.033895)E 11 0.376 0: 0:39.08 3 60 0.00 9.044E+11 0.394 9.033819(+-0.024560)E 11 0.272 0: 0:42.65 4 60 0.00 9.055E+11 0.366 9.041309(+-0.019726)E 11 0.218 0: 0:46.25 5 60 0.00 9.117E+11 0.383 9.059664(+-0.017171)E 11 0.190 0: 0:49.81 6 60 0.00 9.080E+11 0.572 9.061645(+-0.016303)E 11 0.180 0: 0:53.40 7 60 0.00 9.117E+11 0.379 9.071724(+-0.014746)E 11 0.163 0: 0:56.99 8 60 0.00 9.185E+11 1.229 9.073619(+-0.014622)E 11 0.161 0: 1: 0.56 9 60 0.00 9.045E+11 0.426 9.070021(+-0.013671)E 11 0.151 0: 1: 4.16 10 60 0.00 9.113E+11 0.412 9.075021(+-0.012847)E 11 0.142 0: 1: 7.79 11 60 0.00 9.014E+11 0.367 9.066958(+-0.011974)E 11 0.132 0: 1:11.36 12 60 0.00 9.088E+11 0.626 9.067858(+-0.011717)E 11 0.129 0: 1:14.91 13 60 0.00 9.066E+11 0.374 9.067710(+-0.011074)E 11 0.122 0: 1:18.50 14 60 0.00 9.022E+11 0.422 9.064131(+-0.010634)E 11 0.117 0: 1:22.08 15 60 0.00 9.101E+11 0.375 9.067373(+-0.010152)E 11 0.112 0: 1:25.68 16 60 0.00 9.049E+11 0.412 9.066126(+-0.009795)E 11 0.108 0: 1:29.26 17 60 0.00 9.045E+11 0.385 9.064551(+-0.009429)E 11 0.104 0: 1:32.83 18 60 0.00 9.043E+11 0.371 9.062971(+-0.009077)E 11 0.100 0: 1:36.41 19 60 0.00 9.141E+11 0.854 9.064018(+-0.009016)E 11 0.099 0: 1:40.00 ------------------------------------------------------------------------------ ****** END OF BASES ********* << Computing Time Information >> (1) For BASES H: M: Sec Overhead : 0: 0: 0.00 Grid Optim. Step : 0: 0:31.90 Integration Step : 0: 1: 8.08 Go time for all : 0: 1:40.00 (2) Expected event generation time Expected time for 1000 events : 0.05 Sec calculating xpq_dif in grid IFULL/IDIR/IDISDIF check existence of file with pdf"s file 31 0 file found with consistent parameters xpq for DIF done 0 0 1 -31 -30 10150 Now calcualting order alpha_s matrix elements Please be patient, for diffraction it can take long | y | q2 | x_pom | qpm | qqb | qcdc | ccb | bbb | .744E+00 |.538E+00 |.141E-03 |.133E+02 |.000E+00 |.000E+00 |.141E+02 |.000E+00 | QCDGRID: RSUM > RQPM 14.0971327 13.2508831 including 5% safety margin increase pt cutoff PT2CUT program stopped now