# ====================== In 1.95:  whizard.cut1 ================ 

! changed from default in dbd
!>> 1.95  default_mass_cut = 4. (on pairs of m=0 colourless parts)
!>> 1.95   default_q_cut = 4.    (for photon & gluon exchange)
!>> 1.95 ! in addition 1.95 defaults:
!>> 1.95  default_energy_cut = 10. (photon and gluon emission)
!>> 1.95  default_jet_cut = 10.    (pairs of m=0 colored parts)

! attempt to formulate that:
real default_M_cut = 4 GeV
real default_jet_cut = 10 GeV
real default_E_cut = 10 GeV
real default_Q_cut = 4 GeV
alias uncoloured=lepton:neutrino:photon
alias masslessuncoloured=electrons:neutrino
alias masslessquark=quark
cuts =
  all M > default_M_cut [masslessuncoloured,masslessuncoloured]
  and
  all M > default_jet_cut [masslessquark,masslessquark]
  and
!!  let subevt @non_isr_photon = select if index > 2 [A] 
!  let subevt @non_isr_photon = select if Pt > 2 [A] 
!  in 
!  all E > default_E_cut [@non_isr_photon]
!  and
  all M < - default_Q_cut [incoming e1, e1]
  and
  all M < - default_Q_cut [incoming E1, E1]