Placing restrictions

Asked by Sheikh Farah Tabira

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/share/whizard/models/SM.mdl'
| Preloaded model: SM
| Process library 'default_lib': initialized
| Preloaded library: default_lib
| Reading model file '/usr/local/share/whizard/models/SM_hadrons.mdl'
| 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.ufo.mdl' generated
| Reading model file 'z2btype2ufo.ufo.mdl'
| 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.f90,.mod,.lo
  OMEGA[UFO]eeaa_i1.f90
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:
Revision history for this message
Juergen Reuter (j.r.reuter) said :
#1

Dear Tabira,
thanks for the question. As you are not showing your input file, it is hard to be sure what exactly you are aiming at. I assume that a in our process e1, E1 => a, a is an additional Higgs boson, as well as h1? in any case, by using the restriction `3+4 ~ h1` you are forcing the particles 3+4 in your process to originate from a h1. I.e. you have an s-channel process e+ e- -> h1 -> a a. Somehow I think you want to have a Higgs pair production, where only one of the Higgs bosons decays into your Dark Matter particles. In generating factorized event samples with two identical heavy particles, e1, E1 => X X, and specifying an exclusive decay (X -> yy, assuming 100% BR), both X will be decayed with that decay, so in the events you would get e1, E1 => y, y, y, y. If there are more decay processes, (X -> yy and X -> zz) you would get mixtures with the corresponding final states 4y, 4z, 2y2z with corresponding probabilities. If you have a production of two identical particles and let only one decay, you could do something like e1, E1 => X, y, y { $restrictions = "4+5 ~ X" } and then use a different decay in a factorized approach for the second X decay.
Cheers,
    JRR

Revision history for this message
Sheikh Farah Tabira (farahtabira) said :
#2

Apologies for the lack of clarity. I am looking into the direct production of dark matter particles, and they are denoted by 'a' (SARAH convention). And I want to look at the case where they are being produced via the light higgs (h1 ~ 96gev).

the sindarin script,

model = z2btype2ufo (ufo)

Ma = 48 GeV
sqrts = 250 GeV

betaH = 1.41088320
process eezh = e1 , e1bar => a, a {$restrictions = "3+4 ~ h1"}

compile

integrate (eeaa)

beams = e1 , e1bar !=> isr
!beams_pol_density = @(+1), @(+1)
!beams_pol_fraction = 90%, 30%
?fatal_beam_decay = false

isr_q_max = 20 GeV
isr_alpha = 1./137.
epa_q_min = 5 GeV
epa_x_min = 0.01

!cuts = let subevt @meA = select if Index > 2 [A] in
  !all E > 10 GeV [@meA] and all Theta > 10 degree [@meA]

!!isr_mass = Me1
!!isr_handler = true
!!$isr_handler_mode = "recoil"

simulate (eeaa)
show (results)

Revision history for this message
Best Juergen Reuter (j.r.reuter) said :
#3

Thanks for providing the SINDARIN. This input file looks fine, except for the fact that you of course produce the Higgs boson h1 as an s-channel resonance here. That is the reason why the cross section is so small (10^-10 fb) because it is proportional to the square of the electron Yukawa coupling.

Revision history for this message
Sheikh Farah Tabira (farahtabira) said :
#4

Thanks Juergen Reuter, that solved my question.