tHq with MadWeight

Hi,

Here is the patch. I never use MW for single top, but I do not see any problem to use our code in that context.
You might just need to be carefull on the transfer function of the photon since the method should be quite sensitive on it.

Cheers,

Olivier

=== modified file 'madgraph/madweight/write_MadWeight.py'
--- madgraph/madweight/write_MadWeight.py 2014-03-08 18:07:26 +0000
+++ madgraph/madweight/write_MadWeight.py 2017-05-16 16:26:34 +0000
@@ -164,10 +164,12 @@
         data = []

         for part in self.prop_content:
- if part.mass:
+ if part.mass and part.level > 0:
                 mass = self.dict_Fmass[abs(part.pid)]
                 width = self.dict_Fwidth[abs(part.pid)]
                 content = [p.MG for p in get_all_child(part)]
+ if 1 in content or 2 in content:
+ continue
                 content.append(0)
                 initcontent = "\n".join(""" content(%i) = %i""" % (a+1, b) for a,b in enumerate(content))
                 if mass.lower() == 'zero' or width.lower() == 'zero':

Hi Olivier,

Thank you very much for your quick response.

Thanks for fixing the problem with running MadWeight in the latest version of MG5_aMC.

I tried to run MadWeight again by choosing various default options for transfer functions but I still get the same crashes, e.g.:

Fortran runtime error: Index '1' of dimension 1 of array 'e_error' below lower bound of 3

Error termination. Backtrace:
#0 0x114a86729
#1 0x114a873f5
#2 0x114a8779a
#3 0x10eaebeaa
#4 0x10eaec083
#5 0x10ead971d
#6 0x10eb4c003

Would you have any suggestions of what could be tried further in this respect ?

I am now using mg5core2/mg5amcnlo/2.5.5

Many thanks again,

Kirill

Hi,

You still have this with the above patch?
Did you regenerate your directly since applying the above patch?
For me the above patch was fixing the issue.

Cheers,

Olivier

Hi Olivier,

Thank you very much again for your response.

I downloaded the head:

bzr branch lp:~mg5core2/mg5amcnlo/2.5.5

And then regenerated the whole thing by doing:

./bin/mg5_aMC

with:

import model heft
define p = g u c d s b u~ c~ d~ s~ b~
define j = p
define ll = e+ e- mu+ mu- ta+ ta-
define bb = b b~
define vlvl = ve vm vt ve~ vm~ vt~
generate p p > h t j $$ w+ w-, h > a a, t > ll bb vlvl
add process p p > h t~ j $$ w+ w-, h > a a, t~ > ll bb vlvl
output madweight THQMadWeight

And then from the new directory:

bin/madweight.py

I get:

Error termination. Backtrace:
#0 0x109e6f729
#1 0x109e703f5
#2 0x109e7079a
#3 0x103ed2ae3
#4 0x103ed2cbc
#5 0x103ec032d
#6 0x103f32e1d
At line 165 of file permutation_weight.f
Fortran runtime error: Index '1' of dimension 1 of array 'e_error' below lower bound of 3

Example of LHCO events:

# typ eta phi pt jmass ntrk btag had/em dummy dummy
0 0 6
1 0 -0.460059128535 1.56883234104 125.661903614 0.0 1 0 0 0 0
2 0 -0.053990181711 2.73211725586 90.462274526 0.0 1 0 0 0 0
3 1 -0.150253006048 4.2772992407 18.6097556269 0.0 1 0 0 0 0
4 4 0.481140097076 4.40803602669 35.9507387039 4.7 1 2 0 0 0
5 4 -1.97659749898 5.34166438301 125.120496299 0.0 1 0 0 0 0
6 6 -2.6023155368 5.95784122874 29.2478857942 0.0 1 0 0 0 0
# typ eta phi pt jmass ntrk btag had/em dummy dummy
0 1 6
1 0 -1.18708805156 3.38261794543 94.6573981548 0.0 1 0 0 0 0
2 0 -0.909545433148 5.62230546818 49.7487450221 0.0 1 0 0 0 0
3 3 0.588007174048 2.51139499191 16.0872882663 1.777 1 0 0 0 0
4 4 -0.943246567534 0.593037401926 109.110493176 4.7 1 2 0 0 0
5 4 1.81193236175 2.63640507757 29.9130056782 0.0 1 0 0 0 0
6 6 -1.83367958645 4.75452148426 31.8346055217 0.0 1 0 0 0 0

Could you please let me know if there is something else to be done in order to include your recent commit ?

Thanks,

Kirill

Hi again,

I think the problem is that I am not getting your latest commits, I still see the old madgraph/madweight/write_MadWeight.py after doing bzr branch lp:~mg5core2/mg5amcnlo/2.5.5

By any chance, could you please remind me how to get the latest version of the package with your commits ?

Thanks a lot,

Kirill

Hi,

I just push that fix in 2.5.5 branch.
So if you do “bzr pull” you should get it.

Cheers,

Olivier

> On 18 May 2017, at 02:24, Kirill Skovpen <email address hidden> wrote:
>
> Question #632687 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/632687
>
> Status: Answered => Open
>
> Kirill Skovpen is still having a problem:
> Hi Olivier,
>
> Thank you very much again for your response.
>
> I downloaded the head:
>
> bzr branch lp:~mg5core2/mg5amcnlo/2.5.5
>
> And then regenerated the whole thing by doing:
>
> ./bin/mg5_aMC
>
> with:
>
> import model heft
> define p = g u c d s b u~ c~ d~ s~ b~
> define j = p
> define ll = e+ e- mu+ mu- ta+ ta-
> define bb = b b~
> define vlvl = ve vm vt ve~ vm~ vt~
> generate p p > h t j $$ w+ w-, h > a a, t > ll bb vlvl
> add process p p > h t~ j $$ w+ w-, h > a a, t~ > ll bb vlvl
> output madweight THQMadWeight
>
> And then from the new directory:
>
> bin/madweight.py
>
> I get:
>
> Error termination. Backtrace:
> #0 0x109e6f729
> #1 0x109e703f5
> #2 0x109e7079a
> #3 0x103ed2ae3
> #4 0x103ed2cbc
> #5 0x103ec032d
> #6 0x103f32e1d
> At line 165 of file permutation_weight.f
> Fortran runtime error: Index '1' of dimension 1 of array 'e_error' below lower bound of 3
>
> Example of LHCO events:
>
> # typ eta phi pt jmass ntrk btag had/em dummy dummy
> 0 0 6
> 1 0 -0.460059128535 1.56883234104 125.661903614 0.0 1 0 0 0 0
> 2 0 -0.053990181711 2.73211725586 90.462274526 0.0 1 0 0 0 0
> 3 1 -0.150253006048 4.2772992407 18.6097556269 0.0 1 0 0 0 0
> 4 4 0.481140097076 4.40803602669 35.9507387039 4.7 1 2 0 0 0
> 5 4 -1.97659749898 5.34166438301 125.120496299 0.0 1 0 0 0 0
> 6 6 -2.6023155368 5.95784122874 29.2478857942 0.0 1 0 0 0 0
> # typ eta phi pt jmass ntrk btag had/em dummy dummy
> 0 1 6
> 1 0 -1.18708805156 3.38261794543 94.6573981548 0.0 1 0 0 0 0
> 2 0 -0.909545433148 5.62230546818 49.7487450221 0.0 1 0 0 0 0
> 3 3 0.588007174048 2.51139499191 16.0872882663 1.777 1 0 0 0 0
> 4 4 -0.943246567534 0.593037401926 109.110493176 4.7 1 2 0 0 0
> 5 4 1.81193236175 2.63640507757 29.9130056782 0.0 1 0 0 0 0
> 6 6 -1.83367958645 4.75452148426 31.8346055217 0.0 1 0 0 0 0
>
> Could you please let me know if there is something else to be done in
> order to include your recent commit ?
>
> Thanks,
>
> Kirill
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Thanks, Olivier !