Ticket #419: ticket419.sin

model = SM
alias quark = u:d:s:c
alias antiq = U:D:S:C
process qqee = quark:antiq, quark:antiq => e1, E1
process qqeeg = quark:antiq, quark:antiq => e1, E1, g
process qgeeq = quark:antiq, g => e1, E1, quark:antiq
process gqeeq = g, quark:antiq => e1, E1, quark:antiq

ms = 0
mc = 0

compile

sqrts = 1960 GeV
beams = p, pbar => lhapdf

$title = "Pt of Z"
$description = "Pt"
$label = "E"
$physical_unit = "GeV"
$xlabel = ""
$ylabel = ""
histogram Pt_distribution_1 (0, 250, 1)
histogram Pt_distribution_2 (0, 250, 1)

cuts = all M > 80 GeV [combine[e1,E1]]
       and all M < 100 GeV [combine[e1,E1]]
       and all Pt > 10 GeV [g:quark:antiq]
       and all Eta < 5 [g:quark:antiq]

seed = 1

integrate(qqee, qqeeg, qgeeq, gqeeq) { iterations = 2:20000 }

n_events = 40

?rebuild_events=true
$sample = "ticket419"
sample_format = lhef
?write_raw = false

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 
!!!!!! for true the correct histogram is filled
! analysis = if 2>1 then
!!!!!! for false expression both histograms stay empty
analysis = if 1>2 then
                record Pt_distribution_1 (1)
        else
                record Pt_distribution_2 (1)
        endif

simulate(qqee, qqeeg, qgeeq, gqeeq)

compile_analysis { $out_file = "ticket419.dat" }

exec("grep entries ticket419.dat")