Changeset 991 for trunk/source/processes/electromagnetic/standard/include/G4BraggModel.hh

Timestamp:

Apr 17, 2009, 12:17:14 PM (15 years ago)

Author:

garnier

Message:

update

File:

• trunk/source/processes/electromagnetic/standard/include/G4BraggModel.hh (modified) (4 diffs)

trunk/source/processes/electromagnetic/standard/include/G4BraggModel.hh

@@ -24 +24 @@
 // ********************************************************************
 //
-// $Id: G4BraggModel.hh,v 1.13 2009/02/20 12:06:37 vnivanch Exp $
-// GEANT4 tag $Name: geant4-09-02-ref-02 $
+// $Id: G4BraggModel.hh,v 1.12 2008/09/14 17:11:48 vnivanch Exp $
+// GEANT4 tag $Name: geant4-09-02 $
 //
 // -------------------------------------------------------------------
@@ -80 +80 @@
   virtual void Initialise(const G4ParticleDefinition*, const G4DataVector&);
 
-  virtual G4double MinEnergyCut(const G4ParticleDefinition*,
-                                const G4MaterialCutsCouple*);
+  G4double MinEnergyCut(const G4ParticleDefinition*,
+                        const G4MaterialCutsCouple*);
                         
   virtual G4double ComputeCrossSectionPerElectron(
@@ -131 +131 @@
 protected:
 
-  virtual G4double MaxSecondaryEnergy(const G4ParticleDefinition*,
-                                      G4double kinEnergy);
+  G4double MaxSecondaryEnergy(const G4ParticleDefinition*,
+                              G4double kinEnergy);
 
 private:
 
-  inline void SetParticle(const G4ParticleDefinition* p);
+  void SetParticle(const G4ParticleDefinition* p);
 
   G4bool HasMaterial(const G4Material* material);
@@ -182 +182 @@
   G4bool   isInitialised;
 };
+
+//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
+
+inline G4double G4BraggModel::MaxSecondaryEnergy(
+                                            const G4ParticleDefinition* pd,
+                                                  G4double kinEnergy)
+{
+  if(pd != particle) SetParticle(pd);
+  G4double tau  = kinEnergy/mass;
+  G4double tmax = 2.0*electron_mass_c2*tau*(tau + 2.) /
+                  (1. + 2.0*(tau + 1.)*ratio + ratio*ratio);
+  return tmax;
+}
 
 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......