Welcome to Sherpa, Daniel Reichelt. Initialization of framework underway. The local time is Thu Nov 30 10:04:26 2023. Run_Parameter::Init(): Setting memory limit to 503.276 GB. Random::SetSeed(): Seed set to 3032 ----------------------------------------------------------------------------- ----------- 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} = 85.0292 mb at PT_{min} = 3.39111 GeV } Profile_Function_Base::CalculateOMean(4.99237): Results are { k = 23.7135 <\tilde{O}> = 0.210525 } 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.40119e+08 pb +- ( 1.07296e+06 pb = 0.765747 % )  exp. eff: 5.86561e-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 ( 10s elapsed / 11d 23h 13m 9s left ) -> ETA: Tue Dec 12 09:17 XS = 9.99604 pb +- ( 9.99604 pb = 100 % ) Event 2 ( 10s elapsed / 6d 1m 29s left ) -> ETA: Wed Dec 06 10:06 XS = 2.49184e+06 pb +- ( 2.49183e+06 pb = 99.99 % ) Event 3 ( 10s elapsed / 4d 12m 2s left ) -> ETA: Mon Dec 04 10:16 XS = 4.98367e+06 pb +- ( 2.49183e+06 pb = 49.99 % ) Event 4 ( 10s elapsed / 3d 17m 19s left ) -> ETA: Sun Dec 03 10:21 XS = 7.47551e+06 pb +- ( 0 pb = 0 % ) Event 5 ( 10s elapsed / 2d 9h 56m 29s left ) -> ETA: Sat Dec 02 20:01 XS = 9.96735e+06 pb +- ( -nan pb = -2.14748e+07 % ) Event 6 ( 10s elapsed / 2d 33m 42s left ) -> ETA: Sat Dec 02 10:38 XS = 1.24592e+07 pb +- ( -nan pb = -2.14748e+07 % ) Event 7 ( 10s elapsed / 1d 17h 44m 35s left ) -> ETA: Sat Dec 02 03:49 XS = 1.4951e+07 pb +- ( -nan pb = -2.14748e+07 % ) Event 8 ( 10s elapsed / 1d 12h 46m 4s left ) -> ETA: Fri Dec 01 22:50 XS = 1.74429e+07 pb +- ( -nan pb = -2.14748e+07 % ) Event 9 ( 10s elapsed / 1d 8h 52m 2s left ) -> ETA: Fri Dec 01 18:56 XS = 1.99347e+07 pb +- ( -nan pb = -2.14748e+07 % ) Event 10 ( 10s elapsed / 1d 5h 41m 29s left ) -> ETA: Fri Dec 01 15:46 XS = 2.24265e+07 pb +- ( -nan pb = -2.14748e+07 % ) Event 20 ( 11s elapsed / 15h 38m 8s left ) -> ETA: Fri Dec 01 01:42 XS = 3.20444e+07 pb +- ( -nan pb = -2.14748e+07 % ) Event 30 ( 11s elapsed / 10h 47m 1s left ) -> ETA: Thu Nov 30 20:51 XS = 2.68506e+07 pb +- ( -nan pb = -2.14748e+07 % ) Event 40 ( 12s elapsed / 8h 23m 32s left ) -> ETA: Thu Nov 30 18:28 XS = 1.97348e+07 pb +- ( -nan pb = -2.14748e+07 % ) Event 50 ( 12s elapsed / 7h 4m 7s left ) -> ETA: Thu Nov 30 17:08 XS = 1.97768e+07 pb +- ( -nan pb = -2.14748e+07 % ) Event 60 ( 13s elapsed / 6h 5m 36s left ) -> ETA: Thu Nov 30 16:10 XS = 4.03366e+07 pb +- ( -nan pb = -2.14748e+07 % ) Event 70 ( 13s elapsed / 5h 23m 20s left ) -> ETA: Thu Nov 30 15:28 XS = 4.09649e+07 pb +- ( -nan pb = -2.14748e+07 % ) Event 80 ( 14s elapsed / 4h 54m 45s left ) -> ETA: Thu Nov 30 14:59 XS = 5.28421e+07 pb +- ( -nan pb = -2.14748e+07 % ) Event 90 ( 14s elapsed / 4h 30m 52s left ) -> ETA: Thu Nov 30 14:35 XS = 6.47193e+07 pb +- ( -nan pb = -2.14748e+07 % ) Event 100 ( 15s elapsed / 4h 11m 44s left ) -> ETA: Thu Nov 30 14:16 XS = 7.65966e+07 pb +- ( -nan pb = -2.14748e+07 % ) Event 200 ( 20s elapsed / 2h 46m 39s left ) -> ETA: Thu Nov 30 12:51 XS = 1.8389e+08 pb +- ( -nan pb = -2.14748e+07 % ) Event 300 ( 25s elapsed / 2h 19m 11s left ) -> ETA: Thu Nov 30 12:24 XS = 1.5018e+08 pb +- ( -nan pb = -2.14748e+07 % ) Event 400 ( 30s elapsed / 2h 5m 59s left ) -> ETA: Thu Nov 30 12:10 XS = 1.40739e+08 pb +- ( -nan pb = -2.14748e+07 % ) Event 500 ( 36s elapsed / 2h 9s left ) -> ETA: Thu Nov 30 12:05 XS = 1.23105e+08 pb +- ( -nan pb = -2.14748e+07 % ) Decay_Channel::GenerateKinematics(B+ --> Db vmu mu+ ) warning: d\Gamma(x)=2.4192e-13 > max(d\Gamma)=1.859e-13 Event 600 ( 42s elapsed / 1h 56m 34s left ) -> ETA: Thu Nov 30 12:01 XS = 1.14743e+08 pb +- ( -nan pb = -2.14748e+07 % ) Event 700 ( 49s elapsed / 1h 55m 52s left ) -> ETA: Thu Nov 30 12:01 XS = 1.10988e+08 pb +- ( -nan pb = -2.14748e+07 % ) Event 800 ( 55s elapsed / 1h 53m 48s left ) -> ETA: Thu Nov 30 11:59 XS = 1.08732e+08 pb +- ( -nan pb = -2.14748e+07 % ) Decay_Channel::GenerateKinematics(B+ --> D*(2007) D*(2007)b K+ pi ) warning: d\Gamma(x)=1.84251e-06 > max(d\Gamma)=9.231e-07 Event 900 ( 1m 1s elapsed / 1h 52m 8s left ) -> ETA: Thu Nov 30 11:57 XS = 1.12864e+08 pb +- ( -nan pb = -2.14748e+07 % ) Event 1000 ( 1m 6s elapsed / 1h 49m 52s left ) -> ETA: Thu Nov 30 11:55 XS = 1.19709e+08 pb +- ( -nan pb = -2.14748e+07 % ) Event 2000 ( 2m 6s elapsed / 1h 43m 6s left ) -> ETA: Thu Nov 30 11:49 XS = 9.81878e+07 pb +- ( -nan pb = -2.14748e+07 % ) Event 3000 ( 3m 8s elapsed / 1h 41m 20s left ) -> ETA: Thu Nov 30 11:48 XS = 9.15179e+07 pb +- ( -nan pb = -2.14748e+07 % ) Blob_List::FourMomentumConservation(): (0x215fee8) Four Momentum is not conserved. p_{in} = (13000,0,0,0) vs. p_{out} = (12992.5,0.230433,0.208341,-7.45407), diff = (-7.46054,0.230433,0.208341,-7.45407). Event 4000 ( 4m 6s elapsed / 1h 38m 25s left ) -> ETA: Thu Nov 30 11:47 XS = 1.40137e+08 pb +- ( -nan pb = -2.14748e+07 % ) Event 5000 ( 5m 6s elapsed / 1h 36m 57s left ) -> ETA: Thu Nov 30 11:46 XS = 1.57838e+08 pb +- ( -nan pb = -2.14748e+07 % ) Event 6000 ( 6m 7s elapsed / 1h 35m 59s left ) -> ETA: Thu Nov 30 11:46 XS = 1.36883e+08 pb +- ( -nan pb = -2.14748e+07 % ) Event 7000 ( 7m 9s elapsed / 1h 35m 7s left ) -> ETA: Thu Nov 30 11:46 XS = 1.37273e+08 pb +- ( -nan pb = -2.14748e+07 % ) Event 8000 ( 8m 10s elapsed / 1h 33m 58s left ) -> ETA: Thu Nov 30 11:46 XS = 1.28731e+08 pb +- ( -nan pb = -2.14748e+07 % ) Event 9000 ( 9m 12s elapsed / 1h 33m 2s left ) -> ETA: Thu Nov 30 11:46 XS = 1.26807e+08 pb +- ( -nan pb = -2.14748e+07 % ) Event 10000 ( 10m 10s elapsed / 1h 31m 35s left ) -> ETA: Thu Nov 30 11:46 XS = 1.23901e+08 pb +- ( -nan pb = -2.14748e+07 % ) Decay_Channel::GenerateKinematics(B+ --> D*(2007) D*(2010)+b K+ pi+ ) warning: d\Gamma(x)=2.63436e-06 > max(d\Gamma)=8.179e-07 Event 20000 ( 20m 1s elapsed / 1h 20m 5s left ) -> ETA: Thu Nov 30 11:44 XS = 1.29787e+08 pb +- ( -nan pb = -2.14748e+07 % ) Blob_List::FourMomentumConservation(): (0x215fee8) Four Momentum is not conserved. p_{in} = (13000,0,0,0) vs. p_{out} = (13001.2,0.155318,0.0529176,-1.14363), diff = (1.15537,0.155318,0.0529176,-1.14363). Event 30000 ( 30m 5s elapsed / 1h 10m 13s left ) -> ETA: Thu Nov 30 11:44 XS = 1.35173e+08 pb +- ( -nan pb = -2.14748e+07 % ) Decay_Channel::GenerateKinematics(B --> D*(2010)+ D+b K pi ) warning: d\Gamma(x)=3.24319e-06 > max(d\Gamma)=1.273e-06 Event 40000 ( 40m 3s elapsed / 1h 4s left ) -> ETA: Thu Nov 30 11:44 XS = 1.36627e+08 pb +- ( -nan pb = -2.14748e+07 % ) Decay_Channel::GenerateKinematics(B --> D*(2010)+ D*(2010)+b K pi ) warning: d\Gamma(x)=2.57572e-06 > max(d\Gamma)=7.854e-07 Event 50000 ( 50m 8s elapsed / 50m 8s left ) -> ETA: Thu Nov 30 11:44 XS = 1.33726e+08 pb +- ( -nan pb = -2.14748e+07 % ) Event 60000 ( 1h 6s elapsed / 40m 4s left ) -> ETA: Thu Nov 30 11:44 XS = 1.36337e+08 pb +- ( -nan pb = -2.14748e+07 % ) Decay_Channel::GenerateKinematics(B+ --> D*(2010)+ D*(2010)+b K pi+ ) warning: d\Gamma(x)=1.15586e-06 > max(d\Gamma)=7.299e-07 Blob_List::FourMomentumConservation(): (0x215fee8) Four Momentum is not conserved. p_{in} = (13000,0,0,0) vs. p_{out} = (12982.7,-4.08444,3.35397,-16.3813), diff = (-17.2935,-4.08444,3.35397,-16.3813). ...... Cluster_Part::TestDecay fails for ------------------------------------------------------------- Cluster [cluster, 168, 0] ((25.0108,-0.69657,-0.267591,24.6776), mass = 4.00017, y = 2.50233) : l : ss_1b (18.471,-0.959243,0.589417,18.3756) 1.5, kt_max = 0, pt = 1.12586, y = 2.97788 l : ss_1 (6.5398,0.262673,-0.857008,6.30202) 1.5, kt_max = 0, pt = 0.896359, y = 1.99457 Ahadic::Hadronize(): Momentum imbalance: (25.0108,-0.69657,-0.267591,24.6776) Decay_Channel::GenerateKinematics(B+ --> D D*(2010)+b K+ pi+ ) warning: d\Gamma(x)=2.97408e-06 > max(d\Gamma)=1.305e-06 Decay_Channel::GenerateKinematics(B --> D*(2007) D*(2010)+b K+ pi ) warning: d\Gamma(x)=1.71993e-06 > max(d\Gamma)=8.84e-07 Event 70000 ( 1h 10m 12s elapsed / 30m 5s left ) -> ETA: Thu Nov 30 11:44 XS = 1.41298e+08 pb +- ( -nan pb = -2.14748e+07 % ) Decay_Channel::GenerateKinematics(B+ --> D*(2007) D*(2010)+b K+ pi+ ) warning: d\Gamma(x)=3.74164e-06 > max(d\Gamma)=2.63436e-06 Decay_Channel::GenerateKinematics(B+ --> D*(2007)b D*(2010)+ K pi ) warning: d\Gamma(x)=2.77362e-06 > max(d\Gamma)=8.264e-07 Event 80000 ( 1h 19m 57s elapsed / 19m 59s left ) -> ETA: Thu Nov 30 11:44 Event_Handler::MPISync() {  Memory usage increased by 18 MB, now 179 MB.  This might indicate a memory leak! Please monitor this process closely. } XS = 1.37608e+08 pb +- ( -nan pb = -2.14748e+07 % ) Event 90000 ( 1h 29m 41s elapsed / 9m 57s left ) -> ETA: Thu Nov 30 11:44 XS = 1.35352e+08 pb +- ( -nan pb = -2.14748e+07 % ) Decay_Channel::GenerateKinematics(B --> D*(2007)b D+ K pi+b ) warning: d\Gamma(x)=3.98676e-06 > max(d\Gamma)=1.247e-06 Decay_Channel::GenerateKinematics(B+ --> D*(2007)b D*(2010)+ K pi ) warning: d\Gamma(x)=3.72773e-06 > max(d\Gamma)=2.77362e-06 Event 100000 ( 5965 s total ) = 1.44834e+06 evts/day In Event_Handler::Finish : Summarizing the run may take some time. Rivet_Interface::Finish(Analysis_CSS/3032){ } +-------------------------------------------------------------+ | | | Total XS is 1.3358e+08 pb +- ( -nan pb = -2.14748e+07 % ) | | | +-------------------------------------------------------------+ Return_Value::PrintStatistics(): Statistics { Generated events: 100000 New events { From "Beam_Remnants": 3 (100123) -> 0 % From "Jet_Evolution:CSS": 3 (798509) -> 0 % } Retried events { From "Beam_Remnants": 62 (100123) -> 0 % From "Hadronization:Ahadic": 58 (134271) -> 0 % From "Jet_Evolution:CSS": 7075 (798509) -> 0.8 % } Retried phases { From "Hadron_Decay_Handler::RejectExclusiveChannelsFromFragmentation": 16301 (0) -> 16301. } Retried methods { From "Decay_Channel::GenerateKinematics": 37 (42731048) -> 0 % } } Time: 1h 39m 28s on Thu Nov 30 11:43:55 2023 (User: 1h 39m 15s, System: 1s, Children User: 0s, Children System: 0s)