Changeset 961 for trunk/source/processes/electromagnetic/muons/src/G4MuPairProduction.cc

Timestamp:

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

Author:

garnier

Message:

update processes

File:

• trunk/source/processes/electromagnetic/muons/src/G4MuPairProduction.cc (modified) (5 diffs)

trunk/source/processes/electromagnetic/muons/src/G4MuPairProduction.cc

@@ -24 +24 @@
 // ********************************************************************
 //
-// $Id: G4MuPairProduction.cc,v 1.48 2007/05/22 17:35:58 vnivanch Exp $
-// GEANT4 tag $Name: geant4-09-01-patch-02 $
+// $Id: G4MuPairProduction.cc,v 1.52 2009/02/20 14:48:16 vnivanch Exp $
+// GEANT4 tag $Name: geant4-09-02-ref-02 $
 //
 // -------------------------------------------------------------------
@@ -89 +89 @@
     lowestKinEnergy(1.*GeV),
     isInitialised(false)
-{}
+{
+  SetProcessSubType(fPairProdByCharged);
+}
 
 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
@@ -98 +100 @@
 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 
-void G4MuPairProduction::InitialiseEnergyLossProcess(const G4ParticleDefinition* part,
-                                                     const G4ParticleDefinition*)
+G4bool G4MuPairProduction::IsApplicable(const G4ParticleDefinition& p)
+{
+  return (p.GetPDGCharge() != 0.0 && p.GetPDGMass() > 10.0*MeV);
+}
+
+//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
+
+G4double G4MuPairProduction::MinPrimaryEnergy(const G4ParticleDefinition*,
+                                              const G4Material*,
+                                              G4double)
+{
+  return lowestKinEnergy;
+}
+
+//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
+
+void G4MuPairProduction::InitialiseEnergyLossProcess(
+                         const G4ParticleDefinition* part,
+                         const G4ParticleDefinition*)
 {
   if (!isInitialised) {
@@ -110 +129 @@
     G4MuPairProductionModel* em = new G4MuPairProductionModel();
     em->SetLowestKineticEnergy(lowestKinEnergy);
-    G4VEmFluctuationModel* fm = new G4UniversalFluctuation();
-    em->SetLowEnergyLimit(0.1*keV);
-    em->SetHighEnergyLimit(100.0*TeV);
+    G4VEmFluctuationModel* fm = 0;
+    em->SetLowEnergyLimit(MinKinEnergy());
+    em->SetHighEnergyLimit(MaxKinEnergy());
     AddEmModel(1, em, fm);
   }
@@ -120 +139 @@
 
 void G4MuPairProduction::PrintInfo()
-{
-  G4cout << "      Parametrised model "
-         << G4endl;
-}
+{}
 
 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....