Reweighting does not know set command

Asked by Kristin Lohwasser

Dear all,

I am trying to generate ee (or pp) events with aTGCs using the reweighting method. However this does not work, saying:
INFO: Running Reweight
launch
INFO: detected model: sm. Loading...
INFO: generating the square matrix element for reweighting
INFO: generate p p > e+ ve mu- vm~ ;
INFO: Done 6.009
set Dim6 1 10
WARNING: invalid set command Dim6 1 10
set Dim6 2 0
WARNING: invalid set command Dim6 2 0
set Dim6 3 0
WARNING: invalid set command Dim6 3 0

I am using: MG5_aMC_v2_2_3

and am running the following script:
./bin/mg5_aMC script.mg5

#***********script.py

set group_subprocesses Auto
set ignore_six_quark_processes False
set loop_optimized_output True
set gauge unitary
set complex_mass_scheme False
import model sm
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
define l+ = e+ mu+
define l- = e- mu-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
generate p p > e+ ve mu- vm~
output WW_testpp1

launch
reweight=ON
./reweight_card.dat

set lpp1 1
set lpp2 1

set nevents 1000
set iseed 1234

set ebeam2 3500
set ebeam1 3500

the reweighting

#***********reweight_card.dat
launch
set Dim6 1 10
set Dim6 2 0
set Dim6 3 0

I am supecting, there is something wrong in the script file. It would be great, if you could have a look.
(Also: Is it possible, to get e.g. not EFT operators but also the LEP ones, i.e. gZ and gKappa?)

thanks & best
Kristin

Question information

Language:
English Edit question
Status:
Answered
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Last query:
Last reply:
Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#1

Hi,

The problem is that the model that you use is the sm.
The input_card for the reweighting is then a standard model card which do not define Dim6 block.

You either need to generate your events directly in the model of your choice (EWdim6 I guess)
In a future version, (probably 2.3.2) you will be allowed to use two different model but not for the moment.

> (Also: Is it possible, to get e.g. not EFT operators but also the LEP ones, i.e. gZ and gKappa?)

In 1205.4231, Eq. 12-19 you can find the correspondance between those operator.
As stated in the paper, the basis used at LEP should not be use anymore. If you look at the paper who introduced those operators, they include a remark in that sense as well (that it should not be use at higher energy).
So I will not recommend anyone to implement a new model with those operator.

Cheers,

Olivier

 and therefore the model do not know what the dimension 6 operator are.
This is why the bl

Can you help with this problem?

Provide an answer of your own, or ask Kristin Lohwasser for more information if necessary.

To post a message you must log in.