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 9.1.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 //                                                 26 //
 27 /// \file field/field01/field01.cc             <<  27 // $Id: field01.cc,v 1.8 2006/06/29 17:15:22 gunter Exp $
 28 /// \brief Main program of the field/field01 e <<  28 // GEANT4 tag $Name: geant4-09-01-patch-03 $
 29 //                                                 29 //
 30 //....oooOO0OOooo........oooOO0OOooo........oo <<  30 // 
 31 //....oooOO0OOooo........oooOO0OOooo........oo <<  31 // --------------------------------------------------------------
 32                                                <<  32 //      GEANT 4 - TestF01 
 33 #include "F01ActionInitialization.hh"          <<  33 //
 34 #include "F01DetectorConstruction.hh"          <<  34 // --------------------------------------------------------------
 35 #include "F01RunAction.hh"                     <<  35 // Comments
 36 #include "F01SteppingVerbose.hh"               <<  36 //     
                                                   >>  37 //   
                                                   >>  38 // --------------------------------------------------------------
 37                                                    39 
 38 #include "G4RunManagerFactory.hh"              <<  40 #include "G4RunManager.hh"
 39 #include "G4Types.hh"                          << 
 40 #include "G4UImanager.hh"                          41 #include "G4UImanager.hh"
 41                                                <<  42 #include "G4UIterminal.hh"
 42 #include "FTFP_BERT.hh"                        << 
 43                                                << 
 44 #include "G4EmParameters.hh"                   << 
 45 #include "G4HadronicParameters.hh"             << 
 46 #include "G4PhysicsListHelper.hh"              << 
 47 #include "G4StepLimiterPhysics.hh"             << 
 48 #include "G4UIExecutive.hh"                    << 
 49 #include "G4VisExecutive.hh"                   << 
 50 #include "Randomize.hh"                            43 #include "Randomize.hh"
 51                                                    44 
 52 // For Printing statistic from Transporation p <<  45 #ifdef G4VIS_USE
 53 #include "G4CoupledTransportation.hh"          <<  46 #include "G4VisExecutive.hh"
 54 #include "G4Electron.hh"                       <<  47 #endif
 55 #include "G4Transportation.hh"                 << 
 56 #include "G4TransportationParameters.hh"       << 
 57                                                    48 
 58 //....oooOO0OOooo........oooOO0OOooo........oo <<  49 #include "F01DetectorConstruction.hh"
                                                   >>  50 // #include "F01FieldSetup.hh"
                                                   >>  51 #include "F01PhysicsList.hh"
                                                   >>  52 #include "F01PrimaryGeneratorAction.hh"
                                                   >>  53 #include "F01RunAction.hh"
                                                   >>  54 #include "F01EventAction.hh"
                                                   >>  55 #include "F01SteppingAction.hh"
                                                   >>  56 #include "F01SteppingVerbose.hh"
 59                                                    57 
 60 int main(int argc, char** argv)                <<  58 int main(int argc,char** argv) 
 61 {                                                  59 {
 62   // Instantiate G4UIExecutive if there are no << 
 63   //                                           << 
 64   G4UIExecutive* ui = nullptr;                 << 
 65   if (argc == 1) {                             << 
 66     ui = new G4UIExecutive(argc, argv);        << 
 67   }                                            << 
 68                                                    60 
 69   // Setting the application-specific Stepping <<  61   //choose the Random engine
 70   //                                           << 
 71   auto verbosity = new F01SteppingVerbose;     << 
 72                                                    62 
                                                   >>  63   CLHEP::HepRandom::setTheEngine(new CLHEP::RanecuEngine);
                                                   >>  64   
                                                   >>  65   //my Verbose output class
                                                   >>  66 
                                                   >>  67   G4VSteppingVerbose::SetInstance(new F01SteppingVerbose);
                                                   >>  68   
 73   // Construct the default run manager             69   // Construct the default run manager
 74   auto runManager = G4RunManagerFactory::Creat << 
 75                                                    70 
 76   // G4TransportationWithMscType: fDisabled, f <<  71   G4RunManager * runManager = new G4RunManager;
 77   // G4EmParameters::Instance()->SetTransporta <<  72 
                                                   >>  73   // Create class that constructs & messages the magnetic field
                                                   >>  74   //  *** F01FieldSetup* field = new F01FieldSetup() ;
                                                   >>  75   // --> Now done in detector construction
 78                                                    76 
 79   // Set mandatory initialization classes          77   // Set mandatory initialization classes
 80   //                                           << 
 81   // Detector construction                     << 
 82   auto detector = new F01DetectorConstruction( << 
 83   // detector->SetUseFSALstepper();  // Uncomm << 
 84                                                    78 
                                                   >>  79   F01DetectorConstruction* detector;
                                                   >>  80   detector = new F01DetectorConstruction;
 85   runManager->SetUserInitialization(detector);     81   runManager->SetUserInitialization(detector);
                                                   >>  82   runManager->SetUserInitialization(new F01PhysicsList(detector));
                                                   >>  83   
                                                   >>  84 #ifdef G4VIS_USE
 86                                                    85 
 87   // Configure the use of low thresholds for l <<  86   // visualization manager
 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                              << 
109   G4VModularPhysicsList* physicsList = new FTF << 
110   physicsList->RegisterPhysics(new G4StepLimit << 
111   runManager->SetUserInitialization(physicsLis << 
112                                                << 
113   // User action initialization                << 
114   runManager->SetUserInitialization(new F01Act << 
115                                                << 
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                                                    87 
141     G4cout << "field01: Using F01RunAction to  <<  88   G4VisManager* visManager = new G4VisExecutive;
142     runManager->SetUserAction(runAction);      <<  89   visManager->Initialize();
143   }                                            << 
144                                                    90 
145   // Note: this mechanism overwrites the thres <<  91 #endif 
146   //       the call to UseLowLooperThresholds( <<  92  
                                                   >>  93   // Set user action classes
147                                                    94 
148   // Suppress large verbosity from EM & hadron <<  95   runManager->SetUserAction(new F01PrimaryGeneratorAction(detector));
149   G4EmParameters::Instance()->SetVerbose(0);   << 
150   G4HadronicParameters::Instance()->SetVerbose << 
151                                                    96 
152   // Initialize G4 kernel                      <<  97   F01RunAction* runAction = new F01RunAction;
153   //                                           << 
154   runManager->Initialize();                    << 
155                                                    98 
156   // Initialize visualization                  <<  99   runManager->SetUserAction(runAction);
157   //                                           << 
158   // G4VisExecutive can take a verbosity argum << 
159   G4VisManager* visManager = new G4VisExecutiv << 
160   visManager->Initialize();                    << 
161                                                   100 
162   // Get the pointer to the User Interface man << 101   F01EventAction* eventAction = new F01EventAction(runAction);
163   //                                           << 102 
164   G4UImanager* UImanager = G4UImanager::GetUIp << 103   runManager->SetUserAction(eventAction);
165                                                << 104 
166   if (!ui)  // batch mode                      << 105   F01SteppingAction* steppingAction = new F01SteppingAction();
167   {                                            << 106   runManager->SetUserAction(steppingAction);
168     G4String command = "/control/execute ";    << 107   
169     G4String fileName = argv[1];               << 108   // Initialize G4 kernel, physics tables ...
170     UImanager->ApplyCommand(command + fileName << 109 
                                                   >> 110   runManager->Initialize();
                                                   >> 111     
                                                   >> 112   // Get the pointer to the User Interface manager 
                                                   >> 113 
                                                   >> 114   G4UImanager* UI = G4UImanager::GetUIpointer();  
                                                   >> 115  
                                                   >> 116   if (argc==1)   // Define UI terminal for interactive mode  
                                                   >> 117   { 
                                                   >> 118      G4UIsession * session = new G4UIterminal;
                                                   >> 119      session->SessionStart();
                                                   >> 120      delete session;
171   }                                               121   }
172   else {  // interactive mode : define UI sess << 122   else           // Batch mode
173     UImanager->ApplyCommand("/control/execute  << 123   { 
174     if (ui->IsGUI()) UImanager->ApplyCommand(" << 124      G4String command = "/control/execute ";
175     ui->SessionStart();                        << 125      G4String fileName = argv[1];
176     delete ui;                                 << 126      UI->ApplyCommand(command+fileName);
177   }                                               127   }
                                                   >> 128     
                                                   >> 129   // job termination
178                                                   130 
179   // Statistics of tracks killed by G4Transpor << 131 #ifdef G4VIS_USE
180   //  printed in the RunAction's EndOfEvent ac << 
181   // ( Eventually a summary could be provided  << 
182                                                << 
183   delete verbosity;                            << 
184   delete visManager;                              132   delete visManager;
                                                   >> 133 #endif
                                                   >> 134   // delete field;
185   delete runManager;                              135   delete runManager;
186                                                   136 
187   return 0;                                       137   return 0;
188 }                                                 138 }
189                                                   139 
190 //....oooOO0OOooo........oooOO0OOooo........oo << 
191                                                   140