Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/management/include/G4VContinuousDiscreteProcess.hh

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/management/include/G4VContinuousDiscreteProcess.hh (Version 11.3.0) and /processes/management/include/G4VContinuousDiscreteProcess.hh (Version 8.1.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 // G4VContinuousDiscreteProcess                << 
 27 //                                                 26 //
 28 // Class description:                          <<  27 // $Id: G4VContinuousDiscreteProcess.hh,v 1.8 2006/06/29 21:07:44 gunter Exp $
                                                   >>  28 // GEANT4 tag $Name: geant4-08-01-patch-01 $
                                                   >>  29 //
                                                   >>  30 // 
                                                   >>  31 // ------------------------------------------------------------
                                                   >>  32 //  GEANT 4 class header file 
                                                   >>  33 //
                                                   >>  34 //  Abstract class which defines the public behavior of
 29 //                                                 35 // 
 30 // Abstract class which defines the public beh <<  36 // Class Description 
 31 // discrete physics interactions.              <<  37 //  discrete physics interactions.
                                                   >>  38 //
                                                   >>  39 // ------------------------------------------------------------
                                                   >>  40 //   New Physics scheme            8  Mar. 1997  H.Kurahige
                                                   >>  41 // ------------------------------------------------------------
                                                   >>  42 //   fix bugs in GetGPILSelection() 24 Jan. 1998 H.Kurashige
                                                   >>  43 //   modified for new ParticleChange 12 Mar. 1998  H.Kurashige
                                                   >>  44 //   Fixed a bug in PostStepGetPhysicalInteractionLength  
                                                   >>  45 //                                15 Apr. 2002 H.Kurashige 
                                                   >>  46 //
 32                                                    47 
 33 // Authors:                                    <<  48 #ifndef G4VContinuousDiscreteProcess_h
 34 // - 2 December 1995, G.Cosmo - First implemen <<  49 #define G4VContinuousDiscreteProcess_h 1
 35 // - 8 January 1997, H.Kurashige - New Physics << 
 36 // ------------------------------------------- << 
 37 #ifndef G4VContinuousDiscreteProcess_hh        << 
 38 #define G4VContinuousDiscreteProcess_hh 1      << 
 39                                                    50 
 40 #include "globals.hh"                              51 #include "globals.hh"
 41 #include "G4ios.hh"                                52 #include "G4ios.hh"
 42                                                    53 
 43 #include "G4VProcess.hh"                           54 #include "G4VProcess.hh"
 44                                                    55 
 45 class G4VContinuousDiscreteProcess : public G4     56 class G4VContinuousDiscreteProcess : public G4VProcess 
 46 {                                                  57 {
                                                   >>  58   //  Abstract class which defines the public behavior of
                                                   >>  59   //  discrete physics interactions.
 47   public:                                          60   public:     
 48                                                    61 
 49     G4VContinuousDiscreteProcess(const G4Strin <<  62       G4VContinuousDiscreteProcess(const G4String& ,
 50                                  G4ProcessType <<  63            G4ProcessType   aType = fNotDefined );
 51     G4VContinuousDiscreteProcess(G4VContinuous <<  64       G4VContinuousDiscreteProcess(G4VContinuousDiscreteProcess &);
 52                                                    65 
 53     virtual ~G4VContinuousDiscreteProcess();   <<  66       virtual ~G4VContinuousDiscreteProcess();
 54                                                    67 
 55     G4VContinuousDiscreteProcess& operator=(co <<  68   public :// with description
 56                                                <<  69       virtual G4double PostStepGetPhysicalInteractionLength(
 57     virtual G4double PostStepGetPhysicalIntera << 
 58                              const G4Track& tr     70                              const G4Track& track,
 59                              G4double   previo <<  71            G4double   previousStepSize,
 60                              G4ForceCondition* <<  72            G4ForceCondition* condition
 61                             );                 <<  73           );
 62                                                <<  74 
 63     virtual G4VParticleChange* PostStepDoIt(   <<  75       virtual G4VParticleChange* PostStepDoIt(
 64                              const G4Track& ,  <<  76            const G4Track& ,
 65                              const G4Step&     <<  77            const G4Step& 
 66                             );                 <<  78           );
 67                                                    79 
 68     virtual G4double AlongStepGetPhysicalInter <<  80       virtual G4double AlongStepGetPhysicalInteractionLength(
 69                              const G4Track&,       81                              const G4Track&,
 70                              G4double  previou     82                              G4double  previousStepSize,
 71                              G4double  current     83                              G4double  currentMinimumStep,
 72                              G4double& current <<  84            G4double& currentSafety,
 73                              G4GPILSelection*      85                              G4GPILSelection* selection
 74                             );                 <<  86           );
 75                                                    87 
 76     virtual G4VParticleChange* AlongStepDoIt(  <<  88       virtual G4VParticleChange* AlongStepDoIt(
 77                              const G4Track& ,  <<  89            const G4Track& ,
 78                              const G4Step&     <<  90            const G4Step& 
 79                             );                 <<  91           );
 80                                                    92  
 81     virtual G4double AtRestGetPhysicalInteract <<  93      //  no operation in  AtRestDoIt
                                                   >>  94       virtual G4double AtRestGetPhysicalInteractionLength(
 82                              const G4Track& ,      95                              const G4Track& ,
 83                              G4ForceCondition* <<  96            G4ForceCondition* 
 84                             ) { return -1.0; } <<  97           ) { return -1.0; };
 85       // No operation in AtRestDoIt            << 
 86                                                    98 
 87     virtual G4VParticleChange* AtRestDoIt(     <<  99      //  no operation in  AtRestDoIt
 88                              const G4Track& ,  << 100       virtual G4VParticleChange* AtRestDoIt(
 89                              const G4Step&     << 101            const G4Track& ,
 90                             ) { return 0; }    << 102            const G4Step&
 91       // No operation in AtRestDoIt            << 103           ) {return 0;};
 92                                                   104  
 93   protected:                                   << 105  protected:// with description
 94                                                << 106      virtual G4double GetMeanFreePath(const G4Track& aTrack,
 95     virtual G4double GetMeanFreePath(const G4T << 107                              G4double   previousStepSize,
 96                              G4double previous << 108                              G4ForceCondition* condition
 97                              G4ForceCondition* << 109                                                                )=0;
 98       // Calculates from the macroscopic cross << 110       //  Calculates from the macroscopic cross section a mean
 99       // free path, the value is returned in u << 111       //  free path, the value is returned in units of distance.
100                                                   112 
                                                   >> 113   protected:// with description
101     virtual G4double GetContinuousStepLimit(co    114     virtual G4double GetContinuousStepLimit(const G4Track& aTrack,
102                              G4double  previou    115                              G4double  previousStepSize,
103                              G4double  current    116                              G4double  currentMinimumStep,
104                              G4double& current << 117            G4double& currentSafety
                                                   >> 118                                                              )=0;
                                                   >> 119   private:
                                                   >> 120     // this is the returnd value of  G4GPILSelection in 
                                                   >> 121     // the arguments of AlongStepGPIL()
                                                   >> 122     G4GPILSelection  valueGPILSelection;
                                                   >> 123 
                                                   >> 124   protected:// with description
                                                   >> 125     // these two methods are set/get methods for valueGPILSelection
                                                   >> 126     void SetGPILSelection(G4GPILSelection selection)
                                                   >> 127     { valueGPILSelection = selection;};
                                                   >> 128 
                                                   >> 129     G4GPILSelection GetGPILSelection() const{return valueGPILSelection;};
                                                   >> 130 
                                                   >> 131    private:
                                                   >> 132   // hide default constructor and assignment operator as private 
                                                   >> 133       G4VContinuousDiscreteProcess();
                                                   >> 134       G4VContinuousDiscreteProcess & operator=(const G4VContinuousDiscreteProcess &right);
105                                                   135 
106     inline void SetGPILSelection(G4GPILSelecti << 136 };
107       { valueGPILSelection = selection; }      << 137 // -----------------------------------------
                                                   >> 138 //  inlined function members implementation
                                                   >> 139 // -----------------------------------------
                                                   >> 140 #include "G4Step.hh"
                                                   >> 141 #include "G4Track.hh"
                                                   >> 142 #include "G4MaterialTable.hh"
                                                   >> 143 #include "G4VParticleChange.hh"
108                                                   144 
109     inline G4GPILSelection GetGPILSelection()  << 145 inline G4double G4VContinuousDiscreteProcess::PostStepGetPhysicalInteractionLength(
110       { return valueGPILSelection; }           << 146                              const G4Track& track,
                                                   >> 147            G4double   previousStepSize,
                                                   >> 148            G4ForceCondition* condition
                                                   >> 149           )
                                                   >> 150 {
                                                   >> 151   if ( (previousStepSize <=0.0) || (theNumberOfInteractionLengthLeft<=0.0)) {
                                                   >> 152     // beggining of tracking (or just after DoIt of this process)
                                                   >> 153     ResetNumberOfInteractionLengthLeft();
                                                   >> 154   } else {
                                                   >> 155     // subtract NumberOfInteractionLengthLeft 
                                                   >> 156     SubtractNumberOfInteractionLengthLeft(previousStepSize);
                                                   >> 157     if(theNumberOfInteractionLengthLeft<0.)
                                                   >> 158        theNumberOfInteractionLengthLeft=perMillion;        
                                                   >> 159   }
                                                   >> 160 
                                                   >> 161   // condition is set to "Not Forced"
                                                   >> 162   *condition = NotForced;
                                                   >> 163 
                                                   >> 164   // get mean free path
                                                   >> 165   currentInteractionLength = GetMeanFreePath(track, previousStepSize, condition);
                                                   >> 166 
                                                   >> 167 #ifdef G4VERBOSE
                                                   >> 168   if ((currentInteractionLength <=0.0) || (verboseLevel>2)){
                                                   >> 169     G4cout << "G4VContinuousDiscreteProcess::PostStepGetPhysicalInteractionLength ";
                                                   >> 170     G4cout << "[ " << GetProcessName() << "]" <<G4endl;
                                                   >> 171     track.GetDynamicParticle()->DumpInfo();
                                                   >> 172     G4cout << " in Material  " << track.GetMaterial()->GetName() <<G4endl;
                                                   >> 173     G4cout << "MeanFreePath = " << currentInteractionLength/cm << "[cm]" <<G4endl;
                                                   >> 174   }
                                                   >> 175 #endif
111                                                   176 
112   private:                                     << 177   G4double value;
                                                   >> 178   if (currentInteractionLength <DBL_MAX) {
                                                   >> 179     value = theNumberOfInteractionLengthLeft * currentInteractionLength;
                                                   >> 180   } else {
                                                   >> 181     value = DBL_MAX;
                                                   >> 182   }
                                                   >> 183 #ifdef G4VERBOSE
                                                   >> 184   if (verboseLevel>1){
                                                   >> 185     G4cout << "G4VContinuousDiscreteProcess::PostStepGetPhysicalInteractionLength ";
                                                   >> 186     G4cout << "[ " << GetProcessName() << "]" <<G4endl;
                                                   >> 187     track.GetDynamicParticle()->DumpInfo();
                                                   >> 188     G4cout << " in Material  " <<  track.GetMaterial()->GetName() <<G4endl;
                                                   >> 189     G4cout << "InteractionLength= " << value/cm <<"[cm] " <<G4endl;
                                                   >> 190   }
                                                   >> 191 #endif
                                                   >> 192   return value;
                                                   >> 193 }
113                                                   194 
114     G4VContinuousDiscreteProcess();            << 
115       // Hidden default constructor            << 
116                                                   195 
117     G4GPILSelection valueGPILSelection = Candi << 196 inline G4VParticleChange* G4VContinuousDiscreteProcess::PostStepDoIt(
118       // This is the returned value of G4GPILS << 197            const G4Track& ,
119       // the arguments of AlongStepGPIL()      << 198            const G4Step&
120 };                                             << 199           )
                                                   >> 200 { 
                                                   >> 201 //  clear  NumberOfInteractionLengthLeft
                                                   >> 202     ClearNumberOfInteractionLengthLeft();
                                                   >> 203     return pParticleChange;
                                                   >> 204 }
                                                   >> 205 
                                                   >> 206 inline G4VParticleChange* G4VContinuousDiscreteProcess::AlongStepDoIt(
                                                   >> 207            const G4Track& ,
                                                   >> 208            const G4Step&
                                                   >> 209           )
                                                   >> 210 { 
                                                   >> 211 //  clear  NumberOfInteractionLengthLeft
                                                   >> 212     ClearNumberOfInteractionLengthLeft();
                                                   >> 213     return pParticleChange;
                                                   >> 214 }
                                                   >> 215 
                                                   >> 216 inline G4double G4VContinuousDiscreteProcess::AlongStepGetPhysicalInteractionLength(
                                                   >> 217                              const G4Track& track,
                                                   >> 218            G4double previousStepSize,
                                                   >> 219            G4double currentMinimumStep,
                                                   >> 220            G4double& currentSafety,
                                                   >> 221                              G4GPILSelection* selection
                                                   >> 222           )
                                                   >> 223 {
                                                   >> 224   // GPILSelection is set to defaule value of CandidateForSelection
                                                   >> 225   valueGPILSelection = CandidateForSelection;
                                                   >> 226 
                                                   >> 227   // get Step limit proposed by the process
                                                   >> 228   G4double steplength = GetContinuousStepLimit(track,previousStepSize,currentMinimumStep, currentSafety);
                                                   >> 229 
                                                   >> 230   // set return value for G4GPILSelection
                                                   >> 231   *selection = valueGPILSelection;
                                                   >> 232 
                                                   >> 233 #ifdef G4VERBOSE
                                                   >> 234    if (verboseLevel>1){
                                                   >> 235     G4cout << "G4VContinuousDiscreteProcess::AlongStepGetPhysicalInteractionLength ";
                                                   >> 236     G4cout << "[ " << GetProcessName() << "]" <<G4endl;
                                                   >> 237     track.GetDynamicParticle()->DumpInfo();
                                                   >> 238     G4cout << " in Material  " <<  track.GetMaterial()->GetName() <<G4endl;
                                                   >> 239     G4cout << "IntractionLength= " << steplength/cm <<"[cm] " <<G4endl;
                                                   >> 240   }
                                                   >> 241 #endif
                                                   >> 242   return  steplength ;
                                                   >> 243 }
121                                                   244 
122 #endif                                            245 #endif
                                                   >> 246 
                                                   >> 247 
                                                   >> 248 
                                                   >> 249 
                                                   >> 250 
                                                   >> 251 
123                                                   252