madanalysis, fastjet and BSM particles

Asked by Luca on 2019-03-16

Hi,

I have the following problem. I am trying to plot the pt distribution of a BSM particle with code 9000005. I generated a sample of 100 events with mg5 and I showered them with pythia8. Everything seems normal if I use the hempc file and I plot the distribution, but something funny seems to happen if I try to use fastjet. I am using the following commands:

set main.fastsim.package = fastjet
set main.fastsim.ptmin = 20.0
set main.fastsim.radius = 0.4
define ax = 9000005
define invisible = invisible 9000005
import /Users/rottoli/codes/MG5_aMC_v2_6_5/axj/Events/run_05/tag_1_pythia8_events.hepmc.gz as TenTeV
plot PT(ax) 25 750 2000 [logY]
submit

The resulting histogram is empty. I tried to have a look at the associated LHE files to see if the particle appears; with the latest version of MadAnalysis, the lhe file appears to be corrupted (though the particle appears in all the events I can see before the file ends abruptly); with an older version (MA5 1.5), the LHE file is fine, the particle appears, but the cross section is zero

Path to the event file
Nr. of events
Cross section (pb)
Negative wgts (%)
ANALYSIS_0/Output/_onetev/mysample.lhe
100
0.0 @ 0.0%
0.0

and therefore the histograms are all still empty. Do you know what could cause this issue?

Thanks,
Luca

Question information

Language:
English Edit question
Status:
Solved
For:
MadAnalysis 5 Edit question
Assignee:
No assignee Edit question
Solved by:
Benjamin Fuks
Solved:
2019-03-22
Last query:
2019-03-22
Last reply:
2019-03-19
Benjamin Fuks (fuks) said : #1

Hi Luca,

In the reco mode, all invisible particles are collected as missing tranverse momentum, so that there is effectively no more '9000005' particle after jet clustering. There is instead an object with a pdg code 12 that corresponds to the MET. This is why your histogram is empty.

The cross section thingie is an old issue that we should fix. In the meantime, feel free to use set TenTeV.xsection = … to manually set the cross section to the correct one.

For the fact that the file is corrupted, do you mind sharing a small hepmc file that allows me to reproduce the issue? Thanks in advance.

Cheers,

Benjanin

> On 15 Mar 2019, at 20:22 , Luca <email address hidden> wrote:
>
> New question #679240 on MadAnalysis 5:
> https://answers.launchpad.net/madanalysis5/+question/679240
>
> Hi,
>
> I have the following problem. I am trying to plot the pt distribution of a BSM particle with code 9000005. I generated a sample of 100 events with mg5 and I showered them with pythia8. Everything seems normal if I use the hempc file and I plot the distribution, but something funny seems to happen if I try to use fastjet. I am using the following commands:
>
> set main.fastsim.package = fastjet
> set main.fastsim.ptmin = 20.0
> set main.fastsim.radius = 0.4
> define ax = 9000005
> define invisible = invisible 9000005
> import /Users/rottoli/codes/MG5_aMC_v2_6_5/axj/Events/run_05/tag_1_pythia8_events.hepmc.gz as TenTeV
> plot PT(ax) 25 750 2000 [logY]
> submit
>
> The resulting histogram is empty. I tried to have a look at the associated LHE files to see if the particle appears; with the latest version of MadAnalysis, the lhe file appears to be corrupted (though the particle appears in all the events I can see before the file ends abruptly); with an older version (MA5 1.5), the LHE file is fine, the particle appears, but the cross section is zero
>
> Path to the event file
> Nr. of events
> Cross section (pb)
> Negative wgts (%)
> ANALYSIS_0/Output/_onetev/mysample.lhe
> 100
> 0.0 @ 0.0%
> 0.0
>
> and therefore the histograms are all still empty. Do you know what could cause this issue?
>
> Thanks,
> Luca
>
>
>
>
>
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Luca (r8li) said : #2

Hi Benjamin,

Thanks for your reply and for the tip, now I am able to fill the histogram setting the cross section to a non-zero value.

You should find at this address
https://www.dropbox.com/s/f7kxieg52uwawmv/tag_1_pythia8_events.hepmc.gz?dl=0

a small hempc file (100 events, ~10mb) which gives back a corrupted the file with the latest version of MadAnalysis. I am not sure if this is relevant, but note that I am using the latest version of fastjet, which I however installed manually as the 'install fastjet' command was failing on my machine.

This last question is not really relevant to the error, but may I also ask if there's a simple way to use a specific event weight from the hepmc file rather than the first one (which I assume is the default?). I had a quick look to the previous answers but I didn't find anything.

Thanks,
Luca

Best Benjamin Fuks (fuks) said : #3

Dear Luca,

The only issue I have found in the outputetd LHE file is that the final xml closing tag was missing. I have fixed it in ma5 v1.8.2. Yiou can either download this version from launchpad:
  https://code.launchpad.net/~ma5dev/madanalysis5/v1.8_beta

In any case, I have managed to create an LHE file from your hepmc sample without any problem. Could you please try with v1.8.2 (if you want a tarball, please send me your e-mail address directly)?

Concerning the fastjet installation, do you mind reytring in debug mode (./bin/ma5 -d) and telling me what is printed to your screen? I will try to fix it.

Cheers,

Benjamin

> On 17 Mar 2019, at 14:08 , Luca <email address hidden> wrote:
>
> Question #679240 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/679240
>
> Status: Answered => Open
>
> Luca is still having a problem:
> Hi Benjamin,
>
> Thanks for your reply and for the tip, now I am able to fill the
> histogram setting the cross section to a non-zero value.
>
> You should find at this address
> https://www.dropbox.com/s/f7kxieg52uwawmv/tag_1_pythia8_events.hepmc.gz?dl=0
>
> a small hempc file (100 events, ~10mb) which gives back a corrupted the
> file with the latest version of MadAnalysis. I am not sure if this is
> relevant, but note that I am using the latest version of fastjet, which
> I however installed manually as the 'install fastjet' command was
> failing on my machine.
>
> This last question is not really relevant to the error, but may I also
> ask if there's a simple way to use a specific event weight from the
> hepmc file rather than the first one (which I assume is the default?). I
> had a quick look to the previous answers but I didn't find anything.
>
> Thanks,
> Luca
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Luca (r8li) said : #4

Dear Benjamin,

Thanks a lot! The lhe files is now totally fine (and I can also gunzip it, which I couldn't do before).

Regarding fastjet, I am not exactly sure what is happening; I can now compile it without errors (in debug mode) both with 1.7 and with the beta version, and I cannot reproduce the error that I had before. I still have the old compilation log which documents the failure (see below the relevant lines), but I am a bit puzzled. Perhaps I had some environmental variable saved in the shell I was using which caused the problem, but at this point I am not sure if it is worth to investigate further...

Thanks,

Luca

Making all in .
/bin/sh ../../libtool --tag=CXX --mode=compile g++ -DHAVE_CONFIG_H -I. -I../../include/fastjet -O2 -Wall -g -Woverloaded-virtual -DDROP_CGAL -I. -I./siscone -I./../../include -I./siscone -MT libSISConePlugin_la-SISConeBasePlugin.lo -MD -MP -MF .deps/libSISConePlugin_la-SISConeBasePlugin.Tpo -c -o libSISConePlugin_la-SISConeBasePlugin.lo `test -f 'SISConeBasePlugin.cc' || echo './'`SISConeBasePlugin.cc
/bin/sh ../../libtool --tag=CXX --mode=compile g++ -DHAVE_CONFIG_H -I. -I../../include/fastjet -O2 -Wall -g -Woverloaded-virtual -DDROP_CGAL -I. -I./siscone -I./../../include -I./siscone -MT libSISConePlugin_la-SISConePlugin.lo -MD -MP -MF .deps/libSISConePlugin_la-SISConePlugin.Tpo -c -o libSISConePlugin_la-SISConePlugin.lo `test -f 'SISConePlugin.cc' || echo './'`SISConePlugin.cc
/bin/sh ../../libtool --tag=CXX --mode=compile g++ -DHAVE_CONFIG_H -I. -I../../include/fastjet -O2 -Wall -g -Woverloaded-virtual -DDROP_CGAL -I. -I./siscone -I./../../include -I./siscone -MT libSISConePlugin_la-SISConeSphericalPlugin.lo -MD -MP -MF .deps/libSISConePlugin_la-SISConeSphericalPlugin.Tpo -c -o libSISConePlugin_la-SISConeSphericalPlugin.lo `test -f 'SISConeSphericalPlugin.cc' || echo './'`SISConeSphericalPlugin.cc
libtool: compile: g++ -DHAVE_CONFIG_H -I. -I../../include/fastjet -O2 -Wall -g -Woverloaded-virtual -DDROP_CGAL -I. -I./siscone -I./../../include -I./siscone -MT libSISConePlugin_la-SISConePlugin.lo -MD -MP -MF .deps/libSISConePlugin_la-SISConePlugin.Tpo -c SISConePlugin.cc -fno-common -DPIC -o .libs/libSISConePlugin_la-SISConePlugin.o
libtool: compile: g++ -DHAVE_CONFIG_H -I. -I../../include/fastjet -O2 -Wall -g -Woverloaded-virtual -DDROP_CGAL -I. -I./siscone -I./../../include -I./siscone -MT libSISConePlugin_la-SISConeBasePlugin.lo -MD -MP -MF .deps/libSISConePlugin_la-SISConeBasePlugin.Tpo -c SISConeBasePlugin.cc -fno-common -DPIC -o .libs/libSISConePlugin_la-SISConeBasePlugin.o
libtool: compile: g++ -DHAVE_CONFIG_H -I. -I../../include/fastjet -O2 -Wall -g -Woverloaded-virtual -DDROP_CGAL -I. -I./siscone -I./../../include -I./siscone -MT libSISConePlugin_la-SISConeSphericalPlugin.lo -MD -MP -MF .deps/libSISConePlugin_la-SISConeSphericalPlugin.Tpo -c SISConeSphericalPlugin.cc -fno-common -DPIC -o .libs/libSISConePlugin_la-SISConeSphericalPlugin.o
SISConeBasePlugin.cc:12:12: error: no matching member function for call to 'structure_of'
  return a.structure_of<UserScaleBase::StructureType>().ordering_var2()
         ~~^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
./../../include/fastjet/PseudoJet.hh:1006:60: note: candidate template ignored: substitution failure [with TransformerType = fastjet::SISConeBasePlugin::UserScaleBase::StructureType]: ISO C++ specifies that qualified reference to 'StructureType' is a constructor name rather than a type in this context, despite preceding 'typename' keyword
const typename TransformerType::StructureType & PseudoJet::structure_of() const{
      ~~~~~~~~ ^
SISConeBasePlugin.cc:13:12: error: no matching member function for call to 'structure_of'
       > b.structure_of<UserScaleBase::StructureType>().ordering_var2();
         ~~^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
./../../include/fastjet/PseudoJet.hh:1006:60: note: candidate template ignored: substitution failure [with TransformerType = fastjet::SISConeBasePlugin::UserScaleBase::StructureType]: ISO C++ specifies that qualified reference to 'StructureType' is a constructor name rather than a type in this context, despite preceding 'typename' keyword
const typename TransformerType::StructureType & PseudoJet::structure_of() const{
      ~~~~~~~~ ^
2 errors generated.
make[3]: *** [libSISConePlugin_la-SISConeBasePlugin.lo] Error 1
make[3]: *** Waiting for unfinished jobs....
libtool: compile: g++ -DHAVE_CONFIG_H -I. -I../../include/fastjet -O2 -Wall -g -Woverloaded-virtual -DDROP_CGAL -I. -I./siscone -I./../../include -I./siscone -MT libSISConePlugin_la-SISConePlugin.lo -MD -MP -MF .deps/libSISConePlugin_la-SISConePlugin.Tpo -c SISConePlugin.cc -o libSISConePlugin_la-SISConePlugin.o >/dev/null 2>&1
libtool: compile: g++ -DHAVE_CONFIG_H -I. -I../../include/fastjet -O2 -Wall -g -Woverloaded-virtual -DDROP_CGAL -I. -I./siscone -I./../../include -I./siscone -MT libSISConePlugin_la-SISConeSphericalPlugin.lo -MD -MP -MF .deps/libSISConePlugin_la-SISConeSphericalPlugin.Tpo -c SISConeSphericalPlugin.cc -o libSISConePlugin_la-SISConeSphericalPlugin.o >/dev/null 2>&1
mv -f .deps/libSISConePlugin_la-SISConePlugin.Tpo .deps/libSISConePlugin_la-SISConePlugin.Plo
mv -f .deps/libSISConePlugin_la-SISConeSphericalPlugin.Tpo .deps/libSISConePlugin_la-SISConeSphericalPlugin.Plo
make[2]: *** [all-recursive] Error 1
make[1]: *** [all-recursive] Error 1
make: *** [all-recursive] Error 1

Luca (r8li) said : #5

Thanks Benjamin Fuks, that solved my question.

Benjamin Fuks (fuks) said : #6

Hi Luca,

it looks like 2 fastjet installations (an old and a new) have been found and madanalysis5 was mixing them. However, without reproducing exactly the issue, it is hard to tell more. Anyways, I am happy to read that the problem is solved.

Cheers,

Benjamin

> On 21 Mar 2019, at 06:44 , Luca <email address hidden> wrote:
>
> Question #679240 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/679240
>
> Status: Answered => Solved
>
> Luca confirmed that the question is solved:
> Dear Benjamin,
>
> Thanks a lot! The lhe files is now totally fine (and I can also gunzip
> it, which I couldn't do before).
>
> Regarding fastjet, I am not exactly sure what is happening; I can now
> compile it without errors (in debug mode) both with 1.7 and with the
> beta version, and I cannot reproduce the error that I had before. I
> still have the old compilation log which documents the failure (see
> below the relevant lines), but I am a bit puzzled. Perhaps I had some
> environmental variable saved in the shell I was using which caused the
> problem, but at this point I am not sure if it is worth to investigate
> further...
>
> Thanks,
>
> Luca
>
>
> Making all in .
> /bin/sh ../../libtool --tag=CXX --mode=compile g++ -DHAVE_CONFIG_H -I. -I../../include/fastjet -O2 -Wall -g -Woverloaded-virtual -DDROP_CGAL -I. -I./siscone -I./../../include -I./siscone -MT libSISConePlugin_la-SISConeBasePlugin.lo -MD -MP -MF .deps/libSISConePlugin_la-SISConeBasePlugin.Tpo -c -o libSISConePlugin_la-SISConeBasePlugin.lo `test -f 'SISConeBasePlugin.cc' || echo './'`SISConeBasePlugin.cc
> /bin/sh ../../libtool --tag=CXX --mode=compile g++ -DHAVE_CONFIG_H -I. -I../../include/fastjet -O2 -Wall -g -Woverloaded-virtual -DDROP_CGAL -I. -I./siscone -I./../../include -I./siscone -MT libSISConePlugin_la-SISConePlugin.lo -MD -MP -MF .deps/libSISConePlugin_la-SISConePlugin.Tpo -c -o libSISConePlugin_la-SISConePlugin.lo `test -f 'SISConePlugin.cc' || echo './'`SISConePlugin.cc
> /bin/sh ../../libtool --tag=CXX --mode=compile g++ -DHAVE_CONFIG_H -I. -I../../include/fastjet -O2 -Wall -g -Woverloaded-virtual -DDROP_CGAL -I. -I./siscone -I./../../include -I./siscone -MT libSISConePlugin_la-SISConeSphericalPlugin.lo -MD -MP -MF .deps/libSISConePlugin_la-SISConeSphericalPlugin.Tpo -c -o libSISConePlugin_la-SISConeSphericalPlugin.lo `test -f 'SISConeSphericalPlugin.cc' || echo './'`SISConeSphericalPlugin.cc
> libtool: compile: g++ -DHAVE_CONFIG_H -I. -I../../include/fastjet -O2 -Wall -g -Woverloaded-virtual -DDROP_CGAL -I. -I./siscone -I./../../include -I./siscone -MT libSISConePlugin_la-SISConePlugin.lo -MD -MP -MF .deps/libSISConePlugin_la-SISConePlugin.Tpo -c SISConePlugin.cc -fno-common -DPIC -o .libs/libSISConePlugin_la-SISConePlugin.o
> libtool: compile: g++ -DHAVE_CONFIG_H -I. -I../../include/fastjet -O2 -Wall -g -Woverloaded-virtual -DDROP_CGAL -I. -I./siscone -I./../../include -I./siscone -MT libSISConePlugin_la-SISConeBasePlugin.lo -MD -MP -MF .deps/libSISConePlugin_la-SISConeBasePlugin.Tpo -c SISConeBasePlugin.cc -fno-common -DPIC -o .libs/libSISConePlugin_la-SISConeBasePlugin.o
> libtool: compile: g++ -DHAVE_CONFIG_H -I. -I../../include/fastjet -O2 -Wall -g -Woverloaded-virtual -DDROP_CGAL -I. -I./siscone -I./../../include -I./siscone -MT libSISConePlugin_la-SISConeSphericalPlugin.lo -MD -MP -MF .deps/libSISConePlugin_la-SISConeSphericalPlugin.Tpo -c SISConeSphericalPlugin.cc -fno-common -DPIC -o .libs/libSISConePlugin_la-SISConeSphericalPlugin.o
> SISConeBasePlugin.cc:12:12: error: no matching member function for call to 'structure_of'
> return a.structure_of<UserScaleBase::StructureType>().ordering_var2()
> ~~^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
> ./../../include/fastjet/PseudoJet.hh:1006:60: note: candidate template ignored: substitution failure [with TransformerType = fastjet::SISConeBasePlugin::UserScaleBase::StructureType]: ISO C++ specifies that qualified reference to 'StructureType' is a constructor name rather than a type in this context, despite preceding 'typename' keyword
> const typename TransformerType::StructureType & PseudoJet::structure_of() const{
> ~~~~~~~~ ^
> SISConeBasePlugin.cc:13:12: error: no matching member function for call to 'structure_of'
>> b.structure_of<UserScaleBase::StructureType>().ordering_var2();
> ~~^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
> ./../../include/fastjet/PseudoJet.hh:1006:60: note: candidate template ignored: substitution failure [with TransformerType = fastjet::SISConeBasePlugin::UserScaleBase::StructureType]: ISO C++ specifies that qualified reference to 'StructureType' is a constructor name rather than a type in this context, despite preceding 'typename' keyword
> const typename TransformerType::StructureType & PseudoJet::structure_of() const{
> ~~~~~~~~ ^
> 2 errors generated.
> make[3]: *** [libSISConePlugin_la-SISConeBasePlugin.lo] Error 1
> make[3]: *** Waiting for unfinished jobs....
> libtool: compile: g++ -DHAVE_CONFIG_H -I. -I../../include/fastjet -O2 -Wall -g -Woverloaded-virtual -DDROP_CGAL -I. -I./siscone -I./../../include -I./siscone -MT libSISConePlugin_la-SISConePlugin.lo -MD -MP -MF .deps/libSISConePlugin_la-SISConePlugin.Tpo -c SISConePlugin.cc -o libSISConePlugin_la-SISConePlugin.o >/dev/null 2>&1
> libtool: compile: g++ -DHAVE_CONFIG_H -I. -I../../include/fastjet -O2 -Wall -g -Woverloaded-virtual -DDROP_CGAL -I. -I./siscone -I./../../include -I./siscone -MT libSISConePlugin_la-SISConeSphericalPlugin.lo -MD -MP -MF .deps/libSISConePlugin_la-SISConeSphericalPlugin.Tpo -c SISConeSphericalPlugin.cc -o libSISConePlugin_la-SISConeSphericalPlugin.o >/dev/null 2>&1
> mv -f .deps/libSISConePlugin_la-SISConePlugin.Tpo .deps/libSISConePlugin_la-SISConePlugin.Plo
> mv -f .deps/libSISConePlugin_la-SISConeSphericalPlugin.Tpo .deps/libSISConePlugin_la-SISConeSphericalPlugin.Plo
> make[2]: *** [all-recursive] Error 1
> make[1]: *** [all-recursive] Error 1
> make: *** [all-recursive] Error 1
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.