Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/parameterisation/include/G4FastStep.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/parameterisation/include/G4FastStep.hh (Version 11.3.0) and /processes/parameterisation/include/G4FastStep.hh (Version 10.2)


  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: G4FastStep.hh 68056 2013-03-13 14:44:48Z gcosmo $
 27 //                                                 28 //
 28 //                                             <<  29 // 
 29 //--------------------------------------------     30 //---------------------------------------------------------------
 30 //                                                 31 //
 31 //  G4FastStep.hh                                  32 //  G4FastStep.hh
 32 //                                                 33 //
 33 //  Description:                                   34 //  Description:
 34 //    The G4FastStep class insures a friendly      35 //    The G4FastStep class insures a friendly interface
 35 //    to manage the primary/secondaries final  <<  36 //    to manage the primary/secondaries final state for 
 36 //    Fast Simulation Models. This includes fi     37 //    Fast Simulation Models. This includes final states of parent
 37 //    particle (normalized direction of the mo <<  38 //    particle (normalized direction of the momentum, energy, etc) and 
 38 //    secondary particles generated by the par     39 //    secondary particles generated by the parameterisation.
 39 //                                                 40 //
 40 //    The G4FastStep class acts also as the G4     41 //    The G4FastStep class acts also as the G4ParticleChange
 41 //    for the Fast Simulation Process. So it i <<  42 //    for the Fast Simulation Process. So it inherites from 
 42 //    the G4VParticleChange class and redefine <<  43 //    the G4VParticleChange class and redefines the four virtual 
 43 //    methods :                                    44 //    methods :
 44 //                                                 45 //
 45 //     virtual G4Step* UpdateStepForAtRest(G4S     46 //     virtual G4Step* UpdateStepForAtRest(G4Step* Step);
 46 //     virtual G4Step* UpdateStepForAlongStep(     47 //     virtual G4Step* UpdateStepForAlongStep(G4Step* Step);
 47 //     virtual G4Step* UpdateStepForPostStep(G     48 //     virtual G4Step* UpdateStepForPostStep(G4Step* Step);
 48 //     virtual void Initialize(const G4Track&)     49 //     virtual void Initialize(const G4Track&);
 49 //                                                 50 //
 50 //  History:                                       51 //  History:
 51 //    Oct 97: Verderi && MoraDeFreitas - First     52 //    Oct 97: Verderi && MoraDeFreitas - First Implementation.
 52 //    Dec 97: Verderi - ForceSteppingHitInvoca     53 //    Dec 97: Verderi - ForceSteppingHitInvocation(),
 53 //                      Set/GetTotalEnergyDepo     54 //                      Set/GetTotalEnergyDeposited() methods.
 54 //    Apr 98: MoraDeFreitas - G4FastStep becom     55 //    Apr 98: MoraDeFreitas - G4FastStep becomes the G4ParticleChange
 55 //                      for the Fast Simulatio     56 //                      for the Fast Simulation Process.
 56 //    Nov 04: Verderi - Add ProposeXXX methods     57 //    Nov 04: Verderi - Add ProposeXXX methods. SetXXX ones are kept
 57 //                      for backward compatibi <<  58 //                      for backward compatibility. 
 58 //                                                 59 //
 59 //--------------------------------------------     60 //---------------------------------------------------------------
 60                                                    61 
                                                   >>  62 
 61 #ifndef G4FastStep_h                               63 #ifndef G4FastStep_h
 62 #define G4FastStep_h                               64 #define G4FastStep_h
 63                                                    65 
 64 #include "G4ParticleMomentum.hh"               << 
 65 #include "G4ThreeVector.hh"                    << 
 66 #include "G4ios.hh"                            << 
 67 #include "globals.hh"                              66 #include "globals.hh"
                                                   >>  67 #include "G4ios.hh"
                                                   >>  68 #include "G4ThreeVector.hh"
                                                   >>  69 #include "G4ParticleMomentum.hh"
 68 class G4DynamicParticle;                           70 class G4DynamicParticle;
 69 #include "G4FastTrack.hh"                      << 
 70 #include "G4VParticleChange.hh"                    71 #include "G4VParticleChange.hh"
                                                   >>  72 #include "G4FastTrack.hh"
 71                                                    73 
 72 //-------------------------------------------      74 //-------------------------------------------
 73 //                                                 75 //
 74 //        G4FastStep class                         76 //        G4FastStep class
 75 //                                                 77 //
 76 //-------------------------------------------      78 //-------------------------------------------
 77                                                    79 
 78 // Class Description:                              80 // Class Description:
 79 //  The final state of the particles after par <<  81 //  The final state of the particles after parameterisation has to be returned through a G4FastStep 
 80 //  reference. This final state is described a <<  82 //  reference. This final state is described as "requests" the tracking will apply after your 
 81 //  parameterisation has been invoked.             83 //  parameterisation has been invoked.
 82 //                                                 84 //
 83 //  To facilitate the developers work, changes <<  85 //  To facilitate the developers work, changes of position/normalized direction of the 
 84 //  momentum/polarization can be specified in  <<  86 //  momentum/polarization can be specified in the local coordinate system of the envelope or in the 
 85 //  global one.                                    87 //  global one.
 86 //  The default is local system coordinates.       88 //  The default is local system coordinates.
 87 //                                                 89 //
 88                                                    90 
 89 class G4FastStep : public G4VParticleChange    <<  91 class G4FastStep: public G4VParticleChange
 90 {                                                  92 {
 91   public:  // Without description              <<  93 public: // with Description
 92     //------------------------                 <<  94   void KillPrimaryTrack();
 93     // Constructor/Destructor                  <<  95   // Set the kinetic energy of the primary to zero, and set the "fStopAndKill" signal
 94     //------------------------                 <<  96   // used by the stepping.
 95     G4FastStep() = default;                    <<  97 
 96     ~G4FastStep() override = default;          <<  98   // -- Methods used to change the position, normalized direction of 
 97                                                <<  99   // the momentum, time etc... of the primary.
 98     G4FastStep(const G4FastStep& right) = dele << 100   // .. space and time:
 99     G4FastStep& operator=(const G4FastStep& ri << 101   void ProposePrimaryTrackFinalPosition (const G4ThreeVector &, 
100                                                << 102            G4bool localCoordinates = true);
101     // Set the kinetic energy of the primary t << 103   // Set the primary track final position.
102     // used by the stepping.                   << 104   void     SetPrimaryTrackFinalPosition (const G4ThreeVector &, 
103     void KillPrimaryTrack();                   << 105            G4bool localCoordinates = true);
104                                                << 106   // Set the primary track final position -- maintained for backward compatibility.
105     // -- Methods used to change the position, << 107   
106     // the momentum, time etc... of the primar << 108   
107     // .. space and time:                      << 109   void ProposePrimaryTrackFinalTime (G4double);
108                                                << 110   // Set the primary track final time.
109     // Set the primary track final position.   << 111   void    SetPrimaryTrackFinalTime (G4double);
110     void ProposePrimaryTrackFinalPosition(cons << 112   // Set the primary track final time -- maintained for backward compatibility.
111                                                << 113 
112     // Set the primary track final position -- << 114 
113     [[deprecated("use ProposePrimaryTrackFinal << 115   void ProposePrimaryTrackFinalProperTime (G4double);
114     void SetPrimaryTrackFinalPosition(const G4 << 116   // Set the primary final track Proper Time.
115                                                << 117   void     SetPrimaryTrackFinalProperTime (G4double);
116     // Set the primary track final time.       << 118   // Set the primary final track Proper Time -- maintained for backward compatibility.
117     void ProposePrimaryTrackFinalTime(G4double << 119 
118                                                << 120 
119     // Set the primary track final time -- mai << 121   // .. dynamics:
120     [[deprecated("use ProposePrimaryTrackFinal << 122   void ProposePrimaryTrackFinalMomentumDirection (const G4ThreeVector &, 
121     void SetPrimaryTrackFinalTime(G4double);   << 123               G4bool localCoordinates = true);
122                                                << 124   // Be careful: the Track Final Momentum means the normalized direction 
123     // Set the primary final track Proper Time << 125   // of the momentum!
124     void ProposePrimaryTrackFinalProperTime(G4 << 126   void     SetPrimaryTrackFinalMomentum          (const G4ThreeVector &, 
125                                                << 127               G4bool localCoordinates = true);
126     // Set the primary final track Proper Time << 128   // Set the primary track final momentum -- maintained for backward compatibility. Same as ProposePrimaryTrackMomentumDirection(...)
127     [[deprecated("use ProposePrimaryTrackPrope << 129 
128     void SetPrimaryTrackFinalProperTime(G4doub << 130 
129                                                << 131   void ProposePrimaryTrackFinalKineticEnergy (G4double);
130     // .. dynamics:                            << 132   // Set the primary track final kinetic energy.
131                                                << 133   void     SetPrimaryTrackFinalKineticEnergy (G4double);
132     // Be careful: the Track Final Momentum me << 134   // Set the primary track final kinetic energy-- maintained for backward compatibility.
133     // of the momentum!                        << 135 
134     void ProposePrimaryTrackFinalMomentumDirec << 136 
135                                                << 137   void ProposePrimaryTrackFinalKineticEnergyAndDirection(G4double, 
136                                                << 138                const G4ThreeVector &, 
137     // Set the primary track final momentum -- << 139                G4bool localCoordinates 
138     // ProposePrimaryTrackMomentumDirection(.. << 140                = true);
139     [[deprecated("use ProposePrimaryTrackMomen << 141   // Set the primary track final kinetic energy and direction.
140     void SetPrimaryTrackFinalMomentum(const G4 << 142   void     SetPrimaryTrackFinalKineticEnergyAndDirection(G4double, 
141                                                << 143                const G4ThreeVector &, 
142     // Set the primary track final kinetic ene << 144                G4bool localCoordinates 
143     void ProposePrimaryTrackFinalKineticEnergy << 145                = true);
144                                                << 146   // Set the primary track final kinetic energy and direction -- maintained for backward compatibility.
145     // Set the primary track final kinetic ene << 147 
146     [[deprecated("use ProposePrimaryTrackFinal << 148 
147     void SetPrimaryTrackFinalKineticEnergy(G4d << 149 
148                                                << 150   void ProposePrimaryTrackFinalPolarization(const G4ThreeVector &, 
149     // Set the primary track final kinetic ene << 151               G4bool localCoordinates = true);
150     void ProposePrimaryTrackFinalKineticEnergy << 152   // Set the primary track final polarization.
151                                                << 153   void     SetPrimaryTrackFinalPolarization(const G4ThreeVector &, 
152                                                << 154               G4bool localCoordinates = true);
153     // Set the primary track final kinetic ene << 155   // Set the primary track final polarization.
154     // compatibility.                          << 156 
155     [[deprecated("use ProposePrimaryTrackFinal << 157 
156     void SetPrimaryTrackFinalKineticEnergyAndD << 158   void ProposePrimaryTrackPathLength (G4double);
157                                                << 159   // Set the true path length of the primary track during the step.
158                                                << 160   void     SetPrimaryTrackPathLength (G4double);
159     // Set the primary track final polarizatio << 161   // Set the true path length of the primary track during the step -- maintained for backward compatibility.
160     void ProposePrimaryTrackFinalPolarization( << 162 
161                                                << 163   void ProposePrimaryTrackFinalEventBiasingWeight (G4double);
162     // Set the primary track final polarizatio << 164   // Set the weight applied for event biasing mechanism.
163     [[deprecated("use ProposePrimaryTrackFinal << 165   void     SetPrimaryTrackFinalEventBiasingWeight (G4double);
164     void SetPrimaryTrackFinalPolarization(cons << 166   // Set the weight applied for event biasing mechanism -- kept for backward compatibility.
165                                                << 167 
166     // Set the true path length of the primary << 168   // ------------------------------
167     void ProposePrimaryTrackPathLength(G4doubl << 169   // -- Management of secondaries:
168                                                << 170   // ------------------------------
169     // Set the true path length of the primary << 171 
170     // compatibility.                          << 172   // ----------------------------------------------------
171     [[deprecated("use ProposePrimaryTrackPathL << 173   // -- The creation of secondaries is Done in two steps:
172     void SetPrimaryTrackPathLength(G4double);  << 174   // --      1) Give the total number of secondaries
173                                                << 175   // --         that the FastStep returns
174     // Set the weight applied for event biasin << 176   // --         to the tracking using:
175     void ProposePrimaryTrackFinalEventBiasingW << 177   // --         SetNumberOfSecondaryTracks()
176                                                << 178   // --
177     // Set the weight applied for event biasin << 179   // --      2) Invoke the CreateSecondaryTrack() method
178     [[deprecated("use ProposePrimaryTrackFinal << 180   // --         to create one secondary at each time.
179     void SetPrimaryTrackFinalEventBiasingWeigh << 181   // ----------------------------------------------------
180                                                << 182 
181     // ------------------------------          << 183   // -- Total Number of secondaries to be created,
182     // -- Management of secondaries:           << 184   // -- (to be called first)
183     // ------------------------------          << 185   void SetNumberOfSecondaryTracks(G4int);
184                                                << 186   // Set the total number of secondaries that will be created.
185     // --------------------------------------- << 187 
186     // -- The creation of secondaries is Done  << 188   // -- Number of secondaries effectively stored:
187     // --      1) Give the total number of sec << 189   // -- (incremented at each CreateSecondaryTrack()
188     // --         that the FastStep returns    << 190   // -- call)
189     // --         to the tracking using:       << 191   G4int GetNumberOfSecondaryTracks();
190     // --         SetNumberOfSecondaryTracks() << 192   // Returns the number of secondaries effectively stored.
191     // --                                      << 193 
192     // --      2) Invoke the CreateSecondaryTr << 194   // -- Create a secondary: the arguments are:
193     // --         to create one secondary at e << 195   // --     * G4DynamicsParticle: see header file, many constructors exist
194     // --------------------------------------- << 196   // --                           (allow to set particle type + energy + 
195                                                << 197   // -- the normalized direction of momentum...)
196     // Set the total number of secondaries tha << 198   // --     * G4ThreeVector     : Polarization (not in G4ParticleChange constructor)
197     // -- Total Number of secondaries to be cr << 199   // --     * G4ThreeVector     : Position
198     // -- (to be called first)                 << 200   // --     * G4double          : Time
199     void SetNumberOfSecondaryTracks(G4int);    << 201   // --     * G4bool            : says if Position/Momentum are given in the
200                                                << 202   // --                           local coordinate system (true by default)
201     // Returns the number of secondaries effec << 203   // -- Returned value: pointer to the track created.
202     // -- Number of secondaries effectively st << 204   G4Track* CreateSecondaryTrack(const G4DynamicParticle&,
203     // -- (incremented at each CreateSecondary << 205         G4ThreeVector,
204     // -- call)                                << 206         G4ThreeVector,
205     G4int GetNumberOfSecondaryTracks();        << 207         G4double,
206                                                << 208         G4bool localCoordinates=true);
207     // -- Create a secondary: the arguments ar << 209   // Create a secondary. The arguments are:
208     // --     * G4DynamicsParticle: see header << 210   // 
209     // --                           (allow to  << 211   //   G4DynamicsParticle: see the G4DynamicsParticle reference, many constructors exist
210     // -- the normalized direction of momentum << 212   //                       (allow to set particle type + energy + the normalized direction of 
211     // --     * G4ThreeVector     : Polarizati << 213   //                       momentum...);
212     // --     * G4ThreeVector     : Position   << 214   //   G4ThreeVector     : Polarization;
213     // --     * G4double          : Time       << 215   //   G4ThreeVector     : Position;
214     // --     * G4bool            : says if Po << 216   //   G4double          : Time;
215     // --                           local coor << 217   //   G4bool            : says if Position/Momentum are given in the local envelope coordinate 
216     // -- Returned value: pointer to the track << 218   //                       system (true by default).
217     G4Track* CreateSecondaryTrack(const G4Dyna << 219   //
218                                   G4bool local << 220   // Returned value: pointer to the track created.
219                                                << 221   //
220     //-- Create a secondary: the difference wi << 222   
221     //-- is that the Polarization is not given << 223   //-- Create a secondary: the difference with he above declaration
222     //-- in the G4DynamicParticle.             << 224   //-- is that the Polarization is not given and is assumed already set
223     //-- Returned value: pointer to the track  << 225   //-- in the G4DynamicParticle.
224     G4Track* CreateSecondaryTrack(const G4Dyna << 226   //-- Returned value: pointer to the track created
225                                   G4bool local << 227   G4Track* CreateSecondaryTrack(const G4DynamicParticle&,
226                                                << 228         G4ThreeVector,
227     // Returns a pointer on the i-th secondary << 229         G4double,
228     G4Track* GetSecondaryTrack(G4int);         << 230         G4bool localCoordinates=true);
229                                                << 231   //  Create a secondary. The difference with he above declaration is that the Polarization is not 
230     //---------------------------------------- << 232   //  given and is assumed already set in the G4DynamicParticle.
231     //                                         << 233   //
232     //   Total energy deposit in the "fast Ste << 234   //  Returned value: pointer to the track created
233     //   (a default should be provided in futu << 235 
234     //    which can be:                        << 236   
235     //      delta energy of primary -          << 237 
236     //      energy of the secondaries)         << 238   G4Track* GetSecondaryTrack(G4int);
237     //   This allow the user to Store a consis << 239   // Returns a pointer on the i-th secondary track created.
238     //   information in the G4Trajectory.      << 240 
239     //                                         << 241   //------------------------------------------------
240     //---------------------------------------- << 242   //
241     // Set the total energy deposited.         << 243   //   Total energy deposit in the "fast Step"
242     void ProposeTotalEnergyDeposited(G4double  << 244   //   (a default should be provided in future,
243                                                << 245   //    which can be:
244     // Set the total energy deposited -- kept  << 246   //      delta energy of primary -
245     // It should be the delta energy of primar << 247   //      energy of the secondaries)
246     [[deprecated("use ProposeTotalEnergyDeposi << 248   //   This allow the user to Store a consistent
247     void SetTotalEnergyDeposited(G4double anEn << 249   //   information in the G4Trajectory.
248                                                << 250   //
249     // Returns the total energy deposited.     << 251   //------------------------------------------------
250     G4double GetTotalEnergyDeposited() const;  << 252    void ProposeTotalEnergyDeposited(G4double anEnergyPart);
251                                                << 253   // Set the total energy deposited.
252     // Control of the stepping manager Hit inv << 254   void      SetTotalEnergyDeposited(G4double anEnergyPart);
253     //                                         << 255   // Set the total energy deposited -- kept for backward compatibility.
254     // In a usual parameterisation, the contro << 256   // It should be the delta energy of primary less the energy of the secondaries.
255     // responsability in his G4VFastSimulation << 257 
256     //                                         << 258   G4double GetTotalEnergyDeposited() const;
257     // However, in the particular case the G4F << 259   // Returns the total energy deposited.
258     // replacement only (ie replaces all the * << 260 
259     // G4Step), the user can delegate to the G << 261   void ForceSteppingHitInvocation();
260     // the Hit()method of the current sensitiv << 262   // Control of the stepping manager Hit invocation.
261     //                                         << 263   //
262     // By default, the G4SteppingManager is as << 264   // In a usual parameterisation, the control of the hits production is under the user
263     // parameterisation is invoked.            << 265   // responsability in his G4VFastSimulationModel (he generally produces several hits at once.)
264     void ForceSteppingHitInvocation();         << 266   //
265                                                << 267   // However, in the particular case the G4FastSimulation user's model acts as the physics
266     // ======================================= << 268   // replacement only (ie replaces all the ***DoIt() and leads to the construction of a meaningful
267     // Stepping interface.                     << 269   // G4Step), the user can delegate to the G4SteppingManager the responsability to invoke
268     // ======================================= << 270   // the Hit()method of the current sensitive if any.
269     // --- the following methods are for updat << 271   //
270     // Return the pointer to the G4Step after  << 272   // By default, the G4SteppingManager is asked to NOT invoke this Hit() method when parameterisation
271     // by using final state information of the << 273   // is invoked.
272     //                                         << 274   //   
273     // The Fast Simulation Mechanism doesn't c << 275 
274     // state on the AlongDoIt loop, so the def << 276 
275     // virtual G4Step* UpdateStepForAlongStep( << 277 public: // Without description
276                                                << 278   //=======================================================
277     G4Step* UpdateStepForAtRest(G4Step* Step)  << 279   // Implementation section and kernel interfaces.
278     G4Step* UpdateStepForPostStep(G4Step* Step << 280   //=======================================================
279                                                << 281   //------------------------
280     // A Model gives the final state of the pa << 282   // Constructor/Destructor
281     // based on information of G4FastTrack. So << 283   //------------------------
282     // Initialize method is an interface to th << 284   G4FastStep();
283     // G4FastSimulationManager to Initialize t << 285   virtual ~G4FastStep();
284     // G4FastStep.                             << 286   
285                                                << 287   // equal/unequal operator
286     void Initialize(const G4FastTrack&);       << 288   G4bool operator==(const G4FastStep &right) const;
287                                                << 289   G4bool operator!=(const G4FastStep &right) const;
288     // for Debug                               << 290 
289     void DumpInfo() const override;            << 291 protected:
290     G4bool CheckIt(const G4Track&) override;   << 292   // hide copy constructor and assignment operator as protected
291                                                << 293   G4FastStep (const G4FastStep &right);
292   private:                                     << 294   G4FastStep & operator= (const G4FastStep &right);
293     //======================================== << 295 
294     // Private Internal methods (implementatio << 296 public:
295     //======================================== << 297   // ===============================================
296                                                << 298   // Stepping interface.
297     // G4FastStep should never be Initialized  << 299   // ===============================================
298     // but we must define it to avoid compiler << 300   // --- the following methods are for updating G4Step -----   
299     void Initialize(const G4Track&) override;  << 301   // Return the pointer to the G4Step after updating the Step information
300                                                << 302   // by using final state information of the track given by a Model.
301     //  -- Utility functions --                << 303   //
302     //--- methods to keep information of the f << 304   // The Fast Simulation Mechanism doesn't change the track's final 
303     //  IMPORTANT NOTE: Although the name of t << 305   // state on the AlongDoIt loop, so the default one all we need.
304     //   "Change", what it stores (and returns << 306   //virtual G4Step* UpdateStepForAlongStep(G4Step* Step);
305     //   values of the Position, the normalize << 307 
306     //   etc.                                  << 308   G4Step* UpdateStepForAtRest(G4Step* Step);
307                                                << 309   G4Step* UpdateStepForPostStep(G4Step* Step);
308     // Set theMomentumChange vector: it is the << 310 
309     // direction.                              << 311   // A Model gives the final state of the particle 
310     void SetMomentumChange(G4double Px, G4doub << 312   // based on information of G4FastTrack. So the
311     void SetMomentumChange(const G4ThreeVector << 313   // Initialize method is an interface to the 
312                                                << 314   // G4FastSimulationManager to Initialize the 
313     //======================================== << 315   // G4FastStep.
314     // Data members.                           << 316 
315     //======================================== << 317   void Initialize(const G4FastTrack&);
316     //  theMomentumChange is the vector contai << 318 
317     //  direction after the invoked process. T << 319 private:
318     //  of the momentum direction of the parti << 320   //===================================================
319     //  The responsibility to apply the change << 321   // Private Internal methods (implementation).
320     //  which invoked the process.             << 322   //===================================================
321     G4ParticleMomentum theMomentumChange;      << 323 
322                                                << 324   // G4FastStep should never be Initialized in this way
323     //  The changed (final) polarization of a  << 325   // but we must define it to avoid compiler warnings.
324     G4ThreeVector thePolarizationChange;       << 326   void Initialize(const G4Track&);
325                                                << 327 
326     //  The final kinetic energy of the curren << 328   //  -- Utility functions --
327     G4double theEnergyChange = 0.0;            << 329   //--- methods to keep information of the final state--
328                                                << 330   //  IMPORTANT NOTE: Although the name of the class and methods are
329     //  The changed (final) position of a give << 331   //   "Change", what it stores (and returns in get) are the "FINAL" 
330     G4ThreeVector thePositionChange;           << 332   //   values of the Position, the normalized direction of Momentum, 
331                                                << 333   //   etc.
332     //  The changed (final) global time of a g << 334   
333     G4double theTimeChange = 0.0;              << 335   // Set theMomentumChange vector: it is the final unitary momentum 
334                                                << 336   // direction.
335     //  The changed (final) proper time of a g << 337   void SetMomentumChange(G4double Px, G4double Py, G4double Pz);
336     G4double theProperTimeChange = 0.0;        << 338   void SetMomentumChange(const G4ThreeVector& Pfinal);
337                                                << 339   
338     // The reference G4FastTrack               << 340   //=====================================================
339     const G4FastTrack* fFastTrack = nullptr;   << 341   // Data members.
340                                                << 342   //=====================================================
341     // weight for event biasing mechanism:     << 343   //  theMomentumChange is the vector containing the final momentum 
342     G4double theWeightChange = 0.0;            << 344   //  direction after the invoked process. The application of the change
                                                   >> 345   //  of the momentum direction of the particle is not Done here.
                                                   >> 346   //  The responsibility to apply the change is up the entity
                                                   >> 347   //  which invoked the process.
                                                   >> 348   G4ParticleMomentum theMomentumChange;
                                                   >> 349 
                                                   >> 350   //  The changed (final) polarization of a given particle.
                                                   >> 351   G4ThreeVector thePolarizationChange;
                                                   >> 352 
                                                   >> 353   //  The final kinetic energy of the current particle.
                                                   >> 354   G4double theEnergyChange;
                                                   >> 355 
                                                   >> 356   //  The changed (final) position of a given particle.
                                                   >> 357   G4ThreeVector thePositionChange;
                                                   >> 358 
                                                   >> 359   //  The changed (final) global time of a given particle.
                                                   >> 360   G4double theTimeChange;
                                                   >> 361 
                                                   >> 362   //  The changed (final) proper time of a given particle.
                                                   >> 363   G4double theProperTimeChange;
                                                   >> 364 
                                                   >> 365   // The reference G4FastTrack
                                                   >> 366   const G4FastTrack* fFastTrack;
                                                   >> 367 
                                                   >> 368   // weight for event biasing mechanism:
                                                   >> 369   G4double theWeightChange;
                                                   >> 370 
                                                   >> 371 
                                                   >> 372 public:
                                                   >> 373   // for Debug 
                                                   >> 374   void DumpInfo() const;
                                                   >> 375   G4bool CheckIt(const G4Track&);
343 };                                                376 };
344                                                   377 
345 //********************************************    378 //*******************************************************************
346 //                                                379 //
347 //  Inline functions                              380 //  Inline functions
348 //                                                381 //
349 //********************************************    382 //*******************************************************************
350                                                   383 
351 #include "G4FastStep.icc"                         384 #include "G4FastStep.icc"
352                                                   385 
353 #endif                                            386 #endif
354                                                   387