Error in extracting lhe.gz file

Hi Daniele,

I got a similar error. This is weird. It looks like the end of the file is actually missin. The produced file can be analyzed with MA5, with the fraction of it that is legit. We will investigate that more deeply. In the meantime, writing an uncompressed file works like a charm (this consists in a quick temporary way over).

Cheers,

Benj

PS: I know I still need to dig into your other e-mail ;)

Hi Benjamin,

Thanks a lot. Actually the problem I have in analyzing this events is that they correspond to a process like

generate p p > X X [NLO]
add process p p > X X j [NLO]

for which my physical xs is the one in summary.txt times the events after pythia (with FxFx merging) over the generated number of events. I do not however know the number of real events after pythia, hence I don't know the cross-section (if what I said above is correct).

With "writing an uncompressed file works like a charm" you mean writing an ASCII with the fraction of events which are readable?

Thanks again

Daniele

Hi Danieie,

You can type
  set main.outputfile = blablabla.lhe
instead. This works, most probably also in your case.

Cheers,

Benj

> On 3 Oct 2018, at 10:03 , Daniele Barducci <email address hidden> wrote:
>
> Question #674599 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/674599
>
> Daniele Barducci posted a new comment:
> Hi Benjamin,
>
> Thanks a lot. Actually the problem I have in analyzing this events is
> that they correspond to a process like
>
> generate p p > X X [NLO]
> add process p p > X X j [NLO]
>
> for which my physical xs is the one in summary.txt times the events
> after pythia (with FxFx merging) over the generated number of events. I
> do not however know the number of real events after pythia, hence I
> don't know the cross-section (if what I said above is correct).
>
> With "writing an uncompressed file works like a charm" you mean writing
> an ASCII with the fraction of events which are readable?
>
> Thanks again
>
> Daniele
>
> --
> You received this question notification because you are subscribed to
> the question.

Thanks, I'll try this!

Hi Benj,

Sorry for the really naive question, but I'm not sure were to write that string. I am using the built-in MG chain of MG+Pythia8+MA5 with a madanalysis_hadron_card.dat which has inside, among the others these strings,

@MG5aMC inputs = *.hepmc, *.hep, *.stdhep, *.lhco, *.fifo

# Reconstruction using FastJet
@MG5aMC reconstruction_name = BasicReco
@MG5aMC reco_output = lhe

As output I get a file which is called

events_PYTHIA8_0_BasicReco.lhe.gz

I tried to play a bit around with the strings along the lines of your suggestion, but I didn't succeed in getting a pure .lhe output. Am I missing something?

Thanks again

Daniele

Hi Daniele,

please try to add:
  set main.outputfile = blablabla.lhe
in the card.

Cheers,

Benjamin

> On 3 Oct 2018, at 13:13 , Daniele Barducci <email address hidden> wrote:
>
> Question #674599 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/674599
>
> Daniele Barducci posted a new comment:
> Hi Benj,
>
> Sorry for the really naive question, but I'm not sure were to write that
> string. I am using the built-in MG chain of MG+Pythia8+MA5 with a
> madanalysis_hadron_card.dat which has inside, among the others these
> strings,
>
> @MG5aMC inputs = *.hepmc, *.hep, *.stdhep, *.lhco, *.fifo
>
> # Reconstruction using FastJet
> @MG5aMC reconstruction_name = BasicReco
> @MG5aMC reco_output = lhe
>
> As output I get a file which is called
>
> events_PYTHIA8_0_BasicReco.lhe.gz
>
> I tried to play a bit around with the strings along the lines of your
> suggestion, but I didn't succeed in getting a pure .lhe output. Am I
> missing something?
>
> Thanks again
>
> Daniele
>
> --
> You received this question notification because you are subscribed to
> the question.

Hi Benj,

It works. Not many events are lost in a test I made (940 vs 944) but the funny things is that the events looks reordered (I didn't check all of them, just the first and last).

In particular the event cut in the lhe.gz is the second event in the file that you told me how to create, and also there are (tiny) difference in the 4-momenta. I don't know if this will help you, but in case I can send the events

Daniele

Hi Daniele,

We are on it. We don't need the events (at least for now) :)

Cheers,

Benj

> On 3 Oct 2018, at 18:07 , Daniele Barducci <email address hidden> wrote:
>
> Question #674599 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/674599
>
> Daniele Barducci posted a new comment:
> Hi Benj,
>
> It works. Not many events are lost in a test I made (940 vs 944) but the
> funny things is that the events looks reordered (I didn't check all of
> them, just the first and last).
>
> In particular the event cut in the lhe.gz is the second event in the
> file that you told me how to create, and also there are (tiny)
> difference in the 4-momenta. I don't know if this will help you, but in
> case I can send the events
>
> Daniele
>
> --
> You received this question notification because you are subscribed to
> the question.

Perfect, thanks a lot

Cheers

Daniele