Changeset 961 for trunk/source/processes/electromagnetic/highenergy/src/G4hhIonisation.cc

Timestamp:

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

Author:

garnier

Message:

update processes

File:

• trunk/source/processes/electromagnetic/highenergy/src/G4hhIonisation.cc (modified) (6 diffs)

trunk/source/processes/electromagnetic/highenergy/src/G4hhIonisation.cc

@@ -24 +24 @@
 // ********************************************************************
 //
-// $Id: G4hhIonisation.cc,v 1.6 2007/05/22 17:37:30 vnivanch Exp $
-// GEANT4 tag $Name: geant4-09-01-patch-02 $
+// $Id: G4hhIonisation.cc,v 1.9 2009/02/20 16:38:33 vnivanch Exp $
+// GEANT4 tag $Name: geant4-09-02-ref-02 $
 //
 // -------------------------------------------------------------------
@@ -64 +64 @@
     isInitialised(false)
 {
-  minKinEnergy = 0.1*keV;
-  maxKinEnergy = 100.*TeV;
-  SetDEDXBinning(120);
-  SetMinKinEnergy(minKinEnergy);
-  SetMaxKinEnergy(maxKinEnergy);
   SetStepFunction(0.1, 0.1*mm);
   SetVerboseLevel(1);
+  SetProcessSubType(fIonisation);
   mass = 0.0;
   ratio = 0.0;
@@ -79 +75 @@
 G4hhIonisation::~G4hhIonisation()
 {}
+
+//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
+
+G4bool G4hhIonisation::IsApplicable(const G4ParticleDefinition& p)
+{
+  return (p.GetPDGCharge() != 0.0 && p.GetPDGMass() > 100.0*MeV &&
+         !p.IsShortLived());
+}
+
+//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
+
+G4double G4hhIonisation::MinPrimaryEnergy(const G4ParticleDefinition*,
+                                          const G4Material*,
+                                          G4double cut)
+{
+  G4double x = 0.5*cut/electron_mass_c2;
+  G4double y = electron_mass_c2/mass;
+  G4double g = x*y + std::sqrt((1. + x)*(1. + x*y*y));
+  return mass*(g - 1.0);
+}
 
 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
@@ -101 +117 @@
   G4int nm = 1;
 
+  minKinEnergy = MinKinEnergy();
+
   if(eth > minKinEnergy) {
     G4VEmModel* em = new G4BraggNoDeltaModel();
@@ -109 +127 @@
   }
 
-  if(eth < maxKinEnergy) {
+  if(eth < MaxKinEnergy()) {
     G4VEmModel* em1 = new G4BetheBlochNoDeltaModel();
     em1->SetLowEnergyLimit(std::max(eth,minKinEnergy));
-    em1->SetHighEnergyLimit(maxKinEnergy);
+    em1->SetHighEnergyLimit(MaxKinEnergy());
     AddEmModel(nm, em1, flucModel);
   }
 
-  if(verboseLevel>0)
+  if(verboseLevel>1) {
     G4cout << "G4hhIonisation is initialised: Nmodels= " << nm << G4endl;
-
+  }
   isInitialised = true;
 }
@@ -127 +145 @@
 {
   G4cout << "      Delta-ray will not be produced; "
-         << "Bether-Bloch model for E > " << std::max(eth,minKinEnergy)
          << G4endl;
-  if(eth > minKinEnergy) G4cout
-         << "      ICRU49 parametrisation scaled from protons below.";
-  G4cout << G4endl;
 }