whizard is hosted by Hepforge, IPPP Durham

Opened 11 years ago

Closed 11 years ago

#555 closed defect (fixed)

Error message for WHIZARD FSR shower

Reported by: Juergen Reuter Owned by: Juergen Reuter
Priority: P0 Milestone: v2.2.0
Component: core Version: 2.1.1
Severity: minor Keywords: Shower interface
Cc: jan.strube@…

Description (last modified by Juergen Reuter)

Reported by Jan Strube: The following SINDARIN input script

model = SM_CKM

mb = 0
process eebb = e1, E1 => b, B
?ps_fsr_active = true
?ps_isr_active = false
?hadronization_active = true
?ps_use_PYTHIA_shower = false

ps_max_n_flavors = 5
ps_mass_cutoff = 1 GeV

cuts = all M > 10 GeV [b,B]

sqrts = 250 GeV
beams = e1, E1 => isr
beams_pol_density = @(1), @(-1)

integrate(eebb) { iterations = 10:50000 }
simulate(eebb) {
    n_events = 10000
}

leads to the strange error message:

STOP Bug: pdf only implemented for (anti-)proton

(adopting to the new 2.2.0 syntax)

Change History (8)

comment:1 Changed 11 years ago by Juergen Reuter

This is triggered by the usage of a structure function (ISR), and presumably the presence of beam remnants in the event (the shower code tries to abuse them as protons).

comment:2 Changed 11 years ago by Juergen Reuter

In Würzburg, we found out that is the usage of a structure function other than pdf_builtin or lhapdf which triggers this error. In this case, it is ISR. Some shower function tries to do something (probably color correlations) with the beam remnants and stumbles over the fact that in this case this clearly isn't a proton.

comment:3 Changed 11 years ago by Juergen Reuter

Description: modified (diff)

comment:4 Changed 11 years ago by Juergen Reuter

By the way, this is awfully bad convergent (probably the problem WK reported for s-channel and ISR).

comment:5 Changed 11 years ago by Juergen Reuter

I can't reproduce this at the moment, on the other hand, I cannot switch PYTHIA on any more. WTF!?

comment:6 Changed 11 years ago by Juergen Reuter

Owner: changed from kilian to Juergen Reuter
Status: newassigned

comment:7 Changed 11 years ago by Juergen Reuter

Priority: P4P0

This is a weirdo, but taking care of it now... well, soon

comment:8 Changed 11 years ago by Juergen Reuter

Resolution: fixed
Status: assignedclosed

Ok, it works with r5497. Happily closing.

Note: See TracTickets for help on using tickets.