Fork me on GitHub

Opened 3 years ago

Closed 3 years ago

#1507 closed Bug (fixed)

Compatible versions of MadGraph5 and PYTHIA8 that Delphes can read? ERROR: Unsupported block type

Reported by: Matthew Feickert Owned by:
Priority: minor Milestone:
Component: Delphes code Version: Delphes 3
Keywords: Unsupported block type, MadGraph5, PYTHIA Cc:

Description

Hi. Are there known upper limits on compatible versions of MadGraph5 and PYTHIA8 that Delphes can read?

With the following Docker image () I am able to run the example shown in the Delphes README

$ docker pull scailfin/delphes-python-centos:3.5.0
$ docker run --rm -ti -v $PWD:$PWD -w $PWD scailfin/delphes-python-centos:3.5.0
[root@8643ec787938 delphes_output]# curl -sLO http://cp3.irmp.ucl.ac.be/downloads/z_ee.hep.gz
[root@8643ec787938 delphes_output]# gunzip z_ee.hep.gz
[root@8643ec787938 delphes_output]# DelphesSTDHEP /usr/local/venv/cards/delphes_card_ATLAS.tcl delphes_output.root z_ee.hep
** INFO: adding module        ParticlePropagator       ParticlePropagator
** INFO: adding module        Efficiency               ChargedHadronTrackingEfficiency
** INFO: adding module        Efficiency               ElectronTrackingEfficiency
** INFO: adding module        Efficiency               MuonTrackingEfficiency
** INFO: adding module        MomentumSmearing         ChargedHadronMomentumSmearing
** INFO: adding module        MomentumSmearing         ElectronMomentumSmearing
** INFO: adding module        MomentumSmearing         MuonMomentumSmearing
** INFO: adding module        Merger                   TrackMerger
** INFO: adding module        SimpleCalorimeter        ECal
** INFO: adding module        SimpleCalorimeter        HCal
** INFO: adding module        PdgCodeFilter            ElectronFilter
** INFO: adding module        PdgCodeFilter            ChargedHadronFilter
** INFO: adding module        Merger                   Calorimeter
** INFO: adding module        Merger                   EFlowMerger
** INFO: adding module        PdgCodeFilter            EFlowFilter
** INFO: adding module        Efficiency               PhotonEfficiency
** INFO: adding module        Isolation                PhotonIsolation
** INFO: adding module        Efficiency               ElectronEfficiency
** INFO: adding module        Isolation                ElectronIsolation
** INFO: adding module        Efficiency               MuonEfficiency
** INFO: adding module        Isolation                MuonIsolation
** INFO: adding module        Merger                   MissingET
** INFO: adding module        Merger                   ScalarHT
** INFO: adding module        PdgCodeFilter            NeutrinoFilter
** INFO: adding module        FastJetFinder            GenJetFinder
** INFO: adding module        Merger                   GenMissingET
** INFO: adding module        FastJetFinder            FastJetFinder
** INFO: adding module        EnergyScale              JetEnergyScale
** INFO: adding module        JetFlavorAssociation     JetFlavorAssociation
** INFO: adding module        BTagging                 BTagging
** INFO: adding module        TrackCountingTauTagging  TauTagging
** INFO: adding module        UniqueObjectFinder       UniqueObjectFinder
** INFO: adding module        TreeWriter               TreeWriter
** INFO: initializing module  Delphes
** INFO: initializing module  ParticlePropagator
** INFO: initializing module  ChargedHadronTrackingEfficiency
** INFO: initializing module  ElectronTrackingEfficiency
** INFO: initializing module  MuonTrackingEfficiency
** INFO: initializing module  ChargedHadronMomentumSmearing
** INFO: initializing module  ElectronMomentumSmearing
** INFO: initializing module  MuonMomentumSmearing
** INFO: initializing module  TrackMerger
** INFO: initializing module  ECal
** INFO: initializing module  HCal
** INFO: initializing module  Calorimeter
** INFO: initializing module  EFlowMerger
** INFO: initializing module  EFlowFilter
** INFO: initializing module  PhotonEfficiency
** INFO: initializing module  PhotonIsolation
** INFO: initializing module  ElectronFilter
** INFO: initializing module  ElectronEfficiency
** INFO: initializing module  ElectronIsolation
** INFO: initializing module  ChargedHadronFilter
** INFO: initializing module  MuonEfficiency
** INFO: initializing module  MuonIsolation
** INFO: initializing module  MissingET
** INFO: initializing module  NeutrinoFilter
** INFO: initializing module  GenJetFinder
#--------------------------------------------------------------------------
#                         FastJet release 3.3.4
#                 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 GNU GPL v2 or higher.
# 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  GenMissingET
** INFO: initializing module  FastJetFinder
** INFO: initializing module  JetEnergyScale
** INFO: initializing module  JetFlavorAssociation
** INFO: initializing module  BTagging
** INFO: initializing module  TauTagging
** INFO: initializing module  UniqueObjectFinder
** INFO: initializing module  ScalarHT
** INFO: initializing module  TreeWriter
** Reading z_ee.hep
** [################################################################] (100.00%)
** Exiting...
[root@8643ec787938 delphes_output]# root -l delphes_output.root
root [0]
Attaching file delphes_output.root as _file0...
Warning in <TClass::Init>: no dictionary for class LHEFEvent is available
Warning in <TClass::Init>: no dictionary for class Event is available
Warning in <TClass::Init>: no dictionary for class GenParticle is available
Warning in <TClass::Init>: no dictionary for class SortableObject is available
Warning in <TClass::Init>: no dictionary for class Track is available
Warning in <TClass::Init>: no dictionary for class Tower is available
Warning in <TClass::Init>: no dictionary for class Jet is available
Warning in <TClass::Init>: no dictionary for class MissingET is available
Warning in <TClass::Init>: no dictionary for class Electron is available
Warning in <TClass::Init>: no dictionary for class Photon is available
Warning in <TClass::Init>: no dictionary for class Muon is available
Warning in <TClass::Init>: no dictionary for class ScalarHT is available
(TFile *) 0x3ccae20
root [1] _file0->ls()
TFile**		delphes_output.root
 TFile*		delphes_output.root
  KEY: TProcessID	ProcessID0;1	105910a6-d324-11eb-ab85-020011acbeef
  KEY: TTree	Delphes;1	Analysis tree

However, if I try to run on some events that I produced with MadGraph5 v3.1.1 and PYTHIA8 (https://github.com/scailfin/MadGraph5_aMC-NLO/blob/main/docker/centos/Dockerfile) I get an "ERROR: Unsupported block type."

$ docker run --rm -ti -v $PWD:$PWD -w $PWD scailfin/delphes-python-centos:3.5.0
[root@cfe9c43c2bda nevents_10e4_run]# DelphesSTDHEP /usr/local/venv/cards/delphes_card_ATLAS.tcl test_output.root drell-yan_output/Events/run_01/tag_1_pythia8_events.hepmc
** INFO: adding module        ParticlePropagator       ParticlePropagator
** INFO: adding module        Efficiency               ChargedHadronTrackingEfficiency
** INFO: adding module        Efficiency               ElectronTrackingEfficiency
** INFO: adding module        Efficiency               MuonTrackingEfficiency
** INFO: adding module        MomentumSmearing         ChargedHadronMomentumSmearing
** INFO: adding module        MomentumSmearing         ElectronMomentumSmearing
** INFO: adding module        MomentumSmearing         MuonMomentumSmearing
** INFO: adding module        Merger                   TrackMerger
** INFO: adding module        SimpleCalorimeter        ECal
** INFO: adding module        SimpleCalorimeter        HCal
** INFO: adding module        PdgCodeFilter            ElectronFilter
** INFO: adding module        PdgCodeFilter            ChargedHadronFilter
** INFO: adding module        Merger                   Calorimeter
** INFO: adding module        Merger                   EFlowMerger
** INFO: adding module        PdgCodeFilter            EFlowFilter
** INFO: adding module        Efficiency               PhotonEfficiency
** INFO: adding module        Isolation                PhotonIsolation
** INFO: adding module        Efficiency               ElectronEfficiency
** INFO: adding module        Isolation                ElectronIsolation
** INFO: adding module        Efficiency               MuonEfficiency
** INFO: adding module        Isolation                MuonIsolation
** INFO: adding module        Merger                   MissingET
** INFO: adding module        Merger                   ScalarHT
** INFO: adding module        PdgCodeFilter            NeutrinoFilter
** INFO: adding module        FastJetFinder            GenJetFinder
** INFO: adding module        Merger                   GenMissingET
** INFO: adding module        FastJetFinder            FastJetFinder
** INFO: adding module        EnergyScale              JetEnergyScale
** INFO: adding module        JetFlavorAssociation     JetFlavorAssociation
** INFO: adding module        BTagging                 BTagging
** INFO: adding module        TrackCountingTauTagging  TauTagging
** INFO: adding module        UniqueObjectFinder       UniqueObjectFinder
** INFO: adding module        TreeWriter               TreeWriter
** INFO: initializing module  Delphes
** INFO: initializing module  ParticlePropagator
** INFO: initializing module  ChargedHadronTrackingEfficiency
** INFO: initializing module  ElectronTrackingEfficiency
** INFO: initializing module  MuonTrackingEfficiency
** INFO: initializing module  ChargedHadronMomentumSmearing
** INFO: initializing module  ElectronMomentumSmearing
** INFO: initializing module  MuonMomentumSmearing
** INFO: initializing module  TrackMerger
** INFO: initializing module  ECal
** INFO: initializing module  HCal
** INFO: initializing module  Calorimeter
** INFO: initializing module  EFlowMerger
** INFO: initializing module  EFlowFilter
** INFO: initializing module  PhotonEfficiency
** INFO: initializing module  PhotonIsolation
** INFO: initializing module  ElectronFilter
** INFO: initializing module  ElectronEfficiency
** INFO: initializing module  ElectronIsolation
** INFO: initializing module  ChargedHadronFilter
** INFO: initializing module  MuonEfficiency
** INFO: initializing module  MuonIsolation
** INFO: initializing module  MissingET
** INFO: initializing module  NeutrinoFilter
** INFO: initializing module  GenJetFinder
#--------------------------------------------------------------------------
#                         FastJet release 3.3.4
#                 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 GNU GPL v2 or higher.
# 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  GenMissingET
** INFO: initializing module  FastJetFinder
** INFO: initializing module  JetEnergyScale
** INFO: initializing module  JetFlavorAssociation
** INFO: initializing module  BTagging
** INFO: initializing module  TauTagging
** INFO: initializing module  UniqueObjectFinder
** INFO: initializing module  ScalarHT
** INFO: initializing module  TreeWriter
** Reading drell-yan_output/Events/run_01/tag_1_pythia8_events.hepmc
** ERROR: Unsupported block type.

(I'm using these Docker images because I'm under annoying constraints that force me to have a very old version of glibc so I've had to use centos:7 as a base image to build on.)

Are there known bounds that I should avoid?

Change History (3)

comment:1 by Pavel Demin, 3 years ago

In the second example, the format of the tag_1_pythia8_events.hepmc file seems to be HepMC. So, the command to process this file with Delphes should be

DelphesHepMC2 /usr/local/venv/cards/delphes_card_ATLAS.tcl test_output.root drell-yan_output/Events/run_01/tag_1_pythia8_events.hepmc

comment:2 by Matthew Feickert, 3 years ago

Sorry about that --- you're obviously totally right. This is pretty embarrassing, but thanks for pointing out this error!

This

$ docker run --rm -ti -v $PWD:$PWD -w $PWD scailfin/delphes-python-centos:3.5.0
[root@aa8fd8aa3d98 drell-yan]# DelphesHepMC2 /usr/local/venv/cards/delphes_card_ATLAS.tcl test_output.root drell-yan_output/Events/run_01/tag_1_pythia8_events.hepmc

works now as expected.

Thanks again and this ticket can get closed.

comment:3 by Pavel Demin, 3 years ago

Resolution: fixed
Status: newclosed
Note: See TracTickets for help on using tickets.