MadGraph reweighting only works in some specific cases
Dear MadGraph Group,
I generated the events with an external model which contains a heavy higgs ( hh ), and the process is:
import model SMwithHeavyScal
generate p p > z hh, z > l- l+, hh > j j j j
add process p p > z hh, z > j j, hh > e- e+ j j
add process p p > z hh, z > j j, hh > mu- mu+ j j
add process p p > w- hh, w- > j j, hh > e- e+ j j
add process p p > w- hh, w- > j j, hh > mu- mu+ j j
add process p p > w+ hh, w+ > j j, hh > e- e+ j j
add process p p > w+ hh, w+ > j j, hh > mu- mu+ j j
There are two parameters we can change. They are fw and fww that influence the coupling between heavy higgs and vector boson and the decay width of heavy higgs.
I would like to use MadGraph to reweight the generated events with a different setting of fw and fww from parameter point with fw=1200 and fww=2900. I always reweight from broad width point to narrow width point.
The reweight_card is:
change keep_ordering True
launch --rwgt_
set DIM6COEFF 3 1000
set DIM6COEFF 4 2412
set DECAY 254 auto
launch --rwgt_
set DIM6COEFF 3 800
set DIM6COEFF 4 1930
set DECAY 254 auto
launch --rwgt_
set DIM6COEFF 3 -800
set DIM6COEFF 4 -1930
set DECAY 254 auto
launch --rwgt_
set DIM6COEFF 3 800
set DIM6COEFF 4 0
set DECAY 254 auto
I will use (fw,fww) to repreaent parameter points for convinent. The reweighted cross-section and cross-section calculated by MG5 directly for each points is:
point, reweighted cs[pb] , calculated by MG5 directly[pb]
A: (1000, 2412) 0.000109022232265 +- 7.95662029691e-07, 0.0001079
B: (800,1930) 7.19540848184e-05 +- 9.65657147264e-07 , 7.143e-05
C: (800, 0) 5.54344901268e-05 +- 7.29610350547e-06 , 5.804e-05
D: (-800, -1930) 8.27684121487e-05 +- 9.03141792537e-06 , 8.0754e-05
They are obtained from E:(1200,2900).
I found when use reweighting between parameter points in the same radial direction with respect to (0,0) original point, the reweighting works good. For example, from E to A, B
But when reweighting between parameter points in the different radial direction, there is a big difference between the reweighted cs and and the MG5 cs. And the fluctuation of reweighted cs is very large. For example, from E to C, D.
Could you help me to figure out what cases this? A bug?
The model used is here:
https:/
The process is here:
https:/
Thanks,
Yue
Question information
- Language:
- English Edit question
- Status:
- Solved
- Assignee:
- No assignee Edit question
- Solved by:
- Yue Xu
- Solved:
- Last query:
- Last reply: