Add tests for the EventReader plugins

add currentEventNumber function to check which event one is at
Add small input files for the different event readers

DDG4/include/DDG4/Geant4InputAction.h

@@ -71,6 +71,8 @@ namespace DD4hep  {
       const std::string& name()  const   {  return m_name;   }
       /// Flag if direct event access (by event sequence number) is supported (Default: false)
       bool hasDirectAccess() const  {  return m_directAccess; }
+      /// return current Event Number
+      int currentEventNumber() const { return m_currEvent; }
       /// Move to the indicated event number.
       /** For pure sequential access, the default implementation
        *  will skip events one by one.

DDTest/CMakeLists.txt

@@ -27,3 +27,8 @@ dd4hep_add_test_reg ( test_bitfield64          BUILD_EXEC REGEX_FAIL "TEST_FAILE
 dd4hep_add_test_reg ( test_PolarGridRPhi2      BUILD_EXEC REGEX_FAIL "TEST_FAILED" )
 dd4hep_add_test_reg ( test_cellDimensions      BUILD_EXEC REGEX_FAIL "TEST_FAILED" )
 dd4hep_add_test_reg ( test_cellDimensionsRPhi2 BUILD_EXEC REGEX_FAIL "TEST_FAILED" )
+
+if (DD4HEP_USE_GEANT4)
+  dd4hep_add_test_reg ( test_EventReaders BUILD_EXEC REGEX_FAIL "TEST_FAILED"
+    EXEC_ARGS ${CMAKE_CURRENT_SOURCE_DIR} )
+endif()
\ No newline at end of file

DDTest/inputFiles/Muons10GeV.HEPEvt

+    1
+         1        13         0         0        0.00000         0.000000      10.0000000        0.105658
+    1
+         1        13         0         0        0.756007        0.688893       0.0423393        0.105658
+    1
+         1        13         0         0         5.51014         5.08509        -4.48918        0.105658
+    1
+         1        13         0         0        -2.30284         3.54976        0.851954        0.105658
+    1
+         1        13         0         0        0.842587         5.48815         4.16012        0.105658
+    1
+         1        13         0         0         1.23048         5.15334        -3.69239        0.105658
+    1
+         1        13         0         0         2.20157         3.34482         1.73308        0.105658
+    1
+         1        13         0         0         7.56195          3.2313        -2.75249        0.105658
+    1
+         1        13         0         0         2.05324         4.80771        -1.32544        0.105658
+    1
+         1        13         0         0        -3.82285         5.76216         6.65821        0.105658
+    1
+         1        13         0         0         1.74145         2.98833         1.27817        0.105658
+    1
+         1        13         0         0        -3.56344         7.32584         4.57862        0.105658
+    1
+         1        13         0         0        -1.30765         7.52284         2.52024        0.105658
+    1
+         1        13         0         0        -5.23125          2.6905        -3.23875        0.105658
+    1
+         1        13         0         0         7.58875         1.98954         4.39831        0.105658
+    1
+         1        13         0         0        -4.03736         2.50751      -0.0326271        0.105658
+    1
+         1        13         0         0        0.522977        0.507085        0.552055        0.105658
+    1
+         1        13         0         0        -5.78818         5.33183         3.62157        0.105658
+    1
+         1        13         0         0       -0.696683         8.01513        -1.87437        0.105658
+    1
+         1        13         0         0       -0.479692           4.988        -1.47847        0.105658
+    1
+         1        13         0         0        -4.35663         5.80079         1.74009        0.105658
+    1
+         1        13         0         0        -4.43793         2.14182        0.704618        0.105658
+    1
+         1        13         0         0         4.22506         6.46341        -1.03411        0.105658
+    1
+         1        13         0         0          1.3372        0.799973        -1.14639        0.105658
+    1
+         1        13         0         0         1.03267          1.5407         1.65847        0.105658

DDTest/src/test_EventReaders.cc

+#include "DD4hep/DDTest.h"
+
+#include <iostream>
+#include <iomanip>
+#include <vector>
+#include <algorithm>
+#include <exception>
+
+#include "DD4hep/Plugins.h"
+#include "DD4hep/Primitives.h"
+#include "DDG4/Geant4InputAction.h"
+#include "DDG4/Geant4Particle.h"
+
+//using namespace DD4hep::Simulation;
+typedef DD4hep::Simulation::Geant4Particle Particle;
+
+DD4hep::DDTest test = DD4hep::DDTest( "EventReader" ) ;
+
+class TestTuple {
+public:
+  std::string readerType;
+  std::string inputFile;
+  bool skipEOF; //LCIO skipNEvents does not signal end of file??
+  TestTuple( std::string const& rT, std::string const& iF, bool sEOF=false): readerType(rT), inputFile(iF), skipEOF(sEOF) {}
+};
+
+
+int main(int argc, char** argv ){
+
+  if( argc < 2 ) {
+    std::cout << " usage:  test_EventReaders Path/To/InputFiles " << std::endl ;
+    exit(1) ;
+  }
+
+  std::string inputFileFolder = argv[1];
+
+  std::vector<TestTuple> tests;
+  tests.push_back( TestTuple( "LCIOStdHepReader", "bbudsc_3evt.stdhep" ) );
+  tests.push_back( TestTuple( "LCIOFileReader",   "muons.slcio" , /*skipEOF= */ true ) );
+  tests.push_back( TestTuple( "Geant4EventReaderHepEvtShort", "Muons10GeV.HEPEvt" ) );
+  tests.push_back( TestTuple( "Geant4EventReaderHepMC", "g4pythia.hepmc" ) );
+
+
+  try{
+    for(std::vector<TestTuple>::const_iterator it = tests.begin(); it != tests.end(); ++it) {
+      std::string readerType = (*it).readerType;
+      std::string fileName = (*it).inputFile;
+      bool skipEOF =  (*it).skipEOF;
+      //InputFiles are in DDTest/inputFiles, argument is cmake_source directory
+      std::string inputFile = argv[1]+ std::string("/inputFiles/") + fileName;
+      DD4hep::Simulation::Geant4EventReader* thisReader = DD4hep::PluginService::Create<DD4hep::Simulation::Geant4EventReader*>(readerType, inputFile);
+      test( thisReader->currentEventNumber() == 0 , readerType + std::string("Initial Event Number") );
+      thisReader->moveToEvent(1);
+      test( thisReader->currentEventNumber() == 1 , readerType + std::string("Event Number after Skip") );
+      std::vector<Particle*> particles;
+      DD4hep::Simulation::Geant4EventReader::EventReaderStatus sc = thisReader->readParticles(3,particles);
+      std::for_each(particles.begin(),particles.end(),DD4hep::deleteObject<Particle>);
+      test( thisReader->currentEventNumber() == 2 && sc == DD4hep::Simulation::Geant4EventReader::EVENT_READER_OK,
+	    readerType + std::string("Event Number Read") );
+
+      //Reset Reader to check what happens if moving to far in the file
+      if (not skipEOF) {
+	thisReader = DD4hep::PluginService::Create<DD4hep::Simulation::Geant4EventReader*>(readerType, inputFile);
+	sc = thisReader->moveToEvent(1000000);
+	test( sc != DD4hep::Simulation::Geant4EventReader::EVENT_READER_OK , readerType + std::string("EventReader False") );
+      }
+    }
+
+  } catch( std::exception &e ){
+    //} catch( ... ){
+
+    test.log( e.what() );
+    test.error( "exception occurred" );
+  }
+  return 0;
+}
+