Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/medical/dna/clustering/clustering.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 ]

Diff markup

Differences between /examples/extended/medical/dna/clustering/clustering.cc (Version 11.3.0) and /examples/extended/medical/dna/clustering/clustering.cc (Version 10.4.p3)


  1 //                                                  1 //
  2 // *******************************************      2 // ********************************************************************
  3 // * License and Disclaimer                         3 // * License and Disclaimer                                           *
  4 // *                                                4 // *                                                                  *
  5 // * The  Geant4 software  is  copyright of th      5 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
  6 // * the Geant4 Collaboration.  It is provided      6 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
  7 // * conditions of the Geant4 Software License      7 // * conditions of the Geant4 Software License,  included in the file *
  8 // * LICENSE and available at  http://cern.ch/      8 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
  9 // * include a list of copyright holders.           9 // * include a list of copyright holders.                             *
 10 // *                                               10 // *                                                                  *
 11 // * Neither the authors of this software syst     11 // * Neither the authors of this software system, nor their employing *
 12 // * institutes,nor the agencies providing fin     12 // * institutes,nor the agencies providing financial support for this *
 13 // * work  make  any representation or  warran     13 // * work  make  any representation or  warranty, express or implied, *
 14 // * regarding  this  software system or assum     14 // * regarding  this  software system or assume any liability for its *
 15 // * use.  Please see the license in the file      15 // * use.  Please see the license in the file  LICENSE  and URL above *
 16 // * for the full disclaimer and the limitatio     16 // * for the full disclaimer and the limitation of liability.         *
 17 // *                                               17 // *                                                                  *
 18 // * This  code  implementation is the result      18 // * This  code  implementation is the result of  the  scientific and *
 19 // * technical work of the GEANT4 collaboratio     19 // * technical work of the GEANT4 collaboration.                      *
 20 // * By using,  copying,  modifying or  distri     20 // * By using,  copying,  modifying or  distributing the software (or *
 21 // * any work based  on the software)  you  ag     21 // * any work based  on the software)  you  agree  to acknowledge its *
 22 // * use  in  resulting  scientific  publicati     22 // * use  in  resulting  scientific  publications,  and indicate your *
 23 // * acceptance of all terms of the Geant4 Sof     23 // * acceptance of all terms of the Geant4 Software license.          *
 24 // *******************************************     24 // ********************************************************************
 25 //                                                 25 //
 26 // This example is provided by the Geant4-DNA      26 // This example is provided by the Geant4-DNA collaboration
 27 // Any report or published results obtained us     27 // Any report or published results obtained using the Geant4-DNA software
 28 // shall cite the following Geant4-DNA collabo     28 // shall cite the following Geant4-DNA collaboration publication:
 29 // Med. Phys. 37 (2010) 4692-4708                  29 // Med. Phys. 37 (2010) 4692-4708
 30 // The Geant4-DNA web site is available at htt     30 // The Geant4-DNA web site is available at http://geant4-dna.org
 31 //                                                 31 //
                                                   >>  32 // $Id$
 32 //                                                 33 //
 33 /// \file clustering.cc                            34 /// \file clustering.cc
 34 /// \brief clustering example                      35 /// \brief clustering example
 35                                                    36 
 36 #include "ActionInitialization.hh"             << 
 37 #include "CommandLineParser.hh"                << 
 38 #include "DetectorConstruction.hh"                 37 #include "DetectorConstruction.hh"
 39 #include "PhysicsList.hh"                          38 #include "PhysicsList.hh"
                                                   >>  39 #include "ActionInitialization.hh"
                                                   >>  40 
                                                   >>  41 #ifdef G4MULTITHREADED
                                                   >>  42 #include "G4MTRunManager.hh"
                                                   >>  43 #else
                                                   >>  44 #include "G4RunManager.hh"
                                                   >>  45 #endif
 40                                                    46 
 41 #include "G4DNAChemistryManager.hh"                47 #include "G4DNAChemistryManager.hh"
 42 #include "G4RunManagerFactory.hh"              << 
 43 #include "G4UIExecutive.hh"                    << 
 44 #include "G4UImanager.hh"                          48 #include "G4UImanager.hh"
                                                   >>  49 #include "G4UIExecutive.hh"
 45 #include "G4VisExecutive.hh"                       50 #include "G4VisExecutive.hh"
 46                                                    51 
                                                   >>  52 #include "CommandLineParser.hh"
                                                   >>  53 
 47 using namespace G4DNAPARSER;                       54 using namespace G4DNAPARSER;
 48 CommandLineParser* parser(0);                      55 CommandLineParser* parser(0);
 49                                                    56 
 50 //....oooOO0OOooo........oooOO0OOooo........oo     57 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 51                                                    58 
 52 void Parse(int& argc, char** argv);                59 void Parse(int& argc, char** argv);
 53                                                    60 
 54 //....oooOO0OOooo........oooOO0OOooo........oo     61 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 55                                                    62 
 56 int main(int argc, char** argv)                    63 int main(int argc, char** argv)
 57 {                                                  64 {
 58   //////////                                       65   //////////
 59   // Parse options given in commandLine            66   // Parse options given in commandLine
 60   //                                               67   //
 61   Parse(argc, argv);                               68   Parse(argc, argv);
 62                                                    69 
 63   //////////                                       70   //////////
 64   // Construct the run manager according to wh     71   // Construct the run manager according to whether MT is activated or not
 65   //                                               72   //
 66   Command* commandLine(0);                         73   Command* commandLine(0);
 67                                                    74 
 68   auto* runManager = G4RunManagerFactory::Crea <<  75 #ifdef G4MULTITHREADED
 69   if ((commandLine = parser->GetCommandIfActiv <<  76   G4MTRunManager* runManager= new G4MTRunManager;
                                                   >>  77   if ((commandLine = parser->GetCommandIfActive("-mt")))
                                                   >>  78   {
 70     int nThreads = 2;                              79     int nThreads = 2;
 71     if (commandLine->GetOption() == "NMAX") {  <<  80     if(commandLine->GetOption() == "NMAX")
                                                   >>  81     {
 72       nThreads = G4Threading::G4GetNumberOfCor     82       nThreads = G4Threading::G4GetNumberOfCores();
 73     }                                              83     }
 74     else {                                     <<  84     else
                                                   >>  85     {
 75       nThreads = G4UIcommand::ConvertToInt(com     86       nThreads = G4UIcommand::ConvertToInt(commandLine->GetOption());
 76     }                                              87     }
 77     runManager->SetNumberOfThreads(nThreads);  <<  88     G4cout << "===== clustering is started with "
                                                   >>  89     << runManager->GetNumberOfThreads()
                                                   >>  90     << " threads =====" << G4endl;
 78                                                    91 
 79     G4cout << "===== clustering is started wit <<  92     runManager->SetNumberOfThreads(nThreads);
 80            << " threads =====" << G4endl;      << 
 81   }                                                93   }
                                                   >>  94 #else
                                                   >>  95   G4RunManager* runManager = new G4RunManager();
                                                   >>  96 #endif
 82                                                    97 
 83   //////////                                       98   //////////
 84   // Set mandatory user initialization classes     99   // Set mandatory user initialization classes
 85   //                                              100   //
 86   DetectorConstruction* detector = new Detecto    101   DetectorConstruction* detector = new DetectorConstruction;
 87   runManager->SetUserInitialization(new Physic    102   runManager->SetUserInitialization(new PhysicsList);
 88   runManager->SetUserInitialization(detector);    103   runManager->SetUserInitialization(detector);
 89   runManager->SetUserInitialization(new Action    104   runManager->SetUserInitialization(new ActionInitialization());
 90                                                   105 
 91   // Initialize G4 kernel                         106   // Initialize G4 kernel
 92   runManager->Initialize();                       107   runManager->Initialize();
 93                                                   108 
 94   // Initialize visualization                     109   // Initialize visualization
 95   G4VisManager* visManager = new G4VisExecutiv    110   G4VisManager* visManager = new G4VisExecutive;
 96   visManager->Initialize();                       111   visManager->Initialize();
 97                                                   112 
 98   // Get the pointer to the User Interface man    113   // Get the pointer to the User Interface manager
 99   G4UImanager* UImanager = G4UImanager::GetUIp    114   G4UImanager* UImanager = G4UImanager::GetUIpointer();
100   G4UIExecutive* ui(0);                           115   G4UIExecutive* ui(0);
101                                                   116 
102   // interactive mode : define UI session         117   // interactive mode : define UI session
103   if ((commandLine = parser->GetCommandIfActiv << 118   if ((commandLine = parser->GetCommandIfActive("-gui")))
                                                   >> 119   {
104     ui = new G4UIExecutive(argc, argv, command    120     ui = new G4UIExecutive(argc, argv, commandLine->GetOption());
105                                                   121 
106     if (ui->IsGUI()) UImanager->ApplyCommand("    122     if (ui->IsGUI()) UImanager->ApplyCommand("/control/execute gui.mac");
107                                                   123 
108     if (parser->GetCommandIfActive("-novis") =    124     if (parser->GetCommandIfActive("-novis") == 0)
109     // visualization is used by default           125     // visualization is used by default
110     {                                             126     {
111       if ((commandLine = parser->GetCommandIfA    127       if ((commandLine = parser->GetCommandIfActive("-vis")))
112       // select a visualization driver if need    128       // select a visualization driver if needed (e.g. HepFile)
113       {                                           129       {
114         UImanager->ApplyCommand(G4String("/vis << 130         UImanager->ApplyCommand(
                                                   >> 131             G4String("/vis/open ") + commandLine->GetOption());
115       }                                           132       }
116       else                                        133       else
117       // by default OGL is used                   134       // by default OGL is used
118       {                                           135       {
119         UImanager->ApplyCommand("/vis/open OGL    136         UImanager->ApplyCommand("/vis/open OGL 800x600-0+0");
120       }                                           137       }
121       UImanager->ApplyCommand("/control/execut    138       UImanager->ApplyCommand("/control/execute vis.mac");
122     }                                             139     }
123   }                                               140   }
124   else                                            141   else
125   // to be use visualization file (= store the    142   // to be use visualization file (= store the visualization into
126   // an external file:                            143   // an external file:
127   // ASCIITree ;  DAWNFILE ; HepRepFile ; VRML    144   // ASCIITree ;  DAWNFILE ; HepRepFile ; VRML(1,2)FILE ; gMocrenFile ...
128   {                                               145   {
129     if ((commandLine = parser->GetCommandIfAct << 146     if ((commandLine = parser->GetCommandIfActive("-vis")))
130       UImanager->ApplyCommand(G4String("/vis/o << 147     {
                                                   >> 148       UImanager->ApplyCommand(
                                                   >> 149           G4String("/vis/open ") + commandLine->GetOption());
131       UImanager->ApplyCommand("/control/execut    150       UImanager->ApplyCommand("/control/execute vis.mac");
132     }                                             151     }
133   }                                               152   }
134                                                   153 
135   if ((commandLine = parser->GetCommandIfActiv << 154   if ((commandLine = parser->GetCommandIfActive("-mac")))
                                                   >> 155   {
136     G4String command = "/control/execute ";       156     G4String command = "/control/execute ";
137     UImanager->ApplyCommand(command + commandL    157     UImanager->ApplyCommand(command + commandLine->GetOption());
138   }                                               158   }
139   else {                                       << 159   else
                                                   >> 160   {
140     UImanager->ApplyCommand("/control/execute     161     UImanager->ApplyCommand("/control/execute run.in");
141   }                                               162   }
142                                                   163 
143   if ((commandLine = parser->GetCommandIfActiv << 164   if ((commandLine = parser->GetCommandIfActive("-gui")))
                                                   >> 165   {
144     ui->SessionStart();                           166     ui->SessionStart();
145     delete ui;                                    167     delete ui;
146   }                                               168   }
147                                                   169 
148   // Job termination                              170   // Job termination
149   // Free the store: user actions, physics_lis    171   // Free the store: user actions, physics_list and detector_description are
150   // owned and deleted by the run manager, so     172   // owned and deleted by the run manager, so they should not be deleted
151   // in the main() program !                      173   // in the main() program !
152                                                   174 
153   delete visManager;                              175   delete visManager;
154   delete runManager;                              176   delete runManager;
155                                                   177 
156   CommandLineParser::DeleteInstance();            178   CommandLineParser::DeleteInstance();
157                                                   179 
158   return 0;                                       180   return 0;
159 }                                                 181 }
160                                                   182 
161 //....oooOO0OOooo........oooOO0OOooo........oo    183 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
162                                                   184 
163 void Parse(int& argc, char** argv)                185 void Parse(int& argc, char** argv)
164 {                                                 186 {
165   //////////                                      187   //////////
166   // Parse options given in commandLine           188   // Parse options given in commandLine
167   //                                              189   //
168   parser = CommandLineParser::GetParser();        190   parser = CommandLineParser::GetParser();
169                                                   191 
170   parser->AddCommand("-gui", Command::OptionNo << 192   parser->AddCommand(
171                      "Select geant4 UI or just << 193       "-gui", Command::OptionNotCompulsory,
172                                                << 194       "Select geant4 UI or just launch a geant4 terminal session", "qt");
173   parser->AddCommand("-mac", Command::WithOpti << 195 
                                                   >> 196   parser->AddCommand("-mac", Command::WithOption, "Give a mac file to execute",
                                                   >> 197                      "macFile.mac");
                                                   >> 198 
                                                   >> 199 // You cann your own command, as for instance:
                                                   >> 200 //  parser->AddCommand("-seed", 
                                                   >> 201 //                     Command::WithOption,
                                                   >> 202 //                     "Give a seed value in argument to be tested", "seed");
                                                   >> 203 // it is then up to you to manage this option
                                                   >> 204 
                                                   >> 205 #ifdef G4MULTITHREADED
                                                   >> 206   parser->AddCommand("-mt",
                                                   >> 207       Command::WithOption,
                                                   >> 208       "Launch in MT mode (events computed in parallel,"
                                                   >> 209       " NOT RECOMMENDED WITH CHEMISTRY)", "2");
                                                   >> 210 #endif
174                                                   211 
175   // You cann your own command, as for instanc << 212   parser->AddCommand("-vis", Command::WithOption,
176   //  parser->AddCommand("-seed",              << 213                      "Select a visualization driver", "OGL 600x600-0+0");
177   //                     Command::WithOption,  << 
178   //                     "Give a seed value in << 
179   // it is then up to you to manage this optio << 
180                                                << 
181   parser->AddCommand("-mt", Command::WithOptio << 
182                      "Launch in MT mode (event << 
183                      " NOT RECOMMENDED WITH CH << 
184                      "2");                     << 
185                                                   214 
186   parser->AddCommand("-vis", Command::WithOpti << 215   parser->AddCommand("-novis", Command::WithoutOption,
187                      "OGL 600x600-0+0");       << 216                      "Deactivate visualization when using GUI");
188                                                << 
189   parser->AddCommand("-novis", Command::Withou << 
190                                                   217 
191   //////////                                      218   //////////
192   // If -h or --help is given in option : prin    219   // If -h or --help is given in option : print help and exit
193   //                                              220   //
194   if (parser->Parse(argc, argv) != 0)  // help << 221   if (parser->Parse(argc, argv) != 0) // help is being printed
195   {                                               222   {
196     // if you are using ROOT, create a TApplic    223     // if you are using ROOT, create a TApplication in this condition in order
197     // to print the help from ROOT as well        224     // to print the help from ROOT as well
198     CommandLineParser::DeleteInstance();          225     CommandLineParser::DeleteInstance();
199     // Exit is kept                               226     // Exit is kept
200     std::exit(0);                                 227     std::exit(0);
201   }                                               228   }
202                                                   229 
203   ///////////                                     230   ///////////
204   // Kill application if wrong argument in com    231   // Kill application if wrong argument in command line
205   //                                              232   //
206   if (parser->CheckIfNotHandledOptionsExists(a << 233   if (parser->CheckIfNotHandledOptionsExists(argc, argv))
                                                   >> 234   {
207     // if you are using ROOT, you should initi    235     // if you are using ROOT, you should initialise your TApplication
208     // before this condition                      236     // before this condition
209     abort();                                      237     abort();
210   }                                               238   }
211 }                                                 239 }
212                                                   240