]>
Commit | Line | Data |
---|---|---|
26911512 | 1 | |
2 | // Flugg tag | |
3 | ||
4 | /////////////////////////////////////////////////////////////////// | |
5 | // | |
6 | // WrapInit.hh - Sara Vanini | |
7 | // | |
8 | // Wrapper for geometry initialisation. | |
9 | // | |
10 | // modified 12-IV-2000 | |
11 | // modified 24.10.01: by I. Hrivnacova | |
12 | // functions declarations separated from implementation | |
13 | // (moved to Wrappers.hh); | |
14 | // | |
15 | ////////////////////////////////////////////////////////////////// | |
16 | ||
17 | ||
18 | #include "Wrappers.hh" | |
19 | #include "FGeometryInit.hh" | |
20 | #include "G4PhysicalVolumeStore.hh" | |
21 | #include "G4VPhysicalVolume.hh" | |
22 | #include "globals.hh" | |
23 | ||
24 | void jomiwr(const G4int & nge, const G4int& lin, const G4int& lou, | |
25 | G4int& flukaReg) | |
26 | { | |
27 | //flag | |
28 | #ifdef G4GEOMETRY_DEBUG | |
29 | G4cout << "================== JOMIWR =================" << G4endl; | |
30 | #endif | |
31 | ||
32 | ||
33 | //Geoinit Pointer | |
34 | #ifdef G4GEOMETRY_DEBUG | |
35 | G4cout << "\t *==> JOMIWR: Getting FGeometry..." << G4endl; | |
36 | #endif | |
37 | FGeometryInit * ptrGeoInit=FGeometryInit::GetInstance(); | |
38 | ||
39 | //initialize geometry:construct detector and set world volume | |
40 | #ifdef G4GEOMETRY_DEBUG | |
41 | G4cout << "\t *==> JOMIWR: Setting the detector..." << G4endl; | |
42 | #endif | |
7a9edf9c | 43 | // ptrGeoInit->setDetector(); |
26911512 | 44 | #ifdef G4GEOMETRY_DEBUG |
45 | G4cout << "\t *==> JOMIWR: Setting mother volume..." << G4endl; | |
46 | #endif | |
7a9edf9c | 47 | // ptrGeoInit->setMotherVolume(); |
26911512 | 48 | |
49 | //close geometry for optimization | |
50 | #ifdef G4GEOMETRY_DEBUG | |
51 | G4cout << "\t *==> JOMIWR: Closing geometry..." << G4endl; | |
52 | #endif | |
7a9edf9c | 53 | // ptrGeoInit->closeGeometry(); |
26911512 | 54 | |
55 | //initialize wrappers utility histories at the beginning of run and set flag | |
56 | #ifdef G4GEOMETRY_DEBUG | |
57 | G4cout << "\t *==> JOMIWR: InitHistories..." << G4endl; | |
58 | #endif | |
7a9edf9c | 59 | // ptrGeoInit->InitHistories(); |
26911512 | 60 | |
61 | //initialize lattice array | |
62 | #ifdef G4GEOMETRY_DEBUG | |
63 | G4cout << "\t *==> JOMIWR: Init lattice array..." << G4endl; | |
64 | #endif | |
7a9edf9c | 65 | // ptrGeoInit->InitJrLtGeantArray(); |
26911512 | 66 | |
67 | //initialize debug-array | |
68 | #ifdef G4GEOMETRY_DEBUG | |
69 | G4cout << "\t *==> JOMIWR: Init debug array..." << G4endl; | |
70 | #endif | |
7a9edf9c | 71 | // ptrGeoInit->InitHistArray(); |
26911512 | 72 | |
73 | //create Fluka material cards in flukaMat.inp file | |
74 | #ifdef G4GEOMETRY_DEBUG | |
75 | G4cout << "\t *==> JOMIWR: Init fluka materials..." << G4endl; | |
76 | #endif | |
7a9edf9c | 77 | // ptrGeoInit->createFlukaMatFile(); |
26911512 | 78 | |
79 | //returns number of volumes + 1 | |
80 | #ifdef G4GEOMETRY_DEBUG | |
81 | G4cout << "\t *==> JOMIWR: Returning..." << G4endl; | |
82 | #endif | |
83 | G4PhysicalVolumeStore * pVolStore = G4PhysicalVolumeStore::GetInstance(); | |
84 | G4int numVol = G4int(pVolStore->size()); | |
85 | flukaReg = numVol + 1; | |
86 | ||
87 | #ifdef G4GEOMETRY_DEBUG | |
88 | G4cout << "Number of volumes + 1: " << flukaReg << G4endl; | |
89 | G4cout << "================== Out of JOMIWR =================" << G4endl; | |
90 | #endif | |
91 | } | |
92 | ||
93 | ||
94 |