anchor.diff:      rhoOut  = (nBeam*normFac) * exp_mb(-half*(xProj/sigmaX)**2 -half*(yProj/sigmaY)**2)
   
     case default
  -    write(lerr,"(a,i0,a)") "SCATTER> ERROR Type ",scatter_proList(idPro)%proType," for profile '"//&
  -      trim(scatter_proList(idPro)%proName)//"' not understood"
  +
  +    write(lerr,"(a,i0,a)") "SCATTER> ERROR Unknown profile type ",scatter_proList(idPro)%proType," for profile '"//&
  +      trim(scatter_proList(idPro)%proName)//"'"
       call prror
  +
     end select
   
   end subroutine scatter_getDensity
  
positive.code: #endif
  
    case default
      write(lerr,"(a,i0,a)") "SCATTER> ERROR Generator type ",scatter_genList(idGen)%genType," not understood"
      call prror
  
    end select
  
negative.code:       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)
  
  