Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/medical/dna/pdb4dna/pdb4dna.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

  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 // This example is provided by the Geant4-DNA collaboration
 27 // Any report or published results obtained using the Geant4-DNA software
 28 // shall cite the following Geant4-DNA collaboration publication:
 29 // Med. Phys. 37 (2010) 4692-4708
 30 // Delage et al. PDB4DNA: implementation of DNA geometry from the Protein Data
 31 //                  Bank (PDB) description for Geant4-DNA Monte-Carlo
 32 //                  simulations (submitted to Comput. Phys. Commun.)
 33 // The Geant4-DNA web site is available at http://geant4-dna.org
 34 //
 35 //
 36 /// \file pdb4dna.cc
 37 /// \brief Main program of the pdb4dna example
 38 
 39 #include "ActionInitialization.hh"
 40 #include "DetectorConstruction.hh"
 41 #include "PhysicsList.hh"
 42 
 43 #ifdef G4MULTITHREADED
 44 #  include "G4MTRunManager.hh"
 45 #else
 46 #  include "G4RunManager.hh"
 47 #endif
 48 
 49 #include "G4UIExecutive.hh"
 50 #include "G4UImanager.hh"
 51 #include "G4VisExecutive.hh"
 52 #include "Randomize.hh"
 53 #ifdef G4UI_USE_QT
 54 #  include "G4UIQt.hh"
 55 #endif
 56 #include "CommandLineParser.hh"
 57 
 58 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 59 
 60 using namespace G4DNAPARSER;
 61 CommandLineParser* parser(0);
 62 
 63 void Parse(int& argc, char** argv);
 64 
 65 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 66 
 67 int main(int argc, char** argv)
 68 {
 69   //////////
 70   // Parse options given in commandLine
 71   //
 72   Parse(argc, argv);
 73 
 74   //////////
 75   // Construct the run manager according to whether MT is activated or not
 76   //
 77   Command* commandLine(0);
 78 
 79 #ifdef G4MULTITHREADED
 80   G4MTRunManager* runManager = new G4MTRunManager;
 81   if ((commandLine = parser->GetCommandIfActive("-mt"))) {
 82     int nThreads = 2;
 83     if (commandLine->GetOption() == "NMAX") {
 84       nThreads = G4Threading::G4GetNumberOfCores();
 85     }
 86     else {
 87       nThreads = G4UIcommand::ConvertToInt(commandLine->GetOption());
 88     }
 89     G4cout << "===== PDB4DNA is started with " << runManager->GetNumberOfThreads()
 90            << " threads =====" << G4endl;
 91 
 92     runManager->SetNumberOfThreads(nThreads);
 93   }
 94 #else
 95   G4RunManager* runManager = new G4RunManager();
 96 #endif
 97 
 98   // Set user classes
 99   //
100   runManager->SetUserInitialization(new DetectorConstruction());
101   runManager->SetUserInitialization(new PhysicsList);
102   runManager->SetUserInitialization(new ActionInitialization());
103 
104   // Initialize G4 keRnel
105   //
106   // runManager->Initialize();
107 
108   // Initialize visualization
109   G4VisManager* visManager = new G4VisExecutive;
110   visManager->Initialize();
111 
112   // Get the pointer to the User Interface manager
113   G4UImanager* UImanager = G4UImanager::GetUIpointer();
114   G4UIExecutive* ui(0);
115 
116   // interactive mode : define UI session
117   if ((commandLine = parser->GetCommandIfActive("-gui"))) {
118     ui = new G4UIExecutive(argc, argv, commandLine->GetOption());
119 
120     if (parser->GetCommandIfActive("-novis") == 0)
121     // visualization is used by default
122     {
123       if ((commandLine = parser->GetCommandIfActive("-vis")))
124       // select a visualization driver if needed (e.g. HepFile)
125       {
126         UImanager->ApplyCommand(G4String("/vis/open ") + commandLine->GetOption());
127       }
128       else
129       // by default OGL is used
130       {
131         UImanager->ApplyCommand("/vis/open OGL 600x600-0+0");
132       }
133       UImanager->ApplyCommand("/control/execute vis.mac");
134     }
135     if (ui->IsGUI()) UImanager->ApplyCommand("/control/execute gui.mac");
136   }
137   else
138   // to be use visualization file (= store the visualization into
139   // an external file:
140   // ASCIITree ;  DAWNFILE ; HepRepFile ; VRML(1,2)FILE ; gMocrenFile ...
141   {
142     if ((commandLine = parser->GetCommandIfActive("-vis"))) {
143       UImanager->ApplyCommand(G4String("/vis/open ") + commandLine->GetOption());
144       UImanager->ApplyCommand("/control/execute vis.mac");
145     }
146   }
147 
148   if ((commandLine = parser->GetCommandIfActive("-mac"))) {
149     G4String command = "/control/execute ";
150     UImanager->ApplyCommand(command + commandLine->GetOption());
151   }
152   else {
153     UImanager->ApplyCommand("/control/execute init.mac");
154   }
155 
156   if ((commandLine = parser->GetCommandIfActive("-gui"))) {
157 #ifdef G4UI_USE_QT
158     G4UIQt* UIQt = static_cast<G4UIQt*>(UImanager->GetG4UIWindow());
159     if (UIQt) {
160       UIQt->AddViewerTabFromFile("README", "README from " + G4String(argv[0]));
161     }
162 #endif
163     ui->SessionStart();
164     delete ui;
165   }
166 
167   // Job termination
168   delete visManager;
169   delete runManager;
170 
171   return 0;
172 }
173 
174 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
175 
176 void Parse(int& argc, char** argv)
177 {
178   //////////
179   // Parse options given in commandLine
180   //
181   parser = CommandLineParser::GetParser();
182 
183   parser->AddCommand("-gui", Command::OptionNotCompulsory,
184                      "Select geant4 UI or just launch a geant4 terminal session", "qt");
185 
186   parser->AddCommand("-mac", Command::WithOption, "Give a mac file to execute", "pdb4dna.in");
187 
188   // You cann your own command, as for instance:
189   //  parser->AddCommand("-seed",
190   //                     Command::WithOption,
191   //                     "Give a seed value in argument to be tested", "seed");
192   // it is then up to you to manage this option
193 
194 #ifdef G4MULTITHREADED
195   parser->AddCommand("-mt", Command::WithOption, "Launch in MT mode (events computed in parallel)",
196                      "2");
197 #endif
198 
199   parser->AddCommand("-vis", Command::WithOption, "Select a visualization driver",
200                      "OGL 600x600-0+0");
201 
202   parser->AddCommand("-novis", Command::WithoutOption, "Deactivate visualization when using GUI");
203 
204   //////////
205   // If -h or --help is given in option : print help and exit
206   //
207   if (parser->Parse(argc, argv) != 0)  // help is being printed
208   {
209     // if you are using ROOT, create a TApplication in this condition in order
210     // to print the help from ROOT as well
211     CommandLineParser::DeleteInstance();
212     std::exit(0);
213   }
214 
215   ///////////
216   // Kill application if wrong argument in command line
217   //
218   if (parser->CheckIfNotHandledOptionsExists(argc, argv)) {
219     // if you are using ROOT, you should initialise your TApplication
220     // before this condition
221     abort();
222   }
223 }
224