whizard is hosted by Hepforge, IPPP Durham

Opened 15 years ago

Last modified 14 years ago

#351 closed defect

WHIZARD gets stuck in a loop when generating new events after reading in old ones — at Initial Version

Reported by: Christian Speckner Owned by: kilian
Priority: P1 Milestone: v2.0.4
Component: core Version: 2.0.2
Severity: critical Keywords:
Cc:

Description

Occasionally, while trying to generate new events after reading in old ones, WHIZARD gets stuck in a loop searching for a valid event. Unfortunately, I have not been able to build a testcase significantly simpler than the actual input file triggering the issue, so I am quoting the full input below. Note that, while the issue seems to arise erratically, it persists once it has surfaced - cancelling the run and retrying does not help. I have looked at the code but have been unable to identify the issue.

model = SM

alias pr = u:U:d:D:c:C:s:S:b:B:g
process drellyan = pr, pr => pr, pr

sqrts = 14 TeV
beams = p, p => lhapdf {$lhapdf_file = "cteq6l.LHpdf"}

mc = 0
ms = 0
mb = 0
mZ = 91.1875 GeV
mW = 80.45 GeV
wW = 0
wZ = 0

real iae = 137.0359895
GF = sqrt (2.) * pi / iae / mW**2 / (1 - mW**2 / mZ**2) / 2

?alpha_s_is_fixed = false
?alpha_s_from_lhapdf = true
scale = (eval Pt [extract index 1 [pr]]) / 2

cuts = all Pt > 100 GeV [pr] and 
   all abs (Eta) < 2.5 [pr]

integrate (drellyan) {iterations = 5:10000, 10:10000}

histogram ptdist (100 GeV, 4 TeV) {n_bins = 40}

analysis = record ptdist (eval Pt [extract index 1 [pr]])
n_events = 1000000
?unweighted = false
simulate (drellyan) {$event_normalization = "sigma/n" checkpoint = 10000}

$out_file = "test"
write_analysis () {?out_custom = true}
!compile_analysis (ptdist) {$out_file = "test.ps"}

Change History (0)

Note: See TracTickets for help on using tickets.