reweighting lhe.gz file unseccussfully

Asked by ozgun

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/Elements/pp2zajj/FCC_pp_zajj_signal/t8_001_40_2/Events/run_t8_0.01_e90_deneme/unweighted_events.lhe have now the additional weight
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_events.lhe
INFO: Done
INFO:

Also I modify the reweight_card_default.dat like at below.

#************************************************************************
# ENTER YOUR COMMANDS BELOW.
#************************************************************************

change mode NLO # Define type of Reweighting. For LO sample this command
                   # has no effect since only "LO" mode is allowed.

# 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+03</rscale>
<asrwt>0</asrwt>
<pdfrwt beam="1"> 1 21 0.10209771E-01 0.41458180E+03</pdfrwt>
<pdfrwt beam="2"> 1 21 0.28160686E-01 0.41458180E+03</pdfrwt>
<totfact> 0.91965338E+05</totfact>
</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
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

Are you sure that you look at the correct output file?
They are no crash at all, so it should have worked...

Olivier

Revision history for this message
ozgun (ozgunkdz) said (last edit ):
#2

Hi Olivier

My original lhe.gz file is in the Events/run_t8_0.01_e90 directory.

After using the ./bin/madevent reweight run_t8_0.01_e90 -f command, I check it in the same directory and with the same name. If this is wrong, where does it create the new reweighted lhe.gz file? Or do I have to specify it at the beginning when reweighting?

Thanks,

Revision history for this message
Jhoney Bairstow (jhoney12) said :
#3

It appears that after attempting to reweight the coupling FT8 from 0.01 to 0.05, the cross-section value in the lhe file remains the same. The rwgt_2 output suggests that the reweighting process did not result in the expected changes to the cross-section. Further investigation into the reweighting parameters and the command syntax may be needed to achieve the desired adjustments.

Can you help with this problem?

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

To post a message you must log in.