Ticket #658: foo.sin

model = SM_Higgs

process pph = g, g => H

compile

sqrts = 13 TeV
beams = p, p => pdf_builtin

histogram h_rap (-5, 5, 0.2)

analysis = record h_rap (eval Rap [H])

n_events = 1000
?unweighted = false
?update_weight = true

sample_format = lhef

$sample = "pph_original"
simulate (pph)
compile_analysis {$out_file=$sample}

clear (analysis)

$sample = "pph_rescan"
rescan "pph_original" (pph) { alt_setup = { $pdf_builtin_set = "MSTW2008LO" } }
compile_analysis {$out_file = $sample}