Placing restrictions
Hello,
for the process e+e- => XX within the model 2HDMS, I want to make sure whizard only takes into one of the higgses decaying into the dark matter particles so I placed the restrictions as: $restrictions = "3+4 ~ h1" and I got,
| Reading model file '/usr/local/
| Preloaded model: SM
| Process library 'default_lib': initialized
| Preloaded library: default_lib
| Reading model file '/usr/local/
| Reading commands from file 'EEAA.sin'
| Model: Generating model 'z2btype2ufo' from UFO sources
| Model: Searching for UFO sources in working directory
| Model: Found UFO sources for model 'z2btype2ufo'
UFO: warning: symbol 'VVS99' <> name 'VVSpp' in lorentz.py: while legal in UFO, it is unusual and can cause problems!
UFO.Lorentz: stripping unused structure ProjP(2,1)
| Model: Model file 'z2btype2ufo.
| Reading model file 'z2btype2ufo.
| Switching to model 'z2btype2ufo' (generated from UFO source)
z2btype2ufo.Ma => 4.800000000000E+01
sqrts = 2.500000000000E+02
z2btype2ufo.betaH => 1.410883200000E+00
$restrictions = "3+4 ~ h1"
| Process library 'default_lib': recorded process 'eeaa'
| Restoring model 'z2btype2ufo'
| Process library 'default_lib': compiling ...
| Process library 'default_lib': writing makefile
| Process library 'default_lib': writing driver
| Process library 'default_lib': creating source code
RM eeaa_i1.
OMEGA[
UFO: warning: symbol 'VVS99' <> name 'VVSpp' in lorentz.py: while legal in UFO, it is unusual and can cause problems!
UFO.Lorentz: stripping unused structure ProjP(2,1)
[1/1] e1 e1bar -> a a ... allowed. [time: 0.02 secs, total: 0.02 secs, remaining: 0.00 secs]
all processes done. [total time: 0.02 secs]
SUMMARY: 4 fusions, 1 propagators, 2 diagrams
SUMMARY: 3 vertices
SUMMARY: 3 UFO vertices: FFS1, FFS2, SSS1
| Process library 'default_lib': compiling sources
FC default_lib.lo
| Process library 'default_lib': linking
FCLD default_lib.la
| Process library 'default_lib': loading
| Process library 'default_lib': ... success.
| RNG: Initializing TAO random-number generator
| RNG: Setting seed for random-number generator to 14477
| Initializing integration for process eeaa:
| Beam structure: [any particles]
| Beam data (collision):
| e1 (mass = 5.1100000E-04 GeV)
| e1bar (mass = 5.1100000E-04 GeV)
| sqrts = 2.500000000000E+02 GeV
| Phase space: discarding old configuration (process changed)
| Phase space: generating configuration ...
Warning: Intermediate decay of zero-width particle h2 may be possible.
| Phase space: ... success.
| Phase space: writing configuration file 'eeaa.i1.phs'
| -------
| Process [scattering]: 'eeaa'
| Library name = 'default_lib'
| Process index = 1
| Process components:
| 1: 'eezh_i1': e1, e1bar => a, a [omega]
| -------
| Phase space: 2 channels, 2 dimensions
| Phase space: found 2 channels, collected in 2 groves.
| Phase space: Using 4 equivalences between channels.
| Phase space: wood
Warning: No cuts have been defined.
| Starting integration for process 'eezh'
| Integrate: iterations not specified, using default
| Integrate: iterations = 3:1000:"gw", 3:10000:""
| Integrator: 2 chains, 2 channels, 2 dimensions
| Integrator: Using VAMP channel equivalences
| Integrator: 1000 initial calls, 20 bins, stratified = T
| Integrator: VAMP
|======
| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] |
|======
| VAMP: parameter mismatch, discarding grid file 'eezh.m1.vg'
1 784 4.8113539E-10 0.00E+00 0.00 0.00* 100.00
2 784 4.8113539E-10 0.00E+00 0.00 0.00 100.00
3 784 4.8113539E-10 0.00E+00 0.00 0.00 100.00
|------
3 2352 4.8113539E-10 0.00E+00 0.00 0.00 100.00 0.00 3
|------
4 9984 4.8113539E-10 2.77E-19 0.00 0.00 100.00
5 9984 4.8113539E-10 2.77E-19 0.00 0.00 100.00
6 9984 4.8113539E-10 2.77E-19 0.00 0.00 100.00
|------
6 29952 4.8113539E-10 1.60E-19 0.00 0.00 100.00 0.00 3
|======
| Time estimate for generating 10000 events: 0d:00h:00m:00s
?fatal_beam_decay = false
isr_q_max = 2.000000000000E+01
isr_alpha = 7.299270072993E-03
epa_q_min = 5.000000000000E+00
epa_x_min = 1.000000000000E-02
| Starting simulation for process 'eeaa'
| Simulate: using integration grids from file 'eeaa.m1.vg'
| RNG: Initializing TAO random-number generator
| RNG: Setting seed for random-number generator to 14478
| Simulation: requested number of events = 0
| corr. to luminosity [fb-1] = 0.0000E+00
| Events: reading from raw file 'eeaa.evx'
| Events: parameter mismatch, discarding old event set
| Events: closing raw file 'eeaa.evx'
| Events: generating new events
| Events: writing to raw file 'eeaa.evx'
| Events: generating 0 unweighted, unpolarized events ...
| Events: event normalization mode '1'
| ... event sample complete.
| Events: closing raw file 'eeaa.evx'
eezh:
4.8113539E-10 +- 1.60E-19 fb ( 0.000 %)
| There were no errors and 2 warning(s).
| WHIZARD run finished.
I wanted to know, is there anything else I should have taken care of while placing the restriction? Thank you.
Question information
- Language:
- English Edit question
- Status:
- Solved
- For:
- WHIZARD Edit question
- Assignee:
- Juergen Reuter Edit question
- Solved by:
- Juergen Reuter
- Solved:
- Last query:
- Last reply: