Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/management/src/G4VRestContinuousDiscreteProcess.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 ]

  1 //
  2 // ********************************************************************
  3 // * License and Disclaimer                                           *
  4 // *                                                                  *
  5 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
  6 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
  7 // * conditions of the Geant4 Software License,  included in the file *
  8 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
  9 // * include a list of copyright holders.                             *
 10 // *                                                                  *
 11 // * Neither the authors of this software system, nor their employing *
 12 // * institutes,nor the agencies providing financial support for this *
 13 // * work  make  any representation or  warranty, express or implied, *
 14 // * regarding  this  software system or assume any liability for its *
 15 // * use.  Please see the license in the file  LICENSE  and URL above *
 16 // * for the full disclaimer and the limitation of liability.         *
 17 // *                                                                  *
 18 // * This  code  implementation is the result of  the  scientific and *
 19 // * technical work of the GEANT4 collaboration.                      *
 20 // * By using,  copying,  modifying or  distributing the software (or *
 21 // * any work based  on the software)  you  agree  to acknowledge its *
 22 // * use  in  resulting  scientific  publications,  and indicate your *
 23 // * acceptance of all terms of the Geant4 Software license.          *
 24 // ********************************************************************
 25 //
 26 // G4VRestContinuousDiscreteProcess
 27 // 
 28 // Authors:
 29 // - 2 December 1995, G.Cosmo - First implementation, based on object model
 30 // - 8 January 1997, H.Kurashige - New Physics scheme
 31 // --------------------------------------------------------------------
 32 
 33 #include "G4VRestContinuousDiscreteProcess.hh"
 34 #include "G4SystemOfUnits.hh"
 35 
 36 #include "G4Step.hh"
 37 #include "G4Track.hh"
 38 #include "G4MaterialTable.hh"
 39 #include "G4VParticleChange.hh"
 40 
 41 // --------------------------------------------------------------------
 42 G4VRestContinuousDiscreteProcess::G4VRestContinuousDiscreteProcess()
 43   : G4VProcess("No Name Discrete Process")
 44 {
 45   G4Exception("G4VRestContinuousDiscreteProcess::G4VRestContinuousDiscreteProcess()",
 46               "ProcMan102", JustWarning, "Default constructor is called");
 47 }
 48 
 49 // --------------------------------------------------------------------
 50 G4VRestContinuousDiscreteProcess::
 51 G4VRestContinuousDiscreteProcess(const G4String& aName, G4ProcessType aType)
 52   : G4VProcess(aName, aType)
 53 {
 54 }
 55 
 56 // --------------------------------------------------------------------
 57 G4VRestContinuousDiscreteProcess::~G4VRestContinuousDiscreteProcess()
 58 {
 59 }
 60 
 61 // --------------------------------------------------------------------
 62 G4VRestContinuousDiscreteProcess::
 63 G4VRestContinuousDiscreteProcess(G4VRestContinuousDiscreteProcess& right)
 64   : G4VProcess(right),
 65     valueGPILSelection(right.valueGPILSelection)
 66 {
 67 }
 68 
 69 // --------------------------------------------------------------------
 70 G4double G4VRestContinuousDiscreteProcess::
 71 AtRestGetPhysicalInteractionLength( const G4Track& track,
 72                                     G4ForceCondition* condition )
 73 {
 74   // beginning of tracking 
 75   ResetNumberOfInteractionLengthLeft();
 76 
 77   // condition is set to "Not Forced"
 78   *condition = NotForced;
 79 
 80   // get mean life time
 81   currentInteractionLength = GetMeanLifeTime(track, condition);
 82 
 83   G4double time = (currentInteractionLength < DBL_MAX) ?
 84     theNumberOfInteractionLengthLeft * currentInteractionLength : DBL_MAX;
 85  
 86 #ifdef G4VERBOSE
 87   if ((currentInteractionLength <0.0) || (verboseLevel>2))
 88   {
 89     G4double t = (currentInteractionLength < DBL_MAX) ?
 90       currentInteractionLength/ns : DBL_MAX; 
 91     G4cout << "G4VRestContinuousDiscreteProcess::AtRestGetPhysicalInteractionLength() - ";
 92     G4cout << "[ " << GetProcessName() << "]" << G4endl;
 93     track.GetDynamicParticle()->DumpInfo();
 94     G4cout << " in Material  " << track.GetMaterial()->GetName() << G4endl;
 95     G4cout << "MeanLifeTime = " << t << " [ns]"
 96            << G4endl;
 97   }
 98 #endif
 99 
100   return time;
101 }
102 
103 // --------------------------------------------------------------------
104 G4VParticleChange*
105 G4VRestContinuousDiscreteProcess::AtRestDoIt( const G4Track&,
106                                               const G4Step& )
107 {
108   ClearNumberOfInteractionLengthLeft();
109 
110   return pParticleChange;
111 }
112 
113 // --------------------------------------------------------------------
114 G4double G4VRestContinuousDiscreteProcess::
115 AlongStepGetPhysicalInteractionLength( const G4Track& track,
116                                        G4double previousStepSize,
117                                        G4double currentMinimumStep,
118                                        G4double& currentSafety,
119                                        G4GPILSelection* selection )
120 {
121   // GPILSelection is set to defaule value of CandidateForSelection
122   valueGPILSelection = CandidateForSelection;
123 
124   // get Step limit proposed by the process
125   G4double steplength = GetContinuousStepLimit(track, previousStepSize,
126                                      currentMinimumStep, currentSafety);
127 
128   // set return value for G4GPILSelection
129   *selection = valueGPILSelection;
130 
131 #ifdef G4VERBOSE
132   if (verboseLevel>1)
133   {
134     G4cout << "G4VRestContinuousDiscreteProcess::AlongStepGetPhysicalInteractionLength() - ";
135     G4cout << "[ " << GetProcessName() << "]" << G4endl;
136     track.GetDynamicParticle()->DumpInfo();
137     G4cout << " in Material  " <<  track.GetMaterial()->GetName() << G4endl;
138     G4cout << "IntractionLength= " << steplength/cm <<"[cm] " << G4endl;
139   }
140 #endif
141   return steplength;
142 }
143 
144 // --------------------------------------------------------------------
145 G4VParticleChange*
146 G4VRestContinuousDiscreteProcess::AlongStepDoIt( const G4Track& ,
147                                                  const G4Step& )
148 { 
149   return pParticleChange;
150 }
151 
152 // --------------------------------------------------------------------
153 G4double G4VRestContinuousDiscreteProcess::
154 PostStepGetPhysicalInteractionLength( const G4Track& track,
155                                       G4double previousStepSize,
156                                       G4ForceCondition* condition )
157 {
158   if ( (previousStepSize < 0.0) || (theNumberOfInteractionLengthLeft<=0.0))
159   {
160     // beginning of tracking (or just after DoIt() of this process)
161     ResetNumberOfInteractionLengthLeft();
162   }
163   else if ( previousStepSize > 0.0)
164   {
165     // subtract NumberOfInteractionLengthLeft 
166     SubtractNumberOfInteractionLengthLeft(previousStepSize);
167   }
168   else
169   {
170     // zero step
171     // DO NOTHING
172   }
173 
174   // condition is set to "Not Forced"
175   *condition = NotForced;
176 
177   // get mean free path
178   currentInteractionLength = GetMeanFreePath(track,previousStepSize,condition);
179 
180 
181   G4double value;
182   if (currentInteractionLength < DBL_MAX)
183   {
184     value = theNumberOfInteractionLengthLeft * currentInteractionLength;
185   }
186   else
187   {
188     value = DBL_MAX;
189   }
190 #ifdef G4VERBOSE
191   if (verboseLevel>1)
192   {
193     G4cout << "G4VRestContinuousDiscreteProcess::PostStepGetPhysicalInteractionLength() - ";
194     G4cout << "[ " << GetProcessName() << "]" << G4endl;
195     track.GetDynamicParticle()->DumpInfo();
196     G4cout << " in Material  " <<  track.GetMaterial()->GetName() << G4endl;
197     G4cout << "InteractionLength= " << value/cm <<"[cm] " << G4endl;
198   }
199 #endif
200   return value;
201 }
202 
203 // --------------------------------------------------------------------
204 G4VParticleChange*
205 G4VRestContinuousDiscreteProcess::PostStepDoIt( const G4Track& ,
206                                                 const G4Step& )
207 { 
208   ClearNumberOfInteractionLengthLeft();
209 
210   return pParticleChange;
211 }
212