Changes between Initial Version and Version 1 of ExampleCascadeDecays

Timestamp:

Apr 12, 2010, 9:44:21 AM (15 years ago)

Author:

ohl

Comment:

--

ExampleCascadeDecays

@@ +1 @@
+{{{
+model = MSSM
+
+process dec_su_q        = su1 => u, neu2
+process dec_neu_sl2     = neu2 => SE12, e1
+
+process susybg = u,U => SU1, su1
+process full = u, U => SU1, u, e1, SE12
+
+compile
+
+?slha_read_decays = true
+read_slha("sps1ap_decays.slha")
+
+integrate (dec_su_q, dec_neu_sl2) { iterations = 1:1000 }
+
+sqrts = 14000
+beams = p, p => lhapdf
+
+integrate (susybg) { iterations = 5:10000, 2:10000 }
+integrate (full)
+
+n_events = 10000
+
+$title = "Full process"
+$description = 
+  "$p + p \to u+ \bar u \to \bar{\tilde u}_1 + u + \tilde e_{12}^+ + e^-$"
+$xlabel = "$M_{\rm inv}(ue^-)$"
+histogram inv_mass1_full (0,600,20)
+
+$title = "Factorized process with complete spin correlations"
+$description = 
+  "$p + p \to u\bar u \to \bar{\tilde u}_1 + 
+  (\tilde u_1 \to u + (\tilde\chi_2^0 \to \tilde e_{12}^+ + e^-))$"
+histogram inv_mass1 (0,600,20)
+
+$title = "Factorized process with classical spin correlations"
+histogram inv_mass1_dia (0,600,20)
+
+$title = "Factorized process with isotropic decay"
+histogram inv_mass1_iso (0,600,20)
+
+simulate (full) {
+  $sample = "casc_dec_full"
+  analysis = 
+    record inv_mass1_full (eval M / 1 GeV [combine[u,e1]])
+}
+
+unstable neu2 (dec_neu_sl2)
+
+unstable su1 (dec_su_q)
+unstable neu2 (dec_neu_sl2)
+
+simulate (susybg) {
+  $sample = "casc_dec"
+  analysis = 
+    record inv_mass1 (eval M / 1 GeV [combine[u,e1]])
+}
+
+?diagonal_decay = true
+unstable su1 (dec_su_q)
+unstable neu2 (dec_neu_sl2)
+
+simulate (susybg) {
+  $sample = "casc_dec_dia"
+  analysis = 
+    record inv_mass1_dia (eval M / 1 GeV [combine[u,e1]])
+}
+
+?isotropic_decay = true
+unstable su1 (dec_su_q)
+unstable neu2 (dec_neu_sl2)
+
+simulate (susybg) {
+  $sample = "casc_dec_iso"
+  analysis = 
+    record inv_mass1_iso (eval M / 1 GeV [combine[u,e1]])
+}
+
+write_analysis
+$analysis_filename = "casc_dec"
+write_analysis
+}}}