Without beam polarisation the calculated cross section is about 4790 fb.
When I set electron beam polarisation the cross section is exactly twice as large (9580 fb) for both -80% and +80% polarisation. Something has to be wrong in my opinion, as the unpolarised cross section should be the average of -80% and +80% cross sections (the positron beam is assumed not to be polarised). The cross section is the same for any polarisation setting. Even if I set it to 0%:

beams_pol_density = @(-1), @(+1)
beams_pol_fraction = 0%, 0%

it is 9580 fb. But when I remove these lines (comment them out) it is back to 4790 fb...

I use same circe2 file in all cases (standard CLIC file). I tested it with version 2.7.0 and 2.8.2 (same results).

I will be grateful for your reply. Which cross section number is the correct one?

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.

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

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...

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.

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.