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.2.p2)


  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 main.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 << 
 51                                                << 
 52 void Parse(int& argc, char** argv);                57 void Parse(int& argc, char** argv);
 53                                                    58 
 54 //....oooOO0OOooo........oooOO0OOooo........oo << 
 55                                                << 
 56 int main(int argc, char** argv)                    59 int main(int argc, char** argv)
 57 {                                                  60 {
 58   //////////                                       61   //////////
 59   // Parse options given in commandLine            62   // Parse options given in commandLine
 60   //                                               63   //
 61   Parse(argc, argv);                               64   Parse(argc, argv);
 62                                                    65 
 63   //////////                                       66   //////////
 64   // Construct the run manager according to wh     67   // Construct the run manager according to whether MT is activated or not
 65   //                                               68   //
 66   Command* commandLine(0);                         69   Command* commandLine(0);
 67                                                    70 
 68   auto* runManager = G4RunManagerFactory::Crea <<  71 #ifdef G4MULTITHREADED
 69   if ((commandLine = parser->GetCommandIfActiv <<  72   G4MTRunManager* runManager= new G4MTRunManager;
                                                   >>  73   if ((commandLine = parser->GetCommandIfActive("-mt")))
                                                   >>  74   {
 70     int nThreads = 2;                              75     int nThreads = 2;
 71     if (commandLine->GetOption() == "NMAX") {  <<  76     if(commandLine->GetOption() == "NMAX")
                                                   >>  77     {
 72       nThreads = G4Threading::G4GetNumberOfCor     78       nThreads = G4Threading::G4GetNumberOfCores();
 73     }                                              79     }
 74     else {                                     <<  80     else
                                                   >>  81     {
 75       nThreads = G4UIcommand::ConvertToInt(com     82       nThreads = G4UIcommand::ConvertToInt(commandLine->GetOption());
 76     }                                              83     }
 77     runManager->SetNumberOfThreads(nThreads);  <<  84     G4cout << "===== clustering is started with "
                                                   >>  85     << runManager->GetNumberOfThreads()
                                                   >>  86     << " threads =====" << G4endl;
 78                                                    87 
 79     G4cout << "===== clustering is started wit <<  88     runManager->SetNumberOfThreads(nThreads);
 80            << " threads =====" << G4endl;      << 
 81   }                                                89   }
                                                   >>  90 #else
                                                   >>  91   G4RunManager* runManager = new G4RunManager();
                                                   >>  92 #endif
 82                                                    93 
 83   //////////                                       94   //////////
 84   // Set mandatory user initialization classes     95   // Set mandatory user initialization classes
 85   //                                               96   //
 86   DetectorConstruction* detector = new Detecto     97   DetectorConstruction* detector = new DetectorConstruction;
 87   runManager->SetUserInitialization(new Physic     98   runManager->SetUserInitialization(new PhysicsList);
 88   runManager->SetUserInitialization(detector);     99   runManager->SetUserInitialization(detector);
 89   runManager->SetUserInitialization(new Action    100   runManager->SetUserInitialization(new ActionInitialization());
 90                                                   101 
 91   // Initialize G4 kernel                         102   // Initialize G4 kernel
 92   runManager->Initialize();                       103   runManager->Initialize();
 93                                                   104 
 94   // Initialize visualization                     105   // Initialize visualization
 95   G4VisManager* visManager = new G4VisExecutiv    106   G4VisManager* visManager = new G4VisExecutive;
 96   visManager->Initialize();                       107   visManager->Initialize();
 97                                                   108 
 98   // Get the pointer to the User Interface man    109   // Get the pointer to the User Interface manager
 99   G4UImanager* UImanager = G4UImanager::GetUIp    110   G4UImanager* UImanager = G4UImanager::GetUIpointer();
100   G4UIExecutive* ui(0);                           111   G4UIExecutive* ui(0);
101                                                   112 
102   // interactive mode : define UI session         113   // interactive mode : define UI session
103   if ((commandLine = parser->GetCommandIfActiv << 114   if ((commandLine = parser->GetCommandIfActive("-gui")))
                                                   >> 115   {
104     ui = new G4UIExecutive(argc, argv, command    116     ui = new G4UIExecutive(argc, argv, commandLine->GetOption());
105                                                   117 
106     if (ui->IsGUI()) UImanager->ApplyCommand("    118     if (ui->IsGUI()) UImanager->ApplyCommand("/control/execute gui.mac");
107                                                   119 
108     if (parser->GetCommandIfActive("-novis") =    120     if (parser->GetCommandIfActive("-novis") == 0)
109     // visualization is used by default           121     // visualization is used by default
110     {                                             122     {
111       if ((commandLine = parser->GetCommandIfA    123       if ((commandLine = parser->GetCommandIfActive("-vis")))
112       // select a visualization driver if need    124       // select a visualization driver if needed (e.g. HepFile)
113       {                                           125       {
114         UImanager->ApplyCommand(G4String("/vis << 126         UImanager->ApplyCommand(
                                                   >> 127             G4String("/vis/open ") + commandLine->GetOption());
115       }                                           128       }
116       else                                        129       else
117       // by default OGL is used                   130       // by default OGL is used
118       {                                           131       {
119         UImanager->ApplyCommand("/vis/open OGL    132         UImanager->ApplyCommand("/vis/open OGL 800x600-0+0");
120       }                                           133       }
121       UImanager->ApplyCommand("/control/execut    134       UImanager->ApplyCommand("/control/execute vis.mac");
122     }                                             135     }
123   }                                               136   }
124   else                                            137   else
125   // to be use visualization file (= store the    138   // to be use visualization file (= store the visualization into
126   // an external file:                            139   // an external file:
127   // ASCIITree ;  DAWNFILE ; HepRepFile ; VRML    140   // ASCIITree ;  DAWNFILE ; HepRepFile ; VRML(1,2)FILE ; gMocrenFile ...
128   {                                               141   {
129     if ((commandLine = parser->GetCommandIfAct << 142     if ((commandLine = parser->GetCommandIfActive("-vis")))
130       UImanager->ApplyCommand(G4String("/vis/o << 143     {
                                                   >> 144       UImanager->ApplyCommand(
                                                   >> 145           G4String("/vis/open ") + commandLine->GetOption());
131       UImanager->ApplyCommand("/control/execut    146       UImanager->ApplyCommand("/control/execute vis.mac");
132     }                                             147     }
133   }                                               148   }
134                                                   149 
135   if ((commandLine = parser->GetCommandIfActiv << 150   if ((commandLine = parser->GetCommandIfActive("-mac")))
                                                   >> 151   {
136     G4String command = "/control/execute ";       152     G4String command = "/control/execute ";
137     UImanager->ApplyCommand(command + commandL    153     UImanager->ApplyCommand(command + commandLine->GetOption());
138   }                                               154   }
139   else {                                       << 155   else
                                                   >> 156   {
140     UImanager->ApplyCommand("/control/execute     157     UImanager->ApplyCommand("/control/execute run.in");
141   }                                               158   }
142                                                   159 
143   if ((commandLine = parser->GetCommandIfActiv << 160   if ((commandLine = parser->GetCommandIfActive("-gui")))
                                                   >> 161   {
144     ui->SessionStart();                           162     ui->SessionStart();
145     delete ui;                                    163     delete ui;
146   }                                               164   }
147                                                   165 
148   // Job termination                              166   // Job termination
149   // Free the store: user actions, physics_lis    167   // Free the store: user actions, physics_list and detector_description are
150   // owned and deleted by the run manager, so     168   // owned and deleted by the run manager, so they should not be deleted
151   // in the main() program !                      169   // in the main() program !
152                                                   170 
153   delete visManager;                              171   delete visManager;
154   delete runManager;                              172   delete runManager;
155                                                   173 
156   CommandLineParser::DeleteInstance();            174   CommandLineParser::DeleteInstance();
157                                                   175 
158   return 0;                                       176   return 0;
159 }                                                 177 }
160                                                   178 
161 //....oooOO0OOooo........oooOO0OOooo........oo << 
162                                                << 
163 void Parse(int& argc, char** argv)                179 void Parse(int& argc, char** argv)
164 {                                                 180 {
165   //////////                                      181   //////////
166   // Parse options given in commandLine           182   // Parse options given in commandLine
167   //                                              183   //
168   parser = CommandLineParser::GetParser();        184   parser = CommandLineParser::GetParser();
169                                                   185 
170   parser->AddCommand("-gui", Command::OptionNo << 186   parser->AddCommand(
171                      "Select geant4 UI or just << 187       "-gui", Command::OptionNotCompulsory,
                                                   >> 188       "Select geant4 UI or just launch a geant4 terminal session", "qt");
                                                   >> 189 
                                                   >> 190   parser->AddCommand("-mac", Command::WithOption, "Give a mac file to execute",
                                                   >> 191                      "macFile.mac");
                                                   >> 192 
                                                   >> 193 // You cann your own command, as for instance:
                                                   >> 194 //  parser->AddCommand("-seed", 
                                                   >> 195 //                     Command::WithOption,
                                                   >> 196 //                     "Give a seed value in argument to be tested", "seed");
                                                   >> 197 // it is then up to you to manage this option
                                                   >> 198 
                                                   >> 199 #ifdef G4MULTITHREADED
                                                   >> 200   parser->AddCommand("-mt",
                                                   >> 201       Command::WithOption,
                                                   >> 202       "Launch in MT mode (events computed in parallel,"
                                                   >> 203       " NOT RECOMMENDED WITH CHEMISTRY)", "2");
                                                   >> 204 #endif
172                                                   205 
173   parser->AddCommand("-mac", Command::WithOpti << 206   parser->AddCommand("-vis", Command::WithOption,
                                                   >> 207                      "Select a visualization driver", "OGL 600x600-0+0");
174                                                   208 
175   // You cann your own command, as for instanc << 209   parser->AddCommand("-novis", Command::WithoutOption,
176   //  parser->AddCommand("-seed",              << 210                      "Deactivate visualization when using GUI");
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                                                << 
186   parser->AddCommand("-vis", Command::WithOpti << 
187                      "OGL 600x600-0+0");       << 
188                                                << 
189   parser->AddCommand("-novis", Command::Withou << 
190                                                   211 
191   //////////                                      212   //////////
192   // If -h or --help is given in option : prin    213   // If -h or --help is given in option : print help and exit
193   //                                              214   //
194   if (parser->Parse(argc, argv) != 0)  // help << 215   if (parser->Parse(argc, argv) != 0) // help is being printed
195   {                                               216   {
196     // if you are using ROOT, create a TApplic    217     // if you are using ROOT, create a TApplication in this condition in order
197     // to print the help from ROOT as well        218     // to print the help from ROOT as well
198     CommandLineParser::DeleteInstance();          219     CommandLineParser::DeleteInstance();
199     // Exit is kept                            << 
200     std::exit(0);                                 220     std::exit(0);
201   }                                               221   }
202                                                   222 
203   ///////////                                     223   ///////////
204   // Kill application if wrong argument in com    224   // Kill application if wrong argument in command line
205   //                                              225   //
206   if (parser->CheckIfNotHandledOptionsExists(a << 226   if (parser->CheckIfNotHandledOptionsExists(argc, argv))
                                                   >> 227   {
207     // if you are using ROOT, you should initi    228     // if you are using ROOT, you should initialise your TApplication
208     // before this condition                      229     // before this condition
209     abort();                                      230     abort();
210   }                                               231   }
211 }                                                 232 }
212                                                   233