| Geant4 Cross Reference |
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 //
26 // This example is provided by the Geant4-DNA collaboration
27 // Any report or published results obtained using the Geant4-DNA software
28 // shall cite the following Geant4-DNA collaboration publication:
29 // Med. Phys. 37 (2010) 4692-4708
30 // J. Comput. Phys. 274 (2014) 841-882
31 // The Geant4-DNA web site is available at http://geant4-dna.org
32 //
33 //
34 /// \file DetectorConstruction.cc
35 /// \brief Implementation of the DetectorConstruction class
36
37 #include "ActionInitialization.hh"
38
39 #include "DetectorConstruction.hh"
40 #include "PrimaryGeneratorAction.hh"
41 #include "RunAction.hh"
42 #include "SteppingAction.hh"
43 #include "TrackingAction.hh"
44
45 #include "G4DNAChemistryManager.hh"
46 #include "G4RunManager.hh"
47 #include "G4SystemOfUnits.hh"
48 #include "G4Threading.hh"
49 #include "G4UnitsTable.hh"
50
51 // chemistry
52 #include "ITSteppingAction.hh"
53 #include "ITTrackingAction.hh"
54 #include "ITTrackingInteractivity.hh"
55 #include "StackingAction.hh"
56 #include "TimeStepAction.hh"
57
58 #include "G4Scheduler.hh"
59
60 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
61
62 ActionInitialization::ActionInitialization() : G4VUserActionInitialization() {}
63
64 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
65
66 ActionInitialization::~ActionInitialization() {}
67
68 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
69
70 void ActionInitialization::BuildForMaster() const
71 {
72 // In MT mode, to be clearer, the RunAction class for the master thread might
73 // be different than the one used for the workers.
74 // This RunAction will be called before and after starting the
75 // workers.
76 // For more details, please refer to :
77 // https://twiki.cern.ch/twiki/bin/view/Geant4/Geant4MTForApplicationDevelopers
78 //
79 // RunAction* runAction= new RunAction();
80 // SetUserAction(runAction);
81 }
82
83 void ActionInitialization::Build() const
84 {
85 PrimaryGeneratorAction* primGenAction = new PrimaryGeneratorAction;
86 SetUserAction(primGenAction);
87
88 // Set optional user action classes
89 SetUserAction(new RunAction());
90 SetUserAction(new TrackingAction());
91 SetUserAction(new SteppingAction());
92 SetUserAction(new StackingAction());
93
94 // chemistry part
95 if (G4DNAChemistryManager::IsActivated()) {
96 G4Scheduler::Instance()->SetUserAction(new TimeStepAction());
97
98 // Uncomment and set to stop chemistry stage after:
99 // ...given number of time steps
100 // G4Scheduler::Instance()->SetMaxNbSteps(1000);
101
102 // ...OR reaching this time
103 G4Scheduler::Instance()->SetEndTime(100 * nanosecond);
104
105 G4Scheduler::Instance()->SetVerbose(1);
106
107 ITTrackingInteractivity* itInteractivity = new ITTrackingInteractivity();
108 itInteractivity->SetUserAction(new ITSteppingAction);
109 itInteractivity->SetUserAction(new ITTrackingAction);
110 G4Scheduler::Instance()->SetInteractivity(itInteractivity);
111 }
112 /*
113 // To output the pre-chemical stage
114 //
115 G4String fileName ("output");
116
117 if(G4RunManager::GetRunManager()->GetRunManagerType() ==
118 G4RunManager::sequentialRM)
119 {
120 // write initial situation at 1 picosecond
121 G4DNAChemistryManager::Instance()->WriteInto(fileName + ".txt");
122 }
123 else
124 {
125 G4int id = G4Threading::G4GetThreadId();
126
127 G4String fileName_mt = fileName;
128 fileName_mt += G4UIcommand::ConvertToString(id);
129 fileName_mt += ".txt";
130
131 G4cout << "chosen file name : " << fileName_mt << G4endl;
132
133 G4DNAChemistryManager::Instance()->WriteInto(fileName_mt);
134 }
135 */
136 }
137