Welcome to Sherpa, . Initialization of framework underway. The local time is Thu Nov 12 11:31:35 2015. Run_Parameter::Init(): Seed mode '*' Run_Parameter::Init(): Setting memory limit to 9.55738 GB. Random::SetSeed(): Seed set to 1234 ----------------------------------------------------------------------------- ----------- Event generation run with SHERPA started ....... ----------- ----------------------------------------------------------------------------- ................................................ | + ................................................ || | + + ................................... .... | | / + ................. ................ _,_ | .... || +| + + ............................... __.' ,\| ... || / +| + .............................. ( \ \ ... | | | + + \ + ............................. ( \ -/ .... || + | + ........ ................... (~~~~~~~~~## + + + ............................. ~~(! '~~~~~~~ \ + + + + ............................... `~~~QQQQQDb // | + + + + ........................ .......... IDDDDP|| \ + + + + + + .................................... IDDDI|| \ + .................................... IHD HD|| \ + + + + + + + + ................................... IHD ##| :-) + +\ + ......... ............... ......... IHI ## / / + + + + +\ + ................................... IHI/ / / + + + + + ................................... ## | | / / + + + + / + ....................... /TT\ ..... ##/ /// / + + + + + + +/ + ......................./TTT/T\ ... /TT\/\\\ / + + + + + + +/ \ + ....................../TTT/TTTT\...|TT/T\\\/ + ++ + / ----------------------------------------------------------------------------- SHERPA version 2.2.0 (Cho Oyu) Authors: Stefan Hoeche, Frank Krauss, Silvan Kuttimalai, Marek Schoenherr, Holger Schulz, Steffen Schumann, Frank Siegert, Korinna Zapp. Former Authors: Timo Fischer, Tanju Gleisberg, Hendrik Hoeth, Ralf Kuhn, Thomas Laubrich, Andreas Schaelicke, Jan Winter This program uses a lot of genuine and original research work by other people. Users are encouraged to refer to the various original publications. Users are kindly asked to refer to the documentation published under JHEP 02(2009)007 Please visit also our homepage http://sherpa.hepforge.org for news, bugreports, updates and new releases. ----------------------------------------------------------------------------- SVN branch trunk/SHERPA, revision 26832. ===> PHASIC++/Channels has revision 26832M, first seen was 26832 <=== ===> PHASIC++/Main has revision 26832M, first seen was 26832 <=== ===> SHERPA/PerturbativePhysics has revision 26832M, first seen was 26832 <=== Initialization_Handler::SetGlobalVariables(): Set cluster mode 16. Beam_Spectra_Handler : type = Monochromatic*Monochromatic for P+ ((4000,0,0,4000)) and P+ ((4000,0,0,-4000)) PDF set 'NNPDF30NNLO' loaded for beam 1 (P+). PDF set 'NNPDF30NNLO' loaded for beam 2 (P+). Initialized the ISR: (SF)*(SF) One_Running_AlphaS::One_Running_AlphaS() { Setting \alpha_s according to PDF perturbative order 2 \alpha_s(M_Z) = 0.118 } One_Running_AlphaS::One_Running_AlphaS() { Setting \alpha_s according to PDF perturbative order 2 \alpha_s(M_Z) = 0.118 } List of Particle Data IDName kfc MASS[] WIDTH[] STABLE[] MASSIVE[] ACTIVE[] YUKAWA[] d 1 0.01 0 1 0 1 0 u 2 0.005 0 1 0 1 0 s 3 0.2 0 1 0 1 0 c 4 1.42 0 1 0 1 0 b 5 4.8 0 1 0 1 0 t 6 173.21 2 0 1 1 173.21 e- 11 0.000511 0 1 0 1 0 ve 12 0 0 1 0 1 0 mu- 13 0.105 0 1 0 1 0 vmu 14 0 0 1 0 1 0 tau- 15 1.777 2.26735e-12 0 0 1 0 vtau 16 0 0 1 0 1 0 G 21 0 0 1 0 1 0 P 22 0 0 1 0 1 0 Z 23 91.1876 0 0 1 1 91.1876 W+ 24 80.385 0 0 1 1 80.385 h0 25 125 0.00407 0 1 1 125 List of Particle Containers IDName kfc Constituents l 90 {e-,e+,mu-,mu+,tau-,tau+} v 91 {ve,veb,vmu,vmub,vtau,vtaub} j 93 {d,db,u,ub,s,sb,c,cb,b,bb,G} Q 94 {d,db,u,ub,s,sb,c,cb,b,bb} r 99 {d,db,u,ub,s,sb,c,cb,b,bb,G} Initialized the Fragmentation_Handler. Initialized the Soft_Collision_Handler. CS_Shower::CS_Shower(): Set respect Q2 mode 0 CS_Shower::CS_Shower(): Set exclusive cluster mode 1. CS_Shower::CS_Shower(): Set color setter mode 0 CS_Shower::CS_Shower(): Set respect Q2 mode 0 CS_Shower::CS_Shower(): Set exclusive cluster mode 1. CS_Shower::CS_Shower(): Set color setter mode 0 Initialized the Shower_Handler. +-----------------------------------------+ | X X X XXXX XXX XXX XXX | | X X XX XX X X X X X X | | X X X X X XXX X XXX X X XXX XXX | | XXXXX X X X X X X X X X | | X X X X XXXX XXX XXX XXX | +-----------------------------------------+ | please cite: JHEP 0202:044,2002 | +-----------------------------------------+ ME_Generator_Base::SetPSMasses(): Massive PS flavours for Comix: (c,cb,b,bb,e-,e+,mu-,mu+,tau-,tau+) +----------------------------------+ | | | CCC OOO M M I X X | | C O O MM MM I X X | | C O O M M M I X | | C O O M M I X X | | CCC OOO M M I X X | | | +==================================+ | Color dressed Matrix Elements | | http://comix.freacafe.de | | please cite JHEP12(2008)039 | +----------------------------------+ ME_Generator_Base::SetPSMasses(): Massive PS flavours for Amegic: (c,cb,b,bb,e-,e+,mu-,mu+,tau-,tau+) Amegic::Initialize(): Set gauge 1. ME_Generator_Base::SetPSMasses(): Massive PS flavours for Internal: (c,cb,b,bb,e-,e+,mu-,mu+,tau-,tau+) Matrix_Element_Handler::BuildProcesses(): Looking for processes ...................................................................Caution: Using Dummy_QCD_Virtual for (db,u,Z,W+,G) ..Caution: Using Dummy_QCD_Virtual for (db,G,Z,W+,ub) ..Caution: Using Dummy_QCD_Virtual for (u,db,Z,W+,G) ..Caution: Using Dummy_QCD_Virtual for (u,G,Z,W+,d) ..Caution: Using Dummy_QCD_Virtual for (G,db,Z,W+,ub) ..Caution: Using Dummy_QCD_Virtual for (G,u,Z,W+,d) ............................................................................................................. done ( 71 MB, 2s / 1s ). Matrix_Element_Handler::InitializeProcesses(): Performing tests .......................................................................................................... done ( 75 MB, 0s / 0s ). Initialized the Matrix_Element_Handler for the hard processes. Initialized the Beam_Remnant_Handler. Initialized the Soft_Photon_Handler. Read in channels from directory : Results/Amegic/MC_2_3__j__j__Z__W+__j__QCD(BVI) Process_Group::CalculateTotalXSec(): Calculate xs for '2_3__j__j__Z__W+__j__QCD(BVI)' (Amegic) 2_3__j__j__Z__W+__j__QCD(BVI) : 8.9961 pb +- ( 0.103831 pb = 1.15417 % )  exp. eff: 0.193594 % Read in channels from directory : Results/Comix/MC_2_4__j__j__Z__W+__j__j__QCD(RS) Process_Group::CalculateTotalXSec(): Calculate xs for '2_4__j__j__Z__W+__j__j__QCD(RS)' (Comix) 2_4__j__j__Z__W+__j__j__QCD(RS) : 0.0772886 pb +- ( 0.181296 pb = 234.57 % )  exp. eff: 0.00052491 % reduce max for 2_3__j__j__Z__W+__j__QCD(BVI) to 0.474434 ( eps = 0.001 ) reduce max for 2_4__j__j__Z__W+__j__j__QCD(RS) to 0.938379 ( eps = 0.001 ) ---------------------------------------------------------- -- SHERPA generates events with the following structure -- ---------------------------------------------------------- Event generation : Weighted Perturbative : Signal_Processes Perturbative : Hard_Decays Perturbative : Jet_Evolution:CSS Perturbative : Lepton_FS_QED_Corrections:Photons Perturbative : Multiple_Interactions:None Perturbative : Minimum_Bias:Off Hadronization : Beam_Remnants Hadronization : Hadronization:Off Hadronization : Hadron_Decays --------------------------------------------------------- Event_Handler::GenerateEvent(43) { MCatNLO_Process::OneEvent() { BBar_Emission_Generator::GeneratePoint() { in EEG: #1 = 0.3488, #2 = 0.4803, #3 = 0.8785 Born point { (92.06,0,0,92.06) (188.1,0,0,-188.1) (108.4,-11.47,46,-34.56) (98.63,33.83,-45.34,8.123) (73.06,-22.36,-0.6617,-69.55) } selected FI_4[4,5]_0[0] FI_Dipole::GeneratePoint() { vegased : x = 0.3488, z = 0.4803, phi = 0.8785, xmin = 0.02302 transformed : x = 0.2004, z = 0.4767, phi = 5.52 } } BBar_Emission_Generator::GenerateWeight() { Dipole FI_4[4,5]_0[0] { again : x = 0.2004, z = 0.4767, phi = 5.52 recovered : x = 0.3488, z = 0.4803, phi = 0.8785, xmin = 0.02302 } -> w = 3154 ( a = 0.1 ) } Single_Process::Differential(2_3__sb__G__Z__W+__cb__QCD(BVI), RS:0) { Scale_Setter_Base::CalculateScale(2_3__db__G__Z__W+__ub__QCD(BVI)) { Variable_Scale_Setter::Calculate(): Set { \mu_f = 90 \mu_r = 90 \mu_q = 90 } <- 2_3__db__G__Z__W+__ub__QCD(BVI) Coupling_Data::Calculate(0x1dfa960): scl = 90 -> 'Alpha_QCD'{fac=1.002,cpl=0.118} } KP_Terms::Get() { } PDF(fla=sb, xa=0.02302, ta=8100) * PDF(flb=G, xb=0.04701, tb=8100) -> 528.3 type=BVIKP, B=6.96e-11, VI=3.19e-11, KP=2.885e-11, K=1 muR2=8100, muF2=8100, oqcd=2, oew=2, fl1=-3, fl2=21, x1=1, x2=1, x1p=0.8527, x2p=0.5787 wren=(3.537e-11,2.83e-11) wfac=(4.495e-11,-3.045e-11,0,-6.266e-13,0,-1.014e-12,2.98e-11,-1.692e-11,2.005e-11,-2e-11,-0,-4.791e-13,-0,-3.387e-12,1.313e-11,-1.756e-11) Cluster sequence: pdfwgt=528.3, flux=7.22e-06, counter term=0, no cluster steps Calculating DADS terms Single_Process::Differential(2_4__sb__G__Z__W+__G__cb__QCD(S)_RS4_5_0, RS:0) { Scale_Setter_Base::CalculateScale(2_3__db__G__Z__W+__ub__QCD(S)__S4_0) { Variable_Scale_Setter::Calculate(): Set { \mu_f = 90 \mu_r = 90 \mu_q = 90 } <- 2_3__db__G__Z__W+__ub__QCD(S)__S4_0 Coupling_Data::Calculate(0x220a020): scl = 90 -> 'Alpha_QCD'{fac=1.002,cpl=0.118} } 586:m_lastxs=-9.93e-12 PDF(fla=sb, xa=0.1149, ta=8100) * PDF(flb=G, xb=0.04701, tb=8100) -> 17.63 type=none, B=-1.436e-17, VI=0, KP=0, K=0 muR2=8100, muF2=8100, oqcd=2, oew=2, fl1=-3, fl2=21, x1=1, x2=1, x1p=1, x2p=1 Cluster sequence: pdfwgt=17.63, flux=1.447e-06, counter term=0, no cluster steps } DADS_0 = 7.986e-12 DADS: wgt=4.531e-13, x1=0.1149, x2=0.04701, fl1=-3, fl2=21 } 2_3__sb__G__Z__W+__cb__QCD(BVI) ME = 6.8871932525051e-08 , PS = 0.0387272785114049 -> 2.66722251251634e-09 p_lab[0]=Vec4D(92.0649834855704,0,0,92.0649834855705); p_lab[1]=Vec4D(188.055857090885,0,0,-188.055857090885); p_lab[2]=Vec4D(108.428858427042,-11.4677435124342,45.9976285173266,-34.5593399176929); p_lab[3]=Vec4D(98.6280384747466,33.8290725544059,-45.3359439735672,8.12328643532916); p_lab[4]=Vec4D(73.0639436746661,-22.3613290419717,-0.661684543759371,-69.5548201229506); ========================================================== S event Cluster_Algorithm::Cluster() { Combine_Table::CalcJet(mode = 512, nl = 5, complete = 1) { Combine_Table (0x3febf80) 1 [OQCD=1] (up=#) id content flavcut qcd qed mom 0 (0) sb 0 0/0 0/0 (92.06,0,0,92.06) 1 (1) G 0 0/0 0/0 (188.1,0,0,-188.1) 2 (2) Z 0 0/0 0/0 (108.4,-11.47,46,-34.56) 3 (3) W+ 0 0/0 0/0 (98.63,33.83,-45.34,8.123) 4 (4) cb 0 0/0 0/0 (73.06,-22.36,-0.6617,-69.55) --------------- 0&2%1[ sb] CP{kt=0,op=0,x=0,mu=0,k=0,m=0} 0 6 0&2%3[ sb] CP{kt=0,op=0,x=0,mu=0,k=0,m=0} 0 2,6 0&2%4[ sb] CP{kt=0,op=0,x=0,mu=0,k=0,m=0} 0 2,6 0&3%1[ cb] CP{kt=0,op=0,x=0,mu=0,k=0,m=0} 0 3 0&3%2[ cb] CP{kt=0,op=0,x=0,mu=0,k=0,m=0} 0 1,3 0&3%4[ cb] CP{kt=0,op=0,x=0,mu=0,k=0,m=0} 0 1,3 1&4%0[ cb] CP{kt=0,op=0,x=0,mu=0,k=0,m=0} 0 0,1,2 1&4%2[ cb] CP{kt=0,op=0,x=0,mu=0,k=0,m=0} 0 1 1&4%3[ cb] CP{kt=0,op=0,x=0,mu=0,k=0,m=0} 0 2 2&3%0[ W+] CP{kt=0,op=0,x=0,mu=0,k=0,m=0} 0 0,4 2&3%4[ W+] CP{kt=0,op=0,x=0,mu=0,k=0,m=0} 0 0,4 2&4%0[ cb] CP{kt=0,op=0,x=0,mu=0,k=0,m=0} 0 3,5 2&4%1[ cb] CP{kt=0,op=0,x=0,mu=0,k=0,m=0} 0 3 2&4%3[ cb] CP{kt=0,op=0,x=0,mu=0,k=0,m=0} 0 3,5 3&4%0[ sb] CP{kt=0,op=0,x=0,mu=0,k=0,m=0} 0 6,7 3&4%1[ sb] CP{kt=0,op=0,x=0,mu=0,k=0,m=0} 0 6 3&4%2[ sb] CP{kt=0,op=0,x=0,mu=0,k=0,m=0} 0 6,7 k_{T,min} (0,1,2,3,4) -> 0 CS_Cluster_Definitions::KPerp2(kin = -1, mode = 16) { Set fixed n_perp No splitting function (j) } Calculate m_perp(0[sb],2[Z],1[G],sb) -> CP{kt=143.4,op=-1,x=0.2186,mu=143.4,k=1,m=0} CS_Cluster_Definitions::KPerp2(kin = -1, mode = 16) { No splitting function (j) } Calculate m_perp(0[sb],2[Z],3[W+],sb) -> CP{kt=128.3,op=-1,x=0.3744,mu=128.3,k=1,m=0} CS_Cluster_Definitions::KPerp2(kin = -1, mode = 16) { No splitting function (j) } Calculate m_perp(0[sb],2[Z],4[cb],sb) -> CP{kt=72.8,op=-1,x=0.7987,mu=72.8,k=1,m=0} CS_Cluster_Definitions::KPerp2(kin = -1, mode = 16) { Set fixed n_perp No splitting function (j) } Calculate m_perp(0[sb],3[W+],1[G],cb) -> CP{kt=116.9,op=-1,x=0.1796,mu=116.9,k=1,m=0} CS_Cluster_Definitions::KPerp2(kin = -1, mode = 16) { No splitting function (j) } Calculate m_perp(0[sb],3[W+],2[Z],cb) -> CP{kt=102.1,op=-1,x=0.3744,mu=102.1,k=1,m=0} CS_Cluster_Definitions::KPerp2(kin = -1, mode = 16) { No splitting function (j) } Calculate m_perp(0[sb],3[W+],4[cb],cb) -> CP{kt=81.23,op=-1,x=0.6041,mu=81.23,k=1,m=0} CS_Cluster_Definitions::KPerp2(kin = -1, mode = 16) { Set fixed n_perp Kernel weight (A=1/NLO=16) [m=3,c=1] ( x = 0.02829 ) LF_VFF_II|CF_QCD { (1) G (-188.1,-0,-0,188.1) -5.395e-06 (-1,-1) [0|0] (4) cb (73.06,-22.36,-0.6617,-69.55) -9.537e-07 (-1,-1) [0|0] (0) s (-92.06,-0,-0,-92.06) -4.046e-06 (-1,-1) [0|0] } -> w = 0.0006027 (1659) } Calculate m_perp(1[G],4[cb],0[sb],cb) -> CP{kt=36.33,op=40.73,x=0.6017,mu=16.37,k=1,m=0} CS_Cluster_Definitions::KPerp2(kin = -1, mode = 16) { Invalid spectator Z } Calculate m_perp(1[G],4[cb],2[Z],cb) -> CP{kt=36.33,op=-1,x=0.6363,mu=36.33,k=1,m=0} CS_Cluster_Definitions::KPerp2(kin = -1, mode = 16) { Invalid spectator W+ } Calculate m_perp(1[G],4[cb],3[W+],cb) -> CP{kt=36.33,op=-1,x=0.5902,mu=36.33,k=1,m=0} CS_Cluster_Definitions::KPerp2(kin = -1, mode = 17) { No splitting function (i) } Calculate m_perp(2[Z],3[W+],0[sb],W+) -> CP{kt=79.9,op=-1,x=0.3744,mu=79.9,k=0,m=1} CS_Cluster_Definitions::KPerp2(kin = -1, mode = 17) { No splitting function (i) } Calculate m_perp(2[Z],3[W+],4[cb],W+) -> CP{kt=79.75,op=-1,x=1,mu=79.75,k=0,m=1} CS_Cluster_Definitions::KPerp2(kin = -1, mode = 16) { No splitting function (i) } Calculate m_perp(2[Z],4[cb],0[sb],cb) -> CP{kt=72.8,op=-1,x=0.7987,mu=72.8,k=1,m=0} CS_Cluster_Definitions::KPerp2(kin = -1, mode = 16) { No splitting function (i) } Calculate m_perp(2[Z],4[cb],1[G],cb) -> CP{kt=100.5,op=-1,x=0.6363,mu=100.5,k=1,m=0} CS_Cluster_Definitions::KPerp2(kin = -1, mode = 16) { No splitting function (i) } Calculate m_perp(2[Z],4[cb],3[W+],cb) -> CP{kt=80.54,op=-1,x=1,mu=80.54,k=1,m=0} CS_Cluster_Definitions::KPerp2(kin = -1, mode = 16) { No splitting function (i) } Calculate m_perp(3[W+],4[cb],0[sb],sb) -> CP{kt=81.23,op=-1,x=0.6041,mu=81.23,k=1,m=0} CS_Cluster_Definitions::KPerp2(kin = -1, mode = 16) { No splitting function (i) } Calculate m_perp(3[W+],4[cb],1[G],sb) -> CP{kt=128.3,op=-1,x=0.5902,mu=128.3,k=1,m=0} CS_Cluster_Definitions::KPerp2(kin = -1, mode = 16) { No splitting function (i) } Calculate m_perp(3[W+],4[cb],2[Z],sb) -> CP{kt=110.4,op=-1,x=1,mu=110.4,k=1,m=0} Combine_Table (0x3febf80) 1 [OQCD=1] (up=#) idcontent flav cut qcd qed mom 0 (0) sb 0 0/0 0/0 (92.06,0,0,92.06) 1 (1) G 0 0/0 0/0 (188.1,0,0,-188.1) 2 (2) Z 0 0/0 0/0 (108.4,-11.47,46,-34.56) 3 (3) W+ 0 0/0 0/0 (98.63,33.83,-45.34,8.123) 4 (4) cb 0 0/0 0/0 (73.06,-22.36,-0.6617,-69.55) --------------- 0&2%1[sb ] CP{kt=143.4,op=-1,x=0.2186,mu=143.4,k=1,m=0} 0 6 0&2%3[sb ] CP{kt=128.3,op=-1,x=0.3744,mu=128.3,k=1,m=0} 0 2,6 0&2%4[sb ] CP{kt=72.8,op=-1,x=0.7987,mu=72.8,k=1,m=0} 0 2,6 0&3%1[cb ] CP{kt=116.9,op=-1,x=0.1796,mu=116.9,k=1,m=0} 0 3 0&3%2[cb ] CP{kt=102.1,op=-1,x=0.3744,mu=102.1,k=1,m=0} 0 1,3 0&3%4[cb ] CP{kt=81.23,op=-1,x=0.6041,mu=81.23,k=1,m=0} 0 1,3 1&4%0[cb ] CP{kt=36.33,op=40.73,x=0.6017,mu=16.37,k=1,m=0} 0 0,1,2 <- 1&4%2[cb ] CP{kt=36.33,op=-1,x=0.6363,mu=36.33,k=1,m=0} 0 1 1&4%3[cb ] CP{kt=36.33,op=-1,x=0.5902,mu=36.33,k=1,m=0} 0 2 2&3%0[W+ ] CP{kt=79.9,op=-1,x=0.3744,mu=79.9,k=0,m=1} 0 0,4 2&3%4[W+ ] CP{kt=79.75,op=-1,x=1,mu=79.75,k=0,m=1} 0 0,4 2&4%0[cb ] CP{kt=72.8,op=-1,x=0.7987,mu=72.8,k=1,m=0} 0 3,5 2&4%1[cb ] CP{kt=100.5,op=-1,x=0.6363,mu=100.5,k=1,m=0} 0 3 2&4%3[cb ] CP{kt=80.54,op=-1,x=1,mu=80.54,k=1,m=0} 0 3,5 3&4%0[sb ] CP{kt=81.23,op=-1,x=0.6041,mu=81.23,k=1,m=0} 0 6,7 3&4%1[sb ] CP{kt=128.3,op=-1,x=0.5902,mu=128.3,k=1,m=0} 0 6 3&4%2[sb ] CP{kt=110.4,op=-1,x=1,mu=110.4,k=1,m=0} 0 6,7 k_{T,min} (0,1,2,3,4) -> 0 set last k_T = 36.33 for (0,1,2,3,4) from (1,4) set last k_T = 36.33 for (0,1,2,3,4) from (1,4) Combine_Table::CalcJet(mode = 512, nl = 4, complete = 1) { Combine_Table (0x4152d00) 2 [OQCD=0] (up=1) idcontent flav cut qcd qed mom 0 (0) sb 0 0/0 0/0 (92.06,0,0,92.06) 1 (1,4) cb 0 1/3 0/0 (113.2,0,0,-113.2) 2 (2) Z 0 0/0 0/0 (109.3,-23.04,45.66,-31.99) 3 (3) W+ 0 0/0 0/0 (95.89,23.04,-45.66,10.89) --------------- k_{T,min} (0,1,2,3,4) -> 36.33 Default_Core_Scale::Calculate() { (0x3daa1c0): 2 -> 2 { \mu_r = 0, \mu_f = 0, \mu_q = 0, \mu = 0 k_T = 0, z = 0, phi = 0, kin = 0, K = 0 oew = 0, oqcd = 0, nlo = 0, new = (none), ncl = 0, flag = 0 (0) s (-92.06,-0,-0,-92.06) -4.046e-06 (-1,-1) [0|0] (1,4) c (-113.2,-0,-0,113.2) -3.016e-06 (-1,-1) [0|0] (2) Z (109.3,-23.04,45.66,-31.99) 91.19 (-1,-1) [0|0] (3) W+ (95.89,23.04,-45.66,10.89) 80.39 (-1,-1) [0|0] } DY like \mu_f = 204.1 \mu_r = 204.1 \mu_q = 204.1 } } } Cluster_Algorithm::Convert(): { Default_Core_Scale::Calculate() { (0x3e1f510): 2 -> 2 { \mu_r = 90, \mu_f = 90, \mu_q = 90, \mu = 0 k_T = 0, z = 0, phi = 0, kin = 1, K = 0 oew = 2, oqcd = 1, nlo = 0, new = (none), ncl = 0, flag = 0 (0) s (-92.06,-0,-0,-92.06) -4.046e-06 (0,0) [1|0] (1,4) c (-113.2,-0,-0,113.2) -3.016e-06 (0,0) [1|0] <-> (0) (2) Z (109.3,-23.04,45.66,-31.99) 91.19 (0,0) [1|0] (3) W+ (95.89,23.04,-45.66,10.89) 80.39 (0,0) [1|0] } DY like \mu_f = 204.1 \mu_r = 204.1 \mu_q = 204.1 } (0x3e1f510): 2 -> 2 { \mu_r = 90, \mu_f = 90, \mu_q = 204.1, \mu = 204.1 k_T = 204.1, z = 0, phi = 0, kin = 1, K = 0 oew = 2, oqcd = 1, nlo = 0, new = (none), ncl = 0, flag = 0 (0) s (-92.06,-0,-0,-92.06) -4.046e-06 (0,0) [1|4] (1,4) c (-113.2,-0,-0,113.2) -3.016e-06 (0,0) [1|4] <-> (0) (2) Z (109.3,-23.04,45.66,-31.99) 91.19 (0,0) [1|4] (3) W+ (95.89,23.04,-45.66,10.89) 80.39 (0,0) [1|4] } (0x3daa1c0): 2 -> 3 { \mu_r = 90, \mu_f = 90, \mu_q = 90, \mu = 16.37 k_T = 36.33, z = 0, phi = 0, kin = 0, K = 0 oew = 2, oqcd = 2, nlo = 0, new = (4), ncl = 0, flag = 0 (0) s (-92.06,-0,-0,-92.06) -4.046e-06 (0,0) [0|4] (1) G (-188.1,-0,-0,188.1) -5.395e-06 (0,0) [0|4] (2) Z (108.4,-11.47,46,-34.56) 91.19 (0,0) [0|4] (3) W+ (98.63,33.83,-45.34,8.123) 80.39 (0,0) [0|4] (4) cb (73.06,-22.36,-0.6617,-69.55) -9.537e-07 (0,0) [0|4] } } } CS_MCatNLO::GeneratePoint() { CS_MCatNLO::PrepareMCatNLO(): { (0x3daa1c0): 2 -> 3 { \mu_r = 90, \mu_f = 90, \mu_q = 90, \mu = 16.37 k_T = 36.33, z = 0, phi = 0, kin = 0, K = 0 oew = 2, oqcd = 2, nlo = 0, new = (4), ncl = 0, flag = 0 (0) s (-92.06,-0,-0,-92.06) -4.046e-06 (0,0) [0|4] (1) G (-188.1,-0,-0,188.1) -5.395e-06 (0,0) [0|4] (2) Z (108.4,-11.47,46,-34.56) 91.19 (0,0) [0|4] (3) W+ (98.63,33.83,-45.34,8.123) 80.39 (0,0) [0|4] (4) cb (73.06,-22.36,-0.6617,-69.55) -9.537e-07 (0,0) [0|4] } Singlet lists: Singlet parton list from CS_MCatNLO : 0x3ebe170 IS Parton 0x3d62a50 (1)[(0)]: sb : (92.06,0,0,92.06) (0,0)[0,0] k_T start : 36.33 k_T test : 0 k_T veto : 36.33(36.33) x_B : 0.02302 IS Parton 0x3f33ec0 (1)[(1)]: G : (188.1,0,0,-188.1) (0,0)[0,0] k_T start : 36.33 k_T test : 0 k_T veto : 36.33(36.33) x_B : 0.04701 FS Parton 0x39c4ef0 (1)[(2)]: Z : (108.4,-11.47,46,-34.56) (0,0)[0,0] k_T start : 36.33 k_T test : 0 k_T veto : 36.33(36.33) x_B : 1 FS Parton 0x3f2c720 (1)[(3)]: W+ : (98.63,33.83,-45.34,8.123) (0,0)[0,0] k_T start : 36.33 k_T test : 0 k_T veto : 36.33(36.33) x_B : 1 FS Parton 0x4153260 (1)[(4)]: cb : (73.06,-22.36,-0.6617,-69.55) (0,0)[0,0] k_T start : 36.33 k_T test : 0 k_T veto : 36.33(36.33) x_B : 1 ------------------------------------------------------------------------- } filling partner list for sb(0) ... G(1) cb(4) -> 2 dipole(s) filling partner list for G(1) ... sb(0) cb(4) -> 2 dipole(s) filling partner list for Z(2) ... -> 0 dipole(s) filling partner list for W+(3) ... -> 0 dipole(s) filling partner list for cb(4) ... sb(0) G(1) -> 2 dipole(s) before mc@nlo step Singlet parton list from CS_MCatNLO : 0x3ebe170 IS Parton 0x3d62a50 (1)[(0)]: sb : (92.06,0,0,92.06) (0,0)[0,0] k_T start : 36.33 k_T test : 0 k_T veto : 36.33(36.33) x_B : 0.02302 IS Parton 0x3f33ec0 (1)[(1)]: G : (188.1,0,0,-188.1) (0,0)[0,0] k_T start : 36.33 k_T test : 0 k_T veto : 36.33(36.33) x_B : 0.04701 FS Parton 0x39c4ef0 (1)[(2)]: Z : (108.4,-11.47,46,-34.56) (0,0)[0,0] k_T start : 36.33 k_T test : 0 k_T veto : 36.33(36.33) x_B : 1 FS Parton 0x3f2c720 (1)[(3)]: W+ : (98.63,33.83,-45.34,8.123) (0,0)[0,0] k_T start : 36.33 k_T test : 0 k_T veto : 36.33(36.33) x_B : 1 FS Parton 0x4153260 (1)[(4)]: cb : (73.06,-22.36,-0.6617,-69.55) (0,0)[0,0] k_T start : 36.33 k_T test : 0 k_T veto : 36.33(36.33) x_B : 1 ------------------------------------------------------------------------- Emission G -> sb s at kt = 21.56, z = 0.8756, y = 0.005879 for IS Parton 0x3d62a50 (1)[(0)]: sb : (92.06,0,0,92.06) (0,0)[0,0] k_T start : 36.33 k_T test : 21.56 k_T veto : 36.33(36.33) x_B : 0.02302 IS Parton 0x3f33ec0 (1)[(1)]: G : (188.1,0,0,-188.1) (0,0)[0,0] k_T start : 36.33 k_T test : 0 k_T veto : 36.33(36.33) x_B : 0.04701 Shower::MakeKinematics(0) { Singlet::ArrangeColours for sb --> G & s CS_Gamma::TrialWeight() { CS_Gamma::CalculateWeight(0x3f33d90) { (0x3d60bb0): 2 -> 4 { \mu_r = 90, \mu_f = 90, \mu_q = 90, \mu = 16.37 k_T = 21.56, z = 0, phi = 0, kin = 0, K = 0 oew = 2, oqcd = 2, nlo = 0, new = (5), ncl = 0, flag = 0 (0) G (-105.1,-0,-0,-105.1) -4.473e-06 (0,0) [0|4] (1) G (-188.1,-0,-0,188.1) -5.395e-06 (0,0) [0|4] (2) Z (106.9,-10.05,43.49,-33.43) 91.19 (0,0) [0|4] (3) W+ (100.4,35.34,-48.01,8.157) 80.39 (0,0) [0|4] (5) s (13.57,-3.655,6.463,11.36) -4.13e-07 (0,0) [0|4] (4) cb (72.34,-21.64,-1.941,-69) -1.349e-06 (0,0) [0|4] } Set fixed n_perp CS_Gamma::SingleWeight((0),(4)<->(1)) { B config -> (0x3edd8a0): 2 -> 3 { \mu_r = 90, \mu_f = 21.56, \mu_q = 0, \mu = 0 k_T = 0, z = 0, phi = 0, kin = 0, K = 0 oew = 0, oqcd = 0, nlo = 0, new = (none), ncl = 0, flag = 0 (0,4) s (-92.06,-0,-0,-92.06) -4.046e-06 (0,0) [0|0] <-> (1) (1) G (-188.1,-0,-0,188.1) -5.395e-06 (0,0) [0|0] (2) Z (108.4,-11.47,46,-34.56) 91.19 (0,0) [0|0] (3) W+ (98.63,33.83,-45.34,8.123) 80.39 (0,0) [0|0] (5) cb (73.06,-22.36,-0.6617,-69.55) -1.349e-06 (0,0) [0|0] } -> CS{kt=21.56,z=0.8756,phi=0.2707,mode=3,kin=1} ( G -> sb s ) add ( z = 0.8756, y = 0.005879, kt = 21.56 ) { (0) G (-105.1,-0,-0,-105.1) -4.473e-06 (0,0) [0|4] (4) s (13.57,-3.655,6.463,11.36) -4.13e-07 (0,0) [0|4] (1) G (-188.1,-0,-0,188.1) -5.395e-06 (0,0) [0|4] } -> w = 9.64e-06 * 0.00377 -> 3.635e-08 ( S = 1 ) } } Accumulate weights { [(0,5),(1)] -> 3.635e-08 2_3__sb__G__Z__W+__cb [ G -> sb s ] ( \mu_F = 21.56, \mu_R = 90 ) <- active } -> w = 3.635e-08 me / ecss = 3.095e-08 / 3.635e-08 = 0.8515 } w = 0.2838 -> accept 114:p_gamma->Weight()=3 } nem = 1 vs. maxem = 1 after mc@nlo step with 1 emission(s), w = 3 Singlet parton list from CS_MCatNLO : 0x3ebe170 IS Parton 0x4172f70 (1)[(0)]: G : (105.1,0,0,105.1) (0,601)[0,601] k_T start : 21.56 k_T test : 0 k_T veto : 36.33(36.33) x_B : 1 FS Parton 0x3edd730 (1)[(none)]: s : (13.57,-3.655,6.463,11.36) (601,0)[601,0] k_T start : 21.56 k_T test : 0 k_T veto : 36.33(36.33) x_B : 1 IS Parton 0x3f33ec0 (1)[(1)]: G : (188.1,0,0,-188.1) (0,0)[0,0] k_T start : 21.56 k_T test : 0 k_T veto : 36.33(36.33) x_B : 0.04701 FS Parton 0x39c4ef0 (1)[(2)]: Z : (106.9,-10.05,43.49,-33.43) (0,0)[0,0] k_T start : 21.56 k_T test : 0 k_T veto : 36.33(36.33) x_B : 1 FS Parton 0x3f2c720 (1)[(3)]: W+ : (100.4,35.34,-48.01,8.157) (0,0)[0,0] k_T start : 21.56 k_T test : 0 k_T veto : 36.33(36.33) x_B : 1 FS Parton 0x4153260 (1)[(4)]: cb : (72.34,-21.64,-1.941,-69) (0,0)[0,0] k_T start : 21.56 k_T test : 0 k_T veto : 36.33(36.33) x_B : 1 ------------------------------------------------------------------------- } CS_Cluster_Definitions::KPerp2(kin = 1, mode = 0) { Set fixed n_perp Kernel weight (A=1/NLO=0) [m=3,c=1] ( x = 0.02302 ) LF_VFF_II|CF_QCD { (0) G (-105.1,-0,-0,-105.1) -4.473e-06 (0,601) [0|4] (5) s (13.57,-3.655,6.463,11.36) -4.13e-07 (601,0) [0|4] (1) G (-188.1,-0,-0,188.1) -5.395e-06 (0,0) [0|4] } -> w = 0.001766 (566.3) } R selected via Sudakov (0x3f33d90): 2 -> 4 { \mu_r = 90, \mu_f = 90, \mu_q = 90, \mu = 21.56 k_T = 21.56, z = 0, phi = 0, kin = 0, K = 0 oew = 2, oqcd = 3, nlo = 0, new = (5), ncl = 0, flag = 0 (0) G (-105.1,-0,-0,-105.1) -4.473e-06 (0,601) [0|4] (1) G (-188.1,-0,-0,188.1) -5.395e-06 (0,0) [0|4] (2) Z (106.9,-10.05,43.49,-33.43) 91.19 (0,0) [0|4] (3) W+ (100.4,35.34,-48.01,8.157) 80.39 (0,0) [0|4] (5) s (13.57,-3.655,6.463,11.36) -4.13e-07 (601,0) [0|4] (4) cb (72.34,-21.64,-1.941,-69) -1.349e-06 (0,0) [0|4] } ( w = 3 ) } Signal_Processes::FillBlob(1) { Scale_Variations::ComputeVariations(2_4__G__G__Z__W+__s__cb) { } } Event_Handler::IterateEventPhases(): run 'Signal_Processes' -> Success Blob List with 1 elements { Blob [6]( 1, Signal Process , 2 -> 4 @ (0,0,0,0) Incoming particles : [G] 2 G 1 ( -> 1) [( 1.0515e+02, 0.0000e+00, 0.0000e+00, 1.0515e+02), p^2=-2.0009e-11, m= 0.0000e+00] ( 0, 0) [G] 2 G 2 ( -> 1) [( 1.8806e+02, 0.0000e+00, 0.0000e+00,-1.8806e+02), p^2=-2.9104e-11, m= 0.0000e+00] ( 0, 0) Outgoing particles : [H] 1 Z 3 ( 1 -> ) [( 1.0689e+02,-1.0049e+01, 4.3488e+01,-3.3426e+01), p^2= 8.3152e+03, m= 9.1188e+01] ( 0, 0) [H] 1 W+ 4 ( 1 -> ) [( 1.0041e+02, 3.5341e+01,-4.8010e+01, 8.1566e+00), p^2= 6.4617e+03, m= 8.0385e+01] ( 0, 0) [H] 1 s 5 ( 1 -> ) [( 1.3572e+01,-3.6549e+00, 6.4633e+00, 1.1360e+01), p^2=-1.7053e-13, m= 0.0000e+00] ( 0, 0) [H] 1 cb 6 ( 1 -> ) [( 7.2337e+01,-2.1638e+01,-1.9414e+00,-6.8997e+01), p^2=-1.8190e-12, m= 0.0000e+00] ( 0, 0) Data_Container: * Decay_Info ({NULL}) * Enhance (1) * Factorisation_Scale (8100) * MEWeight (1.805e-05) * MEWeightInfo (0x1deb518) * NLOEWType () * NLOQCDType () * Orders ((2,2)) * PDFInfo ( pdf1 = (-3,0.02302,90:0.04173) , pdf2 = (21,0.04701,90:2.281)) * Renormalization_Scale (8100) * ScaleVariations (0 variations) * Trials (1) * Weight (756.5) * Weight_Norm (140.3) } ME_Generator_Base::ShiftMasses(Comix) { Before shift: (0x3eaadf0): 2 -> 4 { \mu_r = 90, \mu_f = 90, \mu_q = 90, \mu = 21.56 k_T = 21.56, z = 0, phi = 0, kin = 0, K = 1.873 oew = 2, oqcd = 3, nlo = 1, new = (5), ncl = 0, flag = 0 (0) G (-105.1,-0,-0,-105.1) -4.473e-06 (0,601) [0|4] (1) G (-188.1,-0,-0,188.1) -5.395e-06 (0,0) [0|4] (2) Z (106.9,-10.05,43.49,-33.43) 91.19 (0,0) [0|4] (3) W+ (100.4,35.34,-48.01,8.157) 80.39 (0,0) [0|4] (5) s (13.57,-3.655,6.463,11.36) -4.13e-07 (601,0) [0|4] (4) cb (72.34,-21.64,-1.941,-69) -1.349e-06 (0,0) [0|4] } After shift: (0x3eaadf0): 2 -> 4 { \mu_r = 90, \mu_f = 90, \mu_q = 90, \mu = 21.56 k_T = 21.56, z = 0, phi = 0, kin = 0, K = 1.873 oew = 2, oqcd = 3, nlo = 1, new = (5), ncl = 0, flag = 0 (0) G (-105.1,0,0,-105.1) 1.349e-06 (0,601) [0|4] (1) G (-188.1,0,0,188.1) 3.815e-06 (0,0) [0|4] (2) Z (106.9,-10.05,43.48,-33.43) 91.19 (0,0) [0|4] (3) W+ (100.4,35.34,-48,8.153) 80.39 (0,0) [0|4] (5) s (13.57,-3.654,6.462,11.36) -1.686e-07 (601,0) [0|4] (4) cb (72.35,-21.63,-1.941,-68.99) 1.42 (0,0) [0|4] } } Cluster_Definitions_Base::ReCluster(0x17e5bf0) { (0x3eaadf0): 2 -> 4 { \mu_r = 90, \mu_f = 90, \mu_q = 90, \mu = 21.56 k_T = 21.56, z = 0, phi = 0, kin = 0, K = 1.873 oew = 2, oqcd = 3, nlo = 1, new = (5), ncl = 0, flag = 0 (0) G (-105.1,0,0,-105.1) 1.349e-06 (0,601) [0|4] (1) G (-188.1,0,0,188.1) 3.815e-06 (0,0) [0|4] (2) Z (106.9,-10.05,43.48,-33.43) 91.19 (0,0) [0|4] (3) W+ (100.4,35.34,-48,8.153) 80.39 (0,0) [0|4] (5) s (13.57,-3.654,6.462,11.36) -1.686e-07 (601,0) [0|4] (4) cb (72.35,-21.63,-1.941,-68.99) 1.42 (0,0) [0|4] } (0x3ebc410): 2 -> 3 { \mu_r = 90, \mu_f = 90, \mu_q = 90, \mu = 16.37 k_T = 36.33, z = 0, phi = 0, kin = 1, K = 1.873 oew = 2, oqcd = 2, nlo = 5, new = (4), ncl = 0, flag = 0 (0,5) s (-92.07,-0,-0,-92.07) 1.349e-06 (0,0) [1|4] <-> (1) (1) G (-188.1,0,0,188.1) 3.815e-06 (0,0) [1|4] (2) Z (108.4,-11.47,45.99,-34.56) 91.19 (0,0) [1|4] (3) W+ (98.62,33.82,-45.33,8.12) 80.39 (0,0) [1|4] (4) cb (73.07,-22.36,-0.6613,-69.55) 1.42 (0,0) [1|4] } (0x3d60bb0): 2 -> 2 { \mu_r = 90, \mu_f = 90, \mu_q = 204.1, \mu = 204.1 k_T = 204.1, z = 0, phi = 0, kin = 1, K = 1.873 oew = 2, oqcd = 1, nlo = 1, new = (none), ncl = 0, flag = 0 (0,5) s (-92.07,-0,-0,-92.07) 1.349e-06 (0,0) [1|4] (1,4) c (-113.2,-0,-0,113.1) 1.42 (0,0) [1|4] <-> (0,5) (2) Z (109.3,-23.04,45.65,-31.97) 91.19 (0,0) [1|4] (3) W+ (95.89,23.04,-45.65,10.9) 80.39 (0,0) [1|4] } } CS_Shower::PrepareStandardShower() { Exception_Handler::SignalHandler: Signal (11) caught.   Cannot continue. Exception_Handler::GenerateStackTrace(..): Generating stack trace  { 0x7f87139dd02d in 'ATOOLS::Exception_Handler::GenerateStackTrace(std::ostream&, bool, std::string const&)' (Exception_Handler.C:318) from '/scratch/sherpa/trunk/lib/SHERPA-MC/libToolsOrg.so.0' 0x7f87139de511 in 'ATOOLS::Exception_Handler::Terminate()' (basic_string.h:2689) from '/scratch/sherpa/trunk/lib/SHERPA-MC/libToolsOrg.so.0' 0x7f87139df0da in 'ATOOLS::Exception_Handler::SignalHandler(int)' (Exception_Handler.C:305) from '/scratch/sherpa/trunk/lib/SHERPA-MC/libToolsOrg.so.0' 0x7f87139df562 in 'ATOOLS::SignalHandler(int)' (Exception_Handler.C:244) from '/scratch/sherpa/trunk/lib/SHERPA-MC/libToolsOrg.so.0' 0x7f872a3b3e56 in 'CSSHOWER::CS_Shower::SetColours(ATOOLS::Cluster_Amplitude*)' (Cluster_Amplitude.H:94) from '/scratch/sherpa/trunk/lib/SHERPA-MC/libCSMain.so.0' 0x7f872a3b459f in 'CSSHOWER::CS_Shower::PrepareStandardShower(ATOOLS::Cluster_Amplitude*)' (CS_Shower.C:284 (discriminator 1)) from '/scratch/sherpa/trunk/lib/SHERPA-MC/libCSMain.so.0' 0x7f872bc87a85 in 'SHERPA::Perturbative_Interface::DefineInitialConditions(ATOOLS::Blob*)' (Perturbative_Interface.C:205) from '/scratch/sherpa/trunk/lib/SHERPA-MC/libSherpaPerturbativePhysics.so.0' 0x7f872120600d in 'SHERPA::Jet_Evolution::AttachShowers(ATOOLS::Blob*, ATOOLS::Blob_List*, SHERPA::Perturbative_Interface*)' (Jet_Evolution.C:156) from '/scratch/sherpa/trunk/lib/SHERPA-MC/libSherpaSingleEvents.so.0' 0x7f8721207155 in 'SHERPA::Jet_Evolution::Treat(ATOOLS::Blob_List*, double&)' (Jet_Evolution.C:110) from '/scratch/sherpa/trunk/lib/SHERPA-MC/libSherpaSingleEvents.so.0' 0x7f87211f2eca in 'SHERPA::Event_Handler::IterateEventPhases(SHERPA::eventtype::code&, double&)' (Event_Handler.C:192) from '/scratch/sherpa/trunk/lib/SHERPA-MC/libSherpaSingleEvents.so.0' 0x7f87211f39c6 in 'SHERPA::Event_Handler::GenerateStandardPerturbativeEvent(SHERPA::eventtype::code&)' (Event_Handler.C:260) from '/scratch/sherpa/trunk/lib/SHERPA-MC/libSherpaSingleEvents.so.0' 0x7f87211f5b0d in 'SHERPA::Event_Handler::GenerateEvent(SHERPA::eventtype::code)' (Event_Handler.C:142) from '/scratch/sherpa/trunk/lib/SHERPA-MC/libSherpaSingleEvents.so.0' 0x7f872c0f0234 in 'SHERPA::Sherpa::GenerateOneEvent(bool)' (Sherpa.C:217) from '/scratch/sherpa/trunk/lib/SHERPA-MC/libSherpaMain.so.0' 0x4090cd in 'main' (Main.C:27) from '/scratch/sherpa/trunk/bin/Sherpa' } Exception_Handler::Terminate(): Pre-crash status saved to '/home/fsiegert/sherpa/tests/libs/trunk/wzj/Status__Thu_Nov_12_11-31-39_2015'. Exception_Handler::ApproveTerminate(): Asking for termination ... ... approved. Exception_Handler::PrepareTerminate(): Preparing termination ... Event 43 ( 0 s total ) = 1.615e+07 evts/day In Event_Handler::Finish : Summarizing the run may take some time. +------------------------------------------------------+ | | | Total XS is 7.19097 pb +- ( 3.16245 pb = 43.97 % ) | | | +------------------------------------------------------+ ... prepared. Exception_Handler::Exit: Exiting Sherpa with code (2)[?25h