Geant4  10.02
exampleB01.cc
Go to the documentation of this file.
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
28 //
29 //
30 // $Id: exampleB01.cc 77475 2013-11-25 09:38:51Z gcosmo $
31 //
32 //
33 // --------------------------------------------------------------
34 // GEANT 4 - exampleB01
35 //
36 // --------------------------------------------------------------
37 // Comments
38 //
39 // This example intends to show how to use importance sampling and scoring
40 // in the mass (tracking) geometry.
41 // A simple geometry consisting of a 180 cm high concrete cylinder
42 // divided into 18 slabs of 10cm each is created.
43 // Importance values are assigned to the 18 concrete slabs in the
44 // detector construction class for simplicity.
45 // Pairs of G4GeometryCell and importance values are stored in
46 // the importance store.
47 // Scoring is carried out by the multifunctional detector (MFD) and
48 // sensitive detectors
49 //
50 // Alex Howard (alexander.howard@cern.ch):
51 // 22/11/13: Migrated to the new MT compliant design which moves the
52 // biasing process to the physicslist constructor - here
53 // via the modular physicslists
54 //
55 
56 // --------------------------------------------------------------
57 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
58 
59 #include <iostream>
60 
61 #include <stdlib.h>
62 
63 #ifdef G4MULTITHREADED
64 #include "G4MTRunManager.hh"
65 #else
66 #include "G4RunManager.hh"
67 #endif
68 
69 #include "G4VPhysicalVolume.hh"
70 #include "G4UImanager.hh"
71 #include "G4GeometryManager.hh"
72 
73 // user classes
75 #include "FTFP_BERT.hh"
76 #include "G4ImportanceBiasing.hh"
77 #include "G4WeightWindowBiasing.hh"
78 
80 // #include "B01PrimaryGeneratorAction.hh"
81 // #include "B01RunAction.hh"
82 
83 // Files specific for biasing and scoring
84 #include "G4GeometrySampler.hh"
85 #include "G4IStore.hh"
86 #include "G4VWeightWindowStore.hh"
88 
89 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
90 
91 int main(int argc, char **argv)
92 {
93  G4int mode = 0;
94  if (argc>1) mode = atoi(argv[1]);
95 
96  G4int numberOfEvents = 100;
97  G4long myseed = 345354;
98 
99 #ifdef G4MULTITHREADED
100  G4MTRunManager * runManager = new G4MTRunManager;
101  G4cout << " Number of cores: " << G4Threading::G4GetNumberOfCores() << G4endl;
102  G4cout << " and using 2 of them! " << G4endl;
103  runManager->SetNumberOfThreads(2);
104  // runManager->SetNumberOfThreads(G4Threading::G4GetNumberOfCores());
105 #else
106  G4RunManager * runManager = new G4RunManager;
107 #endif
108 
109  G4Random::setTheSeed(myseed);
110 
111  G4VWeightWindowAlgorithm *wwAlg = 0; // pointer for WeightWindow (mode>0)
112 
113  // create the detector ---------------------------
115  runManager->SetUserInitialization(detector);
116  G4GeometrySampler mgs(detector->GetWorldVolume(),"neutron");
117 
118  G4VModularPhysicsList* physicsList = new FTFP_BERT;
119  if(mode == 0)
120  {
121  physicsList->RegisterPhysics(new G4ImportanceBiasing(&mgs));
122  }
123  else
124  {
125  wwAlg = new G4WeightWindowAlgorithm(1, // upper limit factor
126  1, // survival factor
127  100); // max. number of splitting
128 
129  physicsList->RegisterPhysics(new G4WeightWindowBiasing
130  (&mgs, wwAlg, onBoundary));
131  // place of action
132  }
133  runManager->SetUserInitialization(physicsList);
134 
135  // Set user action classes through Worker Initialization
136  //
138  runManager->SetUserInitialization(actions);
139 
140  runManager->Initialize();
141 
142  if (mode == 0)
143  {
144  detector->CreateImportanceStore();
145  }
146  else
147  {
148  detector->CreateWeightWindowStore();
149  }
150 
151  // runManager->BeamOn(numberOfEvents);
152 
153  //temporary fix before runManager->BeamOn works...
154  G4UImanager* UImanager = G4UImanager::GetUIpointer();
155  G4String command1 = "/control/cout/setCoutFile threadOut";
156  UImanager->ApplyCommand(command1);
157  G4String command2 = "/run/beamOn " +
158  G4UIcommand::ConvertToString(numberOfEvents);;
159  UImanager->ApplyCommand(command2);
160 
161  // open geometry for clean biasing stores clean-up
162  //
164 
165  if (wwAlg) {
166  delete wwAlg;
167  }
168 
169  mgs.ClearSampling();
170 
171  delete runManager;
172 
173  return 0;
174 }
175 
176 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
virtual void SetUserInitialization(G4VUserDetectorConstruction *userInit)
int main(int argc, char **argv)
Definition: exampleB01.cc:91
G4VPhysicalVolume * GetWorldVolume()
void SetNumberOfThreads(G4int n)
Definition of the B01ActionInitialization class.
long G4long
Definition: G4Types.hh:80
static G4String ConvertToString(G4bool boolVal)
Definition: G4UIcommand.cc:371
Definition of the B01DetectorConstruction class.
int G4int
Definition: G4Types.hh:78
static G4UImanager * GetUIpointer()
Definition: G4UImanager.cc:58
G4GLOB_DLL std::ostream G4cout
G4VWeightWindowStore * CreateWeightWindowStore()
G4int G4GetNumberOfCores()
Definition: G4Threading.cc:127
static G4GeometryManager * GetInstance()
TFTFP_BERT< G4VModularPhysicsList > FTFP_BERT
Definition: FTFP_BERT.hh:63
virtual void Initialize()
#define G4endl
Definition: G4ios.hh:61
void OpenGeometry(G4VPhysicalVolume *vol=0)
G4int ApplyCommand(const char *aCommand)
Definition: G4UImanager.cc:446