Crash in EGamma with electron structure function

[Frank Siegert]

Examples/EGamma crashes if used with an electron structure function, e.g. by switching from PDF_SET=None to PDF_SET_1=None:

Starting the calculation. Lean back and enjoy ... .
CS_Cluster_Definitions::CoreScale(): Parton energies exceed collider energy.
  p_A = (14.7716,0,0,14.7716)
  p_B = (390.074,0,0,-390.074)
SF::CalculateWeight : x out of bounds 1 at 23048, xrange = 1e-06 ... 0.999999

Exception_Handler::SignalHandler: Signal (11) caught. 
Exception_Handler::GenerateStackTrace(..): Generating stack trace 
{
   3078685696 in   <unknown function>
  0xb6d1c155      in 'std::_Rb_tree<ATOOLS::Flavour, ATOOLS::Flavour, std::_Identity<ATOOLS::Flavour>, std::less<ATOOLS::Flavour>, std::allocator<ATOOLS::Flavour> >::find(ATOOLS::Flavour const&) const' (Splitting_Function_Base.C:375)                                                                                                                                                                               
  0xb665c709      in 'PDF::ISR_Handler::Weight(int, ATOOLS::Vec4<double>, ATOOLS::Vec4<double>, double, double, ATOOLS::Flavour, ATOOLS::Flavour)' (stl_set.h:609)
  0xb6449dd8      in 'PHASIC::Single_Process::BeamISRWeight(double const&, int) const' (Single_Process.C:63)
  0xb644a3e6      in 'PHASIC::Single_Process::Differential(std::vector<ATOOLS::Vec4<double>, std::allocator<ATOOLS::Vec4<double> > > const&)' (Single_Process.C:125)
  0xb64b1430      in 'PHASIC::Phase_Space_Handler::CalculateME()' (Phase_Space_Handler.C:855)
  0xb64b396c      in 'PHASIC::Phase_Space_Handler::Differential(PHASIC::Process_Integrator*, PHASIC::psm::code)' (Phase_Space_Handler.C:397)
  0xb64b420a      in 'PHASIC::Phase_Space_Handler::Differential()' (Phase_Space_Handler.C:248)
  0xb64bcf36      in 'PHASIC::Phase_Space_Integrator::Calculate(PHASIC::Phase_Space_Handler*, double, int)' (Phase_Space_Integrator.C:197)
  0xb64b11eb      in 'PHASIC::Phase_Space_Handler::Integrate()' (Phase_Space_Handler.C:181)
  0xb644908e      in 'PHASIC::Single_Process::CalculateTotalXSec(std::string const&, bool)' (Single_Process.C:185)
  0xb733516a      in 'SHERPA::Matrix_Element_Handler::CalculateTotalXSecs()' (Matrix_Element_Handler.C:122)
  0xb73ac52d      in 'SHERPA::Initialization_Handler::CalculateTheHardProcesses()' (Initialization_Handler.C:739)
  0xb78090c0      in 'SHERPA::Sherpa::InitializeTheRun(int, char**)' (Sherpa.C:169)
  0x804947c       in 'main' (Main.C:75)
}   Cannot continue.

Seems like the combination of beam spectrum and ISR doesn't work.

[Stefan Hoeche]

Closing ticket due to inactivity. Please reopen if still relevant.