| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer 3 // * License and Disclaimer *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. 9 // * include a list of copyright holders. *
10 // * 10 // * *
11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. *
17 // * 17 // * *
18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. *
24 // ******************************************* 24 // ********************************************************************
25 // 25 //
26 /// \file biasing/B03/exampleB03.cc <<
27 /// \brief Main program of the biasing/B03 exa <<
28 // <<
29 // <<
30 // 26 //
>> 27 // $Id: exampleB03.cc,v 1.12 2006/06/29 16:34:59 gunter Exp $
>> 28 // GEANT4 tag $Name: geant4-08-03-patch-02 $
31 // 29 //
>> 30 //
32 // ------------------------------------------- 31 // --------------------------------------------------------------
33 // GEANT 4 - exampleB03 32 // GEANT 4 - exampleB03
34 // 33 //
35 // ------------------------------------------- 34 // --------------------------------------------------------------
36 // Comments 35 // Comments
>> 36 // This main may not be used at all. Instead you may use python or
>> 37 // lizard an execute the script B03RunApplication in them.
37 // 38 //
38 // This example intends to show how to use bot << 39 // This main function may used to test B03AppBase::GetB03AppBase().
39 // customized scoring making use of the scorin << 40 //
40 // in a parallel geometry. << 41 // --------------------------------------------------------------
41 // <<
42 // A simple geometry consisting of a 180 cm hi <<
43 // is constructed in the mass geometry. <<
44 // A geometry is constructed in the parallel g <<
45 // in order to assign importance values to sla <<
46 // of width 10cm and for scoring. The parallel <<
47 // overlap the mass world volume and the radi <<
48 // than the radius of the concrete cylinder in <<
49 // Pairs of G4GeometryCell and importance valu <<
50 // the importance store. <<
51 // The scoring uses the primitive scorers via <<
52 // <<
53 // <<
54 //....oooOO0OOooo........oooOO0OOooo........oo <<
55 <<
56 #include "B03ActionInitialization.hh" <<
57 #include "B03DetectorConstruction.hh" <<
58 #include "B03PhysicsList.hh" <<
59 <<
60 #include "G4RunManagerFactory.hh" <<
61 #include "G4SystemOfUnits.hh" <<
62 #include "G4Types.hh" <<
63 #include "G4UImanager.hh" <<
64 #include "G4VPhysicalVolume.hh" <<
65 <<
66 #include <iostream> <<
67 // #include "B03PrimaryGeneratorAction.hh" <<
68 // #include "B03RunAction.hh" <<
69 <<
70 // construction for the parallel geometry <<
71 #include "B03ImportanceDetectorConstruction.hh <<
72 <<
73 // Files specific for biasing and scoring <<
74 // #include "G4Scorer.hh" <<
75 // #include "G4GeometrySampler.hh" <<
76 #include "G4IStore.hh" <<
77 <<
78 //....oooOO0OOooo........oooOO0OOooo........oo <<
79 <<
80 int main(int, char**) <<
81 { <<
82 G4int numberOfEvents = 100; <<
83 <<
84 G4long myseed = 345354; <<
85 <<
86 auto* runManager = G4RunManagerFactory::Crea <<
87 runManager->SetNumberOfThreads(2); <<
88 <<
89 G4Random::setTheSeed(myseed); <<
90 <<
91 // create the detector ---------------- <<
92 B03DetectorConstruction* detector = new B03D <<
93 runManager->SetUserInitialization(detector); <<
94 // ---------------------------------------- <<
95 <<
96 // create a parallel detector <<
97 G4String parallelName("ParallelBiasingWorld" <<
98 B03ImportanceDetectorConstruction* pdet = ne <<
99 detector->RegisterParallelWorld(pdet); <<
100 <<
101 // G4GeometrySampler pgs(pdet->GetWorldVolum <<
102 // B03PhysicsList* physlist = new B03Physics <<
103 // name of first parallel world: <<
104 B03PhysicsList* physlist = new B03PhysicsLis <<
105 // push parallel world to store in case of m <<
106 physlist->AddParallelWorldName(parallelName) <<
107 // physlist->AddParallelWorldName(parallelNa <<
108 // physlist->AddParallelWorldName(sparallelN <<
109 // physlist->AddBiasing(&pgs,parallelName); <<
110 <<
111 runManager->SetUserInitialization(physlist); <<
112 <<
113 // Set user action classes through Worker In <<
114 // <<
115 B03ActionInitialization* actions = new B03Ac <<
116 runManager->SetUserInitialization(actions); <<
117 <<
118 // runManager->SetUserAction(new B03Primar <<
119 // // runManager->SetUserAction(new B03Pr <<
120 // runManager->SetUserAction(new B03RunAct <<
121 <<
122 runManager->Initialize(); <<
123 <<
124 G4VPhysicalVolume& aghostWorld = pdet->GetWo <<
125 G4cout << " ghost world: " << pdet->GetName( <<
126 <<
127 // create an importance <<
128 G4IStore* aIstore = G4IStore::GetInstance(pd <<
129 <<
130 // create a geometry cell for the world volu <<
131 G4GeometryCell gWorldVolumeCell(aghostWorld, <<
132 // set world volume importance to 1 <<
133 aIstore->AddImportanceGeometryCell(1, gWorld <<
134 <<
135 // set importance values and create scorers <<
136 G4int cell(1); <<
137 for (cell = 1; cell <= 18; cell++) { <<
138 G4GeometryCell gCell = pdet->GetGeometryCe <<
139 G4cout << " adding cell: " << cell << " re <<
140 << " name: " << gCell.GetPhysicalVo <<
141 G4double imp = std::pow(2.0, cell - 1); <<
142 // x aIstore.AddImportanceGeometryCell( <<
143 aIstore->AddImportanceGeometryCell(imp, gC <<
144 } <<
145 <<
146 // creating the geometry cell and add both t <<
147 // G4GeometryCell gCell = pdet->GetGeometry <<
148 <<
149 // create importance geometry cell pair for <<
150 // with the same importance as the last conc <<
151 G4GeometryCell gCell = pdet->GetGeometryCell <<
152 // G4double imp = std::pow(2.0,18); <<
153 G4double imp = std::pow(2.0, 17); <<
154 aIstore->AddImportanceGeometryCell(imp, gCel <<
155 <<
156 // temporary fix before runManager->BeamOn w <<
157 G4UImanager* UImanager = G4UImanager::GetUIp <<
158 G4String command1 = "/control/cout/setCoutFi <<
159 UImanager->ApplyCommand(command1); <<
160 <<
161 G4String command2 = "/run/beamOn " + G4UIcom <<
162 ; <<
163 UImanager->ApplyCommand(command2); <<
164 42
165 // runManager->BeamOn(numberOfEvents); << 43 #include "B03AppBase.hh"
>> 44 #include "globals.hh"
166 45
167 // pgs.ClearSampling(); << 46 int main(int argc, char **argv)
>> 47 {
168 48
169 delete runManager; << 49 B03AppBase &base = B03AppBase::GetB03AppBase();
170 50
171 return 0; 51 return 0;
172 } 52 }
173 53
174 //....oooOO0OOooo........oooOO0OOooo........oo <<
175 54