Revert "Merge branch 'cherry-pick-36c54a3b [formerly 5c0031d52e]-21.3' into '21.3'"

This reverts commit ccc04df0, reversing
changes made to a10ac508.

Simulation/ISF/ISF_Fatras/ISF_FatrasToolsG4/ISF_FatrasToolsG4/G4HadIntProcessor.h

@@ -130,25 +130,20 @@ namespace iFatras {
       /** describe deflection parametric/do real deflection */
       bool                                  m_parametricScattering;
 
-      /*
-       * Geant4 engine
-       * - Used only for `G4RunManager::SetVerboseLevel()` -> mutable
-       * - Does not take ownership -> no smart pointer
-      */
-      mutable G4RunManager* m_g4runManager;
-
-      /*
-       * Geant4 processes <PDGcode, process>  TODO : fission, capture
-       * - Used only for "easy access" to processes
-       * - Does not take ownership -> no smart pointers
-      */
+      //!< Geant4 engine
+      mutable G4RunManager*                m_g4runManager;
+      G4VUserPhysicsList*                  m_g4physicsList;
+      G4LayerDetectorConstruction*         m_g4detector;
+
+      //!< Geant4 processes <PDGcode, process>  TODO : fission, capture
       mutable std::map<int, G4VProcess*>   m_g4HadrInelasticProcesses;
       mutable std::map<int, G4VProcess*>   m_g4HadrElasticProcesses;
 
       //!< locally stored Geant4 instances (speeds up processing)
-      mutable std::unique_ptr<G4ThreeVector> m_g4zeroPos;
-      mutable std::unique_ptr<G4Step>        m_g4step;
-      mutable std::unique_ptr<G4StepPoint>   m_g4stepPoint;
+      mutable G4DynamicParticle*           m_g4dynPar;
+      mutable const G4ThreeVector*         m_g4zeroPos;
+      mutable G4Step*                      m_g4step;
+      mutable G4StepPoint*                 m_g4stepPoint;
       mutable std::vector<std::pair<float,std::pair< G4Material*, G4MaterialCutsCouple*> > > m_g4Material;
 
       /** ISF services & Tools */
@@ -158,11 +153,7 @@ namespace iFatras {
       /** MCTruth process code for TruthIncidents created by this tool */
       Barcode::PhysicsProcessCode          m_processCode;
 
-      /*
-       * Random engine
-       * - Initialised by `ServiceHandle<IAtRndmGenSvc>::GetEngine()`
-       * - Does not take ownership -> no smart pointer
-      */
+      /** Random engine  */
       CLHEP::HepRandomEngine*              m_randomEngine;
       std::string                          m_randomEngineName;       //!< Name of the random number stream
 

Simulation/ISF/ISF_Fatras/ISF_FatrasToolsG4/src/G4HadIntProcessor.cxx

@@ -89,10 +89,13 @@ iFatras::G4HadIntProcessor::G4HadIntProcessor(const std::string& t, const std::s
   m_minMomentum(50.0),
   m_cloneParameters(false),
   m_parametricScattering(false),
-  m_g4runManager(nullptr),
-  m_g4zeroPos(nullptr),
-  m_g4step(nullptr),
-  m_g4stepPoint(nullptr),
+  m_g4runManager(0),
+  m_g4physicsList(0),
+  m_g4detector(0),
+  m_g4dynPar(0),
+  m_g4zeroPos(0),
+  m_g4step(0),
+  m_g4stepPoint(0),
   m_particleBroker("ISF_ParticleBrokerSvc", n),
   m_truthRecordSvc("ISF_ValidationTruthService", n), 
   m_processCode(121),
@@ -204,6 +207,12 @@ StatusCode iFatras::G4HadIntProcessor::initialize()
 StatusCode iFatras::G4HadIntProcessor::finalize()
 {
 
+  // clean up locally stored Geant4 instances
+  delete m_g4dynPar;
+  delete m_g4zeroPos;
+  //delete m_g4step; will be deleted via m_g4stepPoint
+  delete m_g4stepPoint;
+
   ATH_MSG_INFO( " ---------- Statistics output -------------------------- " );
   //ATH_MSG_INFO( "                     Minimum energy cut for brem photons : " <<   m_minimumBremPhotonMomentum  );
   //ATH_MSG_INFO( "                     Brem photons (above cut, recorded)  : " <<   m_recordedBremPhotons        );
@@ -337,10 +346,11 @@ bool iFatras::G4HadIntProcessor::initG4RunManager() const {
   initProcessPDG(-321);
 
   // set up locally stored Geant4 instances
-  m_g4zeroPos.reset(new G4ThreeVector(0, 0, 0)),
-  m_g4step.reset(new G4Step()),
-  m_g4stepPoint.reset(new G4StepPoint()),
-  m_g4step->SetPreStepPoint(m_g4stepPoint.get());
+  m_g4dynPar          = new G4DynamicParticle();
+  m_g4zeroPos         = new G4ThreeVector(0, 0, 0);
+  m_g4step            = new G4Step();
+  m_g4stepPoint       = new G4StepPoint();
+  m_g4step->SetPreStepPoint( m_g4stepPoint);
 
   // define the available G4Material
   m_g4Material.clear();
@@ -446,7 +456,9 @@ ISF::ISFParticleVector iFatras::G4HadIntProcessor::getHadState(const ISF::ISFPar
   const G4ThreeVector mom( momentum.x(), momentum.y(), momentum.z() );
   inputPar->SetMomentum( mom);
   // position and timing dummy
-  std::unique_ptr<G4Track> g4track = std::make_unique<G4Track>( inputPar, 0 /* time */, *m_g4zeroPos);
+  G4Track* g4track=new G4Track( inputPar, 0 /* time */, *m_g4zeroPos);
+  //G4TouchableHandle g4touchable(new G4TouchableHistory());     // TODO check memory handling here
+  //g4track->SetTouchableHandle( g4touchable);
 
   // setup up G4Material ---------------------------------------------------------------------------
   std::pair<G4Material*, G4MaterialCutsCouple*> g4mat = retrieveG4Material(ematprop);
@@ -458,10 +470,9 @@ ISF::ISFParticleVector iFatras::G4HadIntProcessor::getHadState(const ISF::ISFPar
 
   // preparing G4Step and G4Track
   m_g4step->DeleteSecondaryVector();
-  g4track->SetStep(m_g4step.get());
+  g4track->SetStep( m_g4step);
 
   // by default, the current process is the inelastic hadr. interaction
-  // Does not take ownership (see definition of m_g4HadrInelasticProcesses) -> no delete or smart pointer
   G4VProcess *process = processIter_inelast!=m_g4HadrInelasticProcesses.end() ? processIter_inelast->second : 0;
  
   // if elastic interactions are enabled and there is a elastic process
@@ -473,18 +484,16 @@ ISF::ISFParticleVector iFatras::G4HadIntProcessor::getHadState(const ISF::ISFPar
     if( rand < 0.5) process = processIter_elast->second;
   }
 
-  if (!process)
-  {
-    ATH_MSG_WARNING( " [ ---- ] Cannot get G4VProcess from map (nullptr)!");
-    return chDef;
-  }
-
   ATH_MSG_VERBOSE ( " [ g4sim ] Computing " << process->GetProcessName() << " process with current particle" );
 
-  // Does not take ownership -> no delete or smart pointer
+  // do the G4VProcess (actually a G4HadronicProcess) ------------------------------------
+  //process->SetVerboseLevel(10);
+  //ATH_MSG_VERBOSE ( "Verbose Level is " << process->GetVerboseLevel() ); 
+  
   G4VParticleChange* g4change = process->PostStepDoIt(*g4track, *m_g4step);
   if (!g4change) {
     ATH_MSG_WARNING( " [ ---- ] Geant4 did not return any hadronic interaction information of particle with pdg=" << pdg );
+    delete g4track;
     return chDef;
   }
 
@@ -566,12 +575,14 @@ ISF::ISFParticleVector iFatras::G4HadIntProcessor::getHadState(const ISF::ISFPar
 
     // free up memory
     g4change->Clear();
+    delete g4track;
     return children;
 
   }
 
   // free up memory
   g4change->Clear();
+  delete g4track;
   return chDef;
 }