Error in <TObjArray::At>: index -32649 out of bounds

[Sergei Chekanov]

Hello,

I have very busy events with 140 pile up events at 14 TeV (using HEPMC). Delphes3 breaks on such events with the attached error.

I should say it works fine when I have 50 and 100 overplayed pileup events, it start to give problems when the number of overplayed events more than 140.

best, Sergei

** 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        EnergySmearing           ElectronEnergySmearing   
** INFO: adding module        MomentumSmearing         MuonMomentumSmearing     
** INFO: adding module        Merger                   TrackMerger              
** INFO: adding module        Calorimeter              Calorimeter              
** INFO: adding module        Merger                   EFlowMerger              
** 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        FastJetFinder            FastJetFinder            
** INFO: adding module        BTagging                 BTagging                 
** INFO: adding module        TauTagging               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  ElectronEnergySmearing   
** INFO: initializing module  MuonMomentumSmearing     
** INFO: initializing module  TrackMerger              
** INFO: initializing module  Calorimeter              
** INFO: initializing module  EFlowMerger              
** INFO: initializing module  PhotonEfficiency         
** INFO: initializing module  PhotonIsolation          
** INFO: initializing module  ElectronEfficiency       
** INFO: initializing module  ElectronIsolation        
** INFO: initializing module  MuonEfficiency           
** INFO: initializing module  MuonIsolation            
** INFO: initializing module  MissingET                
** INFO: initializing module  FastJetFinder            
#--------------------------------------------------------------------------
#                         FastJet release 3.0.3
#                 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 LHC_signal_pileup.hepmc
Error in <TObjArray::At>: index -32649 out of bounds (size: 65536, this: 0x1cea2ab0)

 *** Break *** segmentation violation



===========================================================
There was a crash.
This is the entire stack trace of all threads:
===========================================================
#0  0x00000031772997b5 in waitpid () from /lib64/libc.so.6
#1  0x000000317723c761 in do_system () from /lib64/libc.so.6
#2  0x00002ae27c792722 in TUnixSystem::StackTrace() ()
   from /share/sl5/root_v5.34.Linux-sl5-native/lib/libCore.so
#3  0x00002ae27c78f4aa in TUnixSystem::DispatchSignals(ESignals) ()
   from /share/sl5/root_v5.34.Linux-sl5-native/lib/libCore.so
#4  <signal handler called>
#5  0x00002ae27c6d32d7 in TObject::TObject(TObject const&) ()
   from /share/sl5/root_v5.34.Linux-sl5-native/lib/libCore.so
#6  0x00002ae27ff28529 in TVector3::TVector3(TVector3 const&) ()
   from /share/sl5/root_v5.34.Linux-sl5-native/lib/libPhysics.so
#7  0x00000000004ce9d1 in Calorimeter::Process() ()
#8  0x00002ae27c716f99 in TTask::ExecuteTasks(char const*) ()
   from /share/sl5/root_v5.34.Linux-sl5-native/lib/libCore.so
#9  0x00002ae27c7169eb in TTask::ExecuteTask(char const*) ()
   from /share/sl5/root_v5.34.Linux-sl5-native/lib/libCore.so
#10 0x000000000058d25f in main ()
===========================================================


The lines below might hint at the cause of the crash.
If they do not help you then please submit a bug report at
http://root.cern.ch/bugs. Please post the ENTIRE stack trace
from above as an attachment in addition to anything else
that might help us fixing this issue.
===========================================================
#5  0x00002ae27c6d32d7 in TObject::TObject(TObject const&) ()
   from /share/sl5/root_v5.34.Linux-sl5-native/lib/libCore.so
#6  0x00002ae27ff28529 in TVector3::TVector3(TVector3 const&) ()
   from /share/sl5/root_v5.34.Linux-sl5-native/lib/libPhysics.so
#7  0x00000000004ce9d1 in Calorimeter::Process() ()
#8  0x00002ae27c716f99 in TTask::ExecuteTasks(char const*) ()
   from /share/sl5/root_v5.34.Linux-sl5-native/lib/libCore.so
#9  0x00002ae27c7169eb in TTask::ExecuteTask(char const*) ()
   from /share/sl5/root_v5.34.Linux-sl5-native/lib/libCore.so
#10 0x000000000058d25f in main ()
===========================================================

[Pavel Demin]

Looks like an index overflow in the Calorimeter module. I wrongly assumed that there won't be more than 30k particles in an event.

Could you, please send me a link to one of your HEPMC files? It'll help me to fix this problem.

[Pavel Demin]

I've changed the type of the particle index in modules/Calorimeter.c from Short_t to Int_t.

Could you, please check if the new version works better?

​​http://cp3.irmp.ucl.ac.be/downloads/Delphes-3.0.4.tar.gz

[Pavel Demin]

Delphes 3.0.4 is released now with some additional changes.