Loop induced process in Madgraph+madspin

[Keyhan]

Hello,
I generate p p > h [QCD] on Madgraph and use Madspin for the decay h > a a. It works well.
The madspin_card I used is :

--------------------------------------
import model heft /cmsdata3/pghorbani/MG5aMCv242/PPH/Cards/param_card.dat --bypass_check
set spinmode none
#
# set seed 1
# set Nevents_for_max_weigth 75 # number of events for the estimate of the max. weight
# set BW_cut 15 # cut on how far the particle can be off-shell
 set max_weight_ps_point 400 # number of PS to estimate the maximum for each event
#
# specify the decay for the final state particles
decay t > w+ b, w+ > all all
decay t~ > w- b~, w- > all all
decay w+ > all all
decay w- > all all
decay z > all all
decay h > a a
# running the actual code
launch
----------------------------

But, when I generate p p > h xd xd~ in the same model and use madspin for the decay h > a a
I get the following error:
----
INFO: Extracting the banner ...
INFO: process: p p > h xd xd~
INFO: options:
INFO: detected model: Mono-higgs. Loading...
import model heft /cmsdata3/pghorbani/MG5aMCv242/Signal-A2/Cards/param_card.dat --bypass_check
set spinmode none
INFO: detected model: heft. Loading...
INFO: Restrict model heft with file models/heft/restrict_default.dat .
INFO: Run "set stdout_level DEBUG" before import for more information.
Pass the definition of 'j' and 'p' to 4 flavour scheme.
Kept definitions of multiparticles l- / vl / l+ / vl~ unchanged
Removed obsolete multiparticles all
Defined multiparticle all = g u c d s u~ c~ d~ s~ a ve vm vt e- m- ve~ vm~ vt~ e+ m+ t b t~ b~ z w+ h h1 w- tt- tt+
set max_weight_ps_point 400 # number of PS to estimate the maximum for each event
decay t > w+ b, w+ > all all
decay t~ > w- b~, w- > all all
decay w+ > all all
decay w- > all all
decay z > all all
decay h > a a
launch
INFO: Will use seed 566661156
INFO: generate 10000 decay event for particle h
Command "generate_events run_02" interrupted with error:
AttributeError : 'NoneType' object has no attribute 'get_name'
Please report this bug on https://bugs.launchpad.net/mg5amcnlo
More information is found in '/cmsdata3/pghorbani/MG5aMCv242/Signal-A2/run_02_tag_1_debug.log'.
Please attach this file to your report.
INFO: storring files of previous run
gzipping output file: unweighted_events.lhe
----------------------------------------------------

Please note that the link to the param_card is the one belonging to the process p p > h xd xd~.

I would highly appreciate any help.

Best
Keyhan

[Olivier Mattelaer (olivier-mattelaer)]

Hi,

Can you attach the debug file?
> /cmsdata3/pghorbani/MG5aMCv242/Signal-A2/run_02_tag_1_debug.log
This would be very helpful.

Now I have a theory, the problem is that you are modifying the model between the generation and madspin.
Then madspin has no idea of the particle xd (since it uses the new model) and crashes because of that.
Now obviously in “spinmode=None”, we can/should just ignore those particles.

Cheers,

Olivier

> On Nov 21, 2016, at 12:25, Keyhan <email address hidden> wrote:
>
> Public bug reported:
>
> Hello,
> I generate p p > h [QCD] on Madgraph and use Madspin for the decay h > a a. It works well.
> The madspin_card I used is :
>
> --------------------------------------
> import model heft /cmsdata3/pghorbani/MG5aMCv242/PPH/Cards/param_card.dat --bypass_check
> set spinmode none
> #
> # set seed 1
> # set Nevents_for_max_weigth 75 # number of events for the estimate of the max. weight
> # set BW_cut 15 # cut on how far the particle can be off-shell
> set max_weight_ps_point 400 # number of PS to estimate the maximum for each event
> #
> # specify the decay for the final state particles
> decay t > w+ b, w+ > all all
> decay t~ > w- b~, w- > all all
> decay w+ > all all
> decay w- > all all
> decay z > all all
> decay h > a a
> # running the actual code
> launch
> ----------------------------
>
> But, when I generate p p > h xd xd~ in the same model and use madspin for the decay h > a a
> I get the following error:
> ----
> INFO: Extracting the banner ...
> INFO: process: p p > h xd xd~
> INFO: options:
> INFO: detected model: Mono-higgs. Loading...
> import model heft /cmsdata3/pghorbani/MG5aMCv242/Signal-A2/Cards/param_card.dat --bypass_check
> set spinmode none
> INFO: detected model: heft. Loading...
> INFO: Restrict model heft with file models/heft/restrict_default.dat .
> INFO: Run "set stdout_level DEBUG" before import for more information.
> Pass the definition of 'j' and 'p' to 4 flavour scheme.
> Kept definitions of multiparticles l- / vl / l+ / vl~ unchanged
> Removed obsolete multiparticles all
> Defined multiparticle all = g u c d s u~ c~ d~ s~ a ve vm vt e- m- ve~ vm~ vt~ e+ m+ t b t~ b~ z w+ h h1 w- tt- tt+
> set max_weight_ps_point 400 # number of PS to estimate the maximum for each event
> decay t > w+ b, w+ > all all
> decay t~ > w- b~, w- > all all
> decay w+ > all all
> decay w- > all all
> decay z > all all
> decay h > a a
> launch
> INFO: Will use seed 566661156
> INFO: generate 10000 decay event for particle h
> Command "generate_events run_02" interrupted with error:
> AttributeError : 'NoneType' object has no attribute 'get_name'
> Please report this bug on https://bugs.launchpad.net/mg5amcnlo
> More information is found in '/cmsdata3/pghorbani/MG5aMCv242/Signal-A2/run_02_tag_1_debug.log'.
> Please attach this file to your report.
> INFO: storring files of previous run
> gzipping output file: unweighted_events.lhe
> ----------------------------------------------------
>
> Please note that the link to the param_card is the one belonging to the
> process p p > h xd xd~.
>
> I woul...

[Keyhan (kghorbani)]

Hi,
please find the attached debug file.
Best
Keyhan

[Keyhan (kghorbani)]

Hi,

Meanwhile, I tried to use a model for madspin containing the new particle xd, which was in the model used for generation as well. Again, I find the same error as above.
Best
Keyhan

[Olivier Mattelaer (olivier-mattelaer)]

Hi,

Could you send me the model them such that I can try to reproduce your problem and see how to solve it.
Since this is clearly more complicated that was I thought.

Cheers,

Olivier

> On Nov 23, 2016, at 11:42, Keyhan <email address hidden> wrote:
>
> Hi,
>
> Meanwhile, I tried to use a model for madspin containing the new particle xd, which was in the model used for generation as well. Again, I find the same error as above.
> Best
> Keyhan
>
> --
> You received this bug notification because you are subscribed to
> MadGraph5_aMC@NLO.
> https://bugs.launchpad.net/bugs/1643519
>
> Title:
> Loop induced process in Madgraph+madspin
>
> Status in MadGraph5_aMC@NLO:
> New
>
> Bug description:
> Hello,
> I generate p p > h [QCD] on Madgraph and use Madspin for the decay h > a a. It works well.
> The madspin_card I used is :
>
> --------------------------------------
> import model heft /cmsdata3/pghorbani/MG5aMCv242/PPH/Cards/param_card.dat --bypass_check
> set spinmode none
> #
> # set seed 1
> # set Nevents_for_max_weigth 75 # number of events for the estimate of the max. weight
> # set BW_cut 15 # cut on how far the particle can be off-shell
> set max_weight_ps_point 400 # number of PS to estimate the maximum for each event
> #
> # specify the decay for the final state particles
> decay t > w+ b, w+ > all all
> decay t~ > w- b~, w- > all all
> decay w+ > all all
> decay w- > all all
> decay z > all all
> decay h > a a
> # running the actual code
> launch
> ----------------------------
>
> But, when I generate p p > h xd xd~ in the same model and use madspin for the decay h > a a
> I get the following error:
> ----
> INFO: Extracting the banner ...
> INFO: process: p p > h xd xd~
> INFO: options:
> INFO: detected model: Mono-higgs. Loading...
> import model heft /cmsdata3/pghorbani/MG5aMCv242/Signal-A2/Cards/param_card.dat --bypass_check
> set spinmode none
> INFO: detected model: heft. Loading...
> INFO: Restrict model heft with file models/heft/restrict_default.dat .
> INFO: Run "set stdout_level DEBUG" before import for more information.
> Pass the definition of 'j' and 'p' to 4 flavour scheme.
> Kept definitions of multiparticles l- / vl / l+ / vl~ unchanged
> Removed obsolete multiparticles all
> Defined multiparticle all = g u c d s u~ c~ d~ s~ a ve vm vt e- m- ve~ vm~ vt~ e+ m+ t b t~ b~ z w+ h h1 w- tt- tt+
> set max_weight_ps_point 400 # number of PS to estimate the maximum for each event
> decay t > w+ b, w+ > all all
> decay t~ > w- b~, w- > all all
> decay w+ > all all
> decay w- > all all
> decay z > all all
> decay h > a a
> launch
> INFO: Will use seed 566661156
> INFO: generate 10000 decay event for particle h
> Command "generate_events run_02" interrupted with error:
> AttributeError : 'NoneType' object has no attribute 'get_name'
> Please report this bug on https://bugs.launchpad.net/mg5amcnlo
> More information is found in '/cmsdata3/pghorbani/MG5aMCv242/Signal-A2/run_02_tag_1_debug.log'.
> Please attach this file to your report.
> INFO: storring files of previous run
> gzipping output file: unweighted_events.lhe
> -----------...

[Keyhan (kghorbani)]

Hi again,

I attached the model files we use to generate the process: p p > xd xd~ h [QCD] and h > a a.

Best
Keyhan

[Olivier Mattelaer (olivier-mattelaer)]

Ok I reproduce your origninal problem.
I hope to have a fix in less than 24 hours

Cheers,

Olivier

[Olivier Mattelaer (olivier-mattelaer)]

Hi,

In attachement, you will find a patch.
It is also fixing a second problem when madanalysis5 is not installed correctly.

in that case my madspin_card was

set spinmode None
import model heft PROC_heft_0/Cards/param_card.dat --bypass_check
decay h > a a HIW=1
# running the actual code
launch

Cheers,

Olivier

[Keyhan (kghorbani)]

Hi,

Thank you very much for the quick fix.

Have a nice weekends,
Keyhan