Ticket #460: ticket460.sin~

model = MSSM

alias j = g:u:U:d:D:s:S:c:C
!alias quark = d
alias antiq = D

process base = antiq,antiq => su1:SU1:su2:SU2,su1:SU1:su2:SU2
process base2 = antiq,antiq => sd1:SD1:sd2:SD2,sd1:SD1:sd2:SD2

process dec_su1 = su1 => j,e3:E3,SE31:se31
process dec_su2 = SU1 => j,e3:E3,SE31:se31
process dec_su3 = su2 => j,e3:E3,SE31:se31
process dec_su4 = SU2 => j,e3:E3,SE31:se31

process dec_sd1 = sd1 => j,e3:E3,SE31:se31
process dec_sd2 = SD1 => j,e3:E3,SE31:se31
process dec_sd3 = sd2 => j,e3:E3,SE31:se31
process dec_sd4 = SD2 => j,e3:E3,SE31:se31

compile
?slha_read_decays = true
read_slha("staus.slha")

sqrts = 7000
beams = p, p => lhapdf

mc=0
ms=0

cuts =  all Pt > 20 [j]
        and all abs(Eta) < 2.5 [j]
        and all Dist > 0.5 [j,j]

integrate (base,base2) !{ iterations = 10:100000, 5:200000 }

n_events = 100

unstable su1 (dec_su1)
unstable SU1 (dec_su2)
unstable su2 (dec_su3)
unstable SU2 (dec_su4)

unstable sd1 (dec_sd1)
unstable SD1 (dec_sd2)
unstable sd2 (dec_sd3)
unstable SD2 (dec_sd4)


simulate (base, base2) {
  $sample = "janstau_compl"
  sample_format = hepmc, lhef
  checkpoint = 10       
}