QuantumCosmicDetector/sim.cc at dev · carolFR13/QuantumCosmicDetector · GitHub

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 <iostream>

#include "G4RunManager.hh"
#include "G4MTRunManager.hh"
#include "G4UImanager.hh"
#include "G4VisManager.hh"
#include "G4VisExecutive.hh"
#include "G4UIExecutive.hh"
#include "G4AnalysisManager.hh"

#include "QDPhysicsList.hh"
#include "QDDetectorConstruction.hh"
#include "QDActionInitialization.hh"

int main(int argc, char** argv){

    #ifdef G4MULTITHREADED
        G4MTRunManager *runManager = new G4MTRunManager;
    #else
        G4RunManager *runManager = new G4RunManager;
    #endif

    // Detector Construction
    runManager -> SetUserInitialization(new QDDetectorConstruction());

    // Physics list
    runManager -> SetUserInitialization(new QDPhysicsList());

    // Action Initialization
    runManager -> SetUserInitialization(new QDActionInitialization());

    // Initialize G4 kernel
    runManager -> Initialize();

    // user interface
    G4UIExecutive *ui = 0;

    // if no arguments are provided, start the UI session
    if (argc == 1) {
        ui = new G4UIExecutive(argc, argv);
    }

    G4VisManager* visManager = new G4VisExecutive();
    visManager->Initialize();

    G4UImanager *UImanager = G4UImanager::GetUIpointer();

    if(ui){
        // If we are in interactive mode, we can use the UI
        UImanager->ApplyCommand("/control/execute vis.mac");
        ui -> SessionStart();
    } else {
        // If we are in batch mode, we can execute a macro file
        G4String command = "/control/execute ";
        G4String fileName = argv[1];
        UImanager->ApplyCommand(command + fileName);
    }


    delete ui;
    delete visManager;
    delete runManager;

    return 0;

}