Geant4  10.02.p01
pythia6_decayer.cc
Go to the documentation of this file.
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 // $Id: pythia6_decayer.cc 81443 2014-05-28 14:26:55Z gcosmo $
27 //
30 
31 #include "P6DExtDecayerPhysics.hh"
32 
33 #include "ExG4DetectorConstruction01.hh"
34 #include "ExG4PrimaryGeneratorAction01.hh"
35 #include "ExG4RunAction01.hh"
36 #include "ExG4EventAction01.hh"
37 
38 #include "G4RunManager.hh"
39 #include "G4UImanager.hh"
40 #include "G4ThreeVector.hh"
41 #include "QGSP_BERT.hh"
42 #include "G4SystemOfUnits.hh"
43 
44 #ifdef G4VIS_USE
45 #include "G4VisExecutive.hh"
46 #endif
47 
48 #ifdef G4UI_USE
49 #include "G4UIExecutive.hh"
50 #endif
51 
52 #include "Randomize.hh"
53 
54 
55 int main(int argc,char** argv)
56 {
57  // Choose the Random engine
58  //
59  G4Random::setTheEngine(new CLHEP::RanecuEngine);
60 
61  // Construct the default run manager
62  //
63  G4RunManager * runManager = new G4RunManager;
64 
65  // Set mandatory initialization classes
66  //
68 
69  //
70  G4VModularPhysicsList* physicsList = new QGSP_BERT;
71  physicsList->RegisterPhysics(new P6DExtDecayerPhysics());
72  runManager->SetUserInitialization(physicsList);
73 
74 
75  // Set user action classes
76  //
77  runManager->SetUserAction(
78  new ExG4PrimaryGeneratorAction01("B-", 50.*MeV));
79  // B- meson has not defined decay in Geant4
80  runManager->SetUserAction(new ExG4RunAction01);
81  runManager->SetUserAction(new ExG4EventAction01);
82 
83  runManager->Initialize();
84 
85 #ifdef G4VIS_USE
86  // Initialize visualization
87  G4VisManager* visManager = new G4VisExecutive;
88  // G4VisExecutive can take a verbosity argument - see /vis/verbose guidance.
89  // G4VisManager* visManager = new G4VisExecutive("Quiet");
90  visManager->Initialize();
91 #endif
92 
93  // Get the pointer to the User Interface manager
94  //
95  G4UImanager* UImanager = G4UImanager::GetUIpointer();
96 
97  if (argc!=1) {
98  // batch mode{
99  G4String command = "/control/execute ";
100  G4String fileName = argv[1];
101  UImanager->ApplyCommand(command+fileName);
102  }
103  else { // interactive mode : define UI session
104 #ifdef G4UI_USE
105  G4UIExecutive* ui = new G4UIExecutive(argc, argv);
106 #ifdef G4VIS_USE
107  UImanager->ApplyCommand("/control/execute init_vis.mac");
108 #else
109  UImanager->ApplyCommand("/control/execute init.mac");
110 #endif
111  ui->SessionStart();
112  delete ui;
113 #endif
114  }
115 
116 #ifdef G4VIS_USE
117  delete visManager;
118 #endif
119  // Job termination
120  // Free the store: user actions, physics_list and detector_description are
121  // owned and deleted by the run manager, so they should not
122  // be deleted in the main() program !
123  delete runManager;
124 
125  return 0;
126 }
virtual void SetUserInitialization(G4VUserDetectorConstruction *userInit)
static const double MeV
Definition: G4SIunits.hh:211
TQGSP_BERT< G4VModularPhysicsList > QGSP_BERT
Definition: QGSP_BERT.hh:63
void RegisterPhysics(G4VPhysicsConstructor *)
The builder for external decayer.
int main(int argc, char **argv)
Event action class.
Simple detector construction with only a world volume.
static G4UImanager * GetUIpointer()
Definition: G4UImanager.cc:58
void Initialize()
Run action class.
Definition of the P6DExtDecayerPhysics class.
virtual void Initialize()
The primary generator class with particle gun.
G4int ApplyCommand(const char *aCommand)
Definition: G4UImanager.cc:446
virtual void SetUserAction(G4UserRunAction *userAction)