reweighting lhe.gz file unseccussfully
Dear experts.
I try to reweight the coupling FT8 = 0.01 to FT8 = 0.05 and I use a command at below.
./bin/madevent reweight run_t8_0.01_e90 -f
I see this output on terminal
INFO: Extracting the banner ...
INFO: process: pp pp > z a j j
INFO: options: NP=1
INFO: Running Reweighting
change mode NLO # Define type of Reweighting. For LO sample this command
WARNING: Missing information for NLO type of reweighting. Keeping the mode to LO.
launch
INFO: detected model: SM_LT8_LT9_UFO. Loading...
INFO: generating the square matrix element for reweighting
INFO: generate pp pp > z a j j NP=1 ;
INFO: Done 46.37
INFO: starts to compute weight for events with the following modification to the param_card:
INFO: set param_card anoinputs 19 5e-14 # orig: 1e-14
KEEP:
INFO: Event nb 0 current time: 16h25
INFO: Event nb 10 current time: 16h25
INFO: Event nb 20 0.021s
INFO: Event nb 30 0.03s
INFO: Event nb 40 0.04s
INFO: Event nb 50 0.05s
INFO: Event nb 60 0.06s
INFO: Event nb 70 0.07s
INFO: Event nb 80 0.08s
INFO: Event nb 90 0.09s
INFO: Event nb 100 0.1s
INFO: Event nb 200 0.19s
INFO: Event nb 300 0.29s
INFO: Event nb 400 0.38s
INFO: All event done (nb_event: 500) 0.48s
INFO: Event /Volumes/
INFO: new cross-section is : 49.7201 pb (indicative error: 17.1587 pb)
INFO: Original cross-section: 7.739208 +- 0.1137808 pb (cross-section from sum of weights: 7.739209200000004)
INFO: Computed cross-section:
INFO: rwgt_2 : 49.72009205791405 +- 17.15866617022166 pb
INFO: quit rwgt
quit
INFO: storing files of previous run
gzipping output file: unweighted_
INFO: Done
INFO:
Also I modify the reweight_
#******
# ENTER YOUR COMMANDS BELOW.
#******
change mode NLO # Define type of Reweighting. For LO sample this command
# SPECIFY A PATH OR USE THE SET COMMAND LIKE THIS:
# set sminputs 1 130 # modify 1/alpha_EW
launch
set ANOINPUTS 19 5.00e-14
#set DECAY 23 auto
Then when I check the lhe file, I see former cross-section value like this.
<event>
6 1 +7.7392092e+00 4.14581800e+02 7.81860800e-03 1.04101000e-01
21 -1 0 0 503 502 +0.0000000000e+00 +0.0000000000e+00 +5.1048837771e+02 5.1048837771e+02 0.0000000000e+00 0.0000e+00 1.0000e+00
21 -1 0 0 501 503 -0.0000000000e+00 -0.0000000000e+00 -1.4080348494e+03 1.4080348494e+03 0.0000000000e+00 0.0000e+00 1.0000e+00
23 1 1 2 0 0 -1.6869980018e+02 -6.3069419824e+01 -6.0103487571e+02 6.3403113055e+02 9.1187600000e+01 0.0000e+00 1.0000e+00
22 1 1 2 0 0 +3.4078742388e+01 -8.2016225331e+00 +3.0820325717e+00 3.5187017773e+01 0.0000000000e+00 0.0000e+00 -1.0000e+00
4 1 1 2 501 0 -7.1405943259e+01 -1.9991006297e+02 -6.9543964302e+02 7.2711700509e+02 0.0000000000e+00 0.0000e+00 -1.0000e+00
-4 1 1 2 0 502 +2.0602700105e+02 +2.7118110533e+02 +3.9584601451e+02 5.2218807364e+02 0.0000000000e+00 0.0000e+00 1.0000e+00
<mgrwt>
<rscale> 2 0.40810098E+
<asrwt>0</asrwt>
<pdfrwt beam="1"> 1 21 0.10209771E-01 0.41458180E+
<pdfrwt beam="2"> 1 21 0.28160686E-01 0.41458180E+
<totfact> 0.91965338E+
</mgrwt>
<rwgt>
<wgt id='rwgt_1'> +7.7392092e+00 </wgt>
<wgt id='rwgt_2'> +7.7392092e+00 </wgt>
I see anything in rwgt_2 but it still with same cross-section
Question information
- Language:
- English Edit question
- Status:
- Answered
- Assignee:
- No assignee Edit question
- Last query:
- Last reply:
Can you help with this problem?
Provide an answer of your own, or ask ozgun for more information if necessary.