Problem with beam polarisation in equivalent photon approximation (EPA)

This question was expired because it remained in the 'Open' state without activity for the last 15 days.

Was opened on Christmas. Strange. Nobody reacted^^

Additional hint:
It looks like this behaviour is only observed when circe2 is switched on. If there is no circe2 beam
spectrum, the cross section between the case of beams_pol_fraction = 0%, 0% and the one without
any polarisation setting are identical.

This question was expired because it remained in the 'Open' state without activity for the last 15 days.

We still need some more time.

Filip, I can confirm your findings, using the SINDARIN below for the simple process
gamma e+ -> u dbar nuebar, the unpolarized cross section using 1.4 TeV is 2.5 pb, the
cross section for the polarization combinations are
e-, e+ -+/ ++ 10 pb, while -- and +- are 10^-12 pb. So the average would be the sum of the
four divided by 4, which 20/4 = 5 pb. This is what I get when I do
beams_pol_density = +1,+1
beams_pol_fraction = 0,0
then this cross section corresponds to 5 pb, so the one expected from the average of the
polarized components. The same is true for starting with a forbidden combination of
beams_pol_density = -1,-1
beams_pol_fraction = 0,0

Dear Juergen,

Thank you for checking. All the numbers you give agree with what I get from my simulation.
BUT! The point is that, if you remove the polarisation definition from your script, i.e. remove the lines:
  beams_pol_density = -1,-1
  beams_pol_fraction = 0,0
you will get only 2.5 pb (and not 5pb, as one would expect). Please check!

I would assume that without any polarisation definition the unpolarised cross section should be calculated. But then it should be equal to the cross section for 0%x0% polarisation...

When not using CIRCE2 spectra (monochromatic beam) the cross section for "unpolarised" beams (no polarisation settings) and 0%polarisation are both around 3pb. So it looks like CIRCE2 is doubling the polarised EPA cross section...

  Regards
    Filip

Yes, I also see that not using a beam polariztaion setting, the cross section gets diminuished by a factor of 2. I haven't switched of CIRCE2 yet, but I checked that switching off ISR on the positron hard process leg does not change the picture.

After a more careful inspection (US & WK):

1.
Without CIRCE, one should set
   beams_pol_fraction = 0.0001, 0.0001
to inspect the zero-polarization-fraction case, because otherwise the polarization setup will be recognized as trivial and ignored.

2.
With the "tiny polarization" setting, the problem occurs both with and without CIRCE.

3.
There is no problem for processes that do not involve EPA.

My conclusion: The chaining of the polarization density matrix from the beam with the polarization transfer matrix for EPA is wrong. Most likely this is an actual bug. I'll create a gitlab issue for that and will try to find a solution.