NLO QCD for tth process with 5flav scheme

Asked by SUBHASISH BEHERA on 2018-03-26

Dear Oliver,

My appology for not able to check any relevent question related to my error. I have two questions :

Q1. The process I have is as follow :

import model loop_sm
define p = p b b~
define j = p
generate p p > t t~ h [all=QCD]
output pp2tth_all_qcd_5flav
launch pp2tth_all_qcd_5flav
order NLO
fixed_order OFF
shower HERWIG6
#madspin onshell
#reweight
#madanalysis ON
done
#run_card
set run_card ebeam1 7000
set run_card ebeam2 7000
set run_card nevents 50000
set fixed_ren_scale True
set fixed_fac_scale True
set mur_ref_fixed 235.0
set muf_ref_fixed 235.0
set dynamical_scale_choice 0
set mur_over_ref 2.58
set muf_over_ref 2.58
set run_card parton_shower HERWIG6
#showered_card
set shower_card is_bbar T
done

The error is as follow :

Error detected in "launch auto "
write debug file /home/subhasish153.local/RESEARCH/EVENT_GENERATORS/MADGRAPH5/MG5_aMC_v2_6_1/pp2tth_all_qcd_5flav/run_01_tag_1_debug.log
If you need help with this issue please contact us on https://answers.launchpad.net/mg5amcnlo
aMCatNLOError : An error occurred during the collection of results.
 Please check the .log files inside the directories which failed:
 /home/subhasish153.local/RESEARCH/EVENT_GENERATORS/MADGRAPH5/MG5_aMC_v2_6_1/pp2tth_all_qcd_5flav/SubProcesses/P0_bbx_httx/GF1/log.txt
 /home/subhasish153.local/RESEARCH/EVENT_GENERATORS/MADGRAPH5/MG5_aMC_v2_6_1/pp2tth_all_qcd_5flav/SubProcesses/P0_bbx_httx/GF2/log.txt
 /home/subhasish153.local/RESEARCH/EVENT_GENERATORS/MADGRAPH5/MG5_aMC_v2_6_1/pp2tth_all_qcd_5flav/SubProcesses/P0_bbx_httx/GF3/log.txt
 /home/subhasish153.local/RESEARCH/EVENT_GENERATORS/MADGRAPH5/MG5_aMC_v2_6_1/pp2tth_all_qcd_5flav/SubProcesses/P0_bbx_httx/GF4/log.txt
 /home/subhasish153.local/RESEARCH/EVENT_GENERATORS/MADGRAPH5/MG5_aMC_v2_6_1/pp2tth_all_qcd_5flav/SubProcesses/P0_bxb_httx/GF1/log.txt
 /home/subhasish153.local/RESEARCH/EVENT_GENERATORS/MADGRAPH5/MG5_aMC_v2_6_1/pp2tth_all_qcd_5flav/SubProcesses/P0_bxb_httx/GF2/log.txt
 /home/subhasish153.local/RESEARCH/EVENT_GENERATORS/MADGRAPH5/MG5_aMC_v2_6_1/pp2tth_all_qcd_5flav/SubProcesses/P0_bxb_httx/GF3/log.txt
 /home/subhasish153.local/RESEARCH/EVENT_GENERATORS/MADGRAPH5/MG5_aMC_v2_6_1/pp2tth_all_qcd_5flav/SubProcesses/P0_bxb_httx/GF4/log.txt

Q2. I want to use "madspin" for spin analysis here but it has some issue which will be changed in the upcoming version ! Is it ?
You have also mentioned that we can still download from "bzr" depository which I was unable to use. Can you give a link so that I could download and use without installing any further software in my system ?

Thanks,
Subhasish

Question information

Language:
English Edit question
Status:
Answered
For:
MadGraph5_aMC@NLO Edit question
Assignee:
Valentin Hirschi Edit question
Last query:
2018-03-28
Last reply:
2018-03-28

This question was reopened

Dear Oliver,

I have tried with the model : "loop_sm-no_b_mass"; which works fine now. But I couldn't understand why m_b=0 has solved the problem with 5flav scheme ! Can you explain this issue ?

I also tried to switch on "madspin" now but I had issues as follow (This refers to my Q2 above):

Command "launch auto " interrupted with error:
InvalidCmd : "define" command requires symbols "=" at the second position

I have refer to : https://answers.launchpad.net/mg5amcnlo/+question/661988
but unable to get the latest 2.6.2.

Thanks,
Subhasish

Hi,

I'm not working today,
so I will just give you the following link:
https://bazaar.launchpad.net/~mg5core1/mg5amcnlo/2.6.2/revision/322
where you can find a tar ball with (I hope) the full code.

Cheers,

Olivier

Thanks Oliver,

Now after generating process in the new version : https://bazaar.launchpad.net/~mg5core1/mg5amcnlo/2.6.2/revision/322

I got the following error :

INFO: generate add process p p > t t~ h, (t~ > b~ w- , w- > all all QCD=99), (t > b w+ , w+ > all all QCD=99) --no_warning=duplicate;define pert_QCD = -5 -4 -3 -2 -1 1 2 3 4 5 21;add process add process p p > t t~ h pert_QCD, (t~ > b~ w- , w- > all all QCD=99), (t > b w+ , w+ > all all QCD=99) --no_warning=duplicate;
Command "launch auto " interrupted with error:
InvalidCmd : No particle add in model
quit
INFO:

MG5_aMC>

Thanks,
Subhasish

Dear Oliver,

Please ignore my post at #3.

The problem is solved with the new version6.2.2.

Thanks,
Subhasish

Dear Oliver,

Please ignore my post at #3.

The problem is solved with the new version6.2.2.

Thanks,
Subhasish

Dear Oliver,

I am having the following issue now :

INFO: generate p p > t t~ h, (t~ > b~ w- , w- > j j QCD=99), h > b b~ QCD=99, (t > b w+ , w+ > l+ vl QCD=99) --no_warning=duplicate;define pert_QCD = -5 -4 -3 -2 -1 1 2 3 4 5 21;add process p p > t t~ h pert_QCD, (t~ > b~ w- , w- > j j QCD=99), h > b b~ QCD=99, (t > b w+ , w+ > l+ vl QCD=99) --no_warning=duplicate;
Command "launch auto " interrupted with error:
NoDiagramException : No amplitudes generated from process Process: h > b b~ WEIGHTED=4. Please enter a valid process
Please report this bug on https://bugs.launchpad.net/mg5amcnlo
More information is found in '/home/subhasish153.local/RESEARCH/EVENT_GENERATORS/MADGRAPH5/TEST_MG5_aMC_v2_6_2/pp2tth_all_qcd_5flav/run_01_tag_1_debug.log'.
Please attach this file to your report.
quit
INFO:
MG5_aMC>*************************
A coding exception was thrown and uncaught in a Task.

Full message: TypeError: NetworkError when attempting to fetch resource.
Full stack: loadFile@resource://gre/modules/services-common/blocklist-clients.js:119:30
TaskImpl_run@resource://gre/modules/Task.jsm:319:42
TaskImpl@resource://gre/modules/Task.jsm:277:3
asyncFunction@resource://gre/modules/Task.jsm:252:14
Task_spawn@resource://gre/modules/Task.jsm:166:12
loadDumpFile@resource://gre/modules/services-common/blocklist-clients.js:118:12
syncCollection@resource://gre/modules/services-common/blocklist-clients.js:212:39
TaskImpl_run@resource://gre/modules/Task.jsm:319:42

*************************

MG5_aMC>

Thanks,
Subhasish

Hi,

Not sure if this can be done or not inside that model. (this needs a specific handling or the yukawa mass for the renormalization).
Valentin will confirm if that model can do it or not.

Cheers,

Olivier

Dear Olivier,

The first part of execution is all right but when you pass it to madspin, it causes error.
I think we can think of a mechanism where the madspin should take the value from some other card (may be a parameter setting in the madspin_card itself, but not for all may be only for b- quark (like non zero mass to b-quark) or some NP parameters behave differently at decay level.).

Thanks,
Subhasish

Hi,

This already exists but we never check the impact of changing a mass of a particle when doing so.
But if you look at the madspin wiki page you will find such instructions.

Cheers,

Olivier

Can you help with this problem?

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

To post a message you must log in.