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

Timestamp:

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

Author:

garnier

Message:

update

File:

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

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

@@ -24 +24 @@
 // ********************************************************************
 //
-// $Id: G4PAIModel.hh,v 1.22 2009/02/19 19:17:50 vnivanch Exp $
-// GEANT4 tag $Name: geant4-09-02-ref-02 $
 //
 // -------------------------------------------------------------------
@@ -77 +75 @@
   virtual void Initialise(const G4ParticleDefinition*, const G4DataVector&);
 
-  virtual void InitialiseMe(const G4ParticleDefinition*);
-
-  virtual G4double ComputeDEDXPerVolume(const G4Material*,
+  virtual void InitialiseMe(const G4ParticleDefinition*) {};
+
+  virtual G4double ComputeDEDX(const G4MaterialCutsCouple*,
                                const G4ParticleDefinition*,
                                G4double kineticEnergy,
                                G4double cutEnergy);
 
-  virtual G4double CrossSectionPerVolume(const G4Material*,
+  virtual G4double CrossSection(const G4MaterialCutsCouple*,
                                 const G4ParticleDefinition*,
                                 G4double kineticEnergy,
@@ -120 +118 @@
 
   void SetVerboseLevel(G4int verbose){fVerbose=verbose;};
+
+
 
 protected:
@@ -192 +192 @@
 };
 
+/////////////////////////////////////////////////////////////////////
+
+inline G4double G4PAIModel::MaxSecondaryEnergy( const G4ParticleDefinition* p,
+                                                      G4double kinEnergy) 
+{
+  G4double tmax = kinEnergy;
+  if(p == fElectron) tmax *= 0.5;
+  else if(p != fPositron) { 
+    G4double mass = p->GetPDGMass();
+    G4double ratio= electron_mass_c2/mass;
+    G4double gamma= kinEnergy/mass + 1.0;
+    tmax = 2.0*electron_mass_c2*(gamma*gamma - 1.) /
+                  (1. + 2.0*gamma*ratio + ratio*ratio);
+  }
+  return tmax;
+}
+
+///////////////////////////////////////////////////////////////
+
+inline  void G4PAIModel::DefineForRegion(const G4Region* r) 
+{
+  fPAIRegionVector.push_back(r);
+}
+
 #endif