Ticket #250: casc_res.sin

model = MSSM

!process susyres = u,U -> SU1, u, e1, E1, neu1 { $restrictions = "4+5+6+7~su1&&5+6+7~neu2 && 6+7~se2+"}
process susyres = u,U -> SU1, u, e1, E1, neu1 { $restrictions = "4+5+6+7~su1&&5+6+7~neu2 && 5+7~se2-"}

compile

read_slha("sps1a.in")

wneu2 = 0.006
wse2 = 0.2
wsu1 = 1.8

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

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 (susyres)

write_analysis
$analysis_filename = "casc_rminus"
write_analysis