Error while generating NLO events with polarized e+ e- beams

Asked by Prudhvi Bhattiprolu

Dear Whizard team,

I am trying to generate events (using whizard for showering with POWHEG matching, pythia8 for hadronization) at NLO (in QCD) for e+ e- > t tbar with polarized incoming beams, but I am getting the following error:

******************************************************************************
*** FATAL ERROR: Product of density matrices is empty
*** --------------------------------------------
*** This happens when two density matrices are convoluted
*** but the processes they belong to (e.g., production
*** and decay) do not match. This could happen if the
*** beam specification does not match the hard
*** process. Or it may indicate a WHIZARD bug.
******************************************************************************
******************************************************************************

The input script "ttNLO.sin" along with all the associated ".log" files are attached here: https://www.dropbox.com/sh/j1istyiiwptaohu/AAAlyCLa0taG3HpljPltrTTBa?dl=0 .

Actually, I was able to generate NLO events for the same process using unpolarized beams, and LO events using polarized beams, so I am not sure if this is a bug?

Any help would be greatly appreciated!

Thanks a lot,
Prudhvi

Question information

Language:
English Edit question
Status:
Solved
For:
WHIZARD Edit question
Assignee:
Pia Bredt Edit question
Solved by:
Pia Bredt
Solved:
Last query:
Last reply:
Revision history for this message
Best Pia Bredt (bredtpia) said :
#1

Hi Prudhvi,

The function you tried, i. e. NLO POWHEG matched events in combination with polarized, unfortunately is not supported in WHIZARD yet.
What works is indeed what you tried, i.e. LO & pol & shower and NLO & unpol & POWHEG matching. Another thing which works is NLO fixed-order simulation (without POWHEG matching) in combination with polarized beams.

Best regards,
Pia Bredt

Revision history for this message
Prudhvi Bhattiprolu (prudhvibhattiprolu) said :
#2

Thanks Pia Bredt, that solved my question.