whizard is hosted by Hepforge, IPPP Durham
close Warning: Error with navigation contributor "BrowserModule"

Changes between Version 2 and Version 3 of ExampleCascadeDecays


Ignore:
Timestamp:
Apr 12, 2010, 10:16:11 AM (15 years ago)
Author:
ohl
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • ExampleCascadeDecays

    v2 v3  
     1Select the minimal supersymmetric standard model
    12{{{
    23model = MSSM
    34}}}
     5and consider a decay cascade of up-squarks into neutralinos and selectrons
    46{{{
    57process dec_su_q        = su1 => u, neu2
    68process dec_neu_sl2     = neu2 => SE12, e1
    79}}}
     10where the up-squark pairs are generated from up quarks
    811{{{
    912process susybg = u,U => SU1, su1
     13}}}
     14We will compare the results with the full process including all irreducible backgrounds
     15{{{
    1016process full = u, U => SU1, u, e1, SE12
    1117}}}
     18Call O'Mega and compile the matrix elements
    1219{{{
    1320compile
    1421}}}
     22We use the SLHA convention to choose the SPS1a parameter point
    1523{{{
    1624?slha_read_decays = true
    1725read_slha("sps1ap_decays.slha")
    1826}}}
     27Obtain the total widths (simple integration, no singularities, not too many points required)
    1928{{{
    2029integrate (dec_su_q, dec_neu_sl2) { iterations = 1:1000 }
    2130}}}
     31We're at the full energy LHC
    2232{{{
    2333sqrts = 14000
    2434beams = p, p => lhapdf
    2535}}}
     36Obtain the up-squark production cross section and the cross section for full process
    2637{{{
    2738integrate (susybg) { iterations = 5:10000, 2:10000 }
    2839integrate (full)
    2940}}}
     41Simulate 10K events
    3042{{{
    3143n_events = 10000
    3244}}}
     45Allocate histograms
    3346{{{
    3447$title = "Full process"
     
    5366histogram inv_mass1_iso (0,600,20)
    5467}}}
     68Simulate the full process
    5569{{{
    5670simulate (full) {
     
    6074}
    6175}}}
     76Simulate the cascade decay, including all off-diagonal spin correlations
    6277{{{
    6378unstable neu2 (dec_neu_sl2)
     
    7085}
    7186}}}
     87Simulate the cascade decay, keeing only the classical (diagonal) spin correlations
    7288{{{
    7389?diagonal_decay = true
     
    8096}
    8197}}}
     98Simulate the cascade decay, dropping all spin correlations (i.e. assume isotropic decays)
    8299{{{
    83100?isotropic_decay = true
     
    90107}
    91108}}}
     109Write the results to file
    92110{{{
    93 write_analysis
    94111$analysis_filename = "casc_dec"
    95112write_analysis