PGS for NLO sample

Asked by Arian Abrahantes on 2014-03-26

Dear experts:

I am trying to run a simple NLO event generation of 5M events

generate p p > t t~ [QCD]

I am using the latest version MG5-v2.1.0 (I have pythia-pgs, delphes, ExRootAnalysis packages installed) and I want to have the full simulation at least PGS reconstructed objects since Delphes is not supported yet.

At a first sight I see that pgs card is present in the generated run directory. But I a m not asked at execution time to perform the ful simulation (i.e. NLO+MadSpin+PGS). The default herwig hep is there after a first simulation test -I know this is becasue herwig is the default shower code-, see the final files list I got in the Events folder:

run_01:
alllogs_0.html alllogs_2.html res_0.txt run_01_tag_1_banner.txt summary.txt
alllogs_1.html events.lhe.gz res_1.txt RunMaterial.tar.gz

run_01_decayed_1:
events_HERWIG6_0.hep.gz events.lhe.gz run_01_decayed_1_tag_1_banner.txt

under this MG versin which is the correct approach to make the full simulation of many events, i.e. passing events_HERWIG6_0.hep.gz through PGS.

1- In the run_card 5M events like:

5000000 = nevents ! Number of unweighted events requested
    -1 = req_acc ! Required accuracy (-1=auto determined from nevents)
    500 = nevt_job! Max number of events per job in event generation.
                 ! (-1= no split).
average = event_norm ! Normalize events to sum or average to the X sect.

2- change the shower:

PYTHIA6Q/PYTHIA6PT = parton_shower

I assume this would make a hep pythia compliant (name-wise) may be with identical name to the one used by pgs.

3- Just rename any of the pgs cards available -> pgs_card.dat (yet pgs in MG is tuned to work with a pythia hep file -I found something alike in the run_pgs executable-, right?)

any help is deeply appreciatted,

cheers,

arian

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
marco zaro Edit question
Solved by:
Olivier Mattelaer
Solved:
2014-03-27
Last query:
2014-03-27
Last reply:
2014-03-26

Hi,

PGS will never work for NLO sample, since it will not support weight associate to the event.

Cheers,

Olivier

On Mar 26, 2014, at 7:01 AM, Arian Abrahantes <email address hidden> wrote:

> New question #246069 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/246069
>
> Dear experts:
>
> I am trying to run a simple NLO event generation of 5M events
>
> generate p p > t t~ [QCD]
>
> I am using the latest version MG5-v2.1.0 (I have pythia-pgs, delphes, ExRootAnalysis packages installed) and I want to have the full simulation at least PGS reconstructed objects since Delphes is not supported yet.
>
> At a first sight I see that pgs card is present in the generated run directory. But I a m not asked at execution time to perform the ful simulation (i.e. NLO+MadSpin+PGS). The default herwig hep is there after a first simulation test -I know this is becasue herwig is the default shower code-, see the final files list I got in the Events folder:
>
> run_01:
> alllogs_0.html alllogs_2.html res_0.txt run_01_tag_1_banner.txt summary.txt
> alllogs_1.html events.lhe.gz res_1.txt RunMaterial.tar.gz
>
> run_01_decayed_1:
> events_HERWIG6_0.hep.gz events.lhe.gz run_01_decayed_1_tag_1_banner.txt
>
> under this MG versin which is the correct approach to make the full simulation of many events, i.e. passing events_HERWIG6_0.hep.gz through PGS.
>
> 1- In the run_card 5M events like:
>
> 5000000 = nevents ! Number of unweighted events requested
> -1 = req_acc ! Required accuracy (-1=auto determined from nevents)
> 500 = nevt_job! Max number of events per job in event generation.
> ! (-1= no split).
> average = event_norm ! Normalize events to sum or average to the X sect.
>
> 2- change the shower:
>
> PYTHIA6Q/PYTHIA6PT = parton_shower
>
> I assume this would make a hep pythia compliant (name-wise) may be with identical name to the one used by pgs.
>
> 3- Just rename any of the pgs cards available -> pgs_card.dat (yet pgs in MG is tuned to work with a pythia hep file -I found something alike in the run_pgs executable-, right?)
>
> any help is deeply appreciatted,
>
> cheers,
>
> arian
>
>
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Dear olivier:

Do this mean taht the full simulation (i.e. detector simulation) for NLO
sample is not supported at all in MG5 and analysis is just possible up to
parton level (analysis of parton LHE events)? So to say in my case
plotting the pt of the top or leptons

In my example:

run_01:
alllogs_0.html alllogs_2.html res_0.txt run_01_tag_1_banner.txt
 summary.txt
alllogs_1.html events.lhe.gz res_1.txt RunMaterial.tar.gz

run_01_decayed_1:
events_HERWIG6_0.hep.gz events.lhe.gz run_01_decayed_1_tag_1_banner.txt

which should I use the "hep" or the "lhe" (I guess the "lhe" in the
"decayed" folder)

Secondly the 5M events simulation at NLO is possible by setting the sample
number (nevents) and splitting afterwards in the run_card?

Thanks for your reply and apologyze me for the spam,

cheers,

arian.

On Wed, Mar 26, 2014 at 2:31 PM, Olivier Mattelaer <
<email address hidden>> wrote:

> Your question #246069 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/246069
>
> Status: Open => Answered
>
> Olivier Mattelaer proposed the following answer:
> Hi,
>
> PGS will never work for NLO sample, since it will not support weight
> associate to the event.
>
> Cheers,
>
> Olivier
>
> On Mar 26, 2014, at 7:01 AM, Arian Abrahantes
> <email address hidden> wrote:
>
> > New question #246069 on MadGraph5_aMC@NLO:
> > https://answers.launchpad.net/mg5amcnlo/+question/246069
> >
> > Dear experts:
> >
> > I am trying to run a simple NLO event generation of 5M events
> >
> > generate p p > t t~ [QCD]
> >
> > I am using the latest version MG5-v2.1.0 (I have pythia-pgs, delphes,
> ExRootAnalysis packages installed) and I want to have the full simulation
> at least PGS reconstructed objects since Delphes is not supported yet.
> >
> > At a first sight I see that pgs card is present in the generated run
> directory. But I a m not asked at execution time to perform the ful
> simulation (i.e. NLO+MadSpin+PGS). The default herwig hep is there after a
> first simulation test -I know this is becasue herwig is the default shower
> code-, see the final files list I got in the Events folder:
> >
> > run_01:
> > alllogs_0.html alllogs_2.html res_0.txt run_01_tag_1_banner.txt
> summary.txt
> > alllogs_1.html events.lhe.gz res_1.txt RunMaterial.tar.gz
> >
> > run_01_decayed_1:
> > events_HERWIG6_0.hep.gz events.lhe.gz run_01_decayed_1_tag_1_banner.txt
> >
> > under this MG versin which is the correct approach to make the full
> simulation of many events, i.e. passing events_HERWIG6_0.hep.gz through PGS.
> >
> > 1- In the run_card 5M events like:
> >
> > 5000000 = nevents ! Number of unweighted events requested
> > -1 = req_acc ! Required accuracy (-1=auto determined from nevents)
> > 500 = nevt_job! Max number of events per job in event generation.
> > ! (-1= no split).
> > average = event_norm ! Normalize events to sum or average to the X sect.
> >
> > 2- change the shower:
> >
> > PYTHIA6Q/PYTHIA6PT = parton_shower
> >
> > I assume this would make a hep pythia compliant (name-wise) may be with
> identical name to the one used by pgs.
> >
> > 3- Just rename any of the pgs cards available -> pgs_card.dat (yet pgs
> in MG is tuned to work with a pythia hep file -I found something alike in
> the run_pgs executable-, right?)
> >
> > any help is deeply appreciatted,
> >
> > cheers,
> >
> > arian
> >
> >
> >
> > --
> > You received this question notification because you are an answer
> > contact for MadGraph5_aMC@NLO.
>
> --
> 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/246069/+confirm?answer_id=0
>
> 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/246069
>
> You received this question notification because you asked the question.
>

Hi Arian,

1) @NLO the lhe file is NOT physical in the sense that some NLO term are missing and therefore it is important to pass by the (correct) shower in order to reach the NLO precision.
So those files should not be used.

2) If PGS is not an option Delphes3 is working properly.

Cheers,

Olivier

On Mar 26, 2014, at 2:06 PM, Arian Abrahantes <email address hidden> wrote:

> Question #246069 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/246069
>
> Status: Answered => Open
>
> Arian Abrahantes is still having a problem:
> Dear olivier:
>
> Do this mean taht the full simulation (i.e. detector simulation) for NLO
> sample is not supported at all in MG5 and analysis is just possible up to
> parton level (analysis of parton LHE events)? So to say in my case
> plotting the pt of the top or leptons
>
> In my example:
>
> run_01:
> alllogs_0.html alllogs_2.html res_0.txt run_01_tag_1_banner.txt
> summary.txt
> alllogs_1.html events.lhe.gz res_1.txt RunMaterial.tar.gz
>
> run_01_decayed_1:
> events_HERWIG6_0.hep.gz events.lhe.gz run_01_decayed_1_tag_1_banner.txt
>
> which should I use the "hep" or the "lhe" (I guess the "lhe" in the
> "decayed" folder)
>
> Secondly the 5M events simulation at NLO is possible by setting the sample
> number (nevents) and splitting afterwards in the run_card?
>
> Thanks for your reply and apologyze me for the spam,
>
> cheers,
>
> arian.
>
>
> On Wed, Mar 26, 2014 at 2:31 PM, Olivier Mattelaer <
> <email address hidden>> wrote:
>
>> Your question #246069 on MadGraph5_aMC@NLO changed:
>> https://answers.launchpad.net/mg5amcnlo/+question/246069
>>
>> Status: Open => Answered
>>
>> Olivier Mattelaer proposed the following answer:
>> Hi,
>>
>> PGS will never work for NLO sample, since it will not support weight
>> associate to the event.
>>
>> Cheers,
>>
>> Olivier
>>
>> On Mar 26, 2014, at 7:01 AM, Arian Abrahantes
>> <email address hidden> wrote:
>>
>>> New question #246069 on MadGraph5_aMC@NLO:
>>> https://answers.launchpad.net/mg5amcnlo/+question/246069
>>>
>>> Dear experts:
>>>
>>> I am trying to run a simple NLO event generation of 5M events
>>>
>>> generate p p > t t~ [QCD]
>>>
>>> I am using the latest version MG5-v2.1.0 (I have pythia-pgs, delphes,
>> ExRootAnalysis packages installed) and I want to have the full simulation
>> at least PGS reconstructed objects since Delphes is not supported yet.
>>>
>>> At a first sight I see that pgs card is present in the generated run
>> directory. But I a m not asked at execution time to perform the ful
>> simulation (i.e. NLO+MadSpin+PGS). The default herwig hep is there after a
>> first simulation test -I know this is becasue herwig is the default shower
>> code-, see the final files list I got in the Events folder:
>>>
>>> run_01:
>>> alllogs_0.html alllogs_2.html res_0.txt run_01_tag_1_banner.txt
>> summary.txt
>>> alllogs_1.html events.lhe.gz res_1.txt RunMaterial.tar.gz
>>>
>>> run_01_decayed_1:
>>> events_HERWIG6_0.hep.gz events.lhe.gz run_01_decayed_1_tag_1_banner.txt
>>>
>>> under this MG versin which is the correct approach to make the full
>> simulation of many events, i.e. passing events_HERWIG6_0.hep.gz through PGS.
>>>
>>> 1- In the run_card 5M events like:
>>>
>>> 5000000 = nevents ! Number of unweighted events requested
>>> -1 = req_acc ! Required accuracy (-1=auto determined from nevents)
>>> 500 = nevt_job! Max number of events per job in event generation.
>>> ! (-1= no split).
>>> average = event_norm ! Normalize events to sum or average to the X sect.
>>>
>>> 2- change the shower:
>>>
>>> PYTHIA6Q/PYTHIA6PT = parton_shower
>>>
>>> I assume this would make a hep pythia compliant (name-wise) may be with
>> identical name to the one used by pgs.
>>>
>>> 3- Just rename any of the pgs cards available -> pgs_card.dat (yet pgs
>> in MG is tuned to work with a pythia hep file -I found something alike in
>> the run_pgs executable-, right?)
>>>
>>> any help is deeply appreciatted,
>>>
>>> cheers,
>>>
>>> arian
>>>
>>>
>>>
>>> --
>>> You received this question notification because you are an answer
>>> contact for MadGraph5_aMC@NLO.
>>
>> --
>> 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/246069/+confirm?answer_id=0
>>
>> 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/246069
>>
>> You received this question notification because you asked the question.
>>
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

dear olivier:

I understand your answer however I am following the tutorial and yet I can
not run delphes for the generated sample (as running thein

tutorial aMCatNLO

I never end in choices like usual MadEvent+Pythia+Delphes

the actual dir generated is "my-first" and by going there and typing:

~/CERN/MG/MG5_aMC_v2_1_0/my-first$ ./bin/generate_events

I get again:

INFO: load configuration from
/home/arian/CERN/MG/MG5_aMC_v2_1_0/my-first/Cards/amcatnlo_configuration.txt
INFO: load configuration from
/home/arian/CERN/MG/MG5_aMC_v2_1_0/input/mg5_configuration.txt
INFO: load configuration from
/home/arian/CERN/MG/MG5_aMC_v2_1_0/my-first/Cards/amcatnlo_configuration.txt
Using default text editor "vi". Set another one in
./input/mg5_configuration.txt
Using default eps viewer "evince". Set another one in
./input/mg5_configuration.txt
Using default web browser "firefox". Set another one in
./input/mg5_configuration.txt
launch
INFO: Enter mode value: Go to the related mode
The following switches determine which operations are executed:
 1 Perturbative order of the calculation:
 order=NLO
 2 Fixed order (no event generation and no MC@[N]LO matching):
fixed_order=OFF
 3 Shower the generated events:
shower=ON
 4 Decay particles with the MadSpin module:
 madspin=ON
  Either type the switch number (1 to 4) to change its default setting,
  or set any switch explicitly (e.g. type 'order=LO' at the prompt)
  Type '0', 'auto', 'done' or just press enter when you are done.

deeply sorry to come over this again but I can not find the way to make it
work neither to find a suitable tutorial on the web for the current MG
version.

thanks in advance for your patience,

arian

On Wed, Mar 26, 2014 at 3:56 PM, Olivier Mattelaer <
<email address hidden>> wrote:

> Your question #246069 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/246069
>
> Status: Open => Answered
>
> Olivier Mattelaer proposed the following answer:
> Hi Arian,
>
> 1) @NLO the lhe file is NOT physical in the sense that some NLO term are
> missing and therefore it is important to pass by the (correct) shower in
> order to reach the NLO precision.
> So those files should not be used.
>
> 2) If PGS is not an option Delphes3 is working properly.
>
> Cheers,
>
> Olivier
>
> On Mar 26, 2014, at 2:06 PM, Arian Abrahantes
> <email address hidden> wrote:
>
> > Question #246069 on MadGraph5_aMC@NLO changed:
> > https://answers.launchpad.net/mg5amcnlo/+question/246069
> >
> > Status: Answered => Open
> >
> > Arian Abrahantes is still having a problem:
> > Dear olivier:
> >
> > Do this mean taht the full simulation (i.e. detector simulation) for NLO
> > sample is not supported at all in MG5 and analysis is just possible up to
> > parton level (analysis of parton LHE events)? So to say in my case
> > plotting the pt of the top or leptons
> >
> > In my example:
> >
> > run_01:
> > alllogs_0.html alllogs_2.html res_0.txt run_01_tag_1_banner.txt
> > summary.txt
> > alllogs_1.html events.lhe.gz res_1.txt RunMaterial.tar.gz
> >
> > run_01_decayed_1:
> > events_HERWIG6_0.hep.gz events.lhe.gz run_01_decayed_1_tag_1_banner.txt
> >
> > which should I use the "hep" or the "lhe" (I guess the "lhe" in the
> > "decayed" folder)
> >
> > Secondly the 5M events simulation at NLO is possible by setting the
> sample
> > number (nevents) and splitting afterwards in the run_card?
> >
> > Thanks for your reply and apologyze me for the spam,
> >
> > cheers,
> >
> > arian.
> >
> >
> > On Wed, Mar 26, 2014 at 2:31 PM, Olivier Mattelaer <
> > <email address hidden>> wrote:
> >
> >> Your question #246069 on MadGraph5_aMC@NLO changed:
> >> https://answers.launchpad.net/mg5amcnlo/+question/246069
> >>
> >> Status: Open => Answered
> >>
> >> Olivier Mattelaer proposed the following answer:
> >> Hi,
> >>
> >> PGS will never work for NLO sample, since it will not support weight
> >> associate to the event.
> >>
> >> Cheers,
> >>
> >> Olivier
> >>
> >> On Mar 26, 2014, at 7:01 AM, Arian Abrahantes
> >> <email address hidden> wrote:
> >>
> >>> New question #246069 on MadGraph5_aMC@NLO:
> >>> https://answers.launchpad.net/mg5amcnlo/+question/246069
> >>>
> >>> Dear experts:
> >>>
> >>> I am trying to run a simple NLO event generation of 5M events
> >>>
> >>> generate p p > t t~ [QCD]
> >>>
> >>> I am using the latest version MG5-v2.1.0 (I have pythia-pgs, delphes,
> >> ExRootAnalysis packages installed) and I want to have the full
> simulation
> >> at least PGS reconstructed objects since Delphes is not supported yet.
> >>>
> >>> At a first sight I see that pgs card is present in the generated run
> >> directory. But I a m not asked at execution time to perform the ful
> >> simulation (i.e. NLO+MadSpin+PGS). The default herwig hep is there
> after a
> >> first simulation test -I know this is becasue herwig is the default
> shower
> >> code-, see the final files list I got in the Events folder:
> >>>
> >>> run_01:
> >>> alllogs_0.html alllogs_2.html res_0.txt run_01_tag_1_banner.txt
> >> summary.txt
> >>> alllogs_1.html events.lhe.gz res_1.txt RunMaterial.tar.gz
> >>>
> >>> run_01_decayed_1:
> >>> events_HERWIG6_0.hep.gz events.lhe.gz
> run_01_decayed_1_tag_1_banner.txt
> >>>
> >>> under this MG versin which is the correct approach to make the full
> >> simulation of many events, i.e. passing events_HERWIG6_0.hep.gz through
> PGS.
> >>>
> >>> 1- In the run_card 5M events like:
> >>>
> >>> 5000000 = nevents ! Number of unweighted events requested
> >>> -1 = req_acc ! Required accuracy (-1=auto determined from nevents)
> >>> 500 = nevt_job! Max number of events per job in event generation.
> >>> ! (-1= no split).
> >>> average = event_norm ! Normalize events to sum or average to the X
> sect.
> >>>
> >>> 2- change the shower:
> >>>
> >>> PYTHIA6Q/PYTHIA6PT = parton_shower
> >>>
> >>> I assume this would make a hep pythia compliant (name-wise) may be with
> >> identical name to the one used by pgs.
> >>>
> >>> 3- Just rename any of the pgs cards available -> pgs_card.dat (yet pgs
> >> in MG is tuned to work with a pythia hep file -I found something alike
> in
> >> the run_pgs executable-, right?)
> >>>
> >>> any help is deeply appreciatted,
> >>>
> >>> cheers,
> >>>
> >>> arian
> >>>
> >>>
> >>>
> >>> --
> >>> You received this question notification because you are an answer
> >>> contact for MadGraph5_aMC@NLO.
> >>
> >> --
> >> 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/246069/+confirm?answer_id=0
> >>
> >> 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/246069
> >>
> >> You received this question notification because you asked the question.
> >>
> >
> > --
> > You received this question notification because you are an answer
> > contact for MadGraph5_aMC@NLO.
>
> --
> 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/246069/+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/246069
>
> You received this question notification because you asked the question.
>

Indeed this option is not proposed at that stage but should be in the future.
I know that that command Delphes is present for the nlo command interface but I let Marco comment on that.

Otherwise you can run mannually Delphes on the output of the shower.

Cheers,

Olivier

unfortunately:

~/CERN/MG/MG5_aMC_v2_1_0/my-first$ ../Delphes/DelphesHepMC
Cards/delphes_card.dat myfirstout.root Events/run_01/events_HERWIG6_0.hep

gives:

** INFO: initializing module FastJetFinder
#--------------------------------------------------------------------------
# FastJet release 3.0.6
# M. Cacciari, G.P. Salam and G. Soyez
# A software package for jet finding and analysis at colliders
# http://fastjet.fr
#
# Please cite EPJC72(2012)1896 [arXiv:1111.6097] if you use this package
# for scientific work and optionally PLB641(2006)57 [hep-ph/0512210].
#
# FastJet is provided without warranty under the terms of the GNU GPLv2.
# It uses T. Chan's closest pair algorithm, S. Fortune's Voronoi code
# and 3rd party plugin jet algorithms. See COPYING file for details.
#--------------------------------------------------------------------------
** INFO: initializing module BTagging
** INFO: initializing module TauTagging
** INFO: initializing module UniqueObjectFinder
** INFO: initializing module ScalarHT
** INFO: initializing module TreeWriter
** Reading Events/run_01/events_HERWIG6_0.hep
** ERROR: invalid units format--------------------------------------]
(0.00%)
** [################################################################]
(100.00%)
** Exiting...

it is something awkward here, should I move to an older v2 version

cheers,

arian

On Wed, Mar 26, 2014 at 6:36 PM, Olivier Mattelaer <
<email address hidden>> wrote:

> Your question #246069 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/246069
>
> Status: Open => Answered
>
> Olivier Mattelaer proposed the following answer:
> Indeed this option is not proposed at that stage but should be in the
> future.
> I know that that command Delphes is present for the nlo command interface
> but I let Marco comment on that.
>
> Otherwise you can run mannually Delphes on the output of the shower.
>
> 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/246069/+confirm?answer_id=4
>
> 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/246069
>
> You received this question notification because you asked the question.
>

I would suggest that you contact the Delphes author for that.
They should be able to help you.

Cheers,

Olivier
On Mar 26, 2014, at 6:01 PM, Arian Abrahantes <email address hidden> wrote:

> Question #246069 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/246069
>
> Status: Answered => Open
>
> Arian Abrahantes is still having a problem:
> unfortunately:
>
> ~/CERN/MG/MG5_aMC_v2_1_0/my-first$ ../Delphes/DelphesHepMC
> Cards/delphes_card.dat myfirstout.root Events/run_01/events_HERWIG6_0.hep
>
> gives:
>
> ** INFO: initializing module FastJetFinder
> #--------------------------------------------------------------------------
> # FastJet release 3.0.6
> # M. Cacciari, G.P. Salam and G. Soyez
> # A software package for jet finding and analysis at colliders
> # http://fastjet.fr
> #
> # Please cite EPJC72(2012)1896 [arXiv:1111.6097] if you use this package
> # for scientific work and optionally PLB641(2006)57 [hep-ph/0512210].
> #
> # FastJet is provided without warranty under the terms of the GNU GPLv2.
> # It uses T. Chan's closest pair algorithm, S. Fortune's Voronoi code
> # and 3rd party plugin jet algorithms. See COPYING file for details.
> #--------------------------------------------------------------------------
> ** INFO: initializing module BTagging
> ** INFO: initializing module TauTagging
> ** INFO: initializing module UniqueObjectFinder
> ** INFO: initializing module ScalarHT
> ** INFO: initializing module TreeWriter
> ** Reading Events/run_01/events_HERWIG6_0.hep
> ** ERROR: invalid units format--------------------------------------]
> (0.00%)
> ** [################################################################]
> (100.00%)
> ** Exiting...
>
> it is something awkward here, should I move to an older v2 version
>
> cheers,
>
> arian
>
>
> On Wed, Mar 26, 2014 at 6:36 PM, Olivier Mattelaer <
> <email address hidden>> wrote:
>
>> Your question #246069 on MadGraph5_aMC@NLO changed:
>> https://answers.launchpad.net/mg5amcnlo/+question/246069
>>
>> Status: Open => Answered
>>
>> Olivier Mattelaer proposed the following answer:
>> Indeed this option is not proposed at that stage but should be in the
>> future.
>> I know that that command Delphes is present for the nlo command interface
>> but I let Marco comment on that.
>>
>> Otherwise you can run mannually Delphes on the output of the shower.
>>
>> 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/246069/+confirm?answer_id=4
>>
>> 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/246069
>>
>> You received this question notification because you asked the question.
>>
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

thanks olivier you always replying as fast as usual. Iwill consider the thread solved and post Delphes author's reply as soon as I receive it.

cheers,

arian

Pavel Demin (pavel-demin) said : #9

Dear Arian,

Could you try the following command?

../Delphes/DelphesSTDHEP Cards/delphes_card.dat myfirstout.root Events/run_01/events_HERWIG6_0.hep

If I'm not mistaken MG5/aMC produces STDHEP files not HepMC.

There is another problem with the HERWIG output:

https://answers.launchpad.net/mg5amcnlo/+question/229096

HERWIG outputs a significant fraction of particles with an unphysical vertex position (NaN).

As far as I know, the problem is still present and the solution is to switch to Pythia for showering.

Regards,

Pavel

Dear pavel:

thanks for your reply, the STDHEP executable did run over HERWIG showered
events. I will take a look to possible differences with PYTHIA. Yet, I have
to say that seems to work very well.

cheers,

arian

On Thu, Mar 27, 2014 at 9:31 AM, Pavel Demin <
<email address hidden>> wrote:

> Your question #246069 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/246069
>
> Pavel Demin posted a new comment:
> Dear Arian,
>
>
> Could you try the following command?
>
> ../Delphes/DelphesSTDHEP Cards/delphes_card.dat myfirstout.root
> Events/run_01/events_HERWIG6_0.hep
>
> If I'm not mistaken MG5/aMC produces STDHEP files not HepMC.
>
>
> There is another problem with the HERWIG output:
>
> https://answers.launchpad.net/mg5amcnlo/+question/229096
>
> HERWIG outputs a significant fraction of particles with an unphysical
> vertex position (NaN).
>
> As far as I know, the problem is still present and the solution is to
> switch to Pythia for showering.
>
>
> Regards,
>
> Pavel
>
> --
> You received this question notification because you asked the question.
>

Thanks Olivier Mattelaer, that solved my question.