4 // modified 10/IX/99 for including preStepPoint
5 // modified 28/IX/99 for delating allocated memory
6 // modified 4/X/99 function FreeMemory
7 // modified 2/III/00 base class G4TransportationManager included
8 // modified 20/III/00 PrintHistories() included
11 #ifndef FGeometryInit_h
12 #define FGeometryInit_h 1
14 //#include "g4std/fstream"
15 //#include "g4std/iomanip"
16 //#include "g4rw/cstring.h"
19 #include "G4LogicalVolume.hh"
20 #include "G4PhysicalVolumeStore.hh"
21 #include "G4VPhysicalVolume.hh"
22 #include "G4Material.hh"
23 #include "G4MaterialTable.hh"
24 #include "G4Isotope.hh"
25 #include "G4VUserDetectorConstruction.hh"
26 #include "G4TouchableHistory.hh"
27 #include "G4GeometryManager.hh"
28 #include "G4FieldManager.hh"
29 #include "G4UniformMagField.hh"
30 #include "G4TransportationManager.hh"
39 class FGeometryInit : public G4TransportationManager {
42 ~FGeometryInit(); //destructor
43 static FGeometryInit *GetInstance();
44 inline FluggNavigator *getNavigatorForTracking();
45 inline G4FieldManager * getFieldManager();
46 inline void setDetConstruction(G4VUserDetectorConstruction* detector);
47 inline void setDetector();
48 inline void setMotherVolume();
49 void createFlukaMatFile();
52 void PrintHistories();
54 void DeleteHistories();
55 void UpdateHistories(const G4NavigationHistory *, G4int);
56 inline G4TouchableHistory * GetTouchableHistory();
57 inline G4TouchableHistory * GetOldNavHist();
58 inline G4TouchableHistory * GetTempNavHist();
61 inline void DelHistArray();
62 inline G4int * GetHistArray();
64 void InitJrLtGeantArray();
65 inline G4int * GetJrLtGeantArray();
66 inline G4int GetLttcFlagGeant();
67 void SetLttcFlagGeant(G4int);
68 void PrintJrLtGeant();
73 void SetMediumFromName(const char* volName, int med, int volid);
74 //G4int GetRegionFromName(const char* volName) const;
75 G4int GetMedium(int) const;
76 int CurrentVolID(int ir, int& copyNo);
77 int CurrentVolOffID(int ir, int off, int& copyNo);
78 void Gmtod(double* xm, double* xd, int iflag);
79 void Gdtom(double* xd, double* xm, int iflag);
82 void BuildRegionsMap();
83 void PrintRegionsMap(std::ostream& os);
84 void BuildMaterialTables();
85 FlukaMaterial* BuildFlukaMaterialFromElement(const G4Element* element,
87 FlukaMaterial* BuildFlukaMaterialFromIsotope(const G4Isotope* isotope,
89 FlukaCompound* BuildFlukaCompoundFromMaterial(const G4Material* material);
90 FlukaCompound* BuildFlukaCompoundFromElement(const G4Element* element,
92 void PrintMaterialTables(std::ostream& os);
93 void PrintAssignmat(std::ostream& os);
94 void PrintMagneticField(std::ostream& os);
97 FGeometryInit(); //costructor
100 G4VUserDetectorConstruction * fDetector;
101 G4FieldManager * fFieldManager;
102 G4TransportationManager * fTransportationManager;
103 static FGeometryInit *flagInstance;
104 G4VPhysicalVolume * myTopNode;
105 G4GeometryManager * ptrGeoMan;
107 G4TouchableHistory * ptrTouchHist;
108 G4TouchableHistory * ptrOldNavHist;
109 G4TouchableHistory * ptrTempNavHist;
110 G4int * ptrJrLtGeant;
113 int* fRegionMediumMap;
115 std::map<G4VPhysicalVolume*, int, std::less<G4VPhysicalVolume*> > fRegionVolumeMap;
116 std::map<G4VPhysicalVolume*, int, std::less<G4VPhysicalVolume*> > fMediumVolumeMap;
117 std::map<G4VPhysicalVolume*, int, std::less<G4VPhysicalVolume*> > fVolIdVolumeMap;
119 std::map<G4Material*, FlukaMaterial*, std::less<G4Material*> > G4FlukaMaterialMap;
120 std::map<G4Material*, FlukaCompound*, std::less<G4Material*> > G4FlukaCompoundMap;
121 //G4int NOfMaterials;
124 typedef std::map<G4VPhysicalVolume*, int, std::less<G4VPhysicalVolume*> >::const_iterator RegionIterator;
125 typedef std::vector<G4Material*>::const_iterator MatTableIterator;
128 //Include the file with the inline methods
129 #include "FGeometryInit.icc"