whizard is hosted by Hepforge, IPPP Durham

Opened 11 years ago

Closed 11 years ago

#504 closed defect (fixed)

Resurrect simulate

Reported by: Juergen Reuter Owned by: kilian
Priority: P3 Milestone: v2.2.0
Component: core Version: 2.1.1
Severity: normal Keywords:
Cc:

Description

Gives non-matching density matrices at the moment: 

| Starting simulation for process 'lhapdf_drell_yan'
 Evaluator product
 First interaction
 Interaction: 15
 Outgoing:
 Particle 1
  [momentum undefined]
 mask [fch] = [FFgT]
 source: (8)1
 Particle 2
  [momentum undefined]
 mask [fch] = [FFgT]
 source: (8)2
 State matrix:  norm =  1.000000000000E+00
 [f(2) c(3 )]
   [f(-2) c(-4 )] => ME(1) = ( 0.000000000000E+00, 0.000000000000E+00)
   [f(-2) c(-3 )] => ME(2) = ( 0.000000000000E+00, 0.000000000000E+00)

 Second interaction
 Interaction: 13
 Incoming:
 Particle 1
  [momentum undefined]
 mask [fch] = [FFgT]
 internal links: X => 3 4
 source: (9)1
 Particle 2
  [momentum undefined]
 mask [fch] = [FFgT]
 internal links: X => 3 4
 source: (9)2
 Outgoing:
 Particle 3
  [momentum undefined]
 mask [fch] = [FFgT]
 internal links: 1 2 => X
 source: (9)3
 Particle 4
  [momentum undefined]
 mask [fch] = [FFgT]
 internal links: 1 2 => X
 source: (9)4
 State matrix:  norm =  1.000000000000E+00
 [f(2) c(-1 )]
   [f(-2) c(1 )]
     [f(13)]
       [f(-13)] => ME(1) = ( 0.000000000000E+00, 0.000000000000E+00)
 [f(1) c(-1 )]
   [f(-1) c(1 )]
     [f(13)]
       [f(-13)] => ME(2) = ( 0.000000000000E+00, 0.000000000000E+00)
 [f(3) c(-1 )]
   [f(-3) c(1 )]
     [f(13)]
       [f(-13)] => 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.
******************************************************************************
******************************************************************************

Change History (6)

comment:1 Changed 11 years ago by kilian

Err - example code? Does the ticket title match the description?

comment:2 Changed 11 years ago by Juergen Reuter

Also seems a problem with missing defaults. The following SINDARIN file: 

model = "SM"
process drell_yan = u, U => e2, E2
compile
sqrts = 1000
beams = u, U
integrate (drell_yan)
simulate (drell_yan)

results in 

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
|=============================================================================|
| Starting simulation for process 'drell_yan'
 Evaluator product
 First interaction
 Interaction: 15
 Outgoing:
 Particle 1
  [momentum undefined]
 mask [fch] = [FFgT]
 source: (8)1
 Particle 2
  [momentum undefined]
 mask [fch] = [FFgT]
 source: (8)2
 State matrix:  norm =  1.000000000000E+00
 [f(2) c(3 )]
   [f(-2) c(-4 )] => ME(1) = ( 0.000000000000E+00, 0.000000000000E+00)
   [f(-2) c(-3 )] => ME(2) = ( 0.000000000000E+00, 0.000000000000E+00)

 Second interaction
 Interaction: 13
 Incoming:
 Particle 1
  [momentum undefined]
 mask [fch] = [FFgT]
 internal links: X => 3 4
 source: (9)1
 Particle 2
  [momentum undefined]
 mask [fch] = [FFgT]
 internal links: X => 3 4
 source: (9)2
 Outgoing:
 Particle 3
  [momentum undefined]
 mask [fch] = [FFgT]
 internal links: 1 2 => X
 source: (9)3
 Particle 4
  [momentum undefined]
 mask [fch] = [FFgT]
 internal links: 1 2 => X
 source: (9)4
 State matrix:  norm =  1.000000000000E+00
 [f(2) c(-1 )]
   [f(-2) c(1 )]
     [f(13)]
       [f(-13)] => ME(1) = ( 0.000000000000E+00, 0.000000000000E+00)
 [f(1) c(-1 )]
   [f(-1) c(1 )]
     [f(13)]
       [f(-13)] => ME(2) = ( 0.000000000000E+00, 0.000000000000E+00)
 [f(3) c(-1 )]
   [f(-3) c(1 )]
     [f(13)]
       [f(-13)] => 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.

comment:3 Changed 11 years ago by kilian

Status: new → assigned

OK, this problem I can reproduce.

comment:4 Changed 11 years ago by Juergen Reuter

Another question (I hesitate to open a separate ticket): Why are there these funny numbers for the number of calls, at not full numbers like 800 and 10,000. Numerical noise? Is that a bug or a feature?

comment:5 Changed 11 years ago by kilian

@JR: When redoing the implementation, I decided to display the actual number of calls as returned by VAMP. In all previous versions, the program printed just the requested number of calls. In the 'stratified' mode, VAMP chooses a number of calls that fits the current grids and is less or equal than the requested number.

comment:6 Changed 11 years ago by kilian

Resolution: → fixed
Status: assigned → closed

There was a harmless initialization error for the beam colors, but the actual problem was a convention of previous WHIZARD versions: in the hard interaction, the color assignments of incoming particles are inverted. This convention was missed in the reimplementation.

Fixed in r4558.

Note: See TracTickets for help on using tickets.