There seems to be an inversion in the definition of the sigmaz for the strong beam

[giadarol]

  if(mylhcbeam==1)
  {
  beamx = beam%lhcb2->ex;
  beamy = beam%lhcb2->ey;
  sigz  = beam%lhcb1->sigt;
  sige  = beam%lhcb1->sige;
  };  
  if(mylhcbeam==2)
  {
  beamx = beam%lhcb1->ex;
  beamy = beam%lhcb1->ey;
  sigz  = beam%lhcb2->sigt;
  sige  = beam%lhcb2->sige;
  };