fcal/fcal.cc at main · mpotlog/fcal · 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
#include "G4RunManager.hh"
#include "G4UImanager.hh"
#include "G4UIExecutive.hh"

#include "DetectorConstruction.hh"
#include "PhysicsList.hh"
#include "ActionInitialization.hh"

#include "G4PhysListFactory.hh"
#include "G4SystemOfUnits.hh"
#include "G4VisExecutive.hh"


int main(int argc, char** argv)
{
  G4cout << "Aplicatie pornita..."  << G4endl;

  G4UIExecutive* ui = 0;
  if ( argc == 1 ) {
    ui = new G4UIExecutive(argc, argv);
  }

  // Run manager
  G4RunManager * runManager = new G4RunManager;

  // Mandatory initialization classes
  runManager->SetUserInitialization(new DetectorConstruction);
  runManager->SetUserInitialization(new PhysicsList);

 // //G4PhysListFactory factory;
  //G4VModularPhysicsList *physicsList = factory.GetReferencePhysList("FTFP_BERT");
 // G4VModularPhysicsList *physicsList = factory.GetReferencePhysList("QGSP_BERT");
 //  physicsList->SetDefaultCutValue(0.001*mm);
 // runManager->SetUserInitialization(physicsList);

  // Optional user action classes
   runManager->SetUserInitialization(new ActionInitialization());

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

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

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

  // Process macro or start UI session
  if ( ! ui ) {
    // batch mode
    G4String command = "/control/execute ";
    G4String fileName = argv[1];
    UImanager->ApplyCommand(command+fileName);
  }
  else {
    // interactive mode
    UImanager->ApplyCommand("/control/execute vrml1.mac");
    ui->SessionStart();
    delete ui;
  }

  // Delete the managers
  delete visManager;
  delete runManager;

  G4cout << "Aplicatie incheiata cu succes."  << G4endl;

  return 0;
}