WHIZARD

Beam data: polarization fraction array has wrong dimension

Asked by Yehia Abdelaziz

I have a similar problem as https://answers.launchpad.net/whizard/+question/708455 . I used the solution that you provided {beams = "Eone"} in my case. But the problem is then with the polarization definition.

I have

beams_pol_density= @(-1), @(+1)
beams_pol_fraction = 80%, 30%

It seems that the solution of the mentioned post solves the issue with the beam definition but doesn't provide a solution for their polarization. So I get the following error:

*** FATAL ERROR: Beam data: polarization fraction array has wrong dimension

This is my sindarian commands

model=VLL (ufo)
sqrts = 500 GeV

process ee_e1e1_isr = "e-", "e+" => "Eone", "Eone~"
process Eonedec1 = "Eone" => "mu-" , "X"
process Eonedec2 = "Eone~" => "mu+" , "X"

isr_mass = 511 keV
beams = "e-" , "e+" => circe2 => isr
$circe2_file = "ilc500.circe"
$circe2_design = "ILC"
?circe2_polarized = false

beams_pol_density= @(-1), @(+1)
beams_pol_fraction = 80%, 30%

MVLone = 200 GeV
MX = 100 GeV
integrate(ee_e1e1_isr)
unstable "Eone"(Eonedec1) {beams = "Eone"}
unstable "Eone~"(Eonedec2){beams = "Eone~"}

show(results)

Question information

Language:
English Edit question
Status:
Solved
For:
WHIZARD Edit question
Assignee:
Juergen Reuter Edit question
Solved by:
Yehia Abdelaziz
Solved:
Last query:
Last reply:
Revision history for this message
Juergen Reuter (j.r.reuter) said : 		#1

Hi Yehia,
the solution is written in the issue that you cited: just put the unstable statement _before_ the beams statement of the collider. Then all works fine. In the other case you also need to define local beam polarization settings, along the lines of
unstable "Eone"(Eonedec1) {beams = "Eone" beams_pol_density = @() beams_pol_fraction = 1}
or similar.
Cheers,
    JRR

Revision history for this message
Yehia Abdelaziz (yehia95) said : 		#2

Thank you it now worked. The first option wouldn't work for me. This is why I asked about the second option.