-
Notifications
You must be signed in to change notification settings - Fork 0
/
ICESPICE.cc
67 lines (54 loc) · 1.77 KB
/
ICESPICE.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
#include "G4Types.hh"
#include "G4RunManagerFactory.hh"
#include "G4UImanager.hh"
#include "Randomize.hh"
#include "G4VisExecutive.hh"
#include "G4UIExecutive.hh"
#include "G4AnalysisManager.hh"
#include "ICESPICEDetectorConstruction.hh"
#include "ICESPICEPhysicsList.hh"
#include "ICESPICEActionInitializer.hh"
int main(int argc,char** argv) {
//choose the Random engine
G4Random::setTheEngine(new CLHEP::RanecuEngine);
// Construct the default run manager
//
auto* runManager = G4RunManagerFactory::CreateRunManager();
G4int nThreads = 12;
runManager->SetNumberOfThreads(nThreads);
// set mandatory initialization classes
runManager->SetUserInitialization(new ICESPICEDetectorConstruction);
runManager->SetUserInitialization(new ICESPICEPhysicsList);
runManager->SetUserInitialization(new ICESPICEActionInitializer());
// visualization manager
G4VisManager* visManager = new G4VisExecutive;
visManager->Initialize();
//Initialize G4 kernel
runManager->Initialize();
// get the pointer to the User Interface manager
G4UImanager* UImanager = G4UImanager::GetUIpointer();
if (argc==1) // Define UI session for interactive mode.
{
G4UIExecutive* ui = new G4UIExecutive(argc, argv);
UImanager->ApplyCommand("/control/execute vis.mac");
if (ui->IsGUI()) {
UImanager->ApplyCommand("/control/execute gui.mac");
}
ui->SessionStart();
delete ui;
}
else // Batch mode
{
G4String command = "/control/execute ";
G4String fileName = argv[1];
UImanager->ApplyCommand(command+fileName);
}
// // Save histograms
G4AnalysisManager* man = G4AnalysisManager::Instance();
man->Write();
man->CloseFile();
// job termination
delete visManager;
delete runManager;
return 0;
}