35 #include "RunAction.hh"
37 #include "TrackingAction.hh"
40 #include "Analysis.hh"
42 #include "CommandLineParser.hh"
44 using namespace G4DNAPARSER;
46 void PrintNParticles(std::map<const G4ParticleDefinition*, int>& container);
51 fpTrackingAction(0), fInitialized(0), fDebug(
false)
108 G4cout <<
"===================================" <<
G4endl;
110 G4cout <<
"================ RunAction::BeginMaster" <<
G4endl;
112 G4cout <<
"===================================" <<
G4endl;
122 G4cout <<
"===================================" <<
G4endl;
123 G4cout <<
"================ RunAction::BeginWorker" <<
G4endl;
125 G4cout <<
"===================================" <<
G4endl;
144 G4cout <<
"===================================" <<
G4endl;
145 G4cout <<
"================ RunAction::EndWorker" <<
G4endl;
147 G4cout <<
"===================================" <<
G4endl;
151 if ( nofEvents == 0 )
155 G4cout <<
"================ NO EVENTS TREATED IN THIS RUN ==> Exit"
174 std::map<const G4ParticleDefinition*, int>&
177 G4cout <<
"Number and type of particles created outside region \"Target\" :"
184 std::map<const G4ParticleDefinition*, int>&
187 G4cout <<
"Number and type of particles created in region \"Target\" :"
202 GetUserTrackingAction();
207 exDescrption <<
"fpTrackingAction is a null pointer. "
208 "Has it been correctly initialized ?";
231 G4cout <<
"##### Create analysis manager " <<
" " <<
this <<
G4endl;
235 " analysis manager" <<
G4endl;
245 if(command->GetOption().empty() ==
false)
247 fileName = command->GetOption();
251 fileName =
"microdosimetry";
253 analysisManager->
OpenFile(fileName);
257 analysisManager->
CreateNtuple(
"microdosimetry",
"physics");
290 analysisManager->
Write();
295 G4cout <<
"================ ROOT FILES HAVE BEEN WRITTEN"
304 G4cout <<
"================ Run is = "
306 G4cout <<
"================ Run type is = "
308 G4cout <<
"================ Event processed = "
310 G4cout <<
"================ Nevent = "
316 void PrintNParticles(std::map<const G4ParticleDefinition*, int>& container)
318 std::map<const G4ParticleDefinition*, int>::iterator
it;
319 for(it = container.begin() ;
320 it != container.end(); it ++)
322 G4cout <<
"N " << it->first->GetParticleName() <<
" : "