kamikaze reweighting + parton shower

Asked by senka duric

Dear experts,
I am using Kamikaze reweighting for EWKdim6NLO model. I generate a QCD NLO sample with aMC@NLO (0,1 jet), get the lhe file, and then I run the reweighting that adds additional weights into existing lhe file.
My problem is when I try to run fxfx matching pythia8 on the output. I run pythia outside of MG framework. If I start from lhe file before the reweighting the fxfx pythia8 matching runs fine, matching efficiency is ~50%. But if I start from lhe file after reweighting then the matching efficiency is 0.
I notice that when I do the reweighting there are other small changes in the lhe file on top of just adding extra weights for each event. For example:

before reweighting the line in lhe file:
<event npLO=' -1 ' npNLO=' 0 '>

after reweighting the line becomes:
<event>

Here is a ship of event before reweighitng:
----------------------------------------->
  <init>
   2212 2212 0.65000000E+04 0.65000000E+04 -1 -1 244600 244600 -4 2
 0.13912745E+01 0.45413123E-02 0.26872144E+01 0
 0.89207793E+00 0.31144248E-02 0.26872144E+01 1
  </init>
  <event npLO=' -1 ' npNLO=' 0 '>
  6 0 0.26872144E+01 0.18015476E+03 0.78186083E-02 0.11840000E+00
       -1 -1 0 0 0 501 0.00000000E+00 0.00000000E+00 0.79755649E+02 0.79756332E+02 0.33000000E+00 0.0000E+00 0.9000E+01
        2 -1 0 0 501 0 0.00000000E+00 0.00000000E+00 -.76154331E+03 0.76154338E+03 0.33000000E+00 0.0000E+00 0.9000E+01
       23 2 1 2 0 0 0.28105565E+02 -.47137718E+01 0.47827842E+02 0.10677752E+03 0.91114181E+02 0.0000E+00 0.0000E+00
       24 1 1 2 0 0 -.28105565E+02 0.47137718E+01 -.72961550E+03 0.73452220E+03 0.79824360E+02 0.0000E+00 0.9000E+01
      -11 1 3 3 0 0 0.21155049E+02 0.38765322E+01 0.74693435E+02 0.77728198E+02 0.51099893E-03 0.0000E+00 0.9000E+01
       11 1 3 3 0 0 0.69505159E+01 -.85903040E+01 -.26865593E+02 0.29049322E+02 0.51099893E-03 0.0000E+00 0.9000E+01
#aMCatNLO 1 0 0 1 2 0.11056197E+03 0.00000000E+00 9 0 0 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
  <rwgt>
   <wgt id='1001'> 0.26872E+01 </wgt>
   <wgt id='1002'> 0.27017E+01 </wgt>
   <wgt id='1003'> 0.26468E+01 </wgt>
   <wgt id='1004'> 0.25993E+01 </wgt>
   <wgt id='1005'> 0.26149E+01 </wgt>
   <wgt id='1006'> 0.25586E+01 </wgt>
   <wgt id='1007'> 0.27964E+01 </wgt>
   <wgt id='1008'> 0.28095E+01 </wgt>
   <wgt id='1009'> 0.27564E+01 </wgt>
  </rwgt>
  </event>
<---------------------------------

after reweighting:
----------------------------------------->
<init>
2212 2212 0.65000000E+04 0.65000000E+04 -1 -1 244600 244600 -4 2
 0.13912745E+01 0.45413123E-02 0.26872144E+01 0
 0.89207793E+00 0.31144248E-02 0.26872144E+01 1
</init>
<event>
 6 0 +2.6872144e+00 1.80154760e+02 7.81860830e-03 1.18400000e-01
       -1 -1 0 0 0 501 +0.0000000000e+00 +0.0000000000e+00 +7.9755649000e+01 7.9756332000e+01 3.3000000000e-01 0.0000e+00 9.0000e+00
        2 -1 0 0 501 0 +0.0000000000e+00 +0.0000000000e+00 -7.6154331000e+02 7.6154338000e+02 3.3000000000e-01 0.0000e+00 9.0000e+00
       23 2 1 2 0 0 +2.8105565000e+01 -4.7137718000e+00 +4.7827842000e+01 1.0677752000e+02 9.1114181000e+01 0.0000e+00 0.0000e+00
       24 1 1 2 0 0 -2.8105565000e+01 +4.7137718000e+00 -7.2961550000e+02 7.3452220000e+02 7.9824360000e+01 0.0000e+00 9.0000e+00
      -11 1 3 3 0 0 +2.1155049000e+01 +3.8765322000e+00 +7.4693435000e+01 7.7728198000e+01 5.1099893000e-04 0.0000e+00 9.0000e+00
       11 1 3 3 0 0 +6.9505159000e+00 -8.5903040000e+00 -2.6865593000e+01 2.9049322000e+01 5.1099893000e-04 0.0000e+00 9.0000e+00
#amcatnlo 1 0 0 1 2 0.11056197e+03 0.00000000e+00 9 0 0 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00
<rwgt>
<wgt id='1001'> +2.6872000e+00 </wgt>
<wgt id='1002'> +2.7017000e+00 </wgt>
<wgt id='1003'> +2.6468000e+00 </wgt>
<wgt id='1004'> +2.5993000e+00 </wgt>
<wgt id='1005'> +2.6149000e+00 </wgt>
<wgt id='1006'> +2.5586000e+00 </wgt>
<wgt id='1007'> +2.7964000e+00 </wgt>
<wgt id='1008'> +2.8095000e+00 </wgt>
<wgt id='1009'> +2.7564000e+00 </wgt>
<wgt id='mg_reweight_1'> +2.6840847e+00 </wgt>
<wgt id='mg_reweight_2'> +2.7235047e+00 </wgt>
<wgt id='mg_reweight_3'> +2.6832529e+00 </wgt>
</rwgt>
</event>
<---------------------------------

I did not notice any important difference in the header of lhe...
Is the zero matching efficiency in pythia8 expected and related to these changes in lhe?
Thank you in advance!

Greetings,
Senka

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Solved by:
senka duric
Solved:
Last query:
Last reply:
Revision history for this message
senka duric (senka-duric) said :
#1

Dear experts,
any news on this?

Greetings,
Senka

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#2

Dear Senka,

Sorry I miss this question.
This should be fixed in the 2.4.0.beta version, if I'm correct. (if not it should be in 2.4.1).
Since this is the exact same problem as another bug report concerning madspin this time.

Cheers,

Olivier

Revision history for this message
senka duric (senka-duric) said :
#3

Thank you Olivier. I see 2.4.0.beta is public, so I can try with that. 2.4.1 is not available yet, right?

Greetings,
Senka

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#4

Dear Senka,

It is available as nightly build via the command
bzr branch lp:~maddevelopers/mg5amcnlo/2.4.1

Cheers,

Olivier

> On Apr 21, 2016, at 18:07, senka duric <email address hidden> wrote:
>
> Question #290935 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/290935
>
> Status: Answered => Open
>
> senka duric is still having a problem:
> Thank you Olivier. I see 2.4.0.beta is public, so I can try with that.
> 2.4.1 is not available yet, right?
>
> Greetings,
> Senka
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Revision history for this message
senka duric (senka-duric) said :
#5

Thank you!

Senka