whizard is hosted by Hepforge, IPPP Durham

Opened 14 years ago

Closed 14 years ago

Last modified 14 years ago

#199 closed enhancement (fixed)

Weighted events for decays and in analysis

Reported by: Juergen Reuter Owned by: kilian
Priority: P2 Milestone:
Component: core Version: 2.0.0rc2
Severity: major Keywords: events
Cc:

Description

This is mentioned as a TODO in whizard.nw. Should be (re)enabled in WHIZARD 2.

Change History (9)

comment:1 Changed 14 years ago by Juergen Reuter

Milestone: v2.0.0finalv2.0-rc4
Priority: P3P2

comment:2 Changed 14 years ago by Juergen Reuter

Owner: changed from ALL to Juergen Reuter
Priority: P2P1
Status: newassigned

Will do this, at least for weighted events.

comment:3 Changed 14 years ago by Juergen Reuter

In W1, event weights are normalized between 0 and 1, in W2 they are multiplied by the corresponding cross section. If there is no good physics or computing reason for using the sample function values, there is no good reason for me not to continue the W1 nomenclature. Please, comment!!!!

comment:4 Changed 14 years ago by Juergen Reuter

After discussion with CS I will introduce normalized weights between 0 and 1 in W2 as in W1.

comment:5 Changed 14 years ago by Juergen Reuter

In r2096 I finally activated the weighted events. Now all event formats which write out weights work with weighted events as well (this did not work in W1 !). THe event weights are between 0 and 1, and I abandoned the confusing options of the n_events/n_call interplay in W2. Now n_events is always the number of events that will end up in the event file, both for weighted and unweighted events. For me weighted events seem to be complete, there might be some issues with the analysis when using weighted events, but at least the examples from qedtest.sin are working also with weighted events.

Will try to go to negative weight events as soon as possible.

comment:6 Changed 14 years ago by Juergen Reuter

Priority: P1P2

comment:7 Changed 14 years ago by Juergen Reuter

Owner: changed from Juergen Reuter to kilian
Status: assignednew
Summary: Event generation for weighted and negative weight eventsWeighted events for decays and in analysis

The negative weight events have been separated into ticket #259. This ticket now includes the correct handling of weighted events for decay cascades as well as weighted events in event analyses for e.g. plotting of distributions like differential cross sections.

comment:8 Changed 14 years ago by kilian

Resolution: fixed
Status: newclosed

Thinking over it again, I come to the opinion that decays should be kept unweighted for all modes of event generation. Events weights are used only for the hard interaction, the subsequent decays are selected by rejection. Observables and histograms will come out just right.

This is because decays are treated inclusively when called by the current WHIZARD version. The allowed decay channels are just assumed to saturate the decay width.

It should be possible to treat decays more exclusively, taking into account the full decay width as an independent parameter. But this is the subject of #264.

comment:9 Changed 14 years ago by Juergen Reuter

Milestone: v2.0-rc4

Milestone v2.0-rc4 deleted

Note: See TracTickets for help on using tickets.