Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/field/field01/field01.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/field/field01/field01.cc (Version 11.3.0) and /examples/extended/field/field01/field01.cc (Version 10.3.p1)


  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 // $Id: field01.cc 100688 2016-10-31 11:21:51Z gcosmo $
 26 //                                                 27 //
 27 /// \file field/field01/field01.cc                 28 /// \file field/field01/field01.cc
 28 /// \brief Main program of the field/field01 e     29 /// \brief Main program of the field/field01 example
 29 //                                                 30 //
                                                   >>  31 //
                                                   >>  32 //
                                                   >>  33 //
 30 //....oooOO0OOooo........oooOO0OOooo........oo     34 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 31 //....oooOO0OOooo........oooOO0OOooo........oo     35 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 32                                                    36 
 33 #include "F01ActionInitialization.hh"          <<  37 #ifdef G4MULTITHREADED
 34 #include "F01DetectorConstruction.hh"          <<  38 #include "G4MTRunManager.hh"
 35 #include "F01RunAction.hh"                     <<  39 #else
 36 #include "F01SteppingVerbose.hh"                   40 #include "F01SteppingVerbose.hh"
                                                   >>  41 #include "G4RunManager.hh"
                                                   >>  42 #endif
 37                                                    43 
 38 #include "G4RunManagerFactory.hh"              <<  44 #include "F01DetectorConstruction.hh"
 39 #include "G4Types.hh"                          <<  45 #include "F01ActionInitialization.hh"
 40 #include "G4UImanager.hh"                      << 
 41                                                    46 
                                                   >>  47 #include "G4UImanager.hh"
 42 #include "FTFP_BERT.hh"                            48 #include "FTFP_BERT.hh"
 43                                                << 
 44 #include "G4EmParameters.hh"                   << 
 45 #include "G4HadronicParameters.hh"             << 
 46 #include "G4PhysicsListHelper.hh"              << 
 47 #include "G4StepLimiterPhysics.hh"                 49 #include "G4StepLimiterPhysics.hh"
 48 #include "G4UIExecutive.hh"                    << 
 49 #include "G4VisExecutive.hh"                   << 
 50 #include "Randomize.hh"                            50 #include "Randomize.hh"
 51                                                    51 
 52 // For Printing statistic from Transporation p <<  52 #ifdef G4VIS_USE
 53 #include "G4CoupledTransportation.hh"          <<  53 #include "G4VisExecutive.hh"
 54 #include "G4Electron.hh"                       <<  54 #endif
 55 #include "G4Transportation.hh"                 <<  55 
 56 #include "G4TransportationParameters.hh"       <<  56 #ifdef G4UI_USE
                                                   >>  57 #include "G4UIExecutive.hh"
                                                   >>  58 #endif
 57                                                    59 
 58 //....oooOO0OOooo........oooOO0OOooo........oo     60 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 59                                                    61 
 60 int main(int argc, char** argv)                <<  62 int main(int argc,char** argv)
 61 {                                                  63 {
 62   // Instantiate G4UIExecutive if there are no <<  64   // Choose the Random engine
 63   //                                               65   //
 64   G4UIExecutive* ui = nullptr;                 <<  66   G4Random::setTheEngine(new CLHEP::RanecuEngine);
 65   if (argc == 1) {                             << 
 66     ui = new G4UIExecutive(argc, argv);        << 
 67   }                                            << 
 68                                                << 
 69   // Setting the application-specific Stepping << 
 70   //                                           << 
 71   auto verbosity = new F01SteppingVerbose;     << 
 72                                                    67 
 73   // Construct the default run manager             68   // Construct the default run manager
 74   auto runManager = G4RunManagerFactory::Creat <<  69   //
 75                                                <<  70 #ifdef G4MULTITHREADED
 76   // G4TransportationWithMscType: fDisabled, f <<  71   G4MTRunManager * runManager = new G4MTRunManager;
 77   // G4EmParameters::Instance()->SetTransporta <<  72 #else
                                                   >>  73   G4VSteppingVerbose::SetInstance(new F01SteppingVerbose);
                                                   >>  74   G4RunManager * runManager = new G4RunManager;
                                                   >>  75 #endif
 78                                                    76 
 79   // Set mandatory initialization classes          77   // Set mandatory initialization classes
 80   //                                               78   //
 81   // Detector construction                         79   // Detector construction
 82   auto detector = new F01DetectorConstruction( <<  80   F01DetectorConstruction* detector = new F01DetectorConstruction();
 83   // detector->SetUseFSALstepper();  // Uncomm << 
 84                                                << 
 85   runManager->SetUserInitialization(detector);     81   runManager->SetUserInitialization(detector);
 86                                                << 
 87   // Configure the use of low thresholds for l << 
 88   //  ( appropriate for typical applications u << 
 89   // auto plHelper = G4PhysicsListHelper::GetP << 
 90   // plHelper->UseLowLooperThresholds();       << 
 91   // plHelper->UseHighLooperThresholds();      << 
 92   // Request a set of pre-selected values of t << 
 93   //    particles:                             << 
 94   //       - High for collider HEP application << 
 95   //       - Low  for 'low-E' applications, me << 
 96   // Note: If helper is used select low or hig << 
 97   //       values from TransportationParameter << 
 98                                                << 
 99   // They are currently applied in the followi << 
100   // 1. Transportation Parameters - fine grain << 
101   // 2. Physics List Helper       - impose a f << 
102   // 3. Run Action (F01RunAction) - revise val << 
103   // 4. Tracking Action           - could revi << 
104   //     Note that this also could customise b << 
105   //     to mu-/mu+ , e-/e+ vs others)         << 
106   // If multiple are present, later methods ov << 
107                                                << 
108   // Physics list                                  82   // Physics list
109   G4VModularPhysicsList* physicsList = new FTF     83   G4VModularPhysicsList* physicsList = new FTFP_BERT;
110   physicsList->RegisterPhysics(new G4StepLimit     84   physicsList->RegisterPhysics(new G4StepLimiterPhysics());
111   runManager->SetUserInitialization(physicsLis     85   runManager->SetUserInitialization(physicsList);
112                                                << 
113   // User action initialization                    86   // User action initialization
114   runManager->SetUserInitialization(new F01Act     87   runManager->SetUserInitialization(new F01ActionInitialization(detector));
115                                                    88 
116   G4double warningE = 10.0 * CLHEP::keV;       << 
117   G4double importantE = 0.1 * CLHEP::MeV;      << 
118   G4int numTrials = 30;                        << 
119                                                << 
120   G4bool useTransportParams = true;  // Use th << 
121                                                << 
122   if (useTransportParams) {                    << 
123     auto transportParams = G4TransportationPar << 
124     transportParams->SetWarningEnergy(warningE << 
125     transportParams->SetImportantEnergy(import << 
126     transportParams->SetNumberOfTrials(numTria << 
127     G4cout << "field01: Using G4Transportation << 
128   }                                            << 
129   else {                                       << 
130     // Fine grained control of thresholds for  << 
131     auto runAction = new F01RunAction();       << 
132     runAction->SetWarningEnergy(warningE);     << 
133     // Looping particles with E < 10 keV will  << 
134     //   with warning.                         << 
135     // Looping particles with E > 10 keV will  << 
136     runAction->SetImportantEnergy(importantE); << 
137     runAction->SetNumberOfTrials(numTrials);   << 
138     // Looping particles with E > 0.1 MeV will << 
139     //  30 'tracking' steps, and only be kille << 
140                                                << 
141     G4cout << "field01: Using F01RunAction to  << 
142     runManager->SetUserAction(runAction);      << 
143   }                                            << 
144                                                << 
145   // Note: this mechanism overwrites the thres << 
146   //       the call to UseLowLooperThresholds( << 
147                                                << 
148   // Suppress large verbosity from EM & hadron << 
149   G4EmParameters::Instance()->SetVerbose(0);   << 
150   G4HadronicParameters::Instance()->SetVerbose << 
151                                                << 
152   // Initialize G4 kernel                          89   // Initialize G4 kernel
153   //                                               90   //
154   runManager->Initialize();                        91   runManager->Initialize();
155                                                    92 
                                                   >>  93 #ifdef G4VIS_USE
156   // Initialize visualization                      94   // Initialize visualization
157   //                                               95   //
158   // G4VisExecutive can take a verbosity argum <<  96   G4VisManager* visManager = new G4VisExecutive;
159   G4VisManager* visManager = new G4VisExecutiv <<  97   // G4VisExecutive can take a verbosity argument - see /vis/verbose guidance.
                                                   >>  98   // G4VisManager* visManager = new G4VisExecutive("Quiet");
160   visManager->Initialize();                        99   visManager->Initialize();
                                                   >> 100 #endif
161                                                   101 
162   // Get the pointer to the User Interface man    102   // Get the pointer to the User Interface manager
163   //                                              103   //
164   G4UImanager* UImanager = G4UImanager::GetUIp    104   G4UImanager* UImanager = G4UImanager::GetUIpointer();
165                                                   105 
166   if (!ui)  // batch mode                      << 106   if (argc!=1)   // batch mode
167   {                                            << 107     {
168     G4String command = "/control/execute ";    << 108       G4String command = "/control/execute ";
169     G4String fileName = argv[1];               << 109       G4String fileName = argv[1];
170     UImanager->ApplyCommand(command + fileName << 110       UImanager->ApplyCommand(command+fileName);
171   }                                            << 111     }
172   else {  // interactive mode : define UI sess << 112   else
173     UImanager->ApplyCommand("/control/execute  << 113     {  // interactive mode : define UI session
174     if (ui->IsGUI()) UImanager->ApplyCommand(" << 114 #ifdef G4UI_USE
175     ui->SessionStart();                        << 115      G4UIExecutive* ui = new G4UIExecutive(argc, argv);
176     delete ui;                                 << 116      if (ui->IsGUI())
177   }                                            << 117         UImanager->ApplyCommand("/control/execute gui.mac");
178                                                << 118      ui->SessionStart();
179   // Statistics of tracks killed by G4Transpor << 119      delete ui;
180   //  printed in the RunAction's EndOfEvent ac << 120 #endif
181   // ( Eventually a summary could be provided  << 121     }
                                                   >> 122 
                                                   >> 123   // Job termination
                                                   >> 124   // Free the store: user actions, physics_list and detector_description are
                                                   >> 125   //                 owned and deleted by the run manager, so they should not
                                                   >> 126   //                 be deleted in the main() program !
182                                                   127 
183   delete verbosity;                            << 128 #ifdef G4VIS_USE
184   delete visManager;                              129   delete visManager;
                                                   >> 130 #endif
185   delete runManager;                              131   delete runManager;
186                                                   132 
187   return 0;                                       133   return 0;
188 }                                                 134 }
189                                                   135 
190 //....oooOO0OOooo........oooOO0OOooo........oo    136 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
191                                                   137