Welcome to Sherpa, Daniel Reichelt. Initialization of framework underway. The local time is Thu Nov 30 11:32:10 2023. Run_Parameter::Init(): Setting memory limit to 503.276 GB. Random::SetSeed(): Seed set to 3028 ----------------------------------------------------------------------------- ----------- 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.15 (Cho Oyu) Authors: Enrico Bothmann, 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. ----------------------------------------------------------------------------- Git branch rel-2-2-15, revision 9f988e9b9e17ef2f1eee27b92ac20c73c969aec6. ===> DIRE/Shower has local modifications <=== Beam_Spectra_Handler : type = Monochromatic*Monochromatic for P+ ((6500,0,0,6500)) and P+ ((6500,0,0,-6500)) WARNING: M_TOP=0 replacing with SHERPA top mass: M_TOP=173.21 PDF set 'NNPDF30NNLO' loaded for beam 1 (P+). WARNING: M_TOP=0 replacing with SHERPA top mass: M_TOP=173.21 PDF set 'NNPDF30NNLO' loaded for beam 2 (P+). Initialized the ISR: (SF)*(SF) WARNING: M_TOP=0 replacing with SHERPA top mass: M_TOP=173.21 WARNING: M_TOP=0 replacing with SHERPA top mass: M_TOP=173.21 Standard_Model::FixEWParameters(EW_SCHEME=1){ 288:1./aqed->Default()=128.802 289:csin2thetaW=(0.222928,-0.00110708) 290:cvev=(243.034,-3.75493) } 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} 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} Library_Loader::LoadLibrary(): Failed to load library 'libSherpaHepMCOutput.so'. 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 color setter mode 0 CS_Shower::CS_Shower(): Set respect Q2 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 pilot run mode 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 /mt/home/dreichelt/hepsoftware/src/OpenLoops #################################################################### # ___ ___ ____ _ _ __ __ ___ __ ___ # # / \ |__| |__ |\ | | / \ / \ |__| /__ __/ # # \___/ | |___ | \| |__ \__/ \__/ | __/ /__ # # You are using OpenLoops 2 to evaluate loop amplitudes # # F. Buccioni, J.-N. Lang, J. Lindert, P. Maierhoefer, # # S. Pozzorini, M. Zoller, H. Zhang # # Eur.Phys.J. C79 (2019) no.10, 866 # # arXiv:1907.13071 # #################################################################### Matrix_Element_Handler::BuildProcesses(): Looking for processes ................................ done ( 24 MB, 0s / 0s ). Matrix_Element_Handler::InitializeProcesses(): Performing tests done ( 25 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} = 39.8363 mb at PT_{min} = 4.65296 GeV } Profile_Function_Base::CalculateOMean(2.33893): Results are { k = 7.24375 <\tilde{O}> = 0.32289 } Initialized the Multiple_Interactions_Handler (MI_Handler). Initialized the Soft_Photon_Handler. Variations::InitialiseParametersVector(){ Number of variations: 7 Named variations: MUR0.5_MUF0.5_PDF261000_PSMUR0.5_PSMUF0.5 MUR0.5_MUF1_PDF261000_PSMUR0.5_PSMUF1 MUR1_MUF0.5_PDF261000_PSMUR1_PSMUF0.5 MUR1_MUF1_PDF261000 MUR1_MUF2_PDF261000_PSMUR1_PSMUF2 MUR2_MUF1_PDF261000_PSMUR2_PSMUF1 MUR2_MUF2_PDF261000_PSMUR2_PSMUF2 } Read in channels from directory : Results/Amegic/MC_2_2__j__j__j__j Process_Group::CalculateTotalXSec(): Calculate xs for '2_2__j__j__j__j' (Amegic) 2_2__j__j__j__j : 1.40117e+08 pb +- ( 1.07294e+06 pb = 0.765748 % )  exp. eff: 5.86553e-07 % reduce max for 2_2__j__j__j__j to 0.00494632 ( eps = 0.001 ) ---------------------------------------------------------- -- 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 Unknown : Analysis : Rivet --------------------------------------------------------- Event 1 ( 11s elapsed / 13d 11h 36m 28s left ) -> ETA: Wed Dec 13 23:08 XS = 12783.8 pb +- ( 12783.8 pb = 100 % ) Event 2 ( 11s elapsed / 6d 18h 29m 48s left ) -> ETA: Thu Dec 07 06:02 XS = 25567.6 pb +- ( 25567.6 pb = 100 % ) Event 3 ( 11s elapsed / 4d 12h 47m 34s left ) -> ETA: Tue Dec 05 00:19 XS = 38351.4 pb +- ( 38351.4 pb = 100 % ) Event 4 ( 11s elapsed / 3d 9h 56m 28s left ) -> ETA: Sun Dec 03 21:28 XS = 51135.3 pb +- ( 51135.3 pb = 100 % ) Event 5 ( 11s elapsed / 2d 17h 49m 48s left ) -> ETA: Sun Dec 03 05:22 XS = 63919.1 pb +- ( 63919.1 pb = 100 % ) Event 6 ( 11s elapsed / 2d 6h 59m 48s left ) -> ETA: Sat Dec 02 18:32 XS = 76702.9 pb +- ( 76702.9 pb = 100 % ) Event 7 ( 11s elapsed / 1d 23h 20m 16s left ) -> ETA: Sat Dec 02 10:52 XS = 89486.7 pb +- ( 89486.7 pb = 100 % ) Event 8 ( 11s elapsed / 1d 17h 31m 28s left ) -> ETA: Sat Dec 02 05:03 XS = 102271 pb +- ( 102271 pb = 100 % ) Event 9 ( 12s elapsed / 1d 13h 3m 52s left ) -> ETA: Sat Dec 02 00:36 XS = 115054 pb +- ( 115054 pb = 100 % ) Event 10 ( 12s elapsed / 1d 9h 29m 47s left ) -> ETA: Fri Dec 01 21:02 XS = 127838 pb +- ( 127838 pb = 100 % ) Event 20 ( 12s elapsed / 17h 19m 47s left ) -> ETA: Fri Dec 01 04:52 XS = 2.986e+06 pb +- ( 2.9048e+06 pb = 97.28 % ) Event 30 ( 12s elapsed / 11h 56m 27s left ) -> ETA: Thu Nov 30 23:28 XS = 3.96441e+06 pb +- ( 2.44634e+06 pb = 61.7 % ) Event 40 ( 13s elapsed / 9h 9m 21s left ) -> ETA: Thu Nov 30 20:41 XS = 4.47373e+06 pb +- ( 1.83072e+06 pb = 40.92 % ) Event 50 ( 13s elapsed / 7h 31m 46s left ) -> ETA: Thu Nov 30 19:04 XS = 4.69448e+07 pb +- ( 5.39472e+06 pb = 11.49 % ) Event 60 ( 13s elapsed / 6h 26m 42s left ) -> ETA: Thu Nov 30 17:59 XS = 1.08734e+08 pb +- ( -nan pb = -2.14748e+07 % ) Event 70 ( 14s elapsed / 5h 39m 2s left ) -> ETA: Thu Nov 30 17:11 XS = 1.47943e+08 pb +- ( -nan pb = -2.14748e+07 % ) Event 80 ( 14s elapsed / 5h 4m 57s left ) -> ETA: Thu Nov 30 16:37 XS = 1.48272e+08 pb +- ( -nan pb = -2.14748e+07 % ) Event 90 ( 15s elapsed / 4h 38m 49s left ) -> ETA: Thu Nov 30 16:11 XS = 1.34941e+08 pb +- ( -nan pb = -2.14748e+07 % ) Event 100 ( 15s elapsed / 4h 18m 14s left ) -> ETA: Thu Nov 30 15:50 XS = 1.35025e+08 pb +- ( -nan pb = -2.14748e+07 % ) Event 200 ( 19s elapsed / 2h 39m 50s left ) -> ETA: Thu Nov 30 14:12 XS = 9.45849e+07 pb +- ( -nan pb = -2.14748e+07 % ) Event 300 ( 22s elapsed / 2h 3m 47s left ) -> ETA: Thu Nov 30 13:36 XS = 2.86011e+08 pb +- ( -nan pb = -2.14748e+07 % ) Event 400 ( 25s elapsed / 1h 47m 19s left ) -> ETA: Thu Nov 30 13:19 XS = 2.77827e+08 pb +- ( -nan pb = -2.14748e+07 % ) Event 500 ( 30s elapsed / 1h 41m 11s left ) -> ETA: Thu Nov 30 13:13 XS = 2.641e+08 pb +- ( -nan pb = -2.14748e+07 % ) Event 600 ( 34s elapsed / 1h 35m 38s left ) -> ETA: Thu Nov 30 13:08 XS = 1.9788e+08 pb +- ( -nan pb = -2.14748e+07 % ) Event 700 ( 38s elapsed / 1h 30m 23s left ) -> ETA: Thu Nov 30 13:03 XS = 1.59339e+08 pb +- ( -nan pb = -2.14748e+07 % ) Event 800 ( 42s elapsed / 1h 26m 54s left ) -> ETA: Thu Nov 30 12:59 XS = 1.60371e+08 pb +- ( -nan pb = -2.14748e+07 % ) Event 900 ( 46s elapsed / 1h 24m 40s left ) -> ETA: Thu Nov 30 12:57 XS = 1.55776e+08 pb +- ( -nan pb = -2.14748e+07 % ) Event 1000 ( 49s elapsed / 1h 21m 52s left ) -> ETA: Thu Nov 30 12:54 XS = 1.5657e+08 pb +- ( -nan pb = -2.14748e+07 % ) Decay_Channel::GenerateKinematics(B --> D*(2007)b D*(2010)+ K pi+b ) warning: d\Gamma(x)=1.11779e-06 > max(d\Gamma)=7.635e-07 Event 2000 ( 1m 27s elapsed / 1h 11m 17s left ) -> ETA: Thu Nov 30 12:44 XS = 1.28856e+08 pb +- ( -nan pb = -2.14748e+07 % ) Event 3000 ( 2m 3s elapsed / 1h 6m 23s left ) -> ETA: Thu Nov 30 12:40 XS = 1.37529e+08 pb +- ( -nan pb = -2.14748e+07 % ) Event 4000 ( 2m 40s elapsed / 1h 4m 6s left ) -> ETA: Thu Nov 30 12:38 XS = 1.39898e+08 pb +- ( -nan pb = -2.14748e+07 % ) Decay_Channel::GenerateKinematics(B --> D*(2010)+b pi+ pi pi ) warning: d\Gamma(x)=1.86626e-13 > max(d\Gamma)=1.381e-13 Event 5000 ( 3m 20s elapsed / 1h 3m 26s left ) -> ETA: Thu Nov 30 12:38 XS = 1.31893e+08 pb +- ( -nan pb = -2.14748e+07 % ) Event 6000 ( 3m 59s elapsed / 1h 2m 34s left ) -> ETA: Thu Nov 30 12:38 XS = 1.27244e+08 pb +- ( -nan pb = -2.14748e+07 % ) Event 7000 ( 4m 36s elapsed / 1h 1m 8s left ) -> ETA: Thu Nov 30 12:37 XS = 1.27771e+08 pb +- ( -nan pb = -2.14748e+07 % ) Event 8000 ( 5m 11s elapsed / 59m 42s left ) -> ETA: Thu Nov 30 12:37 XS = 1.47488e+08 pb +- ( -nan pb = -2.14748e+07 % ) Event 9000 ( 5m 48s elapsed / 58m 46s left ) -> ETA: Thu Nov 30 12:36 XS = 1.45814e+08 pb +- ( -nan pb = -2.14748e+07 % ) Event 10000 ( 6m 29s elapsed / 58m 23s left ) -> ETA: Thu Nov 30 12:37 XS = 1.43496e+08 pb +- ( -nan pb = -2.14748e+07 % ) Decay_Channel::GenerateKinematics(B+ --> D*(2010)+ D*(2010)+b K pi+ ) warning: d\Gamma(x)=2.00255e-06 > max(d\Gamma)=7.299e-07 Event 20000 ( 12m 47s elapsed / 51m 9s left ) -> ETA: Thu Nov 30 12:36 XS = 1.36394e+08 pb +- ( -nan pb = -2.14748e+07 % ) Event 30000 ( 18m 59s elapsed / 44m 18s left ) -> ETA: Thu Nov 30 12:35 XS = 1.35488e+08 pb +- ( -nan pb = -2.14748e+07 % ) Decay_Channel::GenerateKinematics(B+ --> D*(2007) D*(2010)+b K+ pi+ ) warning: d\Gamma(x)=2.28281e-06 > max(d\Gamma)=8.179e-07 Decay_Channel::GenerateKinematics(B --> D*(2007) D+b K+ pi ) warning: d\Gamma(x)=4.41019e-06 > max(d\Gamma)=1.354e-06 Decay_Channel::GenerateKinematics(B+ --> D*(2010)+ D*(2010)+b K pi+ ) warning: d\Gamma(x)=5.53693e-06 > max(d\Gamma)=2.00255e-06 Event 40000 ( 25m 7s elapsed / 37m 40s left ) -> ETA: Thu Nov 30 12:34 XS = 1.45563e+08 pb +- ( -nan pb = -2.14748e+07 % ) Decay_Channel::GenerateKinematics(B+ --> K+ omega(782) pi ): Rejected decay kinematics 10000 times. This indicates a wrong maximum. Will accept kinematics. Decay_Channel::GenerateKinematics(B --> D*(2007)b D*(2010)+ K pi+b ) warning: d\Gamma(x)=3.80769e-06 > max(d\Gamma)=1.11779e-06 Decay_Channel::GenerateKinematics(B --> Db D*(2010)+ K pi+b ) warning: d\Gamma(x)=3.53362e-06 > max(d\Gamma)=1.257e-06 Decay_Channel::GenerateKinematics(B+ --> D+ D+b K pi+ ) warning: d\Gamma(x)=1.89477e-06 > max(d\Gamma)=1.145e-06 Decay_Channel::GenerateKinematics(B+ --> D*(2007)b D*(2010)+ K pi ) warning: d\Gamma(x)=1.48964e-06 > max(d\Gamma)=8.264e-07 Event 50000 ( 31m 16s elapsed / 31m 16s left ) -> ETA: Thu Nov 30 12:34 XS = 1.47592e+08 pb +- ( -nan pb = -2.14748e+07 % ) Decay_Channel::GenerateKinematics(B+ --> D*(2007) D*(2007)b K+ pi ) warning: d\Gamma(x)=1.56536e-06 > max(d\Gamma)=9.231e-07 Decay_Channel::GenerateKinematics(B+ --> D*(2007) D*(2007)b K+ pi ) warning: d\Gamma(x)=4.53024e-06 > max(d\Gamma)=1.56536e-06 Decay_Channel::GenerateKinematics(B+ --> D*(2007) D+b K+ pi+ ) warning: d\Gamma(x)=2.45333e-06 > max(d\Gamma)=1.294e-06 Decay_Channel::GenerateKinematics(B --> D*(2007)b D+ K pi+b ) warning: d\Gamma(x)=2.91396e-06 > max(d\Gamma)=1.247e-06 Event 60000 ( 37m 18s elapsed / 24m 52s left ) -> ETA: Thu Nov 30 12:34 XS = 1.43805e+08 pb +- ( -nan pb = -2.14748e+07 % ) Decay_Channel::GenerateKinematics(B --> D*(2007) Db K+ pi+b ) warning: d\Gamma(x)=2.10744e-06 > max(d\Gamma)=1.349e-06 Event 70000 ( 43m 16s elapsed / 18m 32s left ) -> ETA: Thu Nov 30 12:34 XS = 1.44153e+08 pb +- ( -nan pb = -2.14748e+07 % ) Event 80000 ( 49m 2s elapsed / 12m 15s left ) -> ETA: Thu Nov 30 12:33 XS = 1.44561e+08 pb +- ( -nan pb = -2.14748e+07 % ) Decay_Channel::GenerateKinematics(B+ --> D*(2007)b D*(2010)+ K pi ) warning: d\Gamma(x)=2.37468e-06 > max(d\Gamma)=1.48964e-06 Decay_Channel::GenerateKinematics(B+ --> Db D*(2010)+ K pi ) warning: d\Gamma(x)=1.90381e-06 > max(d\Gamma)=1.32e-06 Event 90000 ( 54m 56s elapsed / 6m 6s left ) -> ETA: Thu Nov 30 12:33 XS = 1.42354e+08 pb +- ( -nan pb = -2.14748e+07 % ) ...... Cluster_Part::TestDecay fails for ------------------------------------------------------------- Cluster [cluster, 78, 0] ((5.55932,-0.861997,0.447628,4.08175), mass = 3.64718, y = 0.937816) : l : ud_0b (1.33651,-0.524302,0.68626,0.156047) 1.008, kt_max = 0.38765, pt = 0.863624, y = 0.117292 l : ss_1 (4.22281,-0.337695,-0.238632,3.9257) 1.5, kt_max = 0.447062, pt = 0.413501, y = 1.65575 Ahadic::Hadronize(): Momentum imbalance: (5.55932,-0.861997,0.447628,4.08175) Decay_Channel::GenerateKinematics(B+ --> D*(2010)+b D+ K pi+ ) warning: d\Gamma(x)=3.58316e-06 > max(d\Gamma)=1.29581e-06 Decay_Channel::GenerateKinematics(B --> D D+b K+ pi ) warning: d\Gamma(x)=1.73552e-06 > max(d\Gamma)=1.263e-06 Decay_Channel::GenerateKinematics(B --> D*(2007) D*(2007)b K+ pi+b ) warning: d\Gamma(x)=4.62386e-06 > max(d\Gamma)=8.556e-07 Event 100000 ( 3641 s total ) = 2.37267e+06 evts/day Event_Handler::MPISync() {  Memory usage increased by 19 MB, now 175 MB.  This might indicate a memory leak! Please monitor this process closely. } In Event_Handler::Finish : Summarizing the run may take some time. Rivet_Interface::Finish(Analysis_CSS/3028){ } +--------------------------------------------------------------+ | | | Total XS is 1.43293e+08 pb +- ( -nan pb = -2.14748e+07 % ) | | | +--------------------------------------------------------------+ Return_Value::PrintStatistics(): Statistics { Generated events: 100000 Retried events { From "Beam_Remnants": 43 (100052) -> 0 % From "Hadronization:Ahadic": 9 (127593) -> 0 % From "Jet_Evolution:CSS": 2383 (377288) -> 0.6 % } Retried phases { From "Hadron_Decay_Handler::RejectExclusiveChannelsFromFragmentation": 13093 (0) -> 13093. } Retried methods { From "Decay_Channel::GenerateKinematics": 47 (28009099) -> 0 % } } Time: 1h 44s on Thu Nov 30 12:32:54 2023 (User: 1h 34s, System: 1s, Children User: 0s, Children System: 0s)