Welcome to Sherpa, Frank Siegert. Initialization of framework underway. The local time is Tue Feb 9 09:29:42 2016. Run_Parameter::Init(): Seed mode '*' Run_Parameter::Init(): Setting memory limit to 62.5894 GB. Random::SetSeed(): Seed set to 652413 ----------------------------------------------------------------------------- ----------- 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, 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 25839. 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 2.4952 0 1 1 91.1876 W+ 24 80.385 2.085 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} f 92 {d,db,u,ub,s,sb,c,cb,b,bb,e-,e+,mu-,mu+,tau-,tau+,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} Hadron_Init::Init(): Initializing kf table for hadrons. Initialized the Fragmentation_Handler. Initialized the Soft_Collision_Handler. CS_Shower::CS_Shower(): Set respect Q2 mode 0 CS_Shower::CS_Shower(): Set ME weight mode 0 CS_Shower::CS_Shower(): Set color setter mode 0 CS_Shower::CS_Shower(): Set respect Q2 mode 0 CS_Shower::CS_Shower(): Set ME weight mode 0 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_Generators::ME_Generators(): Try loading 'OpenLoops' from 'libSherpaOpenLoops'. Matrix_Element_Handler::Matrix_Element_Handler(): Set seed mode 3. Matrix_Element_Handler::Matrix_Element_Handler(): Set seed increment to 1 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. Initialising OpenLoops generator from /scratch/p_lhc/OpenLoops ######################################################### # ___ version 1.2.3 # # / \ ___ ____ _ _ | __ __ ___ __ # # | | |__| |__ |\ | | / \ / \ |__| /__ # # \___/ | |___ | \| |___ \__/ \__/ | __/ # # # ######################################################### # You are using OpenLoops to evaluate loop amplitudes # # Authors: # # F. Cascioli, J. Lindert, P. Maierhoefer, S. Pozzorini # # Please cite Phys. Rev. Lett. 108 (2012) 111601 # ######################################################### 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 ME_Generators::LoadGenerator(): Try loading 'gg_yy' from 'libSherpagg_yy'. Library_Loader::LoadLibrary(): Failed to load library 'libSherpagg_yy.so'. ME_Generators::LoadGenerator(): ME generator 'gg_yy' not found. Ignoring it. ................................................................................................................................................................................................................................................................................................................................................................................................................................................................................ done ( 193 MB, 12s / 10s ). Matrix_Element_Handler::InitializeProcesses(): Performing tests .................................................................................................................................................................................................................................................................................... Amplitude::GaugeTest(): Large deviation { -3.92408696729e-15 vs -3.92408689264e-15 => 1.90230162733e-08 } .... done ( 205 MB, 0s / 0s ). Initialized the Matrix_Element_Handler for the hard processes. Initialized the Beam_Remnant_Handler. Hadron_Decay_Map::Read: Initializing HadronDecays.dat. This may take some time. Initialized the Hadron_Decay_Handler, Decay model = Hadrons ME_Generator_Base::SetPSMasses(): Massive PS flavours for Amisic: (c,cb,b,bb,e-,e+,mu-,mu+,tau-,tau+) Simple_Chain::InitializeProcessList(): Init processes ................................................................................................................................................................................ done. Grid_Creator::ReadInGrid(): Reading grid ................................................................................................................................................................................ done. Simple_Chain::CalculateTotal(): Result is { \sigma_{hard} = 35.0108 mb at PT_{min} = 4.16597 GeV } Profile_Function_Base::CalculateOMean(2.23852): Results are { k = 6.68402 <\tilde{O}> = 0.334906 } Initialized the Multiple_Interactions_Handler (MI_Handler). Initialized the Soft_Photon_Handler. Read in channels from directory : .//Results/Internal/MC_2_2__G__G__P__P Single_Process::CalculateTotalXSec(): Calculate xs for '2_2__G__G__P__P' (Internal) 2_2__G__G__P__P : 60.3982 pb +- ( 0.371793 pb = 0.61557 % ) exp. eff: 4.55959 % Read in channels from directory : .//Results/Amegic/MC_2_2__j__j__P__P__QCD(BVI) Process_Group::CalculateTotalXSec(): Calculate xs for '2_2__j__j__P__P__QCD(BVI)' (Amegic) 2_2__j__j__P__P__QCD(BVI) : 148.372 pb +- ( 1.092 pb = 0.735987 % ) exp. eff: 0.835688 % Read in channels from directory : .//Results/Comix/MC_2_3__j__j__P__P__j__QCD(RS) Process_Group::CalculateTotalXSec(): Calculate xs for '2_3__j__j__P__P__j__QCD(RS)' (Comix) 2_3__j__j__P__P__j__QCD(RS) : 185.317 pb +- ( 3.0772 pb = 1.6605 % ) exp. eff: 0.0600374 % reduce max for 2_2__j__j__P__P__QCD(BVI) to 0.232494 ( eps = 0.001 ) reduce max for 2_3__j__j__P__P__j__QCD(RS) to 0.728028 ( eps = 0.001 ) Read in channels from directory : .//Results/Amegic/MC_2_3__j__j__P__P__j__QCD(BVI) Process_Group::CalculateTotalXSec(): Calculate xs for '2_3__j__j__P__P__j__QCD(BVI)' (Amegic) 2_3__j__j__P__P__j__QCD(BVI) : 182.662 pb +- ( 22.2792 pb = 12.197 % ) exp. eff: 0.0133373 % Read in channels from directory : .//Results/Comix/MC_2_4__j__j__P__P__j__j__QCD(RS) Process_Group::CalculateTotalXSec(): Calculate xs for '2_4__j__j__P__P__j__j__QCD(RS)' (Comix) 2_4__j__j__P__P__j__j__QCD(RS) : 650.764 pb +- ( 48.7673 pb = 7.49386 % ) exp. eff: 0.000983882 % reduce max for 2_3__j__j__P__P__j__QCD(BVI) to 0.621321 ( eps = 0.001 ) reduce max for 2_4__j__j__P__P__j__j__QCD(RS) to 0.234216 ( eps = 0.001 ) Read in channels from directory : .//Results/Comix/MC_2_4__j__j__P__P__j__j Process_Group::CalculateTotalXSec(): Calculate xs for '2_4__j__j__P__P__j__j' (Comix) 2_4__j__j__P__P__j__j : 505.694 pb +- ( 101.867 pb = 20.1439 % ) exp. eff: 0.000285122 % reduce max for 2_4__j__j__P__P__j__j to 0.475752 ( eps = 0.001 ) Read in channels from directory : .//Results/Comix/MC_2_5__j__j__P__P__j__j__j Process_Group::CalculateTotalXSec(): Calculate xs for '2_5__j__j__P__P__j__j__j' (Comix) 2_5__j__j__P__P__j__j__j : 538.893 pb +- ( 12.8658 pb = 2.38746 % ) exp. eff: 0.00106009 % reduce max for 2_5__j__j__P__P__j__j__j to 0.839089 ( eps = 0.001 ) Exception_Handler::ReadInStatus(): Reading status from './' { Random::ReadInStatus(): Reading status from './random.dat'. Simple_Chain::ReadInStatus(): Reading status from './MIG_P+P+_8000_NNPDF_2/'. Grid_Creator::ReadInGrid(): Reading grid ................................................................................................................................................................................ done. } ---------------------------------------------------------- -- SHERPA generates events with the following structure -- ---------------------------------------------------------- Perturbative : Signal_Processes Perturbative : Hard_Decays Perturbative : Jet_Evolution:CSS Perturbative : Lepton_FS_QED_Corrections:Photons Perturbative : Multiple_Interactions:Amisic Perturbative : Minimum_Bias:Off Hadronization : Beam_Remnants Hadronization : Hadronization:Ahadic Hadronization : Hadron_Decays Analysis : Rivet --------------------------------------------------------- Event_Handler::GenerateEvent(0) { MCatNLO_Process::OneEvent() { BBar_Emission_Generator::GeneratePoint() { in EEG: #1 = 0.22155, #2 = 0.990825, #3 = 0.671524 Born point { (431.728,0,0,431.728) (1.88859,0,0,-1.88859) (70.3677,-18.7309,8.7491,67.2624) (38.0324,-0.266974,-0.0934536,38.0313) (325.217,18.9978,-8.65565,324.546) } selected IF_1[1,5]_4[4] IF_Dipole::GeneratePoint() { vegased : x = 0.430283, u = 0.976596, phi = 0.748703, xmin = 0.000472148 transformed : x = 0.19595, u = 0.953744, phi = 4.70424 } } } MCatNLO_Process::OneEvent() { BBar_Emission_Generator::GeneratePoint() { in EEG: #1 = 0.356341, #2 = 0.137038, #3 = 0.680821 Born point { (1015.92,0,0,1015.92) (2.33526,0,0,-2.33526) (873.486,-3.65423,3.01483,873.473) (124.685,-7.11788,3.43326,124.434) (20.086,10.7721,-6.44808,15.679) } selected II_0[0,5]_1[1] II_Dipole::GeneratePoint() { vegased : x = 0.503793, v = 0.0920044, phi = 0.728303, xmin = 0.25398 transformed : x = 0.568311, v = 0.00366515, phi = 4.57606 II_Dipole::Construct() { (0):(1015.92,0,0,1015.92) -> (0):(1787.62,0,0,1787.62) (5):(765.151,-5.07068,-0.695583,765.134) (1):(2.33526,0,0,-2.33526) -> (1):(2.33526,0,0,-2.33526) (2):(873.486,-3.65423,3.01483,873.473) -> (2):(876.752,-1.46629,3.31496,876.745) (3):(124.685,-7.11788,3.43326,124.434) -> (3):(121.814,-6.6795,3.49339,121.581) (4):(20.086,10.7721,-6.44808,15.679) -> (4):(26.2338,13.2165,-6.11277,21.8213) } } } } MCatNLO_Process::OneEvent() { BBar_Emission_Generator::GeneratePoint() { in EEG: #1 = 0.554378, #2 = 0.041109, #3 = 0.7447 Born point { (0.693185,0,0,0.693185) (1641.84,0,0,-1641.84) (997.198,-13.1315,16.7153,-996.971) (387.209,0.647027,4.03222,-387.188) (258.121,12.4845,-20.7475,-256.983) } selected II_0[0,4]_1[1] II_Dipole::GeneratePoint() { vegased : x = 0.527202, v = 0.0299543, phi = 0.724697, xmin = 0.000173296 transformed : x = 0.284544, v = 0.000646874, phi = 4.55341 II_Dipole::Construct() { (0):(0.693185,0,0,0.693185) -> (0):(2.43613,0,0,2.43613) (4):(2.80343,-2.68558,-0.430594,0.679306) (1):(1641.84,0,0,-1641.84) -> (1):(1641.84,0,0,-1641.84) (2):(997.198,-13.1315,16.7153,-996.971) -> (2):(986.137,-12.1053,16.8798,-985.919) (3):(387.209,0.647027,4.03222,-387.188) -> (3):(388.064,0.985368,4.08647,-388.042) (4):(258.121,12.4845,-20.7475,-256.983) -> (5):(267.266,13.8055,-20.5357,-266.118) } } } } MCatNLO_Process::OneEvent() { Single_Process::Differential(2_4__u__ub__P__P__G__G__QCD(RS), RS:0x43c9b68) { } } MCatNLO_Process::OneEvent() { BBar_Emission_Generator::GeneratePoint() { in EEG: #1 = 0.987001, #2 = 0.897566, #3 = 0.520969 Born point { (20.1913,0,0,20.1913) (1258.28,0,0,-1258.28) (43.9602,20.7706,9.66198,-37.5197) (1049.14,-108.3,-50.3786,-1042.32) (185.372,87.5295,40.7166,-158.251) } selected II_1[1,5]_0[0] II_Dipole::GeneratePoint() { vegased : x = 0.985916, v = 0.84617, phi = 0.350517, xmin = 0.31457 transformed : x = 0.987668, v = 0.00883238, phi = 2.20236 II_Dipole::Construct() { (1):(1258.28,0,0,-1258.28) -> (1):(1273.99,0,0,-1273.99) (5):(4.63636,1.4393,-1.05287,-4.27968) (0):(20.1913,0,0,20.1913) -> (0):(20.1913,0,0,20.1913) (2):(43.9602,20.7706,9.66198,-37.5197) -> (2):(44.0481,20.633,9.76258,-37.6724) (3):(1049.14,-108.3,-50.3786,-1042.32) -> (3):(1059.75,-109.023,-49.8501,-1052.95) (4):(185.372,87.5295,40.7166,-158.251) -> (4):(185.744,86.9502,41.1404,-158.896) } } } Jet_Finder::Trigger(): '2_3__G__u__P__P__u__QCD(BVI)' Q_cut = 9.92952 { CS_Shower::JetVeto(mode = 0) { (0x68f0ab0): 2 -> 3 { \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 = (), ncl = 0, flag = 0 (0) G (-20.1913,-0,-0,-20.1913) -3.36331e-06 (-1,-1) [0|0] (1) ub (-1258.28,-0,-0,1258.28) 0 (-1,-1) [0|0] (2) P (43.9602,20.7706,9.66198,-37.5197) 0 (-1,-1) [0|0] (3) P (1049.14,-108.3,-50.3786,-1042.32) -2.15792e-05 (-1,-1) [0|0] (4) u (185.372,87.5295,40.7166,-158.251) 3.8147e-06 (-1,-1) [0|0] } noem = (), nospec = () No kernel for G P <-> u (22) No kernel for G P <-> P (21) No kernel for G P <-> u (21) No kernel for G P <-> u (22) No kernel for G P <-> P (21) No kernel for G P <-> u (21) Q_{(0)(4),(1)} = 86.588 No kernel for G u <-> P (21) No kernel for G u <-> P (21) No kernel for u P <-> G (22) No kernel for u P <-> P (21) No kernel for u P <-> u (21) No kernel for u P <-> G (22) No kernel for u P <-> P (21) No kernel for u P <-> u (21) Q_{(1)(4),(0)} = 242.835 No kernel for u u <-> P (21) No kernel for u u <-> P (21) No kernel for P P <-> G (12) No kernel for P P <-> u (12) No kernel for P P <-> u (11) No kernel for P u <-> G (12) No kernel for P u <-> u (12) No kernel for P u <-> P (11) No kernel for P u <-> G (12) No kernel for P u <-> u (12) No kernel for P u <-> P (11) --- Jet veto --- } } -> 1 BBar_Emission_Generator::GenerateWeight() { Dipole II_1[1,5]_0[0] { (2):(44.0481,20.633,9.76258,-37.6724) -> (2):(43.9602,20.7706,9.66198,-37.5197) (3):(1059.75,-109.023,-49.8501,-1052.95) -> (3):(1049.14,-108.3,-50.3786,-1042.32) (4):(185.744,86.9502,41.1404,-158.896) -> (4):(185.372,87.5295,40.7166,-158.251) again : x = 0.987668, v = 0.00883238, phi = 2.20236 recovered : x = 0.985916, v = 0.84617, phi = 0.350517, xmin = 0.31457 } -> w = 33.5084 ( a = 0.1 ) } Single_Process::Differential(2_3__G__u__P__P__u__QCD(BVI), RS:0) { Scale_Setter_Base::CalculateScale(2_3__G__u__P__P__u__QCD(BVI)) { METS_Scale_Setter::CalculateStrict(2_3__G__u__P__P__u__QCD(BVI) from 2_3__G__u__P__P__u__QCD(BVI), mcmode = 0) { Cluster_Algorithm::Cluster() { Combine_Table::CalcJet(mode = 512, nl = 5, complete = 1) { Combine_Table (0xc616760) 1 [OQCD=1] (up=#) idcontent flav cut qcd qed mom 0 (0) G 0 0/0 0/0 (20.1913,0,0,20.1913) 1 (1) u 0 0/0 0/0 (1258.28,0,0,-1258.28) 2 (2) P 0 0/0 0/0 (43.9602,20.7706,9.66198,-37.5197) 3 (3) P 0 0/0 0/0 (1049.14,-108.3,-50.3786,-1042.32) 4 (4) u 0 0/0 0/0 (185.372,87.5295,40.7166,-158.251) --------------- 0&4%1[ub ] CP{kt=0,op=0,x=0,mu=0,k=0,m=0} 0 2,5 0&4%2[ub ] CP{kt=0,op=0,x=0,mu=0,k=0,m=0} 0 5 0&4%3[ub ] CP{kt=0,op=0,x=0,mu=0,k=0,m=0} 0 2 1&2%0[u ] CP{kt=0,op=0,x=0,mu=0,k=0,m=0} 0 1 1&2%3[u ] CP{kt=0,op=0,x=0,mu=0,k=0,m=0} 0 1,2 1&2%4[u ] CP{kt=0,op=0,x=0,mu=0,k=0,m=0} 0 1,2 1&3%0[u ] CP{kt=0,op=0,x=0,mu=0,k=0,m=0} 0 4 1&3%2[u ] CP{kt=0,op=0,x=0,mu=0,k=0,m=0} 0 4,5 1&3%4[u ] CP{kt=0,op=0,x=0,mu=0,k=0,m=0} 0 4,5 2&4%0[u ] CP{kt=0,op=0,x=0,mu=0,k=0,m=0} 0 4 2&4%1[u ] CP{kt=0,op=0,x=0,mu=0,k=0,m=0} 0 3,4 2&4%3[u ] CP{kt=0,op=0,x=0,mu=0,k=0,m=0} 0 3,4 3&4%0[u ] CP{kt=0,op=0,x=0,mu=0,k=0,m=0} 0 1 3&4%1[u ] CP{kt=0,op=0,x=0,mu=0,k=0,m=0} 0 0,1 3&4%2[u ] CP{kt=0,op=0,x=0,mu=0,k=0,m=0} 0 0,1 k_{T,min} (0,1,2,3,4) -> 0 CS_Cluster_Definitions::KPerp2(kin = -1, mode = 16) { Set fixed n_perp Kernel weight (A=0/NLO=16) [m=3,c=1] ( x = 0.000968497 ) LF_VFF_II|CF_QCD { (0) G (-20.1913,-0,-0,-20.1913) -3.36331e-06 (-1,-1) [0|0] (4) u (185.372,87.5295,40.7166,-158.251) 3.8147e-06 (-1,-1) [0|0] (1) ub (-1258.28,-0,-0,1258.28) 0 (-1,-1) [0|0] } -> w = 0.00023822 (4197.8) } Calculate m_perp(0[G],4[u],1[u],ub) -> CP{kt=117.798,op=64.7905,x=0.191864,mu=53.0834,k=1,m=0} CS_Cluster_Definitions::KPerp2(kin = -1, mode = 16) { Invalid spectator P } Calculate m_perp(0[G],4[u],2[P],ub) -> CP{kt=117.798,op=-1,x=1,mu=117.798,k=1,m=0} CS_Cluster_Definitions::KPerp2(kin = -1, mode = 16) { Invalid spectator P } Calculate m_perp(0[G],4[u],3[P],ub) -> CP{kt=117.798,op=-1,x=0.164823,mu=117.798,k=1,m=0} CS_Cluster_Definitions::KPerp2(kin = -1, mode = 16) { Set fixed n_perp Invalid spectator G } Calculate m_perp(1[u],2[P],0[G],u) -> CP{kt=55.7648,op=-1,x=0.808136,mu=55.7648,k=1,m=0} CS_Cluster_Definitions::KPerp2(kin = -1, mode = 16) { Invalid spectator P } Calculate m_perp(1[u],2[P],3[P],u) -> CP{kt=97.3089,op=-1,x=0.415776,mu=97.3089,k=1,m=0} CS_Cluster_Definitions::KPerp2(kin = -1, mode = 16) { No Kernel. Set weight -1. } Calculate m_perp(1[u],2[P],4[u],u) -> CP{kt=0.0154541,op=-1,x=1,mu=1.41421,k=1,m=0} CS_Cluster_Definitions::KPerp2(kin = -1, mode = 16) { Set fixed n_perp Invalid spectator G } Calculate m_perp(1[u],3[P],0[G],u) -> CP{kt=131.022,op=-1,x=1.22463e-08,mu=131.022,k=1,m=0} CS_Cluster_Definitions::KPerp2(kin = -1, mode = 16) { Invalid spectator P } Calculate m_perp(1[u],3[P],2[P],u) -> CP{kt=100.146,op=-1,x=0.415776,mu=100.146,k=1,m=0} CS_Cluster_Definitions::KPerp2(kin = -1, mode = 16) { No Kernel. Set weight -1. } Calculate m_perp(1[u],3[P],4[u],u) -> CP{kt=128.474,op=-1,x=0.0385211,mu=128.474,k=1,m=0} CS_Cluster_Definitions::KPerp2(kin = -1, mode = 16) { Invalid spectator G } Calculate m_perp(2[P],4[u],0[G],u) -> CP{kt=0.0154448,op=-1,x=1,mu=0.0154448,k=1,m=0} CS_Cluster_Definitions::KPerp2(kin = -1, mode = 16) { No Kernel. Set weight -1. } Calculate m_perp(2[P],4[u],1[u],u) -> CP{kt=0.0154541,op=-1,x=1,mu=1,k=1,m=0} CS_Cluster_Definitions::KPerp2(kin = -1, mode = 16) { Invalid spectator P } Calculate m_perp(2[P],4[u],3[P],u) -> CP{kt=0.0154526,op=-1,x=1,mu=0.0154526,k=1,m=0} CS_Cluster_Definitions::KPerp2(kin = -1, mode = 16) { Invalid spectator G } Calculate m_perp(3[P],4[u],0[G],u) -> CP{kt=265.589,op=-1,x=0.164823,mu=265.589,k=1,m=0} CS_Cluster_Definitions::KPerp2(kin = -1, mode = 16) { No Kernel. Set weight -1. } Calculate m_perp(3[P],4[u],1[u],u) -> CP{kt=128.474,op=-1,x=0.0385211,mu=128.474,k=1,m=0} CS_Cluster_Definitions::KPerp2(kin = -1, mode = 16) { Invalid spectator P } Calculate m_perp(3[P],4[u],2[P],u) -> CP{kt=286.578,op=-1,x=1,mu=286.578,k=1,m=0} Combine_Table (0xc616760) 1 [OQCD=1] (up=#) idcontent flav cut qcd qed mom 0 (0) G 0 0/0 0/0 (20.1913,0,0,20.1913) 1 (1) u 0 0/0 0/0 (1258.28,0,0,-1258.28) 2 (2) P 0 0/0 0/0 (43.9602,20.7706,9.66198,-37.5197) 3 (3) P 0 0/0 0/0 (1049.14,-108.3,-50.3786,-1042.32) 4 (4) u 0 0/0 0/0 (185.372,87.5295,40.7166,-158.251) --------------- 0&4%1[ub ] CP{kt=117.798,op=64.7905,x=0.191864,mu=53.0834,k=1,m=0} 0 2,5 <- 0&4%2[ub ] CP{kt=117.798,op=-1,x=1,mu=117.798,k=1,m=0} 0 5 0&4%3[ub ] CP{kt=117.798,op=-1,x=0.164823,mu=117.798,k=1,m=0} 0 2 1&2%0[u ] CP{kt=55.7648,op=-1,x=0.808136,mu=55.7648,k=1,m=0} 0 1 1&2%3[u ] CP{kt=97.3089,op=-1,x=0.415776,mu=97.3089,k=1,m=0} 0 1,2 1&2%4[u ] CP{kt=0.0154541,op=-1,x=1,mu=1.41421,k=1,m=0} 0 1,2 1&3%0[u ] CP{kt=131.022,op=-1,x=1.22463e-08,mu=131.022,k=1,m=0} 0 4 1&3%2[u ] CP{kt=100.146,op=-1,x=0.415776,mu=100.146,k=1,m=0} 0 4,5 1&3%4[u ] CP{kt=128.474,op=-1,x=0.0385211,mu=128.474,k=1,m=0} 0 4,5 2&4%0[u ] CP{kt=0.0154448,op=-1,x=1,mu=0.0154448,k=1,m=0} 0 4 2&4%1[u ] CP{kt=0.0154541,op=-1,x=1,mu=1,k=1,m=0} 0 3,4 2&4%3[u ] CP{kt=0.0154526,op=-1,x=1,mu=0.0154526,k=1,m=0} 0 3,4 3&4%0[u ] CP{kt=265.589,op=-1,x=0.164823,mu=265.589,k=1,m=0} 0 1 3&4%1[u ] CP{kt=128.474,op=-1,x=0.0385211,mu=128.474,k=1,m=0} 0 0,1 3&4%2[u ] CP{kt=286.578,op=-1,x=1,mu=286.578,k=1,m=0} 0 0,1 k_{T,min} (0,1,2,3,4) -> 0 set last k_T = 117.798 for (0,1,2,3,4) from (0,4) set last k_T = 117.798 for (0,1,2,3,4) from (0,4) Combine_Table::CalcJet(mode = 512, nl = 4, complete = 1) { Combine_Table (0xc5fc4b0) 2 [OQCD=0] (up=1) idcontent flav cut qcd qed mom 0 (0,4) ub 0 1/3 0/0 (3.87399,0,0,3.87399) 1 (1) u 0 0/0 0/0 (1258.28,0,0,-1258.28) 2 (2) P 0 0/0 0/0 (346.729,56.4247,26.2474,-341.099) 3 (3) P 0 0/0 0/0 (915.425,-56.4247,-26.2474,-913.308) --------------- k_{T,min} (0,1,2,3,4) -> 117.798 Default_Core_Scale::Calculate() { (0xa5fbeb0): 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 = (), ncl = 0, flag = 0 (0,4) u (-3.87399,-0,-0,-3.87399) -6.46099e-07 (-1,-1) [0|0] (1) ub (-1258.28,-0,-0,1258.28) 0 (-1,-1) [0|0] (2) P (346.729,56.4247,26.2474,-341.099) 0 (-1,-1) [0|0] (3) P (915.425,-56.4247,-26.2474,-913.308) -2.15792e-05 (-1,-1) [0|0] } DY like \mu_f = 139.636 \mu_r = 139.636 \mu_q = 139.636 } } } Cluster_Algorithm::Convert(): { Default_Core_Scale::Calculate() { (0xc616930): 2 -> 2 { \mu_r = 8000, \mu_f = 8000, \mu_q = 8000, \mu = 0 k_T = 0, z = 0, phi = 0, kin = 1, K = 0 oew = 2, oqcd = 1, nlo = 0, new = (), ncl = 0, flag = 0 (0,4) u (-3.87399,-0,-0,-3.87399) -6.46099e-07 (0,0) [1|0] <-> (1) (1) ub (-1258.28,-0,-0,1258.28) 0 (0,0) [1|0] (2) P (346.729,56.4247,26.2474,-341.099) 0 (0,0) [1|0] (3) P (915.425,-56.4247,-26.2474,-913.308) -2.15792e-05 (0,0) [1|0] } DY like \mu_f = 139.636 \mu_r = 139.636 \mu_q = 139.636 } (0xc616930): 2 -> 2 { \mu_r = 8000, \mu_f = 8000, \mu_q = 139.636, \mu = 139.636 k_T = 139.636, z = 0, phi = 0, kin = 1, K = 0 oew = 2, oqcd = 1, nlo = 0, new = (), ncl = 0, flag = 0 (0,4) u (-3.87399,-0,-0,-3.87399) -6.46099e-07 (0,0) [1|6] <-> (1) (1) ub (-1258.28,-0,-0,1258.28) 0 (0,0) [1|6] (2) P (346.729,56.4247,26.2474,-341.099) 0 (0,0) [1|6] (3) P (915.425,-56.4247,-26.2474,-913.308) -2.15792e-05 (0,0) [1|6] } (0xa5fbeb0): 2 -> 3 { \mu_r = 8000, \mu_f = 8000, \mu_q = 8000, \mu = 53.0834 k_T = 117.798, z = 0, phi = 0, kin = 0, K = 0 oew = 2, oqcd = 2, nlo = 0, new = (4), ncl = 0, flag = 0 (0) G (-20.1913,-0,-0,-20.1913) -3.36331e-06 (0,0) [0|6] (1) ub (-1258.28,-0,-0,1258.28) 0 (0,0) [0|6] (2) P (43.9602,20.7706,9.66198,-37.5197) 0 (0,0) [0|6] (3) P (1049.14,-108.3,-50.3786,-1042.32) -2.15792e-05 (0,0) [0|6] (4) u (185.372,87.5295,40.7166,-158.251) 3.8147e-06 (0,0) [0|6] } } } Setting scales { \mu_{2} = 1 * 1 * 53.0834, as = 0.128528, O(QCD) = 1 } -> as = 0.128528 -> 53.0834 Core / QCD scale = 139.636 / 53.0834 METS_Scale_Setter::SetScales(): Set { \mu_f = 139.636 \mu_r = 53.0834 \mu_q = 139.636 \mu_3 = 139.636 \mu_4 = 53.0834 \mu_5 = 117.798 } <- 2_3__G__u__P__P__u__QCD(BVI) } Coupling_Data::Calculate(0x25c04e0): scl = 53.0834 -> 'Alpha_QCD'{fac=1.08922,cpl=0.118} } Setting OpenLoops parameter: alpha = 0.00776384 Setting OpenLoops parameter: alphas = 0.128528 Setting OpenLoops parameter: mu = 53.0834 OpenLoops_Virtual::Calc:("43.9602 (G,u,P,P,u) user=0.05 real=0.14 sys=0.01") Single_Process::NLOCounterTerms(2_3__G__u__P__P__u__QCD(BVI)) { \mu_F = 139.636 -> 1 \mu_R = 53.0834 -> 1 METOOLS::AlphaSCounterTerm(2817.85 -> 2817.85) { } METOOLS::CollinearCounterTerms(Q = 139.636 / 139.636) { } METOOLS::CollinearCounterTerms(Q = 139.636 / 139.636) { } Single_Process::CollinearCounterTerms(Q = 139.636 / 139.636) { } 0 Single_Process::CollinearCounterTerms(Q = 139.636 / 139.636) { } 0 C = 0 } PDF(fla=G, xa=0.00504782, ta=19498.3) * PDF(flb=u, xb=0.31457, tb=19498.3) -> 2467.89 Single_Process::BeamISRWeight(2_3__G__u__P__P__u__QCD(BVI), \mu_F = 139.636, mode = 0, #ampls=1) { * [ PDF(fla=G, xa=0.00504782, ta=13876.4) * PDF(flb=u, xb=0.31457, tb=13876.4) -> 2481.39 / PDF(fla=G, xa=0.00504782, ta=19498.3) * PDF(flb=u, xb=0.31457, tb=19498.3) -> 2467.89 ] = 1.00547 Single_Process::CollinearCounterTerms(Q = 117.798 / 139.636) { \mu_F = 139.636 \mu_R = 53.0834 0.128528/(2\pi) * log(117.798/139.636) = -0.00695781 Beam 0: z = 0.9389, f_{G}(0.00504782,139.636) = 2512.51 { P_{d,G}(0.9389) = {F=1.4254,S=0,I=0}, f_{d}(0.00537632,139.636) = 161.838, f_{d}(0.00504782,139.636) = 176.404 P_{db,G}(0.9389) = {F=1.4254,S=0,I=0}, f_{db}(0.00537632,139.636) = 144.707, f_{db}(0.00504782,139.636) = 158.654 P_{u,G}(0.9389) = {F=1.4254,S=0,I=0}, f_{u}(0.00537632,139.636) = 173.793, f_{u}(0.00504782,139.636) = 188.82 P_{ub,G}(0.9389) = {F=1.4254,S=0,I=0}, f_{ub}(0.00537632,139.636) = 141.207, f_{ub}(0.00504782,139.636) = 155.074 P_{s,G}(0.9389) = {F=1.4254,S=0,I=0}, f_{s}(0.00537632,139.636) = 123.264, f_{s}(0.00504782,139.636) = 135.846 P_{sb,G}(0.9389) = {F=1.4254,S=0,I=0}, f_{sb}(0.00537632,139.636) = 124.009, f_{sb}(0.00504782,139.636) = 136.639 P_{c,G}(0.9389) = {F=1.4254,S=0,I=0}, f_{c}(0.00537632,139.636) = 101.602, f_{c}(0.00504782,139.636) = 112.12 P_{cb,G}(0.9389) = {F=1.4254,S=0,I=0}, f_{cb}(0.00537632,139.636) = 101.692, f_{cb}(0.00504782,139.636) = 112.211 P_{b,G}(0.9389) = {F=1.4254,S=0,I=0}, f_{b}(0.00537632,139.636) = 72.9679, f_{b}(0.00504782,139.636) = 80.8105 P_{bb,G}(0.9389) = {F=1.4254,S=0,I=0}, f_{bb}(0.00537632,139.636) = 72.9921, f_{bb}(0.00504782,139.636) = 80.8349 P_{G,G}(0.9389) = {F=-5.26534,S=98.1995,I=0.0303636}, f_{G}(0.00537632,139.636) = 2269.4, f_{G}(0.00504782,139.636) = 2512.51 } -> 0.0293283 } Single_Process::CollinearCounterTerms(Q = 117.798 / 139.636) { \mu_F = 139.636 \mu_R = 53.0834 0.128528/(2\pi) * log(117.798/139.636) = -0.00695781 Beam 1: z = 0.558641, f_{u}(0.31457,139.636) = 0.982242 { P_{u,u}(0.558641) = {F=-2.07819,S=6.04195,I=1.00722}, f_{u}(0.563099,139.636) = 0.0988966, f_{u}(0.31457,139.636) = 0.982242 P_{G,u}(0.558641) = {F=0.253439,S=0,I=0}, f_{G}(0.563099,139.636) = 0.00706531, f_{G}(0.31457,139.636) = 0.230435 } -> 0.0184845 } * [ PDF(fla=ub, xa=0.000968497, ta=19498.3) * PDF(flb=u, xb=0.31457, tb=19498.3) -> 1634.59 / PDF(fla=ub, xa=0.000968497, ta=13876.4) * PDF(flb=u, xb=0.31457, tb=13876.4) -> 1624.15 ] = 1.00643 Single_Process::CollinearCounterTerms(Q = 139.636 / 117.798) { \mu_F = 139.636 \mu_R = 53.0834 0.128528/(2\pi) * log(139.636/117.798) = 0.00695781 Beam 0: z = 0.852511, f_{ub}(0.000968497,139.636) = 1664.14 { P_{ub,ub}(0.852511) = {F=-2.47001,S=18.0805,I=0.00258391}, f_{ub}(0.00113605,139.636) = 1330.44, f_{ub}(0.000968497,139.636) = 1664.14 P_{G,ub}(0.852511) = {F=0.374264,S=0,I=0}, f_{G}(0.00113605,139.636) = 25458.1, f_{G}(0.000968497,139.636) = 32337.7 } -> 0.0366611 } Single_Process::CollinearCounterTerms(Q = 139.636 / 117.798) { \mu_F = 139.636 \mu_R = 53.0834 0.128528/(2\pi) * log(139.636/117.798) = 0.00695781 Beam 1: z = 0.965392, f_{u}(0.31457,139.636) = 0.982242 { P_{u,u}(0.965392) = {F=-2.62052,S=77.0525,I=1.00722}, f_{u}(0.325847,139.636) = 0.898332, f_{u}(0.31457,139.636) = 0.982242 P_{G,u}(0.965392) = {F=0.466589,S=0,I=0}, f_{G}(0.325847,139.636) = 0.19823, f_{G}(0.31457,139.636) = 0.230435 } -> -0.0238107 } } type=BVIKP|METS, B=0.000209866, VI=-6.90652e-05, KP=0.000287762, RS=0 muR2=2817.85, muF2=19498.3, oqcd=2, oew=2, fl1=21, fl2=2, x1=1, x2=1, x1p=0.0266724, x2p=0.97788 wren=(1.64565e-05,2.91325e-21) wfac=(0,0.00176274,-4.80135e-05,0.00396029,0.00223801,-0.0021789,0,-8.91114e-06,-0,-0.000415803,1.00041e-05,-0.000936583,0.000349214,-0.000346974,-0,-1.40763e-06) Cluster sequence: pdfwgt=2497.34, flux=4.92004e-06, counter term=0.0606632, steps: (19498.3;21,0.00504782|2,0.31457) (13876.4;21,0.00504782|2,0.31457) (19498.3;-2,0.000968497|2,0.31457) Calculating DADS terms Single_Process::Differential(2_4__G__G__P__P__u__ub__QCD(S)_RS1_5_0, RS:0) { } DADS_0 = -0 DADS: wgt=-0, x1=0.00504782, x2=0.318498, fl1=21, fl2=21 } 2_3__G__u__P__P__u__QCD(BVI) ME = 1.03847042959957 , PS = 0.000215984564991741 -> 0.00022429358399385 p_lab[0]=Vec4D(20.1912858878542,0,0,20.1912858878545); p_lab[1]=Vec4D(1258.2805522151,0,0,-1258.2805522151); p_lab[2]=Vec4D(43.9602108368727,20.7705742278003,9.66197727340963,-37.5197225234656); p_lab[3]=Vec4D(1049.13980891077,-108.300053436409,-50.3786108704054,-1042.31829717883); p_lab[4]=Vec4D(185.371818355306,87.5294792086084,40.7166335969958,-158.251246624949); ========================================================== S event Jet_Finder::Trigger(): '2_3__G__u__P__P__u' Q_cut = 9.92952 { CS_Shower::JetVeto(mode = 0) { (0x7129fc0): 2 -> 3 { \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 = (), ncl = 0, flag = 0 (0) G (-20.1913,-0,-0,-20.1913) -3.36331e-06 (-1,-1) [0|0] (1) ub (-1258.28,-0,-0,1258.28) 0 (-1,-1) [0|0] (2) P (43.9602,20.7706,9.66198,-37.5197) 0 (-1,-1) [0|0] (3) P (1049.14,-108.3,-50.3786,-1042.32) -2.15792e-05 (-1,-1) [0|0] (4) u (185.372,87.5295,40.7166,-158.251) 3.8147e-06 (-1,-1) [0|0] } noem = (), nospec = () No kernel for G P <-> u (22) No kernel for G P <-> P (21) No kernel for G P <-> u (21) No kernel for G P <-> u (22) No kernel for G P <-> P (21) No kernel for G P <-> u (21) Q_{(0)(4),(1)} = 86.588 No kernel for G u <-> P (21) No kernel for G u <-> P (21) No kernel for u P <-> G (22) No kernel for u P <-> P (21) No kernel for u P <-> u (21) No kernel for u P <-> G (22) No kernel for u P <-> P (21) No kernel for u P <-> u (21) Q_{(1)(4),(0)} = 242.835 No kernel for u u <-> P (21) No kernel for u u <-> P (21) No kernel for P P <-> G (12) No kernel for P P <-> u (12) No kernel for P P <-> u (11) No kernel for P u <-> G (12) No kernel for P u <-> u (12) No kernel for P u <-> P (11) No kernel for P u <-> G (12) No kernel for P u <-> u (12) No kernel for P u <-> P (11) --- Jet veto --- } } -> 1 Single_Process::Cluster(): Found 1 amplitude(s) ... select 1st. CS_MCatNLO::GeneratePoint() { CS_MCatNLO::PrepareMCatNLO(): { (0xaadcb30): 2 -> 3 { \mu_r = 53.0834, \mu_f = 139.636, \mu_q = 139.636, \mu = 53.0834 k_T = 117.798, z = 0, phi = 0, kin = 0, K = 0 oew = 2, oqcd = 2, nlo = 0, new = (4), ncl = 0, flag = 0 (0) G (-20.1913,-0,-0,-20.1913) -3.36331e-06 (0,0) [0|6] (1) ub (-1258.28,-0,-0,1258.28) 0 (0,0) [0|6] (2) P (43.9602,20.7706,9.66198,-37.5197) 0 (0,0) [0|6] (3) P (1049.14,-108.3,-50.3786,-1042.32) -2.15792e-05 (0,0) [0|6] (4) u (185.372,87.5295,40.7166,-158.251) 3.8147e-06 (0,0) [0|6] } Singlet lists: Singlet parton list from CS_MCatNLO : 0xa825550 IS Parton 0xc62bea0 (1)[(0)]: G : (20.1913,0,0,20.1913) (0,0)[0,0] k_T start : 117.798 k_T test : 0 k_T veto : 9.92952(0) x_B : 0.00504782 IS Parton 0xc62db40 (1)[(1)]: u : (1258.28,0,0,-1258.28) (0,0)[0,0] k_T start : 117.798 k_T test : 0 k_T veto : 9.92952(0) x_B : 0.31457 FS Parton 0xa60a280 (1)[(2)]: P : (43.9602,20.7706,9.66198,-37.5197) (0,0)[0,0] k_T start : 117.798 k_T test : 0 k_T veto : 9.92952(0) x_B : 1 FS Parton 0xbfb5900 (1)[(3)]: P : (1049.14,-108.3,-50.3786,-1042.32) (0,0)[0,0] k_T start : 117.798 k_T test : 0 k_T veto : 9.92952(0) x_B : 1 FS Parton 0xc3f18b0 (1)[(4)]: u : (185.372,87.5295,40.7166,-158.251) (0,0)[0,0] k_T start : 117.798 k_T test : 0 k_T veto : 9.92952(0) x_B : 1 ------------------------------------------------------------------------- } filling partner list for G(0) ... u(1) u(4) -> 2 dipole(s) filling partner list for u(1) ... G(0) u(4) -> 2 dipole(s) filling partner list for P(2) ... -> 0 dipole(s) filling partner list for P(3) ... -> 0 dipole(s) filling partner list for u(4) ... G(0) u(1) -> 2 dipole(s) before mc@nlo step Singlet parton list from CS_MCatNLO : 0xa825550 IS Parton 0xc62bea0 (1)[(0)]: G : (20.1913,0,0,20.1913) (0,0)[0,0] k_T start : 117.798 k_T test : 0 k_T veto : 9.92952(0) x_B : 0.00504782 IS Parton 0xc62db40 (1)[(1)]: u : (1258.28,0,0,-1258.28) (0,0)[0,0] k_T start : 117.798 k_T test : 0 k_T veto : 9.92952(0) x_B : 0.31457 FS Parton 0xa60a280 (1)[(2)]: P : (43.9602,20.7706,9.66198,-37.5197) (0,0)[0,0] k_T start : 117.798 k_T test : 0 k_T veto : 9.92952(0) x_B : 1 FS Parton 0xbfb5900 (1)[(3)]: P : (1049.14,-108.3,-50.3786,-1042.32) (0,0)[0,0] k_T start : 117.798 k_T test : 0 k_T veto : 9.92952(0) x_B : 1 FS Parton 0xc3f18b0 (1)[(4)]: u : (185.372,87.5295,40.7166,-158.251) (0,0)[0,0] k_T start : 117.798 k_T test : 0 k_T veto : 9.92952(0) x_B : 1 ------------------------------------------------------------------------- Emission u -> u G at kt = 50.95, z = 0.819142, y = 0.115693 for IS Parton 0xc62db40 (1)[(1)]: u : (1258.28,0,0,-1258.28) (0,0)[0,0] k_T start : 117.798 k_T test : 50.95 k_T veto : 9.92952(0) x_B : 0.31457 IS Parton 0xc62bea0 (1)[(0)]: G : (20.1913,0,0,20.1913) (0,0)[0,0] k_T start : 117.798 k_T test : 20.2546 k_T veto : 9.92952(0) x_B : 0.00504782 Shower::MakeKinematics(0) { Singlet::ArrangeColours for u --> u & G CS_Gamma::TrialWeight() { CS_Gamma::CalculateWeight(0xaadd010) { (0xc5fc850): 2 -> 4 { \mu_r = 53.0834, \mu_f = 139.636, \mu_q = 139.636, \mu = 53.0834 k_T = 50.95, z = 0, phi = 0, kin = 0, K = 0 oew = 2, oqcd = 2, nlo = 0, new = (5), ncl = 0, flag = 0 (0) G (-20.1913,-0,-0,-20.1913) -3.36331e-06 (0,0) [0|6] (1) ub (-1536.1,-0,-0,1536.1) 0 (0,0) [0|6] (2) P (56.6521,23.5307,8.81335,-50.7748) 6.7435e-07 (0,0) [0|6] (3) P (1158.33,-93.4537,-54.9432,-1153.24) -1.52588e-05 (0,0) [0|6] (5) G (102.436,-29.2327,8.98781,-97.764) -1.90735e-06 (0,0) [0|6] (4) u (238.874,99.1557,37.1421,-214.125) 3.8147e-06 (0,0) [0|6] } Set fixed n_perp CS_Gamma::SingleWeight((1),(4)<->(0)) { B config -> (0xc61d5f0): 2 -> 3 { \mu_r = 53.0834, \mu_f = 50.95, \mu_q = 0, \mu = 0 k_T = 0, z = 0, phi = 0, kin = 0, K = 0 oew = 0, oqcd = 0, nlo = 0, new = (), ncl = 0, flag = 0 (0) G (-20.1913,-0,-0,-20.1913) -3.36331e-06 (0,0) [0|0] (1,4) ub (-1258.28,-0,-0,1258.28) 0 (0,0) [0|0] <-> (0) (2) P (43.9602,20.7706,9.66198,-37.5197) 6.7435e-07 (0,0) [0|0] (3) P (1049.14,-108.3,-50.3786,-1042.32) -2.15792e-05 (0,0) [0|0] (5) u (185.372,87.5295,40.7166,-158.251) 3.8147e-06 (0,0) [0|0] } -> CS{kt=50.95,z=0.819142,phi=2.05791,mode=3,kin=1} ( u -> u G ) zero matrix element } } } } } } } Sherpa: CS_Gamma::TrialWeight throws fatal error: Invalid amplitude for 2_4__G__u__P__P__G__u Exception_Handler::GenerateStackTrace(..): Generating stack trace { 0x2b01df14aebf in 'ATOOLS::Exception_Handler::GenerateStackTrace(std::ostream&, bool, std::string const&)' (Exception_Handler.C:318) from '/scratch/p_lhc/sherpa/trunk/r25839/lib/SHERPA-MC/libToolsOrg.so.0' 0x2b01df14c16e in 'ATOOLS::Exception_Handler::Terminate()' (basic_string.h:536) from '/scratch/p_lhc/sherpa/trunk/r25839/lib/SHERPA-MC/libToolsOrg.so.0' 0x2b01df14cd10 in 'ATOOLS::Terminate()' (Exception_Handler.C:135) from '/scratch/p_lhc/sherpa/trunk/r25839/lib/SHERPA-MC/libToolsOrg.so.0' 0x40a7f8 in 'main' (Main.C:34) from '/scratch/p_lhc/sherpa/trunk/r25839/bin/Sherpa' } Exception_Handler::Terminate(): Pre-crash status saved to '/scratch/p_lhc/projects/diphoton/debug/Status__Tue_Feb__9_09-30-09_2016'. Exception_Handler::ApproveTerminate(): Asking for termination ... ... approved. Exception_Handler::PrepareTerminate(): Preparing termination ... In Event_Handler::Finish : Summarizing the run may take some time. +-------------------------------------------------+ | | | Total XS is 0 pb +- ( 0 pb = -2.14748e+07 % ) | | | +-------------------------------------------------+ ... prepared. Exception_Handler::Exit: Exiting Sherpa with code (1)