Ticket #250: casc_dec.sin

model = MSSM

process dec_su_q        = su1 -> u, neu2
process dec_su_qc       = SU1 -> U, neu1
process dec_neu_sl      = neu2 -> se12, E1
process dec_neu_sl2     = neu2 -> SE12, e1
process dec_sl_l        = se12 -> e1, neu1
process dec_sl_l2       = SE12 -> E1, neu1
 

process susybg = u,U -> SU1, su1

compile

read_slha("sps1a.in")

integrate (dec_su_q,dec_su_qc,dec_neu_sl,dec_sl_l,dec_neu_sl2,dec_sl_l2) { iterations = 1:1000 }


sqrts = 14000
beams = p, p -> lhapdf
integrate (susybg) { iterations = 5:10000, 2:10000 }

unstable su1 (dec_su_q)
unstable SU1 (dec_su_qc)
!unstable neu2 (dec_neu_sl)  !either decay into positron
unstable neu2 (dec_neu_sl2)  !or electron, which then should show spin correlation effects in the invariant mass of e+/e- and u
unstable se12 (dec_sl_l)
unstable SE12 (dec_sl_l2)

n_events = 10000

histogram inv_mass1 (0,1000,10)
histogram inv_mass2 (0,1000,10)

analysis = 
  record inv_mass1 (eval M / 1 GeV [combine[u,e1]])
  and
  record inv_mass2 (eval M / 1 GeV [combine[u,E1]])

simulate (susybg)

write_analysis
$analysis_filename = "casc_dec"
write_analysis