LArSoft  v07_13_02
Liquid Argon Software toolkit - http://larsoft.org/
exampleB02.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: exampleB02.cc 105211 2017-07-14 12:49:11Z gcosmo $
31 //
32 //
33 // --------------------------------------------------------------
34 // GEANT 4 - exampleB02
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 //
48 
49 // --------------------------------------------------------------
50 
51 #include <iostream>
52 
53 #include <stdlib.h>
54 
55 #ifdef G4MULTITHREADED
56 #include "G4MTRunManager.hh"
57 #else
58 #include "G4RunManager.hh"
59 #endif
60 
61 #include "G4VPhysicalVolume.hh"
62 #include "G4UImanager.hh"
63 #include "G4GeometryManager.hh"
64 
65 // user classes
66 #include "B02DetectorConstruction.hh"
67 #include "B02ImportanceDetectorConstruction.hh"
68 #include "FTFP_BERT.hh"
69 #include "G4ImportanceBiasing.hh"
70 #include "G4ParallelWorldPhysics.hh"
71 
72 #include "B02ActionInitialization.hh"
73 // #include "B02PrimaryGeneratorAction.hh"
74 // #include "B02RunAction.hh"
75 
76 // Files specific for biasing and scoring
77 #include "G4GeometrySampler.hh"
78 #include "G4IStore.hh"
79 #include "G4VWeightWindowStore.hh"
80 #include "G4WeightWindowAlgorithm.hh"
81 
82 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
83 
84 int main(int argc, char **argv)
85 {
86  G4int mode = 0;
87  if (argc>1) mode = atoi(argv[1]);
88  if(mode != 0) {
89  G4cout << " mode not used yet, refer to B01 to see WeightWindow technique "
90  << mode << G4endl;
91  }
92 
93  G4int numberOfEvents = 100;
94  G4long myseed = 345354;
95 
96 #ifdef G4MULTITHREADED
97  G4MTRunManager * runManager = new G4MTRunManager;
98  G4cout << " Number of cores: " << G4Threading::G4GetNumberOfCores() << G4endl;
99  G4cout << " but using two! " << G4endl;
100  runManager->SetNumberOfThreads(2);
101  // runManager->SetNumberOfThreads(G4Threading::G4GetNumberOfCores());
102 #else
103  G4RunManager * runManager = new G4RunManager;
104 #endif
105 
106  G4Random::setTheSeed(myseed);
107 
108  // create the detector ---------------------------
109  B02DetectorConstruction* detector = new B02DetectorConstruction();
110  runManager->SetUserInitialization(detector);
111 
112  G4String parallelName("ParallelBiasingWorld");
113  B02ImportanceDetectorConstruction* pdet =
114  new B02ImportanceDetectorConstruction(parallelName);
115  detector->RegisterParallelWorld(pdet);
116 
117  G4GeometrySampler pgs(pdet->GetWorldVolume(),"neutron");
118 
119  pgs.SetParallel(true);
120 
121  G4VModularPhysicsList* physicsList = new FTFP_BERT;
122  physicsList->RegisterPhysics(new G4ImportanceBiasing(&pgs,parallelName));
123  physicsList->RegisterPhysics(new G4ParallelWorldPhysics(parallelName));
124 
125  runManager->SetUserInitialization(physicsList);
126 
127  // Set user action classes through Worker Initialization
128  //
129  B02ActionInitialization* actions = new B02ActionInitialization;
130  runManager->SetUserInitialization(actions);
131 
132  runManager->Initialize();
133 
134  pdet->CreateImportanceStore();
135 
136  //temporary fix before runManager->BeamOn works...
137  G4UImanager* UImanager = G4UImanager::GetUIpointer();
138  G4String command1 = "/control/cout/setCoutFile threadOut";
139  UImanager->ApplyCommand(command1);
140  G4String command2 = "/run/beamOn " +
141  G4UIcommand::ConvertToString(numberOfEvents);
142  UImanager->ApplyCommand(command2);
143 
144  // open geometry for clean biasing stores clean-up
145  //
146  G4GeometryManager::GetInstance()->OpenGeometry();
147 
148  // pgs.ClearSampling();
149 
150  delete runManager;
151 
152  return 0;
153 }
154 
155 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
int main(int argc, char **argv)
Definition: exampleB02.cc:84