Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/advanced/brachytherapy/Brachy.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/advanced/brachytherapy/Brachy.cc (Version 11.3.0) and /examples/advanced/brachytherapy/Brachy.cc (Version 8.1)


  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 //
                                                   >>  27 // $Id: Brachy.cc
                                                   >>  28 // GEANT4 tag $Name: geant4-08-01 $
                                                   >>  29 //
 26 // -------------------------------------------     30 // --------------------------------------------------------------
 27 //                 GEANT 4 - Brachytherapy exa     31 //                 GEANT 4 - Brachytherapy example
 28 // -------------------------------------------     32 // --------------------------------------------------------------
 29 //                                                 33 //
 30 // Code developed currently by:                <<  34 // Code developed by:
 31 //  S.Guatelli & D. Cutajar                    <<  35 // S. Agostinelli, F. Foppiano, S. Garelli , M. Tropeano, S.Guatelli
 32                                                    36 
 33 //                                                 37 //
 34 //    *******************************              38 //    *******************************
 35 //    *                             *              39 //    *                             *
 36 //    *    Brachy.cc                *              40 //    *    Brachy.cc                *
 37 //    *                             *              41 //    *                             *
 38 //    *******************************              42 //    *******************************
 39 //                                                 43 //
 40 //                                                 44 //
 41 #include "G4Types.hh"                          <<  45 
 42 #include "G4RunManagerFactory.hh"              <<  46 #include "G4RunManager.hh"
 43 #include "G4UImanager.hh"                          47 #include "G4UImanager.hh"
 44 #include "G4UIExecutive.hh"                    <<  48 #include "G4UIterminal.hh"
 45 #include "BrachyActionInitialization.hh"       <<  49 #include "G4UItcsh.hh"
                                                   >>  50 #include "BrachyFactoryIr.hh"
                                                   >>  51 #ifdef G4UI_USE_XM
                                                   >>  52 #include "G4UIXm.hh"
                                                   >>  53 #endif
                                                   >>  54 
                                                   >>  55 #ifdef G4VIS_USE
 46 #include "G4VisExecutive.hh"                       56 #include "G4VisExecutive.hh"
                                                   >>  57 #endif
                                                   >>  58 
                                                   >>  59 #include "BrachyEventAction.hh"
 47 #include "BrachyDetectorConstruction.hh"           60 #include "BrachyDetectorConstruction.hh"
 48 #include "BrachyPhysicsList.hh"                    61 #include "BrachyPhysicsList.hh"
                                                   >>  62 #include "BrachyPhantomSD.hh"
 49 #include "BrachyPrimaryGeneratorAction.hh"         63 #include "BrachyPrimaryGeneratorAction.hh"
 50 #include "G4SDManager.hh"                          64 #include "G4SDManager.hh"
 51 #include "Randomize.hh"                        <<  65 #include "BrachyRunAction.hh"
                                                   >>  66 #include "Randomize.hh"  
 52 #include "G4RunManager.hh"                         67 #include "G4RunManager.hh"
 53 #include "G4SDManager.hh"                          68 #include "G4SDManager.hh"
 54 #include "G4UImanager.hh"                          69 #include "G4UImanager.hh"
 55 #include "G4UImessenger.hh"                        70 #include "G4UImessenger.hh"
 56                                                    71 
 57 #include "G4ScoringManager.hh"                 <<  72 #ifdef G4ANALYSIS_USE
 58 #include "G4UIExecutive.hh"                    <<  73 #include "BrachyAnalysisManager.hh"
 59                                                <<  74 #endif
 60 #include "G4ScoringManager.hh"                 << 
 61 #include "BrachyUserScoreWriter.hh"            << 
 62                                                    75 
 63 int main(int argc ,char ** argv)                   76 int main(int argc ,char ** argv)
 64                                                    77 
 65 {                                                  78 {
 66  // Construct the default run manager          <<  79   G4RunManager* pRunManager = new G4RunManager;
 67  //                                            << 
 68   auto* pRunManager = G4RunManagerFactory::Cre << 
 69   G4int nThreads = 4;                          << 
 70   pRunManager->SetNumberOfThreads(nThreads);   << 
 71                                                << 
 72   G4cout << "***********************" << G4end << 
 73   G4cout << "*** Seed: " << G4Random::getTheSe << 
 74   G4cout << "***********************" << G4end << 
 75                                                << 
 76   // Access to the Scoring Manager pointer     << 
 77                                                << 
 78   G4ScoringManager* scoringManager = G4Scoring << 
 79                                                    80 
 80   // Overwrite the default output file with us <<  81   G4String sensitiveDetectorName = "Phantom";
 81   scoringManager->SetScoreWriter(new BrachyUse << 
 82                                                    82 
 83   // Initialize the physics component              83   // Initialize the physics component
 84   pRunManager -> SetUserInitialization(new Bra     84   pRunManager -> SetUserInitialization(new BrachyPhysicsList);
 85                                                    85 
 86   // Initialize the detector component             86   // Initialize the detector component
 87   auto pDetectorConstruction = new  BrachyDete <<  87   BrachyDetectorConstruction  *pDetectorConstruction = new  BrachyDetectorConstruction(sensitiveDetectorName);
 88   pRunManager -> SetUserInitialization(pDetect     88   pRunManager -> SetUserInitialization(pDetectorConstruction);
 89                                                    89 
 90   // User action initialization                <<  90   // Initialize the primary particles
                                                   >>  91   BrachyPrimaryGeneratorAction* primary = new BrachyPrimaryGeneratorAction();
                                                   >>  92   pRunManager -> SetUserAction(primary);
                                                   >>  93 
                                                   >>  94   // Initialize Optional User Action
                                                   >>  95   BrachyEventAction *pEventAction = new BrachyEventAction();
                                                   >>  96   pRunManager -> SetUserAction(pEventAction );
 91                                                    97 
 92  auto actions = new BrachyActionInitialization <<  98   BrachyRunAction *pRunAction = new BrachyRunAction();
 93   pRunManager->SetUserInitialization(actions); <<  99   pRunManager -> SetUserAction(pRunAction);
 94                                                << 
 95   //Initialize G4 kernel                       << 
 96  // pRunManager -> Initialize();               << 
 97                                                   100 
                                                   >> 101   //// Initialize the Visualization component 
                                                   >> 102 #ifdef G4VIS_USE
 98   // Visualization manager                        103   // Visualization manager
 99   G4VisManager* visManager = new G4VisExecutiv    104   G4VisManager* visManager = new G4VisExecutive;
100   visManager->Initialize();                       105   visManager->Initialize();
101                                                << 106 #endif
102   // get the pointer to the User Interface man << 107   
103   G4UImanager* UImanager = G4UImanager::GetUIp << 108   // Output environment variables concerning the analysis:
                                                   >> 109 #ifdef G4ANALYSIS_USE
                                                   >> 110   G4cout << G4endl << G4endl << G4endl 
                                                   >> 111    << " User Environment " << G4endl
                                                   >> 112    << " Using AIDA 3.2.1 analysis " << G4endl;
                                                   >> 113 # else
                                                   >> 114   G4cout << G4endl << G4endl << G4endl 
                                                   >> 115    << " User Environment " << G4endl
                                                   >> 116    << " G4ANALYSIS_USE environment variable not set, NO ANALYSIS " 
                                                   >> 117    << G4endl;
                                                   >> 118 #endif
                                                   >> 119   
                                                   >> 120   // Initialize the interactive session 
                                                   >> 121   G4UIsession* session = 0;
104   if (argc == 1)   // Define UI session for in    122   if (argc == 1)   // Define UI session for interactive mode.
105     {                                             123     {
106       auto ui = new G4UIExecutive(argc, argv); << 124       session = new G4UIterminal();
                                                   >> 125     }
                                                   >> 126 
                                                   >> 127   //Initialize G4 kernel
                                                   >> 128   pRunManager -> Initialize();
                                                   >> 129 
                                                   >> 130   // get the pointer to the User Interface manager 
                                                   >> 131   G4UImanager* UI = G4UImanager::GetUIpointer();  
                                                   >> 132   if (session)   // Define UI session for interactive mode.
                                                   >> 133     { 
107       G4cout << " UI session starts ..." << G4    134       G4cout << " UI session starts ..." << G4endl;
108       UImanager -> ApplyCommand("/control/exec << 135       UI -> ApplyCommand("/control/execute VisualisationMacro.mac");    
109       ui -> SessionStart();                    << 136       session -> SessionStart();
110       delete ui;                               << 137       delete session;
111     }                                             138     }
112   else           // Batch mode                    139   else           // Batch mode
113     {                                          << 140     { 
114       G4String command = "/control/execute ";     141       G4String command = "/control/execute ";
115       G4String fileName = argv[1];                142       G4String fileName = argv[1];
116       UImanager -> ApplyCommand(command+fileNa << 143       UI -> ApplyCommand(command+fileName);
117     }                                          << 144     }  
118                                                   145 
                                                   >> 146   // Close the output file containing histograms and n-tuple with the 
                                                   >> 147   // results of the simulation.
                                                   >> 148   // Finish the analysis.
                                                   >> 149 #ifdef G4ANALYSIS_USE
                                                   >> 150   BrachyAnalysisManager* analysis = BrachyAnalysisManager::getInstance();
                                                   >> 151   analysis -> finish();
                                                   >> 152 #endif
                                                   >> 153   
119   // Job termination                              154   // Job termination
120  // Close the root file                        << 155 #ifdef G4VIS_USE
121                                                << 
122   delete visManager;                              156   delete visManager;
                                                   >> 157 #endif
                                                   >> 158 
123   delete pRunManager;                             159   delete pRunManager;
124                                                   160 
125   return 0;                                       161   return 0;
126 }                                                 162 }
127                                                   163