Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/transportation/include/G4Transportation.icc

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 /processes/transportation/include/G4Transportation.icc (Version 11.3.0) and /processes/transportation/include/G4Transportation.icc (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 // $Id: G4Transportation.icc,v 1.11 2007/11/08 17:59:52 japost Exp $
                                                   >>  28 // GEANT4 tag $Name: geant4-09-01-patch-03 $
 27 //                                                 29 //
                                                   >>  30 
 28 //                                                 31 // 
 29 //                                                 32 //
 30 //        Inline function implementation.          33 //        Inline function implementation.
 31                                                <<  34 //
 32 // ===========================================     35 // =======================================================================
 33 // Created:  9 June 1998, J. Apostolakis           36 // Created:  9 June 1998, J. Apostolakis
 34 // ===========================================     37 // =======================================================================
 35                                                <<  38 //
 36 #include <CLHEP/Units/SystemOfUnits.h>         << 
 37 #include "G4TransportationLogger.hh"           << 
 38                                                    39 
 39 inline void                                        40 inline void
 40 G4Transportation::SetPropagatorInField( G4Prop     41 G4Transportation::SetPropagatorInField( G4PropagatorInField* pFieldPropagator)
 41 {                                                  42 {
 42    fFieldPropagator= pFieldPropagator;             43    fFieldPropagator= pFieldPropagator;
 43 }                                                  44 }
 44                                                    45 
 45 inline G4PropagatorInField* G4Transportation::     46 inline G4PropagatorInField* G4Transportation::GetPropagatorInField()
 46 {                                                  47 {
 47    return fFieldPropagator;                        48    return fFieldPropagator;
 48 }                                                  49 }
 49                                                    50 
                                                   >>  51 inline G4bool G4Transportation::DoesGlobalFieldExist()
                                                   >>  52 {
                                                   >>  53   G4TransportationManager* transportMgr;
                                                   >>  54   transportMgr= G4TransportationManager::GetTransportationManager();
                                                   >>  55 
                                                   >>  56   // fFieldExists= transportMgr->GetFieldManager()->DoesFieldExist();
                                                   >>  57   // return fFieldExists;
                                                   >>  58   return transportMgr->GetFieldManager()->DoesFieldExist();
                                                   >>  59 }
                                                   >>  60 
                                                   >>  61 inline void G4Transportation::SetVerboseLevel( G4int verboseLev )
                                                   >>  62 {
                                                   >>  63   fVerboseLevel= verboseLev;
                                                   >>  64 }
                                                   >>  65 
                                                   >>  66 inline G4int G4Transportation::GetVerboseLevel( ) const
                                                   >>  67 {
                                                   >>  68   return fVerboseLevel;
                                                   >>  69 }
                                                   >>  70 
 50 inline G4double G4Transportation::GetThreshold     71 inline G4double G4Transportation::GetThresholdWarningEnergy() const
 51 {                                                  72 {
 52   return fThreshold_Warning_Energy;                73   return fThreshold_Warning_Energy;
 53 }                                                  74 }
 54                                                    75  
 55 inline G4double G4Transportation::GetThreshold     76 inline G4double G4Transportation::GetThresholdImportantEnergy() const
 56 {                                                  77 { 
 57   return fThreshold_Important_Energy;              78   return fThreshold_Important_Energy;
 58 }                                                  79 } 
 59                                                    80 
 60 inline G4int G4Transportation::GetThresholdTri     81 inline G4int G4Transportation::GetThresholdTrials() const
 61 {                                                  82 {
 62   return fThresholdTrials;                         83   return fThresholdTrials;
 63 }                                                  84 }
 64                                                    85 
 65 inline void G4Transportation::SetThresholdWarn     86 inline void G4Transportation::SetThresholdWarningEnergy( G4double newEnWarn )
 66 {                                                  87 {
 67   fThreshold_Warning_Energy= newEnWarn;            88   fThreshold_Warning_Energy= newEnWarn;
 68   if( fpLogger ) {                             << 
 69      fpLogger->SetThresholdWarningEnergy( newE << 
 70   } else {                                     << 
 71      ReportMissingLogger("SetThresholdWarningE << 
 72   }                                            << 
 73   // Update the information for correct report << 
 74 }                                                  89 }
 75                                                    90 
 76 inline void G4Transportation::SetThresholdImpo     91 inline void G4Transportation::SetThresholdImportantEnergy( G4double newEnImp )
 77 {                                                  92 {
 78   fThreshold_Important_Energy = newEnImp;      <<  93   fThreshold_Important_Energy = newEnImp; 
 79   if( fpLogger ) {                             << 
 80      fpLogger->SetThresholdImportantEnergy( ne << 
 81     // Update the information for correct repo << 
 82   } else {                                     << 
 83      ReportMissingLogger("SetThresholdImportan << 
 84   }                                            << 
 85 }                                                  94 }
 86                                                    95 
 87 inline void G4Transportation::SetThresholdTria     96 inline void G4Transportation::SetThresholdTrials(G4int newMaxTrials )
 88 {                                                  97 {
 89   fThresholdTrials = newMaxTrials;             <<  98   fThresholdTrials = newMaxTrials; 
 90   if( fpLogger ) {                             << 
 91      fpLogger->SetThresholdTrials(newMaxTrials << 
 92      // Update the information for correct rep << 
 93   } else {                                     << 
 94      ReportMissingLogger("SetThresholdTrials") << 
 95   }                                            << 
 96 }                                                  99 }
 97                                                   100 
                                                   >> 101      // Get/Set parameters for killing loopers: 
                                                   >> 102      //   Above 'important' energy a 'looping' particle in field will 
                                                   >> 103      //   *NOT* be abandoned, except after fThresholdTrials attempts.
                                                   >> 104      // Below Warning energy, no verbosity for looping particles is issued
                                                   >> 105 
 98 inline G4double G4Transportation::GetMaxEnergy    106 inline G4double G4Transportation::GetMaxEnergyKilled() const
 99 {                                                 107 {
100   return fMaxEnergyKilled;                        108   return fMaxEnergyKilled; 
101 }                                                 109 }
102                                                   110 
103 inline G4double G4Transportation::GetSumEnergy    111 inline G4double G4Transportation::GetSumEnergyKilled() const
104 {                                                 112 {
105   return fSumEnergyKilled;                        113   return fSumEnergyKilled;
106 }                                                 114 }
107                                                   115 
108 inline void G4Transportation::ResetKilledStati    116 inline void G4Transportation::ResetKilledStatistics(G4int report)
109 {                                                 117 {
110   // Statistics for tracks killed (currently d << 118   if( report ) { 
111                                                << 119     G4cout << " G4Transportation: Statistics for looping particles " << G4endl;
112   if( report )                                 << 120     G4cout << "   Sum of energy of loopers killed: " <<  fSumEnergyKilled << G4endl;
113   {                                            << 121     G4cout << "   Max energy of loopers killed: " <<  fMaxEnergyKilled << G4endl;
114     G4cout << " G4Transportation: Statistics f << 
115            << G4endl;                          << 
116     G4cout << "   Sum of energy of loopers kil << 
117            << G4endl;                          << 
118     G4cout << "   Max energy of loopers killed << 
119            << G4endl;                          << 
120   }                                               122   } 
121                                                   123 
122   fSumEnergyKilled= 0;                            124   fSumEnergyKilled= 0;
123   fMaxEnergyKilled= -1.0*CLHEP::GeV;           << 125   fMaxEnergyKilled= -1.0*GeV;
124 }                                                 126 }
                                                   >> 127      // Statistics for tracks killed (currently due to looping in field)
125                                                   128 
126 inline void                                    << 129 inline void G4Transportation::EnableShortStepOptimisation(G4bool optimiseShortStep)
127 G4Transportation::EnableShortStepOptimisation( << 
128 {                                                 130 { 
129   fShortStepOptimisation=optimiseShortStep;       131   fShortStepOptimisation=optimiseShortStep;
130 }                                                 132 }
131                                                << 
132 inline void                                    << 
133 G4Transportation::PushThresholdsToLogger()     << 
134 {                                              << 
135   if( fpLogger )                               << 
136      fpLogger->SetThresholds( GetThresholdWarn << 
137                               GetThresholdImpo << 
138                               GetThresholdTria << 
139   else                                         << 
140      ReportMissingLogger("PushThresholdsToLogg << 
141 }                                              << 
142                                                << 
143                                                   133