switch on Delphes and top quark decay

Asked by Winston

Hi experts,
when I generate ttbar samples up to NLO(p p > t t~ [QCD]), why I can't switch on Delphes to do detector simulation as in LO.
Another question, how can I make one top quark decay leptonicly and the other do hadronic decay in simulation? Thank you.

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

Hi,

Indeed we do have a switch to run Delphes from the normal launch.
Now it should be possible if you run the advanced mode
"launch -i" or "./bin/aMCnlo"

Then you have the command delphes available.
Now not all output will have this function working.
As an example all fixed order NLO will not work obviously.
And for aMCNLO mode, I can not guarantee that it work for all parton-shower.

>Another question, how can I make one top quark decay leptonicly and the other do hadronic decay in simulation? Thank you.

For this you need to use MadSpin.

Cheers,

Olivier

Revision history for this message
Winston (adampan) said :
#2

Hi Olivier,
thanks for your help. If I wanna make one top leptonicly and the other do hadronic decay, should I assign their decay products specificly in the line "generate p p > ......." when generate sample or just open the MadSpin when launch?

Best regards.
Winston

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

the assignment of the decay will be inside MadSpin

Cheers,

Olivier

Revision history for this message
Winston (adampan) said :
#4

Hi Olivier, I encounter problems when use Madgraph. there are two errors as following, Can you give me instructions? In error 1 it reminds me "set /path_to_lhapdf/bin/lhapdf-config" from the MadGraph5_aMC@NLO python shell, but where(or in which file) should I set the environment? thank you.

error1:
INFO: /home/renqi/Documents/MG5_aMC_v2_6_3_2/bin/pp_ttbar_nlo/Events/run_02_decayed_1/events.lhe.gz
INFO: MadSpin Done
INFO: Preparing MCatNLO run
WARNING: Trying to shower events using the same PDF in the shower as used in the generation of the events using LHAPDF. However, no valid LHAPDF installation found with the needed PDF set. Will use default internal PDF for the shower instead. To use the same set as was used in the event generation install LHAPDF and set the path using"set /path_to_lhapdf/bin/lhapdf-config" from the MadGraph5_aMC@NLO python shell
INFO: Compiling MCatNLO for PYTHIA8...
 USING PYTHIA 8, version >= 8.200
g++ -O -I/home/renqi/Documents/MG5_aMC_v2_6_3_2/HEPTools/pythia8/include \
  -I/home/renqi/Documents/MG5_aMC_v2_6_3_2/bin/pp_ttbar_nlo/MCatNLO/include Pythia82.cc -o Pythia8.exe \
  -L/home/renqi/Documents/MG5_aMC_v2_6_3_2/HEPTools/pythia8/lib -lpythia8 \
  -I/home/renqi/Documents/MG5_aMC_v2_6_3_2/HEPTools/hepmc/include -L/home/renqi/Documents/MG5_aMC_v2_6_3_2/HEPTools/hepmc/lib -lHepMC -L../lib -L/home/renqi/Documents/MG5_aMC_v2_6_3_2/HEPTools/pythia8//lib -L/home/renqi/Documents/MG5_aMC_v2_6_3_2/HEPTools/boost/lib -L/home/renqi/Documents/MG5_aMC_v2_6_3_2/HEPTools/zlib/lib -lpythia8 -lboost_iostreams -lz -ldl -lstdc++ \
/usr/bin/ld: cannot find -lboost_iostreams
collect2: error: ld returned 1 exit status
Makefile:49: recipe for target 'Pythia82' failed
make: *** [Pythia82] Error 1
Pythia8 compilation did not succeed, exiting

error2:
Pythia8 shower jobs: 0 Idle, 1 Running, 3 Done [23m23s]
Pythia8 shower jobs: 0 Idle, 0 Running, 4 Done [24m20s]
Merging results from the split PY8 runs...
WARNING: Install gnuplot to be able to view the plots generated at :

INFO: Pythia8 shower finished after 25m36s.
INFO: prepare delphes run
INFO: Running Delphes
*** Error in `/home/renqi/Documents/MG5_aMC_v2_6_3_2/Delphes/DelphesHepMC': malloc(): memory corruption: 0x0000000001a78160 ***
======= Backtrace: =========
/lib/x86_64-linux-gnu/libc.so.6(+0x777e5)[0x7f49aca867e5]
/lib/x86_64-linux-gnu/libc.so.6(+0x8213e)[0x7f49aca9113e]
/lib/x86_64-linux-gnu/libc.so.6(__libc_malloc+0x54)[0x7f49aca93184]
/usr/lib/x86_64-linux-gnu/libstdc++.so.6(_Znwm+0x18)[0x7f49ad5a2e78]
/home/renqi/Documents/root6/lib/libCore.so(_ZN5TList7NewLinkEP7TObjectRKSt10shared_ptrI8TObjLinkE+0x5e)[0x7f49afd4c59e]
/home/renqi/Documents/root6/lib/libCore.so(_ZN5TList7AddLastEP7TObject+0x67)[0x7f49afd4cfc7]
/home/renqi/Documents/MG5_aMC_v2_6_3_2/Delphes/DelphesHepMC(_ZN10ExRootTask3AddEP5TTask+0x251)[0x4ffd71]
/home/renqi/Documents/MG5_aMC_v2_6_3_2/Delphes/DelphesHepMC(_ZN7Delphes4InitEv+0x3c0)[0x529060]
/home/renqi/Documents/root6/lib/libCore.so(_ZN5TTask11ExecuteTaskEPKc+0x2d8)[0x7f49afd14118]
/home/renqi/Documents/MG5_aMC_v2_6_3_2/Delphes/DelphesHepMC(main+0x3d6)[0x4bdfa6]
/lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0)[0x7f49aca2f830]
/home/renqi/Documents/MG5_aMC_v2_6_3_2/Delphes/DelphesHepMC(_start+0x29)[0x4be7a9]
======= Memory map: ========

At 2020-04-22 15:22:59, "Olivier Mattelaer" <email address hidden> wrote:
>Your question #690083 on MadGraph5_aMC@NLO changed:
>https://answers.launchpad.net/mg5amcnlo/+question/690083
>
> Status: Open => Answered
>
>Olivier Mattelaer proposed the following answer:
>the assignment of the decay will be inside MadSpin
>
>Cheers,
>
>Olivier
>
>--
>If this answers your question, please go to the following page to let us
>know that it is solved:
>https://answers.launchpad.net/mg5amcnlo/+question/690083/+confirm?answer_id=2
>
>If you still need help, you can reply to this email or go to the
>following page to enter your feedback:
>https://answers.launchpad.net/mg5amcnlo/+question/690083
>
>You received this question notification because you asked the question.

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

Hi,

Try to keep one question per thread, otherwise
1) the title of the topic will not be correct (and other people might miss the information)
2) it is difficult for me to track topic that are still open from those which are close.

For your first error.
Not sure why pythia want to link to "-lboost_iostreams" it might be an old bug which is now fixed.
Could you try with a more recent version of MG5aMC and Pythia8?
Otherwise one solution is to install boost (you can run the command "install boost" inside madgraph terminal)

For your question about lhapdf:
you can type that line in the madgraph terminal as well (the terminal created after you run ./bin/mg5_aMC)
Otherwise you can edit the file input/mg5_configuration.txt)

For your second error. This sounds an issue of Delphes on which I can not comment. My suggestion would be to open a ticket for their git repository.

Cheers,

Olivier

Can you help with this problem?

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

To post a message you must log in.