Changeset 961 for trunk/source/processes/electromagnetic/standard/src/G4MultipleScattering71.cc

Timestamp:

Apr 6, 2009, 12:21:12 PM (15 years ago)

Author:

garnier

Message:

update processes

File:

• trunk/source/processes/electromagnetic/standard/src/G4MultipleScattering71.cc (modified) (3 diffs)

trunk/source/processes/electromagnetic/standard/src/G4MultipleScattering71.cc

@@ -24 +24 @@
 // ********************************************************************
 //
-// $Id: G4MultipleScattering71.cc,v 1.4 2006/10/16 15:26:49 vnivanch Exp $
-// GEANT4 tag $Name:  $
+// $Id: G4MultipleScattering71.cc,v 1.5 2008/07/16 11:27:41 vnivanch Exp $
+// GEANT4 tag $Name: geant4-09-02-ref-02 $
 //
 // -----------------------------------------------------------------------------
@@ -76 +76 @@
 
 #include "G4MultipleScattering71.hh"
+#include "G4LossTableManager.hh"
 #include "G4MscModel71.hh"
 
@@ -176 +177 @@
 }
 
+//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
+
+G4double G4MultipleScattering71::GetContinuousStepLimit(
+                                          const G4Track& track,
+                                                G4double,
+                                                G4double currentMinimalStep,
+                                                G4double&)
+{
+  DefineMaterial(track.GetMaterialCutsCouple());
+  const G4MaterialCutsCouple* couple = CurrentMaterialCutsCouple();
+  G4double e = track.GetKineticEnergy();
+  model = dynamic_cast<G4MscModel71*>(SelectModel(e));
+  const G4ParticleDefinition* p = track.GetDefinition();
+  G4double lambda0 = GetLambda(p, e);
+  range =  G4LossTableManager::Instance()->GetRangeFromRestricteDEDX(p,e,couple);
+  if(range < currentMinimalStep) currentMinimalStep = range;
+  truePathLength = TruePathLengthLimit(track,lambda0,currentMinimalStep);
+  //  G4cout << "StepLimit: tpl= " << truePathLength << " lambda0= "
+  //       << lambda0 << " range= " << currentRange
+  //       << " currentMinStep= " << currentMinimalStep << G4endl;
+  if (truePathLength < currentMinimalStep) valueGPILSelectionMSC = CandidateForSelection;
+  geomPathLength = model->GeomPathLength(LambdaTable(),couple,
+           p,e,lambda0,range,truePathLength);
+  if(geomPathLength > lambda0) geomPathLength = lambda0;
+  return geomPathLength;
+}
+
 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......