--- /dev/null
+
+// Flugg tag
+
+///////////////////////////////////////////////////////////////////
+//
+// WrapInit.hh - Sara Vanini
+//
+// Wrapper for geometry initialisation.
+//
+// modified 12-IV-2000
+// modified 24.10.01: by I. Hrivnacova
+// functions declarations separated from implementation
+// (moved to Wrappers.hh);
+//
+//////////////////////////////////////////////////////////////////
+
+
+#include "Wrappers.hh"
+#include "FGeometryInit.hh"
+#include "G4PhysicalVolumeStore.hh"
+#include "G4VPhysicalVolume.hh"
+#include "globals.hh"
+
+void jomiwr(const G4int & nge, const G4int& lin, const G4int& lou,
+ G4int& flukaReg)
+{
+//flag
+#ifdef G4GEOMETRY_DEBUG
+ G4cout << "================== JOMIWR =================" << G4endl;
+#endif
+
+
+ //Geoinit Pointer
+#ifdef G4GEOMETRY_DEBUG
+ G4cout << "\t *==> JOMIWR: Getting FGeometry..." << G4endl;
+#endif
+ FGeometryInit * ptrGeoInit=FGeometryInit::GetInstance();
+
+ //initialize geometry:construct detector and set world volume
+#ifdef G4GEOMETRY_DEBUG
+ G4cout << "\t *==> JOMIWR: Setting the detector..." << G4endl;
+#endif
+ ptrGeoInit->setDetector();
+#ifdef G4GEOMETRY_DEBUG
+ G4cout << "\t *==> JOMIWR: Setting mother volume..." << G4endl;
+#endif
+ ptrGeoInit->setMotherVolume();
+
+ //close geometry for optimization
+#ifdef G4GEOMETRY_DEBUG
+ G4cout << "\t *==> JOMIWR: Closing geometry..." << G4endl;
+#endif
+ ptrGeoInit->closeGeometry();
+
+ //initialize wrappers utility histories at the beginning of run and set flag
+#ifdef G4GEOMETRY_DEBUG
+ G4cout << "\t *==> JOMIWR: InitHistories..." << G4endl;
+#endif
+ ptrGeoInit->InitHistories();
+
+ //initialize lattice array
+#ifdef G4GEOMETRY_DEBUG
+ G4cout << "\t *==> JOMIWR: Init lattice array..." << G4endl;
+#endif
+ ptrGeoInit->InitJrLtGeantArray();
+
+ //initialize debug-array
+#ifdef G4GEOMETRY_DEBUG
+ G4cout << "\t *==> JOMIWR: Init debug array..." << G4endl;
+#endif
+ ptrGeoInit->InitHistArray();
+
+ //create Fluka material cards in flukaMat.inp file
+#ifdef G4GEOMETRY_DEBUG
+ G4cout << "\t *==> JOMIWR: Init fluka materials..." << G4endl;
+#endif
+ ptrGeoInit->createFlukaMatFile();
+
+ //returns number of volumes + 1
+#ifdef G4GEOMETRY_DEBUG
+ G4cout << "\t *==> JOMIWR: Returning..." << G4endl;
+#endif
+ G4PhysicalVolumeStore * pVolStore = G4PhysicalVolumeStore::GetInstance();
+ G4int numVol = G4int(pVolStore->size());
+ flukaReg = numVol + 1;
+
+#ifdef G4GEOMETRY_DEBUG
+ G4cout << "Number of volumes + 1: " << flukaReg << G4endl;
+ G4cout << "================== Out of JOMIWR =================" << G4endl;
+#endif
+}
+
+
+