anchor.diff:        isExact = .true.
   
       else
  +      ! We just let Pythia generate an event assuming both colliding particles
  +      ! are on-momentum, and collide head on.
   
         call pythia_getEvent(evStat, evType, t, theta, dEE, dPP)
   
  
positive.code: 
      call ranecu(rndVals, 2, 0)
  
      pVec(1) = ((rndVals(1)/sqrt(betaX))*c1m3 - (alphaX/betaX)*(xv1(j)-orbX) + orbXP*oidpsv(j)) * ejfv(j)*c1m3
      pVec(2) = ((rndVals(2)/sqrt(betaY))*c1m3 - (alphaY/betaY)*(xv2(j)-orbY) + orbYP*oidpsv(j)) * ejfv(j)*c1m3
      pVec(3) = -sqrt(ejfv(j)**2 - pVec(1)**2 - pVec(2)**2)
  
      call recuut(scatter_seed1,scatter_seed2)
  
negative.code: #ifdef PYTHIA
  10  continue
      if(pythia_useRealBeam) then
        ! We sample the actual SixTrack tracked beam 1, and sample a colliding
        ! particle from beam 2
  
        ! Compute the particle three-momenta in GeV, using division rather
        ! than multiplication for scaling to avoid additional numerical error
        pIn(1) = yv1(j)*ejfv(j)/c1e6
        pIn(2) = yv2(j)*ejfv(j)/c1e6
        pIn(3) = sqrt(ejfv(j)**2 - pIn(1)**2 - pIn(2)**2)/c1e3
  
        call scatter_generateParticle(idPro, iElem, j, pGen)
  
  