No .root output from simulations

Asked by Giacomo on 2019-07-15

Hello everyone.
May be a stupid question but I cannot find the answer.
I had MG5 but broke some links and everything was falling apart so I decided to wipe out my Mac and restart again. I installed Madgraph, Pythia8, ExRootAnalysis, MadEvent and Delphes.
I want to produce some events at detector level so for example
generate p p > e+ e-

I switch on the Pythia8, Delphes and ExRootAnalysis from the menu after launch (note here that I see "ExRoot" not "ExRootAnalysis" don't know why and if it's correct):

/=====================================================================================\
| 1. Choose the shower/hadronization program shower = Pythia8 |
| 2. Choose the detector simulation program detector = Delphes |
| 3. Choose an analysis package (plot/convert) analysis = ExRoot |
| 4. Decay onshell particles madspin = OFF |
| 5. Add weights to events for new hypp. reweight = OFF |
\=====================================================================================/

The process is now generated and everything works fine:

INFO: Running Pythia8 [arXiv:1410.3012]
Splitting .lhe event file for PY8 parallelization...
Submitting Pythia8 jobs...
Pythia8 shower jobs: 1 Idle, 0 Running, 0 Done [0 second]
Pythia8 shower jobs: 0 Idle, 0 Running, 1 Done [16 seconds]
Merging results from the split PY8 runs...
WARNING: Install gnuplot to be able to view the plots generated at :

INFO: Pythia8 shower finished after 17 seconds.
INFO: prepare delphes run
INFO: Running Delphes
INFO: If you are interested in lhco output. please run root2lhco converter.
INFO: or edit bin/internal/run_delphes3 to run the converter automatically.
INFO: delphes done
  === Results Summary for run: run_01 tag: tag_1 ===

     Cross-section : 7.133e+08 +- 5.576e+06 pb
     Nb of events : 100

INFO: storing files of previous run
INFO: Storing Pythia8 files of previous run
INFO: Done
quit
INFO:
more information in /Users/bcoder/MG5_aMC_v2_6_6/prova/index.html

Except that I don't get any .root file neither at Parton level nor particle level or detector level. I tried different combinations but no output as .root Is produced. I know I can run ExRootLHEFConverter and it works fine on the lhe. But before I wiped my computer the conversion was automatic and I always had a .root in my outputs. It may have something to do with settings I think but I really don't know how to move.
Anyone can help?

Question information

Language:
English Edit question
Status:
Answered
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Last query:
2019-07-15
Last reply:
2019-07-17

Hi,

ExRootAnalysis does not run anymore on my machine, so I can not test it.

but I believe that the following patch should do the trick:

=== modified file 'madgraph/interface/madevent_interface.py'
--- madgraph/interface/madevent_interface.py 2019-06-14 19:54:28 +0000
+++ madgraph/interface/madevent_interface.py 2019-07-17 19:38:41 +0000
@@ -2567,7 +2567,7 @@
                 if self.run_card['time_of_flight']>=0:
                     self.exec_cmd("add_time_of_flight --threshold=%s" % self.run_card['time_of_flight'] ,postcmd=False)

- if switch_mode['analysis'] == 'EXROOTANALYSIS':
+ if switch_mode['analysis'] == 'ExRoot':
                     input = pjoin(self.me_dir, 'Events', self.run_name,'unweighted_events.lhe.gz')
                     output = pjoin(self.me_dir, 'Events', self.run_name, 'unweighted_events.root')
                     self.create_root_file(input , output)

Cheers,

Olivier

> On 15 Jul 2019, at 15:22, Giacomo <email address hidden> wrote:
>
> New question #682051 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/682051
>
> Hello everyone.
> May be a stupid question but I cannot find the answer.
> I had MG5 but broke some links and everything was falling apart so I decided to wipe out my Mac and restart again. I installed Madgraph, Pythia8, ExRootAnalysis, MadEvent and Delphes.
> I want to produce some events at detector level so for example
> generate p p > e+ e-
>
> I switch on the Pythia8, Delphes and ExRootAnalysis from the menu after launch (note here that I see "ExRoot" not "ExRootAnalysis" don't know why and if it's correct):
>
> /=====================================================================================\
> | 1. Choose the shower/hadronization program shower = Pythia8 |
> | 2. Choose the detector simulation program detector = Delphes |
> | 3. Choose an analysis package (plot/convert) analysis = ExRoot |
> | 4. Decay onshell particles madspin = OFF |
> | 5. Add weights to events for new hypp. reweight = OFF |
> \=====================================================================================/
>
> The process is now generated and everything works fine:
>
> INFO: Running Pythia8 [arXiv:1410.3012]
> Splitting .lhe event file for PY8 parallelization...
> Submitting Pythia8 jobs...
> Pythia8 shower jobs: 1 Idle, 0 Running, 0 Done [0 second]
> Pythia8 shower jobs: 0 Idle, 0 Running, 1 Done [16 seconds]
> Merging results from the split PY8 runs...
> WARNING: Install gnuplot to be able to view the plots generated at :
>
> INFO: Pythia8 shower finished after 17 seconds.
> INFO: prepare delphes run
> INFO: Running Delphes
> INFO: If you are interested in lhco output. please run root2lhco converter.
> INFO: or edit bin/internal/run_delphes3 to run the converter automatically.
> INFO: delphes done
> === Results Summary for run: run_01 tag: tag_1 ===
>
> Cross-section : 7.133e+08 +- 5.576e+06 pb
> Nb of events : 100
>
> INFO: storing files of previous run
> INFO: Storing Pythia8 files of previous run
> INFO: Done
> quit
> INFO:
> more information in /Users/bcoder/MG5_aMC_v2_6_6/prova/index.html
>
>
> Except that I don't get any .root file neither at Parton level nor particle level or detector level. I tried different combinations but no output as .root Is produced. I know I can run ExRootLHEFConverter and it works fine on the lhe. But before I wiped my computer the conversion was automatic and I always had a .root in my outputs. It may have something to do with settings I think but I really don't know how to move.
> Anyone can help?
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Can you help with this problem?

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

To post a message you must log in.