whizard is hosted by Hepforge, IPPP Durham

Changes between Version 1 and Version 2 of ExampleCascadeDecays


Ignore:
Timestamp:
Apr 12, 2010, 9:52:36 AM (15 years ago)
Author:
ohl
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • ExampleCascadeDecays

    v1 v2  
    11{{{
    22model = MSSM
    3 
     3}}}
     4{{{
    45process dec_su_q        = su1 => u, neu2
    56process dec_neu_sl2     = neu2 => SE12, e1
    6 
     7}}}
     8{{{
    79process susybg = u,U => SU1, su1
    810process full = u, U => SU1, u, e1, SE12
    9 
     11}}}
     12{{{
    1013compile
    11 
     14}}}
     15{{{
    1216?slha_read_decays = true
    1317read_slha("sps1ap_decays.slha")
    14 
     18}}}
     19{{{
    1520integrate (dec_su_q, dec_neu_sl2) { iterations = 1:1000 }
    16 
     21}}}
     22{{{
    1723sqrts = 14000
    1824beams = p, p => lhapdf
    19 
     25}}}
     26{{{
    2027integrate (susybg) { iterations = 5:10000, 2:10000 }
    2128integrate (full)
    22 
     29}}}
     30{{{
    2331n_events = 10000
    24 
     32}}}
     33{{{
    2534$title = "Full process"
    2635$description =
     
    2837$xlabel = "$M_{\rm inv}(ue^-)$"
    2938histogram inv_mass1_full (0,600,20)
    30 
     39}}}
     40{{{
    3141$title = "Factorized process with complete spin correlations"
    3242$description =
     
    3444  (\tilde u_1 \to u + (\tilde\chi_2^0 \to \tilde e_{12}^+ + e^-))$"
    3545histogram inv_mass1 (0,600,20)
    36 
     46}}}
     47{{{
    3748$title = "Factorized process with classical spin correlations"
    3849histogram inv_mass1_dia (0,600,20)
    39 
     50}}}
     51{{{
    4052$title = "Factorized process with isotropic decay"
    4153histogram inv_mass1_iso (0,600,20)
    42 
     54}}}
     55{{{
    4356simulate (full) {
    4457  $sample = "casc_dec_full"
     
    4659    record inv_mass1_full (eval M / 1 GeV [combine[u,e1]])
    4760}
    48 
     61}}}
     62{{{
    4963unstable neu2 (dec_neu_sl2)
    50 
    5164unstable su1 (dec_su_q)
    5265unstable neu2 (dec_neu_sl2)
    53 
    5466simulate (susybg) {
    5567  $sample = "casc_dec"
     
    5769    record inv_mass1 (eval M / 1 GeV [combine[u,e1]])
    5870}
    59 
     71}}}
     72{{{
    6073?diagonal_decay = true
    6174unstable su1 (dec_su_q)
    6275unstable neu2 (dec_neu_sl2)
    63 
    6476simulate (susybg) {
    6577  $sample = "casc_dec_dia"
     
    6779    record inv_mass1_dia (eval M / 1 GeV [combine[u,e1]])
    6880}
    69 
     81}}}
     82{{{
    7083?isotropic_decay = true
    7184unstable su1 (dec_su_q)
    7285unstable neu2 (dec_neu_sl2)
    73 
    7486simulate (susybg) {
    7587  $sample = "casc_dec_iso"
     
    7789    record inv_mass1_iso (eval M / 1 GeV [combine[u,e1]])
    7890}
    79 
     91}}}
     92{{{
    8093write_analysis
    8194$analysis_filename = "casc_dec"