Beams setup does not work

[Juergen Reuter]

Specification of the beams does not work properly (even for pdf_builtin :

[user variable] lepton = PDG(11, 13, 15, -11, -13, -15)
sqrts =  1.000000000000E+03
 Beam data (collision):
   ?  (mass = 0.0000000E+00 GeV)
   ?  (mass = 0.0000000E+00 GeV)
   sqrts = 1.000000000000E+03 GeV
| Beams: ? ?
| Beams: sqrts = 1.000E+03 GeV
| Phase space: generating configuration ...
| Phase space: ... success.
| Phase space: writing configuration file 'drell_yan_i1.phs'
| Phase space: 2 channels, 2 dimensions
| Phase space: wood

[Juergen Reuter]

To make it more precise: it does not work if (1) no beams statement is made (2) directories for PDFs (3) event generation for "partonic" Drell-Yan.

[kilian]

@JR: Could you provide Sindarin code?

[Juergen Reuter]

(2) should be resolved now in r4551. Using this SINDARIN file:

process drell_yan = u, U => e2, E2
compile
sqrts = 1000
integrate (drell_yan)

results in

sqrts =  1.000000000000E+03
 Beam data (collision):
   ?  (mass = 0.0000000E+00 GeV)
   ?  (mass = 0.0000000E+00 GeV)
   sqrts = 1.000000000000E+03 GeV
| Beams: ? ?
| Beams: sqrts = 1.000E+03 GeV
| Phase space: generating configuration ...
| Phase space: ... success.
| Phase space: writing configuration file 'drell_yan_i1.phs'
| Phase space: 2 channels, 2 dimensions
| Phase space: wood
 Evaluator product
 First interaction
 Interaction: 3
 Outgoing:
 Particle 1
 E =   5.000000000000E+02
 P =   0.000000000000E+00  0.000000000000E+00  5.000000000000E+02
 mask [fch] = [FFgT]
 source: (2)1
 Particle 2
 E =   5.000000000000E+02
 P =   0.000000000000E+00  0.000000000000E+00 -5.000000000000E+02
 mask [fch] = [FFgT]
 source: (2)2
 State matrix:  norm =  1.000000000000E+00
 []
   [] => ME(1) = ( 1.000000000000E+00, 0.000000000000E+00)

 Second interaction
 Interaction: 5
 Incoming:
 Particle 1
  [momentum undefined]
 mask [fch] = [FTT]
 internal links: X => 3 4
 source: (4)1
 Particle 2
  [momentum undefined]
 mask [fch] = [FTT]
 internal links: X => 3 4
 source: (4)2
 Outgoing:
 Particle 3
  [momentum undefined]
 mask [fch] = [TTT]
 internal links: 1 2 => X
 source: (4)3
 Particle 4
  [momentum undefined]
 mask [fch] = [TTT]
 internal links: 1 2 => X
 source: (4)4
 State matrix:  norm =  1.000000000000E+00
 [f(2)]
   [f(-2)]
     []
       [] => ME(1) = ( 0.000000000000E+00, 0.000000000000E+00)
 [f(1)]
   [f(-1)]
     []
       [] => ME(2) = ( 0.000000000000E+00, 0.000000000000E+00)
 [f(3)]
   [f(-3)]
     []
       [] => ME(3) = ( 0.000000000000E+00, 0.000000000000E+00)

******************************************************************************
******************************************************************************
*** 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.
******************************************************************************
******************************************************************************
| There were no errors and    1 warning(s).
WHIZARD run aborted.

[kilian]

Really? It works for me:

sqrts =  1.000000000000E+03
 Beam data (collision):
   u     (mass = 0.0000000E+00 GeV)
   ubar  (mass = 0.0000000E+00 GeV)
   sqrts = 1.000000000000E+03 GeV
| Beams: u ubar
| Beams: sqrts = 1.000E+03 GeV
| Phase space: generating configuration ...
| Phase space: ... success.
| Phase space: writing configuration file 'drell_yan_i1.phs'
| Phase space: 2 channels, 2 dimensions
| Phase space: wood
| Starting integration for process 'drell_yan'
| Integrate: iterations not specified, using default
| Integrator: 2 chains, 2 channels, 2 dimensions
| Integrator: VAMP
|=============================================================================|
| It      Calls  Integral[fb]  Error[fb]   Err[%]    Acc  Eff[%]   Chi2 N[It] |
|=============================================================================|
   1        784  1.9309897E+01  3.90E-02    0.20    0.06*   0.00
   2        784  1.9308256E+01  2.79E-02    0.14    0.04*   0.00
   3        784  1.9343524E+01  2.81E-02    0.15    0.04    0.00
|-----------------------------------------------------------------------------|
   3       2352  1.9320559E+01  1.99E-02    0.10    0.05    0.00    0.00   3
|-----------------------------------------------------------------------------|
   4       9936  1.9303810E+01  3.00E-03    0.02    0.02*   0.00
   5       9936  1.9300514E+01  2.99E-03    0.02    0.02*   0.00
   6       9936  1.9303700E+01  3.03E-03    0.02    0.02    0.00
|-----------------------------------------------------------------------------|
   6      29808  1.9302674E+01  1.87E-03    0.01    0.02    0.00    0.00   3
|=============================================================================|
| There were no errors and    1 warning(s).
| WHIZARD run finished.
|=============================================================================|

[Juergen Reuter]

True. Seems that was an old .libs/XXX lingering around. Shouldn't there be a clean mechanism for this maybe? Or are the -r etc. options enough for this. They are not yet working w\ the MD5 sums, right?