Ticket #195: bug22.sin

alias leptons = e1:E1:e2:E2

!process ee1 = e2,E2 -> leptons, leptons, leptons, leptons
!process ee2 = e2,E2 -> E1, E1, e1,e1
process em2 = e1,E1 -> E2, E1, e2,e1
!process emm2 = e2,E2 -> E2, E2, e2,e2
!process eu = e2, E2 -> e1, E1, u, U


compile

sqrts = 500
beams = e1, E1
iterations = 20:5000


! no cuts
seed=0
echo (" NO CUTS ")
integrate(em2)


cuts = all Pt > 10 GeV [e1:E1:e2:E2:u:U]
       and all Dist > 1 [e1:E1:e2:E2:u:U,e1:E1:e2:E2:u:U]

seed=0
echo (" NORMAL CUTS")
integrate(em2)


cuts = all Pt > 10 GeV [e1:E1:e2:E2:u:U]
       and all Dist > 1 [e1:E1:e2:E2:u:U,e1:E1:e2:E2:u:U]
       and all E > 1 GeV [e1:E1:e2:E2:u:U]
       and all 5 degree < Theta < 175 degree [e1:E1:e2:E2:u:U]
       and all 5 degree < Theta < 175 degree [e1:E1:e2:E2:u:U,e1:E1:e2:E2:u:U]
       and all M2 > (10 GeV)^2 [combine [e1:E1:e2:E2:u:U, e1:E1:e2:E2:u:U]]
       and all M2 < -(10 GeV)^2 [combine [incoming leptons, leptons]]
       and all Pt > 10 GeV [combine[e2,E2]]
       and all Pt > 10 GeV [combine[e1,E1]]
       and all E > 10 GeV [combine[e1,E1]]
       and all E > 10 GeV [combine[e2,E2]]
seed=0
echo (" ALL CUTS ")
integrate(em2)

!  ! Automatically generated set of cuts
!  ! Process eeeemm:
!  !    e a-e ->   e a-e  mu a-mu
!  !   32  16 ->   1   2   4    8
!  process eeeemm
!  cut M of    3      within  1.00000E+01  1.00000E+99
!  cut M of   33      within -1.00000E+99 -1.00000E+01
!  cut M of   18      within -1.00000E+99 -1.00000E+01

cuts = all M2 > (10 GeV)^2 [combine[e1,E1]]
       and all M2 < -(10 GeV)^2 [combine[incoming e1, e1]]
       and all M2 < -(10 GeV)^2 [combine[incoming E1, E1]]
seed = 0
echo(" WHIZARD 1.93 standard cuts ")
integrate(em2)