forked from cms-externals/parfullcms
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathParmainApplication.cc
126 lines (107 loc) · 4.28 KB
/
ParmainApplication.cc
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
#include "G4RunManager.hh"
#include "G4MTRunManager.hh"
#include "G4UImanager.hh"
#include "G4UIsession.hh"
#include "G4UIterminal.hh"
#include "G4UItcsh.hh"
#include "G4UIExecutive.hh"
#ifdef G4VIS_USE
#include "G4VisExecutive.hh"
#endif
#include "CLHEP/Random/MTwistEngine.h"
#include "CLHEP/Random/RanluxEngine.h"
#include "MyDetectorConstruction.hh"
#include "MyPrimaryGeneratorAction.hh"
#include "MyEventAction.hh"
#include "MyUserActionInitialization.hh"
// --- Physics Lists ---
#include "G4PhysListFactory.hh"
// ---
//01.25.2009 Xin Dong: This example came from the original sequential
//program FullCMS. The original program is changed here to support parallel
//computing with multiple threads. All events are assigned to each worker
//thread in a round robin fashion. All threads share most detector data
//including physics table and physics vector for some physics processes.
//The master process initializes the data in a regular way. However, worker
//threads initialize thread private data only.
//#include "G4MTParTopC.icc"
//01.25.2009 Xin Dong: Threads share this object.
//MyDetectorConstruction *detector = 0;
int main(int argc,char** argv) {
//CLHEP::RanluxEngine defaultEngine( 1234567, 4 );
CLHEP::MTwistEngine defaultEngine(1234567);
G4Random::setTheEngine( &defaultEngine );
G4long seed = time( NULL );
G4Random::setTheSeed( seed );
G4cout << G4endl
<< " ===================================================== " << G4endl
<< " Initial seed = " << seed << G4endl
<< " ===================================================== " << G4endl
<< G4endl;
#ifdef G4MULTITHREADED
G4MTRunManager* runManager = new G4MTRunManager; // = new G4RunManager;
//G4int nt = 2;
//if ( argc > 1 ) nt = atoi( argv[1] );
runManager->SetNumberOfThreads( 2 );
#else
G4RunManager* runManager = new G4RunManager;
#endif
//01.25.2009 Xin Dong: The master thread and worker threads have different behavior
//in the phase of initialization.
//if (threadRank == 0)
// runManager = new G4RunManager;
//else
// runManager = new G4RunManager(1);
//01.25.2009 Xin Dong: The master thread and worker threads have different behavior.
//in the phase of initialization.
//if (threadRank == 0)
MyDetectorConstruction* detector = new MyDetectorConstruction;
//else
//detector->SlaveMyDetectorConstruction();
runManager->SetUserInitialization(detector);
// --- Physics Lists ---
G4PhysListFactory factory;
runManager->SetUserInitialization( factory.ReferencePhysList() );
//thePL->SetDefaultCutValue( 0.020 *mm ); // 20 microns
runManager->SetUserInitialization( new MyUserActionInitialization );
#ifdef G4VIS_USE
// Visualization manager
//
G4VisManager* visManager = new G4VisExecutive;
visManager->Initialize();
#endif
runManager->Initialize();
// get the pointer to the UI manager and set verbosities
G4UImanager* UI = G4UImanager::GetUIpointer();
if ( argc==1 ) { // Define UI session for interactive mode.
G4UIExecutive* ui = new G4UIExecutive(argc,argv);
//G4UIsession * session = new G4UIterminal(new G4UItcsh);
UI->ApplyCommand("/control/execute vis.mac");
G4cout << "Now, please, apply beamOn command..." << G4endl;
ui->SessionStart();
//session->SessionStart();
delete ui;
} else { // Batch mode
G4String command = "/control/execute ";
G4String fileName = argv[1];
UI->ApplyCommand(command+fileName);
}
G4cout << G4endl
<< " ===================================================== " << G4endl
<< " Final random number = "
<< CLHEP::HepRandom::getTheEngine()->flat() << G4endl
<< " ===================================================== " << G4endl
<< G4endl;
#ifdef G4VIS_USE
delete visManager;
#endif
//01.25.2009 Xin Dong: The master thread and worker threads have different behavior
//in the phase of exit.
//if (threadRank != 0) detector->SlaveDestroy();
//01.25.2009 Xin Dong: The master thread and worker threads have different behavior
//in the phase of exit. It is better that all threads destroy this object. However,
//it results in panics by destroying shared data many times.
//if (threadRank == 0)
// delete runManager;
return 0;
}