Cuts not applied to factorized processes

[dwiesler]

The attached example demonstrates, that cut statements (on final state particles) are not applied, if the process in question is factorized to production and decay. This is easily circumvented by putting an additional if statement in the analysis environment as in the example, but I'm unsure, if this is intended or not. It may confuse the enduser in that given cuts are not applied (in the special case of factorization) and thus it should at least be mentioned in the manual somehow.

[dwiesler]

example sindarin file

[dwiesler]

example pdf showing the missing cut on pt

[Juergen Reuter]

This should maybe be fixed before release 2.0.5.

[Juergen Reuter]

What do we do with this one? Just docu, or is there a possibility to find a solution?

[kilian]

The 'cuts' expression should apply to the hard subprocess, its main use is to make the cross section finite. It would be confusing to apply it to the generated events, again.

I added a 'selection' expression. This is applied only to simulated events (so the matrix element has already been evaluated), where cuts have already been applied. It determines whether an event should be analyzed and/or written to external file. Since the analysis expression can contain conditionals, the rationale for the new feature is mainly to control which events are written to external files.

The 'raw' event files are not affected (all events are there), and events that fail selection will nevertheless count for 'n_events'. The advantage is that one can change the selection without redoing event generation. The disadvantage is a potential waste of resources.

Done in r3196.