pps-timing/main.cc at master · forthommel/pps-timing · 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
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
127
128
129
130
#ifdef G4MULTITHREADED
#include "G4MTRunManager.hh"
#else
#include "G4RunManager.hh"
#endif
#include "G4UImanager.hh"
#include "G4UIterminal.hh"
#include "G4UItcsh.hh"

#include "G4ios.hh"

//#include "LHEP.hh" //FIXME
#include "QGSP_BERT.hh"
#include "LBE.hh"

#include "G4OpticalPhysics.hh"
#include "G4VModularPhysicsList.hh"

#include "globals.hh"

#include "SteppingVerbose.hh"
#include "ActionInitialization.hh"
#include "GeometryConstructor.hh"

#ifdef G4UI_USE
#include "G4UIExecutive.hh"
#endif

#include "Randomize.hh"

#ifdef G4VIS_USE
#include "G4VisExecutive.hh"
#endif

int main(int argc,char** argv)
{
  // Seed the random number generator manually
  //
  //G4long myseed = 345354;
  //CLHEP::HepRandom::setTheSeed(myseed);

  // User Verbose output class
  //
  G4VSteppingVerbose* verbosity = new PPS::SteppingVerbose;
  G4VSteppingVerbose::SetInstance(verbosity);

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

  // UserInitialization classes - mandatory
  //
  PPS::GeometryConstructor *geom_constructor = new PPS::GeometryConstructor;

  runManager->SetUserInitialization(geom_constructor);

  G4VModularPhysicsList* physics;

  //physics = new QuartLPhysicsList;
  //physics = new LHEP;    //22.01 from AK
  physics = new QGSP_BERT;
  //physics = new LBE;
  physics->SetVerboseLevel(-1);
  physics->RegisterPhysics(new G4OpticalPhysics);

  runManager-> SetUserInitialization(physics);

  PPS::ActionInitialization* action = new PPS::ActionInitialization;
  runManager->SetUserInitialization(action);

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

  /*for (G4int i=0; i<20; i++) {
    const G4ThreeVector cell_center = detector->GetCellCenter(1, i);
    G4cout << "--> Cell center for id = " << i << " -> (" << cell_center.x() << ", " << cell_center.y() << ", " << cell_center.z() << ")" << G4endl;
  }*/

  // Get the pointer to the User Interface manager
  //
  G4UImanager* UImanager = G4UImanager::GetUIpointer();

  if (argc==1) { // Define UI session for interactive mode
    G4UIExecutive* ui;
    ui = new G4UIExecutive(argc, argv);
    // Initialize visualization
#ifdef G4VIS_USE
    // G4VisExecutive can take a verbosity argument - see /vis/verbose guidance.
    G4VisManager* visManager = new G4VisExecutive("Quiet");
    visManager->Initialize();
#endif
    UImanager->ApplyCommand("/control/execute test/init_vis.mac");
    ui->SessionStart();
#ifdef G4VIS_USE
  delete visManager;
#endif
  }
  else { // Batch mode
    G4String command = "/control/execute ";
    G4String fileName = argv[1];
    G4UIExecutive* ui;
    //ui = new G4UIExecutive(argc, argv);
#ifdef G4VIS_USE
    // G4VisExecutive can take a verbosity argument - see /vis/verbose guidance.
    G4VisManager* visManager = new G4VisExecutive("Quiet");
    visManager->Initialize();
#endif
    UImanager->ApplyCommand(command+fileName);
    //ui->SessionStart();
#ifdef G4VIS_USE
    delete visManager;
#endif
  }

  // Job termination
  // Free the store: user actions, physics_list and detector_description are
  //                 owned and deleted by the run manager, so they should not
  //                 be deleted in the main() program !

    delete runManager;
    delete verbosity;

    return 0;
}